Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2608.06754

Quantum simulation of bottomonium dynamics in the quark-gluon plasma via the Lindblad equation

Linyuan Wei, Jin Hu, Anping Huang, Yunpeng Liu, Baoyi Chen

Published Aug 7, 2026• 0 citations
arXivHTMLPDF

Abstract

Quantum computing provides a powerful framework for simulating real-time dynamics in open quantum systems, offering key advantages for modeling heavy-quarkonium transport in high-energy nuclear collisions. In this work, we perform quantum simulations of the isotropic next-to-leading-order Lindblad equation for bottomonium in the quark-gluon plasma using a reduced spherical coordinate representation. We discretize operators and wavefunctions, map the physical state onto qubits, and execute time evolution via parameterized quantum gate operations. By extracting the $Υ(1S)$ survival probability, we quantitatively isolate the color-octet contribution, demonstrating that its overall impact is small in the final production of the bottomonium ground state $Υ(1S)$ in the hot QCD medium at temperatures accessible at the Large Hadron Collider. Additionally, we have further optimized the quantum simulation algorithm for the Lindblad equation. The improved algorithm requires only a single ancillary qubit to realize the Lindblad evolution, thereby minimizing the circuit significantly.

DOI
10.48550/arxiv.2608.06754

Related papers

Simulation of Lindblad equations for quarkonium in the quark-gluon plasma

Takahiro Miura, Yukinao Akamatsu, Masayuki Asakawa +1 more

We study the properties of the Lindbladian quantum mechanical evolution of quarkonia with non-Abelian charges (color-singlet and octet) in the quark-gluon plasma. We confirm that heavy quark recoils in the Lindblad equation correctly thermalize quarkonium colorful states within statistical errors from the simulation method. We also demonstrate that the Lindblad equation in the dipole limit can provide an efficient alternative method, which is applicable to a finite time evolution before thermalization and dramatically reduces the numerical cost. Our findings will serve as a foundation for large-scale simulation of quarkonium dynamics in the relativistic heavy-ion collisions.

arXiv:2205.15551•Physical review. D/Physical review. D.•Oct 5, 2022•23 citations

How Quantum Is Bottomonium in the Quark-Gluon Plasma?

Nora Brambilla, Tom Magorsch

We study how much quantum structure bottomonium retains during its evolution in the quark-gluon plasma. Within the open quantum system framework, we simulate the Lindblad equation obtained from potential nonrelativistic QCD in the quantum Brownian regime. From the resulting real-time evolution of the bottomonium density matrix, we compute the Wigner transform and study its negativity as a measure of nonclassicality and a test of a key condition underlying classical Langevin-type approximations. We find that the medium suppresses the negativity, which nevertheless approaches a finite plateau close to that of the $1\mathrm{S}$ state. The overlaps with the excited $2\mathrm{S}$ and $3\mathrm{S}$ states arise from near-canceling positive and negative phase-space contributions, driven by negative regions in the Wigner transforms of the bound-state projectors, while the $1\mathrm{S}$ overlap is insensitive to such contributions. We further show that the position-space coherence is similarly suppressed and plateaus close to the 1S value. We trace the finite residual quantum structure to conditional purification: the medium preferentially dissolves weakly bound and unbound modes with large $\langle r^2\rangle$, driving the surviving singlet $\mathrm{S}$-wave ensemble toward the $1\mathrm{S}$ state. More broadly, these results reveal that classicalization is partial and observable dependent, a feature that is relevant across applications of open quantum system methods in fundamental physics.

hep-phCross-listedarXiv:2610.02313•Oct 5, 2026•0 citations

General Lindblad equation for quarkonium evolution in a quark-gluon plasma

Aoumeur Daddi Hammou, Jean-Paul Blaizot, Pol Bernard Gossiaux +1 more

