Firefly and PC GAMESS-related discussion club

Learn how to ask questions correctly

Re^4: Optimization of conical intersection with penalty function method: how to control the gap?
Dear Alex,

Thanks for your answer! I wonder how much should I increase the B parameter to converge to a conical intersection with a smaller gap (0.0001)? With the default value of B (3.5) the gap is ~0.003.

Best regards,
Evgeniy

On Sat Oct 22 '16 11:23pm, Alex Granovsky wrote
-----------------------------------------------
>Dear Evgeniy,

>SHIFT0 and SHIFT variables are primary used by Lagrange multiplier-
>based code.

>SHIFT can be used with penalty function based code but has very
>little impact on the energy gap when gap is large and shift
>is small. It is intended to trace rather large differences between
>states like 0.1 or 0.05 Hartree.

>Kind regards,
>Alex
>
>
>
>
>On Sat Oct 22 '16 7:39pm, GrEv wrote
>------------------------------------
>>Dear Alex,

>>Thanks very much for your reply.
>>Just a curios question, what does the SHIFT option do?

>>Best regards,
>>Evgeniy

>>On Sat Oct 22 '16 6:00pm, Alex Granovsky wrote
>>----------------------------------------------
>>>Dear Evgeniy,

>>>there is no direct control over the energy gap for penalty-based CoIn
>>>optimization procedures. The gap can be controlled indirectly using
>>>parameters defining a penalty function. The larger is the "strength"
>>>of penalty function the smaller will be the gap. With Firefly, the
>>>strength is controlled by the parameter B of \$mcaver group.

>>>Kind regards,
>>>Alex
>>>
>>>
>>>
>>>
>>>
>>>
>>>On Sat Oct 22 '16 11:50am, GrEv wrote
>>>-------------------------------------
>>>>Hello,

>>>>I am doing optimization of the minimum energy conical intersection using the penalty function method, with the following input (from the manual, p. 176) :

>>>>\$MCAVER TARGET=MIXED PENLTY=2 XGRAD=.F. \$END
>>>>\$STATPT OPTTOL=0.00005 NSTEP=500 method=qa trmax=0.05 dxmax=0.05
>>>> HSSEND=.f. \$END
>>>>\$NUMGRD ORDER=2 DELTA=0.01 ISTATE=1 JSTATE=2 NGRADS=2 PRAXES=1 \$END

>>>>Using this input I found a conical intersection (CoIn), with the final gap between the states of 0.003. I'd like however to have a smaller gap, let's say 0.0001. How this could be done. I've tried to use the SHIFT option in \$MCAVER, setting SHIFT=0.0001. It however didn't do the job. I got a CoIn with almost the same gap as before.
>>>>Many thanks!

>>>>Evgeniy

Sun Oct 23 '16 2:23pm
This message read 429 times