The BFGS algorithm is used for optimization in many computational chemistry methods, including geometry optimization in DMOL 3 . It is an iterative algorithm used to minimize a scalar function, typically the total energy of a system, with respect to atomic positions. The core of the algorithm involves updating an approximation to the inverse Hessian matrix H -1 (which contains second derivative information about the energy with respect to atomic positions).