计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 求助：周期性氧化石墨烯计算SCF难收敛

求助：周期性氧化石墨烯计算SCF难收敛

查看数: 230 | 评论数: 7 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

Yiliuvito

发布时间: 2025-12-13 12:07

正文摘要:

各位老师好！我现在正在使用CP2K计算周期性GO的REPEAT电荷来确定使用Gromacs的动力学参数，但是现在遇到了几何优化SCF不收敛的问题：我的计算量级为PBE-DZVP-MOLOPT-SR-GTH，使用Multiwfn居中结构，通过Gamma点计算 ...

回复 Reply

sobereva 发表于 Post on 2025-12-17 02:42:12

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

可以这样

Yiliuvito 发表于 Post on 2025-12-16 14:35:43

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

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

Yiliuvito 发表于 Post on 2025-12-16 14:34:47

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

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

Yiliuvito 发表于 Post on 2025-12-16 14:34:02

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

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

ultramanm87 发表于 Post on 2025-12-15 13:59:54

用你的结构试了下，只有第一步128步没法收敛，后面都是没什么问题的。结构优化前期偶尔不收敛一般没啥大问题

卡开发发 发表于 Post on 2025-12-14 15:43:41

我了解的情况是氧或羟基可能会导致GO或rGO上的碳不共平面，或许可以试试xTB、神经网络势甚至是差不多的力场参数快速优化下结构再交给CP2K看看会不会好些。

sobereva 发表于 Post on 2025-12-13 21:56:40

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register