pith. sign in

arxiv: 2605.19985 · v1 · pith:G4HHSOBVnew · submitted 2026-05-19 · ✦ hep-ph · hep-ex· nucl-ex· nucl-th

Lindblad-driven quarkonium production in heavy-ion collisions

Pith reviewed 2026-05-20 04:15 UTC · model grok-4.3

classification ✦ hep-ph hep-exnucl-exnucl-th
keywords quarkoniumheavy-ion collisionsLindblad equationsuppressionrecombinationquark-gluon plasmacharmoniumbottomonium
0
0 comments X

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.

The paper aims to show that an open quantum system treatment based on the Lindblad equation can handle both the dissociation and reformation of quarkonium states in the hot medium of heavy-ion collisions. It starts with a complex in-medium potential to find dissociation temperatures and decay widths for states like charmonium and bottomonium. Survival probabilities are then calculated as the system expands according to Bjorken hydrodynamics. Recombination is incorporated by deriving a coalescence model from the same Lindblad framework in the adiabatic limit. A sympathetic reader would care because this offers a consistent, single-framework way to predict how these bound states are affected by the quark-gluon plasma, rather than using separate models for suppression and regeneration.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.19985 by Elena G. Ferreiro, Miguel \'Angel Escobedo, N\'estor Armesto, V\'ictor L\'opez-Pardo.

Figure 1
Figure 1. Figure 1: Real part of the in-medium potential for quarkonia as given by (3). It [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Imaginary part of the in-medium potential for quarkonia as given [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: From this figure, it is evident that the binding en￾ergy of the J/ψ state can be tracked over nearly the entire tem￾perature range, whereas the corresponding signal for the ψ(2S ) state is barely visible. This is due to the fact that the fit parame￾ters of the ψ(2S ) spectral function cannot be reliably extracted for temperatures T ≳ 155 MeV. As observed in [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: In-medium S -wave spectral functions for vector channel charmonium for different temperatures. The dashed gray vertical lines represent the T = 0 bound states J/ψ and ψ(2S ). As the temperature increases, these peaks gradually shift to￾ward lower frequencies and eventually fade away. At the same time, they become broader and less pronounced, indicating a shorter lifetime of the states. In the vacuum, the s… view at source ↗
Figure 5
Figure 5. Figure 5: Decay width for the different species of charmonium at finite temper￾ature: J/ψ (solid purple line) and ψ(2S ) (dash-dot blue line). The ψ(2S ) line is barely visible because the corresponding peak vanishes at T ≳ 155 MeV. The dissociation temperature Td is defined as the tempera￾ture above which the Schrödinger equation no longer admits 3 [PITH_FULL_IMAGE:figures/full_fig_p003_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Binding energy for the different species of bottomonium at finite temperature: Υ(1S ) (solid blue line), Υ(2S ) (dash-dot magenta line), and Υ(3S ) (dotted orange line). The potential barrier is also shown (dashed gray line). Moreover, increasing the temperature leads to a broadening of the peaks, corresponding to an increase in the decay width. The decay width is therefore a monotonically increasing func￾… view at source ↗
Figure 6
Figure 6. Figure 6: In-medium S -wave spectral functions for vector channel bottomo￾nium for different temperatures. The dashed gray vertical lines represent the T = 0 bound states: Υ(1S ), Υ(2S ), Υ(3S ) and Υ(4S ). monium case, the spectral function in the vicinity of each peak can be fitted with a skewed Breit–Wigner distribution, allow￾ing for the extraction of both the binding energy and the decay width. For bottomonium,… view at source ↗
Figure 9
Figure 9. Figure 9: Nuclear modification factor RAA of the J/ψ in Pb–Pb collisions at √ sNN = 5.02 TeV. The dotted curve shows only cold nuclear matter effects (shadowing); the dashed curve adds suppression from the complex Gauss-law potential; the dash-dot curve shows the recombination contribution alone; and the solid curve gives the total result. 5.2. Results for J/ψ [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: Nuclear modification factor RAA of the J/ψ and Υ(1S ) in Pb–Pb collisions at √ sNN = 5.02 TeV. Experimental data from ALICE [25] and CMS [26, 29] are shown for comparison. . ground states J/ψ and Υ(1S ) in Pb-Pb collisions at √ sNN = 5.02 TeV. Our approach proceeds as follows: First, we use the Gauss-law model of Lafferty and Rothkopf [5], which pro￾vides a single-parameter (mD(T)) parametrization of both… view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

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)
  1. [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.
  2. [§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.
  3. [§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)
  1. [§2] Notation for the complex potential should be introduced with explicit real and imaginary parts at first use to avoid ambiguity in later sections.
  2. [Figures 3-5] Figure captions for survival probability plots should include the specific parameter values of the input potential used.

Simulated Author's Rebuttal

3 responses · 0 unresolved

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
  1. 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

  2. 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

  3. 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

0 steps flagged

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

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the Lindblad equation as the correct master equation for the open system, the validity of the adiabatic approximation when adding recombination, and the assumption that a complex in-medium potential fully encodes the medium interactions. No new particles or forces are postulated.

free parameters (1)
  • parameters of the complex in-medium potential
    Real and imaginary parts of the potential are taken as input; their specific functional form and temperature dependence are not derived inside the paper.
axioms (2)
  • domain assumption The Lindblad equation governs the time evolution of the quarkonium density matrix in the quark-gluon plasma.
    Invoked at the outset as the starting point for both dissociation and recombination calculations.
  • ad hoc to paper The adiabatic approximation is valid when deriving the coalescence model from the Lindblad equation.
    Explicitly stated as the condition under which recombination is added.

pith-pipeline@v0.9.0 · 5653 in / 1632 out tokens · 30250 ms · 2026-05-20T04:15:35.079362+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

34 extracted references · 34 canonical work pages

  1. [1]

    Tetsuo Matsui and Helmut Satz.J/ψSuppression by Quark-Gluon Plasma Formation.Phys. Lett. B, 178:416,

  2. [2]

    doi: 10.1016/0370-2693(86)91404-8

  3. [3]

    Heavy-flavour and quarkonium pro- duction in the LHC era: from proton-proton to heavy- ion collisions.Eur

    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. [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. [5]

    Beraudo, J.-P

    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. [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. [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. [8]

    doi: 10.1103/PhysRevD.85.105011

  9. [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. [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

  11. [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

  12. [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. [13]

    doi: 10.1016/j.physletb.2024.138760

  14. [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

  15. [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. [16]

    Lane, and Tung-Mow Yan

    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

  17. [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. [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. [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. [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

  21. [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

  22. [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. [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. [24]

    doi: 10.3390/universe10060244

  25. [25]

    Dover Publica- tions, 1999

    Albert Messiah.Quantum Mechanics. Dover Publica- tions, 1999. Two volumes bound as one

  26. [26]

    Ferreiro

    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. [27]

    Charm-quark fragmentation frac- tions and production cross section at midrapidity in pp collisions at the LHC.Phys

    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. [28]

    doi: 10.1103/PhysRevD.105.L011103

  29. [29]

    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

    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. [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. [31]

    Suppression of ExcitedΥ States Relative to the Ground State in Pb-Pb Collisions at √sNN =5.02 TeV.Phys

    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. [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. [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. [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