Validity of relativistic hydrodynamics beyond local equilibrium

By constructing formal solutions to the Boltzmann moment equations, this paper demonstrates that relativistic hydrodynamics remains effective far from equilibrium not because systems are near equilibrium, but because non-perturbative modes allow it to smoothly interpolate between free-streaming and collective behavior, effectively explaining its success in modeling quark-gluon plasma.

The Gist

Imagine you are trying to predict how a crowd of people moves through a busy train station.

Usually, we have two ways to look at this:

1. The Micro View: You track every single person, their speed, where they are going, and who they bump into. This is incredibly detailed but impossible to calculate for millions of people. In physics, this is like the Boltzmann equation, which tracks individual particles.
2. The Macro View: You ignore the individuals and just look at the "flow" of the crowd. You treat the crowd like a fluid (like water) with properties like pressure and temperature. This is Hydrodynamics.

The Puzzle
For decades, physicists have been puzzled by a specific situation: Quark-Gluon Plasma (QGP). This is a super-hot, super-dense soup of particles created when heavy atoms smash together.

• The Problem: Hydrodynamics is supposed to only work when things are calm and close to "thermal equilibrium" (like a calm lake). But the QGP is created in a violent, chaotic, far-from-equilibrium state (like a tsunami).
• The Surprise: Despite the chaos, hydrodynamics works amazingly well at predicting how this plasma behaves. It's like using a simple "fluid flow" map to predict the movement of a chaotic riot, and the map turns out to be perfect.

The Paper's Solution
This paper, by Reghukrishnan Gangadharan, asks: Why does the simple fluid map work so well when the system is so messy?

The author uses a mathematical tool called the Relaxation Time Approximation (think of it as a simplified rule for how fast particles calm down after a collision) to solve the complex equations exactly. Here is what they found, using some analogies:

1. The "Gradient Series" is a Broken Ladder

Traditionally, physicists tried to fix the hydrodynamic map by adding "corrections" (gradients) to account for the chaos. Imagine trying to climb a ladder to reach the truth.

• The paper shows that this ladder (the mathematical series) is broken. If you keep climbing higher and higher (adding more corrections), the ladder eventually falls apart and gives nonsense answers. It diverges.
• Why? Because the ladder only tries to reach the "calm equilibrium" state. It forgets about the initial chaos.

2. The "Hidden Ghost" (Non-Perturbative Modes)

The paper reveals that the exact solution to the particle equations isn't just the broken ladder. It has two parts:

• Part A: The divergent ladder (the standard hydrodynamic corrections).
• Part B: A "ghost" term that decays exponentially fast. This term carries the memory of the initial conditions (how the system started).

The Analogy: Imagine you throw a stone into a pond.

• The ripples spreading out are the "hydrodynamic" part (the gradient expansion).
• The splash at the moment of impact is the "non-perturbative" part.
• Standard hydrodynamics tries to describe the ripples but ignores the splash. The paper shows that the splash is essential. It fades away quickly, but while it's there, it changes how the ripples behave.

3. The "Smooth Bridge"

The most important discovery is how these two parts interact.
The paper shows that the "ghost" term (the memory of the initial chaos) doesn't just disappear; it effectively renormalizes (rescales) the rules of the fluid.

• Think of the transport coefficients (like viscosity or friction) as the "rules" of the fluid.
• The paper proves that if you take the standard hydrodynamic rules and adjust the numbers (rescale the coefficients) to account for that initial "splash," the simple fluid model suddenly becomes accurate even in the most chaotic, far-from-equilibrium moments.

The Big Picture

