Efficient Implementation of Gaussian Process Regression Accelerated Saddle Point Searches with Application to Molecular Reactions

Abstract

The task of locating first order saddle points on high-dimensional surfaces describing the variation of energy as a function of atomic coordinates is an essential step for identifying the mechanism and estimating the rate of thermally activated events within the harmonic approximation of transition state theory. When combined directly with electronic structure calculations, the number of energy and atomic force evaluations needed for convergence is a primary issue. Here, we describe an efficient implementation of Gaussian process regression (GPR) acceleration of the minimum mode following method where a dimer is used to estimate the lowest eigenmode of the Hessian. A surrogate energy surface is constructed and updated after each electronic structure calculation. The method is applied to a test set of 500 molecular reactions previously generated by Hermez and coworkers [J. Chem. Theory Comput. 18, 6974 (2022)]. An order of magnitude reduction in the number of electronic structure calculations needed to reach the saddle point configurations is obtained by using the GPR compared to the dimer method. Despite the wide range in stiffness of the molecular degrees of freedom, the calculations are carried out using Cartesian coordinates and are found to require similar number of electronic structure calculations as an elaborate internal coordinate method implemented in the Sella software package. The present implementation of the GPR surrogate model in C++ is efficient enough for the wall time of the saddle point searches to be reduced in 3 out of 4 cases even though the calculations are carried out at a low Hartree-Fock level.

14 pages, 4 figures

References in corpus (9)

• NWChem: Past, Present, and Future
• Materials Cloud, a platform for open computational science
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• Nudged elastic band calculations accelerated with Gaussian process regression
• Gaussian Process Regression for Transition State Search
• Global optimization of atomistic structure enhanced by machine learning
• Minimum energy path calculations with Gaussian process regression
• IRA: A shape matching approach for recognition and comparison of generic atomic patterns
• Efficient Implementation of Gaussian Process Regression Accelerated Saddle Point Searches with Application to Molecular Reactions

Cited by in corpus (5)

• Metatensor and metatomic: foundational libraries for interoperable atomistic machine learning
• Efficient Implementation of Gaussian Process Regression Accelerated Saddle Point Searches with Application to Molecular Reactions
• Bayesian Hierarchical Models for Quantitative Estimates for Performance metrics applied to Saddle Search Algorithms
• Two-dimensional RMSD projections for reaction path visualization and validation
• Reproducible Orchestration of Best Practices for Reaction Path Optimization with the Nudged Elastic Band