计算化学公社

标题: 关于MgO晶胞大小对几何优化的影响 [打印本页]
作者
Author: THY    时间: 2025-7-8 22:40
标题: 关于MgO晶胞大小对几何优化的影响
本帖最后由 THY 于 2025-7-9 10:15 编辑

各位老师，本人使用CP2K计算MgO晶体

发现对于原胞的计算，非常容易收敛，仅需一轮SCF；对于xy方向的2x2超胞也非常容易收敛

但是对于xy方向的3x3及以上的超胞，优化出现了问题，每一轮SCF收敛的非常快，但是能量一直处于振荡的过程。通过xyz轨迹文件发现，内部的原子也一直处于运动状态

请问该如何解决？为什么超胞的大小会影响几何优化的收敛

PS：本人尝试过不同超胞，k点，smearing，OT等方法，均未能解决，似乎超胞大于2x2以后，一切就变得不可控了

附件为：原胞，2x2超胞，3x3超胞
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1)

添加k点后的计算：
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1)

作者
Author: sobereva    时间: 2025-7-9 05:05
先把最严重的硬伤解决再说别的：Z方向只有4.2112埃，却没考虑k点。1x1和2x2的计算的X、Y方向也得考虑k点
作者
Author: THY    时间: 2025-7-9 10:13
sobereva 发表于 2025-7-9 05:05
先把最严重的硬伤解决再说别的：Z方向只有4.2112埃，却没考虑k点。1x1和2x2的计算的X、Y方向也得考虑k点

老师，我对于3x3超胞，添加了2x2x5的k点，发现结果还是一样，内部原子依旧处于一直运动的状态，请问该如何解决？我把加k点的文件附在1楼了
作者
Author: sobereva    时间: 2025-7-9 11:25
参考本论坛里CP2K优化不收敛的讨论，首页全文搜索框搜论坛，回答过很多次



欢迎光临 计算化学公社 (http://ccc.keinsci.com/) Powered by Discuz! X3.3