These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
B68_water_model.pdf (281.98 kB)
An Efficient and Accurate Model for Water with an Improved Non-Bonded Potential
Preprints are manuscripts made publicly available before they have been submitted for formal peer review and publication. They might contain new research findings or data. Preprints can be a draft or final version of an author's research but must not have been accepted for publication at the time of submission.
A molecular mechanical model for liquid water is developed that uses a physically-motivated potential to represent Pauli repulsion and dispersion instead of the standard Lennard-Jones potential. The model has three-atomic sites and a virtual site located on the ∠HOH bisector (i.e., a TIP4P-type model). Pauli-repulsive interactions are represented using a Buckingham-type exponential decay potential. Dispersion interactions are represented by both and terms. This higher order dispersion term has been neglected by most force fields. The ForceBalance code was used to define parameters that optimally reproduce the experimental physical properties of liquid water. The resulting model is in good agreement with the experimental density, dielectric constant, enthalpy of vaporization, isothermal compressibility, thermal expansion coefficient, diffusion coefficient, and radial distribution function. A GPU-accelerated implementation of this improved non-bonded potential can be employed in OpenMM without modification by using the CustomNonBondedForce feature. Efficient and automated parameterization of these non-bonded potentials provides a rational strategy to define a new molecular mechanical force field that treats repulsion and dispersion interactions more rigorously without major modifications to existing simulation codes or a substantially larger computational cost.