计算化学公社

标题: 求助在CP2K上进行Mg原子的空位路径的CI-NEB计算难收敛 [打印本页]

作者
Author: Ying_Zhang 时间: yesterday 18:59
标题: 求助在CP2K上进行Mg原子的空位路径的CI-NEB计算难收敛
之前发过帖子，在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，

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

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