Advances in the study of supercooled water

Breakdown of the Stokes-Einstein relation in Stillinger-Weber silicon

We investigate the dynamical properties of liquid and supercooled liquid silicon, modeled using the Stillinger-Weber potential, to examine the validity of the Stokes-Einstein (SE) relation. Toward this end, we examine the relationship among various dynamical quantities, including (i) the macroscopic transport coefficients-self-diffusion coefficient D and viscosity η, (ii) relaxation time τα, and (iii) lengthscale dependent relaxation times τα(q) over a broad range of temperature T, pressure P, and density ρ covering both equilibrium and metastable liquid state points in the phase diagram. Our study shows a weak breakdown in the SE relation involving D and η, and the loci of the breakdown of the SE relation (SEB) are found in the high T liquid phase. The τα, when used as a proxy to η, shows a distinct breakdown in the SE relation, whose loci are found in the supercooled liquid phase. Interestingly, certain parts of the phase diagram show that the loci of onset of slow dynamics lie below the loci of SEB, suggesting a regime that exhibits Arrhenius but non-Fickian behavior. Computation of τα(q) enables us to extract the lengthscale associated with the Fickian to non-Fickian behavior using which we show that the breakdown of the SE relation occurs only below a specific lengthscale at a given temperature. Furthermore, we also compare the SEB loci with other features of the phase behavior, including the loci of compressibility maximum, density maximum, and diffusivity maximum.

Toward mitigating the impact of non-bulk defects on describing water structure in salt aqueous solutions: Characterizing solution density with a network-based structural indicator

Structural indicators, also known as structural descriptors, including order parameters, have been proposed to quantify the structural properties of water to account for its anomalous behaviors. However, these indicators, mainly designed for bulk water, are not naturally transferrable to the vicinity of ions due to disruptions in the immediate neighboring space and a resulting loss of feature completeness. To address these non-bulk defects, we introduced a structural indicator that draws on the concept of clique number from graph theory and the criterion in agglomerative clustering, denoted as the average cluster number. This structural indicator aims to discern intrinsic structural characteristics within the water molecules regardless of the ions occupying the neighboring space, without requiring additional corrections. From molecular dynamics simulation results for neat water and salt aqueous solutions utilizing the TIP4P/2005 water model and the Madrid-2019 force field, we characterized the variations in densities with temperature using this network-based indicator, thereby demonstrating its practical utility. The findings suggest that at lower temperatures, the addition of ions disrupts the intrinsic structure of water molecules, with this effect diminishing as the temperature rises. Cations with larger charge density tend to induce stronger disruptions. This study highlights the importance of mitigating the impact of non-bulk defects before applying the indicators to analyze water's intrinsic structural properties in solutions. By doing so, the relationship between changes in water structure and solution behaviors can be more accurately assessed.

An Alternative Hypothesis on Enhanced Deep Supercooling of Water: Nucleator Inhibition via Bicarbonate Adsorption

Supercooling allows for retarding water crystallization toward negative Celsius temperatures. Previous findings of CO2 molecules shifting into bicarbonate species upon freezing, the latter which naturally adsorb on hydrophobic interfaces, are put in perspective here to interpret earlier published data. Since it has been shown that ice nucleation is unlikely on negatively charged surfaces, I propose that bicarbonates adsorb on most solid particles present in water that act as nucleators, thus retarding freezing and enhancing supercooling. This hypothesis can now explain the deep supercooling observed for sealed and boiled water samples and oil-topped water samples, promoting both more bicarbonate generation and adsorption. Such an explanation opens new directions for access to cryopreservation.

Study of Protein Hydration Water with the V4S Structural Index: Focus on Binding Site Description

V4S, a new structural indicator for water specially designed to be suitable for hydration and nanoconfined contexts, has been recently introduced and preliminarily applied for water in contact with self-assembled monolayers and graphene-like systems. This index enabled an accurate detection of defective high local density water molecules (called HDA-like given their structural resemblance with the high-density amorphous ice, HDA). In the present work, we shall apply this new metric to characterize protein hydration water with particular interest in protein binding sites. As a first result, we shall find that protein hydration water has a higher concentration of HDA-like molecular arrangements compared to the bulk. Significantly, we shall show that the concentration of HDA-like molecules sharply decreases beyond the first hydration layer. Finally, we shall also reveal a highly nonuniform spatial distribution of the V4S values for the first hydration shell on the protein surface, where the higher hydrophobicity inherent to the ligand binding site will be evident from an enrichment in HDA-like molecules as compared to the population exhibited by the global protein surface.

