25.2.6 Keywords for redundant internal coordinates in $redund_inp

With the parameters in $redund_inp the generation of redundant internal coordinates can be modified. All entries have to be made in the control file before invoking the ired option. Important options are:

iprint n

print parameter for debug output: The larger n is, the more output is printed \(n\ge~0,n\leq~5\) (default: 0)

metric n

method for generating and processing of redundant internal coordinates
\(n\ge~-3,n\leq~3,n~\ne~0\) (default: 3)
Values for the metric option:

n = 1

 “Delocalized Coordinates”
The \(\mathbf{BmB^t}\) matrix is diagonalized for the complete set of redundant internal coordinates, matrix \(\mathbf{m}\) is a unit matrix.

n = -3

Delocalized Coordinates obtained with a modified matrix \(\mathbf{m}\), the values of \(\mathbf{m}\) can be defined by user input (see below).

n = -1

 “Hybrid Coordinates”
Natural internal coordinates are defined as in the old iaut option. If a cage remains, delocalized coordinates (as for n=1) are defined for the cage.

n = -2

Very similar to the \(n=1\) option, but for the remaining cage delocalized coordinates with modified matrix \(\mathbf{m}\) are defined as for \(n=-3\).

n = 2

 “Decoupled coordinates”
The redundant coordinates are divided into a sequence of blocks. These are expected to have decreasing average force constants, i.e. stretches, angle coordinates, torsions and “weak” coordinates. The \(\mathbf{BB^{t}}\) matrix is diagonalized for each block separately after the columns of \(\mathbf{B}\) were orthogonalized against the columns of \(\mathbf{B}\) of the the preceding blocks.

n = 3

 “Generalized natural coordinates”
Natural internal coordinates are defined first, for the remaining cage decoupled coordinates are defined.

type r

a positive real number, which is an approximate “force constant”, can be read in for each type of coordinate (see below). The force constants are used for the definition of the matrix \(\mathbf m\) in \(\mathbf{BmB^{t}}\).

Types of internal coordinates for the definition of m

The matrix \(\mathbf{m}\) is assumed to be a diagonal matrix. For each type of coordinate a different value for the force constants \(m_{ii}\) can be read in. Types of coordinates are:

stre

bond stretch (default: 0.5)

invr

inverse bond stretch (default: 0.5)

bend

bond angle (default: 0.2)

outp

Out of plane angle (default: 0.2)

tors

dihedral or “torsional” angle (default: 0.2)

linc

Special angle coordinate for collinear chains, bending of the chain a–b–c in the plane of b–c–d (default: 0.2)

linp

bending of the chain a–b–c perpendicular to the plane of b–c–d
(default: 0.2)

wstr

stretch of a “weak” bond, i.e. the bond is assumed to have a very low force constant, e.g. a “hydrogen bond” or a “van der Waals bond”
(default: 0.05)

winv

inverse stretch of a weak bond (default: 0.05)

wbnd

bond angle involving at least one weak bond (default: 0.02)

wout

Out of plane angle for weak bonds (default: 0.02)

wtor

dihedral angle for weak bonds (default: 0.02)

wlnc

linc coordinate for weak bonds (default: 0.02)

wlnp

linp coordinate for weak bonds (default: 0.02)