Liquid-liquid phase transitions in supercooled water studied by computer simulations of various water models

Liquid-liquid and liquid-vapor coexistence regions of various water models were determined by Monte Carlo (MC) simulations of isotherms of density fluctuation-restricted systems and by Gibbs ensemble MC simulations. All studied water models show multiple liquid-liquid phase transitions in the supercooled region: we observe two transitions of the TIP4P, TIP5P, and SPCE models and three transitions of the ST2 model. The location of these phase transitions with respect to the liquid-vapor coexistence curve and the glass temperature is highly sensitive to the water model and its implementation. We suggest that the apparent thermodynamic singularity of real liquid water in the supercooled region at about 228 K is caused by an approach to the spinodal of the first (lowest density) liquid-liquid phase transition. The well-known density maximum of liquid water at 277 K is related to the second liquid-liquid phase transition, which is located at positive pressures with a critical point close to the maximum. A possible order parameter and the universality class of liquid-liquid phase transitions in one-component fluids are discussed.     

Four phases of amorphous water: Simulations versus experiment

Multiplicity of the liquid-liquid phase transitions in supercooled water, first obtained in computer simulations [Brovchenko et al., J. Chem. Phys. 118, 9473 (2003)], has got strong support from the recent experimental observation of the two phase transitions between amorphous ices [Loerting et al., Phys. Rev. Lett. 96, 025702 (2006)]. These experimental results allow assignment of the four amorphous water phases (I-IV) obtained in simulations to the three kinds of amorphous ices. Water phase I (rho approximately 0.90 gcm(3)) corresponds to the low-density amorphous ice, phase III (rho approximately 1.10 gcm(3)) to the high-density amorphous ice, and phase IV (rho approximately 1.20 gcm(3)) to the very-high-density amorphous ice. Phase II of model water with density rho approximately 1.00 gcm(3) corresponds to the normal-density water. Such assignment is confirmed by the comparison of the structural functions of the amorphous phases of model water and real water. In phases I and II the first and second coordination shells are clearly divided. Phase I consists mainly of the four coordinated tetrahedrally ordered water molecules. Phase II is enriched with molecules, which have tetrahedrally ordered four nearest neighbors and up six molecules in the first coordination shell. Majority of the molecules in phase III still have tetrahedrally ordered four nearest neighbors. Transition from phase III to phase IV is characterized by a noticeable drop of tetrahedral order, and phase IV consists mainly of molecules with highly isotropic angular distribution of the nearest neighbors. Relation between the structures of amorphous water phases, crystalline ices, and liquid water is discussed.     

Dynamics of water confined in single- and double-wall carbon nanotubes

Using high-resolution quasielastic neutron scattering, we investigated the temperature dependence of single-particle dynamics of water confined in single- and double-wall carbon nanotubes with the inner diameters of 14+/-1 and 16+/-3 A, respectively. The temperature dependence of the alpha relaxation time for water in the 14 A nanotubes measured on cooling down from 260 to 190 K exhibits a crossover at 218 K from a Vogel-Fulcher-Tammann law behavior to an Arrhenius law behavior, indicating a fragile-to-strong dynamic transition in the confined water. This transition may be associated with a structural transition from a high-temperature, low-density (<1.02 gcm(3)) liquid to a low-temperature, high-density (>1.14 gcm(3)) liquid found in molecular dynamics simulation at about 200 K. However, no such dynamic transition in the investigated temperature range of 240-195 K was detected for water in the 16 A nanotubes. In the latter case, the dynamics of water simply follows a Vogel-Fulcher-Tammann law. This suggests that the fragile-to-strong crossover for water in the 16 A nanotubes may be shifted to a lower temperature.     

液-液溶解度法与气-液色谱法研究溶剂-溶剂的相互作用

测量了角鲨烷(B)+邻苯二甲酸二丁酯(C)的溶解度相图, 应用HSFH二元系公式推算得相互作用参数Abc; 用气-液色谱测量了若干正构烷烃、环烷烃用作探测溶质(A)在(B+C)混合溶剂中的无限稀活度系数, 应用HSFH三元系公式推算得相互作用参数Abc, 两者符合热力学一致性。     

