Diffusive dynamics during the high-to-low density transition in amorphous ice

Multicomponent dynamics in amorphous ice studied using X-ray photon correlation spectroscopy at elevated pressure and cryogenic temperatures

Knowing the pressure dependence of glass forming liquids is important in various contexts. Here, we study the case of supercooled water, which has at least two different amorphous states with different densities. The pressure dependencies of the two glass transitions are predicted to show opposite behaviour, crossing in the P-T plane at elevated pressure. The experimental identification of the glass transition at elevated pressure and cryo-conditions is technically difficult. Moreover, in the case of amorphous ices, the glass transition is interrupted by crystallization, which makes it even more challenging. We show the feasibility of performing X-ray photon correlation spectroscopy experiments at elevated pressure using a diamond anvil cell at cryogenic temperatures. We observe two dynamic components when approaching the glass transition temperature. For high-density amorphous ice at a pressure of around (0.08 ± 0.02) GPa we determine the glass transition to be at higher temperatures compared to ambient conditions.

Proton diffusion and hydrogen/deuterium exchange in amorphous solid water at temperatures from 114 to 134 K

The reaction coefficient for hydrogen/deuterium (H/D) exchange and the diffusion of hydrated excess protons within amorphous solid water (ASW) are characterized as a function of temperature. For these experiments, water films are deposited on a Pt(111) substrate at 108 K, and reactions with pre-adsorbed hydrogen atoms produce hydrated protons. Upon heating, protons diffuse within the water, and H/D exchange occurs when they encounter D2O probe molecules deposited in the films. The time-dependent concentration of D2O is monitored with infrared spectroscopy, and it indicates the protons diffusion from the substrate and establish an equilibrium distribution prior to significant H/D exchange for temperatures 114 K ≤T≤ 134 K. By controlling the distance between the D2O molecules and the substrate, we probe the distribution of protons within the film. It decays as x-2 for the examined range of x (12-52 nm) due to the electric field that develops between the diffusing protons and their image charges in the metal substrate. This agrees with the theoretical distance scaling for the equilibrated proton concentration in a dielectric near a metal boundary. From the proton concentration and the measured D2O decay rate, a lower bound for the proton diffusion coefficient ranging from 10-20 m2/s at 114 K to 10-18 m2/s at 134 K is estimated. The diffusion coefficient has an activation energy of 0.40 eV, which is comparable to energies reported for molecular translations and rotations of H2O, suggesting they may play a critical role in the proton diffusion mechanism within ASW.

Nuclear quantum effects on glassy water under pressure: Vitrification and pressure-induced transformations

We perform classical molecular dynamics (MD) and path-integral MD (PIMD) simulations of H2O and D2O using the q-TIP4P/F model over a wide range of temperatures and pressures to study the nuclear quantum effects (NQEs) on (i) the vitrification of liquid water upon isobaric cooling at different pressures and (ii) pressure-induced transformations at constant temperature between low-density amorphous and high-density amorphous ice (LDA and HDA) and hexagonal ice Ih and HDA. Upon isobaric cooling, classical and quantum H2O and D2O vitrify into a continuum of intermediate amorphous ices (IA), with densities in-between those of LDA and HDA (depending on pressure). Importantly, the density of the IA varies considerably if NQEs are included (similar conclusions hold for ice Ih at all pressures studied). While the structure of the IA is not very sensitive to NQE, the geometry of the hydrogen-bond (HB) is. NQE leads to longer and less linear HB in LDA, HDA, and ice Ih than found in the classical case. Interestingly, the delocalization of the H/D atoms is non-negligible and identical in LDA, HDA, and ice Ih at all pressures studied. Our isothermal compression/decompression MD/PIMD simulations show that classical and quantum H2O and D2O all exhibit LDA-HDA and ice Ih-HDA transformations, consistent with experiments. The inclusion of NQE leads to a softer HB-network, which lowers slightly the LDA/ice Ih-to-HDA transformation pressures. Interestingly, the HB in HDA is longer and less linear than in LDA, which is counterintuitive given that HDA is ≈25% denser than LDA. Overall, our results show that, while classical computer simulations provide the correct qualitative phenomenology of ice and glassy water, NQEs are necessary for a quantitative description.

Supercritical density fluctuations and structural heterogeneity in supercooled water-glycerol microdroplets

Recent experiments and theoretical studies strongly indicate that water exhibits a liquid-liquid phase transition (LLPT) in the supercooled domain. An open question is how the LLPT of water can affect the properties of aqueous solutions. Here, we study the structural and thermodynamic properties of supercooled glycerol-water microdroplets at dilute conditions (χg = 3.2% glycerol mole fraction). The combination of rapid evaporative cooling with femtosecond X-ray scattering allows us to outrun crystallization and gain access to the deeply supercooled regime down to T = 229.3 K. We find that the density fluctuations of the glycerol-water solution or, equivalently, its isothermal compressibility, κT, increases upon cooling. This is confirmed by molecular dynamics simulations, which indicate that the presence of glycerol shifts the temperature of maximum κT from T = 230 K in pure water down to T = 223 K in the solution. Our findings elucidate the interplay between the complex behavior of water, including its LLPT, and the properties of aqueous solutions at low temperatures, which can have practical consequences in cryogenic biological applications and cryopreservation techniques.

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.

