Stacking disorder in ice I

Interplay between the Formation of Colloidal Clathrate and Cubic Diamond Crystals

Controlling the valency of directional interactions of patchy particles is insufficient for the selective formation of target crystalline structures due to the competition between phases of similar free energy. Examples of such are stacking hybrids of interwoven hexagonal and cubic diamonds with (i) its liquid phase, (ii) arrested glasses, or (iii) clathrates, all depending on the relative patch size, despite being within the one-bond-per-patch regime. Herein, using molecular dynamics simulations, we demonstrate that although tetrahedral patchy particles with narrow patches can assemble into clathrates or stacking hybrids in the bulk, this behavior can be suppressed by the application of external surface potential. Depending on its strength, the selective growth of either cubic diamond crystals or empty sII clathrate cages can be achieved. The formation of a given ordered network depends on the structure of the first adlayer, which is commensurate with the emerging network.

The impact of alcohol and ammonium fluoride on pressure-induced amorphization of cubic structure I clathrate hydrates

We have investigated pressure-induced amorphization (PIA) of an alcohol clathrate hydrate (CH) of cubic structure type I (sI) in the presence of NH4F utilizing dilatometry and x-ray powder diffraction. PIA occurs at 0.98 GPa at 77 K, which is at a much lower pressure than for other CHs of the same structure type. The amorphized CH also shows remarkable resistance against crystallization upon decompression. While amorphized sI CHs could not be recovered previously at all, this is possible in the present case. By contrast to other CHs, the recovery of the amorphized CHs to ambient pressure does not even require a high-pressure annealing step, where recovery without any loss of amorphicity is possible at 120 K and below. Furthermore, PIA is accessible upon compression at unusually high temperatures of up to 140 K, where it reaches the highest degree of amorphicity. Molecular dynamics simulations confirm that polar alcoholic guests, as opposed to non-polar guests, induce cage deformation at lower pressure. The substitution of NH4F into the host-lattice stabilizes the collapsed state more than the crystalline state, thereby enhancing the collapse kinetics and lowering the pressure of collapse.

Structures and Dynamics of Complex Guest Molecules in Confinement, Revealed by Solid-State NMR, Molecular Dynamics, and Calorimetry

This review gives an overview of current trends in the investigation of confined molecules such as water, small and higher alcohols, carbonic acids, ethylene glycol, and non-ionic surfactants, such as polyethylene glycol or Triton-X, as guest molecules in neat and functionalized mesoporous silica materials employing solid-state NMR spectroscopy, supported by calorimetry and molecular dynamics simulations. The combination of steric interactions, hydrogen bonds, and hydrophobic and hydrophilic interactions results in a fascinating phase behavior in the confinement. Combining solid-state NMR and relaxometry, DNP hyperpolarization, molecular dynamics simulations, and general physicochemical techniques, it is possible to monitor these confined molecules and gain deep insights into this phase behavior and the underlying molecular arrangements. In many cases, the competition between hydrogen bonding and electrostatic interactions between polar and non-polar moieties of the guests and the host leads to the formation of ordered structures, despite the cramped surroundings inside the pores.

Observing growth and interfacial dynamics of nanocrystalline ice in thin amorphous ice films

Ice crystals at low temperatures exhibit structural polymorphs including hexagonal ice, cubic ice, or a hetero-crystalline mixture of the two phases. Despite the significant implications of structure-dependent roles of ice, mechanisms behind the growths of each polymorph have been difficult to access quantitatively. Using in-situ cryo-electron microscopy and computational ice-dynamics simulations, we directly observe crystalline ice growth in an amorphous ice film of nanoscale thickness, which exhibits three-dimensional ice nucleation and subsequent two-dimensional ice growth. We reveal that nanoscale ice crystals exhibit polymorph-dependent growth kinetics, while hetero-crystalline ice exhibits anisotropic growth, with accelerated growth occurring at the prismatic planes. Fast-growing facets are associated with low-density interfaces that possess higher surface energy, driving tetrahedral ordering of interfacial H2O molecules and accelerating ice growth. These findings, based on nanoscale observations, improve our understanding on early stages of ice formation and mechanistic roles of the ice interface.

Ice-Water Equilibrium in Nanoscale Confinement

We show that 2D ^{2}H NMR spectra enable valuable insights into the nature of an ice-water equilibrium in nanoscale confinement, which extends over a broad temperature range. In particular, 2D ^{2}H NMR line-shape analysis allows us to determine the timescale on which the coexisting ice and water phases exchange molecules. For D_{2}O in a silica nanopore with a diameter of 5.4 nm, we find that the residence time of a water molecule in either phase is characterized by an NMR exchange time of τ_{X}=5.7  ms at 220 K. Thus, the ice-water equilibrium is highly dynamic, which is an important aspect for an understanding of deeply cooled confined and, possibly, bulk waters.

DFT study of water on graphene: Synergistic effect of multilayer p-doping

Recent experiments related to a study concerning the adsorption of water on graphene have demonstrated the p-doping of graphene, although most of the ab initio calculations predict nearly zero doping. To shed more light on this problem, we have carried out van der Waals density functional theory calculations of water on graphene for both individual water molecules and continuous water layers with coverage ranging from one to eight monolayers. Furthermore, we have paid attention to the influence of the water molecule orientation toward graphene on its doping properties. In this article, we present the results of the band structure and the Bader charge analysis, showing the p-doping of graphene can be synergistically enhanced by putting 4-8 layers of an ice-like water structure on graphene having the water molecules oriented with oxygen atoms toward graphene.