Medium-density amorphous ice unveils shear rate as a new dimension in water's phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice Ih at 77 K [Rosu-Finsen et al., Science 379, 474-478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water's phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. We find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice Ih, LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 106 s-1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~108 s-1 have density, enthalpy, and structure indistinguishable from those of water "instantaneously" quenched from room temperature to 77 K over 10 ps, making them good approximants for the "true glass" of ambient liquid water.

Role of hydrogen-bond coordination defects in the structural relaxation of supercooled water

In this work, we shall study the role of threefold and fivefold coordination defects in the structure and dynamics of the hydrogen bond network of liquid water, with special emphasis on the glassy regime. A significant defect clusterization propensity will be made evident, with a prevalence of mixed pairs, that is, threefold- and fivefold-coordinated defects being first neighbors of each other. This structural analysis will enable us to determine the existence of defective and nondefective regions compatible with the high local density and low local density molecular states of liquid water, respectively. Hydrogen bond coordination defects will also be shown to promote water's structural relaxation, with the undercoordinated ones playing a main role in driving glassy relaxation dynamics. Moreover, we shall show that the three-foldcoordinated molecules together with their first neighbors present at the initial configuration act as markers of the dynamical heterogeneities that would emerge at later times commensurate with the structural relaxation of the supercooled system.

Phase behavior of metastable water from large-scale simulations of a quantitatively accurate model near ambient conditions: The liquid-liquid critical point

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 TC = 186 ± 4 K and PC = 174 ± 14 MPa, at the end of a liquid-liquid phase separation stretching up 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.

A structural determinant of the behavior of water at hydration and nanoconfinement conditions

The molecular nature of the phases that conform the two-liquid scenario is elucidated in this work in the light of a molecular principle governing water structuring, which unveils the relevance of the contraction and reorientation of the second molecular shell to allow for the existence of coordination defects in water's hydrogen bond network. In turn, such principle is shown to also determine the behavior of hydration and nanoconfined water while enabling to define conditions for wettability (quantifying hydrophobicity and predicting drying transitions), thus opening the possibility to unravel the active role of water in central fields of research.

Estimating metastable thermodynamic properties by isochoric extrapolation from stable states

The description of metastable fluids, those in local but not global equilibrium, remains an important problem of thermodynamics, and it is crucial for many industrial applications and all first order phase transitions. One way to estimate their properties is by extrapolation from nearby stable states. This is often done isothermally, in terms of a virial expansion for gases or a Taylor expansion in density for liquids. This work presents evidence that an isochoric expansion of pressure at a given temperature is superior to an isothermal density expansion. Two different isochoric extrapolation strategies are evaluated, one best suited for vapors and one for liquids. Both are exact for important model systems, including the van der Waals equation of state. Moreover, we present a simple method to evaluate all the coefficients of the isochoric expansion directly from a simulation in the canonical ensemble. Using only the properties of stable states, the isochoric extrapolation methods reproduce simulation results with Lennard-Jones potentials, mostly within their uncertainties. The isochoric extrapolation methods are able to predict deeply metastable pressures accurately even from temperatures well above the critical. Isochoric extrapolation also predicts a mechanical stability limit, i.e., the thermodynamic spinodal. For water, the liquid spinodal pressure is predicted to be monotonically decreasing with decreasing temperature, in contrast to the re-entrant behavior predicted by the direct extension of the reference equation of state.

Water at the nanoscale: From filling or dewetting hydrophobic pores and carbon nanotubes to "sliding" on graphene

In this work, we study the effect of nanoconfinement on the hydration properties of model hydrophobic pores and carbon nanotubes, determining their wetting propensity and the conditions for geometrically induced dehydration. By employing a recently introduced water structural index, we aim at two main goals: (1) to accurately quantify the local hydrophobicity and predict the drying transitions in such systems, and (2) to provide a molecular rationalization of the wetting process. In this sense, we will further discuss the number and strength of the interactions required by the water molecules to promote wetting. In the case of graphene-like surfaces, an explanation for their unexpectedly significant hydrophilicity will also be provided. On the one hand, the structural index will show that the net attraction to the dense carbon network that a water molecule experiences through several simultaneous weak interactions is sufficient to give rise to hydrophilic behavior. On the other hand, we will show that an additional effect is also at play: the hydrating water molecule is retained on the surface by a smooth exchange of such simultaneous weak interactions, as if "sliding" on graphene.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