Complexity of confined water vitrification and its glass transition temperature

The ability of vitrification when crossing the glass transition temperature (Tg) of confined and bulk water is crucial for myriad phenomena in diverse fields, ranging from the cryopreservation of organs and food to the development of cryoenzymatic reactions, frost damage to buildings, and atmospheric water. However, determining water's Tg remains a major challenge. Here, we elucidate the glass transition of water by analyzing the calorimetric behavior of nano-confined water across various pore topologies (diameters: 0.3 to 2.5 nm). Our approach involves subjecting confined water to annealing protocols to identify the temperature and time evolution of nonequilibrium glass kinetics. Furthermore, we complement this calorimetric approach with the dynamics of confined water, as seen by broadband dielectric spectroscopy and linear calorimetric measurements, including the fast scanning technique. This study demonstrated that confined water undergoes a glass transition in the temperature range of 170 to 200 K, depending on the confinement size and the interaction with the confinement walls. Moreover, we also show that the thermal event observed at ~136 K must be interpreted as an annealing prepeak, also referred to as the "shadow glass transition." Calorimetric measurements also allow the detection of a specific heat step above 200 K, which is insensitive to annealing and, thereby, interpreted as a true thermodynamic transition. Finally, by connecting our results to bulk water behavior, we offer a comprehensive understanding of confined water vitrification with potential implications for numerous applications.

AI-NERD: Elucidation of relaxation dynamics beyond equilibrium through AI-informed X-ray photon correlation spectroscopy

Understanding and interpreting dynamics of functional materials in situ is a grand challenge in physics and materials science due to the difficulty of experimentally probing materials at varied length and time scales. X-ray photon correlation spectroscopy (XPCS) is uniquely well-suited for characterizing materials dynamics over wide-ranging time scales. However, spatial and temporal heterogeneity in material behavior can make interpretation of experimental XPCS data difficult. In this work, we have developed an unsupervised deep learning (DL) framework for automated classification of relaxation dynamics from experimental data without requiring any prior physical knowledge of the system. We demonstrate how this method can be used to accelerate exploration of large datasets to identify samples of interest, and we apply this approach to directly correlate microscopic dynamics with macroscopic properties of a model system. Importantly, this DL framework is material and process agnostic, marking a concrete step towards autonomous materials discovery.

The laser pump X-ray probe system at LISA P08 PETRA III

Understanding and controlling the structure and function of liquid interfaces is a constant challenge in biology, nanoscience and nanotechnology, with applications ranging from molecular electronics to controlled drug release. X-ray reflectivity and grazing incidence diffraction provide invaluable probes for studying the atomic scale structure at liquid-air interfaces. The new time-resolved laser system at the LISA liquid diffractometer situated at beamline P08 at the PETRA III synchrotron radiation source in Hamburg provides a laser pump with X-ray probe. The femtosecond laser combined with the LISA diffractometer allows unique opportunities to investigate photo-induced structural changes at liquid interfaces on the pico- and nanosecond time scales with pump-probe techniques. A time resolution of 38 ps has been achieved and verified with Bi. First experiments include laser-induced effects on salt solutions and liquid mercury surfaces with static and varied time scales measurements showing the proof of concept for investigations at liquid surfaces.

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.

TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional

We investigate the performance of time-dependent density functional theory (TDDFT) for reproducing high-level reference x-ray absorption spectra of liquid water and water clusters. For this, we apply the integrated absolute difference (IAD) metric, previously used for x-ray emission spectra of liquid water [T. Fransson and L. G. M. Pettersson, J. Chem. Theory Comput. 19, 7333-7342 (2023)], in order to investigate which exchange-correlation (xc) functionals yield TDDFT spectra in best agreement to reference, as well as to investigate the suitability of IAD for x-ray absorption spectroscopy spectrum calculations. Considering highly asymmetric and symmetric six-molecule clusters, it is seen that long-range corrected xc-functionals are required to yield good agreement with the reference coupled cluster (CC) and algebraic-diagrammatic construction spectra, with 100% asymptotic Hartree-Fock exchange resulting in the lowest IADs. The xc-functionals with best agreement to reference have been adopted for larger water clusters, yielding results in line with recently published CC theory, but which still show some discrepancies in the relative intensity of the features compared to experiment.

Unraveling the dynamic slowdown in supercooled water: The role of dynamic disorder in jump motions

