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

dbxhtw

求助各位老师：
分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确定过渡态是否正确的方法吗
感谢各位老师！