So, from what I understand, there is no way to obtain the wave function of the system in the excited state T1, but I could only study it in T0.
You should set both first and second states with weight of 0.5. Then both the states should be able to solve.
]]>! def2-TZVP def2/J
! moread
%moinp "c2o_t.gbw"
%MaxCore 20000
%scf
TolG 1.0e-04
TolE 1.0e-06
maxiter 500
end
%casscf
nel 14
norb 10
mult 3
nroots 3
weights[1] = 1.0,0.0,0.0
end
%pal nprocs 10
end
*xyz 0 3
C 0.000000000 0.000000000 0.000000000
C 0.000000000 -1.369900105 0.000000000
O 0.000000000 1.162700089 0.000000000
I figured I should change the weight of the state, so I would have this: " weights[1] = 0.0,1.0,0.0 "
However, when I get the molden file and check the occupancies of the natural orbitals both files are very similar.
Thanks advance for any help provided.
R.
]]>