🔬 materials science

Kinetics of Stacking Order Evolution During Heterogeneous Ice Formation

Using in situ cryogenic transmission electron microscopy and molecular dynamics simulations, this study reveals that heterogeneous ice formation involves a surface-constrained, symmetry-breaking recrystallization process where cubic-ice embryos transition into hexagonal-ice dendrites via dynamic stacking-disordered layers to minimize free energy.

Original authors: Xudan Huang, Zifeng Yuan, Chon-Hei Lo, Huacong Sun, Lei Liao, Hongbo Han, Wenxi Li, Wenlong Wang, Zhi Xu, Lei Liu, Xuedong Bai, Limei Xu, Enge Wang, Lifen Wang

Published 2026-03-03

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Xudan Huang, Zifeng Yuan, Chon-Hei Lo, Huacong Sun, Lei Liao, Hongbo Han, Wenxi Li, Wenlong Wang, Zhi Xu, Lei Liu, Xuedong Bai, Limei Xu, Enge Wang, Lifen Wang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are watching a snowflake being born, but instead of a fluffy cloud, you are zooming in so close that you can see the individual water molecules dancing and locking into place. That is essentially what this research paper does.

The scientists used a super-powered microscope (called a cryo-TEM) to watch ice form in real-time, right down to the atomic level. They discovered that ice doesn't just "snap" into its final shape; it goes through a dramatic, chaotic transformation, almost like a caterpillar turning into a butterfly, but with a twist.

Here is the story of how ice grows, explained simply:

1. The Two "Outfits" Ice Can Wear

Water molecules love to stack up in neat layers. Usually, they wear a Hexagonal outfit (Ice Ih). This is the standard, stable ice you find in your freezer or in snowflakes. It's like a honeycomb pattern.

But sometimes, especially when they are just starting to form, they try on a Cubic outfit (Ice Ic). This is a slightly different, less stable arrangement. Think of it like a "practice run" or a "rehearsal" before the real show.

2. The "Rehearsal" Phase (The Cubic Core)

When water vapor hits a cold surface in a vacuum, the very first tiny cluster of ice that forms is usually Cubic.

The Analogy: Imagine a group of dancers starting a routine. They begin in a tight, circular huddle (the cubic core). This shape is easy to form quickly, but it's not the final performance. It's a bit wobbly and unstable.

3. The "Chaos Zone" (Stacking Disorder)

As the ice starts to grow outward from that tiny cubic core, things get messy. The molecules don't immediately know which outfit to wear. They flip-flop between the Cubic and Hexagonal patterns.

The Analogy: This is the "rehearsal chaos." The dancers are trying out different moves, some in the circle, some in the honeycomb. The layers are jumbled up. The scientists call this Stacking-Disordered Ice. It's like a traffic jam where cars are switching lanes constantly. This "jumbled" layer acts as a bridge, helping the ice transition from the unstable core to the stable final shape.

4. The "Grand Finale" (The Hexagonal Branches)

Eventually, the ice settles down. The chaotic middle layer gives way to a beautiful, stable Hexagonal structure. The ice grows out into long, needle-like branches (dendrites), just like the arms of a snowflake.

The Analogy: The dancers finally find their rhythm. They break out of the huddle and form a long, elegant line (the hexagonal branch). The "bridge" of chaos disappears, and the ice locks into its perfect, stable honeycomb pattern.

Why Does This Happen?

The paper explains that this happens because of a tug-of-war between Surface and Symmetry.

The Surface: The cold metal wall the ice grows on forces the first few molecules to sit in a specific way (the Cubic core). It's like a mold that shapes the clay at the very beginning.
The Symmetry: As the ice gets bigger, it wants to be as comfortable and low-energy as possible. The Hexagonal shape is the most comfortable "chair" for the water molecules to sit in. So, the ice naturally evolves from the "uncomfortable" Cubic start to the "comfortable" Hexagonal finish.

The Cool Twist: Tricking the Ice

The most exciting part of the discovery is that the scientists found a way to stop the transformation.

The Analogy: If the ice grows in a straight line away from the direction it started, it never gets the chance to "switch outfits." It stays stuck in the Cubic "rehearsal" mode forever.
The Result: They managed to create pure Cubic ice (Ice Ic) that doesn't turn into normal ice. This is like finding a way to keep a caterpillar frozen in time before it becomes a butterfly. This could be a huge breakthrough for making new types of materials in the future.

The Big Picture

This research solves a mystery that has puzzled scientists for decades: Why does ice sometimes start as one shape and end as another?

It turns out that nature uses a "chaotic bridge" (the stacking disorder) to help ice transition smoothly from a metastable start to a stable finish. By understanding this dance, we can better predict how snowflakes form in the sky and even design new materials that need specific atomic structures.

In short: Ice doesn't just appear; it evolves. It starts with a wobbly, cubic rehearsal, stumbles through a chaotic middle section, and finally finds its perfect, stable hexagonal groove.

1. Problem Statement

