求助在CP2K上进行Mg原子的空位路径的CI-NEB计算难收敛

Ying_Zhang

之前发过帖子，在CP2K上跑不了NEB计算，每个点一直在Step0做SCF计算，始终都不进入Step1。最后按照相关帖子，增加插入NEB点数，可以顺利跑起来。但是，跑了300步，计算都没收敛，选取其中一个NEB点，Step 1是，如下：

*******************************************************************************
RMS DISPLACEMENT =   0.03060     [ 0.00010]                               ( NO)
MAX DISPLACEMENT =   0.49128     [ 0.00020]                               ( NO)
RMS FORCE        =   0.00592     [ 0.00030]                               ( NO)
MAX FORCE        =   0.12813     [ 0.00045]                               ( NO)
*******************************************************************************

Step 300时，如下：

*******************************************************************************
RMS DISPLACEMENT =   0.00129     [ 0.00010]                               ( NO)
MAX DISPLACEMENT =   0.02218     [ 0.00020]                               ( NO)
RMS FORCE        =   0.00153     [ 0.00030]                               ( NO)
MAX FORCE        =   0.02945     [ 0.00045]                               ( NO)
*******************************************************************************

体系就是简单的Mg和Bi的化合物，作Mg原子的空位路径的NEB计算，看过渡态能量和空位路径。下列是输入文件及输出文件，和某一个NEB点的Band文件。请问如何修改能快速收敛，且达到收敛限。如果目前CI-NEB计算的结果，没达到收敛限，是否可以那该结果进行dimer计算，因为发现dimer可以使用CG，

大概可以更好地收敛。感谢大家，帮忙解决！