cp2k做cineb的能量需要再算单点确定吗

dbxhtw · 发表于 Post on 2026-7-28 15:26:06

求助各位老师：
分5个插点结构计算了CINEB，初末结构已优化并且在neb过程中固定。跑完之后，显示 MAX FORCE 还没完全收敛
*******************************************************************************
RMS DISPLACEMENT = 0.00000    [ 0.00010]                            (YES)
MAX DISPLACEMENT = 0.00001    [ 0.00020]                            (YES)
RMS FORCE       = 0.00004    [ 0.00030]                            (YES)
MAX FORCE       = 0.00091    [ 0.00045]                            ( NO)
*******************************************************************************

先把cineb_energies里最后的能量作图，发现是这种趋势，于是感到很疑惑，因为插点结构是很合理的，通过metad也计算过自由能是有能垒的

于是我对5个结构的最后一帧进行单点计算
得到的能量趋势是：

和cineb输出的能量完全不同。
能量不准确的可能原因是什么呢？我的体系是有3个孤立的H原子，因此有三个单电子，会是这个原因导致neb的能量不准确吗？
其次我担忧，cineb输出的能量是可靠的吗，我们做cineb之后都要算单点去确认一遍吗，cp2k有类似高斯irc确定过渡态是否正确的方法吗
感谢各位老师！

Uus/pMeC6H4-/キ · 发表于 Post on 2026-7-28 18:32:44

本帖最后由 Uus/pMeC6H4-/キ于 2026-7-28 18:40 编辑

这是旧版CP2K的输出格式吧，2026.2版在NEB输出方面有不少改进，特别是直接有选项可以展示相对能量趋势（无需提取另外画图），欢迎尝试。本质上NEB并不一定需要DFT计算获得能量和受力，没有具体的DFT计算设置信息没法评价能量准不准。

第一个点和第二个点的几何结构有多少差别，拎出来重新单独在相同条件下做几何优化（不管之前到底有没有优化过）会有什么变化？

验证过渡态可以做dimer优化并做振动分析，观察是否有且仅有一个振动虚频、方向对应于反应的成键断键。

另外，参考下述两帖
http://bbs.keinsci.com/thread-59031-1-1.html
http://bbs.keinsci.com/thread-59102-1-1.html

dbxhtw · 发表于 Post on 2026-7-28 19:05:38

Uus/pMeC6H4-/キ发表于 2026-7-28 18:32
这是旧版CP2K的输出格式吧，2026.2版在NEB输出方面有不少改进，特别是直接有选项可以展示相对能量趋势（无 ...

谢谢老师的回复，我去学习一下新版cp2k的neb计算~~
老师，第二个点我去优化了一下，自发的回到了第一个点的结构，说明这个结构确实没有第一个点更稳定

sobereva · 发表于 Post on 2026-7-28 23:04:44

别乱找地方发帖，CP2K的帖子要发到第一性原理版块的CP2K分类里。这次给你移动了，以后注意

dbxhtw · 发表于 Post on 2026-8-1 16:34:21

sobereva 发表于 2026-7-28 23:04
别乱找地方发帖，CP2K的帖子要发到第一性原理版块的CP2K分类里。这次给你移动了，以后注意

好的老师，我之后注意

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

[CP2K] cp2k做cineb的能量需要再算单点确定吗