计算化学公社

标题: 求助：cp2k的restart文件为什么不是最后一步的结果 [打印本页]

作者
Author: wuyeyangguang 时间: 2022-3-10 11:25
标题: 求助：cp2k的restart文件为什么不是最后一步的结果
如题，想用cp2k restarat文件续算，但发现restart文件非实时更新。CP2K每500步会产生一个restart文件，然后从第1000步开始会额外有一个restart文件，此文件会和500步整数倍的restart文件同步更新，并非最后一步的restart文件。
也就是说，假如一个算例跑了2200步，会产生5个restart文件，最新的那个restart文件和第2000步的restart文件是一样的，如果使用该resart文件续算，ener文件内容是不连续的，从2200步接着是2001步，轨迹文件也是相应重复了一小部分，求教如何解决，谢谢。

作者
Author: Eudaimonia 时间: 2022-3-10 14:58
通过更改
&EACH
MD 10 #Frequency of updating last restart file
字段控制restart文件的更新频率，10是每500步输出一次，社长multiwfn产生的输入文件后面都带注释了，如果不是multiwfn产生的文件，为啥不去翻手册呢

另外，你可以手动删掉重复的轨迹

作者
Author: asd9998876 时间: 2022-3-10 21:12
这里也有个问题想要问一下各位老师，重新开始运算时的前几步的能量会明显的高于之前的同一轨迹的能量，这个该怎么取舍呢

作者
Author: wuyeyangguang 时间: 2022-3-11 16:49

Eudaimonia 发表于 2022-3-10 14:58
通过更改
&EACH
MD 10 #Frequency of updating last restart file

谢谢。慌乱的新手不知道还有这种参数来控制，还以为是固定的。

作者
Author: wuyeyangguang 时间: 2022-3-11 16:50

asd9998876 发表于 2022-3-10 21:12
这里也有个问题想要问一下各位老师，重新开始运算时的前几步的能量会明显的高于之前的同一轨迹的能量，这个 ...

这个我也遇到了，同问

作者
Author: chands 时间: 2022-10-3 17:42

Eudaimonia 发表于 2022-3-10 14:58
通过更改
&EACH
MD 10 #Frequency of updating last restart file

请问这个MD的参数与restart文件的更新频率有什么对应关系，我在CP2K官网手册上查不到。

作者
Author: Eudaimonia 时间: 2022-10-3 19:45

chands 发表于 2022-10-3 17:42
请问这个MD的参数与restart文件的更新频率有什么对应关系，我在CP2K官网手册上查不到。

我也不确定有什么对应关系

@jiangjiangna 他可能做过测试，请问有什么对应关系吗？

作者
Author: asd9998876 时间: 2022-10-3 22:06
默认最低是500步一次，10,5,20都是500步，之后的就是1000就是1000步存一次，2000就是2000步存一次，并没有那么多弯弯绕绕

轨迹文件的重复可以手动删去

作者
Author: chands 时间: 2022-10-3 22:33

asd9998876 发表于 2022-10-3 22:06
默认最低是500步一次，10,5,20都是500步，之后的就是1000就是1000步存一次，2000就是2000步存一次，并没有 ...

好的，谢谢！

作者
Author: jiangjiangna 时间: 2022-10-14 15:47
本帖最后由 jiangjiangna 于 2022-10-14 15:50 编辑

我测试过了，后来下载最新版的Multiwfn去生成输入文件，打开输入文件自己去设置的Each,数字填几就是每几步生成一次restart文件，不再遵循旧版本最低500步才更新。@Eudaimonia，@chands

作者
Author: zhangs 时间: 2024-10-19 14:42

Eudaimonia 发表于 2022-3-10 14:58
通过更改
&EACH
MD 10 #Frequency of updating last restart file

请问下最低只能设置500步输出一次吗

作者
Author: 汪杰 时间: 2024-11-29 17:56

jiangjiangna 发表于 2022-10-14 15:47
**** 作者被禁止或删除内容自动屏蔽 ****

请问一下，这计算式分子动力学还是结构优化，我用multiwfn产生输入文件时，似乎并没有这个关键词

作者
Author: Uus/pMeC6H4-/キ 时间: 2024-11-29 18:18
本帖最后由 Uus/pMeC6H4-/キ于 2024-11-29 18:33 编辑

汪杰发表于 2024-11-29 17:56
请问一下，这计算式分子动力学还是结构优化，我用multiwfn产生输入文件时，似乎并没有这个关键词

关键词在&MOTION/&PRINT/&RESTART/&EACH里，是分子动力学就写MD，是结构优化就写GEO_OPT（不变胞）或者CELL_OPT（变胞），后面再写频率。

另外回复这帖：如果先前输入文件里上面关键词的频率写了1，那.restart文件就是每一步实时更新的。

编辑：楼上说每500步更新一次的是&MOTION/&PRINT/&RESTART_HISTORY/&EACH控制的对应特定步数的历史.restart文件，正因为对应特定步数所以才会显得不是实时更新。区分实时和历史两种.restart文件就看有没有加下划线表明步数，没有的对应实时更新的.restart，续算时只用管这个就行。我看在Multiwfn默认的输出文件设置里面已经通过&RESTART_HISTORY OFF关闭历史.restart文件输出了。

作者
Author: 汪杰 时间: 2024-11-29 18:34

Uus/pMeC6H4-/キ发表于 2024-11-29 18:18
关键词在&MOTION/&PRINT/&RESTART/&EACH里，是分子动力学就写MD，是结构优化就写GEO_OPT（不变胞）或者CE ...

好的，谢谢

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