Molecular Mechanics for Chemical Reactions: A Standard Strategy for Using Multiconfiguration Molecular Mechanics for Variational Transition State Theory with Optimized Multidimensional Tunneling
Citations Over TimeTop 10% of 2001 papers
Abstract
Multiconfiguration molecular mechanics (MCMM) is an extension of molecular mechanics to chemically reactive systems. This dual-level method combines molecular mechanics potentials for the reactant and product configurations with electronic structure Hessians at the saddle point and a small number of nonstationary points to model the potential energy surface in the reaction swath region between reactants and products where neither molecular mechanics potential is valid. The resulting semiglobal potential energy surface is used as input for dynamics calculations of tunneling probabilities and variational transition state theory rate constants. In this paper, we present a standard strategy for applying MCMM to calculate rate constants for atom transfer reactions. In particular, we propose a general procedure for determining where to calculate the electronic structure Hessians. We tested this strategy for a diverse test suite of six reactions involving hydrogen-atom transfer. It yields reasonably accurate rate constants as compared to direct dynamics using an uninterpolated full potential energy surface at the same electronic structure level. Furthermore, the rate constants at each of several successively more demanding levels of dynamical theory are also predicted accurately, which indicates that the MCMM potential energy surface accurately predicts many different details of the potential energy surface with a limited number of electronic structure Hessians.
Related Papers
- → CASPT2 Study of the Potential Energy Surface of the HSO2 System(2011)30 cited
- → Saddle point localization of molecular wavefunctions(2016)8 cited
- A Second-Order Saddle Point in the Reaction Coordinate for the Isomerization of the NH 5 Complex: Ab initio Calculations(2004)
- → Second-Order Methods for the Optimization of Molecular Potential Energy Surfaces(1995)4 cited
- → Microdynamics study on reaction of atom H and radical NH (X3Σ)(1999)