When a liquid is rapidly cooled below its melting point without inducing crystallization, its dynamics slow down significantly without noticeable structural changes. Elucidating the origin of this slowdown has been a long-standing challenge. Here, we report a theoretical investigation into the mechanism of the dynamic slowdown in supercooled water, a ubiquitous yet extraordinary substance characterized by various anomalous properties arising from local density fluctuations. Using molecular dynamics simulations, we found that the jump dynamics, which are elementary structural change processes, deviate from Poisson statistics with decreasing temperature. This deviation is attributed to slow variables competing with the jump motions, i.e., dynamic disorder. The present analysis of the dynamic disorder showed that the primary slow variable is the displacement of the fourth nearest oxygen atom of a jumping molecule, which occurs in an environment created by the fluctuations of molecules outside the first hydration shell. As the temperature decreases, the jump dynamics become slow and intermittent. These intermittent dynamics are attributed to the prolonged trapping of jumping molecules within extended and stable low-density domains. As the temperature continues to decrease, the number of slow variables increases due to the increased cooperative motions. Consequently, the jump dynamics proceed in a higher-dimensional space consisting of multiple slow variables, becoming slower and more intermittent. It is then conceivable that with further decreasing temperature, the slowing and intermittency of the jump dynamics intensify, eventually culminating in a glass transition.

In situ aggregation and early states of gelation of gold nanoparticle dispersions

The aggregation and onset of gelation of PEGylated gold nanoparticles dispersed in a glycerol-water mixture is studied by small-angle X-ray scattering and X-ray photon correlation spectroscopy. Tracking structural dynamics with sub-ms time resolution over a total experimental time of 8 hours corresponding to a time windows larger than 108 Brownian times and varying the temperature between 298 K and 266 K we can identify three regimes. First, while cooling to 275 K the particles show Brownian motion that slows down due to the increasing viscosity. Second, upon further cooling the static structure changes significantly, indicated by a broad structure factor peak. We attribute this to the formation of aggregates while the dynamics are still dominated by single-particle diffusion. Finally, the relaxation functions become more and more stretched accompanied by an increased slow down of the dynamics. At the same time the structure changes continuously indicating the onset of gelation. Our observations further suggest that the colloidal aggregation and gelation is characterized first by structural changes with a subsequent slowing down of the systems dynamics. The analysis also reveals that the details of the gelation process and the gel structure strongly depend on the thickness of the PEG-coating of the gold nanoparticles.

High-pressure X-ray photon correlation spectroscopy at fourth-generation synchrotron sources

A new experimental setup combining X-ray photon correlation spectroscopy (XPCS) in the hard X-ray regime and a high-pressure sample environment has been developed to monitor the pressure dependence of the internal motion of complex systems down to the atomic scale in the multi-gigapascal range, from room temperature to 600 K. The high flux of coherent high-energy X-rays at fourth-generation synchrotron sources solves the problems caused by the absorption of diamond anvil cells used to generate high pressure, enabling the measurement of the intermediate scattering function over six orders of magnitude in time, from 10-3 s to 103 s. The constraints posed by the high-pressure generation such as the preservation of X-ray coherence, as well as the sample, pressure and temperature stability, are discussed, and the feasibility of high-pressure XPCS is demonstrated through results obtained on metallic glasses.

Water Potential from Adaptive Force Matching for Ice and Liquid with Revised Dispersion Predicts Supercooled Liquid Anomalies in Good Agreement with Two Independent Experimental Fits

A revised version of the Water potential from Adaptive force matching for Ice and Liquid (WAIL) was developed by using the previous data set for fitting the WAIL model but with a dispersion term calculated using symmetry adapted perturbation theory (SAPT). The model has no adjustable parameters and relies solely on fitting first-principles information. The new model, named revised WAIL (rWAIL), shows improved predictions of most properties of water when compared to the previously published WAIL model. The rWAIL model also compares favorably to other first-principles-derived water models, such as MB-Pol, at only a fraction of the computational cost. The rWAIL model is used to study the properties of supercooled water. The model shows evidence of a liquid-liquid phase transition (LLPT) in the supercooled regimes with the liquid-liquid critical point (LLCP) at 203 K and 90 MPa. This estimate is in good agreement with a recent polynomial fit to the experimental density of water. Also, the fit to the surface tension of supercooled water based on the rWAIL model shows excellent agreement with the corresponding fit to the experimental data. Consistent with previously published molecular dynamics and experimental data, the surface tension of water exhibits exponential growth in the supercooled regime, which is likely a result of the emergence of a low-density liquid form of water. The simulation thus unites two separate experimental fits with one first-principles-based model, lending strong evidence of an LLPT in real water.

A continuum of amorphous ices between low-density and high-density amorphous ice

