Next: Definition of Internal Coordinates
Up: Program Relax
Previous: Optimization of General Coordinates
Contents
Index
In a Newtontype geometry update procedure often only a crude
approximation to the force constant matrix H^{k}
is available. What
can be done then is to update
F^{k} = (H^{k})^{1}
in each iteration using
information about previous coordinates and gradients. This constitutes
the quasiNewton or variable metric methods of which there are a few
variants:
 Murtagh/Sargent (MS):
F^{k} = F^{k1} + 

 Broyden/Fletcher/Goldfarb/Shanno (BFGS):
F^{k} = F^{k1} + 

 Davidon/Fletcher/Powell (DFP):
F^{k} = F^{k1} +  

 combined method (BFGS/DFP): If
S1 < (S  1)S1
and S1 > 0
perform DFP update, otherwise BFGS.
The meaning of the symbols above is as follows:

F^{k} = (H^{k})^{1}
 approximate inverse force constant
matrix in the kth iteration.s
 q^{k}
 general coordinates in the kth iteration.
 G^{k}
 gradients in the kth iteration.

dq^{k1} = q^{k}  q^{k1}

dg^{k1} = g^{k}  g^{k1}

Z^{k1} = dq^{k1}  F^{k1}dG^{k1}

S1 = (dq^{k1})^{†}dg^{k1}

S = 1 + ((dg^{k1})^{†}F^{k1}dG^{k1})/(S1)

An alternative is to use update algorithms for the hessian H^{k}
itself:
Ehrig, Ahlrichs :
Diagonal update for the hessian by means of a least squares fit
H^{k}_{ii} = 

with the new estimate h
for the diagonal elements obtained by
h_{i} = 

and the error d
obtained by the regression
d_{i} = . 

Another alternative is to use DIISlike methods: structure
optimization by direct inversion in the iterative subspace. (See ref.
[33] for the description of the algorithm). The DIIS
procedure can often be applied with good success, using static or
updated force constant matrices.
Any of the algorithms mentioned above may be chosen. Recommended is
the macro option ahlrichs
, which leads to the following actions
(n
is the maximum number of structures to be included for the
update, default is n = 4
):

ncycles < n
:
 geometry update by
inter/extrapolation using the last 2 geometries.

ncycles≥n
:
 diagonal update for the
hessian as described above; DIISlike update for the geometry.

 G < thr
:
 BFGStype update of the
hessian and quasiNewton update of (generalized) coordinates.
References for the algorithms mentioned above:
[34,30,35,33,36,37]
Next: Definition of Internal Coordinates
Up: Program Relax
Previous: Optimization of General Coordinates
Contents
Index
TURBOMOLE M