Surface properties of the polarizable Baranyai-Kiss water model

The water surface properties using the Baranyai-Kiss (BK) model [A. Baranyai and P. T. Kiss, J. Chem. Phys. 133, 144109 (2010)] are studied by molecular dynamics simulation, and compared to popular rigid water potentials, namely to the extended simple point charge (SPC/E) and the transferable interaction potential with 4 points (TIP4P) models. The BK potential is a polarizable model of water with three Gaussian charges. The negative charge is connected to its field-free position by a classical harmonic spring, and mechanical equilibrium is established between this spring force and the force due to the charge distribution of the system. The aim of this study is, on the one hand, to test the surface properties of the new model, and on the other hand, to identify differences between the models listed above. The obtained results reveal that the BK model reproduces very well a number of properties corresponding to liquid-vapor equilibrium, such as the coexisting liquid and vapor densities, saturated vapor pressure or surface tension. Further, this model reproduces excellently the critical point of water even in comparison with a large number of widely used polarizable and nonpolarizable models. The structural properties of the liquid surface of BK water turns out to be very similar to that of the SPC/E model, while the surface of TIP4P water is found to be somewhat less ordered. This finding is related to the fact that the critical temperature of the TIP4P model is lower than that of either SPC/E or BK.     

The melting temperature of the most common models of water

The melting temperature of ice I(h) for several commonly used models of water (SPC, SPC/E,TIP3P,TIP4P, TIP4P/Ew, and TIP5P) is obtained from computer simulations at p = 1 bar. Since the melting temperature of ice I(h) for the TIP4P model is now known [E. Sanz, C. Vega, J. L. F. Abascal, and L. G. MacDowell, Phys. Rev. Lett. 92, 255701 (2004)], it is possible to use the Gibbs-Duhem methodology [D. Kofke, J. Chem. Phys. 98, 4149 (1993)] to evaluate the melting temperature of ice I(h) for other potential models of water. We have found that the melting temperatures of ice I(h) for SPC, SPC/E, TIP3P, TIP4P, TIP4P/Ew, and TIP5P models are T = 190 K, 215 K, 146 K, 232 K, 245 K, and 274 K, respectively. The relative stability of ice I(h) with respect to ice II for these models has also been considered. It turns out that for SPC, SPC/E, TIP3P, and TIP5P the stable phase at the normal melting point is ice II (so that ice I(h) is not a thermodynamically stable phase for these models). For TIP4P and TIP4P/Ew, ice I(h) is the stable solid phase at the standard melting point. The location of the negative charge along the H-O-H bisector appears as a critical factor in the determination of the relative stability between the I(h) and II ice forms. The methodology proposed in this paper can be used to investigate the effect upon a coexistence line due to a change in the potential parameters.     

Relation between the melting temperature and the temperature of maximum density for the most common models of water

Water exhibits a maximum in density at normal pressure at 4 degrees above its melting point. The reproduction of this maximum is a stringent test for potential models used commonly in simulations of water. The relation between the melting temperature and the temperature of maximum density for these potential models is unknown mainly due to our ignorance about the melting temperature of these models. Recently we have determined the melting temperature of ice I(h) for several commonly used models of water (SPC, SPC/E, TIP3P, TIP4P, TIP4P/Ew, and TIP5P). In this work we locate the temperature of maximum density for these models. In this way the relative location of the temperature of maximum density with respect to the melting temperature is established. For SPC, SPC/E, TIP3P, TIP4P, and TIP4P/Ew the maximum in density occurs at about 21-37 K above the melting temperature. In all these models the negative charge is located either on the oxygen itself or on a point along the H-O-H bisector. For the TIP5P and TIP5P-E models the maximum in density occurs at about 11 K above the melting temperature. The location of the negative charge appears as a geometrical crucial factor to the relative position of the temperature of maximum density with respect to the melting temperature.     

Simulated surface tensions of common water models

Initial simulated values of the surface tension for the SPC/E water model have indicated excellent agreement with experiment. More recently, differing values have been obtained which are significantly lower than previous estimates. Here, we attempt to explain the differences between the previous studies and show that a variety of simulation conditions can affect the final surface tension values. Consistent values for the surface tensions of six common fixed charge water models (TIP3P, SPC, SPC/E, TIP4P, TIP5P, and TIP6P) are then determined for four temperatures between 275 and 350 K. The SPC/E and TIP6P models provide the best agreement with experiment.     

Prediction of self-diffusion coefficient and shear viscosity of water and its binary mixtures with methanol and ethanol by molecular simulation

