casscf GIC扫描时能量曲线出现跳动

隔壁老张sir

输入文件：

#p opt=addgic casscf(12,12) maug-cc-pvtz guess(mix,always) SCF(novaracc,noincfock) int=acc2e=12 nosymm SCF=maxcyc=500

0 1
N              -1     0.47544795    0.00000000   -0.27450000
N              -1    -0.47544795    0.00000000    0.27450000
N               0     0.00000000    0.00000000  -10.40000000
N               0     0.00000000    0.00000000   -9.60000000

XC1=XCntr(1,2)
YC1=YCntr(1,2)
ZC1=ZCntr(1,2)
XC2=XCntr(3,4)
YC2=YCntr(3,4)
ZC2=ZCntr(3,4)
F1F2(Frozen)=[(XC1-XC2)^2+(YC1-YC2)^2+(ZC1-ZC2)^2]
D34(StepSize=0.05,NSteps=40)=R(3,4)



rb=1.7~2.6时与文献结果有差异，当rb=2.7是能量突然跳动（增大），竟然就跟文献的结果对上了。这个能量跳动是什么原因？（是因为扫描步长太大了？）跳动前的能量是对的还是之后的是对的？（跳动之前和文献对不上，跳动之后和文献能对上）