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.     

Vapor+liquid equilibrium of water, carbon dioxide, and the binary system, water+carbon dioxide, from molecular simulation

NVT- and NpT-Gibbs ensemble Monte Carlo (GEMC) simulations were applied to describe the vapor–liquid equilibrium of water (between 323 and 573 K), carbon dioxide (between 230 and 290 K) and their binary mixtures (between 348 and 393 K). The properties of supercritical carbon dioxide were determined between 310 and 520 K by NpT-MC simulations. Literature data for the effective pair potentials (for water: the SPC-, SPC/E-, and TIP4P potential models; for carbon dioxide: the EPM2 potential model) were used to describe the properties of the pure substances. The vapor pressures of water and carbon dioxide are calculated. For water, the SPC- and TIP4P models give superior results for the vapor pressure when compared to the SPC/E model. The vapor–liquid equilibrium of the binary mixture, carbon dioxide–water, was predicted using the SPC- as well as the TIP4P model for water and the EPM2 model for carbon dioxide. The interactions between carbon dioxide and water were estimated from the pair potentials of the pure components using common mixing rules without any adjustable binary parameter. Agreement of the predicted data for the compositions of the coexisting phases in vapor–liquid equilibrium and experimental results is observed within the statistical uncertainties of the simulation results in the investigated range of state, i.e. at pressures up to about 20 MPa.     

Determining the bulk viscosity of rigid water models

We use equilibrium molecular dynamics methods to compute the shear and bulk viscosities of the pairwise additive and rigid SPC/E, TIP4P, and TIP4P/2005 water models. For the latter model it was found in a recent study (J. Chem. Phys. 2009, 131, 246101) an excellent agreement with experiment in the prediction of the shear viscosity over a range of different thermodynamic conditions. Here, we examine, for a wide range of temperatures, whether this remarkable accuracy of the TIP4P/2005 model remains in the prediction of the bulk viscosity. Moreover, we examine whether equilibrium molecular dynamics methods provide reasonable accuracy in the calculation of the bulk viscosity, as it was previously found for the shear viscosity (J. Chem. Phys. 2010, 132, 096101). We concluded that, by performing the appropriate data analysis, accurate estimates of the bulk viscosity can be obtained, while, compared to the other simple rigid/pairwise additive water models, the predictions of the TIP4P/2005 model for the bulk viscosity are significantly closer to the experiment.     

Quantum effects in liquid water and ice: model dependence

This paper explores the influence of choice of potential model on the quantum effects observed in liquid water and ice. This study utilizes standard rigid models and a more formal context for the rigid-body centroid molecular dynamics methodology used to perform the quantum simulations is provided. Quantum and classical molecular dynamics simulations are carried out for liquid water and ice Ih at 298 and 220 K, respectively, with the simple point charge/extended and TIP4P-Ew water models. The results obtained for equilibrium and dynamical properties are compared with those recently reported on TIP4P [L. Hernandez de la Pena and P. G. Kusalik, J. Chem. Phys. 121, 5992 (2004); L. Hernandez de la Pena et al., J. Chem. Phys 123, 144506 (2005)]. For the liquid, an energy shift of about 8% and an average molecular uncertainty of about 11 degrees were found independently of the water model. The self-diffusion coefficient consistently increases by more than 50% when going from the classical to the quantum system and quantum dynamics are found to reproduce the experimental isotopic shifts with the models examined. The ice results compare remarkably well with those previously reported for the TIP4P water model; they confirm that quantum effects are considerable and that the quantum mechanical uncertainty and the energy shifts due to quantization are smaller in ice than in liquid water. The relevance of these findings in the context of the construction of water models is briefly discussed.     

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.     

Computational study of structural and dynamical properties of formamide-water mixtures

