Coupled charm and charmonium transport in a strongly coupled quark-gluon plasma

This paper presents a self-consistent coupled transport framework for open and hidden charm in a strongly coupled quark-gluon plasma, utilizing thermodynamic $T$-matrix interactions constrained by lattice QCD to simultaneously describe charm-quark diffusion and charmonium kinetics, thereby successfully reproducing LHC Pb-Pb collision observables.

The Gist

Imagine you are trying to understand what happens inside a super-hot, super-dense soup made of the fundamental building blocks of the universe. This soup is called the Quark-Gluon Plasma (QGP). It's the state of matter that existed just microseconds after the Big Bang, and scientists recreate it today by smashing heavy lead atoms together at nearly the speed of light in giant machines like the Large Hadron Collider (LHC).

In this chaotic soup, there are heavy particles called charm quarks. Sometimes, a charm quark and an anti-charm quark pair up to form a "charm-onium" (like a tiny, heavy atom). The big question scientists are asking is: Do these pairs survive the heat, or do they get torn apart? And if they get torn apart, can they find each other again and reform?

This paper by Kaiyu Fu, Biaogang Wu, and Ralf Rapp presents a new, more accurate way to answer that question. Here is the story of their work, broken down into simple concepts.

1. The Old Problem: The "K-Factor" Crutch

Previously, scientists tried to model this soup using "perturbative" math. Think of this like trying to predict how a swimmer moves through water by only looking at the water when it's calm and still. But the QGP isn't calm; it's a violent, churning storm.

To make their old math work, scientists had to use a "fudge factor" (called a K-factor) to force their predictions to match reality. It was like saying, "Our theory says the swimmer should go this fast, but we know they go faster, so let's just multiply our answer by 1.5." It worked, but it wasn't a true understanding of the physics.

2. The New Approach: The "Universal Interaction"

The authors of this paper say, "Let's stop fudging it." They built a coupled transport framework.

Imagine the QGP as a crowded dance floor.

• Open Charm: A single charm quark is like a dancer trying to move through the crowd. They get bumped, pushed, and slowed down (diffusion).
• Charmonium: A charm-anticharm pair is like a couple holding hands, trying to dance together.

The breakthrough here is that the authors used the same set of rules for both the single dancer and the dancing couple. They didn't invent separate rules for each. They used a sophisticated tool called a T-matrix, which is like a super-accurate map of how particles bump into each other, based on data from "lattice QCD" (which is like a digital simulation of the universe's laws).

3. The "Ghost" Effect: Off-Shell Spectral Functions

This is the most creative part. In the old models, particles were treated like solid billiard balls. But in this hot soup, particles are fuzzy and blurry. They can exist in "off-shell" states—meaning they are temporarily borrowing energy or acting slightly differently than their standard definition.

The authors realized that these "blurry" states are crucial.

• The Analogy: Imagine trying to catch a fish in a river. If you only look for the fish when it's perfectly still, you might miss it. But if you realize the fish is splashing and moving erratically (broad spectral functions), you realize there are more places it could be.
• The Result: Because the particles are "blurry," it's actually easier for them to break apart (dissociate) because there are more ways for the collision to happen. However, it also means they can reform (regenerate) more easily because the "blur" allows them to find each other in more ways than a rigid model would predict.

4. The Dance of Survival and Rebirth

The paper tracks two competing processes over time:

1. Dissociation (Breaking Up): The heat of the soup is so intense it rips the charm-anticharm couples apart.
2. Regeneration (Reuniting): As the soup cools down, the scattered charm quarks wander around (diffusion) and sometimes bump into each other again, reforming the couple.

The authors' new model shows that because the interactions are so strong (the "strongly coupled" part), the couples break apart very quickly. But, because the quarks are constantly bumping into each other, they also reform very quickly.

5. The "Thermalization" Check

A key insight is about thermalization. When the soup is first created, the charm quarks are moving very fast and chaotically. They need time to slow down and "settle in" with the rest of the soup (reach thermal equilibrium).

• The Metaphor: Imagine a group of rowdy teenagers (charm quarks) entering a quiet library (the medium). At first, they are loud and moving fast. It takes time for them to calm down and start walking quietly.
• The Finding: The paper shows that if you assume the quarks calm down instantly, you get the wrong answer. You have to track their journey as they slowly calm down. The model shows that because the quarks haven't fully calmed down yet, the rate at which they reform couples is slightly lower than the "perfect equilibrium" rate. This is a crucial correction that makes the model much more realistic.

