Electrostatically Embedded Multiconfiguration Molecular Mechanics Based on the Combined Density Functional and Molecular Mechanical Method

Masahiro Higashi and Donald G. Truhlar*
Department of Chemistry and Supercomputing Institute, 207 Pleasant Street SE, University of Minnesota, Minneapolis, Minnesota 55455-0431
J. Chem. Theory Comput., 2008, 4 (5), pp 790–803
DOI: 10.1021/ct800004y
Publication Date (Web): March 26, 2008
Copyright © 2008 American Chemical Society
* Corresponding author e-mail: truhlar@umn.edu.

Abstract

We present a new method for generating global or semiglobal potential energy surfaces in the presence of an electrostatic potential; the new method can be used to model chemical reactions in solution or in an enzyme, nanocavity, or other chemical environment. The method extends the multiconfiguration molecular mechanics method so that the energy depends on the electrostatic potential at each atomic center. The charge distribution of the system can also be calculated. We illustrate the method by applying it to the symmetric bimolecular reaction Cl + CH3Cl′ → ClCH3 + Cl′ in aqueous solution, where the potential energy information is obtained by the combined density functional and molecular mechanical method, that is, by the combined quantum mechanical and molecular mechanical method (QM/MM) with the QM level being density functional theory. It is found that we can describe a semiglobal potential energy surface in aqueous solution with electronic structure information obtained entirely in the gas phase, including the linear and quadratic responses to variations in the electrostatic potential distribution. The semiglobal potential energy surface calculated by the present method is in good agreement with that calculated directly without any fitting.

Tools

History

  • Published In Issue May 13, 2008
  • Article ASAPMarch 26, 2008
  • Received: January 4, 2008

Recommend & Share

Related Content

Other ACS content by these authors: