Non-Gaussian hydrodynamic fluctuations in an expanding relativistic fluid

This paper derives analytical evolution equations for two- and three-point velocity correlators in a boost-invariant relativistic fluid using effective field theory, demonstrating that the Landau frame is optimal for studying non-Gaussian fluctuations and revealing that three-point correlations exhibit nonlinear memory effects dependent on two-point dynamics, which are crucial for the search for the QCD critical point.

The Gist

Imagine the universe's most extreme substance, the quark-gluon plasma, as a super-hot, super-dense soup created when heavy atoms smash together. This soup doesn't just sit still; it explodes outward, expanding and cooling down incredibly fast, much like a giant, invisible balloon being blown up at the speed of light.

This paper is about understanding the tiny, chaotic ripples that happen inside this expanding soup.

The Big Picture: Why Ripples Matter

Scientists are trying to find a "critical point" in the laws of physics (the QCD phase diagram). Think of this like finding the exact moment water turns into ice. Near this turning point, the soup doesn't just have tiny, random ripples; it starts having wild, unpredictable, and interconnected waves.

Usually, scientists study these ripples assuming they are simple, "Gaussian" (bell-curve) fluctuations—like gentle waves on a calm lake. But near that critical point, the waves get messy and "non-Gaussian." They start talking to each other in complex ways. The authors of this paper wanted to figure out exactly how these messy, non-Gaussian waves evolve as the soup expands.

The Problem: A Moving Target

To study these waves, you need a reference frame (a way to measure time and space).

