Looking Inside the Widom Region: Non-Equilibrium Stratification in Supercritical CO2

This study demonstrates that supercritical CO2 under non-equilibrium conditions exhibits spontaneous stratification and Brunt-Vaisala oscillations when crossing Widom lines, revealing that the Widom region functions as a dynamic assembly of phase-like behaviors rather than a homogeneous phase.

The Gist

The Big Idea: A Fluid That Isn't Just "One Thing"

Usually, we think of a fluid (like water or air) as a smooth, uniform soup. If you heat it up, it gets less dense; if you cool it down, it gets denser, but it all mixes together nicely.

However, this paper looks at Supercritical Carbon Dioxide (CO2). Think of this state as a "super-fluid" that is squeezed so hard and heated so much that it's neither a gas nor a liquid. It has the density of a liquid but flows like a gas. Scientists usually assume this super-fluid is perfectly smooth and uniform, even when it's not in perfect balance (non-equilibrium).

The Discovery: The researchers found that when you heat this super-fluid from the bottom and cool it from the top, it doesn't stay smooth. Instead, it spontaneously organizes itself into distinct layers, like a multi-layered cake, even though there are no physical walls separating them.

The Experiment: The "Shadow" Trick

To see this invisible layering, the scientists used a technique called Shadowgraphy.

• The Analogy: Imagine holding a flashlight behind a glass of water. If the water is perfectly clear, the light goes straight through. But if there are tiny ripples or density changes in the water, the light bends, creating shadows or patterns on the wall behind it.
• The Setup: They put a thin layer of supercritical CO2 in a special high-pressure cell. They heated the bottom and cooled the top, creating a temperature gradient.
• The Observation: By taking high-speed photos of the shadows cast by the fluid's density fluctuations, they could "see" how the fluid was moving and vibrating.

The Three Scenarios: From Smooth Cake to Layered Cake

The team ran three different experiments, changing the pressure and temperature to see how the fluid behaved.

1. The "Smooth Cake" (Far from the critical point)

• The Setup: They used conditions where the fluid properties change very slowly from top to bottom.
• The Result: The fluid acted like a single, uniform layer. It wiggled and vibrated at one specific rhythm (frequency).
• The Takeaway: When the fluid is "calm" and far from its critical point, it behaves like a simple, homogeneous fluid.

2. The "Two-Layer Cake" (Crossing the Widom Region)

• The Setup: They increased the temperature difference, pushing the fluid into a special zone called the Widom region. In this zone, the fluid's properties (like how much it expands when heated) change very sharply.
• The Result: Suddenly, the fluid stopped acting like one layer. The data showed two distinct rhythms happening at the same time.
• The Analogy: Imagine a choir singing. In the first experiment, everyone sang the same note. In this one, the choir split into two groups: the bottom half sang a low note, and the top half sang a high note. They were singing together, but they were distinct groups.
• The Takeaway: The fluid had spontaneously stratified into two layers with different physical properties, separated by a transition zone.

3. The "Three-Layer Cake" (Near the Critical Point)

• The Setup: They moved even closer to the critical point (the exact spot where liquid and gas become indistinguishable) and applied a temperature gradient.
• The Result: The fluid split into three distinct layers, each vibrating at its own unique frequency.
• The Takeaway: The closer they got to the critical point, the more the fluid broke apart into different "quasi-phases." One layer acted almost like a liquid, another like a gas, and a middle layer acted as a transition between them.

Why Does This Happen? (The "Gravity vs. Heat" Dance)

The paper explains that this layering happens because of a tug-of-war between heat and gravity.

• The Metaphor: Imagine a crowded dance floor.
  • Heat tries to make everyone move randomly and mix together (diffusion).
  • Gravity tries to keep the heavy people (dense fluid) at the bottom and the light people (less dense fluid) at the top.
  • In the Widom region, the fluid is so sensitive that a tiny change in temperature makes a huge change in density.
  • Because the fluid is so sensitive, the "dance" gets complicated. The heat tries to mix the layers, but gravity pulls them apart. The result is that the fluid organizes itself into stable layers where the "dance steps" (vibrations) are different for each layer.

The "Widom Region" Explained Simply

The paper focuses heavily on the Widom region.

• The Analogy: Think of a hill. Usually, a hill has a gentle slope. But the Widom region is like a cliff edge. If you take one step forward (change the temperature slightly), you drop down a huge amount (the fluid properties change drastically).
• The researchers found that when their experiment crossed this "cliff," the fluid couldn't stay uniform. It had to break into layers to handle the sudden changes in its own properties.

What This Means (According to the Paper)

The paper concludes that the common idea of supercritical fluids being a "smooth, continuous phase" is incomplete.

• The Claim: When you apply a temperature gradient (heat from one side, cold from the other), the supercritical fluid is not homogeneous. It naturally develops a structured, layered architecture.
• The Evidence: They proved this by measuring the "vibrations" (oscillations) of the fluid. Just like you can tell if a room has one echo or three different echoes, they could tell the fluid had one, two, or three distinct layers based on the frequencies they detected.

In summary: This paper shows that supercritical CO2, when heated and cooled, doesn't just mix; it organizes itself into a layered cake of different "quasi-phases," driven by the battle between gravity and the fluid's extreme sensitivity to temperature changes.

Technical Summary

Technical Summary: Non-Equilibrium Stratification in Supercritical CO2