A molecular dynamics simulation study of structural and dynamical properties in liquid mixtures of formamide and water is presented. Site-site radial pair distribution functions, local mole fractions, pair energy distributions, and tetrahedral orientational order are the quantities analyzed to investigate the local structure in the simulated mixtures, along with a review of the intermolecular structure in terms of the distribution of hydrogen bonds. Our results indicate that there is a substitution of formamide molecules by water in the hydrogen bonds and a formation of a common hydrogen bond network. By analyzing the extent of tetrahedral order in the liquid as a function of composition, it is observed that whereas the tetrahedral network of liquid water is progressively lost by increasing the formamide concentration, the water structure within the first coordination shell is preserved and somewhat enhanced. The hydrogen-bond mean lifetimes were estimated by performing a time integration of the autocorrelation functions of bond occupation numbers. The lifetimes associated with hydrogen bonds between water, formamide, and interspecies pairs are found to increase with increasing formamide concentration. The lifetimes of the water hydrogen bonds show the largest variations, supporting the picture of an enhancement of the water structure among the nearest neighbors within the first coordination shell. We have used two different force field models for water, SPC/E [J. C. Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and TIP4P/2005 [J. L. F. Abascal and C. Vega, J. Chem. Phys. 123, 234505 (2005)]. Our results for structural and dynamical properties yield very small differences between those models, the TIP4P/2005 predicting a slightly more structured liquid and, consequently, exhibiting a slightly slower translational and librational dynamics.     

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.     

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.     

Simulation of phase separation in alcohol/water mixtures using two-body force field and standard molecular dynamics

Standard molecular dynamics simulations have been carried out on pure alcohols and alcohol/water mixtures. A simple atom-atom force field consisting of Lennard-Jones potentials plus coulombic terms over atomic point charges, but without explicit polarization terms, has been specifically fitted to reproduce several experimental properties of the pure alcohols, and has been used for mixtures by developing combination rules with the TIP3P water model. Densities, enthalpies of vaporization, radial distribution functions, self-diffusion coefficients, and rotational correlation functions of the pure alcohols are well reproduced and compare favorably with those from more sophisticated force fields. Some key aspects of the phase behaviour are correctly reproduced by the molecular dynamics simulation, showing a distinct demixing process for the n-butanol/water mixture as opposed to the stability of the t-butanol/water mixtures. The results demonstrate the ability of a molecular dynamics simulation, even in its standard form and with easily accessible time ranges, but with a carefully optimized force field, to simulate and, to a certain extent, predict the properties of binary mixtures.     

Structure and dynamics of liquid water with different long-range interaction truncation and temperature control methods in molecular dynamics simulations

We have used molecular dynamics simulations to study the physical properties of modified TIP3P water model included in the CHARMM program, using four different methods-the Ewald summation technique, and three different spherical truncation methods-for the treatment of the long-range interactions. Both the structure and dynamics of the liquid water model were affected by the methods used to truncate the long-range interactions. For some of the methods artificial structuring of the model liquid was observed around the cutoff radius. The model liquid properties were also affected by the commonly applied temperature control methods. Four different methods for controlling the temperature of the system were studied, and the effects of these methods on the bulk properties for liquid water were analyzed. The system size was also found to change the dynamics of the model liquid water. Two control simulations with the SPC/E water model were carried out. The self-diffusion coefficient (D), the radial distribution function (g(OO)), the distance dependent Kirkwood G-factor [G(k)(r)] and the intermolecular potential energy (E(pot)) were determined from the different trajectories and compared with the experimental data.     

Structural properties and thermodynamics of water clusters: a Wang-Landau study

The temperature dependence of structural properties and thermodynamic behavior of water clusters has been studied using Wang-Landau sampling. Four potential models, simple point charge/extended (SPC/E), transferable intermolecular potential 3 point (TIP3P), transferable intermolecular potential 4 point (TIP4P), and Gaussian charge polarizable (GCP), are compared for ground states and properties at finite temperatures. Although the hydrogen bond energy and the distance of the nearest-neighbor oxygen pair are significantly different for TIP4P and GCP models, they approach to similar ground state structures and melting transition temperatures in cluster sizes we considered. Comparing with TIP3P, SPC/E model provides properties closer to that of TIP4P and GCP.     

Effects of water model and simulation box size on protein diffusional motions

We performed molecular dynamics simulations of ubiquitin with the distinct water models, TIP3P, SPC, SPC/E, and SPC/Fw, in different system sizes with different box shapes. The translational diffusion constants of pure water linearly depend on the effective box length, which is known as finite size effect, whereas the first and second rotational times of pure water are nearly constant. We then observed that both the overall translation and rotational motions of the protein are linearly correlated to the viscosity of pure water. As expected from the finite size effect in pure water, the translational diffusion of the protein is significantly affected by the system size, and rotational diffusion is nearly size-independent. After correction for the finite size effect, the SPC/E and SPC/Fw models reproduce both the translational and rotational motion of the protein relatively well. Thus, water models that reproduce the experimentally derived diffusional properties of pure water accurately are expected to also be suitable for simulating protein diffusion quantitatively.     

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.     

Dynamical properties of the soft sticky dipole-quadrupole-octupole water model: a molecular dynamics study

The dynamical properties of the soft sticky dipole-quadrupole-octupole (SSDQO) water model using SPC/E moments are calculated utilizing molecular dynamics simulations. This new potential for liquid water describes the water-water interactions by a Lennard-Jones term and a sticky potential, which is an approximate moment expansion with point dipole, quadrupole, and octupole moments, and reproduces radial distribution functions of pure liquid water using the moments of SPC/E [Ichiye and Tan, J. Chem. Phys. 124, 134504 (2006)]. The forces and torques of SSDQO water for the dipole-quadrupole, quadrupole-quadrupole, and dipole-octupole interactions are derived here. The simulations are carried out at 298 K in the microcanonical ensemble employing the Ewald method for the long-range dipole-dipole interactions. Here, various dynamical properties associated with translational and rotational motions of SSDQO water using the moments of SPC/E (SSDQO:SPC/E) water are compared with the results from SPC/E and also experiment. The self-diffusion coefficient of SSDQO:SPC/E water is found to be in excellent agreement with both SPC/E and experiment whereas the single particle orientational relaxation time for dipole vector is better than SPC/E water but it is somewhat smaller than experiment. The dielectric constant of SSDQO:SPC/E is essentially identical to SPC/E, and both are slightly lower than experiment. Also, molecular dynamics simulations of the SSDQO water model are found to be about twice as fast as three-site models such as SPC/E.     

Solvation free energies of amino acid side chain analogs for common molecular mechanics water models

Quantitative free energy computation involves both using a model that is sufficiently faithful to the experimental system under study (accuracy) and establishing statistically meaningful measures of the uncertainties resulting from finite sampling (precision). In order to examine the accuracy of a range of common water models used for protein simulation for their solute/solvent properties, we calculate the free energy of hydration of 15 amino acid side chain analogs derived from the OPLS-AA parameter set with the TIP3P, TIP4P, SPC, SPC/E, TIP3P-MOD, and TIP4P-Ew water models. We achieve a high degree of statistical precision in our simulations, obtaining uncertainties for the free energy of hydration of 0.02-0.06 kcal/mol, equivalent to that obtained in experimental hydration free energy measurements of the same molecules. We find that TIP3P-MOD, a model designed to give improved free energy of hydration for methane, gives uniformly the closest match to experiment; we also find that the ability to accurately model pure water properties does not necessarily predict ability to predict solute/solvent behavior. We also evaluate the free energies of a number of novel modifications of TIP3P designed as a proof of concept that it is possible to obtain much better solute/solvent free energetic behavior without substantially negatively affecting pure water properties. We decrease the average error to zero while reducing the root mean square error below that of any of the published water models, with measured liquid water properties remaining almost constant with respect to our perturbations. This demonstrates there is still both room for improvement within current fixed-charge biomolecular force fields and significant parameter flexibility to make these improvements. Recent research in computational efficiency of free energy methods allows us to perform simulations on a local cluster that previously required large scale distributed computing, performing four times as much computational work in approximately a tenth of the computer time as a similar study a year ago.     

Prediction of thermodynamic, transport and vapor-liquid equilibrium properties of binary mixtures of ethylene glycol and water

Equilibrium and non-equilibrium molecular dynamics and Monte Carlo simulation techniques were applied to predict various thermodynamic, transport and vapor-liquid equilibrium properties of binary mixtures of ethylene glycol and water (EG-W) based on OPLS-AA and SPC/E force fields. The properties predicted include density, vaporization enthalpy, enthalpy of mixing, heat capacities, diffusion coefficients, shear viscosities, thermal conductivities, vapor-liquid coexistence isotherms and isobaric curves, and saturation vapor pressures. Good agreements with experimental data were obtained for most of these properties. Errors are mostly related to inaccuracy found in predictions of pure fluids; a correction to prediction of pure substance can systematically improve prediction for the mixture. This work suggests that OPLS-AA and SPC/E force fields using the common combining rules are transferable for predicting multiple physical properties of EG-W mixtures.     

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.