|
|
本帖最后由 bblovelp 于 2022-5-31 09:51 编辑
各位老师,最近我在使用CP2K的CINEB算法计算过渡态时发现,在计算每个离子步所使用的SCF迭代轮数差异十分大,有的时候只要一轮SCF迭代就可以完成一次离子步计算,有的时候需要5-6轮SCF迭代,甚至会出现10轮以上的SCF迭代才计算完一个离子步(如图,其中红色方框为完成一次离子步,确认了下,每轮SCF迭代都收敛了),不知道造成这个原因是什么。有没有改进的空间。
CP2K计算采用的是OT算法,
&SCF
MAX_SCF 500
EPS_SCF 9.9999999999999995E-07
SCF_GUESS RESTART
&OT T
MINIMIZER DIIS
LINESEARCH 3PNT
PRECONDITIONER FULL_SINGLE_INVERSE
&END OT
&OUTER_SCF T
EPS_SCF 1.0000000000000001E-05
MAX_SCF 100
&END OUTER_SCF
CINEB设置如下:
&BAND
NPROC_REP 16
BAND_TYPE CI-NEB
NUMBER_OF_REPLICA 4
ROTATE_FRAMES F
ALIGN_FRAMES F
K_SPRING 8.0000000000000002E-02
&CONVERGENCE_CONTROL
MAX_DR 1.0000000000000000E-02
MAX_FORCE 1.0000000000000000E-03
RMS_DR 2.0000000000000000E-02
RMS_FORCE 1.0000000000000000E-03
&END CONVERGENCE_CONTROL
&CI_NEB
NSTEPS_IT 5
&END CI_NEB
&OPTIMIZE_BAND
OPT_TYPE DIIS
OPTIMIZE_END_POINTS
&END OPTIMIZE_BAND
&REPLICA
&COORD
~~~~
&END COORD
&END REPLICA
&PROGRAM_RUN_INFO SILENT
&END PROGRAM_RUN_INFO
&CONVERGENCE_INFO SILENT
&END CONVERGENCE_INFO
&END BAND
谢谢给位老师了。
|
|