求助CP2K分子动力学模拟的程序升温和淬火，输入文件应该怎么写，设置哪些关键词呢？

mayya

sobereva 发表于 2024-5-25 04:08
初速度若对应0K，显然所有原子初速度就都为0
0K下体系没有任何动能，精确维持在0K的情况根本没法跑MD

社长，请问我的理解是否正确？如果需要获取从300K升温到3000K的整个过程。
第一步：在300K下跑NVT获取原子的速度。
第二步：将TEMPERATURE设置为3000K，VELOCITY部分改成300K时候的原子速度，CSVR热浴，NPT。
这样是否能获取整个升温过程的数据？我的目的是想要获取不同温度下的分子轨迹用于后续的deepmd训练，这样是否能满足要求？

Uus/pMeC6H4-/キ

mayya 发表于 2024-11-20 10:37
社长，请问我的理解是否正确？如果需要获取从300K升温到3000K的整个过程。
第一步：在300K下跑NVT获取原 ...

这样输出的分子动力学轨迹确实会体现升温，但是升温速度只能靠时间常数间接控制，每个温度下的轨迹长度并不统一，而且升温过程并非平衡的（关于平衡见社长博文http://sobereva.com/627）。

如果“获取不同温度下的分子轨迹用于后续的deepmd训练”要求每个温度都充分平衡且有足够长度的产生相轨迹来当成数据的话，我估计只能设定多个温度分别跑平衡相和产生相了。

顺便留一个其他相关讨论帖的链接