A molecular dynamics study of water nucleation using the TIP4P/2005 model

Extensive molecular dynamics simulations were conducted using the TIP4P/2005 water model of Abascal and Vega [J. Chem. Phys. 123, 234505 (2005)] to investigate its condensation from supersaturated vapor to liquid at 330 K. The mean first passage time method [J. Wedekind, R. Strey, and D. Reguera, J. Chem. Phys. 126, 134103 (2007); L. S. Bartell and D. T. Wu, 125, 194503 (2006)] was used to analyze the influence of finite size effects, thermostats, and charged species on the nucleation dynamics. We find that the Nose?-Hoover thermostat and the one proposed by Bussi et al. [J. Chem. Phys. 126, 014101 (2007)] give essentially the same averages. We identify the maximum thermostat coupling time to guarantee proper thermostating for these simulations. The presence of charged species has a dramatic impact on the dynamics, inducing a marked change towards a pure growth regime, which highlights the importance of ions in the formation of liquid droplets in the atmosphere. It was found a small but noticeable sign preference at intermediate cluster sizes (between 5 and 30 water molecules) corresponding mostly to the formation of the second solvation shell around the ion. The TIP4P/2005 water model predicts that anions induce faster formation of water clusters than cations of the same magnitude of charge.     

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.     

Structural properties of water: comparison of the SPC, SPCE, TIP4P, and TIP5P models of water

Molecular-dynamics simulations were carried out for the SPC, SPCE, TIP4P, and TIP5P models of water at 298 K. From these results we determine the following quantities: the absolute entropy using the two-particle approximation, the mean lifetime of the hydrogen bond, the mean number of hydrogen bonds per molecule, and the mean energy of the hydrogen bond. From the entropy calculations we find that nearly all contributions to the total entropy originates from the orientation effects. Moreover, we determine the contributions to the total entropy which originate from the first, second, and higher solvation shells. It is interesting that the limits between solvation shells are clearly visible. The first solvation shell (0.22 < r < 0.36 nm) contributes approximately 43 J mol K to the total entropy; the second solvation shell (0.36 < r < 0.60 nm) contributes approximately 12 J mol K, while contributions from the third and other solvation shells are very small, approximately 2 J mol K in summary. This indicates that water molecules are strongly ordered up to 0.55-0.6 nm around the central water molecule, and beyond this limit the ordering diminishes. The results of calculations (entropy and hydrogen bonds) are compared with the experimental data for the choosing of the best water model. We find that the SPC and TIP4P models reproduce the best experimental values, and we recommend these models for computer simulations of the aqueous solution of biomolecules.     

A flexible model for water based on TIP4P/2005

A new flexible water model, TIP4P/2005f, is developed. The idea was to add intramolecular degrees of freedom to the successful rigid model TIP4P/2005 in order to try to improve the predictions for some properties, and to enable the calculation of new ones. The new model incorporates flexibility by means of a Morse potential for the bond stretching and a harmonic term for the angle bending. The parameters have been fitted to account for the peaks of the infrared spectrum of liquid water and to produce an averaged geometry close to that of TIP4P/2005. As for the intermolecular interactions, only a small change in the σ parameter of the Lennard-Jones potential has been introduced. The overall predictions are very close to those of TIP4P/2005. This ensures that the new model may be used with the same confidence as its predecessor in studies where a flexible model is advisable.     

A general purpose model for the condensed phases of water: TIP4P/2005

A potential model intended to be a general purpose model for the condensed phases of water is presented. TIP4P/2005 is a rigid four site model which consists of three fixed point charges and one Lennard-Jones center. The parametrization has been based on a fit of the temperature of maximum density (indirectly estimated from the melting point of hexagonal ice), the stability of several ice polymorphs and other commonly used target quantities. The calculated properties include a variety of thermodynamic properties of the liquid and solid phases, the phase diagram involving condensed phases, properties at melting and vaporization, dielectric constant, pair distribution function, and self-diffusion coefficient. These properties cover a temperature range from 123 to 573 K and pressures up to 40,000 bar. The model gives an impressive performance for this variety of properties and thermodynamic conditions. For example, it gives excellent predictions for the densities at 1 bar with a maximum density at 278 K and an averaged difference with experiment of 7 x 10(-4) g/cm3.     

Nonequilibrium versus equilibrium molecular dynamics for calculating the thermal conductivity of nanofluids

Journal of Thermal Analysis and Calorimetry - Both equilibrium and nonequilibrium molecular dynamics (EMD and NEMD, respectively) methods have been used to predict the thermal conductivity of...     

Widom line and the liquid-liquid critical point for the TIP4P/2005 water model

The Widom line and the liquid-liquid critical point of water in the deeply supercooled region are investigated via computer simulation of the TIP4P/2005 model. The Widom line has been calculated as the locus of compressibility maxima. It is quite close to the experimental homogeneous nucleation line and, in the region studied, it is almost parallel to the curve of temperatures of maximum density at fixed pressure. The critical temperature is determined by examining which isotherm has a region with flat slope. An interpolation in the Widom line gives the rest of the critical parameters. The computed critical parameters are T(c)=193 K, p(c)=1350 bar, and ρ(c)=1.012 g/cm(3). Given the performance of the model for the anomalous properties of water and for the properties of ice phases, the calculated critical parameters are probably close to those of real water.     