Low-Temperature Behaviors, Cold Crystallization, and Glass Transition in Poly(vinylpyrrolidone) Aqueous Solution

Intermediate water is a kind of water around the biocompatible polymer, such as poly(vinylpyrrolidone) (PVP), that exhibits the phenomenon of cold crystallization. We investigate the low-temperature behavior of PVP aqueous solution using small- and wide-angle X-ray scattering and total neutron scattering measurements. The ice formation speed of the intermediate water is extremely reduced by confinement in the PVP moiety during the cooling process. However, around the glass transition temperature, the water-rich phase expands and orders the hydrogen-bond network, behaving as ice nuclei. During the heating process, cubic ice is formed first and then fills the water-rich region. After saturation of the cubic ice formation, the ice transforms from the cubic to the hexagonal ice form.

Ice formation and its elimination in cryopreservation of bovine oocytes

Damage from ice and potential toxicity of ice-inhibiting cryoprotective agents (CPAs) are key issues in assisted reproduction using cryopreserved oocytes and embryos. We use synchrotron-based time-resolved x-ray diffraction and tools from protein cryocrystallography to characterize ice formation within bovine oocytes after cooling at rates between ∼1000 °C/min and ∼600,000°C /min and during warming at rates between 20,000 and 150,000 °C /min. Maximum crystalline ice diffraction intensity, maximum ice volume, and maximum ice grain size are always observed during warming. All decrease with increasing CPA concentration, consistent with the decreasing free water fraction. With the cooling rates, warming rates and CPA concentrations of current practice, oocytes may show no ice after cooling but always develop substantial ice fractions on warming, and modestly reducing CPA concentrations causes substantial ice to form during cooling. With much larger cooling and warming rates achieved using cryocrystallography tools, oocytes soaked as in current practice remain essentially ice free during both cooling and warming, and when soaked in half-strength CPA solution oocytes remain ice free after cooling and develop small grain ice during warming. These results clarify the roles of cooling, warming, and CPA concentration in generating ice in oocytes, establish the character of ice formed, and suggest that substantial further improvements in warming rates are feasible. Ice formation can be eliminated as a factor affecting post-thaw oocyte viability and development, allowing other deleterious effects of the cryopreservation cycle to be studied, and osmotic stress and CPA toxicity reduced.

Significance Statement

Cryopreservation of oocytes and embryos is critical in assisted reproduction of humans and domestic animals and in preservation of endangered species. Success rates are limited by damage from crystalline ice, toxicity of cryoprotective agents (CPAs), and damage from osmotic stress. Time-resolved x-ray diffraction of bovine oocytes shows that ice forms much more readily during warming than during cooling, that maximum ice fractions always occur during warming, and that the tools and large CPA concentrations of current protocols can at best only prevent ice formation during cooling. Using tools from cryocrystallography that give dramatically larger cooling and warming rates, ice formation can be completely eliminated and required CPA concentrations substantially reduced, expanding the scope for species-specific optimization of post-thaw reproductive outcomes.

Surface premelting of ice far below the triple point

Premelting of ice, a quasi-liquid layer (QLL) at the surface below the melting temperature, was first postulated by Michael Faraday 160 y ago. Since then, it has been extensively studied theoretically and experimentally through many techniques. Existing work has been performed predominantly on hexagonal ice, at conditions close to the triple point. Whether the same phenomenon can persist at much lower pressure and temperature, where stacking disordered ice sublimates directly into water vapor, remains unclear. Herein, we report direct observations of surface premelting on ice nanocrystals below the sublimation temperature using transmission electron microscopy (TEM). Similar to what has been reported on hexagonal ice, a QLL is found at the solid-vapor interface. It preferentially decorates certain facets, and its thickness increases as the phase transition temperature is approached. In situ TEM reveals strong diffusion of the QLL, while electron energy loss spectroscopy confirms its amorphous nature. More significantly, the premelting observed in this work is thought to be related to the metastable low-density ultraviscous water, instead of ambient liquid water as in the case of hexagonal ice. This opens a route to understand premelting and grassy liquid state, far away from the normal water triple point.

Canyon Diablo lonsdaleite is a nanocomposite containing c/h stacking disordered diamond and diaphite

In 1967, a diamond polymorph was reported from hard, diamond-like grains of the Canyon Diablo iron meteorite and named lonsdaleite. This mineral was defined and identified by powder X-ray diffraction (XRD) features that were indexed with a hexagonal unit cell. Since 1967, several natural and synthetic diamond-like materials with XRD data matching lonsdaleite have been reported and the name lonsdaleite was used interchangeably with hexagonal diamond. Its hexagonal structure was speculated to lead to physical properties superior to cubic diamond, and as such has stimulated attempts to synthesize lonsdaleite. Despite numerous reports, several recent studies have provided alternative explanations for the XRD, transmission electron microscopy and Raman data used to identify lonsdaleite. Here, we show that lonsdaleite from the Canyon Diablo diamond-like grains are a nanocomposite material dominated by subnanometre-scale cubic/hexagonal stacking disordered diamond and diaphite domains. These nanostructured elements are intimately intergrown, giving rise to structural features erroneously associated with h diamond. Our data suggest that the diffuse scattering in XRD and the hexagonal features in transmission electron microscopy images reported from various natural and laboratory-prepared samples that were previously used for lonsdaleite identification, in fact arise from cubic/hexagonal stacking disordered diamond and diaphite domains. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 2)'.

