Density and bond-orientational relaxations in supercooled water

ABSTRACT

Recent computational studies have reported evidence of a metastable liquid–liquid phase transition (LLPT) in molecular models of water under deeply supercooled conditions. A competing hypothesis suggests, however, that non-equilibrium artefacts associated with coarsening of the stable crystal phase have been mistaken for an LLPT in these models. Such artefacts are posited to arise due to a separation of time scales in which density fluctuations in the supercooled liquid relax orders of magnitude faster than those associated with bond-orientational order. Here, we use molecular simulation to investigate the relaxation of density and bond-orientational fluctuations in three molecular models of water (ST2, TIP5P and TIP4P/2005) in the vicinity of their reported LLPT. For each model, we find that density is the slowly relaxing variable under such conditions. We also observe similar behaviour in the coarse-grained mW model of water. Our findings, therefore, challenge the key physical assumption underlying the competing hypothesis.     

利用扩展系综法计算水的自由能

利用扩展系综法得到了正则系综下水的TIP4P模型的自由能值为-21．485±0．035kJ/mol， 并与其他方法所得的结果作了比较.提出了选择该方法中关键参数(平衡因子)的有效方法， 并讨论了该方法的可移植性. 关键词： 自由能 TIP4P水模型 扩展系综 分子动力学模拟 水分子团     

More accurate X-ray scattering data of deeply supercooled bulk liquid water

Deeply supercooled water droplets held containerless in an acoustic levitator are investigated with high-energy X-ray scattering. The temperature dependence of the X-ray structure function is found to be nonlinear. Comparison with two popular computer models reveals that structural changes are predicted too abrupt by the TIP5P-E model, while the rate of change predicted by TIP4P-Ew is in much better agreement with experiment. The abrupt structural changes, predicted by the TIP5P-E model to occur in the temperature range between 260 and 240?K as water approaches the homogeneous nucleation limit, are unrealistic. Both models underestimate the distance between neighbouring oxygen atoms and overestimate the sharpness of the OO distance distribution.     

Transport properties of bulk water at 243–550 K: a Comparative molecular dynamics simulation study using SPC/E,TIP4P,and TIP4P/2005 water models

ABSTRACT

Molecular dynamics simulations of various water models – SPC/E (extended simple point charge), TIP4P (transferable intermolecular potential 4 points), and TIP4P/2005 – have been carried out in the canonical (NVT fixed) ensemble over the range of temperatures 243–550?K with Ewald summation. The transport properties (self-diffusion coefficients D, viscosities η, and thermal conductivities λ) of SPC/E, TIP4P, and TIP4P/2005 water were evaluated at 243–550?K and compared with experimental data. The temperature dependence of transport properties of SPC/E, TIP4P and TIP4P/2005 water was discussed to determine how reliable the models are over this temperature range.     

Temperature dependent anomalous fluctuations in water: shift of ≈1 kbar between experiment and classical force field simulations

Here we report on the temperature dependence of the anomalous behaviour of water in terms of (i) its growth in tetrahedral structures, (ii) instantaneous spatial correlations from small angle x-ray scattering (SAXS) data, (iii) estimates of thermodynamic response functions of isothermal compressibility and (iv) thermal expansion coefficient. Water’s thermal expansion coefficient is estimated for the first time at supercooled conditions from liquid water’s structure factor. We used previously published data from classical force-fields of TIP4P/2005 and iAMOEBA to compare experimental data with molecular dynamics simulations and observe that these force-fields underestimate water’s anomalous behaviour but perform better upon increasing pressure. We demonstrate that the molecular dynamics simulations can describe better the temperature dependent anomalous behaviour of ambient pressure water if simulated at 1?kbar. The deviation in anomalous fluctuations in the simulations is not restricted to ≈228?K but extends all the way to ambient temperatures.     

Dynamical Fracture in Ice Ih as Modeled by TIP4P/Ice and mW Water Potentials

Fracture in ice I_h is simulated with molecular dynamics utilizing two potential fields, TIP4P/Ice and mW, and in different temperature conditions. The simulations produce propagating crack speeds over a large range of fracture energies. Terminal crack speed simulated with TIP4P/Ice potential can reach more than 200 m/s befitting experimental results. On the other hand, for mW potential, crack speed is around 5 m/s. The TIP4P/ice model suggests a brittle ice while mW potential describes a much more ductile material. The computational simulations are designed to permit direct comparison with experiments which can be performed in the hereafter. This comparison could provide a sensitive test to interatomic potentials.     

电场作用下过冷水形核的分子动力学模拟

