Negative thermal expansion in ice I polytypes

✨

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 water as a master of disguise. We all know it as a liquid, but when it freezes, it doesn't just become "ice." It can dress up in different outfits, called polytypes. The most common outfit is Hexagonal Ice (Ice Ih), the stuff you find in your freezer or in snowflakes. But there's a rare, elusive outfit called Cubic Ice (Ice Ic) that scientists have been trying to put on water for decades, only to fail.

This paper is the story of how the authors finally got water to wear the "Cubic" outfit, took its temperature, and discovered a surprising secret about how it behaves.

Here is the breakdown of their adventure:

1. The Great Ice Heist: Catching the Elusive Cubic Ice

For 50 years, scientists tried to make pure Cubic Ice in the lab. Every time they tried, the water would get confused and put on the Hexagonal outfit instead, or a messy mix of both. It was like trying to teach a chameleon to stay blue; it just kept turning green.

The Breakthrough:
The team used a clever trick. They started with a "skeleton" of ice that had been hollowed out by removing gas molecules (like taking the filling out of a sandwich). This hollow structure, called Ice XVII, was the perfect mold. When they warmed it up just right, the water molecules rearranged themselves into the rare Cubic Ice shape. For the first time, they had a big pile of pure, stable Cubic Ice to study.

2. The "Shrinking" Ice: The Negative Thermal Expansion Mystery

Usually, when things get cold, they shrink. When they get warm, they expand. Think of a balloon: heat it up, it gets bigger; cool it down, it gets smaller.

The Twist:
The authors found that both Hexagonal Ice and their new Cubic Ice do the opposite at very low temperatures. As they get colder, they actually expand (get slightly bigger). As they warm up a little, they shrink.

The Analogy:
Imagine a group of people holding hands in a circle.

Normal behavior: If they get excited (warm up), they dance wildly and push the circle outward (expand).
Ice behavior: At very low temperatures, the "dance" is a specific, low-energy wobble. This wobble actually pulls their hands tighter, making the circle shrink. But as they get a tiny bit warmer, the wobble changes, and they accidentally pull the circle inward before finally expanding again at higher temperatures.

This "shrinking when warming up" is called Negative Thermal Expansion, and the paper confirms that Cubic Ice does this just like its common Hexagonal cousin.

3. The Density Race: Who is Lighter?

The team measured the density (how tightly packed the molecules are) of both ice types.

The Result: Cubic Ice is slightly lighter (less dense) than Hexagonal Ice.
The Catch: At the coldest temperatures, Cubic Ice is about 0.06% lighter. But as the temperature rises, this difference disappears. They become almost identical in weight.

4. The Energy Bill: Why Cubic Ice is "Metastable"

This is the most important part of the paper. The scientists used powerful computer simulations to calculate the enthalpy (think of this as the "energy bill" or the cost to keep the ice in that shape).