Microstructure and crystal order during freezing of supercooled water drops

Supercooled water droplets are widely used to study supercooled water1,2, ice nucleation3-5 and droplet freezing6-11. Their freezing in the atmosphere affects the dynamics and climate feedback of clouds12,13 and can accelerate cloud freezing through secondary ice production14-17. Droplet freezing occurs at several timescales and length scales14,18 and is sufficiently stochastic to make it unlikely that two frozen drops are identical. Here we use optical microscopy and X-ray laser diffraction to investigate the freezing of tens of thousands of water microdrops in vacuum after homogeneous ice nucleation around 234-235 K. On the basis of drop images, we developed a seven-stage model of freezing and used it to time the diffraction data. Diffraction from ice crystals showed that long-range crystalline order formed in less than 1 ms after freezing, whereas diffraction from the remaining liquid became similar to that from quasi-liquid layers on premelted ice19,20. The ice had a strained hexagonal crystal structure just after freezing, which is an early metastable state that probably precedes the formation of ice with stacking defects8,9,18. The techniques reported here could help determine the dynamics of freezing in other conditions, such as drop freezing in clouds, or help understand rapid solidification in other materials.

Observing ice structure of micron-sized vapor-deposited ice with an x-ray free-electron laser

The direct observation of the structure of micrometer-sized vapor-deposited ice is performed at Pohang Accelerator Laboratory x-ray free electron laser (PAL-XFEL). The formation of micrometer-sized ice crystals and their structure is important in various fields, including atmospheric science, cryobiology, and astrophysics, but understanding the structure of micrometer-sized ice crystals remains challenging due to the lack of direct observation. Using intense x-ray diffraction from PAL-XFEL, we could observe the structure of micrometer-sized vapor-deposited ice below 150 K with a thickness of 2-50 μm grown in an ultrahigh vacuum chamber. The structure of the ice grown comprises cubic and hexagonal sequences that are randomly arranged to produce a stacking-disordered ice. We observed that ice with a high cubicity of more than 80% was transformed to partially oriented hexagonal ice when the thickness of the ice deposition grew beyond 5 μm. This suggests that precise temperature control and clean deposition conditions allow μm-thick ice films with high cubicity to be grown on hydrophilic Si3N4 membranes. The low influence of impurities could enable in situ diffraction experiments of ice nucleation and growth from interfacial layers to bulk ice.

Pursuing colloidal diamonds

The endeavor to selectively fabricate a cubic diamond is challenging due to the formation of competing phases such as its hexagonal polymorph or others possessing similar free energy. The necessity to achieve this is of paramount importance since the cubic diamond is the only polymorph exhibiting a complete photonic bandgap, making it a promising candidate in view of photonic applications. Herein, we demonstrate that due to the presence of an external field and delicate manipulation of its strength we can attain selectivity in the formation of a cubic diamond in a one-component system comprised of designer tetrahedral patchy particles. The driving force of such a phenomenon is the structure of the first adlayer which is commensurate with the (110) face of the cubic diamond. Moreover, after a successful nucleation event, once the external field is turned off, the structure remains stable, paving an avenue for further post-synthetic treatment.

Tracking cubic ice at molecular resolution

Ice is present everywhere on Earth and has an essential role in several areas, such as cloud physics, climate change and cryopreservation. The role of ice is determined by its formation behaviour and associated structure. However, these are not fully understood1. In particular, there is a long-standing debate about whether water can freeze to form cubic ice-a currently undescribed phase in the phase space of ordinary hexagonal ice2-6. The mainstream view inferred from a collection of laboratory data attributes this divergence to the inability to discern cubic ice from stacking-disordered ice-a mixture of cubic and hexagonal sequences7-11. Using cryogenic transmission electron microscopy combined with low-dose imaging, we show here the preferential nucleation of cubic ice at low-temperature interfaces, resulting in two types of separate crystallization of cubic ice and hexagonal ice from water vapour deposition at 102 K. Moreover, we identify a series of cubic-ice defects, including two types of stacking disorder, revealing the structure evolution dynamics supported by molecular dynamics simulations. The realization of direct, real-space imaging of ice formation and its dynamic behaviour at the molecular level provides an opportunity for ice research at the molecular level using transmission electron microscopy, which may be extended to other hydrogen-bonding crystals.

Ice crystallization under cryogenic cooling in lipid membrane nanoconfined geometry: Time-resolved structural dynamics

