Heterogeneous nucleation of ice on carbon surfaces

Short-Ranged United-Atom Model for Efficient Simulations of Glycerol and Its Aqueous Mixtures

Glycerol, a versatile cryoprotectant, exhibits a complex conformational landscape governed by intra- and intermolecular hydrogen bonds. Capturing its structural and thermodynamic properties in liquid and glass states remains challenging due to discrepancies between NMR, neutron scattering experiments, and all-atom (AA) simulations. While AA simulations are widely used, they overestimate the α-conformation and incur significant computational costs. Coarse-grained (CG) models provide an efficient alternative but have yet to accurately describe glycerol's conformational distribution and thermodynamic behavior. Here, we introduce SR-UA glycerol, a short-ranged united-atom model parametrized to reproduce experimental density, enthalpy of vaporization, conformational distributions from NMR, and radial distribution functions from neutron scattering data. Inspired by the monatomic water (mW) model, SR-UA glycerol employs short-range anisotropic interactions to mimic hydrogen bonding, achieving about 100-fold computational speedup over AA models. The model captures the conformational shift from γγ to αα as glycerol transitions from gas to the liquid phase, emphasizing the role of intermolecular hydrogen bonds in stabilizing open conformations. When combined with mW water, SR-UA glycerol successfully reproduces key features of glycerol-water mixtures, including the decrease in the temperature of maximum density and the dynamical crossover, in agreement with AA simulations across a range of temperatures and concentrations. This work establishes a robust and efficient model to investigate glycerol's behavior in aqueous mixtures, opening the possibility of addressing with molecular simulations the competition between vitrification and crystallization at cryopreservation-relevant conditions.

Ion-Specific Effects under Electric Fields on Ice Nucleation at the Mica Surface

Applying external electric fields to mineral surfaces can have a substantial impact on ice nucleation, influencing both climate and atmospheric systems. While earlier studies have demonstrated that electric fields can enhance ice nucleation on nonmineral surfaces, the mechanisms driving heterogeneous ice nucleation (HIN) on mineral surfaces under electric fields with different surface ions remain unclear. In this study, we investigate the ion-specific effects under electric fields on HIN efficiency using mica surfaces containing various cations. Our findings reveal that an upward electric field significantly boosts HIN of water droplets atop Na-mica surfaces, raising the nucleation temperature by approximately 6 °C. In contrast, mica surfaces with other cations or those exposed to a downward electric field show no change in nucleation temperature or HIN efficiency. Molecular dynamics simulations suggest that Na+ ions detach more easily from the mica surface under an electric field, exposing more of the flat mica lattice and thus possibly promoting ice nucleation. This study offers new insights into the ion-specific effects of electric fields on HIN, providing a deeper understanding of the role of cations and electric fields in ice nucleation processes.

Water Nanofilms Facilitate Ice Crystal Growth across Droplets

A novel mechanism that underlies the peculiar cascading freezing of multiple supercooled droplets on surfaces is reported. The initial ice crystal growth in large droplets is communicated via connected water nanofilms to smaller droplets. Using high-speed imaging, we show that the presence of a water nanofilm with thickness higher than a critical value around 4.5 nm facilitates a freezing front propagation from one droplet to its neighbors on a hydrophilic surface. The freezing propagation rate through water nanofilms is approximately 40% of the recalescence front growth rate within a droplet, both of which are well described by kinetic crystal growth theory.

Fully Atomistic Molecular Dynamics Simulation of Ice Nucleation Near an Antifreeze Protein

Heterogeneous ice nucleation is a widespread phenomenon in nature. Despite extensive research on ice nucleation near biological antifreeze proteins, a probe for ice nucleation and growth processes at the atomic level is still lacking. Herein, we present simulation evidence of the heterogeneous ice nucleation process on the ice-binding surface (IBS) of the Tenebrio molitor antifreeze protein (TmAFP). Our all-atomistic molecular dynamics simulations reveal detailed steps toward precritical nucleus formation from one-dimensional (1D) channel water to a 2D ice nanolayer and, finally, a 3D ice nucleus. Compared with homogeneous ice nucleation under the same supercooling conditions, the IBS of TmAFP can markedly reduce the critical size of the ice embryo and lower the nucleation free energy barrier, thereby favoring ice nucleation. Additionally, through artificial mutation of selected functional groups on the IBS, we gain deeper insights into how the specific functional groups of the IBS affect ice nucleation. We highlight that the carbonyl groups in the backbone play a crucial role by providing fixed locations for channel water. This function is essential for ensuring alignment between the 2D ice nanolayer and the ice lattice structure.

The Role of Geometry on the Ease of Solidification Inside and Out of Cylindrical Nanopores