Pressure effects on the transitions between disordered phases in supercooled liquid silicon

We investigate the pressure effects on the transitions between the disordered phases in supercooled liquid silicon through Monte Carlo simulations and efficient methods to compute free energies. Our calculations, using an environment dependent interatomic potential for Si, indicate that at zero pressure the liquid-liquid phase transition, between the high density liquid and the low density liquid, occurs at a temperature 325K below melting. We found that the liquid-liquid transition temperature decreases with increasing pressure, following the liquid-solid coexistence curve. As pressure increases, the liquid-liquid coexistence curve approaches the region where the glass transition between the low density liquid and the low density amorphous takes place. Above 5 GPa, our calculations show that the liquid-liquid transition is suppressed by the glassy dynamics of the system. We also found that above 5 GPa, the glass transition temperature is lower than that at lower pressures, suggesting that under these conditions the glass transition occurs between the high density liquid and the high density amorphous.     

A global investigation of phase equilibria using the perturbed-chain statistical-associating-fluid-theory approach

The recently developed perturbed-chain statistical-associating-fluid theory (PC-SAFT) is investigated for a wide range of model parameters including the parameter m representing the chain length and the thermodynamic temperature T and pressure p. This approach is based upon the first-order thermodynamic perturbation theory for chain molecules developed by Wertheim [M. S. Wertheim, J. Stat. Phys. 35, 19 (1984); ibid. 42, 459 (1986)] and Chapman et al. [G. Jackson, W. G. Chapman, and K. E. Gubbins, Mol. Phys. 65, 1 (1988); W. G. Chapman, G. Jackson, and K. E. Gubbins, ibid. 65, 1057 (1988)] and includes dispersion interactions via the second-order perturbation theory of Barker and Henderson [J. A. Barker and D. Henderson, J. Chem. Phys. 47, 4714 (1967)]. We systematically study a hierarchy of models which are based on the PC-SAFT approach using analytical model calculations and Monte Carlo simulations. For one-component systems we find that the analytical model in contrast with the simulation results exhibits two phase-separation regions in addition to the common gas-liquid coexistence region: One phase separation occurs at high density and low temperature. The second demixing takes place at low density and high temperature where usually the ideal-gas phase is expected in the phase diagram. These phenomena, which are referred to as "liquid-liquid" and "gas-gas" equilibria, give rise to multiple critical points in one-component systems, as well as to critical end points and equilibria of three fluid phases, which can usually be found in multicomponent mixtures only. Furthermore, it is shown that the liquid-liquid demixing in this model is not a consequence of a "softened" repulsive interaction as assumed in the theoretical derivation of the model. Experimental data for the melt density of polybutadiene with molecular mass Mw=45,000 gmol are correlated here using the PC-SAFT equation. It is shown that the discrepancies in modeling the polymer density at ambient temperature and high pressure can be traced back to the liquid-liquid phase separation predicted by the equation of state at low temperatures. This investigation provides a basis for understanding possible inaccuracies or even unexpected phase behavior which can occur in engineering applications of the PC-SAFT model aiming at predicting properties of macromolecular substances.     

The ice-vapor interface and the melting point of ice I(h) for the polarizable POL3 water model

We use molecular dynamics simulations to determine the melting point of ice I(h) for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys.1992, 97, 2659). Simulations are performed on a slab of ice I(h) with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields the melting point of the POL3 model to be T(m) = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.     

Comparison of selected polarizable and nonpolarizable water models in molecular dynamics simulations of ice I(h)

