A Caltech Library Service

Polarizable charge equilibration model for predicting accurate electrostatic interactions in molecules and solids

Naserifar, Saber and Brooks, Daniel J. and Goddard, William A., III and Cvicek, Vaclav (2017) Polarizable charge equilibration model for predicting accurate electrostatic interactions in molecules and solids. Journal of Chemical Physics, 146 (12). Art. No. 124117. ISSN 0021-9606. doi:10.1063/1.4978891.

[img] PDF - Published Version
See Usage Policy.


Use this Persistent URL to link to this item:


Electrostatic interactions play a critical role in determining the properties, structures, and dynamics of chemical, biochemical, and material systems. These interactions are described well at the level of quantum mechanics (QM) but not so well for the various models used in force field simulations of these systems. We propose and validate a new general methodology, denoted PQEq, to predict rapidly and dynamically the atomic charges and polarization underlying the electrostatic interactions. Here the polarization is described using an atomic sized Gaussian shaped electron density that can polarize away from the core in response to internal and external electric fields, while at the same time adjusting the charge on each core (described as a Gaussian function) so as to achieve a constant chemical potential across all atoms of the system. The parameters for PQEq are derived from experimental atomic properties of all elements up to Nobelium (atomic no. = 102). We validate PQEq by comparing to QM interaction energy as probe dipoles are brought along various directions up to 30 molecules containing H, C, N, O, F, Si, P, S, and Cl atoms. We find that PQEq predicts interaction energies in excellent agreement with QM, much better than other common charge models such as obtained from QM using Mulliken or ESP charges and those from standard force fields (OPLS and AMBER). Since PQEq increases the accuracy of electrostatic interactions and the response to external electric fields, we expect that PQEq will be useful for a large range of applications including ligand docking to proteins, catalytic reactions, electrocatalysis, ferroelectrics, and growth of ceramics and films, where it could be incorporated into standard force fields as OPLS, AMBER, CHARMM, Dreiding, ReaxFF, and UFF.

Item Type:Article
Related URLs:
URLURL TypeDescription
Naserifar, Saber0000-0002-1069-9789
Goddard, William A., III0000-0003-0097-5716
Additional Information:© 2017 AIP Publishing. Published Online: March 2017 Accepted: March 2017. This work was supported by the Joint Center for Artificial Photosynthesis, a DOE Energy Innovation Hub, supported through the Office of Science of the U.S. DOE under Award No. DE-SC0004993. We thank Dr. Qingsong Zhang, Dr. Sergey Zybin, and Dr. Andres Jaramillo-Botero for useful discussions.
Funding AgencyGrant Number
Department of Energy (DOE)DE-SC0004993
Other Numbering System:
Other Numbering System NameOther Numbering System ID
Issue or Number:12
Record Number:CaltechAUTHORS:20170410-141114724
Persistent URL:
Official Citation:Polarizable charge equilibration model for predicting accurate electrostatic interactions in molecules and solids Saber Naserifar, Daniel J. Brooks, William A. Goddard III, and Vaclav Cvicek The Journal of Chemical Physics 146, 124117 (2017);
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:76482
Deposited By: Ruth Sustaita
Deposited On:10 Apr 2017 21:28
Last Modified:15 Nov 2021 17:00

Repository Staff Only: item control page