We investigated the role of a nanoporous particle on the formation of macroscopic solid in the framework of equilibrium thermodynamics and from the free-energy perspective. The model particle has cylindrical pores with equidistant circular openings on the particle surface. We focused on two potentially limiting steps: (i) the solid nucleation from liquid inside a single pore and (ii) the bridging of multiple pores on the particle surface. We examined the nucleation near the liquid-vapor meniscus inside a pore by considering different solid-vapor and solid-pore wall contact angles, as well as the liquid-vapor meniscus angles. For bridging, we quantified the effects of the proximity of neighboring pores and the number of participating pores where we considered two or three pores, placed two different distances apart, and three contact angles of the solid with the particle surface. Except in special cases where an analytical solution could be developed, we determined the equilibrium nucleus and bridge shapes numerically using the Surface Evolver code. The geometry of these equilibrium shapes was the key for correctly calculating the energy barriers. Our results indicate that the meniscus angle can be an important factor in reducing the barrier for nucleation if the internal angles of the solid nucleus satisfy a certain criterion. For the solid growth out of the pores, we found that the barriers were significantly lower in the presence of multiple, closely packed pores compared to the growth from a single pore. This paper is deliberately written with no reference to material properties or a specific process to highlight the generality of geometry-controlled barriers. A direct application where our findings can be particularly valuable is the ice formation in clouds, which is the subject of intensive research in atmospheric sciences for its role in influencing precipitation patterns and hence the climate.

Computational Characterization of the Structure, Energy, Strengths, and Fracture Resistances of Symmetric Tilt Grain Boundaries in Ice

Using an interatomic potential that can capture the tetrahedral configuration of water molecules (H2O) in ice without the need to explicitly track the motion of the O and H atoms, coarse-grained (CG) atomistic simulations are performed here to characterize the structures, energy, cohesive strengths, and fracture resistance of the grain boundaries (GBs) in polycrystalline ice resulting from water freezing. Taking the symmetric tilt grain boundaries (STGBs) with a tilting axis of ⟨0001⟩ as an example, several main findings from our simulations are (i) the GB energy, EGB, exhibits a strong dependence on the GB misorientation angle, θ. The classical Read-Shockley model only predicts the EGB - θ relation reasonably well when θ < 20° or θ > 45° but fails when 20° < θ < 45°; (ii) two "valleys" appear in the EGB-θ landscape. One occurs at θ = 22° for Σ14(2̅31̅0) GB, and the other is at θ = 32° for Σ26(3̅41̅0) GB. These two GBs might be the most common in polycrystalline ice; (iii) all the STGB structures under consideration here are found to be a collection of edge dislocations with a Burgers vector of b = 1/3⟨112̅0⟩. The core structure of this edge dislocation is composed of a pentagon and a heptagon atomic ring. The separation and orientation of the structure units (SUs) at the GB exhibit a strong dependence on θ; (iv) the length of an atomic bond within the SUs, rather than EGB and θ which are often used in the literature, is identified as one controlling parameter that dictates the intrinsic GB cohesive strength; (v) characterization of the fracture resistance of the GB containing an initial crack is beyond the reach of nanoscale atomistic simulations but is feasible in concurrent atomistic-continuum (CAC) simulations that can simultaneously retain the atomic GB structure together with the long-range stress field within one model. The above findings provide researchers with a stepping stone to understand the complex microstructure of polycrystalline ice and its response to external forces from the bottom up. Such knowledge may be consolidated into constitutive rules and then transferred into the higher length scale models, such as cohesive zone finite element models (CZFEMs), for predicting how polycrystalline ice fractures at laboratory and even geophysical length scales.

Effects of quenched disorder on the kinetics and pathways of phase transition in a soft colloidal system

Although impurities are unavoidable in real-world and experimental systems, most numerical studies on nucleation focus on pure (impurity-free) systems. As a result, the role of impurities in phase transitions remains poorly understood, especially for systems with complex free energy landscapes featuring one or more intermediate metastable phases. In this study, we employed Monte Carlo simulations to investigate the effects of static impurities (quenched disorder) of varying length scales and surface morphologies on the crystal nucleation mechanism and kinetics in the Gaussian core model system-a representative model for soft colloidal systems. We first explored how the nucleation free energy barrier and critical cluster size are influenced by the fraction of randomly pinned (or, static) particles (fp) and the size (np) of the pinned region or cluster. Both the nucleation free energy barrier and critical cluster size increase sharply with increasing fp but decrease as np grows for a given fraction of pinned particles, eventually approaching the homogeneous nucleation limit. On examining the impact of impurity's surface morphology on nucleation kinetics, we observed that the nucleation barrier significantly decreases with increasing the impurity (or, seed) size with crystalline surface morphologies with body-centered cubic showing the greatest facilitation. Interestingly, seeds with random surface roughness had little effect on nucleation kinetics. In addition, the polymorphic identity of particles in the final crystalline phase is influenced by both the seed's surface morphology and system size. This study further provides crucial insights into the intricate relationship between surface-induced local structural fluctuations and the selection of the polymorphic identity in the final crystalline phase, which is essential for understanding and controlling crystallization processes in experiments.

Multitwinned Ice Nanocrystals

