Collisional energy loss distribution of a fast parton in a hot or dense QCD medium

This paper computes the complete probability distribution (quenching weight) for the collisional energy loss of an ultrarelativistic parton traversing a quark-gluon plasma by solving a kinetic equation that resums arbitrary elastic scatterings and accounts for stochastic energy exchange, including thermal energy gain, within a framework applicable to various collision systems.

The Gist

The Big Picture: A Fast Runner in a Crowded Room

Imagine a super-fast runner (a "parton," which is a tiny piece of matter like a quark) sprinting through a crowded, hot room filled with people (a "quark-gluon plasma").

In the past, scientists mostly asked one question: "On average, how much distance does the runner lose because of bumping into people?" They calculated a single number, like "the runner slows down by 5 meters per second."

This new paper asks a much more detailed question: "What is the exact probability of the runner losing a specific amount of energy?"

Instead of just giving an average, the authors created a "quenching weight." Think of this as a weather forecast for the runner's energy. Instead of saying "it will rain 2 inches," they are saying, "there is a 10% chance of a drizzle, a 5% chance of a sudden downpour, and a 2% chance the runner actually gets a boost from a tailwind."

The Two Main Surprises

The paper reveals two things that standard "average" calculations miss:

1. The "Tailwind" Effect (Energy Gain)
Usually, we think of running through a crowd as only slowing you down. But because the room is hot and the people are jiggling around (thermal fluctuations), sometimes a person in the crowd might accidentally bump into the runner from behind, giving them a little push.

• The Paper's Claim: The authors calculated the probability that the runner actually gains energy. In their model, the runner can occasionally get a "free ride" from the thermal energy of the medium.

2. The "Rare Giant" Effect (Non-Gaussian Fluctuations)
If you flip a coin a million times, the results usually look like a smooth bell curve (a normal distribution). You rarely get 1,000 heads in a row.
However, in this crowded room, the runner mostly bumps into people gently. But very rarely, they might slam into a giant boulder (a "hard collision").

• The Paper's Claim: Because these rare, hard collisions happen, the energy loss doesn't follow a smooth bell curve. Instead, it follows a "skewed" distribution (like the famous Landau distribution). This means the runner is likely to lose a small amount of energy, but there is a significant chance they will lose a massive amount all at once due to one bad collision. The "average" calculation hides this danger.

How They Did It: The "Recipe"

To get these results, the authors had to mix two different ways of looking at the problem, like blending two types of flour:

1. The "Soft" Flour (HTL): For the gentle, frequent bumps with the crowd, they used a sophisticated mathematical tool called "Hard Thermal Loop" (HTL) resummation. This accounts for the fact that the crowd is a fluid that screens (blocks) some of the interactions.
2. The "Hard" Flour (Kinetic Theory): For the rare, violent crashes, they used standard kinetic theory, which treats the collisions like billiard balls hitting each other.

They created a smooth "seam" to glue these two methods together, ensuring the math works whether the collision is a gentle tap or a hard smash.

The "No-Scattering" Ghost

The paper also highlights a fascinating quirk depending on how long the runner stays in the room:

• In a Cold, Dense Room: If the room is cold and dense, there is a real chance the runner goes through without hitting anyone. The authors call this the "no-scattering" component. It's like a ghost in the distribution—a spike at zero energy loss.
• In a Hot Room: If the room is hot, the runner is guaranteed to interact with the jiggling crowd. The "ghost" disappears, replaced by a smeared-out probability that includes those rare energy boosts.

Why This Matters (According to the Paper)

The authors argue that for small systems (like collisions in smaller particle accelerators or the edges of big explosions), the "average" energy loss is a bad predictor. Because the path is short, the runner doesn't have time to experience enough collisions to smooth out the averages.

In these short trips, the fluctuations (the rare giant hits or the lucky tailwinds) are the most important part of the story. By providing the full probability distribution, this paper gives physicists a more accurate tool to predict what happens to particles in these complex, chaotic environments.

Summary in One Sentence

This paper replaces the simple "average speed loss" of a particle moving through hot matter with a detailed "probability map" that accounts for rare, massive energy crashes and the surprising possibility of the particle actually gaining energy from the heat of the crowd.

