Robust three-body water simulation model.

The most common potentials used in classical simulations of liquid water assume a pairwise additive form. Although these models have been very successful in reproducing many properties of liquid water at ambient conditions, none is able to describe accurately water throughout its complicated phase diagram. The primary reason for this is the neglect of many-body interactions. To this end, a simulation model with explicit three-body interactions was introduced recently [R. Kumar and J. L. Skinner, J. Phys. Chem. B 112, 8311 (2008)]. This model was parameterized to fit the experimental O-O radial distribution function and diffusion constant. Herein we reparameterize the model, fitting to a wider range of experimental properties (diffusion constant, rotational correlation time, density for the liquid, liquid/vapor surface tension, melting point, and the ice Ih density). The robustness of the model is then verified by comparing simulation to experiment for a number of other quantities (enthalpy of vaporization, dielectric constant, Debye relaxation time, temperature of maximum density, and the temperature-dependent second and third virial coefficients), with good agreement.