Multiconfiguration Molecular Mechanics Based on Combined Quantum Mechanical and Molecular Mechanical Calculations
- 23 August 2006
- journal article
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 2  (5) , 1237-1254
- https://doi.org/10.1021/ct600171u
Abstract
The multiconfiguration molecular mechanics (MCMM) method is a general algorithm for generating potential energy surfaces for chemical reactions by fitting high-level electronic structure data with the help of molecular mechanical (MM) potentials. It was previously developed as an extension of standard MM to reactive systems by inclusion of multidimensional resonance interactions between MM configurations corresponding to specific valence bonding patterns, with the resonance matrix element obtained from quantum mechanical (QM) electronic structure calculations. In particular, the resonance matrix element is obtained by multidimensional interpolation employing a finite number of geometries at which electronic-structure calculations of the energy, gradient, and Hessian are carried out. In this paper, we present a strategy for combining MCMM with hybrid quantum mechanical molecular mechanical (QM/MM) methods. In the new scheme, electronic-structure information for obtaining the resonance integral is obtained by means of hybrid QM/MM calculations instead of fully QM calculations. As such, the new strategy can be applied to the studies of very large reactive systems. The new MCMM scheme is tested for two hydrogen-transfer reactions. Very encouraging convergence is obtained for rate constants including tunneling, suggesting that the new MCMM method, called QM/MM-MCMM, is a very general, stable, and efficient procedure for generating potential energy surfaces for large reactive systems. The results are found to converge well with respect to the number of Hessians. The results are also compared to calculations in which the resonance integral data are obtained by pure QM, and this illustrates the sensitivity of reaction rate calculations to the treatment of the QM-MM border. For the smaller of the two systems, comparison is also made to direct dynamics calculations in which the potential energies are computed quantum mechanically on the fly.Keywords
This publication has 149 references indexed in Scilit:
- Coupled Cluster/Molecular Mechanics Method: Implementation and Application to Liquid WaterThe Journal of Physical Chemistry A, 2003
- All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of ProteinsThe Journal of Physical Chemistry B, 1998
- A Quantum Mechanical and Molecular Mechanical Method Based on CM1A Charges: Applications to Solvent Effects on Organic Equilibria and ReactionsThe Journal of Physical Chemistry B, 1998
- Potential Energy Surfaces for Chemical Reactions: An Analytical Representation from Coarse Grained Data with an Application to Proton Transfer in WaterThe Journal of Physical Chemistry B, 1997
- Merck molecular force field. V. Extension of MMFF94 using experimental data, additional computational data, and empirical rulesJournal of Computational Chemistry, 1996
- Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94Journal of Computational Chemistry, 1996
- 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 moleculesComputer Physics Communications, 1995
- An examination of a density functional/molecular mechanical coupled potentialJournal of Computational Chemistry, 1995
- POLYRATE 4: A new version of a computer program for the calculation of chemical reaction rates for polyatomicsComputer Physics Communications, 1992
- THE FIVE-ELECTRON PROBLEM IN QUANTUM MECHANICS AND ITS APPLICATION TO THE HYDROGEN-CHLORINE REACTIONJournal of the American Chemical Society, 1932