The reviewed record of science sign in
Pith

arxiv: 2607.06137 · v1 · pith:AOVJRIXU · submitted 2026-07-07 · hep-ph · nucl-th· quant-ph

Quantum decoherence: a study applied to quarkonium-like bound states in strongly interacting matter

Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel 2026-07-08 15:22 UTCglm-5.2pith:AOVJRIXUrecord.jsonopen to challenge →

classification hep-ph nucl-thquant-ph
keywords quantumdecoherenceevolutioninteractingboundcouplingmattermedium
0
0 comments X

The pith

Expanding plasma slows quarkonium decoherence, viscosity barely matters

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper models a J/psi charmonium bound state as a quantum harmonic oscillator coupled to a thermal bath representing the quark-gluon plasma, then evolves the system's density matrix under a Lindblad master equation to track how quantum coherence is lost. The central move is extending the standard static-bath Lindblad framework to include a time-dependent temperature T(t) that mimics the one-dimensional hydrodynamic expansion of the fireball produced in heavy-ion collisions. The authors find that when the bath cools as the plasma expands, decoherence proceeds more slowly than in a static bath held at the initial temperature. For LHC-like conditions (initial temperature 0.5 GeV), coherence is lost within roughly 2 fm/c, while for RHIC-like conditions (initial temperature 0.3 GeV), coherence persists beyond 5 fm/c. Viscous corrections to the expansion, modeled at eta/s = 1/4pi, produce negligible change in the decoherence timescale compared to ideal hydrodynamics. The paper also tracks the Wigner function in phase space, showing that negative fringes characteristic of quantum superposition are washed out on the same decoherence timescale, driving the system toward a classical probability distribution.

Core claim

The paper's central result is that introducing a realistic time-dependent temperature profile for the expanding QGP into the Lindblad master equation for a J/psi-like harmonic oscillator slows quantum decoherence relative to a static bath, and that this effect depends strongly on the initial temperature (RHIC vs. LHC conditions) but is insensitive to viscous corrections. The mechanism is straightforward: as the plasma expands and cools, the system-reservoir coupling strength gamma(T) decreases, so the rate at which off-diagonal density matrix elements are suppressed drops over time. At LHC temperatures the initial coupling is strong enough that coherence is destroyed in about 2 fm/c, before;

What carries the argument

Lindblad master equation in the optical regime for a damped harmonic oscillator with amplitude-coupling to a bosonic reservoir, extended to include time-dependent temperature T(t) via a multiplicative factor f(t) on the coupling coefficients. The system frequency omega_0 = 0.457 GeV is fixed from the J/psi root-mean-square radius. The dissociation rate gamma(T) is taken from phenomenological charmonium suppression calculations. Temperature evolution follows a Bjorken-like power law T(t) = T_0 * (tau_0/t)^(c_s^2), with c_s^2 = 1/3 for ideal hydrodynamics or a Boltzmann transport result for viscous expansion.

If this is right

  • If decoherence at RHIC temperatures takes longer than 5 fm/c, quarkonium states may retain partial quantum coherence through much of the fireball lifetime, meaning hadronization models that assume fully classical phase-space distributions may need revision for systems produced at lower initial temperatures.
  • The near-equivalence of ideal and viscous expansion for decoherence suggests that the dominant control parameter is the initial temperature, not the detailed hydrodynamic profile, simplifying future calculations.
  • The Wigner function evolution provides a direct bridge between the open-quantum-systems description and coalescence hadronization models, which assume Gaussian Wigner functions; the paper shows this assumption becomes valid on the decoherence timescale.
  • The framework can be extended to other quarkonium states (Upsilon, chi_c) by changing omega_0 and gamma(T), potentially yielding a hierarchy of decoherence times across the quarkonium spectrum.

Where Pith is reading between the lines

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

  • If the oscillator frequency omega_0 were allowed to vary with temperature (as it physically should, since in-medium binding energy changes), the decoherence timescales could shift substantially, particularly near T_c where color screening is strongly temperature-dependent. The paper's quantitative results are thus best read as order-of-magnitude estimates for the optical regime rather than precise
  • The observation that coherence survives beyond 5 fm/c at RHIC conditions raises the question of whether partially coherent quarkonium states could serve as quantum probes of the early-time QGP, carrying phase information from the initial production mechanism through the medium evolution.
  • The insensitivity to viscosity suggests that decoherence is controlled primarily by the integrated thermal exposure (the time-integrated coupling strength) rather than the instantaneous expansion rate, which could be tested by comparing different equations of state at fixed initial temperature.

Load-bearing premise

The oscillator frequency omega_0 representing the J/psi binding energy is held constant throughout the evolution, even as the temperature drops from 0.5 GeV to near the phase transition. In reality, the in-medium binding energy and thus the oscillator frequency should change substantially with temperature, particularly near the deconfinement transition where color screening is strongly temperature-dependent. The authors flag this themselves.

What would settle it

