Heavy-quark transport across the QCD crossover driven by a lattice-constrained in-medium potential

This paper presents a self-consistent framework for heavy-quark transport across the QCD crossover that unifies perturbative and non-perturbative interactions via a lattice-constrained in-medium potential, successfully explaining the medium's extreme opacity near the critical temperature and achieving quantitative agreement with lattice QCD extractions for the spatial diffusion coefficient.

The Gist

Imagine you are trying to understand how a heavy boulder moves through a thick, swirling river. In the world of particle physics, that "boulder" is a heavy quark (like a charm or bottom quark), and the "river" is the Quark-Gluon Plasma (QGP)—a super-hot, super-dense soup of energy created just after the Big Bang or in giant particle colliders.

For decades, scientists had a problem trying to predict how these boulders move through the river. Their old maps (theories) worked well when the river was flowing fast and smoothly (high energy), but they fell apart when the water got turbulent and sticky near the "phase transition" (the point where the soup turns from a solid-like state to a fluid).

Here is a simple breakdown of what this new paper does to fix those maps.

1. The Old Problem: The "Soft vs. Hard" Split

Previously, scientists treated the river in two separate ways:

• The "Hard" part: Fast, direct hits between particles. They used standard math (perturbation theory) for this.
• The "Soft" part: Slow, gentle nudges. They used a different math for this.

To make the math work, they had to draw an arbitrary line in the sand, saying, "Everything faster than this line is Hard; everything slower is Soft." The problem? That line was made up. It didn't exist in nature. When the river got really sticky (near the critical temperature), this split caused the math to break down, predicting that the heavy boulders would glide through the soup too easily, when in reality, the soup was much stickier.

2. The New Solution: A Unified "Super-Map"

The authors of this paper say, "Let's stop drawing lines in the sand." Instead, they built a unified framework that treats the river as one continuous, complex flow.

They did this by creating a new kind of "interaction map" based on Lattice QCD.

• The Analogy: Imagine you want to know how a car drives on a road. Instead of guessing the road conditions, you look at a high-definition satellite map (Lattice data) that shows exactly where the potholes, speed bumps, and smooth asphalt are.
• The Innovation: They took this "satellite map" of the forces inside the plasma and turned it into a single, smooth mathematical formula. This formula works for both fast hits and slow nudges without needing to switch between different rules.

3. The Secret Ingredient: The "String"

The most exciting part of their discovery is what happens near the "critical temperature" (the point where the soup is just about to boil over).

