计算化学公社

标题: 求助，CP2k做MD模拟怎么输出受力信息 [打印本页]

作者
Author: 楚瑶 时间: 2024-3-22 16:17
标题: 求助，CP2k做MD模拟怎么输出受力信息
如题，请问各位老师CP2K跑AIMD的时候怎么输出受力信息。
目前Print_Level设置是Low，然后输出有restart文件、log文件、*-pos-1.xyz和*-1.ener文件。我看有人运行后会输出*-frc-1.xyz，想问下这个在输出文件中怎么设置才会输出呢?
谢谢各位老师

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

(, 下载次数 Times of downloads: 17)

作者
Author: chen0201 时间: 2024-7-15 15:25
本帖最后由 chen0201 于 2024-7-15 15:27 编辑

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

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

作者
Author: sobereva 时间: 2024-7-16 22:24

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

是
我不知道图中对应什么

作者
Author: chen0201 时间: 2024-7-17 11:35

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

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

作者
Author: yxzinc 时间: 2024-7-21 09:47

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

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

作者
Author: chen0201 时间: 2024-7-25 09:42

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

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

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