The paper argues that hydrodynamics works in heavy-ion collisions not because the system is "close to equilibrium" (which it isn't), but because the mathematical structure of hydrodynamics is flexible enough to interpolate (bridge the gap) between two extremes:

1. Free Streaming: Particles flying apart without hitting each other (the initial chaos).
2. Collective Flow: Particles moving together like a fluid (the final state).

By including the "memory" of the initial state into the fluid's rules (the transport coefficients), the theory naturally covers the transition from chaos to order.

In Summary
The paper claims that the "magic" of hydrodynamics in particle physics isn't a coincidence. It's because the theory, when viewed correctly, contains a hidden mechanism that absorbs the chaotic initial conditions into its own parameters. It's not that the system is calm; it's that the fluid model is smart enough to be "calm" even when the underlying particles are "wild," provided you tweak the model's settings to remember where it started.

Technical Summary

1. Problem Statement

Relativistic hydrodynamics is the standard effective theory used to describe the macroscopic evolution of the Quark-Gluon Plasma (QGP) created in heavy-ion collisions. However, a fundamental paradox exists:

• Theoretical Limit: Hydrodynamics is derived under the assumption of local thermal equilibrium and small gradients (small Knudsen number, $Kn \ll 1$).
• Experimental Reality: The QGP is created in a highly anisotropic, far-from-equilibrium state with large gradients.
• The Puzzle: Despite these conditions, viscous hydrodynamics provides quantitatively accurate descriptions of experimental observables (e.g., flow coefficients, particle spectra) even in small systems (p-Pb, p-p) and at very early times.
• Existing Explanations: Previous attempts to resolve this include the concept of hydrodynamic attractors (universal solutions converging before equilibrium) and rescaled transport coefficients. However, attractors fail to explain universality in non-conformal, 3+1D systems, and the mechanism behind rescaled coefficients often lacks a rigorous first-principles derivation.

The paper aims to provide a systematic analytic framework explaining why and how hydrodynamics remains effective far from equilibrium by analyzing the exact solutions of the Boltzmann equation.

2. Methodology

The author employs a rigorous analytic approach based on kinetic theory and singular perturbation theory:

1. Relaxation Time Approximation (RTA): The study utilizes the Boltzmann equation with the Anderson-Witting RTA collision kernel. This allows for exact analytic solutions while retaining the essential physics of relaxation toward equilibrium.
2. Moment Hierarchy: Instead of solving for the full distribution function $f(x,p)$, the author derives evolution equations for the moments of the distribution function ($\rho_{n,l}$), which correspond to macroscopic quantities like energy density and pressure anisotropies.
3. Exact Solution Construction:
   • The moment equations are treated as a linear system of differential equations.
   • Using operator methods (Liouville operators and propagators), the author constructs the exact formal solution for the moments.
   • This solution is decomposed into two distinct parts: a gradient expansion (perturbative) and exponentially decaying non-perturbative modes (transient terms).
4. Dimensional Analysis:
   • 0+1D (Bjorken Flow): A detailed derivation is performed for the boost-invariant 0+1D system to explicitly demonstrate the decomposition and the structure of the evolution equations.
   • 3+1D Generalization: The results are extended to general 3+1D spacetime using a projection operator formalism and interaction picture techniques (analogous to quantum mechanics) to handle non-commuting operators.

3. Key Contributions

A. Decomposition of the Exact Solution

The paper establishes that the exact solution to the Boltzmann moment equations naturally decomposes into:
$$\rho(\tau) = \underbrace{\rho_{\text{grad}}(\tau)}_{\text{Gradient Expansion}} + \underbrace{\rho_{\text{NP}}(\tau)}_{\text{Non-Perturbative}}$$

• Gradient Expansion ($\rho_{\text{grad}}$): Corresponds to the standard Chapman-Enskog expansion (divergent asymptotic series).
• Non-Perturbative Terms ($\rho_{\text{NP}}$): These are exponentially decaying terms ($\sim e^{-\xi}$) that encode initial conditions and free-streaming dynamics. They are non-analytic in the relaxation time parameter.

B. Structural Equivalence of Evolution Equations

A central finding is that the evolution equation for the gradient component ($\pi^G$) is structurally identical to the evolution equation for the full non-equilibrium component ($\pi$).

• The difference lies not in the form of the differential equations, but in the transport coefficients.
• The exact dynamics can be reproduced by a hydrodynamic equation where the transport coefficients (viscosities, relaxation times) are rescaled by factors dependent on the non-perturbative contributions (specifically involving incomplete gamma functions $\gamma(k, \xi)$).

C. Resolution of the "Hydrodynamization" Puzzle

The paper argues that hydrodynamics works far from equilibrium not because the system is close to equilibrium, but because:

1. Hydrodynamic equations inherently interpolate between free-streaming (collisionless) and collective (equilibrium) behaviors.
2. The "non-hydrodynamic" modes (often discarded in standard derivations) are not irrelevant; their effects are systematically absorbed into renormalized transport coefficients.
3. By adjusting these coefficients to account for initial data and free-streaming, standard hydrodynamic equations can accurately describe the system even when $Kn \sim 1$.

4. Key Results

• 0+1D Bjorken Flow:

  • Derived the exact solution for moments $\rho_{n,l}(\tau)$ showing explicit dependence on initial moments $\rho_{n,l}(\tau_0)$ and the damping factor $\xi = \int d\tau/\tau_R$.
  • Showed that the shear ($\pi$) and bulk ($\Pi$) pressure evolution equations derived from the exact solution match the Israel-Stewart form, provided the transport coefficients $\beta_\pi$ and $\beta_\Pi$ are modified to include terms like $e^{-\xi_0} \times (\text{initial anisotropy})$.
  • Demonstrated that the "gradient expansion" is actually a series of terms weighted by incomplete gamma functions, which smoothly transition from the free-streaming regime ($\xi \ll 1$) to the hydrodynamic regime ($\xi \gg 1$).

• 3+1D Generalization:

  • Proved that the decomposition into gradient and non-perturbative parts holds in 3+1D, even when operators do not commute.
  • Formulated a closure scheme where non-hydrodynamic moments are expressed as functions of hydrodynamic variables plus transient corrections.
  • Showed that the "Israel-Stewart" type theories are not merely phenomenological fixes for causality but emerge naturally from the exact moment hierarchy when non-perturbative modes are retained and their effects renormalized into the coefficients.

• Fixed Point Analysis:

  • Interpreted the dynamics in terms of two fixed points: the collisionless fixed point (free streaming) and the attractive fixed point (local equilibrium).
  • The gradient expansion only captures the behavior near the equilibrium fixed point. The full solution interpolates between these two, explaining why hydrodynamics (which includes the interpolation mechanism via modified coefficients) works early on.

5. Significance and Implications

• Theoretical Foundation: This work provides a first-principles derivation for the "renormalized transport coefficients" proposed in previous phenomenological studies. It explains why tuning transport coefficients in hydrodynamic simulations works so well: it effectively accounts for the missing non-perturbative initial data.
• Redefining Hydrodynamics: It reframes relativistic hydrodynamics not as a strict expansion around local equilibrium, but as an effective theory that interpolates between free-streaming and collective dynamics.
• Causality and Convergence: The non-perturbative terms are identified as essential for ensuring causality and the convergence of the theory, rather than being "non-hydrodynamic" noise to be discarded.
• Future Directions: The framework suggests that future hydrodynamic models should focus on systematically including these transient corrections (via rescaled coefficients) rather than just adding higher-order gradient terms, potentially unifying the description of small and large systems.

In summary, the paper resolves the paradox of hydrodynamics' success in far-from-equilibrium regimes by demonstrating that the exact kinetic evolution shares the same structural form as hydrodynamic equations, differing only by transport coefficients that dynamically encode the system's departure from equilibrium and its initial state.