Homogeneous Nucleation of Ice in Transiently-Heated, Supercooled Liquid Water Films

Reversal of crystallization in cryoprotected samples by laser editing

Advances in cryobiology techniques commonly target either the cooling or the warming cycle, while little thought has been given to ≪repair≫ protocols applicable during cold storage. In particular, crystallization is the dominant threat to cryopreserved samples but proceeds from small nuclei that are innocuous if further growth is forestalled. To this end, we propose a laser editing technique that locally heats individual crystals above their melting point by a focused nanosecond pulse, followed by amorphization during rapid resolidification. As a reference, we first apply the approach to ice crystals in cryoprotected solution and use Raman confocal mapping to study the deactivation of crystalline order. Then, we examine dimethyl sulfoxide trihydrate crystals that can germinate at low temperatures in maximally freeze concentrated regions, as commonly produced by equilibrium cooling protocols. We show how to uniquely identify this phase from Raman spectra and evidence retarded growth of laser-edited crystals during warming.

Liquid-liquid transition and ice crystallization in a machine-learned coarse-grained water model

Mounting experimental evidence supports the existence of a liquid-liquid transition (LLT) in high-pressure supercooled water. However, fast crystallization of supercooled water has impeded identification of the LLT line TLL(p) in experiments. While the most accurate all-atom (AA) water models display a LLT, their computational cost limits investigations of its interplay with ice formation. Coarse-grained (CG) models provide over 100-fold computational efficiency gain over AA models, enabling the study of water crystallization, but have not yet shown to have a LLT. Here, we demonstrate that the CG machine-learned water model Machine-Learned Bond-Order Potential (ML-BOP) has a LLT that ends in a critical point at pc = 170 ± 10 MPa and Tc = 181 ± 3 K. The TLL(p) of ML-BOP is almost identical to the one of TIP4P/2005, adding to the similarity in the equation of state of liquid water in both models. Cooling simulations reveal that ice crystallization is fastest at the LLT and its supercritical continuation of maximum heat capacity, supporting a mechanistic relationship between the structural transformation of water to a low-density liquid (LDL) and ice formation. We find no signature of liquid-liquid criticality in the ice crystallization temperatures. ML-BOP replicates the competition between formation of LDL and ice observed in ultrafast experiments of decompression of the high-density liquid (HDL) into the region of stability of LDL. The simulations reveal that crystallization occurs prior to the coarsening of the HDL and LDL domains, obscuring the distinction between the highly metastable first-order LLT and pronounced structural fluctuations along its supercritical continuation.

Nucleation and growth of crystalline ices from amorphous ices

We here review mostly experimental and some computational work devoted to nucleation in amorphous ices. In fact, there are only a handful of studies in which nucleation and growth in amorphous ices are investigated as two separate processes. In most studies, crystallization temperatures T_x or crystallization rates R_JG are accessed for the combined process. Our Review deals with different amorphous ices, namely, vapor-deposited amorphous solid water (ASW) encountered in many astrophysical environments; hyperquenched glassy water (HGW) produced from μm-droplets of liquid water; and low density amorphous (LDA), high density amorphous (HDA), and very high density amorphous (VHDA) ices produced via pressure-induced amorphization of ice I or from high-pressure polymorphs. We cover the pressure range of up to about 6 GPa and the temperature range of up to 270 K, where only the presence of salts allows for the observation of amorphous ices at such high temperatures. In the case of ASW, its microporosity and very high internal surface to volume ratio are the key factors determining its crystallization kinetics. For HGW, the role of interfaces between individual glassy droplets is crucial but mostly neglected in nucleation or crystallization studies. In the case of LDA, HDA, and VHDA, parallel crystallization kinetics to different ice phases is observed, where the fraction of crystallized ices is controlled by the heating rate. A key aspect here is that in different experiments, amorphous ices of different "purities" are obtained, where "purity" here means the "absence of crystalline nuclei." For this reason, "preseeded amorphous ice" and "nuclei-free amorphous ice" should be distinguished carefully, which has not been done properly in most studies. This makes a direct comparison of results obtained in different laboratories very hard, and even results obtained in the same laboratory are affected by very small changes in the preparation protocol. In terms of mechanism, the results are consistent with amorphous ices turning into an ultraviscous, deeply supercooled liquid prior to nucleation. However, especially in preseeded amorphous ices, crystallization from the preexisting nuclei takes place simultaneously. To separate the time scales of crystallization from the time scale of structure relaxation cleanly, the goal needs to be to produce amorphous ices free from crystalline ice nuclei. Such ices have only been produced in very few studies.

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.

Homogeneous Freezing of Water Using Microfluidics

