⚛️ phenomenology

Hot QCD: transport coefficients in strong and weak coupling regimes

This paper evaluates the shear viscosity and heavy quark spatial diffusion coefficient in a Quasi-Particle Model using novel $N_f=3+1$ lattice QCD data, employing both Chapman-Enskog and Green-Kubo formalisms to compare strong and weak coupling regimes.

Original authors: Gabriele Parisi

Published 2026-09-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Gabriele Parisi

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the first moments after the Big Bang, the universe was not made of the solid atoms that build our world today. Instead, it was a seething, super-hot soup of the most fundamental building blocks of matter: quarks and gluons. Under normal conditions, these particles are forever locked together inside protons and neutrons by an incredibly strong force. But when temperatures soar high enough, that force weakens, and the quarks and gluons break free, forming a state of matter known as the quark-gluon plasma. Scientists recreate this primordial state for fleeting instants by smashing heavy atomic nuclei together at nearly the speed of light in massive particle accelerators. The goal is to understand how this plasma behaves, how it flows, and how it cools down to form the matter we see around us. To do this, researchers measure how the plasma resists flowing, a property called viscosity, and how it slows down heavy particles moving through it, a property called diffusion.

A recent study by Gabriele Parisi, a researcher working with advanced theoretical models, has taken a fresh look at these properties. The work focuses on a specific challenge: the quark-gluon plasma is so strongly interacting that standard mathematical tools often fail to describe it accurately. To solve this, the team used a "quasi-particle model." In this approach, the invisible, massless particles that usually make up the plasma are treated as if they have gained weight due to the intense heat and interactions around them. By adjusting the mass of these particles to match data from supercomputer simulations of the strong force, the researchers created a more realistic picture of the plasma's internal environment. They then calculated how this heavy, hot fluid moves and how it drags on the heaviest particles within it, specifically the charm and bottom quarks, which act as probes because they are too heavy to be created by the heat of the collision itself and instead survive the entire event.

The researchers employed two distinct methods to calculate the plasma's resistance to flow, or shear viscosity. The first was a mathematical expansion technique that approximates the behavior of the fluid by looking at small deviations from a calm state. The second was a computer simulation that tracked millions of individual particles as they bounced off one another, using a statistical approach to mimic the chaotic collisions happening inside the plasma. By comparing the results of these two methods, the team found that they agreed well with each other, giving them confidence in their numbers. They discovered that the plasma's resistance to flow changes significantly with temperature. At lower temperatures, just above the point where the plasma forms, the fluid is quite "thick" and resistant. As the temperature rises, the fluid becomes less resistant, behaving more like a gas where particles interact less frequently.

A key part of the study involved looking at how heavy quarks, specifically the charm and bottom types, move through this plasma. Because these particles are so massive, they do not instantly settle into the same speed as the surrounding fluid. Instead, they drift through it, slowed down by friction and jostled by random collisions. The researchers calculated a "diffusion coefficient," which measures how quickly these heavy particles spread out. They found that the bottom quark, being heavier, moves more sluggishly than the charm quark. Furthermore, the study showed that the presence of a fourth type of quark (the charm quark) in the background plasma changes the way the fluid behaves, slightly altering the drag experienced by the heavy probes.

The most significant finding of the work is a comparison between how easily the heavy particles diffuse and how much the fluid resists flowing. In the world of theoretical physics, there are two famous limits for this relationship. One limit, derived from a theory involving extra dimensions of space, suggests that in a perfectly strongly interacting fluid, this ratio should be very small. The other limit, based on simple gas theory, suggests the ratio should be much larger for a weakly interacting gas. The researchers found that their results do not fit neatly into either box. At lower temperatures, the ratio is indeed small, consistent with a strongly interacting system. However, as the temperature increases, the ratio grows larger, moving toward the behavior of a weakly interacting gas. This confirms that the plasma transitions from a strongly coupled state to a weaker one as it gets hotter.