Time-resolved structural investigations of crystallization of water in lipid/protein/salt mesophases at cryogenic temperatures are significant for comprehension of ice nanocrystal nucleation kinetics in lipid membranous systems and can lead to a better understanding of how to experimentally retard the ice formation that obstructs the protein crystal structure determination. Here, we present a time-resolved synchrotron microfocus X-ray diffraction (TR-XRD) study based on ∼40,000 frames that revealed the dynamics of water-to-ice crystallization in a lipid/protein/salt mesophase subjected to cryostream cooling at 100 K. The monoolein/hemoglobin/salt/water system was chosen as a model composition related to protein-loaded lipid cubic phases (LCP) broadly used for the crystallization of proteins. Under confinement in the nanoscale geometry, metastable short-living cubic ice (I_c) rapidly crystallized well before the formation of hexagonal ice (I_h). The detected early nanocrystalline states of water-to-ice transformation in multicomponent systems are relevant to a broad spectrum of technologies and understanding of natural phenomena, including crystallization, physics of water nanoconfinement, and rational design of anti-freezing and cryopreservation systems.

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 H₂O can be a high-energy geophysical material.

Phase diagrams-Why they matter and how to predict them

Understanding the thermodynamic stability and metastability of materials can help us to, for example, gauge whether crystalline polymorphs in pharmaceutical formulations are likely to be durable. It can also help us to design experimental routes to novel phases with potentially interesting properties. In this Perspective, we provide an overview of how thermodynamic phase behavior can be quantified both in computer simulations and machine-learning approaches to determine phase diagrams, as well as combinations of the two. We review the basic workflow of free-energy computations for condensed phases, including some practical implementation advice, ranging from the Frenkel-Ladd approach to thermodynamic integration and to direct-coexistence simulations. We illustrate the applications of such methods on a range of systems from materials chemistry to biological phase separation. Finally, we outline some challenges, questions, and practical applications of phase-diagram determination which we believe are likely to be possible to address in the near future using such state-of-the-art free-energy calculations, which may provide fundamental insight into separation processes using multicomponent solvents.

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.

Accurate crystal structure of ice VI from X-ray diffraction with Hirshfeld atom refinement

Water is an essential chemical compound for living organisms, and twenty of its different crystal solid forms (ices) are known. Still, there are many fundamental problems with these structures such as establishing the correct positions and thermal motions of hydrogen atoms. The list of ice structures is not yet complete as DFT calculations have suggested the existence of additional and - to date - unknown phases. In many ice structures, neither neutron diffraction nor DFT calculations nor X-ray diffraction methods can easily solve the problem of hydrogen atom disorder or accurately determine their anisotropic displacement parameters (ADPs). Here, accurate crystal structures of H₂O, D₂O and mixed (50%H₂O/50%D₂O) ice VI obtained by Hirshfeld atom refinement (HAR) of high-pressure single-crystal synchrotron and laboratory X-ray diffraction data are presented. It was possible to obtain O-H/D bond lengths and ADPs for disordered hydrogen atoms which are in good agreement with the corresponding single-crystal neutron diffraction data. These results show that HAR combined with X-ray diffraction can compete with neutron diffraction in detailed studies of polymorphic forms of ice and crystals of other hydrogen-rich compounds. As neutron diffraction is relatively expensive, requires larger crystals which can be difficult to obtain and access to neutron facilities is restricted, cheaper and more accessible X-ray measurements combined with HAR can facilitate the verification of the existing ice polymorphs and the quest for new ones.

Effect of nitrogen molecules on the growth kinetics at the interface between a (111) plane of cubic ice and water

The molecular-scale growth kinetics of ice from water in the presence of air molecules are still poorly understood, despite their importance for understanding ice particle formation in nature. In this study, a molecular dynamics simulation is conducted to elucidate the molecular-scale growth kinetics at the interface between a (111) plane of cubic ice and water in the presence of N2 molecules. Two potential models of N2 molecules with and without atomic charges are examined. For both models, N2 molecules bind stably to the interface for a period of 1 ns or longer, and the stability of the binding is higher for the charged model than for the noncharged model. Free-energy surfaces of an N2 molecule along the interface and along an ideal (111) plane surface of cubic ice suggest that for both models, the position where an N2 molecule binds stably is different at the interface and on the ideal plane surface, and the stability of the binding is much higher for the interface than for the ideal plane surface. For both models, stacking-disordered ice grows at the interface, and the formation probability of a hexagonal ice layer in the stacking-disordered ice is higher for the charged model than for the uncharged model. The formation probability for the hexagonal ice layer in the stacking-disordered ice depends not only on the stability of binding but also on the positions where N2 molecules bind on the underlying ice, and the number of N2 molecules that bind stably to the underlying ice.

Accurate prediction of ice nucleation from room temperature water

Crystal nucleation is one of the most fundamental processes in the physical sciences and almost always occurs heterogeneously with the aid of a nucleating substrate. No example of nucleation is more ubiquitous and impactful than the formation of ice, vital to fields as diverse as geology, biology, aeronautics, and climate science. However, despite considerable effort, we still cannot predict a priori the efficacy of a nucleating agent. Here we utilize deep learning methods to accurately predict nucleation ability from images of room temperature liquid water-generated from molecular dynamics simulations-on a broad range of substrates. The resulting model, named IcePic, can rapidly and accurately infer nucleation ability, eliminating the requirement for either notoriously expensive simulations or direct experimental measurement. In an online poll, IcePic was found to significantly outperform humans in predicting the ice nucleating efficacy of materials. By analyzing the typical errors made by humans, as well as the application of reverse interpretation methods, physical insights into the role the water contact layer plays in ice nucleation have been obtained. Moving forward, we suggest that IcePic can be used as an easy, cheap, and rapid way to discern the nucleation ability of substrates, also with potential for learning other properties related to interfacial water.