Density, self-diffusion coefficient, and shear viscosity of pure liquid water are predicted for temperatures between 280 and 373 K by molecular dynamics simulation and the Green-Kubo method. Four different rigid nonpolarizable water models are assessed: SPC, SPC/E, TIP4P, and TIP4P/2005. The pressure dependence of the self-diffusion coefficient and the shear viscosity for pure liquid water is also calculated and the anomalous behavior of these properties is qualitatively well predicted. Furthermore, transport properties as well as excess volume and excess enthalpy of aqueous binary mixtures containing methanol or ethanol, based on the SPC/E and TIP4P/2005 water models, are calculated. Under the tested conditions, the TIP4P/2005 model gives the best quantitative and qualitative agreement with experiments for the regarded transport properties. The deviations from experimental data are of 5% to 15% for pure liquid water and 5% to 20% for the water + alcohol mixtures. Moreover, the center of mass power spectrum of water as well as the investigated mixtures are analyzed and the hydrogen-bonding structure is discussed for different states.     

Computer simulation of two new solid phases of water: Ice XIII and ice XIV

NpT Monte Carlo simulations have been performed for two recently discovered solid phases of water which have been denoted as ice XIII and ice XIV C. G. Salzmann et al. [Science311, 1758 (2006)]. Several potential models of water were considered, namely, the traditional SPC/E, TIP4P, and TIP5P and the more recent TIP5P-E, TIP4P-Ew, TIP4P/Ice, and TIP4P/2005 models. Significant differences in density and oxygen-oxygen radial distribution functions are found between the predictions of the SPC/E, TIP5P, and the models of the TIP4P family. The models TIP4P/Ice and TIP4P/2005 provide the best estimates of the density.     

Study of multipole contributions to the structure of water around ions in solution using the soft sticky dipole-quadrupole-octupole (SSDQO) model of water

The solvation of ions in the soft sticky dipole-quadrupole-octupole (SSDQO) model for liquid water is presented here. This new potential energy function for liquid water describes water-water interactions by a Lennard-Jones term plus a sticky potential consisting of an approximate moment expansion with point dipole, quadrupole, and octupole moments. The SSDQO potential energy function using the moments from extended simple point charge (SPC/E), TIP3P, or TIP5P reproduces the pair potential energy functions and radial distribution functions of the respective multipoint model but it is much faster than even the three-point models. Here, the solvation of ions in SSDQO water is studied using ion-water potential energy functions consisting of moment expansions up to the charge-quadrupole term, up to the charge-octupole term, and up to an approximate charge-hexadecapole term using the moments of SPC/E water. The radial distributions from Monte Carlo simulations show the best agreement with the results for ions in SPC/E water for the expansion up to the charge-hexadecapole term. Thus, the best results are obtained when the water-water and ion-water potentials are exact up to the 1r(4) term and also contain an approximate 1r(5) term. Overall, the simplicity, efficiency, and accuracy of the SSDQO potential make it potentially very useful for computer simulations of aqueous solvation.     

Sources of the deficiencies in the popular SPC/E and TIP3P models of water

Motivated by the results of Vega et al. [J. Phys. Condens. Matter 20, 153101 (2008)] about the phase diagram of water, and by the results of Kiss and Baranyai [J. Chem. Phys. 131, 204310 (2009)] about the properties of gas-phase clusters, we carried out a comparative study of the structure modeled by SPC∕E and TIP3P interactions in ambient liquid water. The gas-phase clusters of SPC∕E and TIP3P models show erroneous structures, while TIP4P-type models, either polarizable or not, provide qualitatively correct results. The trimers of SPC∕E and TIP3P are planar in gas phase, contrary to experimental and TIP4P-type models. The aim of this study was to see whether traces of these false geometries characteristic to SPC∕E and TIP3P in gas phase can also be found in the liquid phase. For this purpose we selected trimers formed by adjacent neighbors of water molecules in the liquid and calculated their geometrical features. We determined angles formed by the HO bonds of the molecules with OO vectors and with the normal vector of the OOO plane in the selected trimers. Our results showed that, despite high temperature, the SPC∕E and TIP3P water contains larger number of planar arrangements than other TIP4P-type models. Although structural differences presented in this study are small, they are accurately detectable. These results weaken the reliability of studies obtained by the SPC∕E or TIP3P models even in the liquid phase.     

A comparison of the value of viscosity for several water models using Poiseuille flow in a nano-channel

