Locating Ab Initio Transition States via Approximate Geodesics on Machine Learned Potential Energy Surfaces

Abstract

Efficient and reliable identification and optimization of transition state structures is a longstanding challenge in computational chemistry. Popular chain-of-states methods require hundreds if not thousands of ab initio calculations to generate initial guesses for local quasi-Newton optimizers, with persistent risk of collapse to an alternative stationary point on the potential energy surface (PES). Here, we show that high-quality guess structures for transition state optimization can be obtained by constructing the geodesic path between reactant and product structures on the PES generated by machine learning potentials (MLPs). We present an algorithm for optimization of such geodesic paths, as well as the associated codebase. We demonstrate effectiveness of this approach using the recent eSEN-sm-cons MLP. On average, the highest-energy point along these MLP geodesics requires 30% fewer quasi-Newton optimization steps to converge to the transition state compared to guesses from the fully ab initio frozen string method. Our approach therefore completely eliminates the need for ab initio calculations for generation of transition state guesses and considerably speeds up subsequent structural optimization. Geodesic construction on ML PES thus promises to be a useful approach for efficient computational elucidation of complex chemical reaction networks.

References in corpus (10)

• String Method for the Study of Rare Events
• Improved initial guess for minimum energy path calculations
• Exploration of Reaction Pathways and Chemical Transformation Networks
• Automated Discovery of Elementary Chemical Reaction Steps Using Freezing String and Berny Optimization Methods
• The Open Molecules 2025 (OMol25) Dataset, Evaluations, and Models
• Well-behaved versus ill-behaved density functionals for single bond dissociation: Separating success from disaster functional by functional for stretched H${}_2$
• Learning Smooth and Expressive Interatomic Potentials for Physical Property Prediction
• Jahn-Teller Distortion and Dissociation of CCl${}_4^{+}$ by Transient X-ray Spectroscopy Simultaneously at the Carbon K- and Chlorine L-Edge
• Efficient Transition State Searches by Freezing String Method with Graph Neural Network Potentials
• Incorporation of Internal Coordinates Interpolation into the Freezing String Method