High cubicity of D2O ice inside spherical nanopores of MIL-101(Cr) framework: a neutron diffraction study

Although cubic ice (ice I_c) is considered to be an important phase of water that impacts ice cloud formation in the Earth's upper atmosphere, its properties have not been studied to the same extent as those of hexagonal ice (ice I_h). This is because pristine ice I_c is not formed in simple laboratory conditions. Ice I_c formed in ambient conditions has a stacking disordered array of both hexagonal and cubic-structured hydrogen-bonded water molecules. It is therefore an active area of research to find ways of developing stacking disorder-free pure ice I_c. We demonstrate the evolution of almost pure ice I_c structure within the spherical nanopores of a hydrostable Cr-based metal-organic framework MIL-101(Cr) with an average pore size of 1 nm by low-temperature neutron diffraction study on D₂O. It is observed that at temperatures below 230 K a fraction of liquid D₂O transforms into ice and more than 94% of ice crystals evolved inside the pore are cubic in shape. This is a significantly high fraction of ice I_c formed under simple conditions inside the spherical pores of a Cr-based MOF. It is also observed that upon increasing the temperature, i_ce I_c remains stable until its melting point, without being transformed into ice I_h. This observation is in contrast to our previous observation of ice structure in the 2D cylindrical nanopores of MCM-41, where H₂O ice after creeping out from the cylindrical channel was seen to be dominated by hexagonal shape. In the present study, the D₂O molecules were confined into well-defined spherical nanopores, which hindered the growth of crystals above a certain size, thus minimizing the stacking disordered array. Nanoconfinement of water inside uniform spherical pores is therefore a promising method for the evolution of a significantly large fraction of cubic ice by minimizing the stacking disorder. This finding may open up the possibility of forming ice I_c with 100% cubicity under simple laboratory conditions, which will help in exploring the microphysics of ice cloud formation in the upper atmosphere.

The role of structural order in heterogeneous ice nucleation

The freezing of water into ice is a key process that is still not fully understood. It generally requires an impurity of some description to initiate the heterogeneous nucleation of the ice crystals. The molecular structure, as well as the extent of structural order within the impurity in question, both play an essential role in determining its effectiveness. However, disentangling these two contributions is a challenge for both experiments and simulations. In this work, we have systematically investigated the ice-nucleating ability of the very same compound, cholesterol, from the crystalline (and thus ordered) form to disordered self-assembled monolayers. Leveraging a combination of experiments and simulations, we identify a “sweet spot” in terms of the surface coverage of the monolayers, whereby cholesterol maximises its ability to nucleate ice (which remains inferior to that of crystalline cholesterol) by enhancing the structural order of the interfacial water molecules. These findings have practical implications for the rational design of synthetic ice-nucleating agents.

The freezing of water into ice is still not fully understood. Here, we investigate the role of structural disorder within the biologically relevant impurities that facilitate this fundamental phase transition.

Signatures of sluggish dynamics and local structural ordering during ice nucleation

We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems.

Following the Crystallization of Amorphous Ice after Ultrafast Laser Heating

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

Comparison of the phase transitions of high-pressure phases of ammonium fluoride and ice at ambient pressure

The phase diagrams of water and ammonium fluoride (NH₄F) display some interesting parallels. Several crystalline NH₄F phases have isostructural ice counterparts and one of the famous anomalies of water, the fact that the liquid is denser than ice Ih, is also found for NH₄F. Here, we investigate the phase transitions of the pressure-quenched high-pressure phases of NH₄F upon heating at ambient pressure with x-ray diffraction and calorimetry, and we compare the results with the corresponding ices. NH₄F II transforms to NH₄F Isd, which is a stacking-disordered variant of the stable hexagonal NH₄F Ih polymorph. Heating NH₄F III gives a complex mixture of NH₄F II and NH₄F Isd, while some NH₄F III remains initially. Complete conversion to NH₄F Isd is achieved above ∼220 K. The NH₄F II obtained from NH₄F III persists to much higher temperatures compared to the corresponding pressure-quenched NH₄F II. Quantification of the stacking disorder in NH₄F Isd reveals a more sluggish conversion to NH₄F Ih for NH₄F Isd from NH₄F III. In general, the presence of stress and strain in the samples appears to have pronounced effects on the phase transition temperatures. NH₄F shows a complete lack of amorphous forms at low temperatures either upon low-temperature compression of NH₄F Ih or heating NH₄F III at ambient pressure. The amorphous forms of ice are often used to explain the anomalies of water. It will, therefore, be interesting to explore if liquid NH₄F displays more water-like anomalies despite the apparent lack of amorphous forms at low temperatures.

Heterogeneous Ice Growth in Micron-Sized Water Droplets Due to Spontaneous Freezing

