Global thermodynamics for heat-conducting fluids under weak gravity

This paper develops a global thermodynamic framework for liquid-gas coexistence under weak gravity and heat conduction by constructing a variational free-energy function that decomposes into an effective-gravity configurational term and a residual excess-latent-heat contribution, the latter being essential for recovering fundamental thermodynamic relations and reshaping the free-energy landscape.

The Gist

Imagine you have a tall, clear jar filled with a mixture of water (liquid) and steam (gas). In a normal, still world, gravity pulls the heavy water to the bottom and lets the light steam float to the top. This is the natural order.

Now, imagine you start heating the bottom of the jar while cooling the top. You are forcing heat to flow through the mixture. This paper asks a fascinating question: What happens to the order of the water and steam when you combine gravity with this flow of heat?

The authors, Naoko Nakagawa and Shin-ichi Sasa, use a new way of looking at physics called "global thermodynamics" to solve this puzzle. Here is the story of their findings, explained simply.

1. The Two Forces in a Tug-of-War

Think of the system as a tug-of-war between two invisible teams:

• Team Gravity: This team wants the heavy liquid at the bottom and the light gas at the top.
• Team Heat Flow: This team wants to push the liquid toward the cold side and the gas toward the hot side.

Usually, gravity wins. But if the heat flow is strong enough, it can act like a "fake gravity" that pushes in the opposite direction. The paper introduces a concept called Effective Gravity.

• If the real gravity is stronger, the water stays at the bottom.
• If the heat flow is strong enough, the "Effective Gravity" flips. Suddenly, the water wants to float on top of the steam, defying normal gravity.

2. The "Magic Map" (The Free-Energy Landscape)

To figure out which team wins, the authors created a "magic map" called a Free-Energy Landscape.

• Imagine this map is a hilly terrain.
• The height of the land represents how "uncomfortable" or "expensive" a specific arrangement is.
• The system always wants to roll down to the lowest valley (the most comfortable state).

In a normal jar, there is one deep valley where water is at the bottom. But when you add heat flow, the map changes shape.

• The "Effective Gravity" Part: This part of the map acts like a giant slope. If the slope points one way, the water rolls to the bottom. If the heat flow reverses the slope, the water rolls to the top. This determines the big picture: Which phase is on top?
• The "Residual" Part: This is the tricky part. Even if the big slope tells us where the water goes, there is a tiny, bumpy texture on the ground (the "residual" contribution) that the big slope doesn't show. This texture is caused by the friction of the heat flowing. It doesn't change where the water ends up, but it changes the shape of the hills and valleys around it. It creates strange "metastable" layers right at the boundary where the water meets the steam, making the interface slightly "supercooled" or "superheated."

3. The Surprise: You Can't Just Look at the Bottom of the Valley

The paper makes a very important point about how we measure things.

• If you only look at the lowest point on the map (the final state), you might think the system behaves exactly like a normal gravity system, just with a different gravity strength.
• However, if you want to measure the pressure or the temperature of the system, you cannot just look at that lowest point. You have to look at the shape of the valley walls (the "residual" part).
• Analogy: Imagine a ball sitting in a bowl. If you just look at the ball, you know where it is. But if you want to know how hard the bowl is pushing back on the ball (the pressure), you need to know the curvature of the bowl, not just the ball's position. The "residual" part of the paper is that curvature. Without it, your measurements of pressure and temperature would be wrong.

4. The "Inversion" Experiment

The authors calculated exactly what it would take to see this "Effective Gravity" flip in a real experiment.

• They suggest using a tall, narrow cylinder filled with water and steam.
• By carefully controlling the temperature difference between the top and bottom, and the size of the cylinder, you could reach a "tipping point."
• At this tipping point, the water would suddenly stop sitting at the bottom and start floating on top of the steam, even though gravity is still pulling it down.
• They estimate that water near room temperature is the best candidate for this experiment. The required temperature difference is small (about 0.6 degrees Celsius), and the container size would be manageable (a few centimeters tall).

Summary

In simple terms, this paper shows that when you heat a fluid from one side and cool it from the other, the heat flow acts like a second, invisible gravity.

1. The Big Picture: This "heat gravity" can be strong enough to flip the liquid and gas, making the heavy liquid float.
2. The Fine Print: Even though the big picture is determined by this "heat gravity," the tiny details of the interface (where the liquid meets the gas) are shaped by a leftover "residual" effect of the heat flow.
3. The Measurement: To correctly predict the pressure and other properties of this strange, floating liquid, you must account for both the big "heat gravity" and the tiny "residual" bumps.

The paper provides a mathematical "map" to predict exactly when this flip happens and what the system looks like, suggesting that with a simple jar of water, we could actually see liquid defying gravity due to heat flow.

Technical Summary

Technical Summary: Global Thermodynamics for Heat-Conducting Fluids under Weak Gravity

Problem Statement
The paper addresses the thermodynamic description of liquid-gas coexistence in a fluid subjected to both gravity and a steady heat flow (nonequilibrium steady state). While conventional approaches rely on hydrodynamic balance equations and local equations of state, they require nontrivial interfacial input (e.g., thermal resistance, temperature jumps) that is difficult to model microscopically. The authors aim to determine whether competing macroscopic liquid-gas arrangements (e.g., liquid below gas vs. liquid floating above gas) can be ranked and described using a global thermodynamic balance, specifically within the framework of "global thermodynamics" where the system is described by global variables (total particle number, volume, global temperature) rather than local profiles.

Methodology
The authors construct a variational free-energy function, $F_g$, for a system in a rectangular container with fixed global temperature $\tilde{T}$, volume $V$, particle number $N$, gravity $g$, and imposed temperature difference $\Xi$.