Multitwinned nanocrystals are commonly found in substances that preferentially adopt tetrahedral local arrangements, but not yet in water crystals. Ice nanocrystals are pivotal in cloud microphysics, and their surfaces become increasingly prominent in determining structure as crystal size decreases. Nevertheless, discussions on nanocrystal structures have predominantly centered on ice polymorphs observed in bulk: hexagonal (Ih), cubic (Ic), and stacking-disordered (Isd) ices. Here, we demonstrate, through molecular dynamics (MD) simulations, that decahedral and icosahedral nanocrystals form from liquid water droplets of a few nanometers in size without violating the ice rule. The brute force spontaneous crystallization is conducted using the mW model, and the thermodynamic stability is examined using the TIP4P/Ice model. During the crystallization process, the formation of twin boundaries precedes the emergence of centers exhibiting 5-fold and icosahedral symmetry. The free energy calculation suggests the icosahedron has comparable stability with ice Ih nanocrystal. The frequent occurrence of these unreported ice nanocrystals aligns with the fact that natural polycrystalline snow crystals predominantly display a 70.5-degree angle between the Ih c-axes of adjacent branches. Moreover, we show that the formation of multitwinned ice nanocrystals is enhanced within a fullerene, providing a potential avenue for experimental observations.

Electric Field-Enhanced Ion Rejection Rate in Freeze Desalination

Freeze desalination is an appealing method for seawater desalination through freezing seawater. The percentage of ions in the liquid phase, which is termed ion rejection rate, is a critical factor affecting the performance of freeze desalination. Improving the ion rejection rate is an important topic for freeze desalination. In this work, we investigate the effects of electric fields on the ion rejection rate during the freezing of seawater through molecular dynamics simulations. It is found that the ion rejection rate increases with increasing electric field strength. The enhanced ion rejection rate is due to the reduction of the energy barrier at the ice-water interface caused by the electric field, which affects the orientation of water molecules and ion-water interactions. However, the electric field hinders the ice growth rate, which affects the productivity of freeze desalination. Nevertheless, the finding in this work offers a new idea to improve the ion rejection rate. Practically, a trade-off needs to be found to optimize the overall performance of freeze desalination.

Ice Nucleation Mechanisms on Platinum Surfaces in PEM Fuel Cells: Effects of Surface Morphology and Wettability

Understanding the ice nucleation mechanism in the catalyst layers (CLs) of proton exchange membrane (PEM) fuel cells and inhibiting icing by designing the CLs can optimize the cold start strategies, which can enhance the performance of PEM fuel cells. Herein, mitigating the structural matching and templating effects by adjusting the surface morphology and wettability can inhibit icing on the platinum (Pt) catalyst surface effectively. The Pt(211) surface can inhibit icing because the atomic spacing of (211) crystalline surface is much larger than the characteristic distance of ice crystal, thereby mitigating the structural matching effects. A water overlayer on the Pt surface induced by the strong attraction of Pt can act as a template for ice layers and plays an important role in the icing process. Buckling of water overlayer due to the larger atomic spacing of (211) crystalline surface mitigates the templating effect and inhibits icing. Moreover, the water overlayer on the hydrophobic Pt(211) surface with fewer water molecules also mitigates the templating effect, which makes ice nucleation more difficult than homogeneous nucleation. These findings reveal the ice nucleation mechanisms on the Pt catalyst surface from the molecular level and are valuable for catalyst designs to inhibit icing in CL.

Icing and Adhesion Behaviors on Surfaces with Varied Lattice Constants

Investigating droplet wetting and icing behavior is crucial for comprehending the principles of surface icing and the design of anti-icing surfaces. In this study, we present the evidence from molecular dynamics (MD) simulations that reveal a hitherto unreported behavior of droplet wetting and icing adhesion on surfaces with lattice constants from 2.7 to 4.5 Å. Here, we observe that the contact angles (CA) of droplets on a face-centered cubic (FCC) lattice surface consistently correlate positively with the lattice constant. Further examination of droplet behavior on an idealized crystal surface reveals that hydrophilic surfaces (e.g., CA = 85°) inhibit freezing more effectively than hydrophobic surfaces (e.g., CA = 97°). This finding contradicts the conventional explanation that hydrophobic surfaces reduce heterogeneous nucleation, thereby delaying icing. This study introduces a mechanistic explanation for the promotion of water icing by hydrophobic surfaces and offers a novel design concept for the development of anti-ice surfaces in future applications.

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k, has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k. By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

Multifunctional Laser-Induced Graphene-Based Microfluidic Chip for High-Performance Oocyte Cryopreservation with Low Concentration of Cryoprotectants

Oocyte cryopreservation is essential in the field of assisted reproduction, but due to the large size and poor environmental tolerance of oocytes, cell freezing technology needs further improvement. Here, a Y-shaped microfluidic chip based on 3D graphene is ingeniously devised by combining laser-induced graphene (LIG) technology and fiber etching technology. The prepared LIG/PDMS microfluidic chip can effectively suppress ice crystal size and delay ice crystal freezing time by adjusting surface hydrophobicity. In addition, LIG endows the microfluidic chip with an outstanding photothermal effect, which allows to sharply increase its surface temperature from 25 to 71.8 °C with 10 s of low-power 808 nm laser irradiation (0.4 W cm-2). Notably, the LIG/PDMS microfluidic chip not only replaces the traditional cryopreservation carriers, but also effectively reduces the dosage of cryoprotectants (CPAs) needed in mouse oocyte cryopreservation. Even when the concentration of CPAs is cut in half (final concentration of 7.5% ethylene glycol (EG) and 7.5% dimethyl sulfoxide (DMSO)), the survival rate of oocytes is still as high as 92.4%, significantly higher than the control group's 85.8%. Therefore, this work provides a novel design strategy to construct multifunctional microfluidic chips for high-performance oocytes cryopreservation.

