求助，CP2k做MD模拟怎么输出受力信息

楚瑶 · 发表于 Post on 2024-3-22 16:17:28

如题，请问各位老师CP2K跑AIMD的时候怎么输出受力信息。
目前Print_Level设置是Low，然后输出有restart文件、log文件、*-pos-1.xyz和*-1.ener文件。我看有人运行后会输出*-frc-1.xyz，想问下这个在输出文件中怎么设置才会输出呢?
谢谢各位老师

sobereva · 发表于 Post on 2024-3-22 16:42:30

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）的ppt：

chen0201 · 发表于 Post on 2024-7-15 15:25:28

本帖最后由 chen0201 于 2024-7-15 15:27 编辑

sobereva 发表于 2024-3-22 16:42
参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）的ppt：

请问老师，打印出来的FORCES包含的是每个原子的xyz受力吗？如果需要计算如图的受力情况，是不是需要加起来呢？

sobereva · 发表于 Post on 2024-7-16 22:24:48

chen0201 发表于 2024-7-15 15:25
请问老师，打印出来的FORCES包含的是每个原子的xyz受力吗？如果需要计算如图的受力情况，是不是需要加起 ...

是
我不知道图中对应什么

chen0201 · 发表于 Post on 2024-7-17 11:35:46

sobereva 发表于 2024-7-16 22:24
是
我不知道图中对应什么

谢谢老师，那个图应该就是总的受力随着坐标X的变化而变化的值

yxzinc · 发表于 Post on 2024-7-21 09:47:05

chen0201 发表于 2024-7-15 15:25
请问老师，打印出来的FORCES包含的是每个原子的xyz受力吗？如果需要计算如图的受力情况，是不是需要加起 ...

这是王阳刚老师的文章吧。这个力根本不是原子受力，而是限制性分子动力学所外加的力。他文章的方法部分应该有介绍，具体cp2k操作sob老师应该也教过

chen0201 · 发表于 Post on 2024-7-25 09:42:28

yxzinc 发表于 2024-7-21 09:47
这是王阳刚老师的文章吧。这个力根本不是原子受力，而是限制性分子动力学所外加的力。他文章的方法部分应 ...

嗯嗯，是的，他这个力是平均力，我已经按照sob老师的方法计算了原子受力

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助，CP2k做MD模拟怎么输出受力信息