The crystallization of water ice, specifically the selection of stacking order in close-packed polymorphic materials, remains a fundamental enigma. While hexagonal ice ( $I_h$ ) is the thermodynamically stable phase at standard conditions, cubic ice ( $I_c$ ) and stacking-disordered ice ( $I_{sd}$ ) frequently appear during nucleation and growth.

The Conflict: Existing theories are contradictory. Some suggest $I_c$ embryos form first (Ostwald's step rule), while others propose $I_{sd}$ is stabilized by structural entropy.
The Gap: Direct experimental evidence linking the microscopic evolution of stacking order (cubic $\to$ disordered $\to$ hexagonal) with macroscopic morphology (hemispherical $\to$ dendritic) has been scarce due to limitations in temporal and spatial resolution. The mechanism driving the transition from metastable cubic cores to stable hexagonal dendrites remains unclear.

2. Methodology

The authors employed a multi-modal approach combining advanced experimental imaging with computational simulations:

In Situ Cryogenic Transmission Electron Microscopy (Cryo-TEM):
- Conditions: Experiments were conducted under far-from-equilibrium conditions at low temperature (~~102 K) and low pressure (~~ $10^{-6}$ Pa).
- Setup: A metal sheet (copper or gold) served as a cold wall for vapor deposition. Residual water vapor in the TEM column condensed onto the substrate.
- Imaging: High-resolution TEM and Electron Energy Loss Spectroscopy (EELS) were used to track the nucleation and growth of individual ice crystallites in real-time, capturing molecular-resolution stacking sequences.
Molecular Dynamics (MD) Simulations:
- Simulations of vapor-deposition growth were performed to calculate surface energies and free energy landscapes.
- A machine learning force field for water was utilized to validate trends in surface energy and structural stability.
- Simulations tracked the evolution of "cubicity" (the fraction of cubic layers) as a function of crystallite size.

3. Key Contributions

Direct Visualization of Stacking Evolution: The study provides the first direct experimental observation of the atomistic pathway where ice nucleates as a cubic ( $I_c$ ) hemisphere, transitions through a stacking-disordered ( $I_{sd}$ ) intermediate, and recrystallizes into hexagonal ( $I_h$ ) prismatic dendrites.
Mechanism of Symmetry Breaking: The authors identify a "surface-constrained, symmetry-breaking crystallization preference." They demonstrate that the transition is driven by the minimization of free energy, where the surface constraints of the initial embryo dictate the initial cubic symmetry, which is later broken to form the thermodynamically stable hexagonal phase.
Controlled Synthesis of Metastable Phases: The paper establishes a method to stabilize pure-phase, single-crystalline $I_c$ by growing ice off-axis (away from the <111> stacking direction), preventing the symmetry-breaking transition to $I_h$ .

4. Key Results

Morphological and Structural Coupling:
- Nucleation: Ice initially nucleates as a hemispherical $I_c$ core (~20 nm diameter).
- Transition: As growth proceeds, the structure evolves into stacking-disordered layers ( $I_{sd}$ ) that act as a dynamic, fluctuating bridge.
- Final State: The growth culminates in needle-like or branched $I_h$ dendrites. The branching angle of ~70.5° corresponds to the crystallographic angle between equivalent close-packed {111} planes of the $I_c$ core, confirming coherent epitaxial growth.
Growth Kinetics:
- Growth transitions from isotropic (hemispherical core) to anisotropic (dendritic branches).
- The "cubicity" (percentage of cubic layers) decreases non-linearly from the base to the tip of the crystal.
Energetic Drivers:
- MD simulations reveal that while bulk energy differences between $I_c$ , $I_{sd}$ , and $I_h$ are negligible, surface energy is the dominant factor.
- Hemispherical $I_c$ and branch-like $I_h$ are the most stable morphologies for their respective sizes. The system evolves to minimize total free energy by adopting the $I_h$ structure for the growing branches.
Off-Axis Stabilization:
- When growth occurs off the <111> axis of the $I_c$ core, the coherent epitaxial link to $I_h$ is broken. This results in the formation of large, defect-free, single-crystalline $I_c$ (up to 110 nm), effectively "kinetically trapping" the metastable phase.

5. Significance

Resolving the Ice Polymorph Debate: The study reconciles conflicting theories by showing that $I_c$ embryos and $I_{sd}$ intermediates are real, transient kinetic states driven by surface constraints, rather than just thermodynamic equilibria.
Mechanism of Snowflake Formation: The findings offer a mechanistic foundation for understanding snowflake morphogenesis, explaining how microscopic stacking fluctuations lead to macroscopic dendritic structures.
Material Design Principles: The discovery that surface constraints and symmetry breaking dictate phase selection provides a general strategy for designing advanced materials. It suggests that metastable, high-symmetry phases (like $I_c$ ) can be stabilized in other polymorphic systems through controlled off-axis deposition and surface engineering.
Non-Equilibrium Physics: The work provides a rare atomistic picture of a non-equilibrium vapor-solid phase transition, highlighting the critical interplay between symmetry, surface energy, and free energy minimization.