计算化学公社

标题: 大体系大活性空间CASSCF如何加快收敛 [打印本页]

作者
Author: BinWang 时间: 2026-2-2 10:19
标题: 大体系大活性空间CASSCF如何加快收敛
各位老师：
   大家好！我用orca计算一个84原子的体系，波函数是PySCF在DRMG下收敛好的，现在已经跑了快两天了，还没有收敛。目前还没有收敛的迹象，这是目前输出的能量：
Maximum residual norm  :    0.000475114390

                  ****Iteration 5****
               E    : -20675.029533545010 Eh
Maximum residual norm  :    0.000096184993

                  ****Iteration 6****
               E    : -20675.029572651649 Eh
Maximum residual norm  :    0.000020258312

                  ****Iteration 7****
               E    : -20675.029581228744 Eh
Maximum residual norm  :    0.000004733406

                  ****Iteration 8****
               E    : -20675.029583302465 Eh
Maximum residual norm  :    0.000001203776

                  ****Iteration 9****
               E    : -20675.029583808267 Eh
Maximum residual norm  :    0.000000294018

                  ****Iteration 10****
               E    : -20675.029583932464 Eh
Maximum residual norm  :    0.000000071615

                  ****Iteration 11****
               E    : -20675.029583963817 Eh
Maximum residual norm  :    0.000000018196

                  ****Iteration 12****
               E    : -20675.029583971329 Eh
Maximum residual norm  :    0.000000004319

                  ****Iteration 13****
               E    : -20675.029583973035 Eh

这是我的输入文件：
! X2C TightSCF moread
%pal
  nprocs 32
  end
%maxcore 10000
%moinp "Diradical_CASSCF_CMO_o.gbw"

%casscf
nel 30
norb 28
mult 3
nroots 1
maxiter 800
CIStep ICE
ci
   maxiter 400
   maxdim 20
end
printlevel 4
end
%coords
  Units = angs
  Charge = 1
  Mult = 3
  Coords
  Ir       5.365000022    8.211162982    10.130604986
我想要在这个CASSCF收敛以后，做态平均的CASSCF计算，目前的速度实在是太慢了，各位老师有什么办法加速收敛？另外Maxdim这个值该怎么设置，我是不是设置太小了？

作者
Author: wzkchem5 时间: 2026-2-2 12:07
这个速度是正常的。如果已经做过DMRG计算了，可以看看是否有轨道的占据数特别接近2或者0，去掉这些轨道再试试

作者
Author: BinWang 时间: 2026-2-2 12:55

wzkchem5 发表于 2026-2-2 12:07
这个速度是正常的。如果已经做过DMRG计算了，可以看看是否有轨道的占据数特别接近2或者0，去掉这些轨道再试 ...

非常感谢王老师

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