Unveiling a common phase transition pathway of high-density amorphous ices through time-resolved x-ray scattering

Here, we investigate the hypothesis that despite the existence of at least two high-density amorphous ices, only one high-density liquid state exists in water. We prepared a very-high-density amorphous ice (VHDA) sample and rapidly increased its temperature to around 205 ± 10 K using laser-induced isochoric heating. This temperature falls within the so-called "no-man's land" well above the glass-liquid transition, wherein the IR laser pulse creates a metastable liquid state. Subsequently, this high-density liquid (HDL) state of water decompresses over time, and we examined the time-dependent structural changes using short x-ray pulses from a free electron laser. We observed a liquid-liquid transition to low-density liquid water (LDL) over time scales ranging from 20 ns to 3 μs, consistent with previous experimental results using expanded high-density amorphous ice (eHDA) as the initial state. In addition, the resulting LDL derived both from VHDA and eHDA displays similar density and degree of inhomogeneity. Our observation supports the idea that regardless of the initial annealing states of the high-density amorphous ices, the same HDL and final LDL states are reached at temperatures around 205 K.

Structure and slow dynamics of protein hydration water with cryopreserving DMSO and trehalose upon cooling

We study, through molecular dynamics simulations, three aqueous solutions with one lysozyme protein and three different concentrations of trehalose and dimethyl sulfoxide (DMSO). We analyze the structural and dynamical properties of the protein hydration water upon cooling. We find that trehalose plays a major role in modifying the structure of the network of HBs between water molecules in the hydration layer of the protein. The dynamics of hydration water presents, in addition to the α-relaxation, typical of glass formers, a slower long-time relaxation process, which greatly slows down the dynamics of water, particularly in the systems with trehalose, where it becomes dominant at low temperatures. In all the solutions, we observe, from the behavior of the α-relaxation times, a shift of the Mode Coupling Theory crossover temperature and the fragile-to-strong crossover temperature toward higher values with respect to bulk water. We also observe a strong-to-strong crossover from the temperature behavior of the long-relaxation times. In the aqueous solution with only DMSO, the transition shifts to a lower temperature than in the case with only lysozyme reported in the literature. We observe that the addition of trehalose to the mixture has the opposite effect of restoring the original location of the strong-to-strong crossover. In all the solutions analyzed in this work, the observed temperature of the protein dynamical transition is slightly shifted at lower temperatures than that of the strong-to-strong crossover, but their relative order is the same, showing a correlation between the motion of the protein and that of the hydration water.

Local and global expansivity in water

The supra-molecular structure of a liquid is strongly connected to its dynamics, which in turn control macroscopic properties such as viscosity. Consequently, detailed knowledge about how this structure changes with temperature is essential to understand the thermal evolution of the dynamics ranging from the liquid to the glass. Here, we combine infrared spectroscopy (IR) measurements of the hydrogen (H) bond stretching vibration of water with molecular dynamics simulations and employ a quantitative analysis to extract the inter-molecular H-bond length in a wide temperature range of the liquid. The extracted expansivity of this H-bond differs strongly from that of the average nearest neighbor distance of oxygen atoms obtained through a common conversion of mass density. However, both properties can be connected through a simple model based on a random loose packing of spheres with a variable coordination number, which demonstrates the relevance of supra-molecular arrangement. Furthermore, the exclusion of the expansivity of the inter-molecular H-bonds reveals that the most compact molecular arrangement is formed in the range of ∼316-331K (i.e., above the density maximum) close to the temperature of several pressure-related anomalies, which indicates a characteristic point in the supra-molecular arrangement. These results confirm our earlier approach to deduce inter-molecular H-bond lengths via IR in polyalcohols [Gabriel et al. J. Chem. Phys. 154, 024503 (2021)] quantitatively and open a new alley to investigate the role of inter-molecular expansion as a precursor of molecular fluctuations on a bond-specific level.

Unraveling thermodynamic anomalies of water: A molecular simulation approach to probe the two-state theory with atomistic and coarse-grained water models

