求助在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，

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

kaikai

1，计算尝试开OT
2，你后面使用的计算核数是副本数*进程数达到了384，这个你是使用的超算吗？有那麽多的核数吗？

Ying_Zhang

kaikai 发表于 2026-3-6 09:48
1，计算尝试开OT
2，你后面使用的计算核数是副本数*进程数达到了384，这个你是使用的超算吗？有那麽多的核 ...

哦哦 好的 谢谢老师建议 我会尝试一下OT 是的 是我们所里的集群 最大额度 一个任务可以用20个字节（48核）

Ying_Zhang

kaikai 发表于 2026-3-6 09:48
1，计算尝试开OT
2，你后面使用的计算核数是副本数*进程数达到了384，这个你是使用的超算吗？有那麽多的核 ...

那老师建议拿现在跑了300步的CI-NEB结构 再跑Dimer吗