计算化学公社

标题: CP2K重启NVE系综的AIMD势能出现跳跃 [打印本页]

作者
Author: lmx0413 时间: 2023-7-8 21:50
标题: CP2K重启NVE系综的AIMD势能出现跳跃
各位老师同学，我用CP2K跑了50ps的AIMD，系综是NVE系综，然后再续算了50ps，发现续算的时间点势能出现了跳跃，并且conserved quantity也出现了跳跃。
我在google groups上看到的一些解答说恒温器没有RESTART会出现这样的情况，但是尝试更改RESTART恒温器，但是conserved quantity依然出现跳跃。

下面是我输出的势能随时间变化的图像，部分.ener文件，inp文件和restart文件
这应该怎么解决呢？如有回复，将不胜感激！

作者
Author: sobereva 时间: 2023-7-9 03:50
NVE根本不牵扯恒温器，改这方面没意义
我不知道你是具体怎么续跑的，做了什么操作、对restart文件做了什么修改

&CENTER_COORDINATES T会自动令分子居中，假设之前最后一步分子不是中间位置，重算的时候这个命令导致分子整体有所位移的话，由于基于均匀格点计算的能量并不精确具备平移不变性，有可能这个会造成能量小幅变化。

作者
Author: lmx0413 时间: 2023-7-9 08:27

sobereva 发表于 2023-7-9 03:50
NVE根本不牵扯恒温器，改这方面没意义
我不知道你是具体怎么续跑的，做了什么操作、对restart文件做了什么 ...

好的谢谢老师，我的restart文件没有做任何修改，是直接进行续算的，我只是想在加一段时间的MD

作者
Author: lmx0413 时间: 2023-7-9 08:39

sobereva 发表于 2023-7-9 03:50
NVE根本不牵扯恒温器，改这方面没意义
我不知道你是具体怎么续跑的，做了什么操作、对restart文件做了什么 ...

非常感谢卢老师！！！！我删去这部分后，看起来正常了

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