Multiconfiguration molecular mechanics algorithm for potential energy surfaces of chemical reactions

We present an efficient algorithm for generating semiglobal potential energy surfaces of reactive systems. The method takes as input molecular mechanics force fields for reactants and products and a quadratic expansion of the potential energy surface around a small number of geometries whose locations are determined by an iterative process. These Hessian expansions might come, for example, from ab initio electronic structure calculations, density functional theory, or semiempirical molecular orbital theory. A 2×2 electronic diabatic Hamiltonian matrix is constructed from these data such that, by construction, the lowest eigenvalue of this matrix provides a semiglobal approximation to the lowest electronically adiabatic potential energy surface. The theory is illustrated and tested by applications to rate constant calculations for three gas-phase test reactions, namely, the isomerization of 1,3-cis-pentadiene, OH+CH4→H2O+CH3, and CH2Cl+CH3F→CH3Cl+CH2F.

[1]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[2]  Michael A. Collins,et al.  POLYATOMIC MOLECULAR POTENTIAL ENERGY SURFACES BY INTERPOLATION IN LOCAL INTERNAL COORDINATES , 1998 .

[3]  Donald L Thompson,et al.  Modern Methods for Multidimensional Dynamics Computations in Chemistry , 1998 .

[4]  P. Halvick,et al.  Ab initio quasidiabatic states for the reaction N + CH → NC + H , 1997 .

[5]  W. Thiel,et al.  Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches , 1996 .

[6]  D. Truhlar,et al.  Reaction-path dynamics in curvilinear internal coordinates including torsions , 1996 .

[7]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[8]  D. Truhlar,et al.  A dual‐level Shepard interpolation method for generating potential energy surfaces for dynamics calculations , 1995 .

[9]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[10]  D. Truhlar,et al.  Reaction‐path potential and vibrational frequencies in terms of curvilinear internal coordinates , 1995 .

[11]  Ruifeng Liu,et al.  Molecular mechanics (MM3) calculations on alkyl radicals , 1994, J. Comput. Chem..

[12]  A. Warshel,et al.  Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches , 1993 .

[13]  Àngels González-Lafont,et al.  Direct dynamics calculation of the kinetic isotope effect for an organic hydrogen-transfer reaction, including corner-cutting tunneling in 21 dimensions , 1993 .

[14]  Àngels González-Lafont,et al.  MORATE: a program for direct dynamics calculations of chemical reaction rates by semiempirical molecular orbital theory , 1993 .

[15]  Donald G. Truhlar,et al.  Molecular modeling of the kinetic isotope effect for the [1,5]-sigmatropic rearrangement of cis-1,3-pentadiene , 1993 .

[16]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[17]  Gillian C. Lynch,et al.  Variational Transition-State Theory with Multidimensional, Semiclassical, Ground-State Transmission Coefficients: Applications to Secondary Deuterium Kinetic Isotope Effects in Reactions Involving Methane and Chloromethane , 1992 .

[18]  Donald G. Truhlar,et al.  POLYRATE 4: A new version of a computer program for the calculation of chemical reaction rates for polyatomics , 1992 .

[19]  W. Miller,et al.  Classical trajectory studies of the molecular dissociation dynamics of formaldehyde: H2CO→H2+CO , 1992 .

[20]  Fernando Bernardi,et al.  Simulation of MC-SCF results on covalent organic multi-bond reactions: molecular mechanics with valence bond (MM-VB) , 1992 .

[21]  Arieh Warshel,et al.  Computer Modeling of Chemical Reactions in Enzymes and Solutions , 1991 .

[22]  D. Truhlar,et al.  An experimentally based family of potential energy surfaces for hydride transfer between NAD+ analogs , 1991 .

[23]  Arieh Warshel,et al.  Role of solvent reorganization energies in the catalytic activity of enzymes , 1991 .

[24]  William H. Miller,et al.  An empirical valence bond model for constructing global potential energy surfaces for chemical reactions of polyatomic molecular systems , 1990 .