• The Old View: Scientists thought the forces between particles were like magnets that get weaker the further apart they are (Yukawa potential).
• The New View: The authors found that near the critical temperature, there is a hidden force acting like a rubber band or a string.
  • Metaphor: Imagine the heavy quark is a swimmer. In the old model, the water just pushed back gently. In the new model, the water is full of invisible rubber bands attached to the swimmer. As the swimmer tries to move, these bands stretch and pull back hard.
  • This "string tension" is a non-perturbative effect (it's too complex for simple math) that makes the medium incredibly "opaque" or sticky. It explains why heavy quarks get stuck and slow down much more than anyone expected.

4. What They Found

By using this new "Rubber Band + Magnet" map, they calculated how fast these heavy quarks diffuse (spread out) in the plasma.

• The Result: Their calculations matched perfectly with the latest, most precise computer simulations (Lattice QCD).
• The Surprise: They found that near the critical temperature, the "rubber bands" (strings) are so strong that they dominate the movement. The heavy quarks are essentially dragging a heavy anchor through the soup.
• The Transition: As the temperature gets even higher (the soup gets hotter and thinner), the rubber bands snap (the strings melt), and the quarks start moving more like they would in a normal fluid, governed by the old "magnet" rules.

5. Why This Matters

This paper is a bridge. It connects two worlds that scientists have struggled to unite:

1. Static Snapshots: Computer simulations that show what the forces look like at a frozen moment.
2. Real-Time Motion: How particles actually move and collide in real-time.

The Takeaway:
Think of this paper as upgrading the GPS for heavy quarks. The old GPS told you to take a shortcut that didn't exist, making you think you'd arrive quickly. The new GPS, using the "rubber band" insight, tells you, "No, the road is actually a sticky, rubber-band-filled swamp right here."

This helps us understand the early universe better. It tells us that right after the Big Bang, the universe wasn't just a hot gas; it was a complex, sticky fluid with hidden "strings" holding everything together, making it much harder for heavy particles to move through it than we previously thought.

Technical Summary

1. Problem Statement

Heavy quarks (charm and bottom) serve as crucial probes for understanding the Quark-Gluon Plasma (QGP) created in ultra-relativistic heavy-ion collisions. However, describing their transport properties across the QCD crossover region (near the critical temperature $T_c$) presents significant theoretical challenges:

• Limitations of Perturbative QCD (pQCD): Traditional approaches rely on a "soft-hard" factorization scheme, separating momentum transfers by an arbitrary scale $\sqrt{-t^*}$. Purely perturbative methods severely underestimate interaction strengths near $T_c$ because they fail to capture strong non-perturbative chromoelectric correlations and confinement effects observed in Lattice QCD (LQCD).
• Disconnect between Static and Dynamic Observables: Existing phenomenological models often struggle to systematically map static LQCD potentials (derived from equilibrium simulations) to real-time scattering amplitudes required for transport calculations without invoking arbitrary cutoffs or ad-hoc parameters.
• Inconsistency in Diffusion Coefficients: Purely perturbative calculations yield spatial diffusion coefficients ($2\pi T D_s$) that are significantly larger than those extracted from LQCD, indicating a missing mechanism for strong coupling near the phase transition.

2. Methodology

The authors propose a unified, self-consistent framework that bridges the gap between static lattice potentials and real-time transport dynamics without arbitrary momentum separation scales.

A. Lattice-Constrained In-Medium Potential

• Generalized Gauss-Law Framework: The authors construct an effective heavy-quark potential $V(r, T)$ that combines a short-range Yukawa (screened Coulomb) term and a long-range confining string term.
• Lattice Calibration: The baseline vacuum parameters (coupling $\tilde{\alpha}_s$, string tension $\sigma$, and shift $V_0$) are fixed by averaging lattice data. The temperature-dependent effective screening mass $M_D(T)$ is extracted by fitting lattice data for the static potential.
• Parametrization: $M_D(T)$ is modeled using a robust HTL-inspired functional form that includes next-to-leading-order corrections and non-perturbative magnetic sector contributions, constrained by LQCD data.

B. Unified Scattering Amplitudes

• Fourier Transform: The coordinate-space potential is analytically Fourier-transformed into momentum space to create an effective interaction kernel $\tilde{V}(q, T)$.
• Lorentz Structure Assignment:
  • The Yukawa (Coulombic) component is treated as a vector interaction, ensuring consistency with leading-order pQCD and correct ultraviolet (UV) behavior.
  • The String (Confining) component is treated as a scalar interaction, reflecting the nature of flux-tube formation and spin-independent confinement. This assignment ensures the interference term between the two vanishes due to chiral symmetry (helicity conservation vs. flipping).
• Continuous Description: The interaction kernel serves as an effective $t$-channel propagator. It smoothly transitions from non-perturbative dominance in the infrared (IR) to perturbative pQCD behavior in the UV, eliminating the need for an artificial soft-hard separation scale.

C. Transport Coefficients Calculation

• The framework calculates scattering amplitudes for heavy quarks interacting with light quarks ($Qq \to Qq$) and gluons ($Qg \to Qg$).
• Energy Loss & Diffusion: Using the unified matrix elements, the authors compute the energy loss per unit path length ($-dE/dz$), transverse/longitudinal momentum diffusion coefficients ($\kappa_{T/L}$), and the jet transport coefficient ($\hat{q}$).
• Spatial Diffusion: The spatial diffusion coefficient $D_s$ is derived via the fluctuation-dissipation theorem in the static limit ($v \to 0$).

3. Key Contributions

1. Elimination of Arbitrary Scales: The model removes the reliance on the unphysical soft-hard separation scale ($\sqrt{-t^*}$), providing a continuous description of interactions across the full momentum range.
2. Non-Perturbative String Tension: The work rigorously demonstrates that the long-range confining string tension is indispensable for capturing the extreme opacity of the medium near $T_c$.
3. Unified Vector-Scalar Framework: By assigning distinct Lorentz structures (vector for Coulomb, scalar for string) to the interaction components, the model preserves pQCD consistency at high energies while incorporating lattice-driven enhancements at low energies.
4. Mass Hierarchy Analysis: The framework successfully extends to both charm and bottom quarks, analyzing the transition from mass-independent strong coupling near $T_c$ to mass-dependent kinematics at higher temperatures.

4. Key Results

• Spatial Diffusion Coefficient ($2\pi T D_s$):
  • The model predicts $2\pi T D_s \approx 0.5 \sim 1.7$ near $T_c$.
  • This shows striking quantitative agreement with recent HotQCD (2025) lattice extractions ($2\pi T D_s \approx 1.0^{+0.3}_{-0.2}$ at $T \approx 1.09 T_c$).
  • Purely perturbative (Yukawa-only) models fail to reproduce this, yielding values significantly larger than lattice data.
• Energy Loss and Momentum Broadening:
  • The non-perturbative string contribution significantly enhances energy loss and momentum broadening ($\hat{q}$) in the low-temperature regime ($T \lesssim 1.7 T_c$).
  • As temperature increases, Debye screening melts the flux tubes, and the results smoothly converge to the perturbative baseline.
  • For high-energy quarks ($E \sim 100$ GeV), the results align with JETSCAPE dynamic extractions, though the static potential model slightly underestimates $\hat{q}$ at very high temperatures due to the omission of radiative (inelastic) processes.
• Flavor Dependence:
  • Near $T_c$, the ratio of spatial diffusion for bottom to charm quarks ($D_s^b / D_s^c$) is close to unity, indicating that strong non-perturbative interactions overwhelm kinematic mass differences.
  • As temperature rises, a mass dependence emerges, with the ratio decreasing, consistent with lattice trends.

5. Significance and Outlook

• Theoretical Validation: The work provides a robust dynamical interpretation of the strong heavy-quark coupling near the QCD crossover, validating that static lattice potentials can be directly mapped to real-time transport coefficients.
• Phenomenological Impact: The derived transport coefficients offer a precise baseline for interpreting experimental heavy-flavor observables (e.g., $R_{AA}$ and $v_2$) at RHIC and LHC.
• Future Directions: The authors acknowledge that while the static potential paradigm excellently describes the crossover region, it is limited to elastic scattering. Future work will integrate these coefficients into a Langevin transport framework with gluon radiation (LGR) to account for inelastic processes and chromomagnetic effects, particularly for high-temperature regimes.

In summary, this paper establishes a unified, lattice-constrained framework that resolves long-standing discrepancies between perturbative calculations and lattice data regarding heavy-quark transport, highlighting the critical role of non-perturbative confinement forces in the QGP near the phase transition.