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.

Enhancement and maximum in the isobaric specific-heat capacity measurements of deeply supercooled water using ultrafast calorimetry.

Knowledge of the temperature dependence of the isobaric specific heat (Cp) upon deep supercooling can give insights regarding the anomalous properties of water. If a maximum in Cp exists at a specific temperature, as in the isothermal compressibility, it would further validate the liquid-liquid critical point model that can explain the anomalous increase in thermodynamic response functions. The challenge is that the relevant temperature range falls in the region where ice crystallization becomes rapid, which has previously excluded experiments. Here, we have utilized a methodology of ultrafast calorimetry by determining the temperature jump from femtosecond X-ray pulses after heating with an infrared laser pulse and with a sufficiently long time delay between the pulses to allow measurements at constant pressure. Evaporative cooling of ∼15-µm diameter droplets in vacuum enabled us to reach a temperature down to ∼228 K with a small fraction of the droplets remaining unfrozen. We observed a sharp increase in Cp, from 88 J/mol/K at 244 K to about 218 J/mol/K at 229 K where a maximum is seen. The Cp maximum is at a similar temperature as the maxima of the isothermal compressibility and correlation length. From the Cp measurement, we estimated the excess entropy and self-diffusion coefficient of water and these properties decrease rapidly below 235 K.

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 (q2). We found that q2 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 q2 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.

Anisotropic X-Ray Scattering of Transiently Oriented Water.

We study the structural dynamics of liquid water by time-resolved anisotropic x-ray scattering under the optical Kerr effect condition. In this way, we can separate the anisotropic scattering decay of 160 fs from the delayed temperature increase of ∼0.1 K occurring at 1 ps and quantify transient changes in the O-O pair distribution function. Polarizable molecular dynamics simulations reproduce well the experiment, indicating transient alignment of molecules along the electric field, which shortens the nearest-neighbor distances. In addition, analysis of the simulated water local structure provides evidence that two hypothesized fluctuating water configurations exhibit different polarizability.

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

Water exists in high- and low-density amorphous ice forms (HDA and LDA), which could correspond to the glassy states of high- (HDL) and low-density liquid (LDL) in the metastable part of the phase diagram. However, the nature of both the glass transition and the high-to-low-density transition are debated and new experimental evidence is needed. Here we combine wide-angle X-ray scattering (WAXS) with X-ray photon-correlation spectroscopy (XPCS) in the small-angle X-ray scattering (SAXS) geometry to probe both the structural and dynamical properties during the high-to-low-density transition in amorphous ice at 1 bar. By analyzing the structure factor and the radial distribution function, the coexistence of two structurally distinct domains is observed at T = 125 K. XPCS probes the dynamics in momentum space, which in the SAXS geometry reflects structural relaxation on the nanometer length scale. The dynamics of HDA are characterized by a slow component with a large time constant, arising from viscoelastic relaxation and stress release from nanometer-sized heterogeneities. Above 110 K a faster, strongly temperature-dependent component appears, with momentum transfer dependence pointing toward nanoscale diffusion. This dynamical component slows down after transition into the low-density form at 130 K, but remains diffusive. The diffusive character of both the high- and low-density forms is discussed among different interpretations and the results are most consistent with the hypothesis of a liquid-liquid transition in the ultraviscous regime.

X-ray studies of the transformation from high- to low-density amorphous water.

Here we report about the structural evolution during the conversion from high-density amorphous ices at ambient pressure to the low-density state. Using high-energy X-ray diffraction, we have monitored the transformation by following in reciprocal space the structure factor SOO( Q) and derived in real space the pair distribution function gOO( r). Heating equilibrated high-density amorphous ice (eHDA) at a fast rate (4 K min-1), the transition to the low-density form occurs very rapidly, while domains of both high- and low-density coexist. On the other hand, the transition in the case of unannealed HDA (uHDA) and very-high-density amorphous ice is more complex and of continuous nature. The direct comparison of eHDA and uHDA indicates that the molecular structure of uHDA contains a larger amount of tetrahedral motives. The different crystallization behaviour of the derived low-density amorphous states is interpreted as emanating from increased tetrahedral coordination present in uHDA. This article is part of the theme issue 'The physics and chemistry of ice: scaffolding across scales, from the viability of life to the formation of planets'.

Radial distribution functions of water: Models vs experiments.

