求助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训练”要求每个温度都充分平衡且有足够长度的产生相轨迹来当成数据的话，我估计只能设定多个温度分别跑平衡相和产生相了。

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

fangyl

sobereva 发表于 2023-8-5 17:31
如果输入文件里没直接定义初速度（&VELOCITY），说明初始温度将为800K，控温也控在800K
如果初速度对应 ...

请问社长，用CSVR热浴，NVT，如果初始速度由300k弛豫的温度下得到，restart文件改为inp之后，在600k下跑，每个原子的速度方向会改变吗，还是在原来速度方向是乘了个缩放因子，感谢

sobereva

fangyl 发表于 2024-11-28 10:58
请问社长，用CSVR热浴，NVT，如果初始速度由300k弛豫的温度下得到，restart文件改为inp之后，在600k下跑 ...

方向不会改变

fangyl

sobereva 发表于 2024-11-28 12:06
方向不会改变

谢谢老师