⚛️ phenomenology

Heavy quark thermalization using phonon-heavy quark EFT

This paper utilizes a phonon-heavy quark effective field theory to calculate the heavy quark drag coefficient, determines unknown coupling coefficients by comparing with AdS/CFT results, and subsequently solves the Boltzmann equation to study heavy quark thermalization in the QGP via soft phonon interactions.

Original authors: Andreas Kirchner, Berndt Mueller, Jyotirmoy Roy, Sebastian Waeber

Published 2026-11-03

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Andreas Kirchner, Berndt Mueller, Jyotirmoy Roy, Sebastian Waeber

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 heart of matter, under conditions of extreme heat and density, protons and neutrons melt away into a primordial soup known as the quark-gluon plasma. This state of existence, which filled the universe mere microseconds after the Big Bang, is recreated for fleeting moments in massive particle colliders on Earth. Within this seething environment, heavy quarks—particles with masses far greater than the surrounding light particles—act as unique probes. Because they are created in the very first instants of a collision and are too heavy to be easily disturbed, they survive the entire lifecycle of the plasma, carrying information about its evolution from the initial explosion to the final freeze-out. However, understanding exactly how these heavy particles move and slow down within the plasma has long been a puzzle. At low speeds, the interactions become so intense and the number of particles so high that standard mathematical tools fail, leaving physicists without a clear picture of how these heavy travelers lose energy and eventually settle into equilibrium with the hot medium.

To solve this, a team of researchers has developed a new theoretical framework that treats the hot plasma not as a chaotic collection of individual particles, but as a fluid that supports wave-like vibrations. Just as a solid crystal can transmit sound through vibrations of its atomic lattice, this fluid medium can support collective ripples, which the researchers describe as phonons. By constructing a simplified set of rules—an effective field theory—that governs how a heavy quark interacts with these ripples, the team has created a manageable way to calculate the forces at play. This approach allows them to bypass the overwhelming complexity of direct particle-to-particle collisions at low speeds and instead focus on the heavy quark's interaction with the smooth, collective motion of the medium.

The researchers began by mapping out the fundamental interactions between the heavy quark and the fluid's vibrations. They found that the most significant interactions occur when the quark scatters off these phonons, a process that can be visualized as the quark bumping into ripples in the fluid. A critical discovery in their calculations was that the heavy quark does not simply glide through the medium; it acquires a finite "width," a measure of its instability, because it is constantly absorbing and emitting these thermal ripples. This width is essential for the math to work, preventing the calculations from breaking down. The team observed that as the heavy quark moves faster, particularly if it exceeds the speed of sound within the fluid, it begins to emit phonons in a cone-shaped pattern, similar to a sonic boom. This emission causes the quark's instability to grow rapidly, signaling the limits of their current theory.

With these interactions defined, the team calculated the drag force acting on the heavy quark, which describes how quickly the particle slows down as it plows through the plasma. They found that the drag is stronger in hotter environments and at lower speeds, meaning the quark loses energy more efficiently when the medium is hotter and when it is moving more slowly. However, their initial calculations contained an unknown factor, a coefficient that determines the precise strength of the connection between the quark and the phonons. To pin down this missing number, the researchers compared their results with a different, well-established method of studying extreme physics known as the AdS/CFT correspondence. By matching the way their drag force scaled with temperature and velocity to the predictions from this alternative theory, they were able to determine that the unknown coefficient is approximately 2.2.

Armed with this specific value, the researchers simulated the journey of a heavy quark as it travels through the plasma. They started with a group of quarks all moving in the same direction, representing a beam of particles just after a collision. As time passed in their simulation, the quarks began to interact with the phonons of the medium. These interactions acted like a series of gentle nudges, gradually scattering the quarks in different directions. Over time, the initial tight beam of particles spread out, losing its preferred direction and becoming a uniform, spherical cloud. This process demonstrated the thermalization of the heavy quark, showing how it eventually comes to rest and shares the same temperature and random motion as the surrounding fluid. The study concludes that this phonon-based approach provides a robust and natural way to describe heavy quark thermalization, offering a clear path forward for understanding how matter behaves in the most extreme conditions in the universe.

Technical Summary: Heavy Quark Thermalization using Phonon-Heavy Quark EFT

