⚛️ nuclear theory

Charmonia at Finite Momentum and Spatial Correlators in Quark-Gluon Plasma

Using a thermodynamic T-matrix approach with a separable potential, state-of-the-art heavy-quark selfenergies, and particle-hole excitations, this study successfully reproduces lattice QCD results for charmonium spatial correlators at finite momentum, revealing that zero modes are essential for explaining the increase in Euclidean correlators but not for the suppression of spatial correlators.

Original authors: Thomas Hardin, Ralf Rapp

Published 2026-08-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Thomas Hardin, Ralf Rapp

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

Deep within the heart of matter, where protons and neutrons usually hold tight together, there exists a state of extreme energy known as the quark-gluon plasma. In this fiery soup, which filled the universe just microseconds after the Big Bang and is recreated today in massive particle colliders, the fundamental building blocks of matter called quarks are no longer bound inside particles. Instead, they roam freely, interacting with a sea of other particles and gluons, the carriers of the strong nuclear force. To understand how this plasma behaves, physicists look for specific "probes" that can survive long enough to tell a story about the environment they are traversing. One of the most important of these probes is the charmonium, a particle made of a heavy charm quark and its antimatter partner, the charm antiquark. By studying how these pairs form, break apart, or change their properties while moving through the plasma, scientists can map out the invisible forces at play in this extreme state of matter.

The challenge lies in connecting the theoretical models of these forces with the actual data collected by supercomputers. These computers, running simulations called lattice quantum chromodynamics, can calculate how charmonium particles behave in a static, resting state. However, in a real collision, these particles are often moving at tremendous speeds, carrying significant momentum. When researchers try to translate the static computer data into a picture of moving particles, they run into a mathematical wall. The equations required to describe a moving particle involve a type of calculation that oscillates wildly, making it nearly impossible to get a clear, stable answer. This has left a gap in our understanding: we know what happens to these particles when they are still, but we have struggled to see how they behave when they are flying through the plasma.

A team of researchers at Texas A&M University has now bridged this gap by developing a new way to calculate the behavior of these moving particles. They started by creating a mathematical framework that respects the rules of relativity, ensuring that their description of the particle remains consistent regardless of how fast it is moving. They treated the interaction between the quark and antiquark not as a rigid force, but as a flexible connection that changes depending on the temperature and the speed of the particle. Crucially, they included a specific type of interaction that had often been overlooked in previous attempts: a process where the heavy quark interacts with the surrounding sea of lighter particles in the plasma. This interaction, which they refer to as a "zero mode," acts like a subtle background hum that influences the particle's journey.

Using this refined approach, the team simulated how the charmonium particles would appear to a detector looking at them from different angles and speeds. They compared their results with the existing data from the supercomputer simulations. The findings revealed that the behavior of these moving particles is indeed different from their stationary counterparts. As the particles move faster, their internal structure shifts, and the energy required to hold them together changes. The researchers found that the "zero mode" interaction plays a vital role in explaining why the data from the supercomputers shows a specific increase in signal at certain speeds. Without including this subtle interaction, the theoretical model would fail to match the computer-generated data.

The study also looked at how these particles are distributed in space as they move through the plasma. This is a particularly difficult calculation because the mathematical tools used to describe the spread of the particle become extremely unstable when the particle is moving fast. The team had to develop a careful numerical strategy to handle these instabilities, effectively cutting off the calculation at a point where the results would no longer change significantly. When they finally mapped out the spatial distribution, they found that the presence of the plasma causes a noticeable suppression, or reduction, in the likelihood of finding the particle at certain distances. This suppression becomes more pronounced as the temperature of the plasma rises.

The results offer a clear, qualitative match with the data from the supercomputer simulations, confirming that the new method captures the essential physics of the situation. The researchers found that while the "zero mode" is essential for explaining the behavior of the particles at high speeds, it is the general screening effect of the hot plasma that drives the suppression seen in the spatial distribution. The study does not claim to have solved every mystery of the quark-gluon plasma, but it provides a robust and reliable tool for interpreting future experiments. By successfully modeling how these heavy particles behave while in motion, the work opens the door to a deeper understanding of the strong force under extreme conditions, helping physicists decode the history of the early universe and the nature of matter itself.

Technical Summary: Charmonia at Finite Momentum and Spatial Correlators in Quark-Gluon Plasma

Problem Statement
Spatial correlation functions of quarkonia in the Quark-Gluon Plasma (QGP) can be calculated with high precision using finite-temperature lattice QCD (lQCD). However, interpreting these results in terms of microscopic spectral functions has proven difficult. Two primary challenges hinder this interpretation:

Momentum Dependence: Reliable calculation of the dependence on the bound state's total 3-momentum ( $P$ ) is non-trivial in potential models, which are not inherently Lorentz-covariant in the vacuum. This often leads to spurious momentum dependencies that obscure genuine medium effects.
Fourier Transform Oscillations: The Fourier transform required to move from momentum space to spatial correlators involves a highly oscillating integral for large $P$ , leading to slow convergence and significant numerical instability.

