Observable Dependence of Viscous Corrections in QGP: Heavy Quarks and Dileptons in Chapman--Enskog Theory

This paper presents the first calculation of heavy quark transport and thermal dilepton production in a QGP using second-order viscous corrections derived from the Chapman-Enskog expansion, revealing that these corrections significantly suppress drag forces and enhance early-time dilepton yields while demonstrating that observable modifications depend on the complex interplay between correction magnitude, momentum dependence, and the specific momentum weighting of each observable.

The Gist

Imagine a massive, ultra-hot soup made of the tiniest building blocks of the universe (quarks and gluons). Scientists call this "Quark-Gluon Plasma" (QGP). When heavy atoms smash together in giant particle colliders, they create this soup for a split second. The paper you're asking about tries to understand how this soup behaves when it's not perfectly calm, but rather "wobbly" and flowing with friction (viscosity).

Here is a simple breakdown of what the researchers did and found, using everyday analogies.

The Big Question: How do we measure the "wobble"?

Scientists know this soup expands and cools down very fast. To understand it, they use math to describe how particles move inside it. Usually, they assume the soup is in a perfect, calm state. But in reality, it's messy.

To fix this, scientists add "corrections" to their math to account for the messiness (viscosity). There are two main ways to do this:

1. The "Grad" Method: Think of this like drawing a smooth, simple curve to fit a messy set of dots. It's a standard, easy-to-use approximation.
2. The "Chapman-Enskog" (CE) Method: This is like a more detailed, step-by-step recipe that accounts for the messiness more precisely, looking at it in layers (first-order, then second-order).

The Goal: The authors wanted to see if using this more detailed "CE recipe" (up to the second layer of detail) changes the results compared to the standard "Grad" method. They tested this using two different "probes" (ways of measuring the soup).

Probe 1: The Heavy Quarks (The "Bowling Balls")

Imagine throwing a heavy bowling ball (a heavy quark) into a pool of water (the QGP).

• Drag: How much does the water slow the ball down?
• Diffusion: How much does the ball jitter and bounce around as it moves?

What they found:

• The "Grad" method and the "First-order CE" method gave somewhat similar results.
• The "Second-order CE" method (the super-detailed one) changed things significantly.
  • Drag: It made the water feel thicker to the bowling ball, slowing it down much more than the other methods predicted, especially at moderate speeds.
  • Jitter (Diffusion): It changed how the ball bounced sideways versus forward. The "second-order" math showed a complex pattern where the ball's movement depended heavily on its speed in a way the simpler methods missed.
• The Lesson: The detailed math didn't just add a little extra friction; it fundamentally changed how the heavy ball interacted with the soup, especially because the heavy ball "feels" the soup's particles in a specific speed range where the detailed math matters most.

Probe 2: The Thermal Dileptons (The "Ghost Messengers")

Now, imagine the soup is glowing and emitting light particles (dileptons) that pass right through the soup without getting stuck, like ghosts.

• Because they don't get stuck, they carry a perfect message from the moment they were created all the way to the detector.
• Scientists can look at these "ghosts" to see what the soup was like at different stages of its life (early hot stage vs. later cooling stage).

What they found:

• Early Times: When the soup is hottest and expanding fastest, the detailed "Second-order CE" math predicted a big burst of these "ghosts."
• Later Times: As the soup cools down, the difference between the "Grad" method and the "CE" method shrinks. They start to agree with each other.
• The Twist: Even though the "Grad" method is simpler, at very high speeds (high momentum), it actually predicted more ghosts than the detailed method.
• The Lesson: Just because the "CE" math says the soup is "messier" in the distribution of particles, it doesn't mean the final count of "ghosts" will always be higher. It depends on which part of the soup's speed range the "ghosts" are sensitive to.

The Main Takeaway: It's About the "Match"

The most important discovery in this paper is a concept the authors call "Observable Dependence."

Think of it like this:

• You have a Soup (the QGP).
• You have a Recipe (the math corrections: Grad vs. CE).
• You have a Taste Test (the observable: Heavy Quarks vs. Dileptons).

The paper shows that the Recipe doesn't change the Soup in a way that looks the same to every Taste Test.

• The Heavy Quark (bowling ball) is sensitive to the "middle-speed" particles in the soup. The detailed CE recipe changes the middle-speed particles the most, so the bowling ball feels a huge difference.
• The Dilepton (ghost) is sensitive to a wide range of speeds, including the very fast ones. The detailed CE recipe changes the fast particles differently than the simple Grad recipe, so the ghost count changes in a different pattern.

Conclusion:
You cannot just look at the math and say, "This correction is bigger, so the result must be bigger." You have to look at how the specific thing you are measuring (the probe) interacts with the specific part of the soup the math is changing.