Equilibrium and nonequilibrium molecular dynamics simulations of the thermal conductivity of molten alkali halides

The thermal conductivity of molten NaCl and KCl was calculated through the Evans-Gillan nonequilibrium molecular dynamics (NEMD) algorithm and Green-Kubo equilibrium molecular dynamics (EMD) simulations. The EMD simulations were performed for a "binary" ionic mixture and the NEMD simulations assumed a pure system for reasons discussed in this work. The cross thermoelectric coefficient obtained from Green-Kubo EMD simulations is discussed in terms of the homogeneous thermoelectric power or Seebeck coefficient of these materials. The thermal conductivity obtained from NEMD simulations is found to be in very good agreement with that obtained through Green-Kubo EMD simulations for a binary ionic mixture. This result points to a possible cancellation between the neglected "partial enthalpy" contribution to the heat flux associated with the interdiffusion of one species through the other and that part of the thermal conductivity related to the coupled fluxes of charge and heat in "binary" ionic mixtures.     

Correlations in liquid water for the TIP3P-Ewald, TIP4P-2005, TIP5P-Ewald, and SWM4-NDP models

Water is one of the simplest molecules in existence, but also one of the most important in biological and engineered systems. However, understanding the structure and dynamics of liquid water remains a major scientific challenge. Molecular dynamics simulations of liquid water were performed using the water models TIP3P-Ewald, TIP4P-2005, TIP5P-Ewald, and SWM4-NDP to calculate the radial distribution functions (RDFs), the relative angular distributions, and the excess enthalpies, entropies, and free energies. In addition, lower-order approximations to the entropy were considered, identifying the fourth-order approximation as an excellent estimate of the full entropy. The second-order and third-order approximations are ~20% larger and smaller than the true entropy, respectively. All four models perform very well in predicting the radial distribution functions, with the TIP5P-Ewald model providing the best match to the experimental data. The models also perform well in predicting the excess entropy, enthalpy, and free energy of liquid water. The TIP4P-2005 and SWM4-NDP models are more accurate than the TIP3P-Ewald and TIP5P-Ewald models in this respect. However, the relative angular distribution functions of the four water models reveal notable differences. The TIP5P-Ewald model demonstrates an increased preference for water molecules to act both as tetrahedral hydrogen bond donors and acceptors, whereas the SWM4-NDP model demonstrates an increased preference for water molecules to act as planar hydrogen bond acceptors. These differences are not uncovered by analysis of the RDFs or the commonly employed tetrahedral order parameter. However, they are expected to be very important when considering water molecules around solutes and are thus a key consideration in modelling solvent entropy.     

Molecular dynamics simulations of the thermal conductivity of methane hydrate