Amorphous ices are usually classified as belonging to low-density or high-density amorphous ice (LDA and HDA) with densities ρLDA ≈ 0.94 g/cm3 and ρHDA ≈ 1.15-1.17 g/cm3. However, a recent experiment crushing hexagonal ice (ball-milling) produced a medium-density amorphous ice (MDA, ρMDA ≈ 1.06 g/cm3) adding complexity to our understanding of amorphous ice and the phase diagram of supercooled water. Motivated by the discovery of MDA, we perform computer simulations where amorphous ices are produced by isobaric cooling and isothermal compression/decompression. Our results show that, depending on the pressure employed, isobaric cooling can generate a continuum of amorphous ices with densities that expand in between those of LDA and HDA (briefly, intermediate amorphous ices, IA). In particular, the IA generated at P ≈ 125 MPa has a remarkably similar density and average structure as MDA, implying that MDA is not unique. Using the potential energy landscape formalism, we provide an intuitive qualitative understanding of the nature of LDA, HDA, and the IA generated at different pressures. In this view, LDA and HDA occupy specific and well-separated regions of the PEL; the IA prepared at P = 125 MPa is located in the intermediate region of the PEL that separates LDA and HDA.

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.

Improvement of ultra-small-angle XPCS with the Extremely Brilliant Source

Recent technical developments and the performance of the X-ray photon correlation spectroscopy (XPCS) method over the ultra-small-angle range with the Extremely Brilliant Source (EBS) at the ESRF are described. With higher monochromatic coherent photon flux (∼1012 photons s-1) provided by the EBS and the availability of a fast pixel array detector (EIGER 500K detector operating at 23000 frames s-1), XPCS has become more competitive for probing faster dynamics in relatively dilute suspensions. One of the goals of the present development is to increase the user-friendliness of the method. This is achieved by means of a Python-based graphical user interface that enables online visualization and analysis of the processed data. The improved performance of XPCS on the Time-Resolved Ultra-Small-Angle X-ray Scattering instrument (ID02 beamline) is demonstrated using dilute model colloidal suspensions in several different applications.

Intrinsic Dynamics of Amorphous Ice Revealed by a Heterodyne Signal in X-ray Photon Correlation Spectroscopy Experiments

Unraveling the mechanism of water's glass transition and the interconnection between amorphous ices and liquid water plays an important role in our overall understanding of water. X-ray photon correlation spectroscopy (XPCS) experiments were conducted to study the dynamics and the complex interplay between the hypothesized glass transition in high-density amorphous ice (HDA) and the subsequent transition to low-density amorphous ice (LDA). Our XPCS experiments demonstrate that a heterodyne signal appears in the correlation function. Such a signal is known to originate from the interplay of a static component and a dynamic component. Quantitative analysis was performed on this heterodyne signal to extract the intrinsic dynamics of amorphous ice during the HDA-LDA transition. An angular dependence indicates non-isotropic, heterogeneous dynamics in the sample. Using the Stokes-Einstein relation to extract diffusion coefficients, the data are consistent with the scenario of static LDA islands floating within a diffusive matrix of high-density liquid water.

Lipid vesicle pools studied by passive X-ray microrheology

Vesicle pools can form by attractive interaction in a solution, mediated by proteins or divalent ions such as calcium. The pools, which are alternatively also denoted as vesicle clusters, form by liquid-liquid phase separation (LLPS) from an initially homogeneous solution. Due to the short range liquid-like order of vesicles in the pool or cluster, the vesicle-rich phase can also be regarded as a condensate, and one would like to better understand not only the structure of these systems, but also their dynamics. The diffusion of vesicles, in particular, is expected to change when vesicles are arrested in a pool. Here we investigate whether passive microrheology based on X-ray photon correlation spectroscopy (XPCS) is a suitable tool to study model systems of artificial lipid vesicles exhibiting LLPS, and more generally also other heterogeneous biomolecular fluids. We show that by adding highly scattering tracer particles to the solution, valuable information on the single vesicle as well as collective dynamics can be inferred. While the correlation functions reveal freely diffusing tracer particles in solutions at low CaCl[Formula: see text] concentrations, the relaxation rate [Formula: see text] shows a nonlinear dependence on [Formula: see text] at a higher concentration of around 8 mM CaCl[Formula: see text], characterised by two linear regimes with a broad cross-over. We explain this finding based on arrested diffusion in percolating vesicle clusters.

Calibrating TDDFT Calculations of the X-ray Emission Spectrum of Liquid Water: The Effects of Hartree-Fock Exchange