The authors successfully calculated these effects for the first time using the detailed "Second-order" math. They found that while the math gets more complex, the results are "well-behaved" (they don't break or go crazy), but they do change our understanding of how heavy particles slow down and how light particles are emitted from the hot soup.

Technical Summary

Technical Summary: Observable Dependence of Viscous Corrections in QGP: Heavy Quarks and Dileptons in Chapman–Enskog Theory

Problem Statement
Relativistic heavy-ion collisions produce a Quark-Gluon Plasma (QGP) that behaves as a fluid with low shear viscosity ($\eta/s$). While the evolution of this system is successfully modeled using causal relativistic viscous hydrodynamics, deviations from local thermodynamic equilibrium are encoded in non-equilibrium corrections to particle distribution functions. Two primary methods exist to determine the form of these corrections: Grad's 14-moment approximation and the Chapman-Enskog (CE) expansion. While Grad's method has been widely used to study observables, the CE method—specifically up to second order in gradients—offers a systematic perturbative solution to the Boltzmann equation. A critical gap exists in understanding how these different viscous corrections, particularly second-order CE terms, propagate from the microscopic distribution functions to macroscopic observables. It is unclear whether the hierarchy of corrections observed in the distribution functions translates directly to experimental signals, as different observables may be sensitive to distinct regions of momentum space.

Methodology
The authors calculate heavy quark (HQ) transport coefficients and thermal dilepton production rates using viscous corrections up to second order in gradients derived from the Chapman-Enskog expansion of the Boltzmann transport equation in the relaxation time approximation. These results are compared against those obtained using Grad's 14-moment approximation.

The study employs the following framework:

• Hydrodynamic Evolution: The QGP expansion is modeled using second-order causal relativistic viscous hydrodynamics within a longitudinal boost-invariant Bjorken flow. Temperature and shear stress evolution profiles are obtained by numerically solving the dynamical equations with an initial temperature $T_0 = 0.5$ GeV and $\eta/s = 1/4\pi$.
• Heavy Quark Transport: HQ dynamics are described via the Fokker-Planck approach. The drag force ($A$) and momentum diffusion coefficients (transverse $B_0$ and longitudinal $B_1$) are calculated by integrating the collision kernel over the medium's distribution functions, which include first- and second-order CE corrections ($\delta f^{(1)}$ and $\delta f^{(2)}$) as well as Grad's ansatz ($\delta f^G$).
• Thermal Dilepton Production: The rate of thermal dilepton production is derived analytically up to second order in the CE expansion. The total yield is obtained by convoluting the production rate over the space-time history of the collision.
• Analysis: The authors systematically analyze the momentum dependence of the viscous corrections in the distribution functions and their subsequent impact on the transport and emission kernels. They examine how the interplay between the momentum structure of the corrections and the weighting of the observable kernels determines the final modification of the observables.

Key Contributions and Results

1. First-Order vs. Second-Order CE Corrections: This work constitutes the first study to include second-order viscous corrections for both heavy quark transport and thermal dilepton production.
2. Heavy Quark Transport:
   • Drag Force: The second-order CE corrections substantially suppress the drag force across the momentum range, a more significant effect than the first-order correction. In contrast, first-order corrections tend to enhance drag at high momenta. The suppression is driven by the strong overlap between the CE corrections and the transport kernel in the intermediate momentum region ($q \sim 1-2$ GeV).
   • Diffusion: For transverse diffusion ($B_0$), viscous effects enhance the coefficient at very low momenta and suppress it at higher momenta, with the total CE correction showing the largest modification. For longitudinal diffusion ($B_1$), viscosity suppresses the coefficient at low momenta but enhances it at high momenta ($p > 3$ GeV). Notably, the total CE contribution to $B_1$ remains subleading compared to Grad's ansatz, unlike the drag and transverse diffusion cases.
   • Spatial Diffusion ($D_s$): The inclusion of second-order CE terms leads to a large enhancement of the spatial diffusion coefficient, particularly near the transition temperature $T_c$.
3. Thermal Dileptons:
   • Time Evolution: CE corrections result in an enhanced contribution at early times (high temperature/strong gradients). As the QGP evolves and cools, these corrections decrease and converge toward the first-order CE correction.
   • Comparison with Grad: Unlike the HQ sector, the Grad's ansatz correction becomes dominant at high transverse momenta ($p_T > 1.5$ GeV) due to its quadratic momentum dependence, whereas CE corrections are well-behaved and remain convergent.
4. Observable Dependence: The study demonstrates that the modification of an observable is not governed solely by the magnitude of the viscous correction in the distribution function. Instead, it is determined by the interplay between the momentum dependence of the correction and the momentum weighting of the specific transport or emission kernel. For instance, HQ transport is sensitive to intermediate momenta where CE corrections are large, while dilepton production samples a broader range, including high-momentum tails where Grad's ansatz dominates.

Significance
The paper establishes that the "observable dependence" of viscous corrections is a critical factor in interpreting heavy-ion collision data. The authors conclude that the momentum structure of viscous corrections at the level of the distribution function is not directly translated to observables because different probes are sensitive to distinct regions of momentum space.

By systematically studying the propagation of non-equilibrium effects from microscopic distributions to macroscopic signals, the work provides a unified framework for understanding dissipative dynamics. It highlights that reliable connections between dissipative dynamics and measured QGP signals require accounting for how specific observables sample the non-equilibrium distribution, rather than assuming a universal hierarchy of corrections across all physical quantities. The results indicate that second-order CE corrections are essential for a precise description of heavy quark transport, while their impact on dileptons is time-dependent and distinct from Grad's approximation.