We present a molecular dynamics simulation study in which we determined the melting point of ice I(h) for the polarizable SWM4-NDP water model (Lamoureux et al., Chem. Phys. Lett., 2006, 418, 245-249) and compared the performance of several popular water force fields, both polarizable and nonpolarizable, in terms of melting temperature, stability and orientational structuring of ice. The simulations yield the melting temperature of SWM4-NDP ice as low as T(m) = 185 ± 10 K, despite the quadrupole moment of a SWM4-NDP water molecule being close to the experimental gas phase value. The results thus show that the dependence of T(m) on the molecular quadrupole, observed for the three- and four-site water models, is generally lost if polarization is explicitly included. The study also shows that adding polarizability to a planar three-charge water model increases orientational disorder in hexagonal ice. In addition, analysis of the tetrahedral order in bulk ice reveals a correlation between the pre-existing degree of orientational disorder in ice simulated using different polarizable and nonpolarizable models and the melting temperature of the models. Our findings thus suggest some new considerations regarding the role of polarization forces in a crystalline solid that may guide future development of reliable polarizable water models for ice.     

Phase behavior and local structure of a binary mixture in pores: mean-field lattice model calculations for analyzing neutron scattering data

We investigate the phase behavior of an asymmetric binary liquid A-W mixture confined between two planar homogenous substrates (slit pore). Molecules of species W interact preferentially with the solid walls via a long-range potential. Assuming nearest-neighbor attractions between the liquid molecules, we employ a lattice-gas model and a mean-field approximation for the grand potential. Minimization of this potential yields the density profiles of thermodynamically stable phases for fixed temperature, chemical potentials of both species, pore width and strengths of attraction. This model is used to analyze experimental small-angle neutron-scattering (SANS) data on the microscopic structure of the binary system isobutyric acid (iBA)+heavy water (D2O) inside a mesoscopic porous matrix (controlled-pore glass of about 10 nm mean pore width). Confinement-independent model parameters are adjusted so that the theoretical liquid-liquid coexistence curve in the bulk matches its experimental counterpart. By choosing appropriate values of the pore width and the attraction strength between substrates and water we analyze the effect of confinement on the phase diagram. In addition to a depression of the liquid-liquid critical point we observe surface induced phase transitions as well as water-film adsorption near the walls. The temperature dependence of the structure of water-rich and iBA-rich phases of constant composition are discussed in detail. The theoretical predictions are consistent with results of the SANS study and assist their interpretation.     

The isomerization of dinitrogen tetroxide: O2N-NO2 --> ONO-NO2

The N2O4 isomerization in gas phase has an energy barrier of 31 kcal mol-1 at 298 K. This energy barrier may be reduced due to the interaction of the N2O4 isomers with water or nitric acid clusters adsorbed on surfaces. The Gibbs free energy barrier for this reaction in water medium is estimated to be reduced to 21.1 kcal mol-1 by using the ab initio calculations and the polarizable continuum model (PCM). By using the transition state theory (TST), this model estimates that the N2O4 isomerization may be as fast as 2.0 x 10(-3) s-1 in aqueous phase at room temperature, which confirms the Finlayson-Pitts model for the heterogeneous hydrolysis of NO2 on surfaces. The activation energy of the N2O4 isomerization is about 21 kcal mol-1. The rate coefficient for this reaction is considerably fast, 1.2 x 10(-2) s-1, in aqueous phase at T = 373 K.     

Artificial multiple criticality and phase equilibria: an investigation of the PC-SAFT approach

The perturbed-chain statistical associating fluid theory (PC-SAFT) is studied for a wide range of temperature, T, pressure, p, and (effective) chain length, m, to establish the generic phase diagram of polymers according to this theory. In addition to the expected gas-liquid coexistence, two additional phase separations are found, termed "gas-gas" equilibrium (at very low densities) and "liquid-liquid" equilibrium (at densities where the system is expected to be solid already). These phase separations imply that in one-component polymer systems three critical points occur, as well as equilibria of three fluid phases at triple points. However, Monte Carlo simulations of the corresponding system yield no trace of the gas-gas and liquid-liquid equilibria, and we conclude that the latter are just artefacts of the PC-SAFT approach. Using PC-SAFT to correlate data for polybutadiene melts, we suggest that discrepancies in modelling the polymer density at ambient temperature and high pressure can be related to the presumably artificial liquid-liquid phase separation at lower temperatures. Thus, particular care is needed in engineering applications of the PC-SAFT theory that aims at predicting properties of macromolecular materials.     

