计算化学公社

标题: casscf GIC扫描时能量曲线出现跳动 [打印本页]

作者
Author: 隔壁老张sir 时间: 2023-4-16 16:52
标题: casscf GIC扫描时能量曲线出现跳动
输入文件：

#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前与文献吻合，但rb=1.7~2.6时与文献结果有差异，当rb=2.7时能量突然跳动（增大），竟然就跟文献的结果对上了。想请教各位老师，这个能量跳动是什么原因？（是因为扫描步长太大了吗？我尝试改小步长但是发散了）跳动前的能量是对的还是之后的是对的？（跳动之前和文献对不上，跳动之后和文献能对上）

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