## 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? **
GrEv

grevgen13@gmail.com

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 **317** times