next up previous contents index
Next: Input Data Blocks Needed Up: Format of Keywords and Previous: Keywords for Module Ricc2   Contents   Index

Keywords for Module RELAX

$optimize options

define what kind of nonlinear parameters are to be optimized by relax and specify some control variables for parameter update.

Available options are:

internal on/off

optimize molecular structures in the space of internal coordinates using definitions of internal coordinates given in $intdef as described in Section 4.1 ( default: on).
redundant on/off

optimize molecular structures in redundant internal coordinates using definitions of redundant internal coordinates given in $redundant. For an optimization in redundant internal coordinates option internal has to be switched on too, and option cartesian has to be switched off (default: on).
cartesian on/off

optimize molecular structures in the space of (symmetry-distinct) cartesian coordinates (default: off).
basis on/off suboptions

optimize basis set exponents (default=off).

Available suboptions are:


exponents of uncontracted basis functions will be optimized after conversion into their logarithms (this improves the condition of the approximate force constant matrix obtained by variable metric methods and the behavior of the optimization procedure); scale factors of contracted basis functions will not be affected by the logarithm suboption

ALL basis set exponents will be optimized as scale factors (i.e. contracted blocks and single functions will be treated in the same way); if both suboptions (scale and logarithm) are given the logarithms of the scale factors will be optimized
global on/off

optimize a global scaling factor for all basis set exponents (default: off).

...\$interconversion}!on@{\normalfont\ttfamily on}}
\end{itemize} \end{minipage} }

$coordinateupdate options

define some variables controlling the update of coordinates.

Available options are:

dqmax real

maximum allowed total change for update of coordinates. The maximum change of individual coordinate will be limited to dqmax/2 and the collective change dq will be damped by dqmax/〈dq | dq〉 if 〈dq | dq〉 > dqmaxq
(default: 0.3)
interpolate on/off

calculate geometry update by inter/extrapolation of geometries of the last two cycles (the interpolate option is always switched on by default, but it is only active ANY time if steepest descent update has been chosen, i.e. $forceupdate method=none; otherwise it will only be activated if the DIIS update for the geometry is expected to fail)
statistics on/integer/off

provide a statistics output in each optimization cycle by displaying all (the last integer, default setting by define is 5) subsequent coordinates, gradient and energy values (default: on).
$gdiishistory file=char

the presence of this keyword forces relax to provide informational output about the usage of DIIS for the update of the molecular geometry.
$interconversion options default=off

special input related to the transformation of atomic coordinates between cartesian and internal coordinate spaces (default: off).

Available options are:


maximum number of iterations for the iterative conversion procedure internal → cartesian coordinates (default: 25).

convergence criterion for the coordinate conversion (default: 1.d-10).
on/off options

this switch activates special tasks: transform coordinates/gradients/ hessians between spaces of internal/cartesian coordinates using the definitions of internal coordinates given in $intdef:

available suboptions are:
cartesian -> internal coordinate gradient hessian
cartesian <- internal coordinate
the direction of the transformation is indicated by the direction of the arrow
\textbf{Note}: specification of \...
...optimize}\index{\$optimize@{\normalfont\ttfamily \$optimize}}!
\end{minipage} }
$forceupdate method options

this data group defines both the method for updating the approximate force constant matrix and some control variables needed for the force constant update.

Options for method:
no update (steepest descent)
ms suboptions
Murtagh-Sargent update
dfp suboptions
Davidon-Fletcher-Powell update
bfgs suboptions
Broyden-Fletcher-Goldfarb-Shanno update
dfp-bfgs suboptions
combined (bfgs+dfp) update
schlegel suboptions
Schlegel update
ahlrichs suboptions
Ahlrichs update (macro option)
suboptions if method=ms, dfp, bfgs, schlegel, ahlrichs
number of structures used
maximum number of geometries (= rank of the update procedure, for ahlrichs only)
minimum number of geometries needed to start update
if \emph{method}=\texttt{ms}, \text...
\texttt{maxgeo=2}, \texttt{mingeo=1} as default
\end{minipage} }
additional suboptions if method=ahlrichs
modus= char fmode
for an explanation see suboptions pulay given below e.g. ahlrichs numgeo=7 mingeo=3 maxgeo=4 modus=<g|dg> dynamic
if the macro option \texttt{...
...mensional BFGS (rank $n$) update for the hessian
\end{itemize} \end{minipage} }
pulay suboptions

