Lindblad-driven quarkonium production in heavy-ion collisions
Pith reviewed 2026-05-20 04:15 UTC · model grok-4.3
The pith
The Lindblad equation provides a unified description of quarkonium suppression and recombination in heavy-ion collisions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Starting from the complex-valued in-medium potential, the Lindblad equation is used to derive the dissociation temperature and thermal decay width for each quarkonium state. Survival probabilities are computed for a system undergoing Bjorken expansion. The framework is extended to include recombination from thermalized charm and bottom quarks, deriving a coalescence model from the Lindblad equation under the adiabatic approximation, providing a unified description of suppression and recombination for both charmonium and bottomonium.
What carries the argument
The Lindblad master equation applied to the quarkonium density matrix, using the complex in-medium potential to capture medium effects and deriving dynamics under the adiabatic approximation for recombination.
If this is right
- Survival probabilities for quarkonium states can be calculated consistently during the Bjorken expansion of the quark-gluon plasma.
- A coalescence model for recombination emerges directly from the Lindblad equation in the adiabatic limit.
- The same approach applies to both charmonium and bottomonium systems.
- Thermal decay widths and dissociation temperatures are derived for each state from the potential.
Where Pith is reading between the lines
- If the framework holds, it could allow predictions for quarkonium yields that incorporate both loss and gain mechanisms without ad hoc adjustments.
- Connections to other heavy-ion observables, such as flow or spectra, might be explored by coupling this to hydrodynamic simulations.
- Testing against data at different collision energies could reveal the range of validity of the adiabatic approximation.
Load-bearing premise
The complex-valued in-medium potential is provided as an external input whose real and imaginary parts fully capture the medium interactions with the quarkonium wave function.
What would settle it
If experimental measurements of quarkonium production yields in heavy-ion collisions deviate significantly from the combined predictions of survival probabilities and coalescence rates derived from this Lindblad approach, the framework would be challenged.
Figures
read the original abstract
We study the production of the conventional quarkonium states in ultrarelativistic heavy-ion collisions using an open quantum system framework based on the Lindblad equation. Starting from the complex-valued in-medium potential, we derive the dissociation temperature and thermal decay width for each state, and compute their survival probabilities for a system undergoing Bjorken expansion. We then extend the framework to include recombination from thermalized charm and bottom quarks in the quark-gluon plasma, deriving a coalescence model for quarkonia from the Lindblad equation under the adiabatic approximation. The methodology provides a unified, first-principles-inspired description of suppression and recombination for both charmonium and bottomonium.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops an open quantum system framework based on the Lindblad equation to model quarkonium production in heavy-ion collisions. Starting from an externally supplied complex in-medium potential, it derives dissociation temperatures and thermal decay widths for charmonium and bottomonium states, computes survival probabilities during Bjorken expansion, and extends the approach to a coalescence model for recombination under the adiabatic approximation, aiming for a unified description of suppression and regeneration.
Significance. If the derivations hold, the work provides a coherent dynamical framework linking microscopic Lindblad evolution to macroscopic observables like survival probabilities and coalescence yields for both charm and bottom sectors. This could offer improved consistency between suppression and regeneration mechanisms compared to separate phenomenological treatments, though its predictive power remains tied to the accuracy of the input potential.
major comments (3)
- [Abstract and §2] Abstract and §2: The central claim of a 'unified, first-principles-inspired description' rests on inserting a pre-existing complex potential V(r,T) whose real part determines binding and imaginary part supplies the width. Because this potential is not derived or constrained from the Lindblad master equation itself, the framework's unification is only as robust as the external input; any medium effects outside this V(r,T) fall outside the Lindblad dynamics.
- [§4] §4 (survival probabilities): The integrals for survival under Bjorken expansion inherit all uncertainties from the externally supplied potential and the adiabatic approximation; without a quantitative assessment of how variations in Re[V] and Im[V] propagate to the final yields, the claimed improvement over standard models cannot be evaluated.
- [§5] §5 (coalescence model): The derivation of the coalescence rate from the Lindblad equation under adiabatic approximation assumes the potential remains valid for the recombining quark-antiquark pair; this assumption is load-bearing for the recombination predictions but lacks an explicit check against the regime where the approximation breaks.
minor comments (2)
- [§2] Notation for the complex potential should be introduced with explicit real and imaginary parts at first use to avoid ambiguity in later sections.
- [Figures 3-5] Figure captions for survival probability plots should include the specific parameter values of the input potential used.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below, clarifying the role of the input potential within the Lindblad framework and agreeing to strengthen the discussion of approximations and uncertainties where appropriate.
read point-by-point responses
-
Referee: [Abstract and §2] Abstract and §2: The central claim of a 'unified, first-principles-inspired description' rests on inserting a pre-existing complex potential V(r,T) whose real part determines binding and imaginary part supplies the width. Because this potential is not derived or constrained from the Lindblad master equation itself, the framework's unification is only as robust as the external input; any medium effects outside this V(r,T) fall outside the Lindblad dynamics.
Authors: We agree that the complex in-medium potential V(r,T) is supplied as an external input, consistent with standard effective-theory approaches to quarkonium in the QGP. The Lindblad master equation then provides the unified dynamical evolution: the real part enters the coherent Hamiltonian term governing binding and oscillation, while the imaginary part generates the dissipative jump operators responsible for decoherence and dissociation. This structure allows both suppression (via the survival probability under Bjorken flow) and recombination (via the derived coalescence rate) to be obtained from the same master equation under the adiabatic approximation. The 'first-principles-inspired' qualifier refers to the open-quantum-system derivation rather than a first-principles computation of V(r,T) itself, which remains an active research topic. We have revised the abstract and §2 to make this distinction explicit. revision: partial
-
Referee: [§4] §4 (survival probabilities): The integrals for survival under Bjorken expansion inherit all uncertainties from the externally supplied potential and the adiabatic approximation; without a quantitative assessment of how variations in Re[V] and Im[V] propagate to the final yields, the claimed improvement over standard models cannot be evaluated.
Authors: We acknowledge that a systematic propagation of uncertainties from the input potential would allow a more quantitative comparison with phenomenological models. In the present work we adopted a representative parametrization of V(r,T) from the literature and demonstrated the framework's consistency across charmonium and bottomonium sectors. In the revised manuscript we will add a short sensitivity study showing how variations in the strength of Im[V] affect the thermal widths and integrated survival probabilities, while noting that a full Monte-Carlo error propagation lies beyond the scope of this initial study. revision: yes
-
Referee: [§5] §5 (coalescence model): The derivation of the coalescence rate from the Lindblad equation under adiabatic approximation assumes the potential remains valid for the recombining quark-antiquark pair; this assumption is load-bearing for the recombination predictions but lacks an explicit check against the regime where the approximation breaks.
Authors: The adiabatic approximation is invoked when the relative motion of the recombining pair is slow compared with the medium relaxation timescale, which holds for thermalized heavy quarks at the temperatures and densities relevant to our Bjorken evolution. We have verified that the condition is satisfied for the momentum scales considered. To address the referee's concern we will insert a brief paragraph in §5 that quantifies the adiabaticity criterion (ratio of expansion rate to the inverse relaxation time) and indicates the temperature and momentum range where the approximation remains reliable. revision: yes
Circularity Check
No significant circularity; derivations proceed forward from external potential input via Lindblad dynamics
full rationale
The paper explicitly starts from a complex-valued in-medium potential supplied as an external input and applies the Lindblad equation plus adiabatic approximation to derive dissociation temperatures, thermal decay widths, Bjorken survival probabilities, and a coalescence model. No load-bearing step reduces any output quantity back to a fitted parameter or prior result by construction, nor does any self-citation chain or ansatz smuggling appear in the derivation. The framework remains self-contained once the potential is given, with all subsequent quantities obtained mathematically from that starting point rather than being equivalent to it by definition or statistical forcing.
Axiom & Free-Parameter Ledger
free parameters (1)
- parameters of the complex in-medium potential
axioms (2)
- domain assumption The Lindblad equation governs the time evolution of the quarkonium density matrix in the quark-gluon plasma.
- ad hoc to paper The adiabatic approximation is valid when deriving the coalescence model from the Lindblad equation.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Starting from the complex-valued in-medium potential, we derive the dissociation temperature and thermal decay width... deriving a coalescence model for quarkonia from the Lindblad equation under the adiabatic approximation.
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The imaginary part of the potential arises from Landau damping... ImV(r,mD(T)) = -αs T ϕ(mDr) + ...
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Tetsuo Matsui and Helmut Satz.J/ψSuppression by Quark-Gluon Plasma Formation.Phys. Lett. B, 178:416,
-
[2]
doi: 10.1016/0370-2693(86)91404-8
-
[3]
Anton Andronic et al. Heavy-flavour and quarkonium pro- duction in the LHC era: from proton-proton to heavy- ion collisions.Eur . Phys. J. C, 76(3):107, 2016. doi: 10.1140/epjc/s10052-015-3819-5
-
[4]
Real-time static potential in hot QCD
Mikko Laine, Owe Philipsen, Paul Romatschke, and Markus Tassler. Real-time static potential in hot QCD. JHEP, 03:054, 2007. doi: 10.1088/1126-6708/2007/03/ 054
-
[5]
A. Beraudo, J.-P. Blaizot, and C. Ratti. Real and imaginary-time Q anti-Q correlators in a thermal medium. Nucl. Phys. A, 806:312, 2008. doi: 10.1016/j.nuclphysa. 2008.03.001
-
[6]
Improved Gauss law model and in-medium heavy quarkonium at finite den- sity and velocity.Phys
David Lafferty and Alexander Rothkopf. Improved Gauss law model and in-medium heavy quarkonium at finite den- sity and velocity.Phys. Rev. D, 101(5):056010, 2020. doi: 10.1103/PhysRevD.101.056010
-
[7]
Stochas- tic potential and quantum decoherence of heavy quarko- nium in the quark-gluon plasma.Phys
Yukinao Akamatsu and Alexander Rothkopf. Stochas- tic potential and quantum decoherence of heavy quarko- nium in the quark-gluon plasma.Phys. Rev. D, 85:105011,
-
[8]
doi: 10.1103/PhysRevD.85.105011
-
[9]
Quarkonium suppression in heavy-ion col- lisions: an open quantum system approach.Phys
Nora Brambilla, Miguel Angel Escobedo, Joan Soto, and Antonio Vairo. Quarkonium suppression in heavy-ion col- lisions: an open quantum system approach.Phys. Rev. D, 96(3):034021, 2017. doi: 10.1103/PhysRevD.96.034021
-
[10]
Quan- tum and classical dynamics of heavy quarks in a quark- gluon plasma.JHEP, 06:034, 2018
Jean-Paul Blaizot and Miguel Angel Escobedo. Quan- tum and classical dynamics of heavy quarks in a quark- gluon plasma.JHEP, 06:034, 2018. doi: 10.1007/ JHEP06(2018)034
work page 2018
-
[11]
Heavy quarkonium suppression in a fire- ball.Phys
Nora Brambilla, Miguel Angel Escobedo, Joan Soto, and Antonio Vairo. Heavy quarkonium suppression in a fire- ball.Phys. Rev. D, 97(7):074009, 2018. doi: 10.1103/ PhysRevD.97.074009. 8
work page 2018
-
[12]
Fer- reiro, and Víctor López-Pardo
Néstor Armesto, Miguel Angel Escobedo, Elena G. Fer- reiro, and Víctor López-Pardo. A potential approach to the X(3872) thermal behavior.Phys. Lett. B, 854:138760,
-
[13]
doi: 10.1016/j.physletb.2024.138760
-
[14]
Fer- reiro, and Víctor López-Pardo
Néstor Armesto, Miguel Ángel Escobedo, Elena G. Fer- reiro, and Víctor López-Pardo. Lindblad-driven recombi- nation of the X(3872) tetraquark. 2026
work page 2026
-
[15]
Evidence for X(3872) in Pb- Pb Collisions and Studies of its Prompt Production at√sNN =5.02 TeV.Phys
Albert M Sirunyan et al. Evidence for X(3872) in Pb- Pb Collisions and Studies of its Prompt Production at√sNN =5.02 TeV.Phys. Rev. Lett., 128(3):032001, 2022. doi: 10.1103/PhysRevLett.128.032001
-
[16]
Estia Eichten, Kurt Gottfried, Toichiro Kinoshita, Ken- neth D. Lane, and Tung-Mow Yan. Charmonium: The Model.Phys. Rev. D, 17:3090, 1978. doi: 10.1103/ PhysRevD.17.3090
work page 1978
-
[17]
Determination of the bottom quark mass from the Upsilon(1S) system.JHEP, 06:022, 2001
Antonio Pineda. Determination of the bottom quark mass from the Upsilon(1S) system.JHEP, 06:022, 2001. doi: 10.1088/1126-6708/2001/06/022
-
[18]
Heavy quarkonium in any channel in resummed hot QCD
Yannis Burnier, Mikko Laine, and Mikael Vepsalainen. Heavy quarkonium in any channel in resummed hot QCD. JHEP, 01:043, 2008. doi: 10.1088/1126-6708/2008/01/ 043
-
[19]
James D. Bjorken. Highly Relativistic Nucleus-Nucleus Collisions: The Central Rapidity Region.Phys. Rev. D, 27:140, 1983. doi: 10.1103/PhysRevD.27.140
-
[20]
Gen- eral theory of slow non-Hermitian evolution
Parveen Kumar, Yuval Gefen, and Kyrylo Snizhko. Gen- eral theory of slow non-Hermitian evolution. 2 2025
work page 2025
-
[21]
Loic Grandchamp, Ralf Rapp, and Gerald E. Brown. In medium effects on charmonium production in heavy ion collisions.Phys. Rev. Lett., 92:212301, 2004. doi: 10. 1103/PhysRevLett.92.212301
work page 2004
-
[22]
Char- monium production from nonequilibrium charm and an- ticharm quarks in quark-gluon plasma.Phys
Taesoo Song, Kyong Chol Han, and Che Ming Ko. Char- monium production from nonequilibrium charm and an- ticharm quarks in quark-gluon plasma.Phys. Rev. C, 85: 054905, 2012. doi: 10.1103/PhysRevC.85.054905
-
[23]
Charmonium Transport in Heavy-Ion Collisions at the LHC.Universe, 10(6):244,
Biaogang Wu and Ralf Rapp. Charmonium Transport in Heavy-Ion Collisions at the LHC.Universe, 10(6):244,
-
[24]
doi: 10.3390/universe10060244
-
[25]
Albert Messiah.Quantum Mechanics. Dover Publica- tions, 1999. Two volumes bound as one
work page 1999
-
[26]
Miguel Angel Escobedo and Elena G. Ferreiro. Sim- ple model to include initial-state and hot-medium ef- fects in the computation of quarkonium nuclear modifi- cation factor.Phys. Rev. D, 105(1):014019, 2022. doi: 10.1103/PhysRevD.105.014019
-
[27]
Shreyasi Acharya et al. Charm-quark fragmentation frac- tions and production cross section at midrapidity in pp collisions at the LHC.Phys. Rev. D, 105(1):L011103,
-
[28]
doi: 10.1103/PhysRevD.105.L011103
-
[29]
Shreyasi Acharya et al. Measurement of beauty and charm production in pp collisions at √s=5.02 TeV via non- prompt and prompt D mesons.JHEP, 05:220, 2021. doi: 10.1007/JHEP05(2021)220
-
[30]
Inclusive J/ψproduction at mid- rapidity in pp collisions at √s=5.02 TeV.JHEP, 10: 084, 2019
Shreyasi Acharya et al. Inclusive J/ψproduction at mid- rapidity in pp collisions at √s=5.02 TeV.JHEP, 10: 084, 2019. doi: 10.1007/JHEP10(2019)084
-
[31]
Albert M Sirunyan et al. Suppression of ExcitedΥ States Relative to the Ground State in Pb-Pb Collisions at √sNN =5.02 TeV.Phys. Rev. Lett., 120:142301, 2018. doi: 10.1103/PhysRevLett.120.142301
-
[32]
(Non)thermal aspects of charmonium production and a new look at J/ψsuppression.Phys
Peter Braun-Munzinger and Johanna Stachel. (Non)thermal aspects of charmonium production and a new look at J/ψsuppression.Phys. Lett. B, 490: 196, 2000. doi: 10.1016/S0370-2693(00)00991-6
-
[33]
Variational Monte Carlo calculations of n+H3 scattering
Robert L. Thews, Martin Schroedter, and Johann Rafel- ski. EnhancedJ/ψproduction in deconfined quark matter. Phys. Rev. C, 63:054905, 2001. doi: 10.1103/PhysRevC. 63.054905
-
[34]
Measurement of prompt and nonprompt charmonium suppression in PbPb collisions at 5.02 TeV.Eur
Albert M Sirunyan et al. Measurement of prompt and nonprompt charmonium suppression in PbPb collisions at 5.02 TeV.Eur . Phys. J. C, 78(6):509, 2018. doi: 10.1140/epjc/s10052-018-5950-6. 9
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.