Thermodynamic and dynamic anomalies of water play a crucial role in supporting life on our planet. The two-state theory attributes these anomalies to a dynamic equilibrium between locally favored tetrahedral structures (LFTSs) and disordered normal liquid structures. This theory provides a straightforward, phenomenological explanation for water's unique thermodynamic and dynamic characteristics. To validate this two-state feature, it is critical to unequivocally identify these structural motifs in a dynamically fluctuating disordered liquid. In this study, we employ a recently introduced structural parameter (θavg) that characterizes the local angular order within the first coordination shell to identify these LFTSs through molecular dynamics simulations. We employ both realistic water models with a liquid-liquid critical point (LLCP) and a coarse-grained water model without an LLCP to study water's anomalies in low-pressure regions below 2 kbar. The two-state theory consistently describes water's thermodynamic anomalies in these models, both with and without an LLCP. This suggests that the anomalies predominantly result from the two-state features rather than criticality, particularly within experimentally accessible temperature-pressure regions.

A water structure indicator suitable for generic contexts: Two-liquid behavior at hydration and nanoconfinement conditions and a molecular approach to hydrophobicity and wetting

In a recent work, we have briefly introduced a new structural index for water that, unlike previous indicators, was devised specifically for generic contexts beyond bulk conditions, making it suitable for hydration and nanoconfinement settings. In this work, we shall study this metric in detail, demonstrating its ability to reveal the existence of a fine-tuned interplay between the local structure and energetics in liquid water. This molecular principle enables the establishment of an extended hydrogen bond network, while simultaneously allowing for the existence of network defects by compensating for uncoordinated sites. By studying different water models and different temperatures encompassing both the normal liquid and the supercooled regime, this molecular mechanism will be shown to underlie the two-state behavior of bulk water. In addition, by studying functionalized self-assembled monolayers and diverse graphene-like surfaces, we shall show that this principle is also operative at hydration and nanoconfinement conditions, thus generalizing the validity of the two-liquid scenario of water to these contexts. This approach will allow us to define conditions for wettability, providing an accurate measure of hydrophobicity and a reliable predictor of filling and drying transitions. Hence, it might open the possibility of elucidating the active role of water in the broad fields of biophysics and materials science. As a preliminary step, we shall study the hydration structure and hydrophilicity of graphene-like systems (parallel graphene sheets and carbon nanotubes) as a function of the confinement dimensionality.

Tuning the low-temperature phase behavior of aqueous ionic liquids

Water's anomalous behavior is often explained using a two-liquid model, where two types of water, high-density liquid (HDL) and low-density liquid (LDL), can be separated via a liquid-liquid phase transition (LLPT) at low temperature. Mixtures of water and the ionic liquid hydrazinium trifluoroacetate were suggested to also show an LLPT but with the advantage that there is no rapid ice crystallization hampering its observation. It remains controversial whether these solutions exhibit an LLPT or are instead associated with complex phase separation phenomena. We here show detailed low-temperature calorimetry and diffraction experiments on aqueous solutions containing hydrazinium trifluoroacetate and other similar ionic liquids, all at a solute mole fraction of x = 0.175. Hydrazinium trifluoroacetate, ammonium trifluoroacetate, ethylammonium trifluoroacetate and hydrazinium pentafluoropropionate all boast exothermic transitions unrelated to crystallization as well as remarkable structural changes upon cooling into the glassy state. We propose a model inspired by micelle formation and decomposition in surfactant solutions, which is complemented by MD simulations and allows rationalizing the rich phase behavior of our mixtures during cooling. The fundamental aspect of the model is the hydrophobic nature of fluorinated anions that enables aggregation, which is reversed upon cooling and culminates in the remarkable exothermic first-order transition observed at low temperature. That is, we assign the first-order transition not to an LLPT but to phase-separations similar to the ones when falling below the Krafft temperature. All other solutions merely show simple vitrification behavior. Still, they exhibit distinct differences in liquid fragility, which is decreased continuously with decreasing hydrophobicity of the anions. This might enable the systematic tuning of ionic liquids with the goal of designing aqueous solutions of specific fragility.

Deciphering Density Fluctuations in the Hydration Water of Brownian Nanoparticles via Upconversion Thermometry

