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

Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

We are NATO-free zone

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

Patrick SK Pang
pangsiukwong@gmail.com

Dear Pavel,

For NSTATE, should we be specify the value equal to or greater than the ISTATE value?

What is the state-tracking? What are the negative values used in the state tracking?

Regards,

Patrick

On Mon Dec 10 '12 8:25am, Solntsev Pasha wrote
----------------------------------------------
>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 Patrick.

>>>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.

>>>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

Mon Dec 10 '12 10:38am

This message read 832 times