x 0 , and △x = x -x 0 . Most nonlinear optimization algorithms are based on this local quadratic approximation of the potential energy surface. By differentiating with respect to the coordinates, one could form an approximation for the gradient, which is given by: g(x) = g 0 + H 0 △x