Classical interaction model for the water molecule

The authors propose a new classical model for the water molecule. The geometry of the molecule is built on the rigid TIP5P model and has the experimental gas phase dipole moment of water created by four equal point charges. The model preserves its rigidity but the size of the charges increases or decreases following the electric field created by the rest of the molecules. The polarization is expressed by an electric field dependent nonlinear polarization function. The increasing dipole of the molecule slightly increases the size of the water molecule expressed by the oxygen-centered sigma parameter of the Lennard-Jones interaction. After refining the adjustable parameters, the authors performed Monte Carlo simulations to check the ability of the new model in the ice, liquid, and gas phases. They determined the density and internal energy of several ice polymorphs, liquid water, and gaseous water and calculated the heat capacity, the isothermal compressibility, the isobar heat expansion coefficients, and the dielectric constant of ambient water. They also determined the pair-correlation functions of ambient water and calculated the energy of the water dimer. The accuracy of theirs results was satisfactory.