Accurate modelling and understanding of quarkonium production in ultrarelativistic heavy ion collisions requires a formalism that preserves the quantum properties of microscopic $\bar{Q}Q$ systems while treating the interaction of such pairs with the quark-gluon plasma (QGP). The open quantum system approach has recently emerged as one of the most fruitful schemes to meet such requirements. However, the quantum master equations obtained so far in this context are derived assuming a strict ordering between the QGP temperature $T$ and the $\bar{Q}Q$ energy gaps ($ΔE$) of the quarkonia bound states. This limits their predictive power since, as the QGP expands and cools down, the system traverse all regimes between the quantum Brownian motion (QBM) regime for $T\gtrsim ΔE$ to the quantum optical (QO) regime for $ΔE\gtrsim T$. In this paper, we derive and present a more general non-abelian quantum master equation of the Lindblad type, which does not suffer from these limitations and thus allows to faithfully describe the quantum evolution of the $Q\bar{Q}$ pairs during the whole QGP-evolution. We also provide some illustration of the key quantities governing this equation.

nucl-thCross-listedarXiv:2608.15347•Aug 18, 2026

From the universal Lindblad equation to Boltzmann equations: in-QGP quarkonium dynamics

Aoumeur Daddi Hammou, Pol Bernard Gossiaux

Recently, a set of coupled singlet-octet universal Lindblad equations (ULEs) was derived within the framework of non-relativistic QCD (NRQCD) to describe quarkonium dynamics in the quark-gluon plasma (QGP). These equations provide a unified quantum description spanning the quantum Brownian and quantum optical regimes. In this work, we further develop this framework and establish its connection with semiclassical transport. We first derive the universal Lindblad equations within the potential non-relativistic QCD (pNRQCD) effective field theory and show that they coincide with the small-dipole limit of the NRQCD ULEs. We then derive the semiclassical limit of the NRQCD ULEs, obtaining a set of coupled singlet-octet Boltzmann equations. To our knowledge, this is the first derivation of Boltzmann transport equations directly from the universal Lindblad framework. The resulting equations are valid beyond the small-dipole approximation, allowing the evolution of heavy-quark pairs from compact to widely separated configurations. Taking their small-dipole limit allows a direct comparison with the Boltzmann equations of Yao et al. [Phys. Rev. D 99, 096028 (2019)], which were derived within pNRQCD from the Davies secular equation, relying on the rotating-wave approximation (RWA). While the singlet equations are found to be in near-complete agreement, the octet equation contains an additional collision term describing transitions within the continuum of octet scattering states, which is absent from the RWA-based derivation. Finally, we derive the leading quantum correction to the singlet Boltzmann equation. Our results establish a more general and systematic theoretical foundation for the semiclassical transport description of quarkonium in the QGP.

hep-phCross-listedarXiv:2608.17372•Aug 19, 2026

Thermalization of Bottomonium in the Quark-Gluon Plasma

Nora Brambilla A, Tom Magorsch, Antonio Vairo

We study the approach to equilibrium of bottomonium in the quark-gluon plasma within the open quantum system framework. We perform large-scale simulations of the long-time behavior in three dimensions using the quantum trajectory method to observe the emergence of steady states and determine the timescale of thermalization in position-, angular-momentum-, and color-space. We find that the thermalization timescale increases with decreasing temperature and decreasing coupling to the medium, which is given by transport coefficients of the medium. Additionally, we observe that the steady states exhibit small corrections to the Gibbs state due to medium interactions and show that these corrections diminish for weaker medium coupling and higher temperature. At a temperature of $450\,$MeV, quarkonium relaxes to a state that is approximately thermal, with the most significant correction being a smaller overlap of the $1S$ state with respect to the Gibbs state. We compare these findings with the master equation obtained at leading order in the expansion of the binding energy over the temperature, which we find to have a trivial steady state.

arXiv:2508.11743•Aug 15, 2025•0 citations

Differentiable quantum-trajectory simulation of Lindblad dynamics for QGP transport-coefficient inference

Lukas Heinrich, Tom Magorsch

