Just in case, you can get (print out) all the states involved in the averaging at each step of optimization run by using the option NPFLG(9)=1 in $MCSCF.
On Tue Apr 11 '17 3:38pm, Bernhard Dick wrote
>I am getting confused with how Firefly renumbers states. I am searching for a conical intersection between the states S0 and T1 of a system with an even number of electrons. I encounter problems with convergence of MCSCF unless I include also the next two states in the MCSCF. Therefore, i am using ALDET with
>NSTATE=4 WSTATE(1)=1,1,1,1 PURES=.F.
>to get state-averaged MCSCF. The states are calculated in the order S0, T1, S1, T2
>I start the search for conical intersection with these optimized orbitals, and set ISTATE=1 in $MCSCF and JSTATE=2 in $MCAVER. This should give me the crossing between S0/T1. Then I stop the calculation after NSTEP=100 (convergence is not obtained, but the energies do not change much, and the gap is 2-3 mHartree). Since Firefly at this point does not tell me the engergies of all 4 states, I do an extra MCSCF with the last geometry and orbitals, and it turns out that the crossing obtained is between S1 and T2, i.e. between states 3 and 4!
>I used NTRACK=4 in $MCSCF and the group
> $TRACK TOL=1.4 FREEZE=.T. STICKY=.F. UPDATE=.T. RESET=.F. DELCIV=.F. $END
>according to the manual.
>Finally, I want to calculate MCQDPT2 corrections at the final geometry. But when I add the group
>$MCQDPT EDSHFT=0.02 ISTSYM=1 KSTATE(1)=1,1,1,1 $END
>I do not get the corrections for the 4 states calculated with ALDET, but for the lowest 4 singlet states (although MULT was not set in $CONTRL). How can I do the MP2-correction for the same states that are calculated in the MCSCF?
>with best regards,