Problem Statement
Supercritical fluids are traditionally described as a continuous phase lacking sharp boundaries between liquid-like and gas-like regions, a view that holds under equilibrium conditions. However, real-world fluid systems are rarely in equilibrium, often featuring inhomogeneities in thermodynamic variables. While the concept of "Widom lines" (loci of maxima in response functions like heat capacity and compressibility) describes the supercritical region near the critical point, these lines have rarely been tested under non-equilibrium conditions. Furthermore, standard equilibrium perspectives neglect the role of fluctuations, which become ubiquitous and significant near critical points or under macroscopic gradients. This study addresses the question of whether a fluid out of equilibrium, specifically supercritical carbon dioxide (CO2) subjected to a stabilizing temperature gradient, can exhibit spontaneous stratification and distinct dynamical behaviors that deviate from the homogeneous equilibrium picture.

Methodology
The authors employed a high-sensitivity shadowgraphy setup to investigate non-equilibrium density fluctuations in pure supercritical CO2 layers. The experiments were conducted in a specialized high-pressure optical cell where a vertical, stabilizing temperature gradient was applied using Peltier elements. Three distinct thermodynamic cases were analyzed:

1. Case 1: Far from the Widom region (15 MPa, mean $T \approx 31.1^\circ$C, $\Delta T = 8$ K).
2. Case 2: Crossing the Widom region (15 MPa, mean $T \approx 40^\circ$C, $\Delta T = 44$ K).
3. Case 3: Close to the critical point, crossing the isochore and Widom region (7.7 MPa, mean $T \approx 31.1^\circ$C, $\Delta T = 4$ K).

Density fluctuations were imaged using dynamic shadowgraphy with a superluminescent diode and an s-CMOS camera. The data were processed using the Differential Dynamic Algorithm (DDA) to compute structure functions (SFs) of the fluctuating intensity. These SFs were analyzed to extract the Intermediate Scattering Function (ISF), which characterizes the time evolution of fluctuations. The analysis relied on fluctuating hydrodynamics (FHD) theory, specifically the linearized fluctuating Oberbeck-Boussinesq equations. To interpret the results, the authors fitted the ISFs with models ranging from a single effective layer (for quasi-homogeneous systems) to a sum of multiple damped oscillatory modes (for stratified systems), where each mode corresponds to a distinct sublayer with specific thermophysical properties.

Key Results

• Case 1 (Far from Widom Region): The fluid behaved as a quasi-homogeneous layer. The structure functions exhibited a single dominant decay mode. For wave numbers $q < q_p$ (where $q_p$ is the crossover wave number), the fluctuations showed damped oscillations (Brunt-Väisälä oscillations) resulting from the coupling between thermal and viscous modes. This allowed for the extraction of a single set of thermophysical properties (thermal diffusivity $\alpha_T$, kinematic viscosity $\nu$, and thermal expansion coefficient $\beta_T$) consistent with literature values for a uniform state.
• Case 2 (Crossing Widom Region): As the system crossed the Widom region, the assumption of a single effective base state failed. The structure functions could not be fitted by a single mode but required a superposition of two distinct propagating modes with different frequencies and decay times. This indicated the spontaneous emergence of fluid stratification into two effective layers with different density gradients and thermal expansion coefficients.
• Case 3 (Near Critical Point): In conditions close to the critical point with a pronounced peak in $\beta_T$, the system required a three-layer model to fit the data. Three distinct oscillation frequencies were observed, corresponding to three effective sublayers: a "quasi-liquid" layer at the bottom, two intermediate layers with drastically changing properties, and a transition region. The thermal expansion coefficients derived from the oscillation frequencies ($\beta_{T,1}, \beta_{T,2}, \beta_{T,3}$) matched the theoretical depth-dependent profiles calculated from NIST data.

Key Contributions

1. Experimental Observation of Non-Equilibrium Stratification: The study provides direct experimental evidence that supercritical CO2 under a stabilizing temperature gradient spontaneously stratifies into distinct dynamical layers, separated by transition regions where thermodynamic properties vary sharply.
2. Hydrodynamic Signatures of Widom Lines: The research demonstrates that the Widom region develops structured dynamical heterogeneity under non-equilibrium conditions, evidenced by the emergence of multiple Brunt-Väisälä frequencies.
3. Methodological Advancement: The authors developed a framework to systematically explore a broad range of thermodynamic states within a single experiment. By analyzing the frequency domain of fluctuations, they can separate and quantify the properties of different depths (sublayers) without spatially resolving the fluid, effectively mapping the phase space continuum.
4. Validation of Fluctuating Hydrodynamics: The results confirm that fluctuating hydrodynamics can describe the complex behavior of supercritical fluids in non-equilibrium states, revealing gravity-driven coupling between thermal and viscous modes that leads to hydrodynamic stratification.

Significance and Claims
The paper claims that the classic equilibrium thermodynamic approach is insufficient to describe the dynamical structure of supercritical fluids under non-equilibrium gradients. Instead, the supercritical state in the Widom region should be viewed as a set of stratified hydrodynamic layers exhibiting various dynamics. The authors emphasize that these findings represent "direct evidence of hydrodynamic stratification" and reveal a coupling between thermal and viscous modes with no analogue in equilibrium fluids.

The study modestly notes that while the observed dynamical stratification suggests underlying microscopic structural heterogeneity, the connection between the two remains an open question, as the measurements probe hydrodynamic fluctuations rather than microscopic structural order. The authors suggest that these findings have implications for applications involving supercritical CO2, such as Carbon Capture, Utilisation and Storage (CCUS), energy systems, and planetary sciences where non-equilibrium gradients are ubiquitous. They propose that future experiments under reduced gravity could further elucidate whether this stratification is dominated by temperature gradients or gravity.