[25]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 3. The van der Waals' potentials and crystal data for aliphatic and aromatic hydrocarbons , 1989 .

[26]  George C. Schatz,et al.  The analytical representation of electronic potential-energy surfaces , 1989 .

[27]  J. Stewart.  Optimization of parameters for semiempirical methods II. Applications , 1989 .

[28]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[29]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[30]  B. C. Garrett,et al.  Variational transition state theory and tunneling for a heavy–light–heavy reaction using an ab initio potential energy surface. 37Cl+H(D) 35Cl→H(D) 37Cl+35Cl , 1983 .

[31]  Donald G. Truhlar,et al.  Polyatomic canonical variational theory for chemical reaction rates. Separable‐mode formalism with application to OH+H2→H2O+H , 1982 .

[32]  B. C. Garrett,et al.  Variational Transition State Theory , 1980 .

[33]  A. Warshel,et al.  An empirical valence bond approach for comparing reactions in solutions and in enzymes , 1980 .

[34]  Donald G. Truhlar,et al.  Improved treatment of threshold contributions in variational transition-state theory , 1980 .

[35]  B. C. Garrett,et al.  Generalized transition state theory. Quantum effects for collinear reactions of hydrogen molecules and isotopically substituted hydrogen molecules , 1979 .

[36]  J. L. Kinsey,et al.  Franck–Condon factors in studies of dynamics of chemical reactions. V. Simple construction of quasiadiabatic potential energy surfaces and numerical evaluation of Franck–Condon integrals , 1979 .

[37]  Donald G. Truhlar,et al.  Generalized transition state theory. Bond energy-bond order method for canonical variational calculations with application to hydrogen atom transfer reactions , 1979 .

[38]  J. Hollas,et al.  The Spectra and Structure of Simple Free Radicals , 1972 .

[39]  C. E. Young,et al.  Energy Distribution Among Products of Exothermic Reactions. II. Repulsive, Mixed, and Attractive Energy Release , 1966 .

[40]  N. L. Allinger,et al.  Organic Quantum Chemistry. IX. The Ultraviolet Spectra of Unsaturated Hydrocarbons. Woodward's Rules1,2 , 1965 .

[41]  F. O. Ellison.  A Method of Diatomics in Molecules. I. General Theory and Application to H2O , 1963 .

[42]  SatoShin,et al.  A New Method of Drawing the Potential Energy Surface , 1955 .

[43]  Arieh Warshel,et al.  Computational Approaches to Biochemical Reactivity , 2002 .

[44]  J. Åqvist,et al.  Molecular Dynamics Simulations of Substrate Dephosphorylation by Low Molecular Weight Protein Tyrosine Phosphatase , 1999 .

[45]  D. Truhlar,et al.  Reaction-Path Dynamics in Redundant Internal Coordinates , 1998 .

[46]  D. Truhlar,et al.  Variational transition-state theory and semiclassical tunnelling calculations with interpolated corrections: a new approach to interfacing electronic structure theory and dynamics for organic reactions , 1994 .

[47]  Peter Pulay,et al.  Geometry optimization in redundant internal coordinates , 1992 .

[48]  J. Dodd,et al.  Dipole moments of highly excited vibrational states of HCN , 1988 .

[49]  Mark S. Gordon,et al.  Potential energy surfaces for polyatomic reaction dynamics , 1987 .

[50]  Michael Baer,et al.  Theory of chemical reaction dynamics , 1985 .

[51]  D. Truhlar.  Potential Energy Surfaces and Dynamics Calculations , 1981 .

[52]  Richard B. Bernstein,et al.  Atom - Molecule Collision Theory , 1979 .

[53]  G. Segal.  Semiempirical Methods of Electronic Structure Calculation , 1977 .

[54]  Donald G. Truhlar,et al.  EXACT TUNNELING CALCULATIONS. , 1971 .