计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2761|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

65

帖子

0

威望

756

eV
积分
821

Level 4 (黑子)

本帖最后由 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

谢谢给位老师了。

QQ图片20220531092951.png (71.35 KB, 下载次数 Times of downloads: 9)

QQ图片20220531092951.png

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 15:11 , Processed in 0.196478 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list