We investigate the intricate relationship among temperature, pH, and Brownian velocity in a range of differently sized upconversion nanoparticles (UCNPs) dispersed in water. These UCNPs, acting as nanorulers, offer insights into assessing the relative proportion of high-density and low-density liquid in the surrounding hydration water. The study reveals a size-dependent reduction in the onset temperature of liquid-water fluctuations, indicating an augmented presence of high-density liquid domains at the nanoparticle surfaces. The observed upper-temperature threshold is consistent with a hypothetical phase diagram of water, validating the two-state model. Moreover, an increase in pH disrupts the organization of water molecules, similar to external pressure effects, allowing simulation of the effects of temperature and pressure on hydrogen bonding networks. The findings underscore the significance of the surface of suspended nanoparticles for understanding high- to low-density liquid fluctuations and water behavior at charged interfaces.

The structure of water: A historical perspective

Attempts to understand the molecular structure of water were first made well over a century ago. Looking back at the various attempts, it is illuminating to see how these were conditioned by the state of knowledge of chemistry and physics at the time and the experimental and theoretical tools then available. Progress in the intervening years has been facilitated by not only conceptual and theoretical advances in physics and chemistry but also the development of experimental techniques and instrumentation. Exploitation of powerful computational methods in interpreting what at first sight may seem impenetrable experimental data has led us to the consistent and detailed picture we have today of not only the structure of liquid water itself and how it changes with temperature and pressure but also its interactions with other molecules, in particular those relevant to water's role in important chemical and biological processes. Much remains to be done in the latter areas, but the experimental and computational techniques that now enable us to do what might reasonably be termed "liquid state crystallography" have opened the door to make possible further advances. Consequently, we now have the tools to explore further the role of water in those processes that underpin life itself-the very prospect that inspired Bernal to develop his ideas on the structure of liquids in general and of water in particular.

A DPD model of soft spheres with waterlike anomalies and poly(a)morphism

Core-softened approaches have been employed to understand the behavior of a large variety of systems in soft condensed matter, from biological molecules to colloidal crystals, glassy phases, and water-like anomalies. At the same time, dissipative particle dynamics (DPD) is a powerful tool suitable for studying larger length and time scales. In this sense, we propose a simple model of soft molecules that exhibits a wide range of interesting phenomena: polyamorphism, with three amorphous phases, polymorphysm, including a recently found gyroid phase and a cubic diamond structure, reentrant liquid phase, and density, diffusion, and structural water-like anomalies. Each molecule is constituted by two collapsing beads, representing a harder central core and a softer corona. This induces a competition between distinct conformations that leads to their unique behavior. This provides a basis for the development of more accurate water-like DPD models that can then be parameterized for specific systems and even used to model and understand the self-assembly of colloidal crystals.

Turning an energy-based defect detector into a multi-molecule structural indicator for water

Recent studies have provided conclusive evidence for the existence of a liquid-liquid critical point in numerical models of water. Such a scenario implies the competition between two local molecular arrangements of different densities: a high-density liquid (HDL) and a low-density liquid (LDL). Within this context, the development of accurate structural indicators to properly characterize the two interconverting local structures is demanded. In a previous study, we introduced a reliable energy-based structural descriptor that properly discriminates water molecules into tetrahedrally arranged molecules (T molecules) and distorted molecules (D molecules). The latter constitute defects in terms of hydrogen bond (HB) coordination and have been shown to represent a minority component, even at high temperatures above the melting point. In addition, the D molecules tend to form high-quality HBs with three T molecules and to be surrounded by T and D molecules at further distances. Thus, it became evident that, while the LDL state might consist of a virtually pure T state, the HDL state would comprise mixed molecular arrangements including the D molecules. Such a need to abandon the single-molecule description requires the investigation of the degree of structural information to be incorporated in order to build an appropriate multi-molecule indicator. Hence, in this work, we shall study the effect of the local structural constraints on the water molecules in order to discriminate the different molecular arrangements into two disjoint classes. This will enable us to build a multi-molecule structural indicator for water whose performance will then be investigated within the water's supercooled regime.

Phase Diagram of Aqueous Solutions of LiCl: a Study of Concentration Effects on the Anomalies of Water

