Chapman-Enskog calculation of the shear viscosity of quark-gluon plasma including all $2\leftrightarrow 2$ scatterings at finite temperature

This paper applies the Chapman-Enskog method to calculate the shear viscosity of quark-gluon plasma using perturbative-QCD cross sections with scaled thermal masses, demonstrating that a specific screening parameter ($\kappa=0.4$) aligns results with the AMY framework and revealing the ratio's strong sensitivity to the momentum scale choice, which yields $\eta/s \sim 0.15$ at the phase transition temperature.

The Gist

Imagine a pot of soup so hot and dense that the individual ingredients—quarks and gluons—stop behaving like distinct particles and instead melt into a chaotic, super-hot fluid called Quark-Gluon Plasma (QGP). This is the state of matter that existed just microseconds after the Big Bang.

The scientists in this paper, Okey Ohanaka and Zi-Wei Lin, are trying to figure out how "sticky" or "thick" this cosmic soup is. In physics, this stickiness is called shear viscosity. Think of it like the difference between honey and water: honey has high viscosity (it resists flowing), while water has low viscosity (it flows easily).

Here is the simple breakdown of what they did and what they found:

1. The Problem: Too Many Collisions

To understand how thick this soup is, you have to watch how the particles bump into each other. In this soup, particles are constantly colliding.

• The Old Way: Previous methods (like the "AMY framework") were like using a very complex, high-tech calculator that accounts for every tiny detail of the universe's rules. It's accurate but hard to use for other types of simulations.
• The New Way: The authors used a different mathematical tool called the Chapman-Enskog method. Think of this as a "general recipe" they recently wrote down. This recipe allows them to calculate the thickness of the soup based on any type of collision rule you give them, not just the specific ones used in the old method.

2. The "Screening" Issue: Fixing the Math Glitches

When they tried to use their new recipe with the standard rules of particle physics (perturbative QCD), the math started to break.