采用分子动力学模拟的方法研究了过冷水的结构和形核过程随电场强度变化的情况。对系统微观结构的分析显示在电场强度为0～10~(10)V/m范围内,处于200K的过冷水在5 ns时间内不足以发生形核,只有当电场强度增加到10~(11) V/m或者温度降低到100K时系统才能在5 ns内发生形核。模拟结果表明电场强度的改变对系统中水分子团簇的微观结构有明显影响,电场强度的增加能促进系统的短程有序性。     

How the liquid-liquid transition affects hydrophobic hydration in deeply supercooled water

We determine the phase diagram of liquid supercooled water by extensive computer simulations using the TIP5P-E model. We find that the transformation of water into a low density liquid in the supercooled range strongly enhances the solubility of hydrophobic particles. The transformation of water into a tetrahedrally structured liquid is accompanied by a minimum in the hydration entropy and enthalpy. The corresponding change in sign of the solvation heat capacity indicates a loss of one characteristic signature of hydrophobic hydration. The observed behavior is found to be qualitatively in accordance with the predictions of the information theory model of Garde et al.     

噬菌体T4溶菌酶溶液构象模拟中分子力场和水模型的选择

蛋白质在溶液中可能以不同构象的集合形式存在,不能用单一的静态结构来表示. 分子动力学模拟已成为对溶液中蛋白质构象进行采样的有用工具,但分子力场和水模型的选择是关键问题. 这项工作介绍了噬菌体T4溶菌酶的个例研究. 本文发现,使用经典的AMBER99SB力场和TIP4P水模型,分子动力学模拟不能很好地描述野生型噬菌体T4溶菌酶在微秒时间尺度上的铰链弯曲结构域运动. 其它新型力场和水模型的组合,如被称为RSFF2+的残基特异性力场和离散校正的水模型TIP4P-D,能够对噬菌体T4溶菌酶溶液构象进行合理的采样,与实验数据有良好的一致性. 这项工作为进一步研究噬菌体T4溶菌酶的溶液构象转变提供了分子力场和水模型的参考.     

水滴撞击结冰过程的分子动力学模拟

利用三维分子动力学模拟方法,研究了纳米尺度水滴撞击冷壁面的结冰过程.数值模拟中,统计系统采用微正则系综,势能函数选用TIP4P/ice模型,温度校正使用速度定标法,牛顿运动方程的求解采用文莱特算法,水滴内部结冰过程则通过统计垂直方向水分子温度分布来判定.研究发现,当冷壁面温度降低时,水滴完全结冰的时间减小,但水滴降至壁面温度的时间却增大;同时随着壁面亲水性降低,水滴内部热传递速度减慢(尤其是冷壁面与水滴底端分子层间),水滴内部温度趋于均匀,但水滴完全结冰时间延长.     

水导热系数的分子动力学模拟

本文选用平衡分子动力学模拟方法,应用不同的势能模型对强极性分子水的导热系数进行了模拟计算.结果表明,用分子动力学模拟方法计算水的导热系数时,模拟结果对模型的选取极为敏感.在目前应用较为广泛的几种模型中以TIP4P模型较为适用,所得的导热系数与实验值较为接近.     

Molecular dynamics computation of the work of ion solvation: comparison of two models of water

Clusters of the TIP4P water model with 125 and 256 molecules and with the K⁺ or K^? ions were simulated by the molecular dynamics method. The radial profiles of local density, energy, normal pressure and molecular orientation were obtained. The work of cluster formation was calculated from the normal pressure profile. The results for the TIP4P and SPC/E model of water are discussed in comparison with each other. A divergence in behaviour of water molecules in the two models was discovered.     

超高压水的分子动力学模拟

 采用平衡分子动力学（EMD）方法,模拟研究了温度范围为243~348 K、压强范围为0.1~400 MPa条件下水的热力学性质、结构和动力学性质,模拟结果与实验值吻合较好。模拟结果表明,随着压强的增大,水分子间的氢键作用增强,扩散系数减小;随着温度的升高,水分子间的氢键作用减弱,有序程度下降,扩散系数增大。但在过冷水中,扩散系数随压强的增大有增加的趋势。     

Diffusion coefficient and shear viscosity of rigid water models

We report the diffusion coefficient and viscosity of popular rigid water models: two non-polarizable ones (SPC/E with three sites, and TIP4P/2005 with four sites) and a polarizable one (Dang-Chang, four sites). We exploit the dependence of the diffusion coefficient on the system size (Yeh and Hummer 2004 J. Phys. Chem. B 108 15873) to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang models overestimate the diffusion coefficient and underestimate the viscosity.     

A molecular dynamics simulation study on the cavitation inception of water with dissolved gases