We perform molecular dynamics simulations in order to study thermodynamics and the structure of supercooled aqueous solutions of lithium chloride (LiCl) at concentrations c = 0.678 and 2.034 mol/kg. We model the solvent using the TIP4P/2005 potential and the ions using the Madrid-2019 force field, a force field particularly suited for studying this solution. We find that, for c = 0.678 mol/kg, the behavior of the equation of state, studied in the P-T plane, indicates the presence of a liquid-liquid phase transition, similar to what was previously found for bulk water. We estimate the position of the liquid-liquid critical point to be at T_c ≈ 174 K, P_c ≈ 1775 bar, and ρ_c ≈ 1.065 g/cm³. When the concentration is tripled to c = 2.034 mol/kg, no critical point is observed, indicating its possible disappearance at this concentration. We also study the water-water and water-ions structure in the two solutions, and we find that at the concentrations examined the effect of ions on the water-water structure is not strong, and all the features found in bulk water are preserved. We also calculate the hydration number of the Li and Cl ions, and in line with experiments, we find the value of 4 for Li⁺ and between 5.5 and 6 for Cl^-, confirming the good performances of the Madrid-2019 force field.

Are Neural Network Potentials Trained on Liquid States Transferable to Crystal Nucleation? A Test on Ice Nucleation in the mW Water Model

Neural network potentials (NNPs) are increasingly being used to study processes that happen on long time scales. A typical example is crystal nucleation, which rate is controlled by the occurrence of a rare fluctuation, i.e., the appearance of the critical nucleus. Because the properties of this nucleus are far from those of the bulk crystal, it is yet unclear whether NN potentials trained on equilibrium liquid states can accurately describe nucleation processes. So far, nucleation studies on NNPs have been limited to ab initio models whose nucleation properties are unknown, preventing an accurate comparison. Here we train a NN potential on the mW model of water─a classical three-body potential whose nucleation time scale is accessible in standard simulations. We show that a NNP trained only on a small number of liquid state points can reproduce with great accuracy the nucleation rates and free energy barriers of the original model, computed from both spontaneous and biased trajectories, strongly supporting the use of NNPs to study nucleation events.

Size-Pore-Dependent Methanol Sequestration from Water-Methanol Mixtures by an Embedded Graphene Slit

The separation of liquid mixture components is relevant to many applications-ranging from water purification to biofuel production-and is a growing concern related to the UN Sustainable Development Goals (SDGs), such as "Clean water and Sanitation" and "Affordable and clean energy". One promising technique is using graphene slit-pores as filters, or sponges, because the confinement potentially affects the properties of the mixture components in different ways, favoring their separation. However, no systematic study has shown how the size of a pore changes the thermodynamics of the surrounding mixture. Here, we focus on water-methanol mixtures and explore, using Molecular Dynamics simulations, the effects of a graphene pore, with size ranging from 6.5 to 13 Å, for three compositions: pure water, 90%-10%, and 75%-25% water-methanol. We show that tuning the pore size can change the mixture pressure, density and composition in bulk due to the size-dependent methanol sequestration within the pore. Our results can help in optimizing the graphene pore size for filtering applications.

Correlations between defect propensity and dynamical heterogeneities in supercooled water

A salient feature of supercooled liquids consists in the dramatic dynamical slowdown they undergo as temperature decreases while no significant structural change is evident. These systems also present dynamical heterogeneities (DH): certain molecules, spatially arranged in clusters, relax various orders of magnitude faster than the others. However, again, no static quantity (such as structural or energetic measures) shows strong direct correlations with such fast-moving molecules. In turn, the dynamic propensity approach, an indirect measure that quantifies the tendency of the molecules to move in a given structural configuration, has revealed that dynamical constraints, indeed, originate from the initial structure. Nevertheless, this approach is not able to elicit which structural quantity is, in fact, responsible for such a behavior. In an effort to remove dynamics from its definition in favor of a static quantity, an energy-based propensity has also been developed for supercooled water, but it could only find positive correlations between the lowest-energy and the least-mobile molecules, while no correlations could be found for those more relevant mobile molecules involved in the DH clusters responsible for the system's structural relaxation. Thus, in this work, we shall define a defect propensity measure based on a recently introduced structural index that accurately characterizes water structural defects. We shall show that this defect propensity measure provides positive correlations with dynamic propensity, being also able to account for the fast-moving molecules responsible for the structural relaxation. Moreover, time dependent correlations will show that defect propensity represents an appropriate early-time predictor of the long-time dynamical heterogeneity.

Absorption of pressurized methane in normal and supercooled p-xylene revealed via high-resolution neutron imaging

