⚛️ phenomenology

Heavy quarkonium dynamics beyond the binding energy over temperature expansion

This paper presents a novel non-perturbative framework for simulating bottomonium dynamics in a quark-gluon plasma across the full temperature range by deriving a master equation that avoids the binding-energy-over-temperature expansion and utilizing an optimized positive-sector Lindblad approximation to enable efficient quantum-trajectory simulations.

Original authors: Nora Brambilla, Arthur Yi-Ting Lin, Tom Magorsch, Antonio Vairo

Published 2026-11-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Nora Brambilla, Arthur Yi-Ting Lin, Tom Magorsch, Antonio Vairo

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 extreme heat of a heavy-ion collision, where atomic nuclei smash together at nearly the speed of light, the ordinary rules of matter dissolve. Protons and neutrons melt into a seething, primordial soup known as the quark-gluon plasma. In this state, the fundamental particles of the universe, quarks and gluons, roam freely rather than being locked inside larger particles. To understand how this plasma behaves, physicists often look for specific "messenger" particles that survive the collision. Among the most important are bottomonium states, which are tightly bound pairs of a heavy bottom quark and its antimatter partner. These pairs act like delicate probes; as they travel through the hot plasma, the intense heat can pull them apart, causing them to "melt" or dissolve. By measuring how many of these pairs survive, scientists can infer the temperature and properties of the plasma, much like checking the depth of a river by seeing how far a stone sinks.

For years, researchers have used a specific mathematical framework to describe how these bottomonium pairs evolve inside the plasma. This framework treats the particle pair as a system interacting with a noisy environment, similar to how a boat is buffeted by waves. However, previous calculations relied on a simplifying assumption: that the plasma was so hot that the energy holding the quark pair together was negligible in comparison. This assumption worked well for the scorching early moments of a collision but broke down as the plasma cooled. As the fire of the collision faded, the plasma spent a significant amount of time at lower temperatures where that simplifying assumption was no longer valid. This left a gap in our understanding, preventing scientists from tracking the fate of these particles through the entire life of the collision, from the hottest start to the cooler end.

A team of physicists has now bridged this gap by developing a new, more complete way to calculate the behavior of bottomonium in the quark-gluon plasma. Instead of relying on the old assumption that the plasma is infinitely hot compared to the binding energy of the particles, they derived a master equation that works across the full temperature range. This new approach does not force the complex interactions of the plasma into a few simple numbers. Instead, it connects the particle's evolution directly to the detailed, real-time fluctuations of the electric fields within the plasma. The researchers translated this complex equation into a format that can be solved by powerful computers, allowing them to simulate the journey of bottomonium particles from the moment of the collision until the plasma cools down to about 155 million degrees.

The core of their work involved tackling a mathematical hurdle that arises when removing the "hot plasma" simplification. In their new equation, some terms appeared with negative signs, which created a problem: in quantum mechanics, probabilities must always be positive, and these negative terms threatened to produce impossible, negative probabilities. Rather than discarding the equation, the team found a clever way to rearrange the mathematical pieces. They rotated the components of the equation to minimize the size of these problematic negative terms, effectively isolating them into a small, manageable sector. They then focused their simulations on the dominant positive terms, which form a stable, reliable equation. By treating the small, discarded negative part as a known uncertainty, they created a robust tool that could be run on large-scale computer simulations without breaking the fundamental laws of physics.

Using this new framework, the researchers simulated the survival of three specific bottomonium states, known as the 1S, 2S, and 3S, as they traveled through a cooling plasma. They found that at high temperatures, their results matched well with previous, simpler models, confirming that the old methods were correct in that regime. However, as the temperature dropped below 190 MeV, the new simulations revealed a different picture. The rate at which the particles dissolved slowed down significantly compared to the predictions of the older models. This suggests that at lower temperatures, the bottomonium pairs are more resilient than previously thought. The team also observed that for the most tightly bound state, the 1S, the particles were less likely to recombine or reform after being broken apart, whereas the looser, excited states showed a greater tendency to regenerate.

The study provides the first large-scale quantum simulations of bottomonium dynamics in this lower temperature regime, a region that was previously inaccessible to this type of analysis. The results show that the older, simplified models tend to overestimate how quickly these particles are destroyed as the plasma cools. By offering a more accurate description of the entire temperature range, this work allows physicists to better interpret data from heavy-ion collision experiments. It suggests that the suppression of bottomonium, the measure of how many pairs are lost, is influenced by the cooling history of the plasma in a way that was previously missed. The researchers emphasize that while their method is a significant step forward, the precise behavior of the plasma's electric fields remains a complex input that requires further study. Nevertheless, this new framework opens the door to a more complete understanding of how matter behaves in the most extreme conditions in the universe, connecting the earliest moments of the Big Bang to the data collected in modern particle accelerators.

Technical Summary: Heavy Quarkonium Dynamics Beyond the Binding Energy over Temperature Expansion

Problem Statement
The suppression of heavy quarkonium (specifically bottomonium) in the Quark-Gluon Plasma (QGP) is a key signature for probing the properties of hot QCD matter. While the open quantum system framework, utilizing Potential Non-Relativistic QCD (pNRQCD), has successfully described bottomonium evolution at high temperatures, existing approaches rely on an expansion in the ratio of the binding energy to the temperature ( $E/(\pi T)$ ). This expansion, which leads to Lindblad-type master equations dependent on a finite set of transport coefficients, is valid only when $\pi T \gg E$ . Consequently, these methods fail to accurately describe bottomonium dynamics at lower temperatures, a regime where the QGP spends a significant amount of time during its cooling evolution in heavy-ion collisions. Previous attempts to solve the unexpanded Redfield equation (which includes non-Lindblad terms) have been limited by the neglect of quantum jumps or the use of small lattices.