try to find an optimal linear combination of the coordinates of the numpul previous optimization cycles as specified by modus (see below).

Available suboptions are:

number of geometries to be utilized

maximum number of geometries

minimum number of geometries needed to start update
modus=char fmode

char=SPMlt;g|g>; or SPMlt;g|dq>; or SPMlt;dq|dq>; defines the quantity to be minimized (dq = internal coordinate change).
fmode specifies the force constants to be used (only if char=SPMlt;g|dq>; or SPMlt;dq|dq>;!)
fmode=static: use static force constants
fmode=dynamic: use updated force constants

real defines the threshold for the quantity g * dq/| g| * | dq| which defines the angle between gradient vector and coordinate change (default: 0.1). If pulay is used in connection with a multidimensional BFGS update for the hessian than the default is real=0.0. If $ {\frac{{g\cdot\mathrm{d}q}}{{\vert g\vert\ast\vert\mathrm{d}q\vert}}}$ > -real the pulay update for the geometry is expected to fail and will be ignored. For example:
pulay numpul=4 maxpul=4 minpul=3 modus=<dq|dq> static fail=0.2
options for $forceupdate

update only the diagonal force constants (update for off-diagonals will be suppressed) (only active if method=ms, dfp, bfgs)
offdamp real

this allows to damp off-diagonal force constants by 1/ real (compare offreset, which discards off-diagonals completely). Only values > 1.0 will be accepted. This option is active only within one relax run and will be disabled automatically by relax. This is useful in difficult cases, where the non-diagonal update has lead to too large non-diagonal elements of the hessian.

reset off-diagonal force constants to zero. This option will be active for the current optimization cycle only, i.e. it will be removed by relax after having discarded off-diagonals!

optimization cycle specification of a maximum energy change allowed (given in mHartree) which will be accepted using the actual approximate force constant matrix from $forceapprox; if this energy change will be exceeded, the force constants will be scaled appropriately
(The default: 0.0 means NO action)

scaling factor for the input hessian (default: 1.0).

lower bound for eigenvalues of the approximate hessian (default: 0.005); if any eigenvalue drops below threig, it will be shifted to a reasonable value defined by:
reseig= real
default: texttt0.005.

upper bound for eigenvalues of the hessian; if any eigenvalue exceeds thrbig, it will limited to this value (default: 1000.0).

damp the variable metric update for the hessian by 1/(1 + real) (default: 0.0).
$forceinit option

specify initialization of the (approximate) force constant matrix.

Available options are:


this activates or deactivates initialization; if on has been set, relax will provide an initial force constant matrix as specified by one of the possible initialization options as described below and will store this matrix in data group $forceapprox; after initialization relax resets $forceinit to off!

provide a diagonal force constant matrix with:

available suboptions are:


this will lead to an assignment of diagonal elements (default: 1.0)).

this will lead to an assignment of initial force constant diagonals depending on the coordinate type.

Provide individual defined force constant diagonals for
  • internal coordinates (supplied in $intdef ... fdiag=..)
  • a global scale factor ( $global ... fdiag=..)
This does not work for basis set optimization. For the correct syntax of `fdiag=..' see descriptions of $intdef, $global

read a cartesian (e.g. analytical) hessian from $hessian and use it as a start force constant matrix; if $optimize internal has been set: use its transform in internal coordinate space. If large molecules are to be optimized, it may be necessary (large core memory requirements!) to deactivate the numerical evaluation of the derivative of the B-matrix with respect to cartesian coordinates, which is needed to transform H(cart)→H(int) exactly by specifying no dbdx.
$last SCF energy change = real
$last MP2 energy change = real

These keywords depend on the optimization task to be processed and are updated by the corresponding program (i.g. SCF energy).
$m-matrix options

This data block contains non-default specifications for the m-matrix diagonals. This is of use if some cartesian atomic coordinates shall be kept fixed during optimization.

Available options are:

integer real real real

atomic index followed by diagonal elements of the m-matrix for this atom
$scratch files

The scratch file ftmp allocated by relax can be placed anywhere in your file systems instead of the working directory by referencing its pathname in this data group as follows:

       $scratch files
        relax   ftmp          path/file

The first column specifies the program, the second column the scratch file and the third column the pathname of the file to be used as scratch file.

next up previous contents index
Next: Input Data Blocks Needed Up: Format of Keywords and Previous: Keywords for Module Ricc2   Contents   Index