H -1 k is the approximation to the inverse Hessian matrix at the k -th step. sk is the change in atomic positions. yk is the change in the gradient (force) at each step. The geometry optimization aims to move atoms in such a way that the total force on each atom approaches zero, i.e., it reaches an energy minimum.