The putative liquid-liquid transition is a liquid-solid transition in atomistic models of water

We use numerical simulation to examine the possibility of a reversible liquid-liquid transition in supercooled water and related systems. In particular, for two atomistic models of water, we have computed free energies as functions of multiple order parameters, where one is density and another distinguishes crystal from liquid. For a range of temperatures and pressures, separate free energy basins for liquid and crystal are found, conditions of phase coexistence between these phases are demonstrated, and time scales for equilibration are determined. We find that at no range of temperatures and pressures is there more than a single liquid basin, even at conditions where amorphous behavior is unstable with respect to the crystal. We find a similar result for a related model of silicon. This result excludes the possibility of the proposed liquid-liquid critical point for the models we have studied. Further, we argue that behaviors others have attributed to a liquid-liquid transition in water and related systems are in fact reflections of transitions between liquid and crystal.     

Adsorption, desorption, and clustering of H2O on Pt111

The adsorption, desorption, and clustering behavior of H2O on Pt111 has been investigated by specular He scattering. The data show that water adsorbed on a clean Pt111 surface undergoes a structural transition from a random distribution to clustered islands near 60 K. The initial helium scattering cross sections as a function of temperature are found to be insensitive to the incident H2O flux over a range of 0.005 monolayers (ML)/s-0.55 ML/s indicating that the clustering process is more complex than simple surface diffusion. The coarsening process of an initially random distribution of water deposited at 25 K is found to occur over a broad temperature range, 60相似文献   

A theoretical study of aqueous solvation of K comparing ab initio, polarizable, and fixed-charge models

The hydration of K(+) is studied using a hierarchy of theoretical approaches, including ab initio Born-Oppenheimer molecular dynamics and Car-Parrinello molecular dynamics, a polarizable force field model based on classical Drude oscillators, and a nonpolarizable fixed-charge potential based on the TIP3P water model. While models based more directly on quantum mechanics offer the possibility to account for complex electronic effects, polarizable and fixed-charges force fields allow for simulations of large systems and the calculation of thermodynamic observables with relatively modest computational costs. A particular emphasis is placed on investigating the sensitivity of the polarizable model to reproduce key aspects of aqueous K(+), such as the coordination structure, the bulk hydration free energy, and the self diffusion of K(+). It is generally found that, while the simple functional form of the polarizable Drude model imposes some restrictions on the range of properties that can simultaneously be fitted, the resulting hydration structure for aqueous K(+) agrees well with experiment and with more sophisticated computational models. A counterintuitive result, seen in Car-Parrinello molecular dynamics and in simulations with the Drude polarizable force field, is that the average induced molecular dipole of the water molecules within the first hydration shell around K(+) is slightly smaller than the corresponding value in the bulk. In final analysis, the perspective of K(+) hydration emerging from the various computational models is broadly consistent with experimental data, though at a finer level there remain a number of issues that should be resolved to further our ability in modeling ion hydration accurately.     

The properties of ion-water clusters. II. Solvation structures of Na+, Cl-, and H+ clusters as a function of temperature

Ion-water-cluster properties are investigated both through the multistate empirical valence bond potential and a polarizable model. Equilibrium properties of the ion-water clusters H+(H2O)100, Na+(H2O)100, Na+(H2O)20, and Cl-(H2O)17 in the temperature region 100-450 K are explored using a hybrid parallel basin-hopping and tempering algorithm. The effect of the solid-liquid phase transition in both caloric curves and structural distribution functions is investigated. It is found that sodium and chloride ions largely reside on the surface of water clusters below the cluster melting temperature but are solvated into the interior of the cluster above the melting temperature, while the solvated proton was found to have significant propensity to reside on or near the surface in both the liquid- and solid-state clusters.     

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.     