Methodology
The authors revisit the pNRQCD open-quantum-system description of bottomonium without expanding in $E/(\pi T)$ . The methodology involves several key theoretical and numerical steps:

Derivation of the All-Orders (AO) Master Equation: Starting from the pNRQCD master equation (a Redfield equation), the authors project the evolution of the density matrix onto a spherical basis. This results in a "pseudo-Lindblad" equation containing 12 jump operators. Unlike standard Lindblad equations, six of these operators enter with negative signs (associated with negative eigenvalues of the Kossakowski matrix), meaning the equation does not strictly preserve the positivity of the density matrix.
Handling Positivity Breaking: Recognizing that simulating the full non-Lindblad equation with quantum jumps is computationally prohibitive due to the need for specialized unraveling schemes (e.g., doubled Hilbert spaces), the authors employ an "optimal operator rotation" technique. They rotate the jump operator basis to minimize the norm of the negative contributions.
Positive-Sector Approximation (AO+): Based on the minimization, the authors define a "positive-sector approximation" where the optimized positive terms are retained, and the residual negative sector is discarded. This yields a genuine Lindblad equation (the AO+ equation) that can be efficiently simulated using the standard Monte Carlo wave function algorithm implemented in the open-source code QTraj. The discarded negative sector is treated as a systematic uncertainty associated with the approximations (Born-Markov) used in the derivation.
Chromoelectric Correlator Modeling: Since the AO equation depends on convolutions of Coulombic Green's functions with real-time chromoelectric correlators rather than simple transport coefficients, the authors construct a physical model for the spectral function of the adjoint chromoelectric correlator. This model interpolates between known infrared (constrained by the heavy-quark diffusion coefficient $\hat{\kappa}$ ) and ultraviolet (perturbative) limits.
Numerical Implementation: The authors perform large-scale simulations on a radial lattice ( $L=40 \text{ GeV}^{-1}$ , $N=1500$ points) for $\Upsilon(1S, 2S, 3S)$ states under a Bjorken temperature evolution ( $T_0 = 400$ MeV down to 155 MeV). They compare results from the full AO equation (via effective Hamiltonian evolution), the AO+ Lindblad equation (with and without quantum jumps), and the standard Leading Order (LO) and Next-to-Leading Order (NLO) expanded master equations.

Key Contributions

Framework Extension: The paper establishes a framework for the numerical treatment of the pNRQCD master equation that does not rely on the $E/(\pi T)$ expansion, thereby extending the applicability of open-quantum-system simulations to the low-temperature regime ( $T < 190$ MeV).
Optimal Operator Rotation: The authors demonstrate the application of optimal operator rotations to the bottomonium problem, successfully minimizing the negative sector to enable efficient Lindblad-based simulations (QTraj) for arbitrary chromoelectric spectral functions.
First Large-Scale Low-Temperature Simulations: This work provides the first large-scale quantum-trajectory simulation results for bottomonium in the low-temperature regime within the pNRQCD master-equation framework.
Systematic Uncertainty Quantification: The paper defines the impact of the discarded negative sector as a systematic uncertainty, quantifying the difference between the full pseudo-Lindblad evolution and the positive-sector approximation.

Results

Convergence at High Temperatures: The computations confirm the convergence of the binding-energy-over-temperature expansion at high temperatures, with the AO results aligning well with NLO expanded results.
Low-Temperature Behavior: At low temperatures, the in-medium width calculated with the AO equation decreases rapidly, contrasting with the LO and NLO widths which remain notably larger. This suggests that applying expanded equations outside their validity regime may lead to an overestimation of suppression.
Quantum Jumps and Regeneration: Simulations including quantum jumps show that while the $\Upsilon(1S)$ survival probability is well-described by effective Hamiltonian evolution (neglecting jumps), excited states ( $\Upsilon(2S, 3S)$ ) receive non-negligible corrections from quantum jumps due to regeneration effects.
Positivity Sector Hierarchy: The numerical analysis reveals a clear hierarchy (approximately one order of magnitude) between the optimized positive and negative sector norms, validating the efficacy of the positive-sector approximation even at low temperatures.
Survival Probabilities: In a Bjorken medium, the $\Upsilon(1S)$ survival probability shows slightly stronger suppression in the AO+ case compared to NLO at early times, consistent with width calculations. For excited states, the AO+ equation predicts less suppression at low temperatures compared to NLO.

Significance and Claims
The paper claims to provide the first large-scale quantum-trajectory simulation of bottomonium dynamics that covers the full temperature range relevant to heavy-ion collisions, including the previously inaccessible low-temperature regime below 190 MeV. By removing the restriction that medium dependence be encoded solely in a finite set of transport coefficients, the framework connects the real-time evolution of bottomonium directly to the full real-time chromoelectric correlator. The authors assert that their approach allows for the treatment of the entire temperature evolution of the QGP, offering a more accurate description of bottomonium suppression than expanded methods in the low-temperature limit. They position the difference between the full and positive-sector results as a quantifiable systematic uncertainty arising from the Born-Markov approximations, rather than a fundamental flaw in the method. The work is presented as a necessary step toward accurate phenomenological applications, with future work planned to include hydrodynamics and feeddown.

More like this