🧬 biology

A Thermodynamic Analysis of Enhanced Metastability in Isochoric Supercooled Liquids

This paper provides a general thermodynamic proof that isochoric conditions enhance the metastability of supercooled liquids (where the solid is less dense than the liquid) by reducing the Helmholtz driving force for solidification compared to isobaric conditions, thereby suppressing nucleation rates and offering a new dimensionless stability criterion for comparing materials.

Original authors: Boris Rubinsky

Published 2026-04-30

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Boris Rubinsky

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

The Big Picture: Why a Tight Box Keeps Water Liquid Longer

Imagine you have a bottle of water. If you put it in a freezer, it usually turns into ice. But sometimes, if you are very careful, you can cool the water below freezing without it turning solid. This is called supercooling. It's like a "metastable" state—the water wants to freeze, but it hasn't started yet.

Scientists have noticed something strange: if you put water in a rigid, unchangeable box (where the volume cannot change) and cool it down, it stays liquid much longer and is much harder to freeze than if you put it in a flexible container (like a plastic bag) where the pressure stays constant.

This paper, written by Boris Rubinsky, explains why this happens using the laws of thermodynamics. It proves that the rigid box creates a "self-defense mechanism" that makes it harder for ice to start growing.

The Analogy: The "Bouncy Castle" vs. The "Concrete Room"

To understand the difference, let's compare two scenarios:

1. The Flexible Container (Isobaric / Constant Pressure)

Imagine the water is in a giant, soft, stretchy balloon.

The Ice Problem: When water turns into ice, it expands (it gets bigger).
The Reaction: As soon as a tiny bit of ice tries to form, the balloon stretches out to make room. The pressure inside stays the same.
The Result: The ice feels no resistance. It keeps growing easily because the container just gives way. The "push" to freeze remains strong.

2. The Rigid Container (Isochoric / Constant Volume)

Now, imagine the water is inside a super-strong, unbreakable steel box.

The Ice Problem: Ice still tries to expand.
The Reaction: The steel box cannot stretch. As soon as a tiny speck of ice tries to form, it immediately pushes against the walls. Because the ice is less dense (fluffier) than the water, it tries to take up more space. The box fights back, creating a massive amount of pressure inside.
The Result: This pressure acts like a shield. It changes the rules of the game, making it harder for the ice to grow.

The Three-Step "Self-Limiting" Feedback Loop

The paper explains that the rigid box triggers a three-step chain reaction that stops ice from forming. Think of it as a thermostat that automatically turns down the heat when it gets too hot.

Step 1: The Squeeze (Pressure Rise)
Because ice takes up more space than water, the moment a tiny crystal tries to form in the rigid box, the liquid gets squeezed. The box forces the pressure to shoot up immediately.

Analogy: It's like trying to inflate a balloon inside a locked safe. The moment you blow a little air in, the safe pushes back with huge force.

Step 2: The Melting Point Shift (The Clapeyron Effect)
Here is the magic part. For water (and a few other materials like silicon or bismuth), high pressure makes ice melt at lower temperatures.

Normally, water freezes at 0°C (32°F).
But under the high pressure created by the rigid box, the "freezing point" drops. Maybe it now needs to be -5°C or -10°C to freeze.
Analogy: Imagine a hill where you usually roll a ball down to start a race. High pressure digs a new, deeper hole at the bottom of the hill. Now, the ball has to roll much further down before it can start moving.

Step 3: The Gap Closes (Less "Supercooling")
The temperature of the box is fixed by the freezer. Let's say the freezer is set to -2°C.

In the flexible bag: The freezing point is still 0°C. The gap between the freezer (-2°C) and the freezing point (0°C) is 2 degrees. This is a strong "drive" to freeze.
In the rigid box: The pressure raised the freezing point down to -5°C. Now, the gap between the freezer (-2°C) and the new freezing point (-5°C) is only 3 degrees of "safety" (or rather, the driving force is reduced because the system is effectively "warmer" relative to its new freezing point).
Correction/Clarification based on the paper: The paper argues that the effective driving force is reduced. Because the pressure lowers the melting point, the system feels "less cold" relative to its new melting point. The "push" to freeze is weaker.

The "Inequality" (The Mathematical Proof)