Technical Summary

1. Problem Statement

The paper addresses a significant gap in the phenomenology of jet quenching in heavy-ion and small-system collisions. While the average collisional energy loss ($\Delta E_{coll}$) of a fast parton traversing a Quark-Gluon Plasma (QGP) has been extensively studied (e.g., by Bjorken and subsequent HTL resummation calculations), the probability distribution of this energy loss (the "quenching weight") has not been derived from first principles.

• Motivation: The average loss is often insufficient for realistic applications, particularly in small systems (high-multiplicity $pp$, $pPb$, peripheral $AA$) or for heavy quarks, where path lengths ($L$) are short and fluctuations are significant.
• The Gap: Unlike radiative energy loss, where probabilistic frameworks (quenching weights) exist, collisional loss has traditionally been treated only via its first moment. A full probability distribution $f(x, \Delta)$ is required to account for:
  • Event-by-event fluctuations (straggling).
  • The possibility of energy gain ($\Delta < 0$) via absorption of thermal gluons.
  • Non-Gaussian behavior arising from rare, hard collisions.

2. Methodology

The authors derive the collisional quenching weight $f(x, \Delta)$ using a rigorous kinetic theory approach combined with perturbative QCD (pQCD) and Hard Thermal Loop (HTL) resummation.

A. Kinetic Equation Formulation

The probability density $f(x, \Delta)$ for a parton with initial energy $E_i$ to lose energy $\Delta$ after traveling a distance $x$ is governed by a linear kinetic equation (analogous to the Landau equation for ionization):
$$\frac{\partial f}{\partial x} = \int d\epsilon \, w(\epsilon) \left[ f(x, \Delta - \epsilon) - f(x, \Delta) \right]$$
where $w(\epsilon)$ is the differential scattering rate per unit distance.

• Solution: The equation is solved formally using a Laplace transform, yielding an integral representation involving a Bromwich contour:
  $$f(x, \Delta) = \frac{1}{2\pi i} \int_B d\nu \, \exp\left( \nu \Delta - x I(\nu) \right)$$
  where $I(\nu) = \int d\epsilon \, w(\epsilon) (1 - e^{-\nu \epsilon})$.
• Physical Interpretation: The solution represents a Poisson process of independent elastic scatterings. The term $e^{-x\Gamma}$ (where $\Gamma$ is the total damping rate) represents the probability of "no scattering," while the series expansion accounts for $n$-fold scatterings.

B. Calculation of the Scattering Rate $w(\epsilon)$

The core technical challenge is determining the differential rate $w(\epsilon)$ across all energy transfer scales. The authors split the calculation into two regimes and match them:

1. Soft Domain ($|\epsilon| \lesssim m_D$):

   • Standard perturbation theory fails due to infrared divergences.
   • HTL Resummation: The authors use Hard Thermal Loop (HTL) effective theory to resum soft gluon exchanges.
   • They derive the functions $F_T(\epsilon)$ and $F_L(\epsilon)$ (transverse and longitudinal contributions) using the spectral densities of the HTL gluon propagator.
   • A key innovation is reformulating the $q$-integral using Cauchy's theorem to sum over the poles of the thermal dispersion relations (plasmon modes), providing a compact, numerically stable expression.

2. Hard Domain ($|\epsilon| \gtrsim \Lambda$):

   • For large energy transfers, standard kinetic theory with exact $2 \to 2$ scattering kinematics is used.
   • The rate is calculated using Born-level matrix elements for scattering off thermal gluons and light quarks.

3. Matching:

   • The soft (HTL) and hard (Born) results are matched using a multiplicative scheme to ensure a smooth transition across the intermediate scale.
   • The final rate $w(\epsilon)$ includes the Bose-Einstein distribution $n_B(\epsilon)$, allowing for $\epsilon < 0$ (energy gain).

C. Numerical Evaluation

The authors perform a numerical inversion of the Laplace transform (Eq. 2.7) using the Bromwich contour along the imaginary axis. They handle the oscillatory nature of the integrand using specialized techniques (Levin u-transformation) to ensure convergence.

