Affiliation:
1. HQS Quantum Simulations GmbH, Rintheimer Straße 23, 76131 Karlsruhe, Germany
2. Theoretical Physics, Saarland University, 66123 Saarbrücken, Germany
Abstract
In this work, we propose a multi-level protocol for routine theoretical studies of chemical reaction mechanisms. The initial reaction paths of our investigated systems are sampled using the Nudged Elastic Band (NEB) method driven by a cheap electronic structure method. Forces recalculated at the more accurate electronic structure theory for a set of points on the path are fitted with a machine learning technique (in our case symmetric gradient domain machine learning or sGDML) to produce a semi-local reactive potential energy surface (PES), embracing reactants, products and transition state (TS) regions. This approach has been successfully applied to a unimolecular (Bergman cyclization of enediyne) and a bimolecular (SN2 substitution) reaction. In particular, we demonstrate that with only 50 to 150 energy-force evaluations with the accurate reference methods (here complete-active-space self-consistent field, CASSCF, and coupled-cluster singles and doubles, CCSD) it is possible to construct a semi-local PES giving qualitative agreement for stationary-point geometries, intrinsic reaction coordinates and barriers. Furthermore, we find a qualitative agreement in vibrational frequencies and reaction rate coefficients. The key aspect of the method’s performance is its multi-level nature, which not only saves computational effort but also allows extracting meaningful information along the reaction path, characterized by zero gradients in all but one direction. Agnostic to the nature of the TS and computationally economic, the protocol can be readily automated and routinely used for mechanistic reaction studies.
Funder
German Federal Ministry of Economic Affairs and Climate Action
Reference54 articles.
1. Dykstra, C.E., Frenking, G., Kim, K.S., and Scuseria, G.E. (2005). Chapter 10—Finding minima, transition states, and following reaction pathways on ab initio potential energy surfaces. Theory and Applications of Computational Chemistry, Elsevier.
2. Current Status of Transition-State Theory;Truhlar;J. Phys. Chem.,1996
3. Formulation of the reaction coordinate;Fukui;J. Phys. Chem.,1970
4. Modeling the Kinetics of Bimolecular Reactions;Miller;Chem. Rev.,2006
5. Helgaker, T., Jørgensen, P., and Olsen, J. (2000). Molecular Electronic Structure Theory, John Wiley & Sons, Ltd.