Firefly and PC GAMESS-related discussion club

 
Learn how to ask questions correctly  
 
 

 

Thomas Pijper
thomaspijper@hotmail.com

You could use a different method, such as UDFT or CASSCF. UMP2 is difficult -- as there is no way to calculate UMP2 gradients in 7.1.G, the only possibility is to use the TRUDGE method, which can only be used for small molecules and which isn't very straightforward in its use. TRUDGE is documented in the current manual; see pages 150 and 183-184. Keywords can be found in the older Firefly_input_rev002.pdf manual. In addition, some input examples using TRUDGE are available.


Kind regards,
Thom


On Mon Nov 25 '13 7:30am, Leira Ruth A. F. wrote
------------------------------------------------
>Dear Thom,

>Thanks for your reply.

>Does it mean that there's no way we can optimize geometry for radicals in Firefly 7.1.G?

>Sincerely,
>Leira

>On Thu Nov 21 '13 11:09am, Thomas Pijper wrote
>----------------------------------------------
>>Dear Leira,

>>Analytical gradients are available only for RHF wavefunctions. This is the case for both Firefly 7.1.G and Firefly 8.0.0. However, Firefly 8.0.0 is able to calculate gradients numerically, so you could use this to perform a geometry optimization with your system.

>>MP2 and numerical gradients are described fully in the manual.
>>
>>
>>Kind regards,
>>Thom
>>
>>
>>
>>On Thu Nov 21 '13 7:19am, Leira Ruth A. F. wrote
>>------------------------------------------------
>>>Dear colleagues,

>>>I made an mp2 optimization run with scftyp=uhf for a methyl radical. However, I get this message -  GRADIENT IS AVAILABLE ONLY FOR RHF MP2 RUNS then the calculation terminates. Does anyone know why this happens? I'm using Firefly 7.1.G in a desktop computer.

>>>Below is my input file and attached is the output file.

>>>#NAME: Leira Ruth A. Fulgueras
>>>#open shell optimization for CH3
>>> $CONTRL ICHARG=0 MULT=2 RUNTYP=Optimize SCFTYP=UHF
>>>    EXETYP=CHECK MPLEVL=2 COORD=cart MAXIT=100
>>>    DFTTYP=None ECP=None INTTYP=POPLE LOCAL=NONE
>>>    MOLPLT=.T. PLTORB=.T. AIMPAC=.F. NOSYM=0
>>> $END
>>> $SYSTEM TIMLIM=600 MEMORY=1000000 $END
>>> $BASIS  GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=0
>>>         NFFUNC=0 DIFFSP=.F. DIFFS=.F.
>>> $END
>>> $SCF DAMP=.false. DEM=.false. DIIS=.false. DIRSCF=.true.
>>>    EXTRAP=.true. RSTRCT=.false. SHIFT=.false. SOSCF=.true.
>>> $END
>>> $CUBE   CUBE=.T. MESH=COARSE
>>> $END
>>> $STATPT METHOD=GDIIS NSTEP=200 hssend=.t. nprt=-2 PURIFY=.t.
>>> $END
>>> $FORCE TEMP=298.15
>>> $END
>>> $GUESS GUESS=HUCKEL
>>> $END
>>> $DATA
>>>CH3
>>>C1
>>>C  6.000000 -0.000004 0.000004 -0.000125
>>>H  1.000000 -0.775777 -0.740656 0.000626
>>>H  1.000000 1.029320 -0.301520 0.000626
>>>H  1.000000 -0.253533 1.042166 0.000626
>>> $END
>>>