If allowing omega_0 to vary with temperature in a self-consistent way changed the decoherence timescale by a factor of two or more, the quantitative distinction between RHIC and LHC conditions would need to be re-established.

Figures

Figures reproduced from arXiv: 2607.06137 by Gabriele Coci, Giuseppe Falci, Salvatore Plumari.

Figure 1
Figure 1. Figure 1: FIG. 1: (color online): Population of [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Time evolution of imaginary part of coherence ele [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Time evolution of the real part of coherence ele [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: (color online): Time evolution of the real part of [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Contour plot of Wigner distribution for initial coher [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Time evolution of Wigner function [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Time evolution of momentum distribution [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Average momentum [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: (Color online) Average relative momentum [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Contour plot of initial cat state Eq. (43) with pa [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Contour plots of Wigner function [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Time evolution of ground state population of [PITH_FULL_IMAGE:figures/full_fig_p012_12.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: (color online): Time evolution of the real part of [PITH_FULL_IMAGE:figures/full_fig_p013_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: (color online): Time evolution of the imaginary part [PITH_FULL_IMAGE:figures/full_fig_p013_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16: Contour plots of Wigner function [PITH_FULL_IMAGE:figures/full_fig_p014_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17: (color online) The real part of coherence elements [PITH_FULL_IMAGE:figures/full_fig_p014_17.png] view at source ↗
read the original abstract

We study the quantum decoherence of a bound state interacting with a reservoir of strongly interacting matter within the framework of open quantum systems. The bound state is modeled as a quantum harmonic oscillator whose parameters are tuned to reproduce the root-mean-square radius of $J/\Psi$ particle. The surrounding medium, representing the many degrees of freedom of strongly interacting matter, acts as an environment that induces dissipation and decoherence through system-reservoir coupling. By analyzing the time evolution of the reduced density matrix, we quantify the loss of quantum coherence and its dependence on medium properties. Subsequently, we extend the model by introducing a time dependence in the system-thermal bath coupling, thereby simulating a temperature evolution similar to that occurring during the expansion of a fireball in the central region of heavy-ion collisions. We find that a temperature evolution has a relevant impact on the way the system loses coherence through the coupling with the expanding medium. Finally, we estimate the impact of the time-dependent temperature on the decoherence process, also analyzing a scenario that includes viscous effects without finding a significant change with respect to ideal hydrodynamical evolution.

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 / 7 minor

Summary. The manuscript studies quantum decoherence of a J/ψ-like harmonic oscillator coupled to a thermal reservoir representing the QGP, within the open quantum systems framework. The oscillator frequency ω₀ is fixed from the J/ψ RMS radius, and the dissociation rate γ(T) is taken from the TAMU model. The authors derive a Lindblad master equation in the optical regime, extend it to include a time-dependent temperature T(t) modeling 1D Bjorken expansion, and solve it numerically. The central results are: (1) in a static bath, coherence is lost on timescales of ~5 fm/c at T = 0.3 GeV; (2) with expanding medium, decoherence is slowed relative to the static case; (3) for LHC conditions (T₀ = 0.5 GeV), coherence is lost within ~2 fm/c, while for RHIC conditions (T₀ = 0.3 GeV), it persists beyond ~5 fm/c; (4) viscous corrections (η/s = 1/4π) do not significantly alter the decoherence timescale. The Wigner function evolution from cat-state initial conditions is also presented, showing classicalization and connection to coalescence hadronization models.

Significance. The paper addresses a timely question — how fast quantum coherence of quarkonium is lost in the QGP — using a well-established Lindblad framework in the optical regime. The extension to a time-dependent temperature is a natural and useful step toward realism. The model has a small number of physically motivated inputs (ω₀ from the J/ψ radius, γ(T) from an external dissociation calculation), and the decoherence timescales are genuine outputs rather than fitted parameters. The Wigner-function analysis connecting decoherence to the validity of classical coalescence hadronization is a nice conceptual contribution. The finding that viscous corrections are negligible for decoherence is a concrete, falsifiable result. The work is exploratory rather than definitive, but it provides a reasonable first estimate and a framework that can be refined.

major comments (3)
  1. Sec. VI, and Eqs. (32): The central quantitative claims (coherence lost within ~2 fm/c at LHC, persisting beyond ~5 fm/c at RHIC) depend on the decoherence rate Re(a_n) = (γ/2)·n·(1+2n̄(ω₀,T)), where n̄(ω₀,T) = 1/(e^{ω₀/T}−1). The authors fix ω₀ = 0.457 GeV (vacuum J/ψ value) throughout the evolution, even as T drops from 0.5 GeV toward T_c ≈ 0.16 GeV. In reality, the in-medium binding energy — and thus ω₀(T) — varies substantially with temperature, particularly near T_c where color screening is strongly T-dependent. Since n̄ enters the decoherence rate exponentially, a T-dependent ω₀(T) could shift the quoted timescales significantly: at early times when T₀ = 0.5 GeV, a smaller ω₀(T₀) (weaker binding) would increase n̄ and accelerate decoherence; near T_c, a larger ω₀(T_c) (stronger binding restored) would suppress n̄ exponentially and extend decoherence well beyond ~5 fm/c. The authors
  2. Sec. II.A, Eq. (17): The optical regime requires τ_S ≪ τ_R, i.e., ω₀ ≫ γ(T). With ω₀ = 0.457 GeV, this translates to γ(T) ≪ 0.457 GeV. The authors should verify and state explicitly that this condition is satisfied for the γ(T) values used at all temperatures considered (T = 0.2–0.5 GeV), particularly at the highest T where γ(T) is largest. If γ(T) approaches ω₀ at any point, the rotating-wave approximation underlying Eq. (14) breaks down, and the quantitative results in that regime would be unreliable. A figure or table of γ(T)/ω₀ versus T would address this.
  3. Sec. II.B, Eq. (22)–(23): The quasi-static approximation T(t') ≈ T(t) requires that the bath correlation time τ_E be much shorter than the hydrodynamic timescale τ_hydro. The authors state this hierarchy but do not provide estimates of τ_E. For the optical regime, τ_E ~ 1/T, which at T = 0.3 GeV gives τ_E ~ 0.66 fm, while τ_hydro ~ τ₀ = 0.6 fm — a ratio of only ~1, not a clean separation. The authors should either provide a more careful justification (e.g., estimating τ_E from the spectral density J(ω) rather than from 1/T) or acknowledge that the quasi-static approximation is marginal at the lower end of the temperature range and discuss how this affects the RHIC results in particular. This is load-bearing because the entire time-dependent T(t) framework — and hence the claim that expansion slows decoherence — rests on this approximation.
minor comments (7)
  1. Sec. II.C: The value σ_r = 0.348 fm is attributed to Table 1 of Ref. [19], but the inversion of Eq. (25) to obtain ω₀ = 0.457 GeV should show the intermediate steps (the reduced mass μ = M_c/2 = 0.7 GeV) so the reader can verify the arithmetic.
  2. Eq. (30): The approximation retaining only terms proportional to n ≫ 1 is stated to work well for n ≥ 4, with 25–40% discrepancy for n = 1,2. Since the n = 1 coherence ρ₀₁ is one of the primary observables shown in the figures, the authors should note this limitation when presenting those results.
  3. Fig. 8: The caption states that γ obtained from the exponential fit is 'equivalent to the input dissociation rate.' This is expected by construction (the Lindblad equation is built from γ), so the statement should be framed as a consistency check rather than a prediction.
  4. Sec. VI: The statement 'the bounded system, characterized by frequency ω₀ increases as the temperature decreases' is confusingly worded — ω₀ is held constant throughout. The intended meaning appears to be that the ratio ω₀/T increases. Please rephrase.
  5. Sec. V.B: The connection to coalescence hadronization (paragraph beginning 'On the one hand, Q̄Q pairs...') is interesting but speculative. A brief statement that this is a qualitative observation motivating future work, rather than an established result, would improve precision.
  6. References: Several arXiv-only citations (e.g., Refs. [35], [36], [76], [79]) should be updated to published versions if available.
  7. Typo in Sec. II.C: 'indipendent' should be 'independent' (also appears in Sec. II.A, Eq. (2) description). Similarly, 'reab-sorbed' in Sec. II.B should be 'reabsorbed,' and 'parametrization' is used inconsistently with 'parametrization' elsewhere.

Simulated Author's Rebuttal

3 responses · 2 unresolved

We thank the referee for a careful and constructive report. The comments identify three important physical issues — the temperature dependence of the binding frequency, the validity of the optical regime, and the quasi-static approximation — that we agree must be addressed more explicitly in the manuscript. Below we respond to each point. We find that two of the three require revision (one partial, one full), and for the third we provide the requested quantitative verification, which also necessitates a revision to make the check explicit in the text.

read point-by-point responses
  1. Referee: Major Comment 1 (Sec. VI, Eqs. 32): The central quantitative claims depend on fixing omega_0 = 0.457 GeV (vacuum J/psi value) throughout the evolution, even as T drops from 0.5 GeV toward T_c. In reality, the in-medium binding energy — and thus omega_0(T) — varies substantially with temperature, particularly near T_c where color screening is strongly T-dependent. Since n_bar enters the decoherence rate exponentially, a T-dependent omega_0(T) could shift the quoted timescales significantly.

    Authors: The referee is correct that a temperature-dependent omega_0(T) could quantitatively shift the decoherence timescales, and we agree this limitation must be stated more prominently. In the current manuscript, we acknowledge on p. 13 that 'this difference depends strongly on the way the dissociation coefficient gamma(T) varies with temperature, as well as on our assumption that the frequency omega_0 remains unaffected by the environment,' and we list the introduction of a T-dependent omega_0 as a future direction in the Conclusions. However, we agree this is insufficiently emphasized given its quantitative impact. We will revise the manuscript to include a dedicated discussion of this point in Sec. VI, explicitly noting the direction of the bias the referee identifies: at early times (high T), a smaller in-medium omega_0 would increase n_bar and accelerate decoherence, while near T_c, a larger omega_0 would suppress n_bar exponentially and extend decoherence. We will also add a quantitative estimate of the sensitivity. That said, we note that a fully self-consistent omega_0(T) is not straightforward to implement within the current framework: the Lindblad equation in the optical regime requires omega_0 to define the system Hamiltonian and the Fock basis at all times, and introducing a time-dependent omega_0(t) would require re-deriving the master equation with a time-dependent system frequency, which involves additional non-trivial terms (parametric driving) beyond the current amplitude-coupling Lindblad form. This is a genuine technical limitation of the present framework, not merely a parameter choice. We will therefore implement a partial revision: adding the sensitivity discussion and the caveat, while deferring the full T-dependent omega_0 treatment to a future study,, revision: partial

  2. Referee: Major Comment 2 (Sec. II.A, Eq. 17): The optical regime requires tau_S << tau_R, i.e., omega_0 >> gamma(T). With omega_0 = 0.457 GeV, this translates to gamma(T) << 0.457 GeV. The authors should verify and state explicitly that this condition is satisfied for the gamma(T) values used at all temperatures considered (T = 0.2-0.5 GeV), particularly at the highest T where gamma(T) is largest. A figure or table of gamma(T)/omega_0 versus T would address this.

    Authors: We agree that this verification should be explicit in the manuscript. We have checked the condition gamma(T) << omega_0 for the TAMU model dissociation rates used in this work (from Grandchamp and Rapp, Refs. [18, 58]). At T = 0.5 GeV (the highest temperature considered), the dissociation rate gamma is approximately 0.08-0.1 GeV, giving gamma/omega_0 ~ 0.2. At T = 0.4 GeV, gamma/omega_0 ~ 0.1, and at T = 0.2 GeV the ratio is well below 0.05. The condition gamma(T) << omega_0 is thus satisfied throughout the temperature range, though we acknowledge the ratio is largest at the highest temperatures and the separation is not as dramatic as in typical quantum optics applications. We will add a table or figure of gamma(T)/omega_0 versus T in the revised manuscript, along with an explicit statement that the optical regime condition is satisfied at all temperatures studied, while noting that the margin is narrowest at T ~ 0.5 GeV and results at the highest temperatures should be interpreted with this caveat. We thank the referee for this concrete and actionable suggestion. revision: yes

  3. Referee: Major Comment 3 (Sec. II.B, Eqs. 22-23): The quasi-static approximation T(t') ~ T(t) requires that the bath correlation time tau_E be much shorter than the hydrodynamic timescale tau_hydro. The authors state this hierarchy but do not provide estimates of tau_E. For the optical regime, tau_E ~ 1/T, which at T = 0.3 GeV gives tau_E ~ 0.66 fm, while tau_hydro ~ tau_0 = 0.6 fm — a ratio of only ~1, not a clean separation. The authors should either provide a more careful justification or acknowledge that the quasi-static approximation is marginal at the lower end of the temperature range and discuss how this affects the RHIC results in particular.

    Authors: The referee raises a legitimate concern that we cannot fully resolve within the present framework. The estimate tau_E ~ 1/T is a reasonable order-of-magnitude proxy for the bath correlation time, and at T = 0.3 GeV it indeed gives tau_E ~ 0.66 fm, comparable to tau_hydro = tau_0 = 0.6 fm. A more rigorous estimate would require computing tau_E from the spectral density J(omega) of the reservoir, which in our model is not independently specified — it enters only through the dissociation rate gamma = 2*pi*J(omega_0). Without an explicit form of J(omega) over a range of frequencies, we cannot extract a model-independent tau_E. We therefore agree that the quasi-static approximation is marginal at the lower end of the temperature range, and this directly affects the RHIC results (T_0 = 0.3 GeV), which are precisely the cases where coherence is claimed to persist beyond ~5 fm/c. We will revise the manuscript to: (1) include the explicit estimate tau_E ~ 1/T and the resulting ratio tau_E/tau_hydro at representative temperatures; (2) acknowledge that the separation is marginal at T ~ 0.3 GeV and improves at higher T; (3) add a caveat that the RHIC results, in particular, should be interpreted with this limitation in mind. We note that at LHC temperatures (T_0 = 0.5 GeV), tau_E ~ 0.4 fm and the ratio tau_E/tau_hydro ~ 0.67 is somewhat better, though still not a clean hierarchy. This is a genuine standing limitation of the quasi-static approach that we will be transparent about rather than attempting to dismiss. revision: yes

standing simulated objections not resolved
  • The quasi-static approximation (Major Comment 3) is genuinely marginal at the lower end of the temperature range (T ~ 0.3 GeV), and we cannot provide a more rigorous justification without an explicit spectral density J(omega) beyond the single point J(omega_0). The RHIC results are the most affected, and we can only flag this caveat rather than resolve it.
  • A fully self-consistent treatment of omega_0(T) (Major Comment 1) requires re-deriving the Lindblad master equation with a time-dependent system frequency, which goes beyond the current amplitude-coupling framework. We can discuss the sensitivity but cannot implement the full treatment in this manuscript.

Circularity Check

0 steps flagged

No significant circularity: decoherence timescales are genuine outputs of the Lindblad evolution with externally-provided inputs

full rationale

The paper's derivation chain is self-contained against external benchmarks. The two key input parameters are: (1) ω₀ = 0.457 GeV, fixed from the J/ψ RMS radius (PDG data, Ref. [52]) via Eq. (25), and (2) γ(T), the dissociation rate taken from the TAMU model (Refs. [1, 18] — external to the present authors). The Lindblad master equation (Eq. 14) is a standard textbook derivation in the optical regime (Born-Markov approximation, amplitude coupling). The central quantitative results — coherence lost within ~2 fm/c at LHC conditions and persisting beyond ~5 fm/c at RHIC conditions — are outputs of numerically solving the Lindblad equation with a time-dependent T(t) from Bjorken hydrodynamics (Eq. 44). These timescales are not fitted to data; they emerge from the evolution. The analytic decoherence rate Re(a_n) = (γ/2)n(1+2n̄) in Eq. (32) is derived as a large-n approximation of Eq. (18), not a fit. Figure 8, where the fitted momentum damping rate recovers the input γ, is explicitly framed by the authors as a consistency check of the Fokker-Planck correspondence ('we find to be equivalent to the input dissociation rate from the TAMU model'), not as a prediction. The viscous correction comparison uses an external Boltzmann transport code (Ref. [74], co-authored by S. Plumari) to generate the T(t) profile with η/s = 1/4π; this is a supporting input for one comparison case, not the central derivation. The self-citations to the Catania group's transport code and coalescence work (Refs [69-71, 73-79]) are not load-bearing for the decoherence results — they provide the viscous T(t) profile for a secondary comparison and context for future hadronization studies. No step in the derivation of the decoherence timescales reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

6 free parameters · 5 axioms · 0 invented entities

The paper does not invent new particles, forces, or entities. It uses a harmonic oscillator as a proxy for J/psi, which is a modeling simplification rather than an invented entity. All parameters are either physical inputs (masses, temperatures) or taken from external models (TAMU dissociation rates). The quasi-static Markov approximation is the main ad hoc element.

free parameters (6)
  • omega_0 (system frequency) = 0.457 GeV
    Fixed by matching the harmonic oscillator ground-state width to the J/psi RMS radius (sigma_r = 0.348 fm from Ref. [19]). Not fitted to decoherence data, but a model input chosen to represent the physical system.
  • gamma(T) (dissociation rate) = Temperature-dependent, from TAMU model [18]
    Phenomenologically inferred from quarkonium in-medium dissociation widths calculated from a non-perturbative many-body theory. Used as the Lindblad dissipation coefficient. Not fitted to the paper's own results.
  • alpha (coherent/cat state parameter) = 1 or 2
    Initial condition parameter for the density matrix, chosen to illustrate decoherence dynamics. Not a physical parameter fitted to data.
  • T_0 (initial temperature) = 0.3-0.5 GeV
    Physical input representing RHIC vs. LHC initial conditions, not a fitted parameter.
  • tau_0 (initial proper time) = 0.6 fm
    Set from Bjorken hydrodynamics convention. Standard physical input.
  • c_s^2 (speed of sound squared) = 1/3 (ideal) or from Boltzmann (viscous)
    Controls the cooling law T(t) = T_0 (tau_0/t)^(c_s^2). Standard hydrodynamic input.
axioms (5)
  • domain assumption Born-Markov approximation: the total density matrix factorizes as rho_t ~ rho(t) x rho_R and bath correlations decay fast (tau_E << tau_R).
    Standard OQS assumption invoked in Sec. II.A (Eq. 6) to derive the Lindblad master equation. Validity depends on the coupling strength and bath properties.
  • domain assumption Optical regime hierarchy: tau_S << tau_R, meaning the system frequency omega_0 >> damping rate gamma. Justified by E_B ~ omega_0 >= T for J/psi.
    Invoked in Sec. II.A (Eq. 17) to justify the rotating-wave approximation. The paper checks that E_B ~ omega_0 = 0.457 GeV >= T for the temperatures considered (0.2-0.5 GeV), but at T=0.5 GeV this is marginal.
  • ad hoc to paper Quasi-static approximation for time-dependent T(t): T(t') ~ T(t) during the bath correlation time, justified by tau_E << tau_hydro.
    Invoked in Sec. II.B (Eq. 22-23) to extend the Markov master equation to time-dependent temperature. The condition tau_E << tau_hydro is stated but not quantitatively verified for the QGP parameters used.
  • ad hoc to paper The system frequency omega_0 remains constant and unaffected by the environment throughout the evolution.
    Stated in Sec. VI: 'our assumption that the frequency omega_0 remains unaffected by the environment.' This ignores in-medium modifications of the binding energy, which are known to be significant near Tc.
  • domain assumption The QGP can be modeled as a static or 1D-expanding thermal bath of independent bosonic oscillators with amplitude coupling to the quarkonium.
    The reservoir model (Eq. 2-3) is a standard spin-boson type model. Its applicability to QGP degrees of freedom is a modeling assumption, not derived from QCD.

pith-pipeline@v1.1.0-glm · 25810 in / 3699 out tokens · 451222 ms · 2026-07-08T15:22:26.128777+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

83 extracted references · 83 canonical work pages · 46 internal anchors

  1. [1]

    Heavy Quarks in the Quark-Gluon Plasma

    R. Rapp and H. van Hees (2010), pp. 111–206, 0903.1096

  2. [2]

    Heavy Quarkonium in Extreme Conditions

    A. Rothkopf, Phys. Rept.858, 1 (2020), 1912.02253

  3. [3]

    Dong and V

    X. Dong and V. Greco, Prog. Part. Nucl. Phys.104, 97 (2019)

  4. [4]

    J. Zhao, K. Zhou, S. Chen, and P. Zhuang, Prog. Part. Nucl. Phys.114, 103801 (2020), 2005.08277

  5. [5]

    Matsui and H

    T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986)

  6. [6]

    Gonin et al

    M. Gonin et al. (NA50), Nucl. Phys. A610, 404C (1996)

  7. [7]

    J/psi production in Indium-Indium collisions at 158 GeV/nucleon

    R. Arnaldi et al. (NA60), Conf. Proc. C060726, 430 (2006), 0706.4361

  8. [8]

    J/psi Production vs Centrality, Transverse Momentum, and Rapidity in Au+Au Collisions at sqrt(s_NN) = 200 GeV

    A. Adare et al. (PHENIX), Phys. Rev. Lett.98, 232301 (2007), nucl-ex/0611020

  9. [9]

    B. I. Abelev et al. (STAR), Phys. Rev. C80, 041902 (2009), 0904.0439

  10. [10]

    Suppression of Upsilon Production in d+Au and Au+Au Collisions at sqrt(s_NN) = 200 GeV

    L. Adamczyk et al. (STAR), Phys. Lett. B735, 127 (2014), [Erratum: Phys.Lett.B 743, 537–541 (2015)], 1312.3675

  11. [11]

    J/$\psi$ suppression at forward rapidity in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV

    B. Abelev et al. (ALICE), Phys. Rev. Lett.109, 072301 (2012), 1202.1383

  12. [12]
  13. [13]

    Observation of sequential Upsilon suppression in PbPb collisions

    S. Chatrchyan et al. (CMS), Phys. Rev. Lett.109, 222301 (2012), [Erratum: Phys.Rev.Lett. 120, 199903 (2018)], 1208.2826

  14. [14]

    Heavy-flavour and quarkonium production in the LHC era: from proton-proton to heavy-ion collisions

    A. Andronic et al., Eur. Phys. J. C76, 107 (2016), 1506.03981

  15. [15]

    Real-time static potential in hot QCD

    M. Laine, O. Philipsen, P. Romatschke, and M. Tassler, JHEP03, 054 (2007), hep-ph/0611300

  16. [16]

    On temperature dependence of quarkonium correlators

    P. Petreczky, Eur. Phys. J. C62, 85 (2009), 0810.0258

  17. [17]

    Potential NRQCD: an effective theory for heavy quarkonium

    N. Brambilla, A. Pineda, J. Soto, and A. Vairo, Nucl. Phys. B566, 275 (2000), hep-ph/9907240

  18. [18]

    Charmonium Suppression and Regeneration from SPS to RHIC

    L. Grandchamp and R. Rapp, Nucl. Phys. A709, 415 (2002), hep-ph/0205305

  19. [19]

    T. Song, J. Aichelin, J. Zhao, P. B. Gossiaux, and E. Bratkovskaya, Phys. Rev. C108, 054908 (2023), 2305.10750

  20. [20]

    Stochastic potential and quantum decoherence of heavy quarkonium in the quark-gluon plasma

    Y. Akamatsu and A. Rothkopf, Phys. Rev. D85, 105011 (2012), 1110.1203

  21. [21]

    Real-time quantum dynamics of heavy quark systems at high temperature

    Y. Akamatsu, Phys. Rev. D87, 045016 (2013), 1209.5068. 16

  22. [22]

    Heavy quark master equations in the Lindblad form at high temperatures

    Y. Akamatsu, Phys. Rev. D91, 056002 (2015), 1403.5783

  23. [23]

    Quarkonium suppression in heavy-ion collisions: an open quantum system approach

    N. Brambilla, M. A. Escobedo, J. Soto, and A. Vairo, Phys. Rev. D96, 034021 (2017), 1612.07248

  24. [24]

    Heavy quarkonium suppression in a fireball

    N. Brambilla, M. A. Escobedo, J. Soto, and A. Vairo, Phys. Rev. D97, 074009 (2018), 1711.04515

  25. [25]

    Bottomonium suppression in an open quantum system using the quantum trajectories method

    N. Brambilla, M. A. Escobedo, M. Strickland, A. Vairo, P. Vander Griend, and J. H. Weber, JHEP05, 136 (2021), 2012.01240

  26. [26]

    X. Yao, W. Ke, Y. Xu, S. A. Bass, and B. M¨ uller, JHEP 01, 046 (2021), 2004.06746

  27. [27]

    Quarkonia Theory: From Open Quantum System to Classical Transport

    X. Yao, EPJ Web Conf.339, 01013 (2025), 2505.00098

  28. [28]

    Heavy quark bound states in a quark-gluon plasma: dissociation and recombination

    J.-P. Blaizot, D. De Boni, P. Faccioli, and G. Garberoglio, Nucl. Phys. A946, 49 (2016), 1503.03857

  29. [29]

    The Schr\"odinger-Langevin equation with and without thermal fluctuations

    R. Katz and P. B. Gossiaux, Annals Phys.368, 267 (2016), 1504.08087

  30. [30]

    Comparative Study of Quarkonium Transport in Hot QCD Matter

    A. Andronic et al., Eur. Phys. J. A60, 88 (2024), 2402.04366

  31. [31]

    Breuer and F

    H.-P. Breuer and F. Petruccione,The Theory of Open Quantum Systems(Oxford University Press, Oxford, 2002), ISBN 978-0-19-921390-0

  32. [32]

    J. Rais, H. van Hees, and C. Greiner, Phys. Rev. C111, 054918 (2025)

  33. [33]

    Open Quantum Systems with Kadanoff-Baym Equations

    T. Neidig, J. Rais, M. Bleicher, H. van Hees, and C. Greiner, Phys. Lett. B851, 138589 (2024), 2308.07659

  34. [34]

    W. A. De Jong, M. Metcalf, J. Mulligan, M. P losko´ n, F. Ringer, and X. Yao, Phys. Rev. D104, 051501 (2021), 2010.03571

  35. [35]

    Open quantum system approach to inclusive jet production in heavy-ion collisions

    Y. Mehtar-Tani, F. Ringer, B. Singh, and V. Vaidya, JHEP02, 048 (2026), 2504.00101

  36. [36]

    G. Coci, G. Parisi, S. Plumari, and M. Ruggieri (2025), 2507.04809

  37. [37]

    H. Iida, T. Kunihiro, A. Ohnishi, and T. T. Takahashi (2014), 1410.7309

  38. [38]

    Heavy quarkonia in a medium as a quantum dissipative system: Master-equation approach

    N. Borghini and C. Gombeaud, Eur. Phys. J. C72, 2000 (2012), 1109.4271

  39. [39]

    Quantum and Classical Dynamics of Heavy Quarks in a Quark-Gluon Plasma

    J.-P. Blaizot and M. A. Escobedo, JHEP06, 034 (2018), 1711.10812

  40. [40]

    Quarkonium dynamics in the quantum Brownian regime with non-abelian quantum master equations

    S. Delorme, R. Katz, T. Gousset, P. B. Gossiaux, and J.-P. Blaizot, JHEP06, 060 (2024), 2402.04488

  41. [41]

    Decoherence, the measurement problem, and interpretations of quantum mechanics

    M. Schlosshauer, Rev. Mod. Phys.76, 1267 (2004), quant-ph/0312059

  42. [42]

    Quantum Decoherence

    M. Schlosshauer, Phys. Rept.831, 1 (2019), 1911.06282

  43. [43]

    W. H. Zurek, Rev. Mod. Phys.75, 715 (2003), quant- ph/0105127

  44. [44]

    W. H. Zurek, Phys. Today44N10, 36 (1991), quant- ph/0306072

  45. [45]

    Gorini, A

    V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, J. Math. Phys.17, 821 (1976)

  46. [46]

    Lindblad, Commun

    G. Lindblad, Commun. Math. Phys.48, 119 (1976)

  47. [47]

    Brambilla, M

    N. Brambilla, M. A. Escobedo, A. Islam, M. Strickland, A. Tiwari, A. Vairo, and P. Vander Griend, JHEP08, 303 (2022), 2205.10289

  48. [48]

    Quantum Brownian motion of a heavy quark pair in the quark-gluon plasma

    T. Miura, Y. Akamatsu, M. Asakawa, and A. Rothkopf, Phys. Rev. D101, 034011 (2020), 1908.06293

  49. [49]

    A. O. Caldeira and A. J. Leggett, Physica A121, 587 (1983)

  50. [50]

    A. O. Caldeira and A. J. Leggett, Annals of physics149, 374 (1983)

  51. [51]

    H. J. Carmichael,An Open Systems Approach to Quan- tum Optics, vol. m18 ofLecture Notes in Physics (Springer, Berlin, Heidelberg, 1993), ISBN 978-3-540- 56689-1

  52. [52]

    Navas et al

    S. Navas et al. (Particle Data Group), Phys. Rev. D110, 030001 (2024)

  53. [53]

    Dynamical dissociation of quarkonia by wave function decoherence

    S. Kajimoto, Y. Akamatsu, M. Asakawa, and A. Rothkopf, Phys. Rev. D97, 014003 (2018), 1705.03365

  54. [54]

    Evolution of a quantum harmonic oscillator coupled to a minimal thermal environment

    A. Vidiella-Barranco, Physica A459, 78 (2016), 1605.01050

  55. [55]

    Weiss,Quantum Dissipative Systems(World Scien- tific, Singapore, 2008), ISBN 978-9-812-79179-5

    U. Weiss,Quantum Dissipative Systems(World Scien- tific, Singapore, 2008), ISBN 978-9-812-79179-5

  56. [56]

    C. A. Brasil, F. F. Fanchini, and R. d. J. Napoli- tano, Revista Brasileira de Ensino de F´ ısica35, 01–09 (2013), ISSN 1806-1117, URLhttp://dx.doi.org/10. 1590/S1806-11172013000100003

  57. [57]

    T. Song, J. Aichelin, J. Zhao, P. B. Gossiaux, and E. Bratkovskaya, EPJ Web Conf.316, 04008 (2025), 2409.19280

  58. [58]

    Thermal versus Direct $J/\Psi$ Production in Ultrarelativistic Heavy-Ion Collisions

    L. Grandchamp and R. Rapp, Phys. Lett. B523, 60 (2001), hep-ph/0103124

  59. [59]

    K. E. Cahill and R. J. Glauber, Phys. Rev.177, 1857 (1969)

  60. [60]

    D. F. Walls and G. J. Milburn (1985)

  61. [61]

    Magnus, Commun

    W. Magnus, Commun. Pure Appl. Math.7, 649 (1954)

  62. [62]

    M. O. Scully and M. S. Zubairy,Quantum Optics(Cam- bridge University Press, 1997)

  63. [63]

    Risken,The Fokker-Planck Equation: Methods of So- lution and Applications, vol

    H. Risken,The Fokker-Planck Equation: Methods of So- lution and Applications, vol. 18 ofSpringer Series in Syn- ergetics(Springer-Verlag, Berlin, Heidelberg, 1989), 2nd ed

  64. [64]

    E. P. Wigner, Phys. Rev.40, 749 (1932)

  65. [65]

    Hillery, R

    M. Hillery, R. F. O’Connell, M. O. Scully, and E. P. Wigner, Phys. Rept.106, 121 (1984)

  66. [66]

    W. B. Case, American Journal of Physics76, 937 (2008)

  67. [67]

    Yurke and D

    B. Yurke and D. Stoler, Phys. Rev. Lett.57, 13 (1986)

  68. [68]

    Saito and H

    Y. Saito and H. Hyuga, Journal of the Physical Society of Japan65, 1648 (1996)

  69. [69]

    Charmed Hadrons from Coalescence plus Fragmentation in relativistic nucleus-nucleus collisions at RHIC and LHC

    S. Plumari, V. Minissale, S. K. Das, G. Coci, and V. Greco, Eur. Phys. J. C78, 348 (2018), 1712.00730

  70. [70]

    Minissale, S

    V. Minissale, S. Plumari, Y. Sun, and V. Greco, Eur. Phys. J. C84, 228 (2024), 2305.03687

  71. [71]

    Minissale, V

    V. Minissale, V. Greco, and S. Plumari, Phys. Lett. B 860, 139190 (2025), 2405.19244

  72. [72]

    J. D. Bjorken, Phys. Rev. D27, 140 (1983)

  73. [73]

    Heavy-light flavour correlations of anisotropic flows at LHC energies within event-by-event transport approach

    S. Plumari, G. Coci, V. Minissale, S. K. Das, Y. Sun, and V. Greco, Phys. Lett. B805, 135460 (2020), 1912.09350

  74. [74]

    Nugara, V

    V. Nugara, V. Greco, and S. Plumari, Eur. Phys. J. C 85, 311 (2025), 2409.12123

  75. [75]

    Nugara, S

    V. Nugara, S. Plumari, L. Oliva, and V. Greco, Eur. Phys. J. C84, 861 (2024), 2311.11921

  76. [76]

    Nugara, N

    V. Nugara, N. Borghini, V. Greco, and S. Plumari, Phys. Lett. B872, 140122 (2026), 2509.05495

  77. [77]

    M. L. Sambataro, Y. Sun, V. Minissale, S. Plumari, and V. Greco, Eur. Phys. J. C82, 833 (2022), 2206.03160

  78. [78]

    M. L. Sambataro, V. Minissale, S. Plumari, and V. Greco, Phys. Lett. B849, 138480 (2024), 2304.02953

  79. [79]

    M. L. Sambataro, V. Minissale, S. Plumari, and V. Greco, Phys. Lett. B872, 140049 (2026), 2508.01024

  80. [80]

    Minissale, S

    V. Minissale, S. Plumari, and V. Greco, Phys. Lett. B 821, 136622 (2021), 2012.12001

Showing first 80 references.