Firefly and PC GAMESS-related discussion club


 
Learn how to ask questions correctly  
 
 
We are NATO-free zone
 



Re^3: Why is the first excited state equal to the third excited state? Mistakes in input files

Solntsev Pasha
solntsev@univ.kiev.ua


Dear Patrick.

from manual:

NSTATE =   Number of states to be found (excluding the
          ground state).

ISTATE =   State for which properties and/or gradient will
          be calculated.  Only one state can be chosen.

Optimization is a gradient calculation. If you are interesting in 1st state, set istate=1, if 2nd istate=2 etc.

I also want to draw your attention to symmetry issue.
ISTSYM - symmetry of states of interest. Default is zero, i.e., does not use
any symmetry during calculations. Setting this to the desired index of irrep
(according to PC GAMESS/Firefly numbering) will solve only for the states of
the desired symmetry and exploiting full (including non-abelian) symmetry of
molecule, thus significantly reducing computation time.

and finally you can activate state-tracking:
The state-tracking feature of the PC GAMESS/Firefly' CIS code can be
activated by selecting negative value of istate in the $cis group.
It is intended for geometry optimization of the excited states in the
case of root flipping.

Hope this help.
Best.

Pavel.


On Mon Dec 10 '12 4:28am, Patrick SK Pang wrote
-----------------------------------------------
>Dear Pavel,

>I want to make sure the meanings of the values for istate and nstate.
>For nstate, 6 is that the progam will give the energy results of the ground state, and the 1st, 2nd, 3rd, 4th, 5th and 6th excited states. Right?
>For istate, 1 is the first excited state, but not ground state. For my case, the program will give the optimized geometry of the first exctied state. Right?

>Regards,

>Patrick
>
>
>On Sun Dec 9 '12 7:35pm, Solntsev Pasha wrote
>---------------------------------------------
>>Dear atrick.

>>I both cases you requested only one state. It means, you don't have access to higher states.

>>For 1-T optimization i would recommend you to use this:
>>$CIS nstate=6 istate=1 mult=3 $end
>>for 3-T :
>>$CIS nstate=6 istate=3 mult=3 $end

>>Also, if you have the same symmetry for first and third triplet states, you can use ISTSYM to filter off all states undesired symmetry.

>>Check http://classic.chem.msu.su/gran/gamess/cis.htmlhttp://classic.chem.msu.su/gran/gamess/cis.htmlhttp://classic.chem.msu.su/gran/gamess/cis.htmlhttp://classic.chem.msu.su/gran/gamess/cis.htmlhttp://classic.chem.msu.su/gran/gamess/cis.htmlhttp://classic.chem.msu.su/gran/gamess/cis.htmlhttp://classic.chem.msu.su/gran/gamess/cis.htmlhttp://classic.chem.msu.su/gran/gamess/cis.html for more detailed information.

>>Good luck.

>>Pavel.

>>
>>
>>
>>On Sun Dec 9 '12 5:19pm, Patrick SK Pang wrote
>>----------------------------------------------
>>>Dear all,

>>>I think there are mistakes in my input files. It is because the optimized gemoetry and the energy of the first excited triplet state are equal to those of the third excited triplet state. Could you help me to point out the mistakes of my input files? My input files are as follows.

>>>First excited triplet state:
>>> $CONTRL RUNTYP=Optimize $END
>>> $CONTRL SCFTYP=RHF CITYP=CIS $END
>>> $CONTRL ICHARG=0 MULT=1 $END
>>> $CONTRL NZVAR=33 $END
>>> $CONTRL ICUT=20 ITOL=30 INTTYP=HONDO $END
>>> $SCF NCONV=10 FDIFF=.FALSE. $END
>>> $ZMAT DLC=.T. AUTO=.T. $END
>>> $BASIS EXTFIL=.T. GBASIS=acc-pvdz $END
>>> $CIS NSTATE=1 ISTATE=2 MULT=3 $END
>>> $SYSTEM TIMLIM=10000 MWORDS=100 $END
>>> $P2P P2P=.T. DLB=.T. $END
>>> $STATPT NSTEP=200 METHOD=GDIIS $END
>>> $DATA
>>>Molecule specification
>>>...
>>> $END

>>>Third excited triplet state:
>>> $CONTRL RUNTYP=Optimize $END
>>> $CONTRL SCFTYP=RHF CITYP=CIS $END
>>> $CONTRL ICHARG=0 MULT=1 $END
>>> $CONTRL NZVAR=33 $END
>>> $CONTRL ICUT=20 ITOL=30 INTTYP=HONDO $END
>>> $SCF NCONV=10 FDIFF=.FALSE. $END
>>> $ZMAT DLC=.T. AUTO=.T. $END
>>> $BASIS EXTFIL=.T. GBASIS=acc-pvdz $END
>>> $CIS NSTATE=1 ISTATE=4 MULT=3 $END
>>> $SYSTEM TIMLIM=10000 MWORDS=100 $END
>>> $P2P P2P=.T. DLB=.T. $END
>>> $STATPT NSTEP=200 METHOD=GDIIS $END
>>> $DATA
>>>Molecule specification
>>>...
>>> $END

>>>Thanks!

>>>Regards,

>>>Patrick


[ Previous ] [ Next ] [ Index ]           Mon Dec 10 '12 8:25am
[ Reply ] [ Edit ] [ Delete ]           This message read 1020 times