Heavy Quark Transport is Non-Gaussian Beyond Leading Log

This paper demonstrates that heavy quark transport in weakly coupled non-Abelian plasmas is intrinsically non-Gaussian with asymmetric exponential tails beyond the leading logarithm, a robust feature shared with strongly coupled holographic plasmas that is essential for accurate equilibration dynamics.

The Gist

The Big Picture: A Heavy Quark in a Hot Soup

Imagine the universe just after the Big Bang, or the center of a heavy-ion collision at a particle accelerator like the LHC. It's a super-hot, super-dense soup of energy called Quark-Gluon Plasma (QGP).

Inside this soup, we drop in a Heavy Quark (like a charm or bottom quark). Think of this heavy quark as a giant, bowling ball trying to roll through a crowd of tiny, hyper-active ping-pong balls (the lighter particles in the plasma).

For decades, physicists have tried to predict how that bowling ball moves. The standard way to describe this movement is using a "Gaussian" model (also known as Brownian motion).

• The Gaussian Idea: Imagine the bowling ball gets hit by the ping-pong balls. Most hits are tiny nudges from the side. Occasionally, a slightly harder hit comes. The math says these hits follow a "Bell Curve." Most of the time, the ball moves a little bit; very rarely, it moves a lot. The path is smooth and predictable, like a drunk person stumbling in a straight line.

The Problem: The "Bell Curve" Lie

This paper says: The Bell Curve is wrong.

While the Bell Curve works for small nudges, it fails to describe what happens when the heavy quark moves fast or when we look at the rare but massive hits.

The authors discovered that the distribution of hits isn't a smooth Bell Curve. Instead, it looks like a Bell Curve with long, lopsided "tails."

• The Core: In the middle, yes, it looks like a Bell Curve (many small nudges).
• The Tails: But on the sides, the curve doesn't drop off quickly like a Bell Curve. Instead, it has exponential tails. This means there is a much higher chance of the heavy quark getting hit by a "monster" ping-pong ball that sends it flying in a specific direction than the old math predicted.

The Analogy: The Drunk vs. The Rollerblader

The Old View (Gaussian/Langevin):
Imagine a drunk person walking home. They stumble left and right randomly. If you watch them for a long time, their path looks like a smooth, predictable drift. If they stumble, it's usually a small step. This is the "Gaussian" view.

The New View (Non-Gaussian):
Now imagine that same drunk person, but every now and then, a giant (a high-energy particle) runs into them and knocks them 20 feet to the left.

• The "Bell Curve" says giant knocks are impossible.
• The "Exponential Tail" says giant knocks happen more often than you think, and they are asymmetric. The giant might knock the person hard to the left, but rarely to the right.

This asymmetry is crucial. If you only use the "Bell Curve" math, you can't explain how the person eventually stops moving and settles down (equilibrates). The "giant knocks" are actually the secret sauce that helps the system reach balance.

How They Found This

The scientists did a massive calculation. They had to combine two different ways of looking at the physics:

1. The "Hard" View: Looking at the high-energy, fast collisions (like the giant knocks).
2. The "Soft" View: Looking at the low-energy, gentle nudges (the ping-pong balls).

They had to be very careful not to double-count the area where these two views overlap. Once they cleaned up the math, they isolated the "fixed-order" kernel (the core rulebook of how the quark moves).

When they looked at the math for the "tails" of the distribution, they found something surprising: The math showed exponential tails.

The "Universal" Discovery

Here is the most mind-blowing part.

• This "Exponential Tail" structure was previously found in Strongly Coupled systems (using a theory called String Theory/Holography).
• This paper proves it also exists in Weakly Coupled systems (standard particle physics).

The Metaphor:
Imagine you are studying how water flows.

• In one lab, you study water flowing through a thick, sticky honey (Strong Coupling). You find the flow has a weird, jagged pattern.
• In another lab, you study water flowing through a thin, watery stream (Weak Coupling). You expect it to be smooth.
• The Discovery: You find that both the honey and the stream have the exact same jagged pattern.