Interestingly, the study also revealed that the relationship between the heavy particle's movement and the fluid's flow is more complex than some previous theories suggested. The researchers found that the simple idea that the distance a heavy particle travels is directly linked to the distance a light particle travels in the fluid is not accurate, especially for the heaviest particles. This means that the standard formulas used to estimate these properties might need adjustment. The work provides a more nuanced view of the quark-gluon plasma, showing that it is a dynamic system that changes its character depending on how hot it is. By refining the models used to describe this state of matter, the study helps physicists better interpret the data coming from particle colliders, bringing us closer to understanding the fundamental nature of the universe's earliest moments.

Technical Summary: Transport Coefficients in Strong and Weak Coupling Regimes within a Quasi-Particle Model

Problem Statement
The study of the Quark-Gluon Plasma (QGP) created in ultra-relativistic heavy-ion collisions requires a theoretical framework capable of describing both perturbative and non-perturbative regimes. While perturbative QCD (pQCD) is valid at very high energies due to asymptotic freedom, experimental evidence from RHIC and LHC suggests that the QGP near the critical temperature ( $T_c \sim 155$ MeV) is a strongly coupled medium where the coupling constant is too large for standard perturbative treatments. Furthermore, the hadronization mechanism and the transition to a deconfined state remain open questions in the Standard Model. A specific challenge lies in accurately calculating transport coefficients—specifically the shear viscosity to entropy density ratio ( $\eta/s$ ) and the spatial diffusion coefficient of heavy quarks ( $D_s$ )—which characterize the medium's response and the interaction of heavy probes (charm and bottom quarks) with the bulk. Existing approaches often rely on massless parton approximations or simplified scattering models that may not capture the non-perturbative dynamics encoded in lattice QCD (lQCD) thermodynamics.

Methodology
This work employs a Quasi-Particle Model (QPM) to describe the QGP. In this effective theory, the strong interactions of the medium are encoded in temperature-dependent effective masses for quarks and gluons, rather than in a running coupling constant alone. The coupling $g(T)$ is determined by fitting the QPM energy density to lQCD data for $N_f = 3$ and $N_f = 3 + 1$ (including a charm quark in the bulk).

The methodology is divided into three main components:

Scattering Matrices with Massive Partons:
The author calculates tree-level scattering matrices ( $\mathcal{M}$ ) for quark-quark, quark-gluon, and gluon-gluon interactions. Crucially, unlike standard pQCD, the gluon propagator is modified to include a non-zero mass ( $M_g$ ), and the polarization sum is adjusted accordingly. These matrices are used to derive differential cross-sections for the transport calculations.
Heavy Quark Dynamics (Fokker-Planck Approach):
The motion of heavy quarks (charm and bottom) through the medium is treated as a Brownian motion, justified by the small momentum transfer approximation. The relativistic Boltzmann transport equation is expanded to derive the Fokker-Planck equation. From this, the drag coefficient $A(p, T)$ and the spatial diffusion coefficient $D_s$ are extracted. The calculations are performed for both $N_f=3$ and $N_f=3+1$ bulk scenarios.
Shear Viscosity ( $\eta$ ) Evaluation:
Two distinct formalisms are employed to calculate $\eta$ and compare their validity:
- Chapman-Enskog (CE) Expansion: A dual-component (binary mixture of quarks and gluons) CE formalism is developed. This involves expanding the distribution function around local equilibrium and solving for the first-order correction. The author explicitly derives the formulas for a binary mixture, accounting for different masses and cross-sections, and compare them against simple averaging methods.
- Green-Kubo (GK) Formalism: A relativistic transport code is developed using the test-particle method and a stochastic implementation of the collision integral. The shear viscosity is derived from the time correlation function of the shear stress tensor, $\langle \pi_{xy}(t)\pi_{xy}(0) \rangle$ , integrated over time.

Key Contributions