The structure and dynamics of liquid water continue to be debated, with insight provided by, among others, X-ray emission spectroscopy (XES), which shows a split in the high-energy 1b1 feature. This split is yet to be reproduced by theory, and it remains unclear if these difficulties are related to inaccuracies in dynamics simulations, spectrum calculations, or both. We investigate the performance of different methods for calculating XES of liquid water, focusing on the ability of time-dependent density functional theory (TDDFT) to reproduce reference spectra obtained by high-level coupled cluster and algebraic-diagrammatic construction scheme calculations. A metric for evaluating the agreement between theoretical spectra termed the integrated absolute difference (IAD), which considers the integral of shifted difference spectra, is introduced and used to investigate the performance of different exchange-correlation functionals. We find that computed spectra of symmetric and asymmetric model water structures are strongly and differently influenced by the amount of Hartree-Fock exchange, with best agreement to reference spectra for ∼40-50%. Lower percentages tend to yield high density of contributing states, resulting in too broad features. The method introduced here is useful also for other spectrum calculations, in particular where the performance for ensembles of structures are evaluated.

Challenges of Electron Correlation Microscopy on Amorphous Silicon and Amorphous Germanium

Electron correlation microscopy experiments were conducted on amorphous germanium (a-Ge) and amorphous silicon (a-Si) with the goal to study self-diffusion. For this purpose, a series of tilted dark-field images were acquired during in situ heating of the samples in a transmission electron microscope. These experiments show that the measurements are greatly affected by artefacts. Contamination, crystallization, electron beam-induced sputtering, and macroscopic bending of the samples pose major obstacles to the measurements. Other, more subtle experimental artefacts could occur in addition to these which makes interpretations regarding the structural dynamics nearly impossible. The data were nonetheless evaluated to see if some useful information could be extracted. One such result is that the distribution of the characteristic times τKWW, which were obtained from stretched exponential fits to the intensity autocorrelation data, is spatially heterogeneous. This spatial heterogeneity is assumed to be caused by a potential nonergodicity of the materials, the artefacts or an inhomogeneous amorphous structure. Further data processing shows that the characteristic times τKWW are moreover temperature independent, especially for the a-Ge data. It is concluded that the structural rearrangements over time are primarily electron beam-driven and that diffusive dynamics are too slow to be measured at the chosen, experimentally accessible annealing temperatures.

Pressure-induced nonmonotonic cross-over of steady relaxation dynamics in a metallic glass

Relaxation dynamics, as a key to understand glass formation and glassy properties, remains an elusive and challenging issue in condensed matter physics. In this work, in situ high-pressure synchrotron high-energy X-ray photon correlation spectroscopy has been developed to probe the atomic-scale relaxation dynamics of a cerium-based metallic glass during compression. Although the sample density continuously increases, the collective atomic motion initially slows down as generally expected and then counterintuitively accelerates with further compression (density increase), showing an unusual nonmonotonic pressure-induced steady relaxation dynamics cross-over at ~3 GPa. Furthermore, by combining in situ high-pressure synchrotron X-ray diffraction, the relaxation dynamics anomaly is evidenced to closely correlate with the dramatic changes in local atomic structures during compression, rather than monotonically scaling with either sample density or overall stress level. These findings could provide insight into relaxation dynamics and their relationship with local atomic structures of glasses.

Coherent X-ray Scattering Reveals Nanoscale Fluctuations in Hydrated Proteins

Hydrated proteins undergo a transition in the deeply supercooled regime, which is attributed to rapid changes in hydration water and protein structural dynamics. Here, we investigate the nanoscale stress-relaxation in hydrated lysozyme proteins stimulated and probed by X-ray Photon Correlation Spectroscopy (XPCS). This approach allows us to access the nanoscale dynamics in the deeply supercooled regime (T = 180 K), which is typically not accessible through equilibrium methods. The observed stimulated dynamic response is attributed to collective stress-relaxation as the system transitions from a jammed granular state to an elastically driven regime. The relaxation time constants exhibit Arrhenius temperature dependence upon cooling with a minimum in the Kohlrausch-Williams-Watts exponent at T = 227 K. The observed minimum is attributed to an increase in dynamical heterogeneity, which coincides with enhanced fluctuations observed in the two-time correlation functions and a maximum in the dynamic susceptibility quantified by the normalized variance χ_T. The amplification of fluctuations is consistent with previous studies of hydrated proteins, which indicate the key role of density and enthalpy fluctuations in hydration water. Our study provides new insights into X-ray stimulated stress-relaxation and the underlying mechanisms behind spatiotemporal fluctuations in biological granular materials.

Dynamic Heterogeneity and Kovacs' Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. We used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water's structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

Temperature Dependence of Hydrogen Bond Networks of Liquid Water: Thermodynamic Properties and Structural Heterogeneity from Topological Descriptors

Topological analyses of hydrogen bond networks were performed based on the complex network and island statistics of liquid water at different temperatures. The influence of temperature on the liquid water structures and the topological properties of the hydrogen bond networks was investigated by Metropolis Monte Carlo simulations with the TIP4P/2005 potential model. The bilinear behavior of the second peak in the radial distribution function with the temperature was properly reproduced by these simulations. The average connectivity also displayed a bilinear behavior consistent with being a local descriptor. The semiglobal average path length (or geodesic distance) descriptor showed an unprecedented trimodal distribution, whose areas were dependent on the temperature. Considering equilibrium between these three sets of networks, standard enthalpy and entropy of equilibrium were determined for the first time, providing new insights into the structural heterogeneities of liquid water with interesting perspectives for modeling these quantitative properties of hydrogen bond networks.

