Accounting for electronic polarization in nonpolarizable force fields

Abstract

The issues of electronic polarizability in molecular dynamics simulations are discussed. We argue that the charges of ionized groups in proteins, and charges of ions in conventional non-polarizable force fields such as CHARMM, AMBER, GROMOS, etc should be scaled by a factor about 0.7. Our model explains why a neglect of electronic solvation energy, which typically amounts to about a half of total solvation energy, in non-polarizable simulations with un-scaled charges can produce a correct result; however, the correct solvation energy of ions does not guarantee the correctness of ion-ion pair interactions in many non-polarizable simulations. The inclusion of electronic screening for charged moieties is shown to result in significant changes in protein dynamics and can give rise to new qualitative results compared with the traditional non-polarizable force field simulations. The model also explains the striking difference between the value of water dipole $\mu$~3D reported in recent ab initio and experimental studies with the value ${\mu }_{eff}$~2.3D typically used in the empirical potentials, such as TIP3P or SPC/E. It is shown that the effective dipole of water can be understood as a scaled value ${\mu }_{eff}=\mu /\sqrt{{\epsilon }_{el}}$, where ${\epsilon }_{el}$=1.78 is the electronic (high-frequency) dielectric constant of water. This simple theoretical framework provides important insights into the nature of the effective parameters, which is crucial when the computational models of liquid water are used for simulations in different environments, such as proteins, or for interaction with solutes.

53 pages, 6 figures, 3 tables

Cited by in corpus (30)

• A force field of Li${}^{+}$ , Na${}^{+}$ , K${}^{+}$, Mg${}^{2+}$, Ca${}^{2+}$, Cl${}^{-}$, and SO${}_4^{2-}$ in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions
• ESPResSo 4.0 -- An Extensible Software Package for Simulating Soft Matter Systems
• Simulations of room temperature ionic liquids: From polarizable to coarse-grained force fields
• Polarizable Mean-Field Model of Water for Biological Simulations with Amber and Charmm force fields
• The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: extension to the ions F${}^{-}$, Br${}^{-}$, I${}^{-}$, Rb${}^{+}$, Cs${}^{+}$
• Accurate biomolecular simulations account for electronic polarization
• Polarizable molecular interactions in condensed phase and their equivalent nonpolarizable models
• Scaled charges for ions: an improvement but not the final word for modeling electrolytes in water
• Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces , One Property
• Structural and Transport Properties of Li/S Battery Electrolytes: Role of the Polysulfide Species
• Giant thermoelectric response of nanofluidic systems driven by water excess enthalpy
• Molecular modeling of aqueous electrolytes at interfaces: effects of long-range dispersion forces and of ionic charge rescaling
• Improved general-purpose five-point model for water: TIP5P/2018
• Local ice-like structure at the liquid water surface
• Tensorial properties via the neuroevolution potential framework: Fast simulation of infrared and Raman spectra
• In Silico Seawater
• Adsorption isotherm and mechanism of $\mathrm{C}{\mathrm{a}}^{2+}$ binding to polyelectrolyte
• Impact of Charged Surfaces on the Structure and Dynamics of Polymer Electrolytes: Insights from Atomistic Simulations
• Polymer electrolytes in strong external electric fields: Modification of structure and dynamics
• Determining the atomic charge of calcium ion requires the information of its coordination geometry in an EF-hand motif
• On the transferability of three water models developed by adaptive force matching
• Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions
• How important is the dielectric constant in water modeling? Evaluation of the performance of the TIP4P/$\epsilon$ force field and its compatibility with the Joung-Cheatham NaCl model
• Modelling structure and ionic diffusion in a class of ionic liquid crystal-based solid electrolytes
• Polarizability Plays a Decisive Role in Modulating Association Between Molecular Cations and Anions
• Static field-gradient polarizabilities of small atoms and molecules in finite temperature
• First Observation of Electron Transfer Mediated Decay in Aqueous Solutions: A Novel Probe of Ion Pairing
• Effect of Succinonitrile on Ion Transport in PEO-based Lithium Ion Battery Electrolytes
• Charge Transfer in Classical Molecular Dynamics Simulations of Met-enkephalin: Improving Traditional Force Field with Data Driven Models
• Molecular dynamics study of perchloric acid using the extended Madrid-2019 force field