⚛️ phenomenology

How Quantum Is Bottomonium in the Quark-Gluon Plasma?

Using an open quantum system framework to simulate bottomonium evolution in the quark-gluon plasma, the study reveals that while the medium suppresses quantum nonclassicality, a finite residual quantum structure persists due to the preferential dissolution of weakly bound states, demonstrating that classicalization is partial and observable-dependent.

Original authors: Nora Brambilla, Tom Magorsch

Published 2026-10-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Nora Brambilla, Tom Magorsch

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 a heavy-ion collision, where protons and neutrons are smashed together with such violence that they melt into a primordial soup, a unique form of matter emerges. This is the quark-gluon plasma, a state of existence that filled the universe mere microseconds after the Big Bang. In this seething, super-hot fluid, the fundamental building blocks of matter—quarks and the gluons that bind them—are no longer confined inside individual particles but roam freely. To understand how this plasma behaves, physicists often study heavy particles called bottomonium, which are pairs of a heavy bottom quark and its antimatter twin, orbiting one another. As these pairs travel through the plasma, they interact with the surrounding heat and chaos, and scientists have long debated a fundamental question: does this interaction strip away their quantum nature, turning them into simple, classical objects that follow predictable paths like tiny billiard balls, or do they retain a ghostly, wave-like complexity that defies classical description?

For decades, the standard approach to modeling these interactions relied on classical physics, treating the particles as if they were subject to random kicks and friction, much like pollen grains jiggling in water. This method, known as a Langevin equation, is computationally simple and has been widely used to predict how many of these heavy particles survive a collision. However, a more rigorous, fully quantum mechanical framework has recently emerged, describing the system as an "open quantum system" where the particle is constantly exchanging information with its environment. This new framework suggests that the particle's evolution is governed by a complex set of rules that preserve its quantum identity. The critical question remained: is the classical approximation good enough, or does the quantum nature of the particle survive the journey through the plasma in a way that matters for what we actually observe?

In a recent study, researchers Nora Brambilla and Tom Magorsch set out to answer this by simulating the journey of bottomonium through a realistic model of the quark-gluon plasma. They did not rely on the simplified classical equations but instead used a sophisticated quantum master equation derived from the fundamental theory of strong interactions. They tracked the evolution of the particle's state over time, watching how its internal structure changed as it moved from the moment of creation until it cooled down. Their goal was to measure a specific property called "negativity," which serves as a strict test for quantumness. In the quantum world, a particle can exist in a superposition of states, and this leads to regions in its mathematical description where the probability of finding it becomes negative—a concept that is impossible in classical physics, where probabilities must always be zero or positive. If this negativity disappears, the particle has effectively become classical. If it remains, the particle is still behaving in a distinctly quantum way.

The simulation revealed a nuanced picture that challenges the idea of a complete transformation. As the bottomonium entered the hot plasma, the surrounding medium acted as a filter, rapidly suppressing the quantum negativity and the coherence between different positions. The researchers observed that the initial, highly complex quantum state quickly lost its wild fluctuations, and the negativity dropped sharply. However, contrary to the expectation that the particle would eventually become a purely classical object, the negativity did not vanish entirely. Instead, it settled at a stable, finite level, hovering around eighteen percent of its maximum possible value. This residual quantum structure was not random; it was closely tied to the specific state of the particle. The simulation showed that the plasma preferentially dissolved the loosely bound, larger versions of the particle, while the tightly bound, compact versions survived. This process, which the authors call "conditional purification," meant that the surviving ensemble of particles became increasingly dominated by the most stable, ground-state configuration.

The study further explored how this remaining quantum structure affects what scientists can actually measure in experiments. They examined the probability of finding the particle in specific excited states, known as the 2S and 3S states, compared to the ground state. They discovered that for the excited states, the calculation of survival probability relies heavily on the delicate cancellation between positive and negative contributions in the quantum description. In a classical world, these negative regions would not exist, and the calculation would fail to reproduce the correct result. This means that for the excited states, the quantum nature of the particle is not just a theoretical curiosity but a practical necessity for accurate predictions. The ground state, however, was found to be largely insensitive to these quantum fluctuations, meaning that classical approximations might still work well for the most stable particles but would likely fail for the more fragile, excited ones.

Ultimately, the work demonstrates that the transition from quantum to classical is not a simple on-off switch but a partial and selective process. The plasma does not erase the quantum nature of bottomonium; it reshapes it. By dissolving the weakly bound components, the medium forces the surviving particles into a state that retains a specific, finite amount of quantum complexity. This residual structure is enough to invalidate a purely classical description for certain observables, particularly those involving excited states. The findings suggest that while classical models can capture the broad behavior of the most stable particles, they miss the subtle, essential details required to understand the full spectrum of heavy quarkonium in the early universe. The study concludes that the quantum world does not simply fade away in the heat of a collision; instead, it endures in a refined form, leaving a measurable imprint on the particles that survive the inferno.

Technical Summary: How Quantum Is Bottomonium in the Quark-Gluon Plasma?