Medium-density amorphous ice

Amorphous ices govern a range of cosmological processes and are potentially key materials for explaining the anomalies of liquid water. A substantial density gap between low-density and high-density amorphous ice with liquid water in the middle is a cornerstone of our current understanding of water. However, we show that ball milling "ordinary" ice Ih at low temperature gives a structurally distinct medium-density amorphous ice (MDA) within this density gap. These results raise the possibility that MDA is the true glassy state of liquid water or alternatively a heavily sheared crystalline state. Notably, the compression of MDA at low temperature leads to a sharp increase of its recrystallization enthalpy, highlighting that H2O can be a high-energy geophysical material.

Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice

Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale < 100 ns and widely separated from ice formation, which begins at times >1 μs.

Single and multi-pulse based X-ray photon correlation spectroscopy

The ability of pulsed nature of synchrotron radiation opens up the possibility of studying microsecond dynamics in complex materials via speckle-based techniques. Here, we present the study of measuring the dynamics of a colloidal system by combining single and multiple X-ray pulses of a storage ring. In addition, we apply speckle correlation techniques at various pulse patterns to collect correlation functions from nanoseconds to milliseconds. The obtained sample dynamics from all correlation techniques at different pulse patterns are in very good agreement with the expected dynamics of Brownian motions of silica nanoparticles in water. Our study will pave the way for future pulsed X-ray investigations at various synchrotron X-ray sources using individual X-ray pulse patterns.

Free energy calculations and unbiased molecular dynamics targeting the liquid-liquid transition in water no man's land

The existence of a first-order phase transition between a low-density liquid (LDL) and a high-density liquid (HDL) form of supercooled water has been a central and highly debated issue of physics and chemistry for the last three decades. We present a computational study that allows us to determine the free-energy landscapes of supercooled water over a wide range of pressure and temperature conditions using the TIP4P/2005 force field. Our approach combines topology-based structural transformation coordinates, state-of-the-art free-energy calculation methods, and extensive unbiased molecular dynamics. All our diverse simulations cannot detect any barrier within the investigated timescales and system size, for a discontinuous transition between the LDL and HDL forms throughout the so-called "no man's land," until the onset of the solid, non-diffusive amorphous forms.

Water adsorption and dynamics on graphene and other 2D materials: Computational and experimental advances

The interaction of water and surfaces, at molecular level, is of critical importance for understanding processes such as corrosion, friction, catalysis and mass transport. The significant literature on interactions with single crystal metal surfaces should not obscure unknowns in the unique behaviour of ice and the complex relationships between adsorption, diffusion and long-range inter-molecular interactions. Even less is known about the atomic-scale behaviour of water on novel, non-metallic interfaces, in particular on graphene and other 2D materials. In this manuscript, we review recent progress in the characterisation of water adsorption on 2D materials, with a focus on the nano-material graphene and graphitic nanostructures; materials which are of paramount importance for separation technologies, electrochemistry and catalysis, to name a few. The adsorption of water on graphene has also become one of the benchmark systems for modern computational methods, in particular dispersion-corrected density functional theory (DFT). We then review recent experimental and theoretical advances in studying the single-molecular motion of water at surfaces, with a special emphasis on scattering approaches as they allow an unparalleled window of observation to water surface motion, including diffusion, vibration and self-assembly.

Using coherent X-rays to follow dynamics in amorphous ices

Amorphous solid water plays an important role in our overall understanding of water's phase diagram. X-ray scattering is an important tool for characterising the different states of water, and modern storage ring and XFEL facilities have opened up new pathways to simultaneously study structure and dynamics. Here, X-ray photon correlation spectroscopy (XPCS) was used to study the dynamics of high-density amorphous (HDA) ice upon heating. We follow the structural transition from HDA to low-density amorphous (LDA) ice, by using wide-angle X-ray scattering (WAXS), for different heating rates. We used a new type of sample preparation, which allowed us to study μm-sized ice layers rather than powdered bulk samples. The study focuses on the non-equilibrium dynamics during fast heating, spontaneous transformation and crystallization. Performing the XPCS study at ultra-small angle (USAXS) geometry allows us to characterize the transition dynamics at length scales ranging from 60 nm-800 nm. For the HDA-LDA transition we observe a clear separation in three dynamical regimes, which show different dynamical crossovers at different length scales. The crystallization from LDA, instead, is observed to appear homogenously throughout the studied length scales.

Cryomicroscopy in situ: what is the smallest molecule that can be directly identified without labels in a cell?

