计算化学公社

标题: 采用CP2K的CINEB计算过渡态计算时完成一次离子步所需要SCF迭代次数的疑问。 [打印本页]

作者
Author: bblovelp 时间: 2022-5-31 09:46
标题: 采用CP2K的CINEB计算过渡态计算时完成一次离子步所需要SCF迭代次数的疑问。
本帖最后由 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

谢谢给位老师了。

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