Molecular dynamics simulations have been performed to study the liquid–vapor equilibrium of water as a function of temperature. The orthobaric densities and the surface tension of water are reported for temperatures from 316 K until 573 K. The extended simple point charge (SPC/E) interaction potential for water molecules is used with full Ewald summation. The normal and tangential components of the pressure tensor were calculated and are presented at 328 K. The nature of the long-range contribution to the surface tension has been studied in detail. At 328 K the calculated surface tension is 66.0±3.0 mN m−1 in comparison with the experimental value of 67 mN m−1. The simulated surface tensions between 316 K and 573 K are in good agreement with experiment. The orthobaric densities are in better agreement with experimental values than those obtained from the Gibbs ensemble calculation for the SPC model of water.

The gas-liquid surface of a system of Lennard-Jones (12, 6 ) molecules has been simulated by Monte Carlo and by Molecular Dynamic methods at temperatures which span most of the liquid range. For systems of 255 molecules the two methods lead to similar results and this agreement confirms that the density profile, as a function of height, falls monotonically from the density of the bulk liquid to that of the gas. The thickness of the surface layer is sensitive to the surface area, and appears to approach its thermodynamic limit for surface areas of 4000~ for a system of 4080 molecules. The density profile can be represented by a hyperbolic tangent of an appropriately scaled height. The thickness of the surface is of the order of two molecular diameters at temperatures near the triple point and increases rapidly as the critical point is approached. The computed surfacetens ions agree well with those calculated by statistical perturbation theory.

The liquid-vapor phase equilibria and surface tension of the TIP4P/2005 water model is obtained by using the Ewald summation method to determine the long range Lennard-Jones and electrostatic interactions. The method is implemented in a straightforward manner into standard simulation programs. The computational cost of using Ewald sums in dispersion interactions of water is estimated in direct simulation of interfaces. The results of this work at 300 K show a dramatic change in surface tension with an oscillatory behavior for surface areas smaller than 5x5sigma(2), where sigma is the Lennard-Jones oxygen diameter. The amplitude of such oscillations substantially decreases with temperature. Finite size effects are less important on coexisting densities. Phase equilibria and interfacial properties can be determined using a small number of water molecules; their fluctuations are around the same size of simulation error at all temperatures, even in systems where the interfaces are separated a few molecular diameters only. The difference in surface tension of this work compared to the results of other authors is not significant (on the contrary, there is a good agreement). What should be stressed is the different and more consistent approach to obtain the surface tension using the Ewald sums for dispersion interactions. There are two relevant aspects at the interface: An adsorption of water molecules is observed at small surface areas and its thickness systematically increases with system size.

The crystal formation of NaCl in water is studied by extensive molecular dynamics simulations. Ionic solutions at room temperature and various concentrations are studied using the SPC/E and TIP4P/2005 water models and seven force fields of NaCl. Most force fields of pure NaCl fail to reproduce the experimental density of the crystal, and in solution some favor dissociation at saturated conditions, while others favor crystal formation at low concentration. A new force field of NaCl is proposed, which reproduces the experimental phase diagram in the solid, liquid, and vapor regions. This force field overestimates the solubility of NaCl in water at saturation conditions when used with standard Lorentz-Berthelot combining rules for the ion-water pair potentials. It is shown that precipitation of ions is driven by the short range interaction between Cl-H pairs, a term which is generally missing in the simulation of ionic solutions. The effects of intramolecular flexibility of water on the solubility of NaCl ions are analyzed and is found to be small compared to rigid models. A flexible water model, extending the rigid SPC/E, is proposed, which incorporates Lennard-Jones interactions centered on the hydrogen atoms. This force field gives liquid-vapor coexisting densities and surface tensions in better agreement with experimental data than the rigid SPC/E model. The Cl-H, Na-O, and Cl-O pair distribution functions of the rigid and flexible models agree well with experiment. The predicted concentration dependence of the electric conductivity is in fair agreement with available experimental data.

A four-site rigid water model is presented, whose parameters are fitted to reproduce the experimental static dielectric constant at 298 K, the maximum density of liquid water and the equation of state at low pressures. The model has a positive charge on each of the three atomic nuclei and a negative charge located at the bisector of the HOH bending angle. This charge distribution allows increasing the molecular dipole moment relative to four-site models with only three charges and improves the liquid dielectric constant at different temperatures. Several other properties of the liquid and of ice Ih resulting from numerical simulations with the model are in good agreement with experimental values over a wide range of temperatures and pressures. Moreover, the model yields the minimum density of supercooled water at 190 K and the minimum thermal compressibility at 310 K, close to the experimental values. A discussion is presented on the structural changes of liquid water in the supercooled region where the derivative of density with respect to temperature is a maximum.