Supercooling of liquids, as described by the Enskog-Vlasov kinetic equation

This paper utilizes an Enskog-Vlasov kinetic model to demonstrate that isochoric cooling allows liquids to reach lower supercooled temperatures than isobaric cooling, while predicting that surface tension diverges at the spinodal temperature due to the emergence of an infinite oscillatory region as the liquid approaches instability.

The Gist

Imagine you have a cup of hot coffee. If you leave it alone, it cools down slowly until it reaches room temperature. But what if you could cool it down really fast, or in a very specific way, so that it gets colder than the freezing point without turning into ice? This is called supercooling. It's like a liquid holding its breath, refusing to turn solid even though it's cold enough to do so.

This paper by E. S. Benilov is like a sophisticated weather forecast for liquids, but instead of predicting rain, it predicts exactly when a liquid will finally "snap" and turn into a solid crystal. The author uses a complex mathematical tool called the Enskog–Vlasov (EV) equation to simulate this process.

Here is a breakdown of the paper's main discoveries using simple analogies:

1. The Tool: A "Crowded Dance Floor" Model

To understand how liquids behave, the author combines two ideas:

• The Bumper Cars (Enskog): Imagine molecules as bumper cars in a very crowded room. They constantly bump into each other. The model accounts for how crowded the room is.
• The Invisible Magnet (Vlasov): Now imagine those bumper cars also have a weak, invisible magnet that pulls them together from a distance. This represents the "van der Waals" force that holds liquids together.

By mixing these two ideas, the author created a simulation that tracks how these "magnetic bumper cars" behave when the room gets very cold.

2. The Big Discovery: The "Spinodal" Breaking Point

The paper calculates a specific temperature called the Spinodal Temperature ($T_s$).

• The Analogy: Think of a liquid as a ball sitting in a valley. As you cool it, the valley gets steeper. At a certain point, the valley disappears, and the ball has nowhere to stay but to roll down into a new shape (a solid crystal).
• The Finding: The paper finds that how you cool the liquid matters. If you cool it while keeping the volume fixed (like in a rigid, unchangeable box), you can get it colder than if you cool it while keeping the pressure fixed (like in a flexible balloon). The "rigid box" method allows the liquid to stay liquid at lower temperatures before it snaps into a solid.

3. The Surface Tension Singularity: The "Shaking Edge"

One of the most striking results concerns the surface tension (the "skin" on the surface of the liquid).

• The Analogy: Imagine the surface of the liquid is a trampoline. As the liquid gets closer to its breaking point ($T_s$), the trampoline starts to vibrate violently.
• The Result: The paper shows that as the liquid approaches this breaking point, a strange "wavy" region appears just under the surface. These waves get bigger and bigger.
• The Singularity: At the exact moment the liquid is about to turn solid, these waves stop fading away and stretch out forever. Because the "skin" of the liquid is trying to contain these infinite waves, the surface tension shoots up to infinity. It's like the surface is screaming, "I can't hold this anymore!"

The author argues this isn't just a math trick; it's a real physical phenomenon. If a liquid is about to crystallize, it starts "radiating" these waves, and the surface tension must diverge (go to infinity) to accommodate them.

4. Testing the Theory: Argon and Water

The author tested this model on several fluids, including Argon (a noble gas) and Water.

• Argon: The model predicts that Argon can be supercooled down to about 40 Kelvin (very cold!) before it spontaneously turns to crystal. This matches reasonably well with experiments, though the experiments had some extra gases mixed in that complicated things.
• Water: The model predicts water can be supercooled to about 250 Kelvin (just below freezing). This is close to what scientists see in experiments, though the model isn't perfect for water because water molecules are complex and rotate, while this model treats them as simple spheres.
• The "No-Man's Land": The paper draws a map showing a "No-Man's Land" region. If you try to cool a liquid into this zone, it becomes unstable and instantly crystallizes. You cannot have a stable liquid there.

5. Why This Matters (According to the Paper)

The author emphasizes that this model is different from older theories.

• Old Way: Some theories try to guess the "microscopic" details of how a crystal starts forming, which is hard to measure and often leads to guessing games.
• This Way: The EV model uses big, easy-to-measure facts (like the temperature where water boils or freezes) to calibrate the math. It doesn't need to guess the tiny details; it just uses the known "personality" of the fluid to predict its breaking point.

Summary

In short, this paper uses a mathematical model of "magnetic bumper cars" to show that:

1. Liquids have a hard limit on how cold they can get before turning solid.
2. How you cool them (in a box vs. a balloon) changes that limit.
3. Right before they turn solid, the liquid's surface starts vibrating wildly, causing the surface tension to theoretically go to infinity.
4. This behavior is a fundamental physical rule that likely applies to all liquids, not just the ones the author calculated.

The paper is a theoretical exploration of the "tipping point" where a liquid loses its patience and becomes a solid.

Technical Summary

Problem Statement
The paper addresses the phenomenon of supercooling in liquids, specifically investigating the thermodynamic limits of stability and the behavior of the liquid-vapor interface as a fluid approaches the spinodal temperature ($T_s$). While classical nucleation theory (CNT) has been widely used to study these phenomena, the author notes a lack of consensus regarding its assumptions in the regime of strong supercooling, where key parameters become ill-defined. Furthermore, existing kinetic models often struggle to bridge the gap between dilute gas dynamics and dense fluid behavior, particularly near phase transitions. The central problem is to determine the spinodal temperature at which a liquid becomes unstable to small perturbations and transitions to a solid, and to analyze the associated behavior of surface tension and interfacial structure using a kinetic approach that incorporates both dense-fluid collisions and long-range intermolecular forces.