The author uses math to prove a simple rule:
The force trying to freeze water in a rigid box is always weaker than the force trying to freeze water in a flexible bag.

He calls this a "thermodynamic inequality." Because the force is weaker, the ice crystals have a much harder time starting. In fact, the paper says that because the rate of freezing depends on the square of this force, even a small reduction in the force leads to a massive drop in how fast ice forms.

The "Stability Number" (A Scorecard for Materials)

The paper introduces a special number, called $\Pi_0$ (Pi-zero).

Think of this as a "Stability Score" for any material.
It is calculated using only basic properties of the material: how much it expands when freezing, how hard it is to compress, and how much energy is needed to melt it.
If a material has a high score, it will be very stable in a rigid box (hard to freeze). If it has a low score, the rigid box won't help much.
This score works for water, but also for materials like silicon, gallium, and bismuth.

What This Means for "Pre-Critical" Fluctuations

The paper also talks about what happens before a full ice crystal forms.

Imagine tiny, microscopic ice specks forming and disappearing constantly due to heat.
In a flexible bag, the "pull" to grow is strong, so these specks are more likely to survive and get bigger.
In a rigid box, the "pull" is weaker. The tiny specks are more likely to melt back into water before they can grow into a dangerous ice crystal.
Analogy: It's like a hill. In the flexible bag, the hill is steep, so a rolling ball (ice speck) speeds up easily. In the rigid box, the hill is flat. The ball rolls a bit, gets tired, and rolls back down.

Summary of Claims

Rigid boxes make supercooling more stable than flexible containers for water and similar materials.
The reason is physics, not magic: The rigid box creates pressure, which lowers the freezing point, which reduces the "drive" to freeze.
It works for any material where the solid takes up more space than the liquid (like water, silicon, bismuth).
It is a universal rule: It doesn't matter how big the box is or what shape it is; the physics holds true.
The result is a "Stability Number" that scientists can use to predict how well a material will behave in a rigid container.

The paper does not claim to solve all freezing problems or promise that this will instantly cure diseases. It strictly provides the thermodynamic proof for why the rigid container method works better than the traditional method, explaining the "why" behind the experimental results.

1. Problem Statement

Supercooling (maintaining a liquid below its equilibrium freezing point without solidifying) is a critical strategy for preserving biological tissues and food, as it avoids the damage caused by ice crystal formation. However, supercooled liquids are metastable and prone to spontaneous nucleation.

While experimental evidence has consistently shown that isochoric (constant-volume) confinement significantly enhances supercooling stability compared to isobaric (constant-pressure) conditions, a general thermodynamic explanation for this phenomenon has been lacking. Previous studies relied on empirical observations or specific system geometries. The central question addressed is: Why does isochoric confinement intrinsically suppress nucleation relative to isobaric conditions for substances where the solid phase is less dense than the liquid phase ( $v_s > v_l$ )?

2. Methodology

The author employs a rigorous thermodynamic approach combining Helmholtz free energy analysis with first-order perturbation theory.

Thermodynamic Potentials: The analysis contrasts the Gibbs free energy ( $G$ ), the natural potential for isobaric systems, with the Helmholtz free energy ( $F$ ), the natural potential for isochoric systems.
Perturbation Parameter: Instead of tracking individual nucleus radii (as in Classical Nucleation Theory), the analysis uses the nascent solid mass fraction ( $\xi \ll 1$ ) as the perturbation variable.
Key Assumptions:
- The system is a rigid container (constant volume $V$ ).
- The substance undergoes a phase transition where the solid is less dense than the liquid ( $v_s > v_l$ ), resulting in a negative Clapeyron slope ( $dT_m/dP < 0$ ). Examples include water/ice, silicon, gallium, and bismuth.
- The derivation focuses on the initial stage of solidification where $\xi$ is infinitesimal.
Analytical Framework:
1. Helmholtz Convexity: Utilizes the convexity of the Helmholtz free energy surface with respect to volume ( $\partial^2 F / \partial V^2 > 0$ ).
2. Volume Conservation: Analyzes the pressure rise ( $\Delta P$ ) required to maintain constant volume when a less-dense solid forms.
3. Clapeyron Relation: Applies the negative slope of the melting curve to determine how the pressure rise affects the equilibrium melting temperature ( $T_m$ ).
4. Driving Force Comparison: Compares the magnitude of the Helmholtz driving force ( $|\Delta f_{isoV}|$ ) in the isochoric system against the Gibbs driving force ( $|\Delta g_{isoP}|$ ) in the isobaric system.

