and is computed by re-evaluating the energy of each configuration generated during a simulation using the perturbed force field parameters. Although more expensive than computing either the forward or backwards derivative, the central difference method should give a more accurate estimate of the gradient at the minima, maxima and transition points. By default a value of ℎ = × 10 -4 is used.