Problem Statement
The evolution of heavy quarkonium (specifically bottomonium) in the quark-gluon plasma (QGP) is frequently modeled using classical transport equations, such as rate, Boltzmann, or Langevin equations. However, recent advances combining non-relativistic effective field theories (EFTs) with the open quantum system framework have enabled full quantum mechanical descriptions of this non-equilibrium evolution. This duality raises a fundamental question: Is a quantum mechanical treatment strictly necessary, or can the evolution be accurately approximated by classical transport equations? Specifically, the transition to a classical description (classicalization) relies on the assumption that the Wigner transform of the system's density matrix becomes non-negative, allowing it to be interpreted as a classical probability distribution. This paper investigates the extent to which bottomonium retains quantum structure (nonclassicality) during its evolution in the QGP and how this structure impacts physical observables.

Methodology
The authors utilize potential non-relativistic QCD (pNRQCD) to describe bottomonium in a strongly coupled medium. Within the "quantum Brownian regime" (where the binding energy $E$ is smaller than the thermal scale $\pi T$ ), they employ a Lindblad equation derived from first principles to model the real-time evolution of the bottomonium density matrix $\rho(t)$ .

Key methodological steps include:

Simulation Setup: Using the open-source code QTraj, the authors simulate the evolution of a bottomonium state in a Bjorken-expanding medium, cooling from $T_0 = 425$ MeV to 190 MeV. The simulation employs leading-order Coulombic singlet and octet potentials and Monte Carlo wave-function methods to reconstruct the full density matrix.
Conditional Density Matrix: To isolate the quantum structure of the surviving ensemble from the depletion of the system, the authors define a conditional density matrix $\rho_{c|\ell_{max}}$ restricted to the color-singlet, low-angular-momentum ( $\ell=0$ ) sector.
Quantifying Nonclassicality:
- Wigner Negativity: They compute the Wigner transform of the conditional density matrix and define a negativity metric $\eta$ as the fraction of the phase-space volume where the distribution is negative.
- Position-Space Coherence: They calculate the relative entropy of coherence $C_r$ in the position basis to measure off-diagonal elements in the density matrix.
Observable Sensitivity Analysis: The authors decompose the overlap integrals (survival probabilities) $\text{Tr}[\rho \Pi_{nS}]$ into contributions from positive and negative regions of the Wigner distributions of both the evolving state and the bound-state projectors ( $\Pi_{nS}$ ). This allows them to assess how sensitive specific observables (1S, 2S, 3S survival) are to the negative regions of the Wigner function.

Key Results

Partial Classicalization: The medium suppresses both Wigner negativity and position-space coherence. However, this suppression is incomplete; both quantities decay rapidly initially but plateau at finite, non-zero values.
Residual Quantum Structure: The late-time plateau of the negativity ( $\approx 18\%$ ) and coherence is close to the intrinsic values of the 1S bound state itself. This indicates that the surviving ensemble does not become fully classical but rather converges toward the quantum structure of the ground state.
Conditional Purification: The finite residual quantum structure is attributed to "conditional purification." The medium-induced dipole interaction preferentially dissolves weakly bound and unbound modes (which have large mean squared radii $\langle r^2 \rangle$ ). Consequently, the surviving singlet S-wave ensemble becomes increasingly dominated by the tightly bound 1S state. Since the 1S state inherently possesses Wigner negativity and coherence, the ensemble retains these features.
Observable Dependence:
- 1S State: The survival probability of the 1S state is largely insensitive to the negative regions of the Wigner distribution.
- Excited States (2S, 3S): The overlaps for excited states are highly sensitive to negative regions. The physical yield for these states arises from a near-cancellation between large positive contributions ( $O_{++}$ ) and large negative contributions ( $O_{+-}$ ) in the phase-space integral.
- Implication for Classical Approximations: Classical transport schemes that assume a non-negative Wigner distribution (or use positive acceptance criteria for bound states) inherently miss the $O_{+-}$ cancellation term. While this may not significantly affect 1S predictions, it can lead to $O(1)$ relative errors in the yields of excited states (2S, 3S).

Significance and Claims
The paper claims that the classicalization of bottomonium in the QGP is partial and observable-dependent.

It challenges the universal applicability of classical Langevin-type approximations for all quarkonium observables. While the 1S survival probability may be robustly described by classical methods, the yields of excited states are fundamentally tied to quantum interference effects arising from negative Wigner regions.
The study demonstrates that the "classicality" of a system is not an intrinsic property of the state alone but depends on the interplay between the state, the operations (dissipation), and the specific measurements (projectors) being performed.
The authors suggest that double ratios of nuclear modification factors, such as $R_{AA}[\Upsilon(3S)]/R_{AA}[\Upsilon(1S)]$ , could serve as experimental handles to probe this residual quantum structure, as they are sensitive to the differential impact of quantum effects on excited versus ground states.

The work does not claim that classical models are entirely invalid but highlights their limitations in describing excited-state yields where quantum interference is critical. It establishes that the "conditional purification" mechanism drives the system toward a state that retains the quantum characteristics of the 1S ground state, preventing a total transition to a classical probability distribution.

Technical Summary: How Quantum Is Bottomonium in the Quark-Gluon Plasma?

More like this