计算化学公社

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

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

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

作者
Author: e-Lis 时间: 2025-5-13 10:16

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

谢谢您，我去试一下哇

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