计算化学公社

标题: CP2K跑AIMD控温控不住如何解决 [打印本页]

作者
Author: 1372133091 时间: 2025-3-31 15:40
标题: CP2K跑AIMD控温控不住如何解决
各位老师好！最近在计算Li7La3Zr2O12这个材料的质子迁移率，参考了公社的几篇帖子之后决定使用CP2K进行分子动力学计算，但是在计算过程中无法将温度稳定在预设温度，例如预设800K实际上只能稳定在500K. 想求助大家有遇到此类情况的解决办法。
输入输出文件已附。

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-3-31 16:04
本帖最后由 Uus/pMeC6H4-/キ于 2025-3-31 16:07 编辑

热浴不建议用Nose-Hoover，改成CSVR吧。而且总模拟时长才1 ps很可能不够平衡的，最好延长并参考http://sobereva.com/627判断何时才真正平衡。另外鉴于氢离子很轻、温度又是800 K高温，1 fs的步长可能偏大了，为使动力学行为更准确最好减小些。再有，UKS和PLUS_U的设置为啥留着呢？

（编辑：我还好奇一点，质子是单独出现在这材料里的么，没有与阴离子杂质或材料的阳离子缺陷之类一起出现使结构保持电中性？以往类似文献研究用的也是带净电荷的周期性体系么，如果是的话对周期镜像之间的静电排斥有什么考虑？）

作者
Author: ljh123 时间: 2025-3-31 17:15

Uus/pMeC6H4-/キ发表于 2025-3-31 16:04
热浴不建议用Nose-Hoover，改成CSVR吧。而且总模拟时长才1 ps很可能不够平衡的，最好延长并参考http://sobe ...

确实很多文献进行AIMD模拟时都没有考虑电荷保持中性，
http://bbs.keinsci.com/thread-33831-1-1.html 中的5楼sob老师也讲到这点

作者
Author: 1372133091 时间: 2025-3-31 22:06

Uus/pMeC6H4-/キ发表于 2025-3-31 16:04
热浴不建议用Nose-Hoover，改成CSVR吧。而且总模拟时长才1 ps很可能不够平衡的，最好延长并参考http://sobe ...

好的，谢谢。因为是想通过AIMD模拟质子在固态电解质中的扩散系数，进而计算出其扩散活化能，因此就这么设置了

作者
Author: sobereva 时间: 2025-4-1 00:24
如今几乎没有任何理由用Nose-Hoover
想控温效果更显著就用较小的TIMECON，但即便用很小的50 fs，1 ps模拟时间长度也完全不够

作者
Author: 1372133091 时间: 2025-4-1 08:55

sobereva 发表于 2025-4-1 00:24
如今几乎没有任何理由用Nose-Hoover
想控温效果更显著就用较小的TIMECON，但即便用很小的50 fs，1 ps模拟 ...

好的，谢谢sob老师

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