Liquid-liquid phase transition in deeply supercooled Stillinger-Weber silicon

The existence of a phase transition between two distinct liquid phases in single-component network-forming liquids (e.g. water, silica, silicon) has elicited considerable scientific interest. The challenge, both for experiments and simulations, is that the liquid-liquid phase transition (LLPT) occurs under deeply supercooled conditions, where crystallization occurs very rapidly. Thus, early evidence from numerical equation of state studies was challenged with the argument that slow spontaneous crystallization had been misinterpreted as evidence of a second liquid state. Rigorous free-energy calculations have subsequently confirmed the existence of a LLPT in some models of water, and exciting new experimental evidence has since supported these computational results. Similar results have so far not been found for silicon. Here, we present results from free-energy calculations performed for silicon modeled with the classical, empirical Stillinger-Weber-potential. Through a careful study employing state-of-the-art constrained simulation protocols and numerous checks for thermodynamic consistency, we find that there are two distinct metastable liquid states and a phase transition. Our results resolve a long-standing debate concerning the existence of a liquid-liquid transition in supercooled liquid silicon and address key questions regarding the nature of the phase transition and the associated critical point.

Correlated Local Fluctuations in the Hydrogen Bond Network of Liquid Water

The hypothesis that liquid water can separate into two phases in the supercooled state has been supported by recent experimental and theoretical studies. However, whether such structural inhomogeneity extends to ambient conditions is under intense debate. Due to the dynamic nature of the hydrogen bond network of liquid water, exploring its structure requires detailed insight into the collective motion of neighboring water molecules, a missing link that has not been examined so far. Here, highly sensitive quantum mechanical calculations detect that the time evolution of nearby hydrogen bonds is strongly correlated, revealing a direct mechanism for the appearance of short-range structural fluctuations in the hydrogen bond network of liquid water for the first time. This correlated dynamics is found to be closely connected to the static structural picture. The distortions from the tetrahedral structure do not occur independently but are correlated due to the preference of nearby donors and acceptors to be in similar environments. The existence of such cooperative fluctuations is further supported by the temperature dependence of the local structural evolution and explained by conventional analysis of localized orbitals. It was found that such correlated structural fluctuations are only observed on a short length scale in simulations at ambient conditions. The correlations of the nearby hydrogen bond pairs of liquid water unveiled here are expected to offer a new insight into connecting the dynamics of individual water molecules and the local structure of the hydrogen bond network.

Structural changes across thermodynamic maxima in supercooled liquid tellurium: A water-like scenario

Liquid polymorphism is an intriguing phenomenon that has been found in a few single-component systems, the most famous being water. By supercooling liquid Te to more than 130 K below its melting point and performing simultaneous small-angle and wide-angle X-ray scattering measurements, we observe clear maxima in its thermodynamic response functions around 615 K, suggesting the possible existence of liquid polymorphism. A close look at the underlying structural evolution shows the development of intermediate-range order upon cooling, most strongly around the thermodynamic maxima, which we attribute to bond-orientational ordering. The striking similarities between our results and those of water, despite the lack of hydrogen-bonding and tetrahedrality in Te, indicate that water-like anomalies may be a general phenomenon among liquid systems with competing bond- and density-ordering.

High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization

Cluster mixture models for liquid water at higher pressures suggest the need for water clusters of higher coordination and density than those commonly based on tetrahedral H-bonding motifs. We show here how proton-ordered water clusters of increased coordination and density can assemble from a starting cyclic tetramer or twisted bicyclic (Möbius-like) heptamer to form extended Aufbau sequences of stable two-, three-, and four-coordinate “windowpane” motifs. Such windowpane clusters exhibit sharply reduced (~90°) bond angles that differ appreciably from the tetrahedral angles of idealized crystalline ice I_h. Computed free energy and natural resonance theory (NRT) bond orders provide quantitative descriptors for the relative stabilities of clusters and strengths of individual coordinative linkages. The unity and consistency of NRT description is demonstrated to extend from familiar supra-integer bonds of the molecular regime to the near-zero bond orders of the weakest linkages in the present H-bond clusters. Our results serve to confirm that H-bonding exemplifies resonance–covalent (fractional) bonding in the sub-integer range and to further discount the dichotomous conceptions of “electrostatics” for intermolecular bonding vs. “covalency” for intramolecular bonding that still pervade much of freshman-level pedagogy and force-field methodology.

Liquid-Liquid Criticality in the WAIL Water Model

The hypothesis that the anomalous behavior of liquid water is related to the existence of a second critical point in deeply supercooled states has long been the subject of intense debate. Recent, sophisticated experiments designed to observe the transformation between the two subcritical liquids on nano- and microsecond time scales, along with demanding numerical simulations based on classical (rigid) models parametrized to reproduce thermodynamic properties of water, have provided support to this hypothesis. A stronger numerical proof requires demonstrating that the critical point, which occurs at temperatures and pressures far from those at which the models were optimized, is robust with respect to model parameterization, specifically with respect to incorporating additional physical effects. Here we show that a liquid-liquid critical point can be rigorously located also in the WAIL model of water [J. Chem. Phys. 137, 014510 (2012)], a model parameterized using ab-initio calculations only. The model incorporates two features not present in many previously-studied water models: it is both flexible and polarizable, properties which can significantly influence the phase behavior of water. The observation of the critical point in a model in which the water-water interaction is estimated using only quantum ab-initio calculations provides strong support to the viewpoint according to which the existence of two distinct liquids is a robust feature in the free energy landscape of supercooled water.