1. Variational Formulation: They define a coarse-grained description where the system is divided into a lower region ($\ell$) and an upper region ($u$) by a formal dividing position $x=l$. Each region is assumed to be in a single phase (liquid or gas). The free energy is expressed as a function of the variational coordinates: the number of particles ($N_\ell$) and volume ($V_\ell$) in the lower region.
2. Decomposition: The central methodological step is the decomposition of the variational free energy into two distinct parts:
   • Effective Gravity Component ($F_{g,eff}$): This term has the same configurational form as the equilibrium weak-gravity free energy, but with the physical gravity $g$ replaced by an effective gravity $g_{eff}$.
   • Residual Component ($F_{g,res}$): This term represents a genuine nonequilibrium contribution proportional to an "excess latent heat" ($\hat{q}_{ex}$). It vanishes when heat conduction is absent ($\Xi=0$).
3. Linear Response Regime: The analysis is conducted in the linear response regime where the dimensionless parameter $\epsilon = \max(|\Xi|/\tilde{T}, |mgL|/k_B\tilde{T}) \ll 1$. Profiles are approximated as piecewise linear, and quantities are expanded to leading order in $\epsilon$.
4. Stability and Landscape Analysis: The authors analyze the Hessian matrix of the free energy to determine the local stability of stationary solutions (heterogeneous liquid-gas coexistence vs. homogeneous single-phase states). They also map the free-energy landscape to understand the barrier geometry and saddle points.

Key Contributions and Results

1. Effective Gravity and Configurational Ordering:
   The paper derives an effective gravity $g_{eff} = g + \frac{\bar{v}}{L} \frac{dp_s}{d\tilde{T}} \Xi$, where $\bar{v}$ is the specific volume and $dp_s/d\tilde{T}$ is the slope of the saturation pressure.

   • The ordering of the two separated liquid-gas configurations (Liquid-Gas vs. Gas-Liquid) is determined solely by the sign of $g_{eff}$.
   • If $g_{eff} > 0$, the configuration with liquid at the bottom and gas at the top is favored. If $g_{eff} < 0$, the configuration with liquid floating above gas is favored.
   • The transition between these configurations at $g_{eff}=0$ is a first-order configurational transition in the sense of global thermodynamics, characterized by a cusp in the minimized free energy.

2. Role of the Residual Contribution:
   While $F_{g,res}$ does not determine which separated configuration has the lower free energy (it is $O(\epsilon^2)$ at the stationary points of the heterogeneous states), it is essential for the thermodynamic consistency of the theory.

   • The fundamental thermodynamic relation (the differential of the free energy) cannot be reconstructed from the effective gravity term alone.
   • The derivatives of the residual term are required to recover the correct laboratory observables, such as the spatially averaged pressure $\bar{p}$ and the global chemical potential $\tilde{\mu}_g$.
   • The residual term accounts for local nonequilibrium signatures, such as the discontinuity of the chemical potential at the interface and the shift of the interface temperature from the equilibrium saturation temperature (leading to supercooled gas or superheated liquid layers near the interface).

3. Free-Energy Landscape and Stability:

   • The authors map the free-energy landscape in the space of variational parameters. They show that the effective gravity controls the global hierarchy of minima (the "valleys" corresponding to the stable configurations).
   • The residual contribution reshapes the local barrier geometry, specifically the "saddle ridge" and "valley" structure near homogeneous states.
   • Stability analysis reveals that while the homogeneous state remains locally stable outside the spinodal region (determined by isothermal compressibility), the heat flow can deform the local barrier shape, potentially causing the disappearance of near-homogeneous stationary roots (a saddle-node bifurcation in the landscape geometry) without changing the global stability of the separated phases.

4. Experimental Feasibility:
   Using NIST saturation-property values for water ($H_2O$), carbon dioxide ($CO_2$), sulfur hexafluoride ($SF_6$), and helium ($He$), the authors estimate experimental scales for detecting the effective-gravity inversion ($g_{eff}=0$).

   • They identify a "window" defined by the cell radius, height, and temperature difference that satisfies constraints against capillary effects and Rayleigh-number-driven convection.
   • Water near room temperature is identified as the most accessible candidate, allowing for a temperature difference of order 1 K with centimeter-scale dimensions, whereas other fluids require much smaller temperature differences or present more challenging scales.

Significance and Claims
The paper claims to provide a rigorous global thermodynamic framework that extends equilibrium statistical mechanics to heat-conducting steady states under gravity. Its primary significance lies in:

• Decoupling Ordering from Thermodynamics: Demonstrating that the macroscopic ordering of phases can be described by an effective field ($g_{eff}$) similar to equilibrium, while the thermodynamic relations (derivatives) require a separate residual term to account for nonequilibrium fluxes.
• Resolving the Interface Problem: By treating the interface as a sharp boundary in a global variational principle, the theory encodes nonequilibrium mismatches (like temperature jumps) into the global free energy without needing a detailed microscopic model of interfacial transport.
• Predictive Power: The theory predicts a specific condition ($g_{eff}=0$) for the inversion of liquid-gas stratification, offering a concrete target for experimental verification.
• Landscape Geometry: It reveals that nonequilibrium driving (heat flow) can qualitatively alter the local geometry of the free-energy landscape (barrier shapes and valley existence) even when the global ordering rule remains simple.

The authors conclude that this framework isolates the "heat-flow-induced configurational bias" from complex hydrodynamic and interfacial mechanisms, providing a baseline for understanding phase coexistence in driven systems. They note that while the theory assumes a sharp interface and does not include surface free energy terms, the residual contribution suggests that apparent interfacial costs in nonequilibrium systems may contain both global and microscopic components.