We study parameter estimation for the transport coefficients of the quark-gluon plasma by differentiating open-quantum-system-based Monte Carlo simulations of quarkonium suppression. The underlying simulator requires solving a Lindblad equation in a large Hilbert space, which makes parameter estimation computationally expensive. We approach the problem using gradient-based optimization. Specifically, we apply the score-function gradient estimator to differentiate through discrete jump sampling in the Monte Carlo wave-function algorithm used to solve the Lindblad equation. The resulting stochastic gradient estimator exhibits sufficiently low variance and can still be estimated in an embarrassingly parallel manner, enabling efficient scaling of the simulations. We implement this gradient estimator in the existing open-source quarkonium suppression code QTraj. To demonstrate its utility for parameter estimation, we infer the two transport coefficients $\hatκ$ and $\hatγ$ using gradient-based optimization on synthetic nuclear modification factor data.

physics.comp-phCross-listedarXiv:2601.14399•Jan 22, 2026•0 citations
Quantum computing provides a powerful framework for simulating real-time dynamics in open quantum systems, offering key advantages for modeling heavy-quarkonium transport in high-energy nuclear collisions. In this work, we perform quantum simulations of the isotropic next-to-leading-order Lindblad equation for bottomonium in the quark-gluon plasma using a reduced spherical coordinate representation. We discretize operators and wavefunctions, map the physical state onto qubits, and execute time evolution via parameterized quantum gate operations. By extracting the Υ(1S)Υ(1S)Υ(1S) survival probability, we quantitatively isolate the color-octet contribution, demonstrating that its overall impact is small in the final production of the bottomonium ground state Υ(1S)Υ(1S)Υ(1S) in the hot QCD medium at temperatures accessible at the Large Hadron Collider. Additionally, we have further optimized the quantum simulation algorithm for the Lindblad equation. The improved algorithm requires only a single ancillary qubit to realize the Lindblad evolution, thereby minimizing the circuit significantly.
We study the properties of the Lindbladian quantum mechanical evolution of quarkonia with non-Abelian charges (color-singlet and octet) in the quark-gluon plasma. We confirm that heavy quark recoils in the Lindblad equation correctly thermalize quarkonium colorful states within statistical errors from the simulation method. We also demonstrate that the Lindblad equation in the dipole limit can provide an efficient alternative method, which is applicable to a finite time evolution before thermalization and dramatically reduces the numerical cost. Our findings will serve as a foundation for large-scale simulation of quarkonium dynamics in the relativistic heavy-ion collisions.
We study how much quantum structure bottomonium retains during its evolution in the quark-gluon plasma. Within the open quantum system framework, we simulate the Lindblad equation obtained from potential nonrelativistic QCD in the quantum Brownian regime. From the resulting real-time evolution of the bottomonium density matrix, we compute the Wigner transform and study its negativity as a measure of nonclassicality and a test of a key condition underlying classical Langevin-type approximations. We find that the medium suppresses the negativity, which nevertheless approaches a finite plateau close to that of the 1S1\mathrm{S}1S state. The overlaps with the excited 2S2\mathrm{S}2S and 3S3\mathrm{S}3S states arise from near-canceling positive and negative phase-space contributions, driven by negative regions in the Wigner transforms of the bound-state projectors, while the 1S1\mathrm{S}1S overlap is insensitive to such contributions. We further show that the position-space coherence is similarly suppressed and plateaus close to the 1S value. We trace the finite residual quantum structure to conditional purification: the medium preferentially dissolves weakly bound and unbound modes with large ⟨r2⟩\langle r^2\rangle⟨r2⟩, driving the surviving singlet S\mathrm{S}S-wave ensemble toward the 1S1\mathrm{S}1S state. More broadly, these results reveal that classicalization is partial and observable dependent, a feature that is relevant across applications of open quantum system methods in fundamental physics.
Accurate modelling and understanding of quarkonium production in ultrarelativistic heavy ion collisions requires a formalism that preserves the quantum properties of microscopic QˉQ\bar{Q}QQˉ​Q systems while treating the interaction of such pairs with the quark-gluon plasma (QGP). The open quantum system approach has recently emerged as one of the most fruitful schemes to meet such requirements. However, the quantum master equations obtained so far in this context are derived assuming a strict ordering between the QGP temperature TTT and the QˉQ\bar{Q}QQˉ​Q energy gaps (ΔEΔEΔE) of the quarkonia bound states. This limits their predictive power since, as the QGP expands and cools down, the system traverse all regimes between the quantum Brownian motion (QBM) regime for T≳ΔET\gtrsim ΔET≳ΔE to the quantum optical (QO) regime for ΔE≳TΔE\gtrsim TΔE≳T. In this paper, we derive and present a more general non-abelian quantum master equation of the Lindblad type, which does not suffer from these limitations and thus allows to faithfully describe the quantum evolution of the QQˉQ\bar{Q}QQˉ​ pairs during the whole QGP-evolution. We also provide some illustration of the key quantities governing this equation.
Recently, a set of coupled singlet-octet universal Lindblad equations (ULEs) was derived within the framework of non-relativistic QCD (NRQCD) to describe quarkonium dynamics in the quark-gluon plasma (QGP). These equations provide a unified quantum description spanning the quantum Brownian and quantum optical regimes. In this work, we further develop this framework and establish its connection with semiclassical transport. We first derive the universal Lindblad equations within the potential non-relativistic QCD (pNRQCD) effective field theory and show that they coincide with the small-dipole limit of the NRQCD ULEs. We then derive the semiclassical limit of the NRQCD ULEs, obtaining a set of coupled singlet-octet Boltzmann equations. To our knowledge, this is the first derivation of Boltzmann transport equations directly from the universal Lindblad framework. The resulting equations are valid beyond the small-dipole approximation, allowing the evolution of heavy-quark pairs from compact to widely separated configurations. Taking their small-dipole limit allows a direct comparison with the Boltzmann equations of Yao et al. [Phys. Rev. D 99, 096028 (2019)], which were derived within pNRQCD from the Davies secular equation, relying on the rotating-wave approximation (RWA). While the singlet equations are found to be in near-complete agreement, the octet equation contains an additional collision term describing transitions within the continuum of octet scattering states, which is absent from the RWA-based derivation. Finally, we derive the leading quantum correction to the singlet Boltzmann equation. Our results establish a more general and systematic theoretical foundation for the semiclassical transport description of quarkonium in the QGP.
We study the approach to equilibrium of bottomonium in the quark-gluon plasma within the open quantum system framework. We perform large-scale simulations of the long-time behavior in three dimensions using the quantum trajectory method to observe the emergence of steady states and determine the timescale of thermalization in position-, angular-momentum-, and color-space. We find that the thermalization timescale increases with decreasing temperature and decreasing coupling to the medium, which is given by transport coefficients of the medium. Additionally, we observe that the steady states exhibit small corrections to the Gibbs state due to medium interactions and show that these corrections diminish for weaker medium coupling and higher temperature. At a temperature of 450 450\,450MeV, quarkonium relaxes to a state that is approximately thermal, with the most significant correction being a smaller overlap of the 1S1S1S state with respect to the Gibbs state. We compare these findings with the master equation obtained at leading order in the expansion of the binding energy over the temperature, which we find to have a trivial steady state.
We study parameter estimation for the transport coefficients of the quark-gluon plasma by differentiating open-quantum-system-based Monte Carlo simulations of quarkonium suppression. The underlying simulator requires solving a Lindblad equation in a large Hilbert space, which makes parameter estimation computationally expensive. We approach the problem using gradient-based optimization. Specifically, we apply the score-function gradient estimator to differentiate through discrete jump sampling in the Monte Carlo wave-function algorithm used to solve the Lindblad equation. The resulting stochastic gradient estimator exhibits sufficiently low variance and can still be estimated in an embarrassingly parallel manner, enabling efficient scaling of the simulations. We implement this gradient estimator in the existing open-source quarkonium suppression code QTraj. To demonstrate its utility for parameter estimation, we infer the two transport coefficients κ^\hatκκ^ and γ^\hatγγ^​ using gradient-based optimization on synthetic nuclear modification factor data.