3. Key Contributions

The paper provides three major theoretical advancements:

General Thermodynamic Proof: It derives a strict inequality proving that for any substance with $v_s > v_l$ , the thermodynamic driving force for solidification under constant volume is strictly smaller than under constant pressure:
$|\Delta f_{isoV}| < |\Delta g_{isoP}|$
This holds regardless of system geometry, container size, or surface chemistry.
Mechanism Identification (Self-Limiting Feedback): The paper elucidates a three-step self-limiting feedback loop unique to isochoric conditions:
- Step 1 (Pressure Rise): Due to the convexity of $F(V)$ , the formation of a small amount of less-dense solid ( $\xi$ ) in a rigid container forces an increase in system pressure ( $\Delta P > 0$ ) to conserve volume.
- Step 2 (Melting Point Depression): Because the substance has a negative Clapeyron slope ( $dT_m/dP < 0$ ), this pressure rise lowers the equilibrium melting temperature ( $T_m$ ).
- Step 3 (Reduced Driving Force): Since the external temperature ( $T$ ) is fixed, the reduction in $T_m$ reduces the effective supercooling gap ( $\Delta T_{isoV} = T_m(P) - T$ ), thereby reducing the thermodynamic driving force for further solidification.
Dimensionless Stability Number ( $\Pi_0$ ): The author introduces a new dimensionless parameter to quantify the intrinsic stability of a material under isochoric conditions:
$\Pi_0 = \frac{(v_s - v_l)^2}{L_f \cdot \kappa_l \cdot v_l}$
Where $L_f$ is the latent heat of fusion and $\kappa_l$ is the liquid isothermal compressibility. This number allows for the comparison of metastability across different materials using only bulk thermodynamic data.

4. Results

Inequality Derivation: The analysis yields the relationship:
$\Delta T_{isoV} = \Delta T_{isoP} - \Pi_0 \cdot T_m \cdot \xi + O(\xi^2)$
This confirms that the effective supercooling in an isochoric system is always less than in an isobaric system for the same imposed temperature.
Nucleation Rate Suppression: According to Classical Nucleation Theory (CNT), the nucleation rate $J$ scales exponentially with the inverse square of the driving force ( $J \propto \exp(-1/\Delta g^2)$ ). Even a modest linear reduction in the driving force ( $|\Delta f_{isoV}| < |\Delta g_{isoP}|$ ) results in a disproportionately large exponential suppression of the nucleation rate.
Pre-Critical Fluctuations: The analysis extends beyond the critical barrier height. It suggests that isochoric conditions reduce the thermodynamic bias favoring the growth of pre-critical fluctuations. This increases the probability that thermal fluctuations will dissolve rather than grow, further stabilizing the liquid state.
Universality: The results apply to water, silicon, gallium, bismuth, germanium, and antimony—materials where the solid phase is less dense than the liquid.

5. Significance

Theoretical Foundation: This work provides the first rigorous, geometry-independent thermodynamic proof explaining the enhanced stability of isochoric supercooling. It moves the explanation from empirical observation to fundamental thermodynamic principles.
Application in Cryopreservation: The findings validate and optimize the use of isochoric containers for preserving biological tissues (e.g., cardiac microtissues, pig livers) and organs, explaining why they withstand mechanical stress and thermal cycling better than isobaric systems.
Materials Science: The framework is applicable to semiconductor processing (silicon solidification) and precision casting (gallium, bismuth), offering a new criterion for controlling phase transitions in rigid molds.
Design Criterion: The dimensionless number $\Pi_0$ offers a practical tool for engineers and scientists to predict and compare the stability of different materials under isochoric confinement without needing complex simulations or specific geometric data.

In conclusion, the paper establishes that isochoric supercooling is intrinsically more stable because the rigid volume constraint creates a self-limiting thermodynamic feedback loop that reduces the driving force for solidification, a mechanism absent in isobaric systems.