Supercooling of liquids leads to peculiarities which are scarcely studied under high-pressure conditions. Here, we report the surface tension, solubility, diffusivity, and partial molar volume for normal and supercooled liquid solutions of methane with p-xylene. Liquid bodies of perdeuterated p-xylene (p-C₈D₁₀), and, for comparison, o-xylene (o-C₈D₁₀), were exposed to pressurized methane (CH₄, up to 101 bar) at temperatures ranging 7.0-30.0 °C and observed at high spatial resolution (pixel size 20.3 μm) using a non-tactile neutron imaging method. Supercooling led to the increase of diffusivity and partial molar volume of methane. Solubility and surface tension were insensitive to supercooling, the latter substantially depended on methane pressure. Overall, neutron imaging enabled to reveal and quantify multiple phenomena occurring in supercooled liquid p-xylene solutions of methane under pressures relevant to the freeze-out in the production of liquefied natural gas.

Self-Diffusion in Confined Water: A Comparison between the Dynamics of Supercooled Water in Hydrophobic Carbon Nanotubes and Hydrophilic Porous Silica

Confined liquids are model systems for the study of the metastable supercooled state, especially for bulk water, in which the onset of crystallization below 230 K hinders the application of experimental techniques. Nevertheless, in addition to suppressing crystallization, confinement at the nanoscale drastically alters the properties of water. Evidently, the behavior of confined water depends critically on the nature of the confining environment and the interactions of confined water molecules with the confining matrix. A comparative study of the dynamics of water under hydrophobic and hydrophilic confinement could therefore help to clarify the underlying interactions. As we demonstrate in this work using a few representative results from the relevant literature, the accurate assessment of the translational mobility of water molecules, especially in the supercooled state, can unmistakably distinguish between the hydrophilic and hydrophobic nature of the confining environments. Among the numerous experimental methods currently available, we selected nuclear magnetic resonance (NMR) in a field gradient, which directly measures the macroscopic translational self-diffusion coefficient, and quasi-elastic neutron scattering (QENS), which can determine the microscopic translational dynamics of the water molecules. Dielectric relaxation, which probes the re-orientational degrees of freedom, are also discussed.

Manifestations of the structural origin of supercooled water’s anomalies in the heterogeneous relaxation on the potential energy landscape

Liquid water is well-known for its intriguing thermodynamic anomalies in the supercooled state. The phenomenological two-state models—based on the assumption of the existence of two types of competing local states (or, structures) in liquid water—have been extremely successful in describing water’s thermodynamic anomalies. However, the precise structural features of these competing local states in liquid water still remain elusive. Here, we have employed a predefined structural order parameter-free approach to unambiguously identify two types of competing local states—entropically and energetically favored—with significantly different structural and energetic features in the TIP4P/2005 liquid water. This identification is based on the heterogeneous structural relaxation of the system in the potential energy landscape (PEL) during the steepest-descent energy minimization. This heterogeneous relaxation is characterized using order parameters inspired by the spin-glass transition in frustrated magnetic systems. We have further established a direct relationship between the population fluctuation of the two states and the anomalous behavior of the heat capacity in supercooled water. The composition-dependent spatial distribution of the entropically favored local states shows an interesting crossover from a spanning network-like single cluster to the spatially delocalized clusters in the close vicinity of the Widom line. Additionally, this study establishes a direct relationship between the topographic features of the PEL and the water’s thermodynamic anomalies in the supercooled state and provides alternate markers (in addition to the locus of maxima of thermodynamic response functions) for the Widom line in the phase plane.

Water's Unusual Thermodynamics in the Realm of Physical Chemistry

While it is known since the early work by Edsall, Frank and Evans, Kauzmann, and others that the thermodynamics of solvation of nonpolar solutes in water is unusual and has implications for the thermodynamics of protein folding, only recently have its connections with the unusual temperature dependence of the density of solvent water been illuminated. Such density behavior is, in turn, one of the manifestations of a nonstandard thermodynamic pattern contemplating a second, liquid-liquid critical point at conditions of temperature and pressure at which water exists as a deeply supercooled liquid. Recent experimental and computational work unambiguously points toward the existence of such a critical point, thereby providing concrete answers to the questions posed by the 1976 pioneering experiments by Speedy and Angell and the associated "liquid-liquid transition hypothesis" posited in 1992 by Stanley and co-workers. Challenges of this phenomenology to the branch of Statistical Mechanics remain.

Homogeneous ice nucleation rates for mW and TIP4P/ICE models through Lattice Mold calculations

