We present a new algorithm for the optimization of molecular structures to saddle points on the potential energy surface using a redundant internal coordinate system. This algorithm automates the procedure of defining the internal coordinate system, including the handling of linear bending angles, e.g. through the addition of dummy atoms. Additionally, the algorithm supports constrained optimization using the null-space sequential quadratic programming formalism. Our algorithm determines the direction of the reaction coordinate through iterative diagonalization of the Hessian matrix, and does not require evaluation of the full Hessian matrix. Geometry optimization steps are chosen using the restricted step partitioned rational function optimization method, and displacements are realized using a high-performance geodesic stepping algorithm. This results in a robust and efficient optimization algorithm suitable for use in automated frameworks. We have implemented our algorithm in Sella, an open source software package designed to optimize atomic systems to saddle point structures. We also introduce a new benchmark test comprising 500 molecular structures that approximate saddle point geometries and show that our saddle point optimization algorithm outperforms the algorithms implemented in several leading electronic structure theory packages.
New benchmark comparison of Sella vs Q-Chem, Psi4, and NWChem. Rearranged some sections (mainly, 2.1 becomes Appendix A and 2.1.1 is promoted to section 2.1). Clarified some math in section 2.2.
Optimization benchmark input files