Methodology
The study employs the Enskog–Vlasov (EV) kinetic equation, a model proposed by Grmela (1971) that extends the classical Boltzmann equation. The EV equation differs from the Boltzmann equation in two primary ways:

1. Collision Integral: It utilizes Enskog's collision integral, which accounts for the finite size of molecules and high-density effects, rather than the integral for dilute fluids.
2. Long-Range Forces: It includes a Vlasov-style term describing the van der Waals force as a collective mean-field force, analogous to electromagnetic forces in plasmas.

The model is applied to spherical molecules (ignoring rotational degrees of freedom) and calibrated using macroscopic parameters: the critical point and the triple point. The calibration involves determining adjustable parameters, including the molecular diameter ($D$) and coefficients ($c_l$) in the high-density expansion, as well as the van der Waals and Korteweg constants ($a$ and $K$), which characterize the strength of the attractive force and surface tension, respectively.

The analysis proceeds in three stages:

• Thermodynamics: Deriving the equation of state, pressure, and chemical potential from the EV energy and entropy integrals.
• Interface Analysis: Solving a nonlinear integral equation for the density profile of a flat liquid-vapor interface to calculate surface tension ($\gamma$) and interfacial structure.
• Stability Analysis: Analyzing "frozen waves" (harmonic disturbances with zero growth rate) to determine the stability boundary (spinodal curve) in the density-temperature plane. This identifies the onset of instability for both infinite-wavelength (fluid-fluid) and finite-wavelength (fluid-solid) perturbations.

The results are illustrated for Argon, Nitrogen, Oxygen, Fluorine, Carbon Monoxide, and Water. For polyatomic fluids, the non-rotational model is treated as a phenomenological approximation, with the caveat that heat capacities remain inaccurate, though qualitative trends are expected to persist.

Key Contributions and Results

1. Spinodal Temperature and Cooling Trajectories: The paper calculates the spinodal temperature ($T_s$) for various fluids. A key finding is that the attainable $T_s$ depends on the cooling trajectory. Isochoric (constant volume) cooling allows a liquid to reach a lower temperature than isobaric (constant pressure) cooling. For isobaric and interfacial cooling, the instability is triggered by finite-wavelength waves (leading to crystallization), whereas isochoric cooling is destabilized by infinite-wavelength waves.
2. Surface Tension Singularity: The study demonstrates that the surface tension of a supercooled liquid-vapor interface diverges as the temperature approaches $T_s$. This singularity is physically interpreted as the emergence of an oscillatory region on the liquid side of the interface. As $T \to T_s$, the liquid approaches instability, and the interface begins to "radiate" waves that are evanescent (decaying) at higher temperatures but become non-decaying (extending to infinity) at $T_s$.
3. Stability Portraits: The paper constructs stability portraits in the density-temperature plane, distinguishing between regions of fluid stability, fluid-crystal instability, and a "no-man's land" where no stable states exist. These portraits reveal that the distance between the triple point and the liquid-crystal separation curve correlates with the compressibility of the liquid (e.g., Argon is highly compressible, while Water is nearly incompressible).
4. Comparison with Experiments: The theoretical spinodal temperatures are compared with experimental data and molecular dynamics simulations.
   • For Argon, the model predicts a spinodal temperature of 40.5 K for interfacial cooling, which is lower than experimental observations of amorphous argon crystallization (20–24 K). The author attributes this discrepancy to the presence of residual gases and non-equilibrium fluxes in the experiments, which alter the effective spinodal curve.
   • For Water, the model predicts a spinodal temperature of 250 K, slightly higher than experimental nucleation limits (235–243 K). The discrepancy is attributed to the limitations of the non-rotational model for water, particularly its inability to fully capture hydrogen bonding and molecular asymmetry.
   • The predicted divergence of surface tension near $T_s$ is consistent with experimental reports of a "second inflection point" in water's surface tension, though the paper notes that direct observation of the singularity is difficult due to scale resolution requirements.

Significance and Claims
The paper claims that the divergence of surface tension at the spinodal point is a robust physical phenomenon with a clear interpretation: the interface becomes unstable as the liquid approaches the crystallization threshold, causing the interfacial profile to develop long-range oscillations. The author argues that because this effect has a clear physical basis, it should occur regardless of the specific model or approximations used, provided the model adequately describes crystallization.

The work positions the Enskog–Vlasov equation as a tool capable of revealing underlying trends in supercooling using only macroscopic calibration parameters, avoiding the microscopic fitting issues often associated with classical nucleation theory. While the author acknowledges that the current non-rotational model yields only qualitative results for polyatomic fluids (like water) and that the model requires generalization for rotating, asymmetric molecules to achieve quantitative accuracy, the study establishes a framework for understanding the thermodynamic limits of supercooled liquids and the nature of the liquid-solid transition. The paper concludes that the EV model, despite its computational complexity involving high-dimensional integrals, provides a consistent kinetic description of gas, liquid, and solid phases and their transitions.