We study the temperature behavior of the first four peaks of the oxygen-oxygen radial distribution function of water, simulated by the TIP4P/2005, MB-pol, TIP5P, and SPC/E models and compare to experimental X-ray diffraction data, including a new measurement which extends down to 235 K [H. Pathak et al., J. Chem. Phys. 150, 224506 (2019)]. We find the overall best agreement using the MB-pol and TIP4P/2005 models. We observe, upon cooling, a minimum in the position of the second shell simulated with TIP4P/2005 and SPC/E potentials, located close to the temperature of maximum density. We also calculated the two-body entropy and the contributions coming from the first, second, and outer shells to this quantity. We show that, even if the main contribution comes from the first shell, the contribution of the second shell can become important at low temperature. While real water appears to be less ordered at short distance than obtained by any of the potentials, the different water potentials show more or less order compared to the experiments depending on the considered length-scale.

Intermediate range O-O correlations in supercooled water down to 235 K.

Wide angle x-ray scattering of supercooled water down to 234.8 K was studied using high energy x rays at the European Synchrotron Radiation Facility. The oxygen-oxygen pair distribution function (PDF) was calculated from the scattering pattern out to the 5th peak at an intermolecular distance, r ≈ 11 Å. We observe that the 4th peak and the 5th peak in the PDF increase in height upon supercooling. We also observe that the 4th peak position (r4) shifts to shorter distances upon supercooling consistent with previous studies, but we see a more rapid change at the lowest temperature. The running oxygen-oxygen coordination number is calculated for 5 different temperatures, and an isosbestic point at riso = 3.31 ± 0.05 Å was found corresponding to a coordination number of 4.39 ± 0.15. The comparison of the PDF of the coldest water with that of amorphous ice shows distinct differences. We propose that there are 5-member pentamer rings in low density liquid-like structures giving rise to the sharp correlations at r ≈ 9 Å and r ≈ 11 Å.

Apparent power-law behavior of water's isothermal compressibility and correlation length upon supercooling.

The isothermal compressibility and correlation length of supercooled water obtained from small-angle X-ray scattering (SAXS) were analyzed by fits based on an apparent power-law in the temperature range from 280 K down to the temperature of maximum compressibility at 229 K. Although the increase in thermodynamic response functions is not towards a critical point, it is still possible to obtain an apparent power law all the way to the maximum values with best-fit exponents of γ = 0.40 ± 0.01 for the isothermal compressibility and ν = 0.26 ± 0.03 for the correlation length. The ratio between these exponents is close to a value of ≈0.5, as expected for a critical point, indicating the proximity of a potential second critical point. Comparison of γ obtained from experiment with molecular dynamics simulations on the iAMOEBA water model shows that it would be located at pressures in the neighborhood of 1 kbar. The high value and sharpness of the compressibility maximum observed in the experiment are not reproduced by any of the existing classical water models, thus inviting further development of simulation models of water.

Temperature-Independent Nuclear Quantum Effects on the Structure of Water.

Nuclear quantum effects (NQEs) have a significant influence on the hydrogen bonds in water and aqueous solutions and have thus been the topic of extensive studies. However, the microscopic origin and the corresponding temperature dependence of NQEs have been elusive and still remain the subject of ongoing discussion. Previous x-ray scattering investigations indicate that NQEs on the structure of water exhibit significant temperature dependence [Phys. Rev. Lett. 94, 047801 (2005)PRLTAO0031-900710.1103/PhysRevLett.94.047801]. Here, by performing wide-angle x-ray scattering of H_{2}O and D_{2}O droplets at temperatures from 275 K down to 240 K, we determine the temperature dependence of NQEs on the structure of water down to the deeply supercooled regime. The data reveal that the magnitude of NQEs on the structure of water is temperature independent, as the structure factor of D_{2}O is similar to H_{2}O if the temperature is shifted by a constant 5 K, valid from ambient conditions to the deeply supercooled regime. Analysis of the accelerated growth of tetrahedral structures in supercooled H_{2}O and D_{2}O also shows similar behavior with a clear 5 K shift. The results indicate a constant compensation between NQEs delocalizing the proton in the librational motion away from the bond and in the OH stretch vibrational modes along the bond. This is consistent with the fact that only the vibrational ground state is populated at ambient and supercooled conditions.

Co-condensation of nonane and D2O in a supersonic nozzle.