Electron cryomicroscopy (cryoEM) has made great strides in the last decade, such that the atomic structure of most biological macromolecules can, at least in principle, be determined. Major technological advances - in electron imaging hardware, data analysis software, and cryogenic specimen preparation technology - continue at pace and contribute to the exponential growth in the number of atomic structures determined by cryoEM. It is now conceivable that within the next decade we will have structures for hundreds of thousands of unique protein and nucleic acid molecular complexes. But the answers to many important questions in biology would become obvious if we could identify these structures precisely inside cells with quantifiable error. In the context of an abundance of known structures, it is appropriate to consider the current state of electron cryomicroscopy for frozen specimens prepared directly from cells, and try to answer to the question of the title, both now and in the foreseeable future.

Multiple insights call for revision of modern thermodynamic models to account for structural fluctuations in water

Modern thermodynamic models incorporate the concept of association (hydrogen bonding) and they can describe very satisfactorily many properties of water containing mixtures. They have not been successful in representing water's anomalous properties and this work provides a possible explanation. We have analyzed and interpreted recent experimental data, molecular simulation results, and two-state theory approaches and compared against the predictions from thermodynamic models. We show that the dominance of the tetrahedral structure implemented in modern thermodynamic models may be the reason for their failure for describing water systems. While this study does not prove the two-state theories for water, it indicates that a high level of tetrahedral structure of water is not in agreement with water's anomalous properties when used in thermodynamic models.

Infrared Spectroscopy on Equilibrated High-Density Amorphous Ice

High-density (HDA) and low-density amorphous ices (LDA) are believed to be counterparts of the high- and low-density liquid phases of water, respectively. In order to better understand how the vibrational modes change during the transition between the two solid states, we present infrared spectroscopy measurements, following the change of the decoupled OD-stretch (vOD) (∼2460 cm-1) and OH-combinational mode (vOH + v2, vOH + 2vR) (∼5000 cm-1). We observe a redshift from HDA to LDA, accompanied with a drastic decrease of the bandwidth. The hydrogen bonds are stronger in LDA, which is caused by a change in the coordination number and number of water molecules interstitial between the first and second hydration shell. The unusually broad uncoupled OD band also clearly distinguishes HDA from other crystalline high-pressure phases, while the shape and position of the in situ prepared LDA are comparable to those of vapor-deposited amorphous ice.

pyXPCSviewer: an open-source interactive tool for X-ray photon correlation spectroscopy visualization and analysis

pyXPCSviewer, a Python-based graphical user interface that is deployed at beamline 8-ID-I of the Advanced Photon Source for interactive visualization of XPCS results, is introduced. pyXPCSviewer parses rich X-ray photon correlation spectroscopy (XPCS) results into independent PyQt widgets that are both interactive and easy to maintain. pyXPCSviewer is open-source and is open to customization by the XPCS community for ingestion of diversified data structures and inclusion of novel XPCS techniques, both of which are growing demands particularly with the dawn of near-diffraction-limited synchrotron sources and their dedicated XPCS beamlines.

Structure of Sulfuric Acid Solutions Using Pair Distribution Function Analysis

Solvation and mesoscale ordering of sulfuric acid and other strong acid solutions leads to suppressed freezing points and strong rheological changes with varying concentration. While the solid-state structures are well-understood, studies focused on the evolving solvation structure in the solution phase have probed a limited concentration range (∼1-6 M). This study applies a total scattering approach in both the wide-angle X-ray scattering (WAXS) and pair distribution function (PDF) regimes to elucidate the evolving solvation structure over its full range of acid concentration (0-18 M). The emergence of a prepeak in the WAXS regime at intermediate concentrations indicates a transition from noninteracting sulfate molecules in the dilute limit to sterically limited sulfates at concentrations near its deep eutectic point. Fits to the PDF data quantify this trend, showing a transition from octahedrally hydrated sulfates up to 6-7 M concentrations, followed by gradual dehydration, and eventually reaching a solution structure similar to that of water-in-salt electrolyte systems at high acid concentrations.

Advancing Chemical Separations: Unraveling the Structure and Dynamics of Phase Splitting in Liquid-Liquid Extraction

