计算化学公社

标题: 求助：cp2k跑aimd温度降不下来是为什么？ [打印本页]

作者
Author: ZQing 时间: 2025-3-27 13:12
标题: 求助：cp2k跑aimd温度降不下来是为什么？
本帖最后由 ZQing 于 2025-3-27 13:14 编辑

各位老师好，我在cp2k在跑一个含有NaHCO3、MgCO3、CaCO3的成分的体系，初始结构是用packmol搭建的，NPT系综，温度为900K，压力为200MPa,目前是打算先跑一个eq平衡相，跑2ps,步长1fs,基组用的是DZVP-MOLOPT-SR-GTH，赝势是GTH-PBE，CUTOFF 是300，收敛精度是1.0E-05，具体inp文件放在下面了，但是现在跑了快两天了，大概才跑了500多步，温度也一直在2000K上下，没有降下来，想请教一下各位老师是什么原因呢？

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-3-27 14:25
（这是模拟混合碳酸盐熔盐么，感觉高压下不该这么稀疏啊，盒子尺寸怎么选的）
记得看输出文件检查每一步SCF收敛情况，参考
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

Na和Mg是典型的在此基组/赝势组合下对截断能要求极高的元素，现在CUTOFF才300、收敛限才1.0E-5，这些原子的能量和受力都计算不够准，动力学行为可能有偏差。另外900 K高温下步长1 fs可能也过大了。有关packmol构建初始结构与温度控制的问题，见这帖与这帖；看起来盒子里部分区域阳离子聚集、部分区域阴离子聚集，这样静电排斥明显也会使结构不稳定。

作者
Author: ZQing 时间: 2025-3-27 20:28

Uus/pMeC6H4-/キ发表于 2025-3-27 14:25
（这是模拟混合碳酸盐熔盐么，感觉高压下不该这么稀疏啊，盒子尺寸怎么选的）
记得看输出文件检查每一步SC ...

感谢大佬的解答

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