We study the unary and binary nucleation and growth of nonane-D2O nanodroplets in a supersonic nozzle. Fourier Transform Infrared spectroscopy measurements provide the overall composition of the droplets and Small Angle X-ray Scattering experiments measure the size and number density of the droplets. The unary nucleation rates Jmax of nonane, 9.4 × 10(15) < Jmax /cm(-3) s(-1) < 2.0 × 10(16), and those of D2O, 2.4 × 10(17) < Jmax /cm(-3) s(-1) < 4.1 × 10(17), measured here agree well with previous results. In most of the binary condensation experiments new particle formation is dominated by D2O, but the observed nucleation rates are decreased by up to a factor of 6 relative to the rates measured for pure D2O, an effect that can be partly explained by non-isothermal nucleation theory. The subsequent condensation of D2O is inhibited both by the increased temperature of the binary droplets relative to the pure D2O droplets, and because the binary droplet surface is expected to be comprised largely of nonane. For the one case where nonane appears to initiate condensation, we find that the nucleation rate is about 50% higher than that observed for pure nonane at comparable pv0, consistent with significant particle formation driven by D2O.

The structure of D2O-nonane nanodroplets.

We study the internal structure of nanometer-sized D2O-nonane aerosol droplets formed in supersonic nozzle expansions using a variety of experimental techniques including small angle X-ray scattering (SAXS). By fitting the SAXS spectra to a wide range of droplet structure models, we find that the experimental results are inconsistent with mixed droplets that form aqueous core-organic shell structures, but are quite consistent with spherically asymmetric lens-on-sphere structures. The structure that agrees best with the SAXS data and Fourier transform infra-red spectroscopy measurements is that of a nonane lens on a sphere of D2O with a contact angle in the range of 40°-120°.

Experimental evidence for surface freezing in supercooled n-alkane nanodroplets.

Intermediate chain length (16 ≤i≤ 50) n-alkanes are known to surface freeze at temperatures that are up to three degrees higher than the equilibrium melting point [B. M. Ocko et al., Phys. Rev. E, 1997, 55, 3164-3182]. Our recent experimental results suggest that highly supercooled nanodroplets of n-octane and n-nonane also surface freeze, and subsequently bulk crystallization occurs. The data yield surface and bulk nucleation rates on the order of ~10(15) cm(-2) s(-1) and ~10(22) cm(-3) s(-1), respectively, at temperatures between 180 K and 200 K. Molecular dynamics simulations at the united atom level were used to follow the freezing of a supercooled n-octane drop and show that an ordered monolayer develops on the surface of the droplet almost immediately, and the rest of the droplet then freezes in a layer-by-layer manner.

Freezing of heavy water (D2O) nanodroplets.

We follow the freezing of heavy water (D2O) nanodroplets formed in a supersonic nozzle apparatus using position resolved pressure trace measurements, Fourier transform infrared spectroscopy, and small-angle X-ray scattering. For these 3-9 nm radii droplets, freezing starts between 223 and 225 K, at volume based ice nucleation rates Jice,V on the order of 10(23) cm(-3) s(-1) or surface based ice nucleation rates Jice,S on the order of 10(16) cm(-2) s(-1). The temperatures corresponding to the onset of D2O ice nucleation are higher than those reported for H2O by Manka et al. [Manka, A.; Pathak, H.; Tanimura, S.; Wölk, J.; Strey, R.; Wyslouzil, B. E. Phys. Chem. Chem. Phys.2012, 14, 4505]. Although the values of Jice,S scale somewhat better with droplet size than values of Jice,V, the data are not accurate enough to state that nucleation is surface initiated. Finally, using current estimates of the thermophysical properties of D2O and the theoretical framework presented by Murray et al. [Murray, B. J.; Broadley, S. L.; Wilson, T. W.; Bull, S. J.; Wills, R. H.; Christenson, H. K.; Murray, E. J. Phys. Chem. Chem. Phys.2010, 12, 10380], we find that the theoretical ice nucleation rates are within 3 orders of magnitude of the measured rates over an ∼15 K temperature range.

Binary nucleation rates for ethanol/water mixtures in supersonic Laval nozzles: analyses by the first and second nucleation theorems.

