Up: IBRION-tag, NFREE-tag
N.B. This document is no longer maintained, please visit our wiki.
A conjugate-gradient algorithm
(a simple discussion of this algorithm can be found for instance in )
is used to relax the ions into their instantaneous
In the first step ions (and cell shape) are changed along the direction of the
steepest descent (i.e. the direction of the calculated forces and stress tensor).
The conjugate gradient method requires a line minimization, which
is performed in several steps:
i) first a trial step into the search direction (scaled gradients) is done,
with the length of the trial step controlled
by the POTIM parameter (section 6.23). Then the energy and
the forces are recalculated.
The approximate minimum of the total energy is calculated from a cubic
(or quadratic) interpolation
taking into account the change of the total energy and the change of the forces
(3 pieces of information), then a
corrector step to the approximate minimum is performed.
iii) After the corrector step the forces and energy are recalculated and
it is checked whether the forces
contain a significant component parallel to the previous search
direction. If this is the
case, the line minimization is improved by further corrector steps
using a variant of Brent's algorithm.
In the first ionic step the forces are calculated for the initial configuration
read from POSCAR , the second step is a trial (or predictor step),
the third step is a corrector step. If the line minimization is sufficiently
accurate in this step,
the next trial step is performed.
||corrector step, i.e. positions corresponding to anticipated minimum
N.B. Requests for support are to be addressed to: email@example.com