On the Role of Starchy Grains in Ice Nucleation Processes

Little is known about the role of starchy food on climate change processes like ice nucleation. Here, we investigate the ice nucleation efficiency (INE) of eight different starchy food materials, namely, corn (CO), potato (PO), barley (BA), brown rice (BR), white rice (WR), oats (OA), wheat (WH), and sweet potato (SP), in immersion freezing mode under mixed-phase cloud conditions. Notably, among all these food materials, PO and BA exhibit the highest ice nucleation efficiency with ice nucleation temperatures as high as -4.3 °C (T50 ∼ -7.0 ± 0.5 °C) and -6.5 °C (T50 ∼ -7.2 ± 0.2 °C), respectively. We also explore the effect of environmentally relevant physicochemical conditions on ice nucleation efficiency, including different pH, temperature, UV/O3/NOx exposure, and various cocontaminants. The change in shape, size, surface properties, hydrophobicity, and crystallinity of materials accounted for the altered INE. The increase in shape, size, and hydrophobicity of the sample generally reduces the INE, whereas an increase in crystallinity enhances the INE of the sample under our experimental conditions. The results suggest that environmentally relevant concentrations slightly alter INE, indicating their role as catalysts in environmental matrices. The outcome of studies on the ice nucleation properties of these food-containing aerosols might help in the physicochemical understanding of other biomolecule-induced ice nucleation, which is still an underdeveloped research area.

Using machine learning with atomistic surface and local water features to predict heterogeneous ice nucleation

Heterogeneous ice nucleation (HIN) has applications in climate science, nanotechnology, and cryopreservation. Ice nucleation on the earth's surface or in the atmosphere usually occurs heterogeneously involving foreign substrates, known as ice nucleating particles (INPs). Experiments identify good INPs but lack sufficient microscopic resolution to answer the basic question: What makes a good INP? We employ molecular dynamics (MD) simulations in combination with machine learning (ML) to address this question. Often, the large amount of computational cost required to cross the nucleation barrier and observe HIN in MD simulations is a practical limitation. We use information obtained from short MD simulations of atomistic surface and water models to predict the likelihood of HIN. We consider 153 atomistic substrates with some surfaces differing in elemental composition and others only in terms of lattice parameters, surface morphology, or surface charges. A range of water features near the surface (local) are extracted from short MD simulations over a time interval (≤300 ns) where ice nucleation has not initiated. Three ML classification models, Random Forest (RF), support vector machine, and Gaussian process classification are considered, and the accuracies achieved by all three approaches lie within their statistical uncertainties. Including local water features is essential for accurate prediction. The accuracy of our best RF classification model obtained including both surface and local water features is 0.89 ± 0.05. A similar accuracy can be achieved including only local water features, suggesting that the important surface properties are largely captured by the local water features. Some important features identified by ML analysis are local icelike structures, water density and polarization profiles perpendicular to the surface, and the two-dimensional lattice match to ice. We expect that this work, with its strong focus on realistic surface models, will serve as a guide to the identification or design of substrates that can promote or discourage ice nucleation.

Heterogeneous Ice Nucleation in Model Crystalline Porous Organic Polymers: Influence of Pore Size on Immersion Freezing

Heterogeneous ice nucleation activity is affected by aerosol particle composition, crystallinity, pore size, and surface area. However, these surface properties are not well understood, regarding how they act to promote ice nucleation and growth to form ice clouds. Therefore, synthesized materials for which surface properties can be tuned were examined in immersion freezing mode in this study. To establish the relationship between particle surface properties and efficiency of ice nucleation, materials, here, covalent organic frameworks (COFs), with different pore diameters and degrees of crystallinity (ordering), were characterized. Results showed that out of all the highly crystalline COFs, the sample with a pore diameter between 2 and 3 nm exhibited the most efficient ice nucleation activity. We posit that the highly crystalline structures with ordered pores have an optimal pore diameter where the ice nucleation activity is maximized and that the not highly crystalline structures with nonordered pores have more sites for ice nucleation. The results were compared and discussed in the context of other synthesized porous particle systems. Such studies give insight into how material features impact ice nucleation activity.

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.

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

RSeeds: Rigid Seeding Method for Studying Heterogeneous Crystal Nucleation