Structural and thermodynamic properties of different phases of supercooled liquid water

Computer simulation results are reported for a realistic polarizable potential model of water in the supercooled region. Three states, corresponding to the low density amorphous ice, high density amorphous ice, and very high density amorphous ice phases are chosen for the analyses. These states are located close to the liquid-liquid coexistence lines already shown to exist for the considered model. Thermodynamic and structural quantities are calculated, in order to characterize the properties of the three phases. The results point out the increasing relevance of the interstitial neighbors, which clearly appear in going from the low to the very high density amorphous phases. The interstitial neighbors are found to be, at the same time, also distant neighbors along the hydrogen bonded network of the molecules. The role of these interstitial neighbors has been discussed in connection with the interpretation of recent neutron scattering measurements. The structural properties of the systems are characterized by looking at the angular distribution of neighboring molecules, volume and face area distribution of the Voronoi polyhedra, and order parameters. The cumulative analysis of all the corresponding results confirms the assumption that a close similarity between the structural arrangement of molecules in the three explored amorphous phases and that of the ice polymorphs I(h), III, and VI exists.     

Thermodynamics of the liquid states of Langmuir monolayers

The liquid states and the liquid-liquid equilibrium of surfactant molecules forming an interphase between air and water have been considered using Monte Carlo computer simulations. Specifically, the expanded and compressed liquid phases observed for surfactant molecules were characterized as a function of pressure and temperature. Simple modified beadlike potentials were implemented in order to describe the interparticle forces between the hydrophobic and hydrophilic portions of surfactant molecules at the air/water interface. A simulation box was defined such that the monolayer was exposed to an externally applied lateral pressure in a modified isothermal-isobaric ensemble, whereas the water bath was modeled in a canonical ensemble. The simulation resembles the experimental setup used to measure lateral pressure (Pi) versus area isotherms obtained with Langmuir troughs. The applied lateral pressure-surface area phase diagram clearly showed the coexistence of the expanded and compressed liquid phases within certain temperature and pressure ranges. Distribution functions of distances and enthalpies for the monolayer were computed to clearly identify each liquid phase and the coexistence region.     

Temperature Effect on Liquid-Liquid Phase Equilibria of Aqueous Solutions of Fluorobenzene in Ethanol at the Pressure 101.2 kPa

Solubilities and liquid-liquid equilibrium (LLE) data for (water + ethanol + fluorobenzene) solutions are determined at T = (288.15, 298.15, 308.15) K and atmospheric pressure. The titration method is used to construct the binodal curves. All measured solubilities and liquid-liquid equilibrium data are reproduced quantitatively by empirical equations and the NRTL activity coefficient model. The correlated data are consistent to better than 0.3% in phase composition. The studied temperatures have a small effect on the phase equilibrium properties of the investigated systems. The influence of water concentration on the distribution coefficient of ethanol in the coexistent liquid phases and on the fluorobenzene selectivity is important. Practically, fluorobenzene may well be used to produce absolute alcohol.     

Liquid-liquid coexistence surface for lysozyme: role of salt type and salt concentration

The liquid-liquid phase separation curves for lysozyme in a salt solution are known to depend on salt type and salt concentration. For the case of monovalent cations, the cloud point temperature typically increases with increasing salt concentration, for fixed lysozyme concentration. For the case of divalent cations, however, a maximum in the cloud point temperature is observed that has been interpreted as being due to ion binding to the protein surface and subsequent water structuring. In this paper, we use a simple square well model due to Grigsby et al. (Biophys. Chem. 2001, 91, 231-243), whose well depth depends on salt type and salt concentration, to determine the phase coexistence surfaces from experimental data. The surfaces are shown as a function of temperature, salt concentration, and protein concentration for two typical salts, NaCl and MgCl2. These surfaces are calculated using the results of a single standard Monte Carlo simulation and a simple scaling argument and are in reasonably good agreement with known experimental results.