求助大佬们有关NiCCO负离子的VDE计算的问题

QJTX

我要算NiCCO负离子的VDEs，首先用高斯先在B3LYP下确认二重态负离子和三重态中性的基态结构，然后计算得到他的VDE值。然后用CCSD(T1)诊断发现它的多组态特征明显。我就拿着高斯得到的负离子结构去ORCA算三重态中性激发态能量，用了casscf caspt2和nevpt2都算了，结果和我的光电子能谱图实验都对不上。为了确认它的活性空间，我又用了MOKIT，但是我用负离子的结构，电荷和多重度用的中性，MOKIT算就不收敛了。但是用负离子结构和负离子多重度mokit可以得到的活性空间，但是这个活性空间和我预期判断的差很多。我不知道我计算VDEs的流程有没有问题，这个已经困扰了我很久了。求大佬帮助。这是我的负离子结构。有没有大佬能帮我看看我活性空间选择或者计算流程和方法上哪里有问题
%mem=20GB
%nprocshared=5
#p CASSCF/def2tzvpp
mokit{}
-1 2
Ni   -0.971214000000     -0.719264000000      0.000000000000
C     1.277902000000      0.922282000000      0.000000000000
C     0.000000000000      0.774109000000      0.000000000000
O     2.440824000000      1.245131000000      0.000000000000

我用ORCA算激发态能量的时候是这样写的，这个活性空间是我自己判断的，也不知道对不对
!DEF2-TZVPP def2/J def2-tzvpp/C Conv NoPop VeryTightSCF
%maxcore 100000   
%pal nprocs 5
end

%casscf
        nel      18               
         norb     13               
        Maxiter    300
        trafostep ri
        nroots 7,7
        mult 3,1
        PTMethod sc_nevpt2
        PTSettings  
                D4Step efficient
        QDType QD_VanVleck
        end
end

*xyz 0 3                           
Ni   -0.971214000000     -0.719264000000      0.000000000000
C     1.277902000000      0.922282000000      0.000000000000
C     0.000000000000      0.774109000000      0.000000000000
O     2.440824000000      1.245131000000      0.000000000000
*

QJTX

求救