The homogeneous freezing of water is important in the formation of ice in clouds, but there remains a great deal of variability in the representation of the homogeneous freezing of water in the literature. The development of new instrumentation, such as droplet microfluidic platforms, may help to constrain our understanding of the kinetics of homogeneous freezing via the analysis of monodisperse, size-selected water droplets in temporally and spatially controlled environments. Here, we evaluate droplet freezing data obtained using the Lab-on-a-Chip Nucleation by Immersed Particle Instrument (LOC-NIPI), in which droplets are generated and frozen in continuous flow. This high-throughput method was used to analyse over 16,000 water droplets (86 μm diameter) across three experimental runs, generating data with high precision and reproducibility that has largely been unrepresented in the microfluidic literature. Using this data, a new LOC-NIPI parameterisation of the volume nucleation rate coefficient (JV(T)) was determined in the temperature region of -35.1 to -36.9 °C, covering a greater JV(T) compared to most other microfluidic techniques thanks to the number of droplets analysed. Comparison to recent theory suggests inconsistencies in the theoretical representation, further implying that microfluidics could be used to inform on changes to parameterisations. By applying classical nucleation theory (CNT) to our JV(T) data, we have gone a step further than other microfluidic homogeneous freezing examples by calculating the stacking-disordered ice-supercooled water interfacial energy, estimated to be 22.5 ± 0.7 mJ m-2, again finding inconsistencies when compared to theoretical predictions. Further, we briefly review and compile all available microfluidic homogeneous freezing data in the literature, finding that the LOC-NIPI and other microfluidically generated data compare well with commonly used non-microfluidic datasets, but have generally been obtained with greater ease and with higher numbers of monodisperse droplets.

Reversible structural transformations in supercooled liquid water from 135 to 245 K

A fundamental understanding of the unusual properties of water remains elusive because of the limited data at the temperatures and pressures needed to decide among competing theories. We investigated the structural transformations of transiently heated supercooled water films, which evolved for several nanoseconds per pulse during fast laser heating before quenching to 70 kelvin (K). Water’s structure relaxed from its initial configuration to a steady-state configuration before appreciable crystallization. Over the full temperature range investigated, all structural changes were reversible and reproducible by a linear combination of high- and low-temperature structural motifs. The fraction of the liquid with the high-temperature motif decreased rapidly as the temperature decreased from 245 to 190 K, consistent with the predictions of two-state “mixture” models for supercooled water in the supercritical regime.

Using natural honey as an anti-oxidant and thermodynamically efficient cryoprotectant in embryo vitrification

Embryo cryopreservation is a common practice in reproductive biology and infertility treatments. Despite major improvements over years, the cryoprotectant solutions are still a major source of concern, mostly due to their chemical toxicity and suboptimal protection against cryoinjuries. In this work, we introduced natural honey as a non-permeating cryoprotectant to replace traditionally used sucrose in embryo vitrification. The proposed media were compared with conventional ones by evaluating vitrified/warmed mouse embryos based on their re-expansion, hatching rate and transcription pattern of selected genes involved in heat-shock response, apoptosis and oxidative stress. Despite the similar high re-expansion rate, molecular fingerprint of the cryopreservation is remarkably reduced when honey is used instead of sucrose. The biological response of the proposed media was explained from a fundamental point of view using antioxidant analysis, DSC and GC techniques. It was found that the proposed honey-based medium is less thermodynamically prone to ice formation, which along with its antioxidant capacity can control the production of oxygen radicals and minimize the stress-induced transcriptional response. Furthermore, this work tries to correlate the physico-chemical properties of the vitrification solutions with the cellular and molecular aspects of the cryopreservation and proposes the application of natural cryoprotectants in cryobiology.

Temperature Dependence of Homogeneous Nucleation in Ice

Ice nucleation is a process of great relevance in physics, chemistry, technology, and environmental sciences; much theoretical effort has been devoted to its understanding, but it still remains a topic of intense research. We shed light on this phenomenon by performing atomistic based simulations. Using metadynamics and a carefully designed set of collective variables, reversible transitions between water and ice are able to be simulated. We find that water freezes into a stacking disordered structure with the all-atom transferable intermolecular potential with 4 points/ice (TIP4P/ice) model, and the features of the critical nucleus of nucleation at the microscopic level are revealed. We have also estimated the ice nucleation rates along with other nucleation parameters at different undercoolings. Our results are in agreement with recent experimental and other theoretical works, and they confirm that nucleation is preceded by a large increase in tetrahedrally coordinated water molecules.

Homogeneous ice nucleation rates and crystallization kinetics in transiently-heated, supercooled water films from 188 K to 230 K

The crystallization kinetics of transiently heated, nanoscale water films were investigated for 188 K < T_pulse < 230 K, where T_pulse is the maximum temperature obtained during a heat pulse. The water films, which had thicknesses ranging from approximately 15-30 nm, were adsorbed on a Pt(111) single crystal and heated with ∼10 ns laser pulses, which produced heating and cooling rates of ∼10⁹-10¹⁰ K/s in the adsorbed water films. Because the ice growth rates have been measured independently, the ice nucleation rates could be determined by modeling the observed crystallization kinetics. The experiments show that the nucleation rate goes through a maximum at T = 216 K ± 4 K, and the rate at the maximum is 10^29±1 m^-3 s^-1. The maximum nucleation rate reported here for flat, thin water films is consistent with recent measurements of the nucleation rate in nanometer-sized water drops at comparable temperatures. However, the nucleation rate drops rapidly at lower temperatures, which is different from the nearly temperature-independent rates observed for the nanometer-sized drops. At T ∼ 189 K, the nucleation rate for the current experiments is a factor of ∼10^4-5 smaller than the rate at the maximum. The nucleation rate also decreases for T_pulse > 220 K, but the transiently heated water films are not very sensitive to the smaller nucleation rates at higher temperatures. The crystallization kinetics are consistent with a "classical" nucleation and growth mechanism indicating that there is an energetic barrier for deeply supercooled water to convert to ice.