高手们好,请教一个casscf报错问题,屡次出现,不知如何是好
计算一个CrCl6_3-的Oh对称性离
molpro.xml
(42.6 KB, 下载次数 Times of downloads: 5)
子的casscf(实际计算中为D2h),优化和hf过程都没问题,但是在multi的时候出现了问题。
在hf计算中显示,三个单电子分别占据在三个简并轨道B1g/B2g/B3g中,LUMO为一组二重简并轨道Ag/Ag。
那么在一个(3,5)活性空间中,不改变自旋方向情况下共有10种态,分别为一个B1g*B2g*B3g=Ag,三个Ag*B1g*B2g=B3g,三个Ag*B1g*B3g=B2g,三个Ag*B2g*B3g=B1g
在molpro中Ag,B1g,B2g,B3g对应的不可约表示分别为1,4,6,7
但是仍然出错,找不到原因。
所以来请高手帮帮忙。
输入文件如下
***,crcl6
gprint,orbital=30
symmetry,x,y,z
angstrom
geometry=zmat
r(i)=2.4702767791
basis=cc-pvdz
{hf
wf,spin=3,charge=-3
save,2100.2
}
{multi
occ,38,22,22,9,11,9,9,0
closed,36,22,22,8,11,8,8,0
frozen,36,22,22,8,11,8,8,0
start,2100.2
wf,129,4,3
state,3
wf,129,6,3
state,4
wf,129,7,3
state,3
save,2140.2
}
e_cas(3,5)_trip=energy
put,molden,cas_35_triplet.molden
输出文件结尾
Micro-Iterations solved by the two-step solver
ITER. MIC NCI NEG ENERGY(VAR) ENERGY(PROJ) ENERGY CHANGE GRAD(0) GRAD(ORB) GRAD(CI) STEP TIME
MORE STATES THAN CONFIGURATIONS SPECIFIED: 0 3
谢谢大家
|