Water freezing is the most common liquid-to-crystal phase transition on Earth, however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240K. We employ two different water models, mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic non-polarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih, Ic and a stacking mixture of ice Ih/Ic; reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner which contributes to benchmarking the freezing behaviour of two popular water models.

Salt Enrichment and Dynamics in the Interface of Supercooled Aqueous Droplets

The interconversion reaction of NaCl between the contact-ion pair (CIP) and the solvent-separated ion pair (SSIP) as well as the free-ion state in cold droplets has not yet been investigated. We report direct computational evidence that the lower is the temperature, the closer to the surface the ion interconversion reaction takes place. In supercooled droplets the enrichment of the subsurface in salt becomes more evident. The stability of the SSIP relative to the CIP increases as the ion-pairing is transferred toward the droplet's outer layers. In the free-ion state, where the ions diffuse independently in the solution, the number density of Cl^- shows a broad maximum in the interior in addition to the well-known maximum in the surface. In the study of the reaction dynamics, we find a weak coupling between the interionic NaCl distance reaction coordinate and the solvent degrees of freedom, which contrasts with the diffusive crossing of the free energy barrier found in bulk solution modeling. The H₂O self-diffusion coefficient is found to be at least an order of magnitude larger than that in the bulk solution. We propose to exploit the enhanced surface ion concentration at low temperature to eliminate salts from droplets in native mass spectrometry ionization methods.

Maximum in density of electrolyte solutions: Learning about ion-water interactions and testing the Madrid-2019 force field

In this work, we studied the effect of Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation. All the salts studied in this work shift the TMD of the solution to lower temperatures and flatten the density vs temperature curves (when compared to pure water) with increasing salt concentration. The shift in the TMD depends strongly on the nature of the electrolyte. In order to explore this dependence, we have evaluated the contribution of each ion to the shift in the TMD concluding that Na⁺, Ca²⁺, and SO₄ ^2- seem to induce the largest changes among the studied ions. The volume of the system has been analyzed for salts with the same anion and different cations. These curves provide insight into the effect of different ions upon the structure of water. We claim that the TMD of electrolyte solutions entails interesting physics regarding ion-water and water-water interactions and should, therefore, be considered as a test property when developing force fields for electrolytes. This matter has been rather unnoticed for almost a century now and we believe it is time to revisit it.

Oxygen NMR of high-density and low-density amorphous ice

Using oxygen-17 as a nuclear probe, spin relaxometry was applied to study the high-density and low-density states of amorphous ice, covering temperatures below and somewhat above their glass transitions. These findings are put in perspective with results from deuteron nuclear magnetic resonance and with calculations based on dielectrically detected correlation times. This comparison reveals the presence of a wide distribution of correlation times. Furthermore, oxygen-17 central-transition echo spectra were recorded for wide ranges of temperature and pulse spacing. The spectra cannot be described by a single set of quadrupolar parameters, suggesting a distribution of H–O–H opening angles that is broader for high-density than for low-density amorphous ice. Simulations of the pulse separation dependent spin-echo spectra for various scenarios demonstrate that a small-step frequency diffusion process, assigned to the presence of homonuclear oxygen–oxygen interactions, determines the shape evolution of the pulse-separation-dependent spectra.

Is It Possible to Follow the Structural Evolution of Water in "No-Man's Land" Using a Pulsed-Heating Procedure?

The anomalous increase in compressibility and heat capacity of supercooled water has been attributed to its structural transformation of into a four-coordinated liquid. Experiments revealed that κ_T and C_p peak at T_W^thermo ≈ 229 K [Kim et al. Science 2017, 358, 1589; Pathak et al. Proc. Natl. Acad. Sci. 2021, 118, e2018379118]. Recently, a pulsed heating procedure (PHP) was employed to interrogate the structure of water, reporting a steep increase in tetrahedrality around T_W^PHP = 210 ± 3 K [Kringle et al. Science 2020, 369, 1490]. This discrepancy questions whether water structure and thermodynamics are decoupled, or if the shift in T_W is an artifact of PHP. Here we implement PHP in molecular simulations. We find that the stationary states captured at the bottom of the pulse are not representative of the thermalized liquid or its inherent structure. Our analysis reveals a temperature-dependent distortion that shifts T_W^PHP to ∼20 K below T_W^thermo. We conclude that 2 orders of magnitude faster rates are required to sample water's inherent structure with PHP.