The promotion/prevention mechanism of dissolved gases on cavitation inception is essential for many high-tech industries and research. In the present study, large-scale molecular dynamics simulations are performed to investigate the effects of water cavitation caused by different gas types by using nitrogen and oxygen gases with TIP4P/2005 water. The cavitation inception behaviour is analyzed via Mean First Passage Time method. Water with dissolved gases has a higher nucleation rate and is easier to cavitate than pure water. At the same gas concentration, the cavitation of water with nitrogen is promoted to a greater extent than that with oxygen. The number and energy of hydrogen bond (HB) are further calculated by the Acceptor-Hydrogen-Donor method to explain this promotion mechanism. The number and energy of HB in water with gases decrease compared with those in pure water. The introduction of gases weakens the HB network and promotes cavitation inception because of weaker interactions between gas and water molecules. A model is developed to describe the relationship between nucleation rate and HB energy. Gas molecules assemble on the surface of bubbles during water cavitation, which may decrease the free energy of bubble surface, maintain the existing bubble, and contribute to the growth process.     

Water and other tetrahedral liquids: order, anomalies and solvation

In order to understand the common features of tetrahedral liquids with water-like anomalies, the relationship between local order and anomalies has been studied using molecular dynamics simulations for three categories of such liquids: (a)?atomistic rigid-body models for water (TIP4P, TIP4P/2005, mTIP3P, SPC/E), (b)?ionic melts, BeF(2) (TRIM model) and SiO(2) (BKS potential) and (c)?Stillinger-Weber liquids parametrized to model water (mW) and silicon. Rigid-body, atomistic models for water and the Stillinger-Weber liquids show a strong correlation between tetrahedral and pair correlation order and the temperature for the onset of the density anomaly is close to the melting temperature. In contrast, the ionic melts show weaker and more variable degrees of correlation between tetrahedral and pair correlation metrics, and the onset temperature for the density anomaly is more than twice the melting temperature. In the case of water, the relationship between water-like anomalies and solvation is studied by examining the hydration of spherical solutes (Na(+), Cl(-), Ar) in water models with different temperature regimes of anomalies (SPC/E, TIP4P and mTIP3P). For both ionic and nonpolar solutes, the local structure and energy of water molecules is essentially the same as in bulk water beyond the second-neighbour shell. The local order and binding energy of water molecules are not perturbed by the presence of a hydrophobic solute. In the case of ionic solutes, the perturbation is largely localized within the first hydration shell. The binding energies for the ions are strongly dependent on the water models and clearly indicate that the geometry of the partial charge distributions, and the associated multipole moments, play an important role. However the anomalous behaviour of the water network has been found to be unimportant for polar solvation.     

Effect of short and long range forces on the thermodynamic properties of water. A simple short range reference system

Three realistic potential models of water, the non-polarizable ST2 and TIP4P models, and the polarizable TIP4P/P model, were used in computer simulations to study the effect of the range of intermolecular interactions on the thermodynamic properties of water. Following the results of recent studies, a short range system is constructed to the full pair potential u(1,2) in such a way that a perturbation expansion can be formulated in powers of the dipole-dipole interaction only. Computations of low density properties and computer simulations performed for several densities on three subcritical and one supercritical isotherms show that the short range reference not only reproduces the structure but approximates also the internal energy and pressure of water surprisingly well. Differences in the internal energy between the full and short range water do not exceed 5% for all models used over the entire range of the thermodynamic conditions considered.     

Molecular dynamics simulations of liquid water using various long-range electrostatics techniques

Water is one of the most extensively studied molecules, owing to its crucial role in biological processes. The water molecule is both highly polar and highly polarizable. Properties of water computed from molecular simulations are therefore critically dependent on both the intermolecular potential and the method for computing long-range electrostatic corrections. In this paper, the effects of the potential and the long-range electrostatic corrections are quantified for liquid water from 260 to 400?K. Simulations were carried out for a system of 256 molecules in the NVT ensemble. Thermodynamic, structural, dynamical, hydrogen bonding and dielectric properties have been computed for the flexible SPC and rigid SPC, SPC/E, TIP4P, TIP4P-Ew and TIP4P-FQ potentials, using the Lekner, Ewald and reaction field techniques to handle long-range electrostatics. The Lekner method gave the best overall agreement with experimental data, while the reaction field approach produced poorer results. Some measurable differences were found between the Lekner and Ewald techniques. For dielectric properties, the performance of the TIP4P-FQ model was superior relative to other potentials. For 256 molecules, the computational speeds of the Ewald and reaction field methods were found to be 2.5 to 3 times and 3.5 to 5 times faster than the Lekner technique, respectively.