计算化学公社

标题: 求助：周期性氧化石墨烯计算SCF难收敛 [打印本页]

作者
Author: Yiliuvito 时间: 2025-12-13 12:07
标题: 求助：周期性氧化石墨烯计算SCF难收敛
各位老师好！我现在正在使用CP2K计算周期性GO的REPEAT电荷来确定使用Gromacs的动力学参数，但是现在遇到了几何优化SCF不收敛的问题：我的计算量级为PBE-DZVP-MOLOPT-SR-GTH，使用Multiwfn居中结构，通过Gamma点计算，开启色散矫正。一开始是冻结Z方向的变胞优化，后来选择只优化结构，并分别选择开启OT或Smear，最终结果都难以实现SCF收敛，请问各位老师是我的GO结构有问题么（如果是是不是因为周期性正好截断在了原子上）？请问该如何操作能够实现几何优化的目的呢？感谢各位老师的指导！并附上我最后一次计算的输入文件。

(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 4)

作者
Author: sobereva 时间: 2025-12-13 21:56
初始结构比较差，O-H不可能正好垂直于石墨烯
你给的inp文件里Z方向盒子太小
完全不需要设冻结。
其它方面按下文说的解决
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

作者
Author: 卡开发发 时间: 2025-12-14 15:43
我了解的情况是氧或羟基可能会导致GO或rGO上的碳不共平面，或许可以试试xTB、神经网络势甚至是差不多的力场参数快速优化下结构再交给CP2K看看会不会好些。

作者
Author: ultramanm87 时间: 2025-12-15 13:59
用你的结构试了下，只有第一步128步没法收敛，后面都是没什么问题的。结构优化前期偶尔不收敛一般没啥大问题

作者
Author: Yiliuvito 时间: 2025-12-16 14:34

sobereva 发表于 2025-12-13 21:56
初始结构比较差，O-H不可能正好垂直于石墨烯
你给的inp文件里Z方向盒子太小
完全不需要设冻结。

谢谢社长的帮忙！通过您和其他老师的帮助已经成功计算了原子电荷了，我还想请问：如果是我想要把MD模拟的规模扩大，也就是把GO的片层结构拓展到10nm*10nm，我可以不可以直接使用Multiwfn拓展晶胞然后把计算好的REPEAT电荷手动再分配到已经更新大小的GO结构中呢？还是我还需要计算一个更大的GO的原子电荷？再次感谢社长

作者
Author: Yiliuvito 时间: 2025-12-16 14:34

卡开发发发表于 2025-12-14 15:43
我了解的情况是氧或羟基可能会导致GO或rGO上的碳不共平面，或许可以试试xTB、神经网络势甚至是差不多的力场 ...

谢谢您的回复，我简单先用xTB进行优化，然后再使用PBE进行计算成功获得了REPEAT电荷了，十分感谢！

作者
Author: Yiliuvito 时间: 2025-12-16 14:35

ultramanm87 发表于 2025-12-15 13:59
用你的结构试了下，只有第一步128步没法收敛，后面都是没什么问题的。结构优化前期偶尔不收敛一般没啥大问 ...

谢谢您的帮忙，之前的优化收敛确实步数设的太小，结合楼上所说并提升了SCF的步数，也成功收敛了，十分感谢

作者
Author: sobereva 时间: 2025-12-17 02:42

Yiliuvito 发表于 2025-12-16 14:34
谢谢社长的帮忙！通过您和其他老师的帮助已经成功计算了原子电荷了，我还想请问：如果是我想要把MD模拟的 ...

可以这样

欢迎光临计算化学公社 (http://bbs.keinsci.com/)