Firefly and PC GAMESS-related discussion club

 
Learn how to ask questions correctly  
 
 

 

Panwang Zhou
pwzhou@gmail.com

I want to calculate vertical excitation energy between the ground state and the first excited state with XMCQDPT2, first I do a state-averaged MCSCF on the optimized S0 geometries, with the following settings:

$MCSCF CISTEP=ALDET FULLNR=.F. SOSCF=.T. MAXIT=100 $END
$DET NCORE=72 NACT=12 NELS=12 NSTATE=6 WSTATE(1)=1,1 $END

From this discussion club and the GAMESS-US manual I understand why the NSTATE should be larger than 2 as determinants do not correspond to any spin number and the triplet will also be obtained when the CISTEP=ALDET. The WSTATE(1)=1,1 will exclude the triplet and only do averaging on the first and second singlet.

Then I using the optimized MCSCF MOs to do the XMCQDPT2 calculation with the following settings:

$MCSCF CISTEP=ALDET FULLNR=.F. SOSCF=.T. MAXIT=100 $END
$DET NCORE=72 NACT=12 NELS=12 NSTATE=6 WSTATE(1)=1,1 $END
$XMCQDPT NSTATE=10 EDSHFT=0.02 THRGEN=1D-12 INORB=2 $END
$XMCQDPT WSTATE(1)=1,1,-0 AVECOE(1)=1,1,-0 $END

Is this correct for calculating the vertical excitation energy between S0 and S1?

I set NSTATE=10 because sanya has said "The size of the effective Hamiltonian given by NSTATE has serious effect on the transition energy (the energy grows with NSTATE and levels off at NSTATE~10)."

However, I'm not clear the relation between the NSTATE, WSTATE and AVECOE in XMCQDPT2, could you comment on this? In the firefly manual I could not find the explain for WSTATE
in XMCQDPT.

Additional two small questions:
From this dicussion club I know there are two ways to avoid CASCI integral transformation step in XMCQDPT calculation.
(1) First do the MCSCF calculation with the following settings:
$mcscf iforb=.t. sd=.t. $end
$contrl wide=.t. $end

The do the XMCQDPT with:
$guess guess=moread norb=the_overall_number_of_mos $end
$xmcqdpt inorb=2 iforb=0 $end

My question is: when the iforb is true in $mcscf, can the obained MOs be used for next MCSCF calculation, where the geometry is close to the previous one.

(2) Do the XMCQDPT wih the following settings:
$mcscf iforb=.t. $end
$xmcqdpt iforb(1)=-1,1,1 $end

My question is:  what is the value should be set for inorb in this condition?


Best Regards,
Panwang Zhou