Geometry optimization speedup through a geodesic approach to internal coordinates.

We present a new geodesic-based method for geometry optimization in a basis of redundant internal coordinates. Our method updates the molecular geometry by following the geodesic generated by a displacement vector on the internal coordinate manifold, which dramatically reduces the number of steps required to reach convergence. Our method can be implemented in any existing optimization code, requiring only implementation of derivatives of the Wilson B-matrix and the ability to solve an ordinary differential equation.