The Finding: Hexagonal Ice always has a lower energy bill than Cubic Ice.
The Metaphor: Imagine two houses. House A (Hexagonal) is built on solid rock with a low mortgage. House B (Cubic) is built on a slight hill with a slightly higher mortgage.
- House B can stand there for a long time (it's metastable). It won't collapse immediately.
- But if you give it a little push (like warming it up), it will naturally roll down the hill into House A because that's the more stable, cheaper option.

The paper proves that Cubic Ice is never the "winner." It is always a temporary guest. It exists, but it is always trying to turn into the common Hexagonal Ice.

5. Why Should We Care?

You might ask, "Who cares about a weird type of ice?"

Space Travel: This is crucial for understanding the universe. Many moons and asteroids (like Europa or Titan) are covered in ice. Scientists use telescopes (like the James Webb Space Telescope) to look at these distant worlds. By knowing exactly how Cubic Ice behaves, expands, and shrinks, astronomers can better guess what those icy worlds are made of and how they formed.
The Atmosphere: High up in Earth's atmosphere, where it's freezing cold, Cubic Ice might form in clouds. Understanding its properties helps us predict weather and climate patterns.

The Bottom Line

The authors successfully created a pure sample of the rare Cubic Ice, proved that it behaves strangely by shrinking when it warms up (just like normal ice), and confirmed that it is a "temporary" form of ice that is always trying to turn into the common kind. It's a victory for understanding the hidden complexity of something as simple as frozen water.

1. Problem Statement

The study addresses the fundamental thermodynamic and structural properties of Ice Ic (cubic ice), a metastable polytype of water ice. While Ice Ih (hexagonal ice) is the stable form at ambient pressure and ubiquitous in nature, Ice Ic has historically been difficult to synthesize in a pure form. Previous attempts typically resulted in Ice Isd (stacking-disordered ice), a mixture of cubic and hexagonal sequences, rather than pure cubic ice.

The core scientific questions are:

Can the density and thermal expansion behavior of pure Ice Ic be accurately measured across its entire metastable thermodynamic range?
How does the thermal expansion of Ice Ic compare to Ice Ih, particularly regarding the phenomenon of Negative Thermal Expansion (NTE) observed at low temperatures?
What is the precise enthalpy difference ( $\Delta H$ ) between the two polytypes, and does this confirm the metastable nature of Ice Ic relative to Ice Ih?

2. Methodology

The researchers employed a multi-faceted approach combining advanced sample synthesis, high-resolution neutron diffraction, and quantum molecular dynamics simulations.

Sample Synthesis:
- Precursor: The team utilized Ice XVII (a low-density, hydrogen-filled hydrate of D $_2$ O) as a precursor.
- Process: D $_2$ O ice XVII was synthesized by pressurizing D $_2$ O with H $_2$ gas at 4.3 kbar and 256 K, followed by quenching to 77 K.
- Transformation: The sample was annealed under dynamic vacuum. Heating from 120 K to 160 K (at 0.04 K/min) removed the hydrogen, leaving pure D $_2$ O cubic ice (Ice Ic).
- Conversion: To enable direct comparison, the Ice Ic sample was subsequently heated to transform it into Ice Ih within the same experimental cell.
Neutron Diffraction Experiments:
- Instrument: High-Resolution Powder Diffractometer (HRPD) at the ISIS Neutron and Muon Source (UK).
- Conditions: Measurements were taken in the temperature range of 10 K to 240 K.
- Technique: Time-of-flight neutron diffraction was used to measure lattice parameters with high precision ( $\Delta d/d \approx 6 \times 10^{-4}$ ). Rietveld refinement was applied to extract structural parameters and calculate density ( $\rho$ ).
- Key Observation: The transformation from cubic to hexagonal symmetry was monitored in real-time between 180 K and 220 K.
Computational Simulations:
- Method: Path-Integral Molecular Dynamics (PIMD) simulations were performed to account for nuclear quantum effects (crucial for light molecules like water).
- Potential: The MB-pol(2023) many-body potential was used for high accuracy in modeling water interactions.
- Ensemble: Constant-volume (NVT) ensemble simulations were run for both H $_2$ O and D $_2$ O in the range of 60 K to 205 K to calculate molar enthalpy ( $H$ ).

3. Key Contributions

First Direct Density Measurement of Pure Ice Ic: The study successfully measured the density of pure cubic ice (D $_2$ O) across its full metastable range, a feat previously impossible due to the inability to synthesize pure samples.
Quantification of Negative Thermal Expansion (NTE): The authors confirmed that Ice Ic exhibits NTE at low temperatures, similar to Ice Ih, and provided precise data on the temperature dependence of this phenomenon.
Thermodynamic Validation of Metastability: By combining experimental density data with PIMD simulations, the study provided a rigorous quantitative proof that Ice Ic is metastable relative to Ice Ih at all measured temperatures, resolving long-standing questions about their relative stability.
Enthalpy Jump Calculation: The study calculated the enthalpy difference ( $\Delta H$ ) between the two phases, finding values consistent with recent calorimetric data for H $_2$ O.

4. Key Results

Density and Thermal Expansion:
- Density Difference: At the base temperature (10 K), Ice Ic is 0.06% less dense than Ice Ih. This difference diminishes as temperature increases, nearly vanishing near the metastability limit.
- NTE Behavior: Both polytypes exhibit negative thermal expansion at low temperatures. The minimum in the linear molecular volume expansivity occurs at 33 K for Ice Ic and 29 K for Ice Ih. Above ~60 K, both exhibit positive thermal expansion.
- Mechanism: The NTE is attributed to low-energy phonons causing bond shortening, which overcomes the effect of higher-energy phonons at low temperatures.
Enthalpy and Stability:
- Enthalpy Jump ( $\Delta H$ ): At $T = 205.286$ K, the enthalpy difference $\Delta H = H_{Ih} - H_{Ic}$ was calculated as $-24.1 \pm 2.1$ J/mol for D $_2$ O and $-31.1 \pm 2.0$ J/mol for H $_2$ O.
- Validation: The H $_2$ O simulation result aligns well with experimental calorimetric data ( $\Delta H \approx -37.7$ J/mol at 226 K).
- Conclusion: Since $H_{Ic} > H_{Ih}$ across the measured range, Ice Ic is confirmed to be metastable relative to Ice Ih. No temperature was found where Ice Ic possesses lower enthalpy.
Entropy and Gibbs Free Energy:
- While direct entropy measurements were not available, the authors estimated the entropy jump ( $\Delta S$ ) to be approximately $-0.10$ J/K/mol for D $_2$ O.
- The condition $\Delta H < T\Delta S$ holds for $T < 205$ K, reinforcing the metastability of Ice Ic.

5. Significance

Fundamental Physics: This work provides a complete thermodynamic characterization of Ice Ic, clarifying the subtle energy differences between the two primary Ice I polytypes. It confirms that the "cubicity" of ice is a structural feature of a metastable state, not a distinct stable phase.
Astrophysics and Planetary Science: The findings are critical for modeling icy bodies in the Solar System (e.g., Titan, Enceladus, Europa, Ceres). Understanding the density and stability of Ice Ic helps explain phenomena like cryo-volcanism and the behavior of ice reservoirs under low-temperature, high-vacuum conditions.
Remote Sensing: The precise density and structural data enable better interpretation of data from advanced space instrumentation, such as the James Webb Space Telescope (JWST), allowing for the remote sensing of ice properties (density and cubicity) on distant celestial bodies.
Material Science: The study highlights the potential of ice as a functional material with tunable porous behaviors, opening avenues for applications in low-density materials and hydrate-based technologies.