Heterogeneous nucleation is the dominant form of liquid-to-solid transition in nature. Although molecular simulations are most uniquely suited to studying nucleation, the waiting time to observe even a single nucleation event can easily exceed the current computational capabilities. Therefore, there exists an imminent need for methods that enable computationally fast and feasible studies of heterogeneous nucleation. Seeding is a technique that has proven to be successful at dramatically expanding the range of computationally accessible nucleation rates in simulation studies of homogeneous crystal nucleation. In this article, we introduce a new seeding method for heterogeneous nucleation called Rigid Seeding (RSeeds). Crystalline seeds are treated as pseudorigid bodies and simulated on a surface with metastable liquid above its melting temperature. This allows the seeds to adapt to the surface and identify favorable seed-surface configurations, which is necessary for reliable predictions of crystal polymorphs that form and the corresponding heterogeneous nucleation rates. We demonstrate and validate RSeeds for heterogeneous ice nucleation on a flexible self-assembled monolayer surface, a mineral surface based on kaolinite, and two model surfaces. RSeeds predicts the correct ice polymorph, exposed crystal plane, and rotation on the surface. RSeeds is semiquantitative and can be used to estimate the critical nucleus size and nucleation rate when combined with classical nucleation theory. We demonstrate that RSeeds can be used to evaluate nucleation rates spanning many orders of magnitude.

Effects of Hydrophobicity Treatment of Gas Diffusion Layers on Ice Crystallization in Polymer Electrolyte Fuel Cells

The unassisted cold-start capability of polymer electrolyte fuel cells (PEFCs) remains challenging for large-scale automotive applications. Various studies have shown that the freezing of produced water at the cathode catalyst layer (CL) and gas diffusion layer (GDL) interface blocks the oxidant gas and leads to a cold-start failure. However, the impact of GDL properties, including substrate, size, and hydrophobicity, on the freezing behavior of supercooled water is yet to be thoroughly investigated. We use differential scanning calorimetry to perform non-isothermal calorimetric measurements on untreated and waterproofed GDLs (Toray TGP-H-060, Freudenberg H23). By conducting a large number of experiments (>100) for each type of GDL, we obtained the corresponding distribution of onset freezing temperature (T_onset) and found noticeable sample-to-sample variations in both untreated and waterproofed GDLs. Furthermore, ice crystallization is affected by GDL wettability, coating load, coating distribution, and GDL size, whereas the impact of the GDL substrate and saturation level is not apparent. The T_onset distribution allows for predicting the capability of PEFC freeze-start and the freezing probability of residual water at a given subzero temperature. Our work paves the way for GDL modifications toward the improved cold-start capability of PEFC by identifying and avoiding the features that systematically trigger the freezing of supercooled water with high probability.

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.

Poly(vinyl alcohol) Molecular Bottlebrushes Nucleate Ice

Ice binding proteins (IBP) have evolved to limit the growth of ice but also to promote ice formation by ice-nucleating proteins (INPs). IBPs, which modulate these seemingly distinct processes, often have high sequence similarities, and molecular size/assembly is hypothesized to be a crucial determinant. There are only a few synthetic materials that reproduce INP function, and rational design of ice nucleators has not been achieved due to outstanding questions about the mechanisms of ice binding. Poly(vinyl alcohol) (PVA) is a water-soluble synthetic polymer well known to effectively block ice recrystallization, by binding to ice. Here, we report the synthesis of a polymeric ice nucleator, which mimics the dense assembly of IBPs, using confined ice-binding polymers in a high-molar-mass molecular bottlebrush. Poly(vinyl alcohol)-based molecular bottlebrushes with different side-chain densities were synthesized via a combination of ring-opening metathesis polymerization (ROMP) and reversible addition-fragmentation chain-transfer (RAFT) polymerization, using "grafting-to" and "grafting-through" approaches. The facile preparation of the PVA bottlebrushes was performed via selective hydrolysis of the acetate of the poly(vinyl acetate) (PVAc) side chains of the PVAc bottlebrush precursors. Ice-binding polymer side-chain density was shown to be crucial for nucleation activity, with less dense brushes resulting in colder nucleation than denser brushes. This bio-inspired approach provides a synthetic framework for probing heterogeneous ice nucleation and a route toward defined synthetic nucleators for biotechnological applications.

The dependence of soot particle ice nucleation ability on its volatile content

Aviation soot can affect contrail and cirrus cloud formation and impact climate. A product of incomplete combustion, soot particles, are fractal and hydrophobic aggregates comprising carbonaceous spheres with complex physicochemical properties. In the cirrus cloud regime, the surface wettability and pore abundance of soot particles are important determinants for their ice nucleation ability via pore condensation and freezing. In the atmosphere, soot particles can undergo various ageing processes which modify their surface chemistry and porosity, thus acting as ice nucleating particles with varying abilities as a function of ageing. In this study, size-selected soot particles were treated by thermal denuding at 573 K in a pure nitrogen (N₂) or synthetic air (N₂ + O₂) flow and then exposed to varying relative humidity conditions at a fixed temperature in the range from 218 to 243 K, to investigate the role of volatile content in the ice nucleation ability. Both organic-lean and organic-rich propane (C₃H₈) flame soot particles, as well as two types of commercially available carbon black soot particles with high and low surface wettability, were tested. The size and mass distribution of soot aerosol were monitored during the ice nucleation experiments. Bulk soot samples also prepared in pure N₂ or synthetic air environments at 573 K were characterised by thermogravimetric analysis, Fourier transform infrared spectroscopy and dynamic vapour sorption measurements, to reveal the relation between denuding volatile content, associated soot particle property modifications and the ice nucleation ability. Our study shows that thermal denuding induces a change in soot particle porosity playing a dominant role in regulating its ice nucleation via the pore condensation and freezing mechanism. The enrichment in mesopore (2-50 nm) availability may enhance soot ice nucleation. The presence of O₂ in the thermal denuding process may introduce new active sites on soot particles for water interaction and increase soot surface wettability. However, these active sites only facilitate soot ice nucleation when mesopore structures are available. We conclude that a change in volatile content modifies both morphological properties and surface chemistry for soot particles, but porosity change plays the dominant role in regulating soot particle ice nucleation ability.

