Heavy quarkonium dynamics beyond the binding energy over temperature expansion
Abstract
We revisit the pNRQCD open-quantum-system description of bottomonium in a quark-gluon plasma, writing the full master equation in a spherical basis and establishing a framework for its numerical treatment. The equation does not rely on an expansion in the ratio of the binding energy over the temperature, while retaining the Born, Markov, and dilute approximations. It depends on suitable convolutions of Coulombic Green's functions with real-time chromoelectric correlators, rather than a finite set of transport coefficients. It contains jump-operator contributions that enter with both positive and negative signs, so it is not of a Lindblad type. We utilize optimal operator rotations to minimize the negative sector and define a positive-sector approximation that retains only the optimized positive terms, obtaining in this way a Lindblad equation that can be simulated efficiently with the open-source code QTraj for arbitrary chromoelectric spectral functions. We treat the discarded sector as a systematic uncertainty associated with positivity-breaking effects of the pseudo-Lindblad generator. Positivity breaking effects can be attributed to the approximations made in the derivation of the quarkonium master equation, such as the Born--Markov approximations. Using a spectral model constrained by known infrared and ultraviolet limits, we compute the in-medium width and Bjorken-evolved survival probabilities of the Υ(1S,2S,3S). The computations show the convergence of the binding-energy-over-temperature expansion at high temperatures and provide the first large-scale quantum-trajectory simulation results for bottomonium in the low temperature regime. We can thus follow the bottomonium evolution over the full temperature range relevant in heavy-ion collisions.
Related
- From the universal Lindblad equation to Boltzmann equations: in-QGP quarkonium dynamics
- Temperature evolution of spin dynamics in extended Kitaev models within an Abrikosov fermion theory beyond mean field
- High-temperature expansion of the inverse-square Ising chain in a magnetic field: exact coefficients and finite-size reconstruction
- Low-temperature thermal expansion of the α and ω phases of titanium
- First-principles transition-state tensorial cluster expansion of vacancy diffusion in Ta-W beyond the kinetically-resolved activation approximation
BibTeX
@misc{brambilla2026heavy,
title = {{Heavy quarkonium dynamics beyond the binding energy over temperature expansion}},
author = {Nora Brambilla and Arthur Yi-Ting Lin and Tom Magorsch and Antonio Vairo},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.00194},
eprint = {2611.00194},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.00194},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Heavy quarkonium dynamics beyond the binding energy over temperature expansion
AU - Nora Brambilla
AU - Arthur Yi-Ting Lin
AU - Tom Magorsch
AU - Antonio Vairo
PY - 2026
DA - 2026/11/01
PB - arXiv
DO - 10.48550/arxiv.2611.00194
UR - https://arxiv.org/abs/2611.00194
AB - We revisit the pNRQCD open-quantum-system description of bottomonium in a quark-gluon plasma, writing the full master equation in a spherical basis and establishing a framework for its numerical treatment. The equation does not rely on an expansion in the ratio of the binding energy over the temperature, while retaining the Born, Markov, and dilute approximations. It depends on suitable convolutions of Coulombic Green's functions with real-time chromoelectric correlators, rather than a finite set of transport coefficients. It contains jump-operator contributions that enter with both positive and negative signs, so it is not of a Lindblad type. We utilize optimal operator rotations to minimize the negative sector and define a positive-sector approximation that retains only the optimized positive terms, obtaining in this way a Lindblad equation that can be simulated efficiently with the open-source code QTraj for arbitrary chromoelectric spectral functions. We treat the discarded sector as a systematic uncertainty associated with positivity-breaking effects of the pseudo-Lindblad generator. Positivity breaking effects can be attributed to the approximations made in the derivation of the quarkonium master equation, such as the Born--Markov approximations. Using a spectral model constrained by known infrared and ultraviolet limits, we compute the in-medium width and Bjorken-evolved survival probabilities of the Υ(1S,2S,3S). The computations show the convergence of the binding-energy-over-temperature expansion at high temperatures and provide the first large-scale quantum-trajectory simulation results for bottomonium in the low temperature regime. We can thus follow the bottomonium evolution over the full temperature range relevant in heavy-ion collisions.
ER -