计算化学公社

标题: CP2K做CI-NEB能量如何监控？ [打印本页]

作者
Author: THY 时间: 2025-4-23 17:19
标题: CP2K做CI-NEB能量如何监控？
本帖最后由 THY 于 2025-4-23 17:25 编辑

请问使用CP2K做CI-NEB的时候，如何监测能量收敛过程，从而画出能量收敛曲线，考虑到每个CI-NEB插点数目不同，有没有什么办法可以忽视插点数目，而直接读取整个NEB的能量变化呢？

作者
Author: sobereva 时间: 2025-4-23 18:41
如果你指整条NEB的总能量，用grep把输出文件里的能量那行提出来就完了

作者
Author: 丁越 时间: 2025-4-24 08:50
需要安装gnuplot，运行：

plot for [i=2:9] 'NEB-1.ener' u 1:i w lp t ''.i

复制代码

这里假设你有8个NEB的副本，for循环的变量i范围便设置为了2~9，输出的图中监控了每个副本的能量变化。

作者
Author: THY 时间: 2025-4-25 10:56

sobereva 发表于 2025-4-23 18:41
如果你指整条NEB的总能量，用grep把输出文件里的能量那行提出来就完了

谢谢老师！

作者
Author: THY 时间: 2025-4-25 10:56

丁越发表于 2025-4-24 08:50
需要安装gnuplot，运行：

这里假设你有8个NEB的副本，for循环的变量i范围便设置为了2~9，输出的图中监控 ...

谢谢老师！

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