Liquid water's isothermal compressibility and isobaric heat capacity, and the magnitude of its thermal expansion coefficient, increase sharply on cooling below the equilibrium freezing point. Many experimental, theoretical and computational studies have sought to understand the molecular origin and implications of this anomalous behaviour. Of the different theoretical scenarios put forward, one posits the existence of a first-order phase transition that involves two forms of liquid water and terminates at a critical point located at deeply supercooled conditions. Some experimental evidence is consistent with this hypothesis, but no definitive proof of a liquid-liquid transition in water has been obtained to date: rapid ice crystallization has so far prevented decisive measurements on deeply supercooled water, although this challenge has been overcome recently. Computer simulations are therefore crucial for exploring water's structure and behaviour in this regime, and have shown that some water models exhibit liquid-liquid transitions and others do not. However, recent work has argued that the liquid-liquid transition has been mistakenly interpreted, and is in fact a liquid-crystal transition in all atomistic models of water. Here we show, by studying the liquid-liquid transition in the ST2 model of water with the use of six advanced sampling methods to compute the free-energy surface, that two metastable liquid phases and a stable crystal phase exist at the same deeply supercooled thermodynamic condition, and that the transition between the two liquids satisfies the thermodynamic criteria of a first-order transition. We follow the rearrangement of water's coordination shell and topological ring structure along a thermodynamically reversible path from the low-density liquid to cubic ice. We also show that the system fluctuates freely between the two liquid phases rather than crystallizing. These findings provide unambiguous evidence for a liquid-liquid transition in the ST2 model of water, and point to the separation of time scales between crystallization and relaxation as being crucial for enabling it.

No abstract

We use the weighted histogram analysis method [S. Kumar, D. Bouzida, R. H. Swendsen, P. A. Kollman, and J. M. Rosenberg, J. Comput. Chem. 13, 1011 (1992)] to calculate the free energy surface of the ST2 model of water as a function of density and bond-orientational order. We perform our calculations at deeply supercooled conditions (T = 228.6 K, P = 2.2 kbar; T = 235 K, P = 2.2 kbar) and focus our attention on the region of bond-orientational order that is relevant to disordered phases. We find a first-order transition between a low-density liquid (LDL, ρ ≈ 0.9 g/cc) and a high-density liquid (HDL, ρ ≈ 1.15 g/cc), confirming our earlier sampling of the free energy surface of this model as a function of density [Y. Liu, A. Z. Panagiotopoulos, and P. G. Debenedetti, J. Chem. Phys. 131, 104508 (2009)]. We demonstrate the disappearance of the LDL basin at high pressure and of the HDL basin at low pressure, in agreement with independent simulations of the system's equation of state. Consistency between directly computed and reweighted free energies, as well as between free energy surfaces computed using different thermodynamic starting conditions, confirms proper equilibrium sampling. Diffusion and structural relaxation calculations demonstrate that equilibration of the LDL phase, which exhibits slow dynamics, is attained in the course of the simulations. Repeated flipping between the LDL and HDL phases in the course of long molecular dynamics runs provides further evidence of a phase transition. We use the Ewald summation with vacuum boundary conditions to calculate long-ranged Coulombic interactions and show that conducting boundary conditions lead to unphysical behavior at low temperatures.

Molecular modeling at the electronic and atomistic levels plays an important and complementary role to experimental studies of confinement effects. Theory and atomistic simulation can provide fundamental understanding, determine the limits of well known macroscopic laws such as Kelvin's equation, provide predictions for systems that are difficult to study via experiment (e.g. adsorption of highly toxic gases), and can be used to gain detailed molecular level information that may not be accessible in the laboratory (e.g. the local structure and composition of confined phases). We describe the most important and useful methods that are based firmly on quantum mechanics and statistical mechanics, including ab intio and classical density functional theories, and Monte Carlo and molecular dynamics simulation. We discuss their strengths and limitations. We then describe examples of applications of these methods to adsorption and equilibrium properties, including testing the Kelvin equation, determination of pore size distributions and capillary phenomena. Applications to self and transport diffusion, including single-file and anomalous diffusion, and viscous flow in nanoporous materials are described. The use of these methods to understand confinement effects on chemical reactions in heterogeneous media is treated, including effects on reaction equilibria, rates and mechanism. Finally we discuss the current status of molecular modeling in this area, and the outlook and future research needs for the next few years. The treatment is suitable for the general technical reader.

Thermodynamic properties of the ST2 model for supercooled liquid water exhibit anomalies similar to those observed in real water. A possible explanation of these anomalies is the existence of a metastable, liquid-liquid transition terminated by a critical point. This phenomenon, whose possible existence in real water is the subject of much current experimental work, has been unambiguously demonstrated for this particular model by most recent simulations. In this work, we reproduce the anomalies of two versions of the ST2 model with an equation of state describing water as a non-ideal "mixture" of two different types of local molecular order. We show that the liquid-liquid transition in the ST2 water is energy-driven. This is in contrast to another popular model, mW, in which nonideality in mixing of two alternative local molecular orders is entropy-driven, and is not sufficiently strong to induce a liquid-liquid transition.