Problem Statement
The description of heavy quarks within the Quark-Gluon Plasma (QGP), particularly at low momenta, presents significant challenges when approached from first principles. In this kinematic regime, the large coupling strength and high gluon occupation numbers render standard perturbative QCD descriptions insufficient. While heavy quarks serve as excellent probes for the full evolution of heavy ion collisions due to their early production, a rigorous theoretical framework for their interactions with the medium at low momentum transfers has remained elusive.

Methodology
To address this, the authors develop an Effective Field Theory (EFT) that describes heavy quark interactions not with individual perturbative gluons, but with the collective excitations of the medium, modeled as phonons. The methodology proceeds through the following steps:

EFT Construction: The authors construct the EFT based on fluid dynamics, utilizing comoving coordinates $\phi^I(t, \vec{x})$ to track fluid volume elements. The theory is built upon symmetries including translations, rotations, and volume-preserving diffeomorphisms. The fluid Lagrangian is defined by an arbitrary function $F(X)$ , where $X = J_\mu J^\mu$ , encoding the equation of state.
Interaction Terms: Small fluctuations around equilibrium are introduced as phonon fields $\pi^I$ . The interaction Lagrangian between the heavy quark field $\psi$ and the fluid is constructed by promoting parameters (such as the heavy quark mass) to functions of the dimensionless variable $X$ . This yields interaction terms involving the heavy quark current and the fluid current $J^\mu$ , as well as couplings to gluon fields.
Scale Determination: The UV scale of the theory, $\Lambda_h$ , is estimated kinematically rather than via analytic matching to non-perturbative QCD. Constraints from sound attenuation lengths in both perturbative and strong coupling limits suggest $\Lambda_h \approx 5T$ .
Scattering and Width Calculation: The authors calculate the leading-order $2 \to 2$ heavy quark-phonon scattering cross-section. A critical component of this calculation is the inclusion of a quasi-particle width for the heavy quark propagator. This width arises from the absorption of thermal phonons and, for supersonic quarks, the emission of phonons (Mach cone effect). This width regulates divergences in the scattering matrix elements.
Drag Coefficient and Matching: Using the calculated cross-section, the heavy quark drag coefficient $A$ is derived within the Langevin equation framework. Since the absolute magnitude depends on an undetermined Wilson coefficient, the authors match the EFT results to AdS/CFT calculations. Specifically, they align the scaling behavior of the drag force with respect to temperature ( $T$ ) and velocity ( $v$ ) in the low-momentum limit, where AdS/CFT predicts $F \propto T^2 v$ .
Thermalization Simulation: With the Wilson coefficient fixed, the authors solve the Boltzmann equation for the heavy quark distribution function $f(t, \vec{p})$ to simulate the thermalization process.

Key Results

Scattering Properties: The heavy quark-phonon scattering cross-section is found to be finite in the collinear limit (due to longer interaction times) but vanishes for head-on collisions. The cross-section exhibits a resonance peak linked to the Mach condition of the quark moving through the medium, scaling inversely with the speed of sound.
Quasi-Particle Width: The heavy quark width is non-zero at all velocities but grows exponentially at large velocities (supersonic regime) due to Mach cone emission, signaling a breakdown of the EFT's unitarity at these extremes.
Wilson Coefficient: By matching the EFT drag force scaling to the AdS/CFT prediction ( $T^2 v$ ), the authors determine the previously unknown Wilson coefficient for the quark-phonon coupling to be approximately $\tilde{M}'(1) \approx 2.2$ .
Thermalization Dynamics: Simulations of the Boltzmann equation demonstrate that an initially anisotropic heavy quark distribution (moving along the z-axis) evolves into an isotropic, fully thermalized distribution through interactions with soft phonons. The results confirm that heavy quarks lose more energy in hotter media and at lower momenta.

Significance and Claims
The paper claims to provide a systematic, power-counted effective field theory description of heavy quark thermalization in the QGP via phonon interactions, bridging the gap where perturbative QCD fails. The primary contribution is the derivation of the phonon-heavy quark interaction vertices and the subsequent calculation of the drag coefficient. By successfully matching the EFT scaling to AdS/CFT results, the work establishes a method to fix the coupling strength of the effective theory without relying on non-perturbative matching. Finally, the framework successfully models the transition of heavy quarks from a non-equilibrium state to thermal equilibrium, offering a perturbative description of a process traditionally dominated by strong coupling considerations.

The authors note that future work will aim to extend this framework by including viscosities and applying it to the study of gluon momentum broadening.

Technical Summary: Heavy Quark Thermalization using Phonon-Heavy Quark EFT

More like this