We performed pressure trace measurements and small angle x-ray scattering measurements to determine the vapor-liquid nucleation rates of EtOH/H2O mixtures including pure EtOH and pure H2O in two supersonic Laval nozzles with different expansion rates. The nucleation rates varied from 0.9 × 10(17) to 16 × 10(17) cm(-3) s(-1) over the temperature range of 210 K to 230 K, EtOH activity range of 0 to 11.6, and H2O activity range of 0 to 124. The first and second nucleation theorems were applied to the nucleation rates to estimate the sizes, compositions, and excess energies of the critical clusters. The critical clusters contained from 4 to 15 molecules for pure H2O and EtOH/H2O clusters, and from 16 to 23 molecules for pure EtOH clusters. Comparing the excess energies of the pure H2O critical clusters with the results of a quantum-chemistry calculation suggested that the pre-factor of the theoretical nucleation rate is almost constant regardless of the monomer concentration. One possible explanation for this result is that cooling of the critical clusters limits the nucleation rate under the highly supersaturated conditions. The results of the analyses also yielded the relation between the surface energy and the composition of the critical clusters, where the latter are predicted to consist only of surface molecules. Applying this relationship to the EtOH/H2O bulk liquid mixtures, we estimated the EtOH mole fraction in the surface layer and found it is higher than that derived from the surface tension based on the Gibbs adsorption equation when the EtOH mole fraction in the liquid is higher than about 0.2 mol/mol. This discrepancy was attributed to the existence of the EtOH depletion layer just below the surface layer of the liquid.

Melting domain size and recrystallization dynamics of ice revealed by time-resolved x-ray scattering.

The phase transition between water and ice is ubiquitous and one of the most important phenomena in nature. Here, we performed time-resolved x-ray scattering experiments capturing the melting and recrystallization dynamics of ice. The ultrafast heating of ice I is induced by an IR laser pulse and probed with an intense x-ray pulse which provided us with direct structural information on different length scales. From the wide-angle x-ray scattering (WAXS) patterns, the molten fraction, as well as the corresponding temperature at each delay, were determined. The small-angle x-ray scattering (SAXS) patterns, together with the information extracted from the WAXS analysis, provided the time-dependent change of the size and the number of liquid domains. The results show partial melting (~13%) and superheating of ice occurring at around 20 ns. After 100 ns, the average size of the liquid domains grows from about 2.5 nm to 4.5 nm by the coalescence of approximately six adjacent domains. Subsequently, we capture the recrystallization of the liquid domains, which occurs on microsecond timescales due to the cooling by heat dissipation and results to a decrease of the average liquid domain size.

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.

Freezing water in no-man's land.

We report homogeneous ice nucleation rates between 202 K and 215 K, thereby reducing the measurement gap that previously existed between 203 K and 228 K. These temperatures are significantly below the homogenous freezing limit, T(H)≈ 235 K for bulk water, and well within no-man's land. The ice nucleation rates are determined by characterizing nanodroplets with radii between 3.2 and 5.8 nm produced in a supersonic nozzle using three techniques: (1) pressure trace measurements to determine the properties of the flow as well as the temperature and velocity of the droplets, (2) small angle X-ray scattering (SAXS) to measure the size and number density of the droplets, and (3) Fourier Transform Infrared (FTIR) spectroscopy to follow the liquid to solid phase transition. Assuming that nucleation occurs throughout the droplet volume, the measured ice nucleation rates J(ice,V) are on the order of 10(23) cm(-3) s(-1), and agree well with published values near 203 K.

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 (Isd), 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 (Ih) was also observed, suggesting that within our timeframe, Isd starts annealing into Ih. 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.

Experimental observation of the liquid-liquid transition in bulk supercooled water under pressure.

We prepared bulk samples of supercooled liquid water under pressure by isochoric heating of high-density amorphous ice to temperatures of 205 ± 10 kelvin, using an infrared femtosecond laser. Because the sample density is preserved during the ultrafast heating, we could estimate an initial internal pressure of 2.5 to 3.5 kilobar in the high-density liquid phase. After heating, the sample expanded rapidly, and we captured the resulting decompression process with femtosecond x-ray laser pulses at different pump-probe delay times. A discontinuous structural change occurred in which low-density liquid domains appeared and grew on time scales between 20 nanoseconds to 3 microseconds, whereas crystallization occurs on time scales of 3 to 50 microseconds. The dynamics of the two processes being separated by more than one order of magnitude provides support for a liquid-liquid transition in bulk supercooled water.