Please, keep in mind that CASSCF is NOT a black-box method. You should know in advance what are you going to see. In your example, HN3 particle (zero charge, singlet state) dissociates into N2 (zero charge, singlet state?) and HN (zero charge, singlet state, too?). Is this mechanism really possible? Does hydrogen azide decompose through a monomolecular mechanism? Are the first-stage products really N2 and HN? If so, are you confident in the charge and/or multiplicity of the resulting particles? Probably, you should try ALDET instead of GUGA to include states of different multiplicity.
Anyway, you should start at the very beginning, i.e., with a hydrogen azide molecule, whose ground state is well described by a single close-shell reference (or by some relatively simple and understandable combination of configurations), and whose lowest excited states mostly consist of single excitations from the ground state. Even if you want to average over 3 lowest states (S0, S1, and S2), set NSTATE=5 or even more and add NTRACK=3 to $MCSCF group. This will help you detect root flipping if it occurs.
On Thu Aug 16 '12 6:05pm, Serzhan wrote
>On Thu Aug 16 '12 2:13pm, sanya wrote
>>The problems in QDPT frequently arise from errors or problems in CASSCF. From your outputs, I see that root flipping occurs: in 136CASSCF-10-10PT34 State 1 was closed-shell (all occupancies are close to 2 and 0), while in 135CASSCF-10-10PT34 I see fractional occupancies of orbital 5 (1.43) and 6 (0.57). Can you show us your CASSCF inputs and outputs?
>>In addition, you'd better set $XMCQDPT group as follows:
>> $XMCQDPT INORB=1
>>! (instead of kstate(1))
>>! (additional parameters: you need 3 states, so 3 1's)
>> wstate(1)=1,1,1,-0 avecoe(1)=1,1,1,-0
>>I think this improvement alone will not help you, unless you find out what is wrong with your CASSCF.
>Thank you for your reply.
>I attached the CASSCF inputs and outputs with this post. It seems the MOs in active space remain intact, while the CI roots differ at points of interest. So that can be the cause of the jump. However, the root flipping, you've mentioned, is due to the physics of the problem. So It will be great if the energy profile could be kept smooth, while root flipping is present.
>Thanks for your help,