Understanding how ice nucleates and grows into larger crystals is of crucial importance for many research fields. The purpose of this study was to shed light on the phase and structure of ice once a nucleus is formed inside a metastable water droplet. Wide-angle X-ray scattering (WAXS) was performed on micron-sized droplets evaporatively cooled to temperatures where homogeneous nucleation occurs. We found that for our weak hits ice grows more cubic compared to the strong hits that are completely hexagonal. Due to efficient heat removal caused by evaporation, we propose that the cubicity of ice at the vicinity of the droplet’s surface is higher than for ice formed within the bulk of the droplet. Moreover, the Bragg peaks were classified based on their geometrical shapes and positions in reciprocal space, which showed that ice grows heterogeneously with a significant population of peaks indicative of truncation rods and crystal defects. Frequent occurrences of the (100) reflection with extended in-planar structure suggested that large planar ice crystals form at the droplet surface, then fracture into smaller domains to accommodate to the curvature of the droplets. Planar faulting due to misaligned domains would explain the increased cubicity close to the droplet surface.

Patterning Configuration of Surface Hydrophilicity by Graphene Nanosheet towards the Inhibition of Ice Nucleation and Growth

Freezing of liquid water occurs in many natural phenomena and affects countless human activities. The freezing process mainly involves ice nucleation and continuous growth, which are determined by the energy and structure fluctuation in supercooled water. Herein, considering the surface hydrophilicity and crystal structure differences between metal and graphene, we proposed a kind of surface configuration design, which was realized by graphene nanosheets being alternately anchored on a metal substrate. Ice nucleation and growth were investigated by molecular dynamics simulations. The surface configuration could induce ice nucleation to occur preferentially on the metal substrate where the surface hydrophilicity was higher than the lateral graphene nanosheet. However, ice nucleation could be delayed to a certain extent under the hindering effect of the interfacial water layer formed by the high surface hydrophilicity of the metal substrate. Furthermore, the graphene nanosheets restricted lateral expansion of the ice nucleus at the clearance, leading to the formation of a curved surface of the ice nucleus as it grew. As a result, ice growth was suppressed effectively due to the Gibbs–Thomson effect, and the growth rate decreased by 71.08% compared to the pure metal surface. Meanwhile, boundary misorientation between ice crystals was an important issue, which also prejudiced the growth of the ice crystal. The present results reveal the microscopic details of ice nucleation and growth inhibition of the special surface configuration and provide guidelines for the rational design of an anti-icing surface.

Transformation process of ice crystallized from a glassy dilute trehalose aqueous solution

Crystal growth of ice Isd occurring after crystallization of a glassy dilute trehalose aqueous solution is slower than that of ice Isd in a dilute glycerol solution and pure ice Isd, and ice Isd in trehalose aqueous solution survives to ∼230 K.

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.

Freezing of few nanometers water droplets

Water-ice transformation of few nm nanodroplets plays a critical role in nature including climate change, microphysics of clouds, survival mechanism of animals in cold environments, and a broad spectrum of technologies. In most of these scenarios, water-ice transformation occurs in a heterogenous mode where nanodroplets are in contact with another medium. Despite computational efforts, experimental probing of this transformation at few nm scales remains unresolved. Here, we report direct probing of water-ice transformation down to 2 nm scale and the length-scale dependence of transformation temperature through two independent metrologies. The transformation temperature shows a sharp length dependence in nanodroplets smaller than 10 nm and for 2 nm droplet, this temperature falls below the homogenous bulk nucleation limit. Contrary to nucleation on curved rigid solid surfaces, ice formation on soft interfaces (omnipresent in nature) can deform the interface leading to suppression of ice nucleation. For soft interfaces, ice nucleation temperature depends on surface modulus. Considering the interfacial deformation, the findings are in good agreement with predictions of classical nucleation theory. This understanding contributes to a greater knowledge of natural phenomena and rational design of anti-icing systems for aviation, wind energy and infrastructures and even cryopreservation systems.

Desorption-induced evolution of cubic and hexagonal ices in an ultrahigh vacuum and cryogenic temperatures

Reflection absorption infrared spectroscopic investigations of multilayer films of acetonitrile (ACN) and water in an ultrahigh vacuum under isothermal conditions showed the emergence of cubic (ice I_c) and hexagonal (ice I_h) ices depending on the composition of the film. The experiments were conducted with a mixed film of 300 monolayers in thickness and the ACN : H₂O monolayer ratios were varied from 1 : 5 to 5 : 1. Mixed films were deposited at 10 K and warmed to 130-135 K, where ACN desorbed subsequently and IR spectral evolution was monitored continuously. While the emergence of ice I_c at 130 K has been reported, the occurrence of ice I_h at this temperature was seen for the first time. Detailed investigations showed that ice I_h can form at 125 K as well. Crystallization kinetics and activation energy (E_a) for the emergence of ice I_h were evaluated using the Avrami equation.

Effect of Graphene on Ice Polymorph

Recently, ice with stacking disorder structure, consisting of random sequences of cubic ice (Ic) and hexagonal ice (Ih) layers, was reported to be more stable than pure Ih/Ic. Due to a much lower free energy barrier of heterogeneous nucleation, in practice, the freezing process of water is controlled by heterogeneous nucleation triggered by an external medium. Therefore, we carry out molecular dynamic simulations to explore how ice polymorphism depends on the lattice structure of the crystalline substrates on which the ice is grown, focusing on the primary source of atmospheric aerosols, carbon materials. It turns out that, during the nucleation stage, the polymorph of ice nuclei is strongly affected by graphene substrates. For ice nucleation on graphene, we find Ih is the dominant polymorph. This can be attributed to structural similarities between graphene and basal face of Ih. Our results also suggest that the substrate only affects the polymorph of ice close to the graphene surface, with the preference for Ih diminishing as the ice layer grows.

