19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

42 GEOMETRY OPTIMIZATION (OPTG) 317<br />

42.2.16 Reaction path following options (OPTION)<br />

OPTION,key=param;<br />

where key can be<br />

IDIR<br />

STPTOL<br />

SLMAX<br />

If starting at a transition state (or near a transition state) determine<br />

where to take the first step. If IDIR=0 is chosen, the first step will be<br />

towards the transition state. This is the default. If IDIR=1 is given<br />

in the input the first optimization step will be along the ”transition<br />

vector” i.e. the hessian eigenvector to the smallest eigenvalue which<br />

points down towards the minimum. If using a larger IDIR parameter,<br />

the first step will be larger; if using a negative value, the first step will<br />

be in the opposite direction.<br />

If using an updated hessian matrix, this parameter determines what<br />

update to take. If the step size between two subsequent points on<br />

which the steepest decent lines are puzzled together is smaller than<br />

stptol (i.e. if we are close to a minimum) the BFGS update is used,<br />

otherwise it is Powell update. The default value of stptol is 1.d − 6.<br />

This option is only valid with the old version of the reaction path following<br />

algorithm (i.e. METHOD,SRSTEEP). In this algorithm slmax<br />

determines the frequency of the recalculation of the numerical hessian.<br />

If the total step size of the last steps exceeds slmax the hessian<br />

will be recalculated, otherwise it will be updated. By default slmax is<br />

two times the maximum step size of the optimization step steplength<br />

(see STEP section 42.2.12). If you are using METHOD,QSD, the<br />

SLMAX option is obsolete and the NUMHES command (see above)<br />

should be used instead.<br />

42.2.17 Optimizing energy variables (VARIABLE)<br />

VARIABLE,name;<br />

Defines a variable name which holds the energy value to be optimized in using finite differences.<br />

By default, this is ENERGY(1) as set by the most recent program. Other variables which can<br />

be used are<br />

ENERGY(i) holds last energy for state i.<br />

ENERGR(i) holds last reference energy for state i.<br />

ENERGD(i) holds last Davidson corrected energy for state i.<br />

ENERGP(i) holds last Pople corrected energy for state i.<br />

ENERGC<br />

ENERGT(1)<br />

ENERGT(2)<br />

ENERGT(3)<br />

holds CCSD (QCI, BCCD) energy in CCSD(T) [QCI(T), BCCD(T)]<br />

calculations (single state optimization).<br />

holds CCSD(T) energy in CCSD(T) calculations (single state)<br />

holds CCSD[T] energy in CCSD(T) calculations (single state).<br />

holds CCSD-T energy in CCSD(T) calculations (single state).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!