Benchmarking Electronic Structure Methods for Accurate Fixed-Charge Electrostatic Models
The accuracy of classical molecular mechanics (MM) force fields used for condensed phase molecular simulations depends strongly on the accuracy of modeling nonbonded interactions between atoms, such as electrostatic interactions. Some popular fixed-charge MM force fields use partial atomic charges derived from gas phase electronic structure calculations using the Hartree-Fock method with the relatively small 6-31G* basis set (HF/6-31G*). It is generally believed that HF/6-31G* generates fortuitously overpolarized electron distributions, as would be expected in the higher dielectric environment of the condensed phase. Using a benchmark set of 47 molecules we show that HF/6-31G* overpolarizes molecules by just under 10% on average with respect to experimental gas phase dipole moments. The overpolarization of this method/basis set combination varies significantly though and, in some cases, even leads to molecular dipole moments that are lower than experimental gas phase measurements. We further demonstrate that using computationally inexpensive density functional theory (DFT) methods, together with appropriate augmented basis sets and a continuum solvent model, can yield molecular dipole moments that are both more strongly and more uniformly overpolarized. These data suggest that these methods – or ones similar to them – should be adopted for the derivation of accurate partial atomic charges for next-generation MM force fields.
Read the published paper
in Journal of Chemical Information and Modeling
CALIFORNIA STATE UNIVERSITY LOS ANGELES PROPOSES TO ESTABLISH A NASA DATA INTENSIVE RESEARCH AND EDUCATION CENTER IN STEM DIRECT STEM TO CONDUCT COLLABORATIVE RESEARCH AND TRAINING IN AREAS OF HYDROLOGY AND CLIMATE CHANGE COMPUTATIONAL PHYSICS AND CLOUD C
National Aeronautics and Space AdministrationFind out more...
Efficient, Polarizable Charge Model for Molecular Simulation
FWF Austrian Science FundFind out more...