While recent progress has been made for the light-flavor sector using schematic ansätze, a microscopic model calculation for quarkonium spatial correlators that addresses both challenges simultaneously has been lacking.

Methodology
The authors employ a thermodynamic T-matrix approach to construct a Lorentz-covariant scattering framework. The methodology proceeds as follows:

Vacuum Formulation: The study begins with the covariant four-dimensional Bethe-Salpeter (BS) equation for the two-body scattering amplitude. To maintain Lorentz covariance while simplifying the interaction, a separable potential approximation is utilized ( $V(p', p) = C v(p) v(p')$ ). This potential is regularized with a momentum cutoff $\Lambda$ .
Dimensional Reduction: The BS equation is reduced to three dimensions using the Blankenbecler-Sugar scheme. This preserves the correct cut structure of the two-body propagator and ensures the T-matrix depends only on the invariant mass squared ( $s = E^2 - \vec{P}^2$ ), thereby avoiding spurious momentum dependencies.
Medium Modifications: The system is embedded in the QGP by incorporating two principal medium effects:
1. Heavy-Quark Selfenergies: The heavy quark and antiquark propagators are dressed with selfenergies ( $\Sigma$ ) derived from a "strongly coupled scenario" (SCS) constrained by lQCD free energies and correlators. These selfenergies include large imaginary parts (collisional widths $\sim 0.5$ GeV) characteristic of the sQGP.
2. Potential Screening: The interaction potential is screened via a Debye-like modification to the vertex function and a temperature-dependent coupling constant $C(T)$ . These parameters are constrained by matching the model's Euclidean correlator ratios (ECRs) at zero momentum to lQCD data.
Spectral Function Construction: The in-medium spectral functions are calculated from the T-matrix. Crucially, the authors include two distinct contributions:
1. Unitarity Cut: The standard two-particle ( $Q\bar{Q}$ ) continuum and bound-state contributions.
2. Zero Mode (Particle-Hole Excitation): A contribution arising from heavy-light scattering (particle-hole excitations), which is encoded in the heavy-quark selfenergies. This term is non-zero only for finite total momentum $\vec{P}$ .
Correlator Evaluation:
- Euclidean Correlators: Calculated via a Laplace transform of the spectral functions.
- Spatial Correlators: Calculated via a Fourier transform. To overcome the numerical challenge of the highly oscillating integrand, the authors implement a momentum cutoff ( $P_{cut} = 250$ GeV) above which the spectral function is approximated as momentum-independent, a strategy verified to yield stable results.

Key Results

T-Matrix Behavior: As the total momentum $P$ increases, the bound-state peaks in the T-matrix shift toward higher invariant masses (indicating a loss of binding energy due to enhanced screening) and exhibit a narrowing at high momenta ( $P \gtrsim 10$ GeV) for lower temperatures, while narrowing occurs at lower momenta for higher temperatures.
Euclidean Correlator Ratios (ECRs): The model successfully reproduces the lQCD observation that Euclidean correlators at large Euclidean time ( $\tau$ ) are significantly enhanced as momentum increases. The authors identify the zero-mode contribution as the primary driver of this enhancement (accounting for a 10–20% increase), with the unitarity cut playing a secondary role.
Spatial Correlators (SCRs): The calculated spatial correlators show a substantial suppression with increasing temperature.
- The unitarity cut drives the primary suppression, showing a mild temperature dependence that becomes more pronounced at large spatial separations ( $z$ ).
- The zero mode exhibits a strong, opposite temperature dependence (producing a negative contribution that enhances the correlator at low $z$ relative to the unitarity cut's suppression) but remains subleading in absolute magnitude compared to the unitarity cut.
- The model qualitatively reproduces the suppression trends found in unquenched lQCD data, though it overpredicts the suppression at the lowest temperature ( $T=194$ MeV).

Significance and Claims
The paper claims to present the first microscopic model calculation for the spatial correlators of quarkonia that simultaneously addresses the challenges of Lorentz-covariant momentum dependence and the numerical difficulties of the Fourier transform.

The authors emphasize that their approach provides a qualitative (and in some cases semi-quantitative) reproduction of lQCD findings:

The zero mode is identified as an essential mechanism for the momentum-dependent enhancement of Euclidean correlators, a feature previously difficult to capture in potential models.
The suppression of spatial correlators is largely driven by the in-medium bound-state properties (mass shifts and broadening) and the unitarity cut, with the zero mode playing a secondary, albeit distinct, role characterized by a negative sign and opposite temperature dependence.
The work bridges the gap between lQCD data and transport calculations by providing a framework that yields spectral functions with explicit momentum dependence, a necessary input for kinetic theory in heavy-ion collisions.

The authors remain modest regarding quantitative precision, noting that the exact origin of the spatial correlator suppression requires further investigation, potentially involving more realistic interaction potentials or updated heavy-quark selfenergies. They suggest that the observed uniform upward mass shift of the bound state with momentum is a reasonable candidate for the observed behavior.

Technical Summary: Charmonia at Finite Momentum and Spatial Correlators in Quark-Gluon Plasma

More like this