Water's unique anomalies are vital in various applications and biological processes, yet the molecular mechanisms behind these anomalies remain debated, particularly in the metastable liquid phase under supercooling and stretching conditions. Experimental challenges in these conditions have led to simulations suggesting a liquid-liquid phase transition between low-density and high-density water phases, culminating in a liquid-liquid critical point (LLCP). However, these simulations are limited by computational expense, small system sizes, and reliability of water models. Using the FS model, we improve accuracy in predicting water's density and response functions across a broad range of temperatures and pressures. The FS model avoid by design first-order phase transitions towards crystalline phases, allowing thorough exploration of the metastable phase diagram. We employ advanced numerical techniques to bypass dynamical slowing down and perform finite-size scaling on systems significantly larger than those used in previous analyses. Our study extrapolates thermodynamic behavior in the infinite-system limit, accurately demonstrating the existence of the LLCP in the 3D Ising universality class at TC = 186 +/- 4 K and PC = 174 +/- 14 MPa, following a liquid-liquid phase separation below 200 MPa. These predictions align with recent experimental data and more sophisticated models, highlighting that hydrogen bond cooperativity governs the LLCP and the origin of water anomalies. Moreover, we observe that the hydrogen bond network exhibits substantial cooperative fluctuations at scales larger than 10 nm, even at temperatures relevant to biopreservation. These findings have significant implications for fields such as nanotechnology and biophysics, offering new insights into water's behavior under varied conditions.

The molecular mechanisms of water's unique anomalies are still debated upon. Experimental challenges have led to simulations suggesting a liquid-liquid (LL) phase transition, culminating in the supercooled region's LL critical point (LLCP). Computational expense, small system sizes, and the reliability of water models often limit these simulations. We adopt the CVF model, which is reliable, transferable, scalable, and efficient across a wide range of temperatures and pressures around ambient conditions. By leveraging the timescale separation between fast hydrogen bonds and slow molecular coordinates, the model allows a thorough exploration of the metastable phase diagram of liquid water. Using advanced numerical techniques to bypass dynamical slowing down, we perform finite-size scaling on larger systems than those used in previous analyses. Our study extrapolates thermodynamic behavior in the infinite-system limit, demonstrating the existence of the LLCP in the 3D Ising universality class in the low-temperature, low-pressure side of the line of temperatures of maximum density, specifically at T-C = 186 +/- 4 K and P-C = 174 +/- 14 MPa, at the end of a liquid-liquid phase separation stretching up to similar to 200 MPa. These predictions align with recent experimental data and sophisticated models, highlighting that hydrogen bond cooperativity governs the LLCP and the origin of water anomalies. We also observe substantial cooperative fluctuations in the hydrogen bond network at scales larger than 10 nm, even at temperatures relevant to biopreservation. These findings have significant implications for nanotechnology and biophysics, providing new insights into water's behavior under varied conditions.

Water modeling is a challenging problem. Its anomalies are difficult to reproduce, promoting the proliferation of a large number of computational models, among which researchers select the most appropriate for the property they study. In this chapter, we introduce a coarse-grained model introduced by Franzese and Stanley (FS) that accounts for the many-body interactions of water. We review mean-field calculations and Monte Carlo simulations on water monolayers for a wide range of pressures and temperatures, including extreme conditions. The results show the presence of two dynamic crossovers and explain the origin of diffusion anomalies. Moreover, the model shows that all the different scenarios, proposed in the last decades as alternative explanations of the experimental anomalies of water, can be related by the fine-tuning of the many-body (cooperative) interaction. Once this parameter is set from the experiments, the FS model predicts a phase transition between two liquids with different densities and energies in the supercooled water region, ending in a liquid-liquid critical point. From this critical point stems a liquid-liquid Widom line, i.e., the locus of maxima of the water correlation length, that in the FS model can be directly calculated. The results are consistent with the extrapolations from experiments. Furthermore, they agree with those from atomistic models but make predictions over a much wider thermodynamic region, allowing for a better interpretation of the available experimental data. All these findings provide a coherent picture of the properties of water and confirm the validity of the FS model that has proved to be useful for large-scale simulations of biological systems.

Many-body interactions can play a relevant role in water properties. Here we study by Monte Carlo simulations a coarse-grained model for bulk water that includes many-body interactions associated to water cooperativity. The model is efficient and allows us to equilibrate water at extreme low temperatures in a wide range of pressures. We find the line of temperatures of maximum density at constant pressure and, at low temperature and high pressure, a phase transition between high-density liquid and low-density liquid phases. The liquid-liquid phase transition ends in a critical point. In the supercritical liquid region we find for each thermodynamic response function a locus of weak maxima at high temperature and a locus of strong maxima at low temperature, with both loci converging at the liquid-liquid critical point where they diverge. Comparing our results with previous works for the phase diagram of a many-body water monolayer, we observe that the weak maxima of the specific heat are less evident in bulk and appear only at negative pressures, while we find those of compressibility and thermal expansion coefficient also at positive pressures in bulk. However, the strong maxima of compressibility and thermal expansion coefficient are very sharp for the bulk case. Our results clarify fundamental properties of bulk water, possibly difficult to detect with atomistic models not accounting for many-body interactions.