Novel Application of lQCD Data: For the first time, the evaluation of transport coefficients within a QPM framework utilizes novel lQCD data for $N_f = 3 + 1$ , allowing for the study of charm quark diffusion within a bulk that includes dynamical charm.
Dual-Component CE Formalism: The thesis develops and implements a rigorous two-component Chapman-Enskog formalism for a quark-gluon mixture. This avoids the simplifications of treating the medium as a single effective fluid or using simple harmonic averages of single-component viscosities.
Stochastic Transport Code: A dedicated transport code was developed to solve the relativistic Boltzmann equation stochastically, enabling the direct calculation of $\eta$ via the Green-Kubo relation for massive partons with pQCD-like scattering matrices.
Systematic Comparison: The work provides a systematic comparison between the CE and GK approaches, validating the transport code against simpler cases (constant cross-sections, single-component gases) before applying it to the full physical case.

Results

Shear Viscosity ( $\eta/s$ ):
- The results from the Green-Kubo method and the two-component Chapman-Enskog approximation show good agreement.
- The study highlights that the anisotropy of the scattering matrices (derived from pQCD) significantly increases $\eta/s$ compared to calculations assuming isotropic cross-sections with the same total magnitude.
- The two-component CE formalism is shown to be superior to simple averaging techniques (such as harmonic means of single-component viscosities), which tend to underestimate the viscosity of the mixture.
- The calculated $\eta/s$ values are consistent with other lQCD and phenomenological (PHSD) results.
Spatial Diffusion Coefficient ( $D_s$ ):
- The dimensionless quantity $2\pi T D_s$ exhibits a minimum near the critical temperature $T_c$ and increases at higher temperatures.
- The diffusion of bottom quarks is significantly slower (lower $D_s$ ) than that of charm quarks, consistent with the inverse dependence on mass.
- The $N_f = 3 + 1$ bulk yields slightly higher $D_s$ values than the $N_f = 3$ case, attributed to a smaller effective coupling constant $g(T)$ in the former at high temperatures.
- The results lie between the quenched lQCD results and those assuming infinite heavy quark mass, aligning with current lattice expectations.
The Ratio $(2\pi T D_s) / (4\pi \eta/s)$ :
- This ratio is analyzed as a discriminator between strong and weak coupling regimes. Theoretical limits suggest a value of $\approx 1$ for strong coupling (AdS/CFT) and $\approx 2.5$ for weak coupling (perturbative kinetic theory).
- The QPM results show that this ratio is an increasing function of temperature, transitioning from values near the strong coupling limit at low $T$ toward the weak coupling limit at high $T$ .
- However, the study finds that for heavy quarks (especially bottom) at low temperatures, the ratio can fall below 1, violating the conjectured AdS/CFT lower bound.
- The analysis reveals that the perturbative relation $D_s \approx \lambda_{tr}$ (where $\lambda_{tr}$ is the momentum transfer mean free path) is a poor approximation, particularly for heavy quarks and at low temperatures, leading to deviations from the expected weak-coupling limit of 2.5.

Significance and Claims
The paper claims that the Quasi-Particle Model, when coupled with a rigorous two-component transport formalism, successfully describes the transition from a strongly coupled to a weakly coupled regime in Hot QCD. The agreement between the microscopic Chapman-Enskog expansion and the stochastic Green-Kubo simulations validates the use of the two-component CE formalism for binary mixtures.

The author modestly concludes that while the ratio $(2\pi T D_s) / (4\pi \eta/s)$ serves as a useful indicator of coupling strength, the theoretical limits derived from AdS/CFT and perturbative kinetic theory are not strict bounds in the realistic, non-conformal QPM framework. Specifically, the violation of the $D_s \approx \lambda_{tr}$ approximation suggests that the relationship between diffusion and mean free path is more complex than previously assumed in simple kinetic estimates. The work emphasizes that future experimental precision at RHIC and LHC, particularly in the low-momentum region, will be necessary to further constrain these transport coefficients and potentially extend the model to include momentum-dependent interactions.

More like this