CP2K测算晶体动力学，630步后跑散了，想问一下是什么设置的不合理导致的呢

e-Lis · 发表于 Post on 2025-5-7 20:37:31

CP2K测算已有晶体的分子动力学，之前在vasp中跑过，前面还挺符合预期的，但是630步后跑散了，前几天培训初学，想问一下是什么设置的不合理导致的呢？是截断能不合理吗，还是系综不合适，因为后续原子全部跑到盒子外边，彻底裂解了。

Uus/pMeC6H4-/キ · 发表于 Post on 2025-5-7 21:41:53

当前体系盒子太小了，一是不考虑非gamma点的k点会使周期性计算不准确，二是分子动力学行为会受限于虚假的周期性而有失真实，至少A和C方向应该扩胞2倍。此外2 fs步长对这种含有氢原子的体系偏大了，考虑减小到1 fs甚至0.5 fs。截断能300 Ry也最好增大些。

e-Lis · 发表于 Post on 2025-5-13 10:16:27

Uus/pMeC6H4-/キ发表于 2025-5-7 21:41
当前体系盒子太小了，一是不考虑非gamma点的k点会使周期性计算不准确，二是分子动力学行为会受限于虚假的周 ...

谢谢您，我去试一下哇

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

[CP2K] CP2K测算晶体动力学，630步后跑散了，想问一下是什么设置的不合理导致的呢