Nuclear quantum effects on the dynamics and glass behavior of a monatomic liquid with two liquid states

We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid-liquid phase transition and liquid-liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., Tg LDL(P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., Tg HDL(P) increases upon compression. (iii) NQE shift both Tg LDL(P) and Tg HDL(P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA-HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA-HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA-HDA transformation.

Low- and High-Density Unknown Waters at Ice-Water Interfaces

Experimental confirmation of liquid polymorphs of water, high-density liquid (HDL) and low-density liquid (LDL), is desired for understanding not only the liquid state of matter but also the origin of the mysterious properties of water. However, this remains challenging because the liquid-liquid critical point of water lies in experimentally inaccessible supercooling conditions known as "no man's land". Here, we show by in situ optical microscopy that droplets and layers of low- and high-density unknown waters (LDUW and HDUW) appear macroscopically depending upon ice polymorphs at non-equilibrium interfaces between water and ices under experimentally accessible (de)pressurization conditions. These unknown waters were found to have characteristic velocities (about 20 and 100 m/s for LDUW and HDUW, respectively) different from water (about 40 m/s) and quasi-liquid layers (QLLs) (about 2 and 0.2 m/s for droplet and layer forms of QLLs, respectively). Our discoveries provide insight on liquid polymorphism of water.

Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

For the past 50 years, researchers have sought molecular models that can accurately reproduce water's microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ∼208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. This work suggests that MB-pol may provide predictive capability for studies of liquid water's physical properties across broad ranges of thermodynamic states, including the so-called water's "no man's land" which is difficult to probe experimentally.

Evidence of a liquid-liquid phase transition in H[Formula: see text]O and D[Formula: see text]O from path-integral molecular dynamics simulations

We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H[Formula: see text]O and D[Formula: see text]O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density [Formula: see text], isothermal compressibility [Formula: see text], and self-diffusion coefficients D(T) of H[Formula: see text]O and D[Formula: see text]O are in excellent agreement with available experimental data; the isobaric heat capacity [Formula: see text] obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H[Formula: see text]O and D[Formula: see text]O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H[Formula: see text]O and D[Formula: see text]O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H[Formula: see text]O, from PIMD simulations, is located at [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D[Formula: see text]O is estimated to be [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of [Formula: see text] for D[Formula: see text]O and, particularly, H[Formula: see text]O suggest that improved water models are needed for the study of supercooled water.

Density Functional Theory of Water with the Machine-Learned DM21 Functional

The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn-Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. We find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.

Finite-size scaling and thermodynamics of model supercooled liquids: long-range concentration fluctuations and the role of attractive interactions

We compute partial structure factors, Kirkwood-Buff integrals (KBIs) and chemical potentials of model supercooled liquids with and without attractive interactions. We aim at investigating whether relatively small differences in the tail of the radial distribution functions result in contrasting thermodynamic properties. Our results suggest that the attractive potential favours the nucleation of long-range structures. Indeed, upon decreasing temperature, Bathia-Thornton structure factors display anomalous behaviour in the k→0 limit. KBIs extrapolated to the thermodynamic limit confirm this picture, and excess coordination numbers identify the anomaly with long-range concentration fluctuations. By contrast, the purely repulsive system remains perfectly miscible for the same temperature interval and only reveals qualitatively similar concentration fluctuations in the crystalline state. Furthermore, differences in both isothermal compressibilities and chemical potentials show that thermodynamics is not entirely governed by the short-range repulsive part of the interaction potential, emphasising the nonperturbative role of attractive interactions. Finally, at higher density, where both systems display nearly identical dynamical properties and repulsive interactions become dominant, the anomaly disappears, and both systems also exhibit similar thermodynamic properties.

Performance of the new biological small- and wide-angle X-ray scattering beamline 13A at the Taiwan Photon Source

A new endstation for biological small- and wide-angle X-ray scattering is detailed, which provides development opportunities for studying correlated local and global structures of biomolecules in solution.

Recent developments in the instrumentation and data analysis of synchrotron small-angle X-ray scattering (SAXS) on biomolecules in solution have made biological SAXS (BioSAXS) a mature and popular tool in structural biology. This article reports on an advanced endstation developed at beamline 13A of the 3.0 GeV Taiwan Photon Source for biological small- and wide-angle X-ray scattering (SAXS–WAXS or SWAXS). The endstation features an in-vacuum SWAXS detection system comprising two mobile area detectors (Eiger X 9M/1M) and an online size-exclusion chromatography system incorporating several optical probes including a UV–Vis absorption spectrometer and refractometer. The instrumentation and automation allow simultaneous SAXS–WAXS data collection and data reduction for high-throughput biomolecular conformation and composition determinations. The performance of the endstation is illustrated with the SWAXS data collected for several model proteins in solution, covering a scattering vector magnitude q across three orders of magnitude. The crystal-model fittings to the data in the q range ∼0.005–2.0 Å⁻¹ indicate high similarity of the solution structures of the proteins to their crystalline forms, except for some subtle hydration-dependent local details. These results open up new horizons of SWAXS in studying correlated local and global structures of biomolecules in solution.

How do interfaces alter the dynamics of supercooled water?

The structure of liquid water in the proximity of an interface can deviate significantly from that of bulk water, with surface-induced structural perturbations typically converging to bulk values at about ∼1 nm from the interface. While these structural changes are well established it is, in contrast, less clear how an interface perturbs the dynamics of water molecules within the liquid. Here, through an extensive set of molecular dynamics simulations of supercooled bulk and interfacial water films and nano-droplets, we observe the formation of persistent, spatially extended dynamical domains in which the average mobility varies as a function of the distance from the interface. This is in stark contrast with the dynamical heterogeneity observed in bulk water, where these domains average out spatially over time. We also find that the dynamical response of water to an interface depends critically on the nature of the interface and on the choice of interface definition. Overall these results reveal a richness in the dynamics of interfacial water that opens up the prospect of tuning the dynamical response of water through specific modifications of the interface structure or confining material.

Following the Crystallization of Amorphous Ice after Ultrafast Laser Heating

Using time-resolved wide-angle X-ray scattering, we investigated the early stages (10 μs–1 ms) of crystallization of supercooled water, obtained by the ultrafast heating of high- and low-density amorphous ice (HDA and LDA) up to a temperature T = 205 K ± 10 K. We have determined that the crystallizing phase is stacking disordered ice (I_sd), with a maximum cubicity of χ = 0.6, in agreement with predictions from molecular dynamics simulations at similar temperatures. However, we note that a growing small portion of hexagonal ice (I_h) was also observed, suggesting that within our timeframe, I_sd starts annealing into I_h. The onset of crystallization, in both amorphous ice, occurs at a similar temperature, but the observed final crystalline fraction in the LDA sample is considerably lower than that in the HDA sample. We attribute this discrepancy to the thickness difference between the two samples.

Structural characteristics of low-density environments in liquid water

The existence of two structural forms in liquid water has been a point of discussion for a long time. A phase transition between these two forms of liquid water has been proposed based on evidence from molecular simulations, and experiments have also been very recently able to track the proposed transition of the low-density liquid form to the high-density liquid form. We propose to use the average angle an oxygen atom makes with its neighbors to describe the structural environment of a water molecule. The distribution of this order parameter is observed to have two peaks with one peak at ∼109.5^{∘}, corresponding to the internal angle of a regular tetrahedron, indicating tetrahedral arrangement. The other peak corresponds to an environment with a tighter arrangement of neighboring molecules. The distribution of O-O-O angles is decomposed into two skewed distributions to estimate the fractions of the two liquid forms in water. A good similarity is observed between the temperature and pressure trends of fractions of locally favored tetrahedral structure (LFTS) form estimated using the new order parameter and the reports in the literature, over a range of temperatures and pressures. We also compare the structural environments indicated by different order parameters and find that the order parameter proposed in this paper captures the structure of first solvation shell of the LFTS accurately.

Isotope effects on the structural transformation and relaxation of deeply supercooled water

We have examined the structure of supercooled liquid D2O as a function of temperature between 185 and 255 K using pulsed laser heating to rapidly heat and cool the sample on a nanosecond timescale. The liquid structure can be represented as a linear combination of two structural motifs, with a transition between them described by a logistic function centered at 218 K with a width of 10 K. The relaxation to a metastable state, which occurred prior to crystallization, exhibited nonexponential kinetics with a rate that was dependent on the initial structural configuration. When the temperature is scaled by the temperature of maximum density, which is an isostructural point of the isotopologues, the structural transition and the non-equilibrium relaxation kinetics of D2O agree remarkably well with those for H2O.

Phase transition of supercooled water confined in cooperative two-state domain

The question of 'what is the structure of water?' has been regarded as one of the major scientific conundrums in condensed-matter physics due to the complex phase behavior and condensed structure of supercooled water. Great effort has been made so far using both theoretical analysis based on various mathematical models and computer simulations such as molecular dynamics and first-principle. However, these theoretical and simulation studies often do not have strong evidences of condensed-matter physics to support. In this study, a cooperative domain model is formulated to describe the dynamic phase transition of supercooled water between supercooled water and amorphous ice, both of which are composed of low- and high-density liquid water. Free volume theory is initially employed to identify the working principle of dynamic phase transition and its connection to glass transition in the supercooled water. Then a cooperative two-state model is developed to characterize the dynamic anomalies of supercooled water, including density, viscosity and self-diffusion coefficient. Finally, the proposed model is verified using the experimental results reported in literature.

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.

Folded network and structural transition in molten tin

The fundamental relationships between the structure and properties of liquids are far from being well understood. For instance, the structural origins of many liquid anomalies still remain unclear, but liquid-liquid transitions (LLT) are believed to hold a key. However, experimental demonstrations of LLTs have been rather challenging. Here, we report experimental and theoretical evidence of a second-order-like LLT in molten tin, one which favors a percolating covalent bond network at high temperatures. The observed structural transition originates from the fluctuating metallic/covalent behavior of atomic bonding, and consequently a new paradigm of liquid structure emerges. The liquid structure, described in the form of a folded network, bridges two well-established structural models for disordered systems, i.e., the random packing of hard-spheres and a continuous random network, offering a large structural midground for liquids and glasses. Our findings provide an unparalleled physical picture of the atomic arrangement for a plethora of liquids, shedding light on the thermodynamic and dynamic anomalies of liquids but also entailing far-reaching implications for studying liquid polyamorphism and dynamical transitions in liquids.

Unraveling the structural origin of liquid anomalies remains a challenging topic. Xu et al. propose a folded-network structural model for molten tin and provide insights into the observed second-order-like structural transition.

Comparison of experimental times in T1-D and D-T2 correlation experiments in single-sided NMR

Diffusion-relaxation correlation experiments in nuclear magnetic resonance are a powerful technique for the characterization of fluid dynamics in confined geometries or soft matter, in which relaxation may be either spin-spin (T₂) or spin-lattice (T₁). The general approach is to acquire a set of bidimensional data in which diffusion is codified by the evolution of the magnetization with either Hahn or stimulated echoes (STE) in the presence of a constant magnetic field gradient. While T₂ is codified by a Carr-Purcell-Meiboom-Gil (CPMG) sequence, T₁ is either encoded by saturation or inversion-recovery methods. In this work, we analyse the measurement time of diffusion-relaxation times in single-sided NMR and show that T₁-D acquisition is always shorter than D-T₂. Depending on the hardware characteristics, this time reduction can be up to an order of magnitude. We present analytical calculations and examples in model porous media saturated with water and in a dairy product.

Temperature-Dependent Liquid Water Structure for Individual Micron-Sized, Supercooled Aqueous Droplets with Inclusions

Herein, we measure the water structure for individual micron-sized droplets of water, salt water, and water containing biologically and marine relevant atmospheric inclusions as a function of temperature. Individual droplets, formed on a hydrophobic substrate, are analyzed with micro-Raman spectroscopy. Analysis of the Raman spectra in the O-H stretching region shows that the equilibrium of partially and fully hydrogen-bonding water interactions change as temperature decreases up until there is a phase transition to form ice. Using these temperature-dependent measurements, the thermodynamic parameters for the interchange between partially and fully hydrogen-bonded water (PHW ⇄ FHW) for different supercooled droplets (water, salt water, and water containing biologically and marine relevant atmospheric inclusions) have been determined.

The pervasive impact of critical fluctuations in liquid-liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Overall, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

Anomalous temperature dependence of the experimental x-ray structure factor of supercooled water

The structural changes of water upon deep supercooling were studied through wide-angle x-ray scattering at SwissFEL. The experimental setup had a momentum transfer range of 4.5 Å^-1, which covered the principal doublet of the x-ray structure factor of water. The oxygen-oxygen structure factor was obtained for temperatures down to 228.5 ± 0.6 K. Similar to previous studies, the second diffraction peak increased strongly in amplitude as the structural change accelerated toward a local tetrahedral structure upon deep supercooling. We also observed an anomalous trend for the second peak position of the oxygen-oxygen structure factor (q₂). We found that q₂ exhibits an unprecedented positive partial derivative with respect to temperature for temperatures below 236 K. Based on Fourier inversion of our experimental data combined with reference data, we propose that the anomalous q₂ shift originates from that a repeat spacing in the tetrahedral network, associated with all peaks in the oxygen-oxygen pair-correlation function, gives rise to a less dense local ordering that resembles that of low-density amorphous ice. The findings are consistent with that liquid water consists of a pentamer-based hydrogen-bonded network with low density upon deep supercooling.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

Exploring the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes

The breakdown of Stokes-Einstein relation in liquid water is one of the many anomalies that take place upon cooling and indicates the decoupling of diffusion and viscosity. It is hypothesized that these anomalies manifest due to the appearance of nanometer-scale spatial fluctuations, which become increasingly pronounced in the supercooled regime. Here, we explore the validity of the Stokes-Einstein relation in supercooled water using nanomolecular probes. We capture the diffusive dynamics of the probes using dynamic light scattering and target dynamics at different length scales by varying the probe size, from ≈100 nm silica spheres to molecular-sized polyhydroxylated fullerenes (≈1 nm). We find that all the studied probes, independent of size, display similar diffusive dynamics with an Arrhenius activation energy of ≈23 kJ mol^-1. Analysis of the diffusion coefficient further indicates that the probes, independent of their size, experience similar dynamic environment, which coincides with the macroscopic viscosity, while single water molecules effectively experience a comparatively lower viscosity. Finally, we conclude that our results indicate that the Stokes-Einstein relation is preserved for diffusion of probes in supercooled water T ≥ 260 K with size as small as ≈1 nm.

Structure and dynamics of nanoconfined water and aqueous solutions

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Minimal Microscopic Model for Liquid Polyamorphism and Waterlike Anomalies

Liquid polyamorphism is the intriguing possibility for a single component substance to exist in multiple liquid phases. We propose a minimal model for this phenomenon. Starting with a binary lattice model with critical azeotropy and liquid-liquid demixing, we allow interconversion of the two species, turning the system into a single-component fluid with two states differing in energy and entropy. Unveiling the phase diagram of the noninterconverting binary mixture gives unprecedented insight on the phase behaviors accessible to the interconverting fluid, such as a liquid-liquid transition with a critical point, or a singularity-free scenario, exhibiting thermodynamic anomalies without polyamorphism. The model provides a unified theoretical framework to describe supercooled water and a variety of polyamorphic liquids with waterlike anomalies.

Phase classification using neural networks: application to supercooled, polymorphic core-softened mixtures

Characterization of phases of soft matter systems is a challenge faced in many physical chemical problems. For polymorphic fluids it is an even greater challenge. Specifically, glass forming fluids, as water, can have, besides solid polymorphism, more than one liquid and glassy phases, and even a liquid-liquid critical point. In this sense, we apply a neural network algorithm to analyze the phase behavior of a mixture of core-softened fluids that interact through the continuous-shouldered well (CSW) potential, which have liquid polymorphism and liquid-liquid critical points, similar to water. We also apply the neural network to mixtures of CSW fluids and core-softened alcohols models. We combine and expand methods based on bond-orientational order parameters to study mixtures, applied to mixtures of hardcore fluids and to supercooled water, to include longer range coordination shells. With this, the trained neural network was able to properly predict the crystalline solid phases, the fluid phases and the amorphous phase for the pure CSW and CSW-alcohols mixtures with high efficiency. More than this, information about the phase populations, obtained from the network approach, can help verify if the phase transition is continuous or discontinuous, and also to interpret how the metastable amorphous region spreads along the stable high density fluid phase. These findings help to understand the behavior of supercooled polymorphic fluids and extend the comprehension of how amphiphilic solutes affect the phases behavior.

From Critical Point to Critical Point: The Two-States Model Describes Liquid Water Self-Diffusion from 623 to 126 K

Liquid’s behaviour, when close to critical points, is of extreme importance both for fundamental research and industrial applications. A detailed knowledge of the structural–dynamical correlations in their proximity is still today a target to reach. Liquid water anomalies are ascribed to the presence of a second liquid–liquid critical point, which seems to be located in the very deep supercooled regime, even below 200 K and at pressure around 2 kbar. In this work, the thermal behaviour of the self-diffusion coefficient for liquid water is analyzed, in terms of a two-states model, for the first time in a very wide thermal region (126 K < T < 623 K), including those of the two critical points. Further, the corresponding configurational entropy and isobaric-specific heat have been evaluated within the same interval. The two liquid states correspond to high and low-density water local structures that play a primary role on water dynamical behavior over 500 K.

Segregation on the nanoscale coupled to liquid water polyamorphism in supercooled aqueous ionic-liquid solution

The most intriguing hypothesis explaining many water anomalies is a metastable liquid-liquid phase transition (LLPT) at high pressure and low temperatures, experimentally hidden by homogeneous nucleation. Recent infrared spectroscopic experiments showed that upon addition of hydrazinium trifluoroacetate to water, the supercooled ionic solution undergoes a sharp, reversible LLPT at ambient pressure, possible offspring of that in pure water. Here, we calculate the temperature-dependent signature of the OH-stretching band, reporting on the low/high density phase of water, in neat water and in the same experimentally investigated ionic solution. The comparison between the infrared signature of the pure liquid and that of the ionic solution can be achieved only computationally, providing insight into the nature of the experimentally observed phase transition and allowing us to investigate the effects of ionic compounds on the high to low density supercooled liquid water transition. We show that the experimentally observed crossover behavior in the ionic solution can be reproduced only if the phase transition between the low- and high-density liquid states of water is coupled to a mixing-unmixing transition between the water component and the ions: at low temperatures, water and ions are separated and the water component is a low density liquid. At high temperatures, water and ions get mixed and the water component is a high-density liquid. The separation at low temperatures into ion-rich and ion-poor regions allows unveiling the polyamorphic nature of liquid water, leading to a crossover behavior resembling that observed in supercooled neat water under high pressure.

Wide-angle X-ray scattering and molecular dynamics simulations of supercooled protein hydration water

Understanding the mechanism responsible for the protein low-temperature crossover observed at T≈ 220 K can help us improve current cryopreservation technologies. This crossover is associated with changes in the dynamics of the system, such as in the mean-squared displacement, whereas experimental evidence of structural changes is sparse. Here we investigate hydrated lysozyme proteins by using a combination of wide-angle X-ray scattering and molecular dynamics (MD) simulations. Experimentally we suppress crystallization by accurate control of the protein hydration level, which allows access to temperatures down to T = 175 K. The experimental data indicate that the scattering intensity peak at Q = 1.54 Å-1, attributed to interatomic distances, exhibits temperature-dependent changes upon cooling. In the MD simulations it is possible to decompose the water and protein contributions and we observe that, while the protein component is nearly temperature independent, the hydration water peak shifts in a fashion similar to that of bulk water. The observed trends are analysed by using the water-water and water-protein radial distribution functions, which indicate changes in the local probability density of hydration water.

Reversal of the Temperature Dependence of Hydrophobic Hydration in Supercooled Water

Using simulations and theory, we examine the enthalpy and entropy of hydrophobic hydration which exhibit minima in supercooled water, contrasting with the monotonically increasing temperature dependence traditionally ascribed to these properties. The enthalpy/entropy minima are marked by a negative to positive sign change in the heat capacity at a size-dependent reversal temperature. A Gaussian fluctuation theory accurately captures the reversal temperature, tracing it to water's distinct thermal expansivity and compressibility influenced by its metastable liquid-liquid critical point.

Effects of Resonant Electromagnetic Fields on Biofilm Formation in Pseudomonas aeruginosa

The global rise of antimicrobial resistance (AMR) constitutes a future health threat and dictates a need to explore alternative and non-chemical approaches. The aim of this study was to explore the use of weak resonant electromagnetic fields as a method to disrupt biofilm formation of a pathogenic bacterium in cystic fibrosis patients. We developed a bioresonance laboratory setup able to distinguish between changes in planktonic growth and changes in biofilm formation and showed that certain resonant frequencies were able to affect biofilm formation without affecting planktonic growth. In addition, we show that the ambient day-to-day magnetic field affects biofilm formation in a non-consistent manner. Overall, we conclude that our assay is suitable for studying the potential of resonant magnetic fields as a treatment and prevention strategy to prevent biofilm infections, and that certain resonant frequencies may be used as future medical applications to combat antimicrobial resistance.

Dynamical crossover and its connection to the Widom line in supercooled TIP4P/Ice water

We perform molecular dynamics simulations with the TIP4P/Ice water model to characterize the relationship between dynamics and thermodynamics of liquid water in the supercooled region. We calculate the relevant properties of the phase diagram, and we find that TIP4P/Ice presents a retracing line of density maxima, similar to what was previously found for atomistic water models and models of other tetrahedral liquids. For this model, a liquid-liquid critical point between a high-density liquid and a low-density liquid was recently found. We compute the lines of the maxima of isothermal compressibility and the minima of the coefficient of thermal expansion in the one phase region, and we show that these lines point to the liquid-liquid critical point while collapsing on the Widom line. This line is the line of the maxima of correlation length that emanates from a second order critical point in the one phase region. Supercooled water was found to follow mode coupling theory and to undergo a transition from a fragile to a strong behavior right at the crossing of the Widom line. We find here that this phenomenology also happens for TIP4P/Ice. Our results appear, therefore, to be a general characteristic of supercooled water, which does not depend on the interaction potential used, and they reinforce the idea that the dynamical crossover from a region where the relaxation mechanism is dominated by cage relaxation to a region where cages are frozen and hopping dominates is correlated in water to a phase transition between a high-density liquid and a low-density liquid.

Core-softened water-alcohol mixtures: the solute-size effects

Water is the most anomalous material on Earth, with a long list of thermodynamic, dynamic and structural behaviors that deviate from what is expected. Recent studies have indicated that these anomalies may be related to a competition between two liquids, which means that water has a potential liquid-liquid phase transition (LLPT) that ends at a liquid-liquid critical point (LLCP). In a recent study [J. Mol. Liq., 2020, 320, 114420], using molecular dynamics simulations and a core-softened potential approach, we have shown that adding a simple solute such as methanol can "kill" the density-anomalous behavior as the LLCP is suppressed by spontaneous crystallization in a hexagonal close packing (HCP) crystal near the LLPT. Now, we extend this work to realize how longer-chain alcohols will affect the complex behavior of water-alcohol mixtures in the supercooled regime. Besides core-softened (CS) methanol, ethanol and 1-propanol were added to a system of identical particles that interact through the continuous shouldered well (CSW) potential. We observed that the density anomaly gradually decreases its extension in phase diagrams until it disappears with the growth of the non-polar chain and the alcohol concentration, different from the liquid-liquid phase transition (and the LLCP), which remained present in all analyzed mixtures, according to Nature, 2001, 409, 692. For our model, the longer non-polar chains and higher concentrations gradually impact the competition between the scales in the CS potential, leading to a gradual disappearing of the anomalies until the TMD total disappearance is observed when the first coordination shell structure is also affected: short-range ordering is favored, leading to less competition between short- and long-range ordering and, consequently, to the extinction of anomalies. Also, the non-polar chain size and concentration have an effect on the solid phases, favoring the hexagonal close packed (HCP) solid and the amorphous solid phase over the body-centered cubic (BCC) crystal. These findings help to elucidate the behavior of water solutions in the supercooled regime and indicate that the LLCP can be observed in systems without density-anomalous behavior.

A high-accuracy machine-learning water model for exploring water nanocluster structures

Water, the most important molecule on the Earth, possesses many essential and unique physical properties that are far from completely understood, partly due to serious difficulties in identifying the precise microscopic structures of water. Hence, identifying the structures of water nanoclusters is a fundamental and challenging issue for studies on the relationship between the macroscopic physical properties of water and its microscopic structures. For large-scale simulations (at the level of nm and ns) of water nanoclusters, a calculation method with simultaneous accuracy at the level of quantum chemistry and efficiency at the level of an empirical potential method is in great demand. Herein, a machine-learning (ML) water model was utilized to explore the microscopic structural features at different length scales for water nanoclusters with a size up to several nm. The ML water model can be employed to efficiently predict the structures of water nanoclusters with a similar accuracy to that of density functional theory and with substantially lower computational resource demands. To validate the low-lying structure search results with experimental spectral results, an ML water model combined with velocity autocorrelation function analysis was used to simulate the vibrational spectra of water nanoclusters with up to thousands of water molecules. By comparing the simulated and experimentally recorded vibrational spectra, the atomic structures determined by a simulation based on the ML water model are all verified. To demonstrate its ability to represent water's structural evolution at large length and time scales, the ML water model was employed to model the structural evolution during the crystal-liquid transition, and the phase transition temperatures of water clusters with different sizes were precisely predicted. The ML water model provides an efficient theoretical calculation tool for exploring the structures and physical properties of water and their relationships, especially for clusters with relatively large sizes and processes with relatively long durations.

Hydrophilic and Hydrophobic Effects on the Structure and Themodynamic Properties of Confined Water: Water in Solutions

NMR spectroscopy is used in the temperature range 180–350 K to study the local order and transport properties of pure liquid water (bulk and confined) and its solutions with glycerol and methanol at different molar fractions. We focused our interest on the hydrophobic effects (HE), i.e., the competition between hydrophilic and hydrophobic interactions. Nowadays, compared to hydrophilicity, little is known about hydrophobicity. Therefore, the main purpose of this study is to gain new information about hydrophobicity. As the liquid water properties are dominated by polymorphism (two coexisting liquid phases of high and low density) due to hydrogen bond interactions (HB), creating (especially in the supercooled regime) the tetrahedral networking, we focused our interest to the HE of these structures. We measured the relaxation times (T1 and T2) and the self-diffusion (DS). From these times, we took advantage of the NMR property to follow the behaviors of each molecular component (the hydrophilic and hydrophobic groups) separately. In contrast, DS is studied in terms of the Adam–Gibbs model by obtaining the configurational entropy (Sconf) and the specific heat contributions (CP,conf). We find that, for the HE, all of the studied quantities behave differently. For water–glycerol, the HB interaction is dominant for all conditions; water–methanol, two different T-regions above and below 265 K are observable, dominated by hydrophobicity and hydrophilicity, respectively. Below this temperature, where the LDL phase and the HB network develops and grows, with the times and CP,conf change behaviors leading to maxima and minima. Above it, the HB becomes weak and less stable, the HDL dominates, and hydrophobicity determines the solution.

From Femtoseconds to Hours—Measuring Dynamics over 18 Orders of Magnitude with Coherent X-rays

X-ray photon correlation spectroscopy (XPCS) enables the study of sample dynamics between micrometer and atomic length scales. As a coherent scattering technique, it benefits from the increased brilliance of the next-generation synchrotron radiation and Free-Electron Laser (FEL) sources. In this article, we will introduce the XPCS concepts and review the latest developments of XPCS with special attention on the extension of accessible time scales to sub-μs and the application of XPCS at FELs. Furthermore, we will discuss future opportunities of XPCS and the related technique X-ray speckle visibility spectroscopy (XSVS) at new X-ray sources. Due to its particular signal-to-noise ratio, the time scales accessible by XPCS scale with the square of the coherent flux, allowing to dramatically extend its applications. This will soon enable studies over more than 18 orders of magnitude in time by XPCS and XSVS.

Mesostructuring in Liquid-Liquid Extraction Organic Phases Originating from Critical Points

Organic phase structure plays an important role in solute extraction energetics and phase behavior of liquid-liquid extraction (LLE) systems. For a binary extractant (amphiphile)/solvent mixture of relevance to LLE, we find that the organic phase mesostructuring is consistent with extractant concentration fluctuations as the compositional isotherm traverses the Widom line above its liquid-liquid critical point. This reveals a different mechanism for the well-documented heterogeneities in LLE organic phases that are typically attributed to micellization.

Liquid-liquid phase separation in an inhomogeneous ternary colloid-polymer mixture

Suspended colloids are often considered as models for molecules, which are sufficiently big so that they can be observed directly in (light) microscopes and for which the effective interaction among each other can be tailored. The Asakura-Oosawa model of ideal colloid-polymer mixtures captures the idea of tuning the interaction between the colloids via a potential, which possesses a range set by the size of the polymers and an attractive strength characterized by the (reservoir) number density of the polymers, which plays the role of an inverse temperature. The celebrated Asakura-Oosawa depletion potential allows one to recreate the bulk phase diagram of a simple fluid by employing a colloid-polymer mixture. This has been verified in theory, by computer simulations, and via experiments. Here, we study the phase behavior of a confined colloid-polymer mixture with two polymer species. The sizes and densities are chosen such that the resulting bulk phase diagram exhibits a second stable critical point within the framework of the classical density functional theory. Our results suggest that a suitably tuned colloid-polymer mixture can be an interesting model system to study fluids with two critical points.

Manifestations of metastable criticality in the long-range structure of model water glasses

Much attention has been devoted to water’s metastable phase behavior, including polyamorphism (multiple amorphous solid phases), and the hypothesized liquid-liquid transition and associated critical point. However, the possible relationship between these phenomena remains incompletely understood. Using molecular dynamics simulations of the realistic TIP4P/2005 model, we found a striking signature of the liquid-liquid critical point in the structure of water glasses, manifested as a pronounced increase in long-range density fluctuations at pressures proximate to the critical pressure. By contrast, these signatures were absent in glasses of two model systems that lack a critical point. We also characterized the departure from equilibrium upon vitrification via the non-equilibrium index; water-like systems exhibited a strong pressure dependence in this metric, whereas simple liquids did not. These results reflect a surprising relationship between the metastable equilibrium phenomenon of liquid-liquid criticality and the non-equilibrium structure of glassy water, with implications for our understanding of water phase behavior and glass physics. Our calculations suggest a possible experimental route to probing the existence of the liquid-liquid transition in water and other fluids.

The subtle connections between water’s supercooled liquid and glassy states are difficult to characterize. Gartner et al. suggest with MD simulations that the long-range structure of glassy water may reflect signatures of water’s debated second critical point in the supercooled liquid.

Two Liquid-Liquid Phase Transitions in Confined Water Nanofilms

The stories behind supercooled bulk and confined water can be different. Bulk water has a metastable liquid-liquid phase transition at deeply supercooled conditions, but the existence of such a phenomenon in confined water is in question. Herein we show simulation results of first-order phase transitions between high- and low-density liquid (HDL and LDL) in confined water in both positive and negative pressures. A mid-density state between these two local states appears, which lets the transition show the hysteresis loop with transiently stable intermediate states. On the basis of Landau theory that we have adapted for mixing of moieties with high- and low-density states, we explain the phase transitions with the order parameter-dependent free energy change which is governed by second- to higher-order interactions among those moieties.

Universal Nucleation Behavior of Sheared Systems

Using molecular simulations and a modified classical nucleation theory, we study the nucleation, under flow, of a variety of liquids: different water models, Lennard-Jones, and hard sphere colloids. Our approach enables us to analyze a wide range of shear rates inaccessible to brute-force simulations. Our results reveal that the variation of the nucleation rate with shear is universal. A simplified version of the theory successfully captures the nonmonotonic temperature dependence of the nucleation behavior, which is shown to originate from the violation of the Stokes-Einstein relation.

Structural and topological changes across the liquid-liquid transition in water

It has recently been shown that the TIP4P/Ice model of water can be studied numerically in metastable equilibrium at and below its liquid-liquid critical temperature. We report here simulations along a subcritical isotherm, for which two liquid states with the same pressure and temperature but different density can be equilibrated. This allows for a clear visualization of the structural changes taking place across the transition. We specifically focus on how the topological properties of the H-bond network change across the liquid-liquid transition. Our results demonstrate that the structure of the high-density liquid, characterized by the existence of interstitial molecules and commonly explained in terms of the collapse of the second neighbor shell, actually originates from the folding back of long rings, bringing pairs of molecules separated by several hydrogen-bonds close by in space.

Structural relaxation and crystallization in supercooled water from 170 to 260 K

The origin of water's anomalous properties has been debated for decades. Resolution of the problem is hindered by a lack of experimental data in a crucial region of temperatures, T, and pressures where supercooled water rapidly crystallizes-a region often referred to as "no man's land." A recently developed technique where water is heated and cooled at rates greater than 10⁹ K/s now enables experiments in this region. Here, it is used to investigate the structural relaxation and crystallization of deeply supercooled water for 170 K < T < 260 K. Water's relaxation toward a new equilibrium structure depends on its initial structure with hyperquenched glassy water (HQW) typically relaxing more quickly than low-density amorphous solid water (LDA). For HQW and T > 230 K, simple exponential relaxation kinetics is observed. For HQW at lower temperatures, increasingly nonexponential relaxation is observed, which is consistent with the dynamics expected on a rough potential energy landscape. For LDA, approximately exponential relaxation is observed for T > 230 K and T < 200 K, with nonexponential relaxation only at intermediate temperatures. At all temperatures, water's structure can be reproduced by a linear combination of two, local structural motifs, and we show that a simple model accounts for the complex kinetics within this context. The relaxation time, τ _rel , is always shorter than the crystallization time, τ _xtal For HQW, the ratio, τ _xtal /τ _rel , goes through a minimum at ∼198 K where the ratio is about 60.