“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

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

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

[复制链接 Copy URL]

65

帖子

0

威望

758

eV
积分
823

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: 16)

QQ图片20220531092951.png

本版积分规则 Credits rule

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

GMT+8, 2026-4-18 21:03 , Processed in 0.208373 second(s), 23 queries , Gzip On.

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