3. Key Contributions

1. First-Principles Derivation: This is the first derivation of the full collisional quenching weight $f(x, \Delta)$ from QCD first principles, bridging the gap between average loss calculations and probabilistic jet-quenching models.
2. Inclusion of Energy Gain: The formalism naturally incorporates the possibility of the fast parton gaining energy from the medium ($\Delta < 0$), a feature absent in standard "energy loss only" models.
3. Finite Path Length Effects: The study explicitly treats finite path lengths $x$, revealing regimes where the distribution is far from the asymptotic Landau limit.
4. Handling Infrared Divergences: The authors demonstrate that while the total scattering rate $\Gamma$ diverges in the HTL framework at finite temperature (due to unscreened magnetic modes), the probability distribution $f(x, \Delta)$ remains infrared safe and well-defined. The divergence of $\Gamma$ merely suppresses the "no-scattering" $\delta(\Delta)$ term, replacing it with a singular but integrable distribution at $\Delta=0$.

4. Key Results

A. The Landau Limit ($x \to \infty$)

For sufficiently large path lengths, the distribution converges to the universal Landau distribution, regardless of the specific medium parameters ($T, \mu$). This confirms that the long-distance behavior is dominated by the $1/\epsilon^2$ tail of the scattering rate (rare hard collisions), consistent with the Central Limit Theorem's failure for power-law distributions.

B. Finite Path Length Regimes

• Cold Dense Medium ($T=0, \mu \neq 0$):
  • The distribution contains a distinct Dirac-$\delta$ component at $\Delta=0$ (probability of no scattering), proportional to $e^{-x\Gamma_0}$.
  • For small $x$, the distribution is dominated by single and double scattering terms.
  • As $x$ increases, the distribution broadens and approaches the Landau shape.
• Thermal Medium ($T > 0$):
  • Smearing of the $\delta$-term: Thermal fluctuations "smear" the zero-scattering peak into a singular distribution $f(x, \Delta) \sim |\Delta|^{-1 + \alpha x T}$ near $\Delta=0$.
  • Energy Gain: The distribution extends to $\Delta < 0$. For small $x$ and $T$, the probability of energy gain is significant.
  • Critical Transition: There is a transition in the behavior at $\Delta=0$ depending on the dimensionless parameter $\bar{x} \alpha_s C_F T$. If $x$ is small enough, $f(x, \Delta)$ is singular at $\Delta=0$; if $x$ is large, it becomes finite.

C. Non-Gaussian Behavior

The distributions are markedly non-Gaussian.

• Large $\Delta$: Dominated by rare, hard scatterings (power-law tails).
• Small $\Delta$: Dominated by soft interactions and thermal fluctuations.
• Most Probable vs. Average: The most probable energy loss ($\Delta_p$) scales logarithmically with $x$, while the average energy loss ($\langle \Delta \rangle$) is much larger due to the heavy tail of the distribution. This highlights the inadequacy of using only the average loss for phenomenology.

5. Significance and Implications

• Phenomenology of Small Systems: The results are crucial for interpreting data from high-multiplicity $pp$ and $pPb$ collisions, where path lengths are short, and fluctuations are not averaged out. Standard quenching models based on average loss may fail in these regimes.
• Heavy Quark vs. Light Parton: While derived for a heavy quark (to simplify tagging), the results apply to light partons in the high-energy limit ($E_i \to \infty$), with minor modifications for exchange diagrams.
• Theoretical Robustness: The work demonstrates that collisional energy loss can be treated consistently within pQCD even in the presence of infrared divergences in the total rate, provided one focuses on the probability distribution.
• Future Applications: The authors provide a numerical code and the functional form of the quenching weight, enabling more accurate calculations of nuclear modification factors ($R_{AA}$) and heavy-flavor suppression that account for stochastic fluctuations and energy gain.

In summary, this paper provides a foundational theoretical framework for the stochastic nature of collisional energy loss, moving beyond mean-field approximations to a complete probabilistic description essential for precision QCD phenomenology in both large and small collision systems.