This suggests that this "Gaussian Core + Exponential Tail" structure is a fundamental law of nature for how heavy particles move through hot plasma, regardless of whether the plasma is "sticky" or "thin," or whether it follows the rules of supersymmetry.

Why Does This Matter?

1. Fixing the Math: For a long time, physicists were stuck because the "Bell Curve" math violated a fundamental rule of physics (the Einstein relation) when the quark was moving fast. The new "Non-Gaussian" math fixes this. It shows that the "giant knocks" (the tails) are necessary to keep the laws of physics consistent.
2. Better Experiments: When scientists at the LHC or RHIC look at data from heavy-ion collisions, they use models to interpret what they see. If they use the old "Bell Curve" models, they might be misinterpreting the data.
3. Real-World Application: This helps us understand the "viscosity" and "temperature" of the early universe. By understanding how these heavy quarks actually move (with their exponential tails), we can build better models of the Quark-Gluon Plasma, which is the state of matter that existed microseconds after the Big Bang.

Summary

The paper proves that heavy particles moving through hot plasma don't just take small, random steps. They occasionally get hit by massive, asymmetric "kicks" that are far more common than previously thought. This "Non-Gaussian" behavior is a universal feature of the universe, appearing in both weak and strong interactions, and it is essential for the system to reach a state of balance.

Technical Summary

1. Problem Statement

Heavy quarks (HQs) serve as critical probes for the Quark-Gluon Plasma (QGP) created in ultra-relativistic heavy-ion collisions. Traditionally, their transport is modeled using Langevin or Fokker-Planck equations, which assume Gaussian momentum transfer statistics. This assumption relies on the Einstein relation ($\kappa_L = 2MT\eta_D$) between the drag coefficient ($\eta_D$) and the longitudinal diffusion coefficient ($\kappa_L$).

However, it has long been known that in weakly coupled plasmas, this Einstein relation is violated at non-zero velocities once calculations go beyond the leading logarithmic (leading-log) accuracy. Similarly, violations occur in strongly coupled $\mathcal{N}=4$ Super Yang-Mills (SYM) theory. This creates a theoretical obstruction: first-principles calculations (perturbative QCD or lattice QCD) cannot be consistently mapped to phenomenological Langevin descriptions if the underlying momentum transfer is not Gaussian. While recent work established that heavy quark evolution is governed by a Kolmogorov equation with a non-Gaussian kernel in strongly coupled theories, it remained unknown whether this non-Gaussian structure exists in weakly coupled QCD.

2. Methodology

The authors address this by computing the leading-order momentum transfer kernel for relativistic heavy quarks in weakly coupled non-Abelian plasmas. Their methodology involves several rigorous steps:

• Framework: They work in the heavy mass limit ($M \gg T$) using Heavy Quark Effective Theory (HQET), neglecting $1/M$ corrections. The momentum transfer distribution $P(k; v)$ is derived from the Fourier transform of a Wilson loop expectation value $\langle W_v \rangle$ on the Schwinger-Keldysh contour.
• Matching Hard and Soft Physics: To obtain a consistent leading-order result ($O(g^4)$), they match the perturbative "hard" momentum transfer ($p \sim T$) with the Hard Thermal Loop (HTL) resummed "soft" physics ($p \sim gT$). They employ a subtraction scheme to remove double-counting:
  $$\text{Result} = (\text{Pert.}) + (\text{HTL}) - (\text{HTL, unresummed})$$
• Isolating the Fixed-Order Kernel: A major technical challenge is that the standard HTL-resummed kernel contains UV-divergent terms when taking derivatives (cumulants) at $n \geq 3$, which obscures the fixed-order structure. The authors resolve this by introducing Heaviside step functions to separate momentum regimes ($|p| > T$ and $|p| < T$). They explicitly isolate the strict fixed-order kernel ($S_{FO}$) by removing the HTL-resummed contributions in the hard region and carefully expanding the soft region integrals to extract the correct $g^4 \ln(1/g)$ terms.
• Legendre Transform: They compute the longitudinal evolution kernel $K(x, v)$ (the cumulant generating function) and perform a Legendre transform to obtain the momentum transfer probability distribution $\tilde{S}(C_L; v)$.
• Analytic Continuation: They analyze the convergence properties of the cumulant series to determine the behavior of the distribution tails.