• The Glitch: In the real world, particles have a "personal space" (thermal mass) that stops them from getting infinitely close. In the math, if you don't account for this, the numbers can go crazy—becoming negative (which is impossible for a collision rate) or infinitely large.
• The Fix: The authors added a "screening" filter to the math. Imagine putting a safety net under a trapeze artist. They adjusted the math so the particles couldn't get too close, preventing the numbers from crashing.
• The Tuning Knob ($\kappa$): They found that using the standard safety net (where the net is exactly the size of the particle's personal space) made their results too high compared to the old, trusted methods. So, they introduced a "tuning knob" called $\kappa$. By turning this knob down to 0.4, they made their new, simpler recipe match the results of the complex, trusted old method perfectly.

3. The "Speed Limit" Choice ($Q$)

In their calculations, they had to choose a "speed limit" for how fast the particles are moving when they collide. This is called the momentum scale ($Q$).

• They found that this choice is like choosing the zoom level on a camera. If you zoom in too much or too little, the picture of the viscosity changes drastically.
• They discovered that choosing a specific zoom level ($Q = 3T$, where $T$ is temperature) gives a very specific result: at the moment the universe cooled down enough for normal matter to form (the phase transition), the plasma was surprisingly thin.
• The Result: The ratio of stickiness to disorder (viscosity/entropy) was about 0.15. This is very close to the theoretical "perfect fluid" limit (0.08), meaning this cosmic soup flows almost as easily as possible.

4. Why the "Extra Fixes" Didn't Matter Much

The authors had to add extra mathematical "patches" to make sure the collision numbers were always positive and finite (not infinite).

• The Surprise: They expected these patches to change the final result a lot. However, they found the patches barely changed the final viscosity.
• The Reason: The "stickiness" of the soup is mostly determined by collisions where particles hit each other with moderate energy. The patches mostly fixed the math for collisions where particles barely touched (very low energy). Since those low-energy collisions don't contribute much to the overall "stickiness," fixing them didn't change the final answer.

Summary

The paper provides a new, flexible "recipe" (the Chapman-Enskog method) for calculating how thick the early universe's soup was. They fixed some mathematical glitches by adding a safety net and a tuning knob. They found that with the right settings, their simple recipe matches the complex, trusted methods, and it suggests that the early universe's plasma was an incredibly smooth, low-viscosity fluid. This new recipe is now ready to be used by other scientists to simulate how this plasma behaves in computer models.

Technical Summary

1. Problem Statement

The shear viscosity ($\eta$) of the Quark-Gluon Plasma (QGP) is a critical transport coefficient for understanding the hydrodynamic evolution of heavy-ion collisions. While the shear viscosity has been calculated using the AMY (Arnold-Moore-Yaffe) framework, which employs perturbative QCD (pQCD) with parton self-energies and includes both $2 \leftrightarrow 2$ and $1 \leftrightarrow 2$ processes, there is a need for a general analytical expression that relates shear viscosity directly to arbitrary parton scattering cross sections.

This is particularly important for parton transport models (e.g., ZPC, AMPT, BAMPS) and QCD effective kinetic theory, where specific cross-section inputs are used. Existing analytical expressions often lack generality or do not account for the specific screening mechanisms (thermal masses vs. self-energies) used in different frameworks. Furthermore, standard pQCD cross sections at finite temperature often diverge or become unphysical (negative) at low energies, requiring regularization that must be carefully mapped to transport coefficients.

2. Methodology

The authors employ the Chapman-Enskog (CE) method to derive and apply a general analytical expression for shear viscosity.

• General Analytical Expression: They utilize a previously derived formula for the shear viscosity of a massless multi-species quark-gluon gas in chemical equilibrium with Boltzmann statistics. This expression accounts for all $2 \leftrightarrow 2$ elastic and inelastic scatterings (including $gg \leftrightarrow gg$, $gq \leftrightarrow gq$, $qq \leftrightarrow qq$, $gg \leftrightarrow q\bar{q}$, $q\bar{q} \leftrightarrow gg$, etc.).
• Cross-Section Regularization:
  • They start with leading-order pQCD matrix elements.
  • Screening: To regulate infrared divergences, they screen the propagators using scaled thermal masses ($\mu_D = \sqrt{\kappa}m_D$ and $\mu_F = \sqrt{\kappa}m_F$) rather than the momentum-dependent self-energies used in the AMY framework.
  • Correction of Unphysical Values: The original matrix elements were found to yield negative cross-sections or divergences at low center-of-mass energy ($\sqrt{s} \to 0$). The authors introduced "minimal" and "extra" screenings (modifying terms like $1/t$ and $1/s$) to ensure all total and transport cross-sections are finite and non-negative for all collision energies. This is crucial for numerical stability in transport models.
• Parameter Tuning:
  • Screening Coefficient ($\kappa$): A multiplicative factor $\kappa$ is introduced to the thermal masses to compensate for the difference between using thermal masses (simpler approach) and self-energies (AMY approach).
  • Momentum Scale ($Q$): The strong coupling constant $\alpha_s(Q^2)$ is evaluated at a momentum scale $Q$. The study investigates the sensitivity of results to the choice of $Q$ (defaulting to $Q = \sqrt{\langle -t \rangle} = 3T$).

3. Key Contributions

• General Mapping: The paper provides a rigorous mapping between specific, finite-temperature parton cross-sections and the resulting shear viscosity using the CE method.
• Resolution of Divergences: They successfully regularized pQCD matrix elements to produce finite, non-negative cross-sections suitable for transport simulations, demonstrating that these modifications have negligible impact on the final viscosity calculation.
• Calibration of Screening: They determined that a screening coefficient of $\kappa = 0.4$ is required to make the CE results (using thermal masses) quantitatively match the leading-order AMY results (using self-energies).
• Explicit Fit Functions: The authors provide explicit fitting functions for $\eta g^4/T^3$ and $\eta/s$ as functions of $m_D/T$ and the number of quark flavors ($N_f$), facilitating easy implementation in kinetic models.

4. Key Results

• Agreement with AMY:
  • At $\kappa = 1$ (standard thermal mass screening), the CE results for $\eta g^4/T^3$ are qualitatively similar but quantitatively higher (by $\sim 50\%$) than AMY results.
  • At $\kappa = 0.4$, the CE results match the AMY results (for Boltzmann statistics and $2 \leftrightarrow 2$ processes) very well. This confirms that the discrepancy arises from the difference between thermal mass screening and self-energy screening.
• Shear Viscosity to Entropy Ratio ($\eta/s$):
  • The ratio $\eta/s$ is highly sensitive to the choice of the momentum scale $Q$.
  • Using the default $Q = 3T$, the authors find $\eta/s \approx 0.15$ at the QCD phase transition temperature ($T_c \approx 156$ MeV) for $N_f = 3$. This value is slightly below twice the conjectured AdS/CFT lower bound of $1/4\pi \approx 0.08$.
  • A smaller $Q$ leads to a larger coupling $g$, larger cross-sections, and consequently a lower $\eta/s$.
• Flavor Dependence ($N_f$):
  • $\eta/s$ initially increases when moving from pure gluon gas ($N_f=0$) to $N_f=1$, then decreases as $N_f$ increases further. This is due to the interplay between the dominance of the large $gg$ elastic cross-section and the increasing number of scattering channels and coupling strength as $N_f$ grows.
• Impact of Inelastic Scattering:
  • Including inelastic $2 \leftrightarrow 2$ processes ($gg \leftrightarrow q\bar{q}$, etc.) reduces the shear viscosity by 4% to 8% compared to elastic-only calculations, with the effect increasing with $N_f$.
• Impact of Extra Screenings:
  • The extra screenings added to fix negative/divergent cross-sections at low energies have a negligible effect (<1%) on the calculated shear viscosity. This is because the weighting functions in the CE integral (PDF$_\eta$) are suppressed at low center-of-mass energies ($v = \sqrt{s}/T$), meaning low-energy collisions contribute little to the total viscosity.

5. Significance

This work bridges the gap between microscopic parton transport models and macroscopic hydrodynamic properties.

1. Model Validation: It allows transport modelers to input specific cross-sections and directly calculate the resulting shear viscosity, enabling a direct comparison with experimental flow data to constrain the QGP equation of state and transport coefficients.
2. Theoretical Consistency: By calibrating the screening coefficient $\kappa$, the study reconciles the simpler thermal mass approach with the more rigorous AMY self-energy approach, providing a robust method for finite-temperature calculations.
3. Practical Utility: The provided fit functions for $\eta/s(T, N_f)$ serve as a ready-to-use tool for QCD effective kinetic theory and hydrodynamic simulations, particularly in the temperature range relevant to current heavy-ion collision experiments ($T \sim T_c$ to $4T_c$).

In summary, the paper establishes a robust, analytically derived framework for calculating QGP shear viscosity from first-principles cross-sections, resolving technical issues with divergences, and quantifying the sensitivity of viscosity to screening parameters and momentum scales.