计算化学公社

标题: 在CP2K中计算AIMD时，进行结构优化，SCF自洽场收敛后，但结构优化一直不收敛 [打印本页]

作者
Author: 方不惑 时间: 7 day ago
标题: 在CP2K中计算AIMD时，进行结构优化，SCF自洽场收敛后，但结构优化一直不收敛
各位老师好，我在使用CP2K作AIMD应用于电解液的反应模拟，一开始我未进行结构优化直接开始模拟，发现其SCF始终不收敛，即便填写了忽略SCF不收敛，继续作AIMD，但模拟的结构十分诡异；我就寻思先使用PBE-DZVP-MOLOPT-SR进行结构优化，然后再去跑AIMD，但在结构优化的过程中，发现其SCF能够做到收敛，但其结构优化始终无法做到很好的收敛，尝试增大截短能，修改BFGS的设置，但始终无法奏效，请教一下各位老师这种情况该如何解决呢

作者
Author: Uus/pMeC6H4-/キ 时间: 7 day ago
9.5埃的盒子尺寸偏小了，初始结构的周期镜像之间有很多不合理的近距离接触。不知道是用什么软件构建的结构，哪怕是研究某个压强下特定密度/浓度的电解液、对盒子尺寸和分子数有固定要求，填充分子时也应该留一点空间而不是完全顶着盒子来，并且产生的结构应该放到如VMD里显示出周期镜像来检查。

作者
Author: 方不惑 时间: 7 day ago

Uus/pMeC6H4-/キ发表于 2026-4-6 17:29
9.5埃的盒子尺寸偏小了，初始结构的周期镜像之间有很多不合理的近距离接触。不知道是用什么软件构建的结构 ...

老师，我是用的是packmol建立的模型，建的是10*10*35，但是出的盒子总是是正方体，老师我这边去修改一下盒子大小再去作提交

作者
Author: sobereva 时间: 6 day ago

方不惑发表于 2026-4-6 18:41
老师，我是用的是packmol建立的模型，建的是10*10*35，但是出的盒子总是是正方体，老师我这边去修改一下 ...

100% packmol输入文件没写对。packmol当前版本建模时可以考虑周期性，记得这么做，免得和镜像盒子里的原子有不合理接触。
另外，最短边长10埃也偏小，会给动力学造成不可忽视的人为虚假性，最好再增加50%

作者
Author: 方不惑 时间: 6 day ago

sobereva 发表于 2026-4-7 03:24
100% packmol输入文件没写对。packmol当前版本建模时可以考虑周期性，记得这么做，免得和镜像盒子里的原 ...

好的，谢谢sob老师，方才尝试了一下，终于能稳稳的动起来了，十分感谢

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