The viscosity-temperature relation is determined for the water models SPC/E, TIP4P, TIP4P/Ew, and TIP4P/2005 by considering Poiseuille flow inside a nano-channel using molecular dynamics. The viscosity is determined by fitting the resulting velocity profile (away from the walls) to the continuum solution for a Newtonian fluid and then compared to experimental values. The results show that the TIP4P/2005 model gives the best prediction of the viscosity for the complete range of temperatures for liquid water, and thus it is the preferred water model of these considered here for simulations where the magnitude of viscosity is crucial. On the other hand, with the TIP4P model, the viscosity is severely underpredicted, and overall the model performed worst, whereas the SPC/E and TIP4P/Ew models perform moderately.     

The melting point of ice Ih for common water models calculated from direct coexistence of the solid-liquid interface

In this work we present an implementation for the calculation of the melting point of ice I(h) from direct coexistence of the solid-liquid interface. We use molecular dynamics simulations of boxes containing liquid water and ice in contact. The implementation is based on the analysis of the evolution of the total energy along NpT simulations at different temperatures. We report the calculation of the melting point of ice I(h) at 1 bar for seven water models: SPC/E, TIP4P, TIP4P-Ew, TIP4P/ice, TIP4P/2005, TIP5P, and TIP5P-E. The results for the melting temperature from the direct coexistence simulations of this work are in agreement (within the statistical uncertainty) with those obtained previously by us from free energy calculations. By taking into account the results of this work and those of our free energy calculations, recommended values of the melting point of ice I(h) at 1 bar for the above mentioned water models are provided.     

Higher-order virial coefficients of water models

We use the Mayer sampling method, with both direct and overlap sampling, to calculate and compare classical virial coefficients up to B6 for various water models (SPC, SPC/E, MSPC/E, TIP3P, and TIP4P). The precision of the computed values ranges from 0.1% for B2 to an average of 25% for B6. When expressed in a form scaled by the critical properties, the values of the coefficients for SPC water are observed to greatly exceed the magnitude of corresponding coefficients for the simple Lennard-Jones model. We examine the coefficients in the context of the equation of state and the Joule-Thomson coefficient. Comparisons of these properties are made both to established molecular simulation data for each respective model and to real water. For all models, the virial series up to B5 describes the equation of state along the saturated vapor line better than the series that includes B6. At supercritical temperatures, however, the sixth-order series often describes pressure-volume-temperature behavior better than the fifth-order series. For example, the sixth-order virial equation of state for SPC/E water predicts the 673 K isotherm within 8% of published molecular simulation values up to a density of 9 mol/L (roughly half the critical density of SPC/E water).     

Monte Carlo calculation of the thermodynamic properties of water

The Monte Carlo method and parallel computing are used to calculate the thermodynamic properties of water (density, heat capacity, compressibility, thermal expansion coefficient, and static dielectric constant) in a wide range of temperatures (from 70 K to 530 K) at constant (atmospheric) pressure. Four groups of computational experiments are carried out, each for its own model of the water molecule: TIP3P (Jorgensen et al., 1983), SPC/E (Berendsen et al., 1987), TIP4P/2005 (Abascal&Vega, 2005), and TIP5P-E (Rick, 2004). An additional calculation based on the replica exchange method is conducted for the TIP4P/2005 model. A comparison of the calculated properties of water with experimental data suggests that the TIP4P/2005 model can provide highly realistic computer simulation results for water and aqueous solutions.     

Molecular dynamics simulation of liquid water and ice nanoclusters using a new effective HFD‐like model

We have determined a new two‐body interaction potential of water by the inversion of viscosity collision integrals of water vapor and fitted to achieve the Hartree–fock dispersion‐like (HFD‐like) potential function. The calculated two‐body potential generates the thermal conductivity, viscosity, and self‐diffusion coefficient of water vapor in an excellent accordance with experimental data at wide temperature ranges. We have also used a new many‐body potential as a function of temperature and density with the HFD‐like pair‐potential of water to improve the two‐body properties better than the SPC, SPC/E, TIP3P, and TIP4P models. We have also used the new corrected potential to simulate the configurational energy and the melting temperatures of the (H₂O)₅₀₀, (H₂O)₈₆₄, (H₂O)₂₀₄₈, and (H₂O)₆₉₁₂ ice nanoclusters in good agreement with the previous simulation data using the TIP4P model. The extrapolated melting point at the bulk limit is also in better agreement with the experimental bulk data. The self‐diffusion coefficients for the ice nanoclusters also simulated at different temperatures. © 2017 Wiley Periodicals, Inc.     

Hydration thermodynamic properties of amino acid analogues: a systematic comparison of biomolecular force fields and water models