Janus regulation of ice growth by hyperbranched polyglycerols generating dynamic hydrogen bonding

In this study, a new phenomenon describing the Janus effect on ice growth by hyperbranched polyglycerols, which can align the surrounding water molecules, has been identified. Even with an identical polyglycerol, we not only induced to inhibit ice growth and recrystallization, but also to promote the growth rate of ice that is more than twice that of pure water. By investigating the polymer architecture and population, we found that the stark difference in the generation of quasi-structured H2O molecules at the ice/water interface played a crucial role in the outcome of these opposite effects. Inhibition activity was induced when polymers at nearly fixed loci formed steady hydrogen bonding with the ice surface. However, the formation-and-dissociation dynamics of the interfacial hydrogen bonds, originating from and maintained by migrating polymers, resulted in an enhanced quasi-liquid layer that facilitated ice growth. Such ice growth activity is a unique property unseen in natural antifreeze proteins or their mimetic materials.

Partition Constant of Binary Mixtures for the Equilibrium between a Bulk and a Confined Phase

It is well-known that the thermodynamic, kinetic and structural properties of fluids, and in particular of water and its solutions, can be drastically affected in nanospaces. A possible consequence of nanoscale confinement of a solution is the partial segregation of its components. Thereby, confinement in nanoporous materials (NPM) has been proposed as a means for the separation of mixtures. In fact, separation science can take great advantage of NPM due to the tunability of their properties as a function of nanostructure, morphology, pore size, and surface chemistry. Alcohol-water mixtures are in this context among the most relevant systems. However, a quantitative thermodynamic description allowing for the prediction of the segregation capabilities as a function of the material-solution characteristics is missing. In the present study we attempt to fill this vacancy, by contributing a thermodynamic treatment for the calculation of the partition coefficient in confinement. Combining the multilayer adsorption model for binary mixtures with the Young equation, we conclude that the liquid-vapor surface tension and the contact angle of the pure substances can be used to predict the separation ability of a particular material for a given mixture to a semiquantitative extent. Moreover, we develop a Kelvin-type equation that relates the partition coefficient to the radius of the pore, the contact angle, and the liquid-vapor surface tensions of the constituents. To assess the validity of our thermodynamic formulation, coarse grained molecular dynamics simulations were performed on models of alcohol-water mixtures confined in cylindrical pores. To this end, a coarse-grained amphiphilic molecule was parametrized to be used in conjunction with the mW potential for water. This amphiphilic model reproduces some of the properties of methanol such as enthalpy of vaporization and liquid-vapor surface tension, and the minimum of the excess enthalpy for the aqueous solution. The partition coefficient turns out to be highly dependent on the molar fraction, on the interaction between the components and the confining matrix, and on the radius of the pore. A remarkable agreement between the theory and the simulations is found for pores of radius larger than 15 Å.

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.

Homogeneous ice nucleation rates for mW and TIP4P/ICE models through Lattice Mold calculations

Water freezing is the most common liquid-to-crystal phase transition on Earth, however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240K. We employ two different water models, mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic non-polarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih, Ic and a stacking mixture of ice Ih/Ic; reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner which contributes to benchmarking the freezing behaviour of two popular water models.

Guiding epitaxial crystallization of amorphous solids at the nanoscale: interfaces, stress, and precrystalline order

The crystallization of amorphous solids impacts fields ranging from inorganic crystal growth to biophysics. Promoting or inhibiting nanoscale epitaxial crystallization and selecting its final products underpins applications in cryopreservation, semiconductor devices, oxide electronics, quantum electronics, structural and functional ceramics, and advanced glasses. As precursors for crystallization, amorphous solids are distinguished from liquids and gases by the comparatively long relaxation times for perturbations of the mechanical stress and for variations in composition or bonding. These factors allow experimentally controllable parameters to influence crystallization processes and to drive materials towards specific outcomes. For example, amorphous precursors can be employed to form crystalline phases, such as polymorphs of Al2O3, VO2, and other complex oxides, that are not readily accessible via crystallization from a liquid or through vapor-phase epitaxy. Crystallization of amorphous solids can further be guided to produce a desired polymorph, nanoscale shape, microstructure, and orientation of the resulting crystals. These effects can enable advances in applications in electronics, magnetic devices, optics, and catalysis. Directions for the future development of the chemical physics of crystallization from amorphous solids can be drawn from the impact of structurally complex and non-equilibrium atomic arrangements in liquids and the atomic-scale structure of liquid-solid interfaces.

