Physical Chemistry

Sella, an open-source automation-friendly molecular saddle point optimizer

Authors

Abstract

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.

Version notes

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.

Content

Thumbnail image of main.pdf

Supplementary material

Thumbnail image of sella_si.tar.xz
Optimization benchmark input files
Contains the xyz files for all 500 initial guess structures for the saddle point benchmark in this work. Additionally, contains a Python script that can be used to perform the saddle point optimizations for these structures using Sella. Requires NWChem.