计算化学公社

标题: 熔盐体系进行AIMD的NPT过程时于8000步SCF未收敛 [打印本页]

作者
Author: 以玉名诗 时间: 2026-3-30 15:42
标题: 熔盐体系进行AIMD的NPT过程时于8000步SCF未收敛
本帖最后由以玉名诗于 2026-3-30 15:47 编辑

各位老师好！我的体系是熔盐体系，参考相似（没有完全相同的体系）的体系密度计算了的理论盒子体积，之后使用packmol构建盒子时构建失败。故将盒子边长稍扩大1埃，这才成功得到了混合好的pdb文件。

于Multiwfn创建cp2k的输入文件之后，完成了5000步的NVT系综进行预平衡过程，成功将温度稳定在目标温度。之后使用.restart中的原子坐标构建了输入文件，为30000步的NPT动力学过程，目标是在该温度下得到合适的盒子体积，之后再进行产生相的计算。

起初体系运算还算顺利，但在4000步之后进行每一步所需的时间越来越长，最终该体系于8000步时体系SCF未收敛，计算程序失败。

参考了Sob老师的帖子（http://sobereva.com/665），检查了初始结构、UKS、净电荷和自旋多重度等方面，应当是没问题的。动力学方面EPS_SCF设置的1E-5，也不是很严格。同时步长也是参考文献中的为1 fs。目前仍无法找到解决办法，故请教各位老师我该在哪里进行修改，相关的输入文件和输出文件已上传。

作者
Author: sobereva 时间: 2026-3-31 09:29
开OT的时候绝对不要用诸如MAX_SCF 100这么大的数值，否则outerSCF没法充分起到效果
并且尝试IRAC
这些下文都说了
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

作者
Author: 以玉名诗 时间: 2026-3-31 16:34

sobereva 发表于 2026-3-31 09:29
开OT的时候绝对不要用诸如MAX_SCF 100这么大的数值，否则outerSCF没法充分起到效果
并且尝试IRAC
这些下 ...

谢谢老师值点！

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