Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liquid phases---low-density liquid (LDL) and high-density liquid (HDL)---deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stability of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature and constant number of molecules N for N<=729 and times up to 1000 ns. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N=512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs tiny crystallites grow and then melt within 1000 ns. Only for N<=343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystallization after crystallites reach an estimated critical size of about 70+/-10 molecules.

Pressure Effects in Supercooled Water: Comparison between a 2D Model of Water and Experiments for Surface Water on a Protein

Experiments in bulk water confirm the existence of two local arrangements of water molecules with different densities, but, because of inevitable freezing at low temperature $T$, can not ascertain whether the two arrangements separate in two phases. To avoid the freezing, new experiments measure the dynamics of water at low $T$ on the surface of proteins, finding a crossover from a non-Arrhenius regime at high $T$ to a regime that is approximately Arrhenius at low $T$. Motivated by these experiments, Kumar et al. [Phys. Rev. Lett. 100, 105701 (2008)] investigated, by Monte Carlo simulations and mean field calculations, the relation of the dynamic crossover with the coexistence of two liquid phases in a cell model for water and predict that: (i) the dynamic crossover is isochronic, i.e. the value of the crossover time $tau_{rm L}$ is approximately independent of pressure $P$; (ii) the Arrhenius activation energy $E_{rm A}(P)$ of the low-$T$ regime decreases upon increasing $P$; (iii) the temperature $T^*(P)$ at which $tau$ reaches a fixed macroscopic time $tau^*geq tau_{rm L}$ decreases upon increasing $P$; in particular, this is true also for the crossover temperature $T_{rm L}(P)$ at which $tau=tau_{rm L}$. Here, we compare these predictions with recent quasi elastic neutron scattering (QENS) experiments performed by X.-Q. Chu {it et al.} on hydrated proteins at different values of $P$. We find that the experiments are consistent with these three predictions.