Fcc vs. hcp competition in colloidal hard-sphere nucleation: on their relative stability, interfacial free energy and nucleation rate

Hard-sphere crystallization has been widely investigated over the last six decades by means of colloidal suspensions and numerical methods. However, some aspects of its nucleation behaviour are still under debate. Here, we provide a detailed computational characterisation of the polymorphic nucleation competition between the face-centered cubic (fcc) and the hexagonal-close packed (hcp) hard-sphere crystal phases. By means of several state-of-the-art simulation techniques, we evaluate the melting pressure, chemical potential difference, interfacial free energy and nucleation rate of these two polymorphs, as well as of a random stacking mixture of both crystals. Our results highlight that, despite the fact that both polymorphs have very similar stability, the interfacial free energy of the hcp phase could be marginally higher than that of the fcc solid, which in consequence, mildly decreases its propensity to nucleate from the liquid compared to the fcc phase. Moreover, we analyse the abundance of each polymorph in grown crystals from different types of inserted nuclei: fcc, hcp and stacking disordered fcc/hcp seeds, as well as from those spontaneously emerged from brute force simulations. We find that post-critical crystals fundamentally grow maintaining the polymorphic structure of the critical nucleus, at least until moderately large sizes, since the only crystallographic orientation that allows stacking close-packed disorder is the fcc (111) plane, or equivalently the hcp (0001) one. Taken together, our results contribute with one more piece to the intricate puzzle of colloidal hard-sphere crystallization.

Experimental evidence for glass polymorphism in vitrified water droplets

The question of whether a first-order liquid-to-liquid transition is at the origin of water’s anomalous properties has been controversial since the pioneering experiments by Mishima et al. in 1985 and molecular simulations by Poole et al. in 1992. Since then, experiments aimed at shedding light on this question have been performed using amorphous ices made from crystalline ice, fueling criticism about their crystal-like nature. In the present study, we avoid crystalline ice at any time of the experiment yet still observe a first-order glass-to-glass transition in vitrified liquid droplets. This makes the strong case for glass polymorphism and the direct thermodynamic connection to the liquid-to-liquid transition at higher temperatures, dismissing the criticism voiced for three decades.

The nature of amorphous ices has been debated for more than 35 years. In essence, the question is whether they are related to ice polymorphs or to liquids. The fact that amorphous ices are traditionally prepared from crystalline ice via pressure-induced amorphization has made a clear distinction tricky. In this work, we vitrify liquid droplets through cooling at ≥10⁶ K ⋅ s⁻¹ and pressurize the glassy deposit. We observe a first order–like densification upon pressurization and recover a high-density glass. The two glasses resemble low- and high-density amorphous ice in terms of both structure and thermal properties. Vitrified water shows all features that have been reported for amorphous ices made from crystalline ice. The only difference is that the hyperquenched and pressurized deposit shows slightly different crystallization kinetics to ice I upon heating at ambient pressure. This implies a thermodynamically continuous connection of amorphous ices with liquids, not crystals.

Ice in biomolecular cryocrystallography

Diffraction data acquired from cryocooled protein crystals often include diffraction from ice. Analysis of ice diffraction from crystals of three proteins shows that the ice formed within solvent cavities during rapid cooling is comprised of a stacking-disordered mixture of hexagonal and cubic planes, with the cubic plane fraction increasing with increasing cryoprotectant concentration and increasing cooling rate. Building on the work of Thorn and coworkers [Thorn et al. (2017), Acta Cryst. D73, 729-727], a revised metric is defined for detecting ice from deposited protein structure-factor data, and this metric is validated using full-frame diffraction data from the Integrated Resource for Reproducibility in Macromolecular Crystallography. Using this revised metric and improved algorithms, an analysis of structure-factor data from a random sample of 89 827 PDB entries collected at cryogenic temperatures indicates that roughly 16% show evidence of ice contamination, and that this fraction increases with increasing solvent content and maximum solvent-cavity size. By examining the ice diffraction-peak positions at which structure-factor perturbations are observed, it is found that roughly 25% of crystals exhibit ice with primarily hexagonal character, indicating that inadequate cooling rates and/or cryoprotectant concentrations were used, while the remaining 75% show ice with a stacking-disordered or cubic character.

Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation

Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.

Routes to cubic ice through heterogeneous nucleation

