Molecular dynamics study of water clusters, liquid, and liquid–vapor interface of water with many-body potentials

The molecular dynamics computer simulation technique is used to develop a rigid, four-site polarizable model for water. The suggested model reasonably describes the important properties of water clusters, the thermodynamic and structural properties of the liquid and the liquid/vapor interface of water. The minimum energy configurations and the binding energies for these clusters are in reasonable agreement with accurate electronic structure calculations. The model predicts that the water trimer, tetramer, and pentamer have cyclic planar minimum energy structures. A prismlike structure is predicted to be lowest in energy for the water hexamer, and a cagelike structure is the second lowest in energy, with an energy of about 0.2 kcal/mol higher than the prismlike structure. The results are consistent with recent quantum Monte Carlo simulations as well as electronic structure calculations. The computed thermodynamic properties for the model, at room temperature, including the liquid density, the enthalpy of vaporization, as well as the diffusion coefficient, are in excellent agreement with experimental values. Structural properties of liquid water, such as the radial distribution functions, neutron, and x-ray scattering intensities, were calculated and critically evaluated against the experimental measurements. In all cases, we found the agreement between the observed data and the computed properties to be quite reasonable. The computed density profile of the water’s liquid/vapor interface shows that the interface is not sharp at a microscopic level and has a thickness of 3.2 Å at 298 K. These results are consistent with those reported in earlier work on the same systems. The calculated surface tension at room temperature is in reasonable agreement with the corresponding experimental data. As expected, the computed average dipole moments of water molecules near the interface are close to their gas phase values, while water molecules far from the interface have dipole moments corresponding to their bulk values.