Icephobic/anti-icing properties of superhydrophobic surfaces

In the winter, icing on solid surfaces is a typical occurrence that may create a slew of hassles and even tragedies. Anti-icing surfaces are one of the effective solutions for this kind of problem. The roughness of a superhydrophobic surface traps air and weakens the contact between the solid surface and liquid water, allowing water droplets to be removed before freezing. At present, the conventional anti-icing methods including mechanical or thermal technology are not only surface structure unfriendly but also have the obsessions of low efficiency, high energy consumption and high manufacturing costs. Hence, developing a way to remove ice by just modifying the surface shape or chemical composition with a low surface energy is extremely desirable. Numerous attempts have been made to investigate the evolution of ice nucleation and icing on superhydrophobic surfaces under the direction of the ice nucleation hypothesis. In this paper, the research progress of ice nucleation in recent years is reviewed from theoretical and application. The icephobic surfaces are described using the wettability and classical nucleation theories. The benefits and drawbacks of anti-icing superhydrophobic surface are summarized, as well as deicing methods. Finally, several applications of ice phobic materials are illustrated, and some problems and challenges in the research field are discussed. We believed that this review will be useful in guiding future water freezing initiatives.

Influences of Oxidation Degree and Size on the Ice Nucleation Efficiency of Graphene Oxide

Figuring out the influences of carbonaceous particle properties on ice nucleation is important to atmospheric science, but it is still a challenge, especially for experimental investigations due to the coupling effect of multiple properties. Here we separately investigate the effects of oxidation degree and size, two typical and debated factors, on ice nucleation efficiency by choosing graphene oxide (GO) as the model. The results show that with the decrease of oxidation degree, ice nucleation efficiency increases through decreasing the ice nucleation free energy barrier (ΔGheter*) on GO surface. Interestingly, although the chosen GO sizes are sufficiently large compared with the sizes of critical ice nuclei, the increase of GO size leads to the increase of ΔGheter* and thus the decrease of ice nucleation efficiency, unlike the general thought that ΔGheter* is not affected by the particle size any more when the size of particle increases to several times that of the critical ice nucleus.

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.

Studying Ice with Environmental Scanning Electron Microscopy

Scanning electron microscopy (SEM) is a powerful imaging technique able to obtain astonishing images of the micro- and the nano-world. Unfortunately, the technique has been limited to vacuum conditions for many years. In the last decades, the ability to introduce water vapor into the SEM chamber and still collect the electrons by the detector, combined with the temperature control of the sample, has enabled the study of ice at nanoscale. Astounding images of hexagonal ice crystals suddenly became real. Since these first images were produced, several studies have been focusing their interest on using SEM to study ice nucleation, morphology, thaw, etc. In this paper, we want to review the different investigations devoted to this goal that have been conducted in recent years in the literature and the kind of information, beyond images, that was obtained. We focus our attention on studies trying to clarify the mechanisms of ice nucleation and those devoted to the study of ice dynamics. We also discuss these findings to elucidate the present and future of SEM applied to this field.

Impact of surface nanostructure and wettability on interfacial ice physics

Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface-ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations.

Unraveling Adhesion Strength between Gas Hydrate and Solid Surfaces

Natural gas hydrate is a promising future energy source, but it also poses a huge threat to oil and gas production due to its ability to deposit within and block pipelines. Understanding the atomistic mechanisms of adhesion between the hydrate and solid surfaces and elucidating its underlying key determining factors can shed light on the fundamentals of novel antihydrate materials design. In this study, large-scale molecular simulations are employed to investigate the hydrate adhesion on solid surfaces, especially with focuses on the atomistic structures of intermediate layer and their influences on the adhesion. The results show that the structure of the intermediate layer formed between hydrate and solid surface is a competitive equilibrium of induced growth from both sides, and is regulated by the content of guest molecules. By comparing the fracture behaviors of the hydrate-solid surface system with different intermediate structures, it is found that both the lattice areal density of water structure and the adsorption of guest molecules on the interface together determine the adhesion strength. Based on the analysis of the adhesion strength distribution, we have also revealed the origins of the drastic difference in adhesion among different water structures such as ice and hydrate. Our simulation indicates that ice-adhesion strength is approximately five times that of lowest hydrate adhesion strength. This finding is surprisingly consistent with the available experimental results.

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.

Effect of interfacial dipole on heterogeneous ice nucleation

In this work, we performed molecular dynamics simulations of ice nucleation on a rigid surface model of cubic zinc blende structure with different surface dipole strength and orientation. Our results show that, although substrates are excellently lattice-matched to cubic ice, ice nucleation merely occurred as the interfacial water molecules (IWs) show identical or similar orientations to that of water molecules in cubic ice. Free energy landscapes revealed that, as substrates have non-suitable dipole strength/orientation, there exist large free energy barriers for rotating dipole IWs to the right orientation to trigger ice formation. This study stresses that, beyond the traditional view of lattice match and the similarity of lattice length between the substrate and new-formed crystal, the similarity between molecular orientations of interfacial component and component in the specific new-formed crystalline face is also critical for promoting ice nucleation.

