On Wed Oct 10 '12 10:12pm, Alex Granovsky wrote
>sorry for large delay on my side.
>>Q1. ok, it sounds very interesting. So, if an order of states from
>CAS-CI part is 1,1,1,0,1 we can setup kstate(1)=1,1,1,0,1 and
>basically remove "undesired" state from PT2 part.
>> Does it make sense?
>There are some situations when this for sure makes (some) sense.
>For instance, assume one is dealing with diatomic molecule.
>The full symmetry group is either Cinfv or Dinfh
>which is a non-abelian symmetry group. The effective Hamiltonian
>is thus block-diagonal and (provided one is interested in states of
>some particular symmetry type) one can remove other states from PT treatment.
>> Is this equivalent to setting of the AVECOE and WSTATE (XMCQDPT)
>arrays to be AVECOE(1)=1,1,1,0,1 and WSTATE(1)=1,1,1,0,1. If not, what
>is the best strategy, to use kstate or avecoe/wstate? Or, maybe, we
>can just use nstate/wstate/avecoe combination and forget about kstate.
>No you in general cannot. Both avecoe and wstate applies to renumbered
>CI states i.e. to state numbers generated after renumbering procedure.
>In general, kstate and avecoe/wstate are more or less independent.
>As to NSTATE/kstate, kstate provides the finer control.
>>Q2. What if those additional eigenvectors contribute to some low
>lying states? Do we need repeat CASSCF part once again and include
>those extra states?
>The answer depends on the magnitude of contribution. If it is large
>the answer can be "yes, it would be better to do that". As far
>as I remember this self-consistent procedure of the selection of
>the states in state-averaging was once discussed a year ago or so
>on the forum, you may find my older comments on this.
>>>>Q3: Is it good idea to do SA-CASSCF/XMCQDPT via singlets and triplets
>>>simultaneously in case of high spin-orbit coupling. Actually, i am
>>>going to check singlets only and singlets-triplets, but maybe it make
>>>I'd suggest you to use to use SA-CASSCF averaged over both singlet
>>>and triplet states to generate the same set of MOs to be used in PT2
>>>computations. Firefly v. 8.0.0 has some new features specifically
>>>tailored for these types of jobs which I'll describe in my next post
>>>to this thread.
>The couple of advises are:
>1. Use aldet CASCI code with PURES=.f. to generate orbitals with SA-CASSCF.
$mcscf iforb=.t. $end $xmcqdpt iforb(1)=-1,1,1 $end
>3. Perform two sets of QDPT calculations, one for singlets:
$xmcqdpt mult=1 $end
>and second for triplets:
$xmcqdpt mult=3 $end
>All the best,