3. Key Contributions

• Derivation of the Fixed-Order Kernel: The paper provides the first complete, consistent calculation of the strict fixed-order momentum transfer kernel in weakly coupled QCD, free from the UV divergences that plagued previous attempts to extract higher-order cumulants.
• Proof of Non-Gaussianity: They demonstrate that beyond leading-log accuracy, the longitudinal momentum transfer distribution is intrinsically non-Gaussian. It consists of a Gaussian core (governed by the Central Limit Theorem) and asymmetric exponential tails.
• Universality of Structure: The authors show that this specific structure (Gaussian core + exponential tails) is identical to that found in strongly coupled $\mathcal{N}=4$ SYM. This suggests the feature is robust and independent of coupling strength, conformality, or supersymmetry.
• Resolution of the Einstein Relation Violation: They clarify that the violation of the Einstein relation is not a failure of the theory but a consequence of the non-Gaussian tails. The generalized equilibration condition (derived in Ref. [30]) is satisfied only when these tails are included.

4. Key Results

• Distribution Shape: The longitudinal momentum transfer probability $P(k_L)$ exhibits:
  • A Gaussian core near $k_L = 0$.
  • Asymmetric exponential tails for large $|k_L|$. Specifically:
    $$P(k_L) \propto e^{-R(v)k_L} \quad \text{as } k_L \to \infty$$
    $$P(k_L) \propto e^{+[R(v)+v/T]k_L} \quad \text{as } k_L \to -\infty$$
    Here, $R(v)$ is a "volatility exponent" determined by the radius of convergence of the cumulant series.
• Role of Hard Physics: The exponential tails are driven by hard physics ($|p| \sim T$), not soft HTL physics. The high-order cumulants ($n \geq 3$) are dominated by the perturbative hard contribution, which dictates the radius of convergence.
• Coupling Independence: The leading contribution to the volatility exponent $R(v)$ is independent of the coupling constant $g$ (or 't Hooft coupling $\lambda$). Coupling dependence only appears at subleading orders.
• Numerical Verification: Figure 1 in the paper shows that when normalized by the first two cumulants (drag and width), the distributions for weakly coupled QCD, weakly coupled $\mathcal{N}=4$ SYM, and strongly coupled $\mathcal{N}=4$ SYM all collapse onto the same universal non-Gaussian curve, deviating significantly from a Gaussian parabola, especially at high velocities ($v=0.95$).

5. Significance

• Theoretical Consistency: This work resolves a long-standing tension between first-principles QCD calculations and phenomenological transport models. It proves that a purely Gaussian (Langevin) description is insufficient for relativistic heavy quarks in QGP.
• Phenomenological Impact: Since non-Gaussian dynamics can alter heavy-flavor observables and the transport coefficients extracted from them, current phenomenological models based on Gaussian truncations may be systematically biased.
• New Framework: The authors propose a more realistic description of HQ transport characterized by three parameters: the drag $\eta_D(v)$, the diffusion $\kappa_L(v)$, and the volatility exponent $R(v)$, which encodes the probability of rare, large momentum kicks.
• Future Directions: The results suggest that future analyses of heavy-ion collision data (e.g., from LHC Run 3 and the proposed ALICE 3 upgrade) should utilize non-Gaussian kernels to accurately extract transport coefficients and compare them with lattice QCD determinations.

In summary, the paper establishes that heavy quark transport in QGP is fundamentally non-Gaussian due to asymmetric exponential tails arising from hard scattering processes, a feature that is universal across both weak and strong coupling regimes.