Nonequilibrium molecular dynamics simulations with the nonpolarizable SPC/E (Berendsen et al., J. Phys. Chem. 1987, 91, 6269) and the polarizable COS/G2 (Yu and van Gunsteren, J. Chem. Phys. 2004, 121, 9549) force fields have been employed to calculate the thermal conductivity and other associated properties of methane hydrate over a temperature range from 30 to 260 K. The calculated results are compared to experimental data over this same range. The values of the thermal conductivity calculated with the COS/G2 model are closer to the experimental values than are those calculated with the nonpolarizable SPC/E model. The calculations match the temperature trend in the experimental data at temperatures below 50 K; however, they exhibit a slight decrease in thermal conductivity at higher temperatures in comparison to an opposite trend in the experimental data. The calculated thermal conductivity values are found to be relatively insensitive to the occupancy of the cages except at low (T相似文献   

A potential model for the study of ices and amorphous water: TIP4P/Ice

The ability of several water models to predict the properties of ices is discussed. The emphasis is put on the results for the densities and the coexistence curves between the different ice forms. It is concluded that none of the most commonly used rigid models is satisfactory. A new model specifically designed to cope with solid-phase properties is proposed. The parameters have been obtained by fitting the equation of state and selected points of the melting lines and of the coexistence lines involving different ice forms. The phase diagram is then calculated for the new potential. The predicted melting temperature of hexagonal ice (Ih) at 1 bar is 272.2 K. This excellent value does not imply a deterioration of the rest of the properties. In fact, the predictions for both the densities and the coexistence curves are better than for TIP4P, which previously yielded the best estimations of the ice properties.     

Dynamics of TIP5P and TIP4P/ice potentials

The dynamics of a thin film of ice Ih deposited on MgO (001) is studied through molecular dynamics simulations performed with two new potential models of ice. This system is chosen because it is possible to compare the results of the simulations to incoherent neutron quasielastic scattering experiments performed few years ago and to previous molecular dynamics simulations using the TIP4P potential model. The present simulations are performed to determine the evolution of the translational and orientational order parameters of the ice film upon temperature increase in the 250-280 K range. They are also used to calculate the translational and orientational diffusion coefficients of the water molecules in the supported film as a function of the temperature. When using the TIP5P potential, the present results show a better agreement with experimental data than those calculated with the TIP4P potential, especially regarding the temperature above which significant changes are obtained in the dynamics of the water film. Similar conclusions are obtained when using the TIP4P/ice potential, although this latter potential clearly underestimates the translational diffusion coefficients.     

Transport coefficients of the TIP4P-2005 water model

A detailed understanding of the dynamics of liquid water at molecular level is of fundamental importance as well as have applications in many branches of science and technology. In this work, the diffusion of the TIP4P-2005 model of water is systematically investigated in liquid phase in the temperature range 210-310 K. The translational and rotational diffusions, as well as correlations between them, are examined. The effects of system size and shape are also probed in this study. The results suggest the presence of a temperature of dynamical arrest of molecular translations in the range of 150-180 K and of molecular rotations in the range of 80-130 K, depending on specific direction. A substantial change in the preferred directions of translations and rotations relative to the molecular coordinate system is observed slightly below (≈15 K) the melting temperature of the model. It is shown that there is a correlation between translational and rotational molecular motions essential for diffusion in the liquid. The presence of hydrodynamic size effects is confirmed and quantified; it is also shown that using a non-cubic simulation box for a liquid system leads to an anisotropic splitting in the diffusion tensor. The findings of this study enhance our general understanding of models of water, specifically the TIP4P-2005 model, as well as provide evidences of the direct connection between thermodynamics of liquid water and dynamics of its molecules.     

Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

Influence of water-protein hydrogen bonding on the stability of Trp-cage miniprotein. A comparison between the TIP3P and TIP4P-Ew water models

We report extensive replica exchange molecular dynamics (REMD) simulations on the folding/unfolding equilibrium of Trp-cage miniprotein using the Amber ff99SB all atom forcefield and TIP3P and TIP4P-Ew explicit water solvent models. REMD simulation-lengths in the 500 ns to the microsecond regime per replica are required to adequately sample the folding/unfolding equilibrium. We observe that this equilibrium is significantly affected by the choice of the water model. Compared with experimental data, simulations using the TIP3P solvent describe the stability of the Trp-cage quite realistically, providing a melting point which is just a few Kelvins above the experimental transition temperature of 317 K. The TIP4P-Ew model shifts the equilibrium towards the unfolded state and lowers the free energy of unfolding by about 3 kJ mol(-1) at 280 K, demonstrating the need to fine-tune the protein-forcefield depending on the chosen water model. We report evidence that the main difference between the two water models is mostly due to the different solvation of polar groups of the peptide. The unfolded state of the Trp-cage is stabilized by an increasing number of hydrogen bonds, destabilizing the α-helical part of the molecule and opening the R-D salt bridge. By reweighting the strength of solvent-peptide hydrogen bonds by adding a hydrogen bond square well potential, we can fully recover the effect of the different water models and estimate the shift in population as due to a difference in hydrogen bond-strength of about 0.4 kJ mol(-1) per hydrogen bond.     

Chemometric study of liquid water simulations. I. The parameters of the TIP4P model potential

The multivariate chemometric techniques two level factorial design (TLFD) and principal component analysis (PCA) were used to investigate the TIP4P model potential behavior with respect to perturbations on all intermolecular interaction parameters. The effects of these perturbations were calculated for the enthalpy of vaporization, the density, the first maximum of the radial distribution functions of the O-H and O-O pairs, and the second maximum of the radial distribution function of the O-H pair obtained from Monte Carlo simulations of liquid water at 25 degrees C. The principal effects were quantified and rationalized in terms of the pair-wise interaction potential of the TIP4P model. They also corroborate previously published sensitivity analysis results using molecular dynamics and other model potentials. In addition, significant interaction effects between some parameters of the TIP4P model potential were observed and quantified, which hardly could be obtained without such a statistic approach. These interaction effects are very regular and systematic, and their behavior has not been encountered in other chemometric studies and cannot be rationalized in terms of the functional form of the pair-wise potential.     

Molecular dynamics simulations of vapor/liquid coexistence using the nonpolarizable water models

The surface tension, vapor-liquid equilibrium densities, and equilibrium pressure for common water models were calculated using molecular dynamics simulations over temperatures ranging from the melting to the critical points. The TIP4P/2005 and TIP4P-i models produced better values for the surface tension than the other water models. We also examined the correlation of the data to scaling temperatures based on the critical and melting temperatures. The reduced temperature (T/T(c)) gives consistent equilibrium densities and pressure, and the shifted temperature T + (T(c, exp) - T(c, sim)) gives consistent surface tension among all models considered in this study. The modified fixed charge model which has the same Lennard-Jones parameters as the TIP4P-FQ model but uses an adjustable molecular dipole moment is also simulated to find the differences in the vapor-liquid coexistence properties between fixed and fluctuating charge models. The TIP4P-FQ model (2.72 Debye) gives the best estimate of the experimental surface tension. The equilibrium vapor density and pressure are unaffected by changes in the dipole moment as well as the surface tension and liquid density.