The freezing of water into ice is one of the most important processes in the physical sciences. However, it is still not understood at the molecular level. In particular, the crystallization of cubic ice ([Formula: see text])-rather than the traditional hexagonal polytype ([Formula: see text])-has become an increasingly debated topic. Although evidence for [Formula: see text] is thought to date back almost 400 y, it is only in the last year that pure [Formula: see text] has been made in the laboratory, and these processes involved high-pressure ice phases. Since this demonstrates that pure [Formula: see text] can form, the question naturally arises if [Formula: see text] can be made from liquid water. With this in mind, we have performed a high-throughput computational screening study involving molecular dynamics simulations of nucleation on over 1,100 model substrates. From these simulations, we find that 1) many different substrates can promote the formation of pristine [Formula: see text]; 2) [Formula: see text] can be selectively nucleated for even the mildest supercooling; 3) the water contact layer's resemblance to a face of ice is the key factor determining the polytype selectivity and nucleation temperature, independent of which polytype is promoted; and 4) substrate lattice match to ice is not indicative of the polytype obtained. Through this study, we have deepened understanding of the interplay of heterogeneous nucleation and ice I polytypism and suggest routes to [Formula: see text] More broadly, the substrate design methodology presented here combined with the insight gained can be used to understand and control polymorphism and stacking disorder in materials in general.

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.

Brief Overview of Ice Nucleation

The nucleation of ice is vital in cloud physics and impacts on a broad range of matters from the cryopreservation of food, tissues, organs, and stem cells to the prevention of icing on aircraft wings, bridge cables, wind turbines, and other structures. Ice nucleation thus has broad implications in medicine, food engineering, mineralogy, biology, and other fields. Nowadays, the growing threat of global warming has led to intense research activities on the feasibility of artificially modifying clouds to shift the Earth’s radiation balance. For these reasons, nucleation of ice has been extensively studied over many decades and rightfully so. It is thus not quite possible to cover the whole subject of ice nucleation in a single review. Rather, this feature article provides a brief overview of ice nucleation that focuses on several major outstanding fundamental issues. The author’s wish is to aid early researchers in ice nucleation and those who wish to get into the field of ice nucleation from other disciplines by concisely summarizing the outstanding issues in this important field. Two unresolved challenges stood out from the review, namely the lack of a molecular-level picture of ice nucleation at an interface and the limitations of classical nucleation theory.

Identification of synthesisable crystalline phases of water – a prototype for the challenges of computational materials design

We discuss the identification of experimentally realisable crystalline phases of water to outline and contextualise some of the diverse building blocks of a computational materials design process.

X-ray scattering study on the crystalline and semi-crystalline structure of water/PEG mixtures in their eutectic phase diagram

Mixtures of water and PEG exhibit a well known eutectic phase diagram. While the thermodynamic properties like eutectic and liquidus temperatures as well as the eutectic concentration are intensely investigated almost nothing is known about the structural properties of water and PEG in the different regions of the phase diagram. Therefore, we report on a combined DSC, SAXS and WAXS study over the full range of polymer water compositions in order to elucidate the crystalline and semi-crystalline structure. Throughout the whole phase diagram no signatures of a mixed-crystalline phase of PEG and water can be found. Below the eutectic temperature, both components demix microscopically into hexagonal ice and crystalline PEG with its well known crystalline structure. In the region between eutectic and liquidus temperature, the solid component is composed of a single phase of either pure semi-crystalline PEG (PEG rich side of the phase diagram) or pure ice (water rich side). The semi-crystalline structure of PEG, in contrast, is changed by the presence of water. Its long spacing d_ac increases due to the incorporation of water molecules in the amorphous regions, while the formation of crystalline regions seems to be enhanced, resulting in an almost unaffected crystallinity.

Slow Crystal Growth of Cubic Ice with Stacking Faults in a Glassy Dilute Glycerol Aqueous Solution

Control of ice formation is an important issue as catastrophic ice growth influences our life activities and many industrial systems. We prepared a homogeneous glass of a dilute glycerol aqueous solution by a pressure liquid cooling vitrification method and examined the effect of solute on the ice formation of solvent water using a powder X-ray diffraction method. The solvent water immediately after the crystallization is composed of nanosized pure cubic ice (ice Ic). The crystal growth of ice Ic with stacking faults is much slower than that of pure water. The presence of glycerol molecules dispersing homogeneously may hinder crystal growth. The macroscopic segregation occurs rapidly during the transformation from stacking disordered ice to hexagonal ice. The results suggest that ice formation can be controlled by changing the solute type and concentration. This study has implications for thawing technology in cryobiology and frozen food engineering.

Crystallization in Confinement

Many crystallization processes of great importance, including frost heave, biomineralization, the synthesis of nanomaterials, and scale formation, occur in small volumes rather than bulk solution. Here, the influence of confinement on crystallization processes is described, drawing together information from fields as diverse as bioinspired mineralization, templating, pharmaceuticals, colloidal crystallization, and geochemistry. Experiments are principally conducted within confining systems that offer well-defined environments, varying from droplets in microfluidic devices, to cylindrical pores in filtration membranes, to nanoporous glasses and carbon nanotubes. Dramatic effects are observed, including a stabilization of metastable polymorphs, a depression of freezing points, and the formation of crystals with preferred orientations, modified morphologies, and even structures not seen in bulk. Confinement is also shown to influence crystallization processes over length scales ranging from the atomic to hundreds of micrometers, and to originate from a wide range of mechanisms. The development of an enhanced understanding of the influence of confinement on crystal nucleation and growth will not only provide superior insight into crystallization processes in many real-world environments, but will also enable this phenomenon to be used to control crystallization in applications including nanomaterial synthesis, heavy metal remediation, and the prevention of weathering.