We present an extensive study on hydration thermodynamic properties of analogues of 13 amino acid side chains at 298 K and 1 atm. The hydration free energies DeltaG, entropies DeltaS, enthalpies DeltaH, and heat capacities Deltac(P)() were determined for 10 combinations of force fields and water models. The statistical sampling was extended such that precisions of 0.3, 0.8, 0.8 kJ/mol and 25 J/(mol K) were reached for DeltaG, TDeltaS, DeltaH, and Deltac(P)(), respectively. The three force fields used in this study are AMBER99, GROMOS 53A6, and OPLS-AA; the five water models are SPC, SPC/E, TIP3P, TIP4P, and TIP4P-Ew. We found that the choice of water model strongly influences the accuracy of the calculated hydration entropies, enthalpies, and heat capacities, while differences in accuracy between the force fields are small. On the basis of an analysis of the hydrophobic analogues of the amino acid side chains, we discuss what properties of the water models are responsible for the observed discrepancies between computed and experimental values. The SPC/E water model performs best with all three biomolecular force fields.     

Hydration free energies of monovalent ions in transferable intermolecular potential four point fluctuating charge water: an assessment of simulation methodology and force field performance and transferability

Hydration free energies of nonpolarizable monovalent atomic ions in transferable intermolecular potential four point fluctuating charge (TIP4P-FQ) are computed using several commonly employed ion-water force fields including two complete model sets recently developed for use with the simple water model with four sites and Drude polarizability and TIP4P water models. A simulation methodology is presented which incorporates a number of finite-system free energy corrections within the context of constant pressure molecular dynamics simulations employing the Ewald method and periodic boundary conditions. The agreement of the computed free energies and solvation structures with previously reported results for these models in finite droplet systems indicates good transferability of ion force fields from these water models to TIP4Q-FQ even when ion polarizability is neglected. To assess the performance of the ion models in TIP4P-FQ, we compare with consensus values for single-ion hydration free energies arising from recently improved cluster-pair estimates and a reevaluation of commonly cited, experimentally derived single-ion hydration free energies; we couple the observed consistency of these energies with a justification of the cluster-pair approximation in assigning single-ion hydration free energies to advocate the use of these consensus energies as a benchmark set in the parametrization of future ion force fields.     

The importance of polarizability in the modeling of solubility: quantifying the effect of solute polarizability on the solubility of small nonpolar solutes in popular models of water

In recent work by Paschek [J. Chem. Phys. 120, 6674 (2004)] and others [see H. Docherty et al., J. Chem. Phys. 125, 074510 (2006) for a review] it has been suggested that, when coupled to a simple Lennard-Jones model for various small nonpolar solute molecules, the most common models of water (e.g., SPC/E and TIP4P) fail to reproduce quantitatively the solubility of small nonpolar solute molecules in water due in part to failing to account for polarization of the solute molecule. Given the importance of such systems as test-case prototype models of the solubility of proteins and biomolecules, in this work, we investigate the impact of using a polarizable solute model with the SPC/E, TIP3P, TIP4P, TIP4P-Ew, and TIP4P/2005 rigid water models. Specifically we consider Ne, Ar, Kr, Xe, and methane as solutes. In all cases we observe that the use of a polarizable solute improves agreement between experiment and simulations, with the best agreement seen for the largest solutes, Kr, CH(4), and Xe and the modern reparametrizations of the TIP4P model, i.e., the TIP4P-Ew and TIP4P/2005 models.     

Free energy of liquid water from a computer simulation via cell theory

A method to calculate the free energy of water from computer simulation is presented. Based on cell theory, it approximates the potential energy surface sampled in the simulation by an anisotropic six-dimensional harmonic potential to model the three hindered translations and three hindered rotations of a single rigid water molecule. The potential is parametrized from the magnitude of the forces and torques measured in the simulation. The entropy of these six harmonic oscillators is calculated and summed with a conformational term to give the total entropy. Combining this with the simulation enthalpy yields the free energy. The six water models examined are TIP3P, SPC, TIP4P, SPC/E, TIP5P, and TIP4P-Ew. The results reproduce experiment well: free energies for all models are within 1.6 kJ mol(-1) and entropies are within 3.6 J K(-1) mol(-1). Approximately two-thirds of the entropy comes from translation, a third from rotation, and 5% from conformation. Vibrational frequencies match those in the experimental infrared spectrum and assist in their assignment. Intermolecular quantum effects are found to be small, with free energies for the classical oscillator lying 0.5-0.7 kJ mol(-1) higher than in the quantum case. Molecular displacements and vibrational and zero point energies are also calculated. Altogether, these results validate the harmonic oscillator as a quantitative model for the liquid state.