Liquid-liquid extraction (LLE), the go-to process for a variety of chemical separations, is limited by spontaneous organic phase splitting upon sufficient solute loading, called third phase formation. In this study we explore the applicability of critical phenomena theory to gain insight into this deleterious phase behavior with the goal of improving separations efficiency and minimizing waste. A series of samples representative of rare earth purification were constructed to include each of one light and one heavy lanthanide (cerium and lutetium) paired with one of two common malonamide extractants (DMDOHEMA and DMDBTDMA). The resulting postextraction organic phases are chemically complex and often form rich hierarchical structures whose statics and dynamics near the critical point were probed herein with small-angle X-ray scattering and high-speed X-ray photon correlation spectroscopy. Despite their different extraction behaviors, all samples show remarkably similar critical behavior with exponents well described by classical critical point theory consistent with the 3D Ising model, where the critical behavior is characterized by fluctuations with a single diverging length scale. This unexpected result indicates a significant reduction in relevant chemical parameters at the critical point, indicating that the underlying behavior of phase transitions in LLE rely on far fewer variables than are generally assumed. The obtained scalar order parameter is attributed to the extractant fraction of the extractant/diluent mixture, revealing that other solution components and their respective concentrations simply shift the critical temperature but do not affect the nature of the critical fluctuations. These findings point to an opportunity to drastically simplify studies of liquid-liquid phase separation and phase diagram development in general while providing insights into LLE process improvement.

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.

Demonstration of 3D photon correlation spectroscopy in the hard X-ray regime

Three-dimensional photon correlation spectroscopy (3D PCS) is a well-known technique developed to suppress multiple scattering contributions in correlation functions, which are inevitably involved when an optical laser is employed to investigate dynamics in a turbid system. Here, we demonstrate a proof-of-principle study of 3D PCS in the hard X-ray regime. We employ an X-ray optical cross-correlator to measure the dynamics of silica colloidal nanoparticles dispersed in polypropylene glycol. The obtained cross correlation functions show very good agreement with auto-correlation measurements. This demonstration provides the foundation for X-ray speckle-based studies of very densely packed soft matter systems.

Long-Range Structures of Amorphous Solid Water

High-energy X-ray diffraction (XRD) and Fourier transform infrared spectroscopy (FTIR) of amorphous solid water (ASW) were studied during vapor deposition and the heating process. From the diffraction patterns, the oxygen–oxygen pair distribution functions (PDFs) were calculated up to the eighth coordination shell and an r = 23 Å. The PDF of ASW obtained both during vapor deposition at 80 K as well as the subsequent heating are consistent with that of low-density amorphous ice. The formation and temperature-induced collapse of micropores were observed in the XRD data and in the FTIR measurements, more specifically, in the OH stretch and the dangling mode. Above 140 K, ASW crystallizes into a stacking disordered ice, I_sd. It is observed that the fourth, fifth, and sixth peaks in the PDF, corresponding to structural arrangements between 8 and 12 Å, are the most sensitive to the onset of crystallization.

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.

High-Pressure Nonequilibrium Dynamics on Second-to-Microsecond Time Scales: Application of Time-Resolved X-ray Diffraction and Dynamic Compression in Ice

The study of nonequilibrium transition dynamics on structural transformation from the second to microsecond regime, a time scale between static and shock compression, is an emerging field of high-pressure research. There are ample opportunities to uncover novel physical phenomena within this time regime. Herein, we briefly review the development and application of a dynamic compression technique based on a diamond anvil cell (DAC) and time-resolved X-ray diffraction (TRXRD) for the study of time-, pressure-, and temperature-dependent structural dynamics. Applications of the techniques are illustrated with our recent investigations on the mechanisms of the interconversions between different high-pressure ice polymorphs. These examples demonstrate that a combination of dynamic compression and TRXRD is a versatile approach capable of providing information on the kinetics and thermodynamic nature associated with structural transformations. Future improvement of rapid compression and TRXRD techniques and potentially interesting research topics in this area are suggested.

Molecular Properties and Chemical Transformations Near Interfaces

The properties of bulk water and aqueous solutions are known to change in the vicinity of an interface and/or in a confined environment, including the thermodynamics of ion selectivity at interfaces, transition states and pathways of chemical reactions, and nucleation events and phase growth. Here we describe joint progress in identifying unifying concepts about how air, liquid, and solid interfaces can alter molecular properties and chemical reactivity compared to bulk water and multicomponent solutions. We also discuss progress made in interfacial chemistry through advancements in new theory, molecular simulation, and experiments.

Interplay between Kinetics and Dynamics of Liquid-Liquid Phase Separation in a Protein Solution Revealed by Coherent X-ray Spectroscopy

Microscopic dynamics of complex fluids in the early stage of spinodal decomposition (SD) is strongly intertwined with the kinetics of structural evolution, which makes a quantitative characterization challenging. In this work, we use X-ray photon correlation spectroscopy to study the dynamics and kinetics of a protein solution undergoing liquid-liquid phase separation (LLPS). We demonstrate that in the early stage of SD, the kinetics relaxation is up to 40 times slower than the dynamics and thus can be decoupled. The microscopic dynamics can be well described by hyper-diffusive ballistic motions with a relaxation time exponentially growing with time in the early stage followed by a power-law increase with fluctuations. These experimental results are further supported by simulations based on the Cahn-Hilliard equation. The established framework is applicable to other condensed matter and biological systems undergoing phase transitions and may also inspire further theoretical work.

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.

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.