6. The Results: Matching the Real World

The authors tested their new model against real data from the LHC (collisions of lead atoms).

• The Outcome: Their model, which uses no "fudge factors" and treats the particles as fuzzy, interacting entities, successfully predicted how many charm-onium particles survive in different collision scenarios.
• The Pattern: In the most violent collisions (central collisions), almost all original couples are destroyed, but many new ones are born from the scattered pieces. In less violent collisions (peripheral), more of the original couples survive. The model captures this "U-shape" pattern perfectly.

Summary

In short, this paper is like upgrading a weather forecast.

• Old Model: Used simple rules and added a "correction factor" to guess the storm's path.
• New Model: Uses a super-complex, physics-based simulation that accounts for the chaotic, fuzzy nature of the storm itself.

By treating the heavy particles and the soup they swim in with a consistent, non-perturbative (no fudging) approach, the authors have given us a clearer, more honest picture of how matter behaves at the extreme temperatures of the early universe. They showed that even in the hottest, most chaotic environment imaginable, particles can break apart and come back together, governed by a delicate balance of forces that we are finally learning to calculate accurately.

Technical Summary

Here is a detailed technical summary of the paper "Coupled charm and charmonium transport in a strongly coupled quark-gluon plasma" by Fu, Wu, and Rapp.

1. Problem Statement

The production of charmonium (bound states of charm-anticharm quarks, such as $J/\psi$, $\chi_c$, and $\psi(2S)$) in ultra-relativistic heavy-ion collisions (URHICs) serves as a primary probe for the Quark-Gluon Plasma (QGP). The observed yields result from a complex interplay between suppression (dissociation due to color screening and scattering) and regeneration (recombination of $c\bar{c}$ pairs).

Current theoretical challenges include:

• Inconsistency in Transport Coefficients: Many models treat open charm (charm quarks) and hidden charm (charmonia) separately, often relying on perturbative QCD interactions with phenomenological $K$-factors to fit data. This fails to capture the strongly coupled nature of the QGP (sQGP), evidenced by small heavy-quark diffusion coefficients and non-perturbative features in lattice QCD.
• Off-Shell Effects: The QGP exhibits broad in-medium spectral functions. Standard perturbative approaches often neglect these off-shell effects, which significantly alter dissociation rates and phase space availability.
• Non-Equilibrium Dynamics: Regeneration depends on the momentum distribution of charm quarks, which are initially far from thermal equilibrium. Previous models often assumed instantaneous thermalization or used simplified relaxation-time approximations, failing to accurately capture the time-dependent evolution of charm quarks and their impact on regeneration.

2. Methodology

The authors develop a coupled charm-charmonium transport framework that self-consistently treats the diffusion of open charm and the kinetics of charmonium within a strongly coupled QGP.

A. Theoretical Foundation: Thermodynamic T-Matrix

The core of the framework relies on non-perturbative T-matrix interactions constrained by recent Wilson-line correlators (WLCs) computed in Lattice QCD.

• Unified Interaction: The same underlying heavy-light scattering amplitude is used to calculate transport coefficients for both open charm and charmonium.
• Off-Shell Spectral Functions: The model incorporates broad in-medium spectral functions for charm quarks and thermal partons. This allows for off-shell scattering, enhancing dissociation rates and modifying the available phase space.

B. Charm-Quark Transport (Open Heavy Flavor)

• Langevin Dynamics: The evolution of the charm-quark phase-space distribution is simulated using relativistic Langevin equations.
• Transport Coefficients: The friction (drag) and momentum diffusion coefficients are derived directly from the T-matrix formalism, including off-shell effects. This replaces the need for phenomenological $K$-factors.
• Time Evolution: The simulation tracks the gradual thermalization of charm quarks from their initial hard-scattering distributions to the medium's thermal equilibrium.

C. Charmonium Transport (Hidden Heavy Flavor)

• Boltzmann Equation: The time evolution of charmonium yields is governed by a Boltzmann equation with dissociation ($\alpha$) and regeneration ($\beta$) terms.
• Non-Perturbative Rates:
  • Dissociation: Calculated using a quasi-free approximation where a thermal parton scatters off a bound state, treating the charm quark as "half-off-shell."
  • Regeneration: Derived using the principle of detailed balance. Crucially, the authors implement an off-shell extension for the regeneration rate. They replace the thermal equilibrium distribution of charm quarks with the time-dependent, non-equilibrium distributions ($f_c(p,t)$) obtained from the Langevin simulations.