Ice Nucleation Activity of Perfluorinated Organic Acids

Perfluorinated acids (PFAs) are widely used synthetic chemical compounds, highly resistant to environmental degradation. The widespread PFA contamination in remote regions such as the High Arctic implies currently not understood long-range atmospheric transport pathways. Here, we report that perfluorooctanoic acid (PFOA) initiates heterogeneous ice nucleation at temperatures as high as -16 °C. In contrast, the eight-carbon octanoic acid, perfluorooctanesulfonic acid, and deprotonated PFOA showed poor ice nucleating capabilities. The ice nucleation ability of PFOA correlates with the formation of a PFOA monolayer at the air-water interface, suggesting a mechanism in which the aligned hydroxyl groups of the carboxylic acid moieties provide a lattice matching to ice. The ice nucleation capabilities of fluorinated compounds like PFOA might be relevant for cloud glaciation in the atmosphere and the removal of these persistent pollutants by wet deposition.

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.

Supercooling as a potentially improved storage option for commercial kimchi

The supercooling degree (SD), which refers to the difference between the ice nucleation temperature and freezing point of kimchi, varies depending on the type of kimchi, manufacturer, recipe, and manufacturing season. The aim of this study is to investigate the major influencing factors for the supercooled storage of kimchi and to analyze the possibility of supercooled storage for commercial kimchi. Pearson correlation analysis determined that, in commercial kimchi manufactured between March and July 2018, the SD of kimchi correlated to the number of aerobic bacteria (P < 0.01), however, was not associated with lactic acid bacteria. Moreover, the ice nucleation temperature of saline solution inoculated with aerobic bacteria was reduced from -3.03 ± 0.04 to -6.18 ± 0.11 °C by 10 kGy gamma ray sterilization. Meanwhile, the ice nucleation temperatures of 1.8 kg of commercial red cabbage kimchi and 500 g of white cabbage kimchi manufactured in February 2020 were -3.93 ± 0.06 °C and -3.57 ± 0.06 °C, respectively, and they could be stored at -2.5 °C for 12 weeks without freezing. Additionally, supercooled storage of kimchi at -2.5 °C caused a fermentation delay effect compared to control storage at 1 °C, considering the acidity and amount of lactic acid bacteria. Therefore, if the number of aerobic bacteria is controlled during the manufacturing process of kimchi, supercooled storage at temperatures below -2.5 °C may extend the shelf life of kimchi. PRACTICAL APPLICATION: We have shown that aerobic bacteria are the key influencing factor for ice nucleation of kimchi during supercooled storage. Aside from the initial sterilization process, fermentation of kimchi can also be delayed by lowering the storage temperature below -2.5 °C. Moreover, the method of direct cool refrigeration may have an industrial-level application.

Ice nucleation activity of iron oxides via immersion freezing and an examination of the high ice nucleation activity of FeO

Heterogeneous ice nucleation is a common process in the atmosphere, but relatively little is known about the role of different surface characteristics on the promotion of ice nucleation. We have used a series of iron oxides as a model system to study the role of lattice mismatch and defects induced by milling on ice nucleation activity. The iron oxides include wüstite (FeO), hematite (Fe₂O₃), magnetite (Fe₃O₄), and goethite (FeOOH). The iron oxides were characterized by X-ray diffraction (XRD) and Brunauer-Emmett-Teller (BET) surface area measurements. The immersion freezing experiments were performed using an environmental chamber. Wüstite (FeO) had the highest ice nucleation activity, which we attribute to its low lattice mismatch with hexagonal ice and the exposure of Fe-OH after milling. A comparison study of MnO and wüstite (FeO) with milled and sieved samples for each suggests that physical defects alone result in only a slight increase in ice nucleation activity. Despite differences in the molecular formula and surface groups, hematite (Fe₂O₃), magnetite (Fe₃O₄), and goethite (FeOOH) had similar ice nucleation activities, which may be attributed to their high lattice mismatch to hexagonal ice. This study provides further insight into the characteristics of a good heterogeneous ice nucleus and, more generally, helps to elucidate the interactions between aerosol particles and ice particles in clouds.

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.

Atomic Steps Induce the Aligned Growth of Ice Crystals on Graphite Surfaces

Heterogeneous ice nucleation on atmospheric aerosols strongly affects the earth's climate, and at the microscopic level, surface-irregularity-induced ice crystallization behaviors are common but crucial. Because of the lack of visual evidence and effective experimental methods, the mechanism of atomic-structure-dependent ice formation on aerosol surfaces is poorly understood. Here we chose highly oriented pyrolytic graphite (HOPG) to represent soot (a primary aerosol), and environmental scanning electron microscopy (ESEM) was performed for in situ observations of ice formation. We found that hexagonal ice crystals show an aligned growth pattern via a two-stage pathway with one a axis coinciding with the direction of atomic step edges on the HOPG surface. Additionally, the ice crystals grow at a noticeably higher speed along this direction. This study reveals the role of atomic surface defects in heterogeneous ice nucleation and may pave the way to control icing-related processes in practical applications.