• The Old Way: Imagine trying to measure the waves on a boat that is rocking wildly. If you try to measure the "average" direction of the boat while it's rocking, your measurements get messy. In physics terms, if you try to define time based on the fluctuating fluid itself, you get mathematical "glitches" (singularities) when looking at complex, three-way interactions between waves.
• The New Way (The Paper's Solution): The authors suggest using a "steady hand" approach. Instead of measuring time based on the wobbly, fluctuating fluid, they measure it based on the average flow of the soup.
  • Analogy: Imagine a river. The water is churning with eddies and whirlpools (fluctuations). If you try to time the river's flow by watching a single leaf spinning in a whirlpool, your clock is useless. But if you time the flow based on the average current of the whole river, your clock is steady.
  • In the specific case of this expanding soup (Bjorken flow), this "average current" frame turns out to be the same as the "density frame." This makes the math much cleaner and avoids the glitches.

The Discovery: Memory and Non-Linear Dance

The authors derived equations to track how these ripples change over time. They focused on two types of ripples:

1. Two-point correlations: How one ripple relates to another (like two waves interacting).
2. Three-point correlations: How three ripples interact together (the messy, non-Gaussian stuff).

Key Finding 1: The "Memory" Effect
The paper shows that the three-way interactions (the messy waves) don't just happen in a vacuum. They are driven by the history of the two-way interactions.

• Analogy: Imagine a dance floor. The "two-point" dancers are the main couples. The "three-point" dancers are the wild, chaotic group dances. The paper shows that the group dances are remembering how the couples danced earlier. The chaotic behavior isn't random; it's a direct, non-linear consequence of the simpler waves evolving over time. Because the soup is expanding so fast, the waves don't have time to settle down, so they carry this "memory" of their past movements.

Key Finding 2: The Universal Fade-Out
Even though these waves get wild and complex in the middle of the expansion, the authors found a pattern for how they eventually calm down.

• Analogy: No matter how wild the party gets in the middle of the night, if you wait long enough, everyone eventually goes home. The paper shows that as the soup expands and cools, these complex three-way ripples eventually die out and return to zero (equilibrium), following a specific, predictable mathematical rule (a power law).

The Method: A New Toolkit

To do this, the authors used a sophisticated mathematical toolkit called Effective Field Theory (EFT).

• Analogy: Think of EFT as a "recipe book" for fluids. Instead of trying to track every single atom (which is impossible), the recipe tells you how to mix the ingredients (energy, pressure, viscosity) to predict the behavior of the whole pot. The authors wrote a new "recipe" specifically for these messy, non-Gaussian ripples in an expanding universe.

Summary

In short, this paper provides a new, cleaner way to calculate how complex, chaotic waves behave in an expanding, cooling fluid.

1. They fixed a mathematical problem by choosing the right "viewpoint" (the average flow).
2. They showed that complex waves are driven by the history of simpler waves (memory).
3. They provided exact formulas (solutions) for how these waves evolve and eventually fade away as the fluid expands.

This work is a crucial step for scientists who want to use these "ripples" as a detector to find the elusive QCD critical point in future experiments, helping them distinguish between normal chaos and the special chaos that happens near a phase transition.

Technical Summary

1. Problem Statement

The paper addresses the need for a quantitative understanding of non-Gaussian hydrodynamic fluctuations in relativistic fluids undergoing rapid expansion, specifically in the context of heavy-ion collisions and the search for the QCD critical point.

• Context: While the QCD critical point is expected to enhance thermal fluctuations, the quark-gluon plasma (QGP) created in heavy-ion collisions expands and cools rapidly. This dynamical evolution prevents many fluctuation modes from reaching thermal equilibrium before "freeze-out" (hadronization).
• The Gap: Existing theoretical predictions often assume thermal equilibrium. However, non-equilibrium effects, such as memory (where the system's history influences current fluctuations), are crucial for accurate predictions of critical point signatures.
• Specific Challenge: Deriving deterministic evolution equations for three-point (and higher) correlators is technically difficult because the fluid velocity itself fluctuates. In standard fluctuating frames (like the Landau frame), this leads to ambiguities and singular time derivatives of stochastic noise in higher-order correlators.

2. Methodology

The authors employ the Effective Field Theory (EFT) framework for fluctuating hydrodynamics, specifically utilizing the Schwinger-Keldysh formalism.

• Framework: They use an effective action principle that incorporates the dynamical Kubo-Martin-Schwinger (KMS) symmetry to ensure the fluctuation-dissipation theorem is satisfied locally.
• Choice of Frame (Key Innovation): To resolve the ambiguity of time evolution in fluctuating frames, the authors work in the Average Landau Frame.
  • In this frame, the time direction is defined by the average fluid velocity ($\bar{u}^\mu$) rather than the fluctuating local velocity.
  • For Bjorken flow (boost-invariant expansion), this frame coincides with the Density Frame.
  • Benefit: This choice absorbs the time-like components of the noise into the definition of energy and momentum fluxes, ensuring that the remaining stochastic noise is purely spatial. This eliminates singular time derivatives of noise in the evolution equations for higher-order correlators.
• Derivation Steps:
  1. Construct the effective Lagrangian up to cubic order in fluctuations (energy density $\delta\epsilon$ and momentum density $\delta\pi$) in the Density Frame.
  2. Derive stochastic equations of motion (Langevin-like) and convert them into deterministic Schwinger-Dyson equations for equal-time correlation functions.
  3. Perform a Wigner transform to move to momentum space, separating slow (relaxational) and fast (oscillatory) modes.
  4. Solve the resulting hierarchy of differential equations analytically for the Bjorken background.

3. Key Contributions

• Resolution of Frame Ambiguity: The paper rigorously demonstrates that the Average Landau (Density) frame is the superior choice for studying non-Gaussian velocity fluctuations in relativistic systems, avoiding the mathematical singularities associated with time derivatives of noise in other frames.
• Analytical Solutions for Three-Point Functions: The authors derive and solve the evolution equations for three-point correlators ($G_3$) in a Bjorken expanding background. This is the first analytical treatment of non-Gaussian dynamics in this specific expanding geometry.
• Hierarchy of Scales: They explicitly show how the evolution of three-point functions is driven by the non-equilibrium dynamics of two-point functions. The three-point correlators act as a "memory" of the system's departure from equilibrium.
• Mode Decoupling: They identify that in the Bjorken background, only transverse momentum modes are "slow" (relaxational) for three-point functions, while longitudinal modes oscillate rapidly and average out.

4. Key Results

• Evolution Equations: The paper provides explicit evolution equations (Eqs. 3.11, 3.12) showing that the time derivative of an $N$-point function depends linearly on itself but nonlinearly on lower-order functions (specifically products of two-point functions).
  • $\partial_\tau G_3 \sim \text{Linear}(G_3) + \text{Nonlinear}(G_2 \times G_2)$.
• Analytical Solutions:
  • Two-point functions: Solutions show exponential decay of initial conditions and a universal late-time approach to equilibrium governed by a power law dependent on the dimensionless parameter $q^2 \eta(\tau) \tau$.
  • Three-point functions: The solutions consist of an initial condition term (decaying) and a dynamical source term driven by the product of two-point functions.
• Memory Effects: The three-point correlators are dynamically generated by the system falling out of equilibrium. Even if initialized at equilibrium values (zero), the expansion causes the two-point functions to deviate from equilibrium, which in turn sources non-zero three-point functions.
• Late-Time Behavior:
  • Three-point correlators eventually relax to zero (as required by isotropy in equilibrium).
  • The relaxation follows a universal power-law scaling determined by the effective relaxation rate $q^2 \eta(\tau) \tau$.
  • The transition to this universal behavior occurs when $\sqrt{q^2 \eta \tau} \sim O(1)$.
• Numerical Illustration: The authors provide numerical plots (Figs. 1 & 2) showing that longer wavelength modes (smaller $q$) relax more slowly and exhibit larger magnitudes of non-Gaussianity compared to short-wavelength modes.

5. Significance and Outlook

• Phenomenological Application: These analytical results provide the necessary input for the Maximum-Entropy Freeze-Out framework. This allows theorists to convert hydrodynamic fluctuations into hadronic fluctuations (observable in experiments) with controlled non-Gaussian corrections.
• QCD Critical Point Search: The findings are critical for the Beam Energy Scan program at RHIC. Understanding how non-Gaussian fluctuations evolve in a non-equilibrium, expanding medium is essential for distinguishing signals of the QCD critical point from background hydrodynamic noise.
• Future Directions: The paper sets the stage for extending these methods to include conserved baryon charges, more general expanding backgrounds, and critical dynamics near the phase transition.

In summary, this work bridges the gap between microscopic fluctuation theory and macroscopic hydrodynamic evolution, providing the first analytical toolkit to quantify how non-Gaussian fluctuations evolve and relax in the expanding fireball of a heavy-ion collision.