• Equilibrium Limit: The framework rigorously defines the equilibrium limit ($f^{eq}_\Psi$) as the ratio of regeneration to dissociation rates. This limit is constrained by lattice QCD data on charm-quark number susceptibility to determine the effective charm mass and fugacity ($\gamma_c$).

D. Medium Evolution

The transport equations are coupled to a schematic, isotropic, and isentropic fireball model representing the expansion of the medium in Pb-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV (LHC). The equation of state switches from Lattice QCD (partonic) to a Hadron Resonance Gas (hadronic) at $T_H \approx 170$ MeV.

3. Key Contributions

1. First Self-Consistent Coupling: This is the first work to use the same underlying heavy-light interaction and off-shell spectral functions to simultaneously evaluate charm-quark diffusion coefficients and charmonium dissociation/regeneration rates.
2. Off-Shell Regeneration: The authors construct a regeneration rate that explicitly accounts for the non-equilibrium momentum distributions of charm quarks and their off-shell spectral properties, ensuring the correct detailed balance and equilibrium limit are recovered only when charm quarks fully thermalize.
3. Lattice-Constrained Parameters: The transport coefficients and equilibrium limits are directly constrained by state-of-the-art Lattice QCD results (WLCs and susceptibility), removing the need for phenomenological tuning of reaction rates.
4. Canonical Ensemble Treatment: The model includes a sophisticated treatment of the canonical ensemble for charm number conservation, switching between canonical and grand-canonical ensembles based on the number of $c\bar{c}$ pairs and correlation volume.

4. Results

The framework was applied to Pb-Pb collisions at the LHC, comparing results with ALICE data for $J/\psi$, $\chi_c$, and $\psi(2S)$.

• Dissociation Rates: The non-perturbative rates are substantially larger than perturbative calculations, particularly for excited states ($\chi_c, \psi(2S)$). They exhibit a decreasing trend with momentum, contrasting with the increasing trend seen in leading-order perturbative calculations.
• Time Evolution:
  • In central collisions, primordial charmonia are almost entirely suppressed early in the QGP phase due to high dissociation rates.
  • Regenerated charmonia gradually approach their equilibrium limits. However, due to incomplete thermalization of charm quarks (especially in peripheral collisions), the actual yield remains below the statistical equilibrium limit.
  • The $\psi(2S)$ reaches equilibrium much faster than $J/\psi$ due to its larger dissociation rate.
• Centrality Dependence ($R_{AA}$):
  • The model reproduces the characteristic "dip" in $J/\psi$ $R_{AA}$ at intermediate centralities ($N_{part} \approx 100-150$) and the rise in central collisions due to regeneration.
  • It captures the stronger suppression of excited states ($\chi_c, \psi(2S)$) compared to $J/\psi$.
• Transverse Momentum ($p_T$) Spectra:
  • The model successfully describes the low-$p_T$ enhancement due to regeneration and the high-$p_T$ dominance of primordial production.
  • Improvement over previous work: By using off-equilibrium charm spectra instead of a blast-wave approximation, the model significantly reduces the previous overestimation of yields at low $p_T$ in peripheral collisions.
  • Discrepancies: The model slightly underestimates data in the intermediate $p_T$ region ($\sim 5$ GeV) for semi-central collisions, suggesting a need for space-momentum correlations (SMCs) and more realistic hydrodynamic evolution.

5. Significance

This work represents a major step forward in the quantitative understanding of heavy-flavor dynamics in the QGP. By unifying open and hidden charm transport under a single non-perturbative framework constrained by Lattice QCD, the authors demonstrate that:

• The strongly coupled nature of the QGP is essential for describing charmonium suppression and regeneration.
• Off-shell effects and non-equilibrium dynamics are critical for determining the final yields, particularly the regeneration component.
• The approach provides a robust, parameter-free (regarding reaction rates) baseline for interpreting experimental data, reducing the reliance on phenomenological $K$-factors.

The study highlights that while the current schematic fireball model captures the main trends, future improvements involving viscous hydrodynamics and space-momentum correlations are necessary to fully resolve discrepancies in the intermediate $p_T$ region. This framework sets a new standard for precision heavy-flavor phenomenology at the LHC and future colliders.