求助，核心空穴激发态的计算结果与文献相差很多

Samson.W

各位老师好，

我是一名分子物理专业的学生，是计算化学领域的新手，我在研究中需要计算N2O分子被X-ray单光子电离后形成的空穴离子激发态的能量，在光电离过程中一个电子被从1sigma轨道，也就是最内层的核心轨道，或者说是从氧原子的点位上被电离到了分子连续轨道上。因此我需要利用molpro计算最内层的核心空穴激发态（single core hole state）的能量。N2O分子可以被近似为C2v点群，有3个核心轨道，它们都是A1对称性的。它的轨道在molpro中的表示是：1.1，2.1，3.1，4.1，5.1，6.1，1.2，1.3，7.1，2.2，2.3， 我选择了4个额外的virtual orbitals : 3.2，3.3，8.1，9.1， 这15个分子轨道构成了我的完全活性空间，我使用"two-step"策略，先冻结核心轨道，优化价轨道，然后在冻结价轨道来优化核心轨道，并使用了state-averaging CASSCF方法。但是我最终得出的能量为-169.0014283 a.u.，而文献中的能量为-163.844022 a.u. 它们相差特别多。老师们能帮我看看我的molpro输入文件写的有没有问题吗。谢谢老师！！！

附件是输出文件，下面是我的输入文件：

***,N2O_12222222222_eq
memory,500,m
print,orbitals,civector,basis

Symmetry,x,y !use C2v symmetry
angstrom,
basis = aug-cc-pVQZ;
geometry={
N
N2, N, dist_NN
O, N2, dist_NO, N, 180.0
}
dist_NN=1.13164231
dist_NO=1.18840001

{hf;closed,0;wf,22,1,0;save,2100.2}
{multi;closed,0;occ,9,3,3,0;start,2100.2;wf,22,1,0;orbital,2140.2}

{multi;config;closed,0;occ,9,3,3,0;start,2140.2;freeze,1.1,2.1,3.1;wf,21,1,1;restrict,1,1,1.1;state,10;accuracy,econv=1.d-5;orbital,2141.2;}
{multi;config;closed,0;occ,9,3,3,0;start,2141.2;freeze,4.1,5.1,6.1,1.2,1.3,7.1,2.2,2.3,8.1,9.1;wf,21,1,1;restrict,1,1,1.1;state,10;accuracy,econv=1.d-5;orbital,2142.2;}
{multi;config;closed,0;occ,9,3,3,0;start,2142.2;freeze,1.1,2.1,3.1;wf,21,1,1;restrict,1,1,1.1;state,10;accuracy,econv=1.d-5;orbital,2143.2;}
{multi;config;closed,0;occ,9,3,3,0;start,2143.2;freeze,4.1,5.1,6.1,1.2,1.3,7.1,2.2,2.3,8.1,9.1;wf,21,1,1;restrict,1,1,1.1;state,10;accuracy,econv=1.d-5;orbital,2144.2;}
cas1=energy(1),cas2=energy(2),cas3=energy(3),cas4=energy(4),cas5=energy(5),cas6=energy(6),cas7=energy(7),cas8=energy(8),cas9=energy(9),cas10=energy(10);

{table,cas1;
save,N2O_12222222222.dat}
---