Recognition: 2 theorem links
· Lean TheoremQuarkonium p_{rm T} spectra in heavy--ion collisions at LHC energies within a hydrodynamic core--corona framework
Pith reviewed 2026-05-11 01:45 UTC · model grok-4.3
The pith
An analytical hydrodynamic model with core-corona separation reproduces quarkonium transverse momentum spectra in lead-lead collisions at LHC energies.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Within the hydrodynamic core-corona framework, the transverse momentum spectra of charmonium (J/ψ, ψ(2S)) and bottomonium (Υ(nS)) states are calculated by integrating the thermal distribution over a constant-temperature freeze-out hypersurface in a cylindrically symmetric, boost-invariant medium with Hubble-like transverse flow; the core-corona separation incorporates thermal production in the dense interior and non-thermal contributions from the corona, yielding spectra and ratios that match ALICE and CMS data across a wide pT range while indicating additional hard processes at high pT for J/ψ and sequential suppression for bottomonia.
What carries the argument
Analytical relativistic hydrodynamics under cylindrical symmetry with boost-invariant longitudinal expansion and Hubble-like transverse flow, evaluated on a constant-temperature freeze-out hypersurface via the Cooper-Frye formalism and combined with core-corona separation of thermal and non-thermal quarkonium contributions.
If this is right
- The model accounts for high-pT deviations in J/ψ spectra by invoking additional hard production mechanisms outside the thermal core.
- Bottomonium spectra and their yield ratios follow the expected sequential suppression ordered by binding energy.
- The ψ(2S)/J/ψ ratio is reproduced across the measured pT range, confirming the role of the core-corona distinction.
- The same framework applies uniformly to both charmonium and bottomonium sectors at LHC energies.
Where Pith is reading between the lines
- If the analytical flow profile suffices here, similar closed-form hydrodynamics could be tested on other heavy-flavor observables without full numerical codes.
- The success at 5.02 TeV suggests the assumed symmetry captures the dominant expansion features for quarkonium decoupling, but deviations at lower energies would test the Hubble-like flow assumption.
- Extending the constant-temperature freeze-out to a temperature-dependent surface could quantify how sensitive the spectra are to the precise decoupling condition.
Load-bearing premise
The expanding medium is taken to obey cylindrical symmetry, boost-invariant longitudinal flow, and Hubble-like transverse flow, with all quarkonia decoupling at one fixed temperature on the freeze-out surface.
What would settle it
A new measurement of the pT spectra or ψ(2S)/J/ψ ratio in a different centrality bin or at a lower collision energy that shows large deviations from the model's predictions while the assumed flow symmetry remains intact would falsify the claim.
read the original abstract
We present a systematic study of the transverse momentum ($p_{\rm T}$) spectra of charmonium (J/$\psi$, $\psi(2S)$) and bottomonium ($\Upsilon(nS)$) states in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV within an analytical relativistic hydrodynamics framework. The medium evolution is described assuming cylindrical symmetry with boost-invariant longitudinal expansion and Hubble-like transverse flow. Quarkonium spectra are evaluated using the Cooper-Frye formalism on a constant-temperature freeze-out hypersurface, supplemented by a core--corona approach to include both thermal and non-thermal contributions. The model describes the measurements from ALICE and CMS over a broad $p_{\rm T}$ range. For charmonium, both the spectra and the $\psi(2S)$/J/$\psi$ ratio are well reproduced, while deviations at high $p_{\rm T}$ for J/$\psi$ indicate additional hard production mechanisms. In the bottomonium sector, the $\Upsilon(nS)$ spectra and their yield ratios are successfully described, consistent with the expected sequential suppression pattern. These results demonstrate that an analytical hydrodynamic approach combined with a core-corona framework provides a unified and transparent description of quarkonium production in heavy--ion collisions at LHC energies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a systematic study of the transverse momentum spectra of charmonium (J/ψ, ψ(2S)) and bottomonium (Υ(nS)) states in Pb-Pb collisions at √s_NN = 5.02 TeV. It employs an analytical relativistic hydrodynamics framework assuming cylindrical symmetry, boost-invariant longitudinal expansion, and Hubble-like transverse flow. Spectra are evaluated with the Cooper-Frye formalism on a constant-temperature freeze-out hypersurface, combined with a core-corona approach to include thermal and non-thermal contributions. The authors report that the model describes ALICE and CMS data for spectra and yield ratios over a broad p_T range, with deviations at high p_T for J/ψ and consistency with sequential suppression for bottomonia.
Significance. If the central results hold, the work offers a transparent, computationally simple framework that unifies thermal and hard production mechanisms for quarkonia, potentially useful as a benchmark or for quick estimates in heavy-ion phenomenology. The analytical hydro plus core-corona split makes the separation of contributions explicit. However, the significance is reduced by the reliance on three adjustable parameters and the idealized flow assumptions, which may limit the extent to which the agreement constitutes a genuine test of the hydrodynamic picture rather than a fit.
major comments (2)
- The hydrodynamic evolution assumes cylindrical symmetry with boost-invariant longitudinal expansion and Hubble-like transverse flow (v_r ∝ r), with spectra obtained via Cooper-Frye on a single constant-T hypersurface. This rigidly fixes the radial boost and effective temperature. For states whose mass is comparable to the flow velocity scale, mismatches with realistic viscous evolution, longitudinal gradients, or fluctuating initial conditions would directly affect the low-to-intermediate p_T shape once the core-corona split is applied. The manuscript should quantify the sensitivity of the spectra to the flow profile (e.g., by varying the Hubble parameter or comparing to numerical hydro) to establish that the thermal component is physically motivated rather than absorbed into the corona weights.
- The abstract states that the model describes the data for most observables while noting deviations at high p_T for J/ψ. Because the Hubble flow parameter, freeze-out temperature, and core-corona separation weights are free parameters typically tuned to spectra, the reported agreement may largely reflect the fit rather than independent predictive power. The paper should report quantitative fit metrics (e.g., χ² per degree of freedom) and show results with the corona component turned off to demonstrate that the analytical hydro component supplies the correct shape on its own.
minor comments (1)
- The abstract mentions deviations at high p_T for J/ψ but does not quantify the p_T threshold or the size of the discrepancy; adding this information would clarify the range of validity.
Simulated Author's Rebuttal
We thank the referee for the constructive comments and the positive overall assessment of our work. We address each major comment point by point below, indicating the revisions we will implement in the manuscript.
read point-by-point responses
-
Referee: The hydrodynamic evolution assumes cylindrical symmetry with boost-invariant longitudinal expansion and Hubble-like transverse flow (v_r ∝ r), with spectra obtained via Cooper-Frye on a single constant-T hypersurface. This rigidly fixes the radial boost and effective temperature. For states whose mass is comparable to the flow velocity scale, mismatches with realistic viscous evolution, longitudinal gradients, or fluctuating initial conditions would directly affect the low-to-intermediate p_T shape once the core-corona split is applied. The manuscript should quantify the sensitivity of the spectra to the flow profile (e.g., by varying the Hubble parameter or comparing to numerical hydro) to establish that the thermal component is physically motivated rather than absorbed into the corona weights.
Authors: We agree that a quantitative assessment of sensitivity to the flow profile strengthens the physical interpretation. In the revised manuscript we will add a dedicated subsection (or appendix) presenting results obtained by varying the Hubble parameter by ±20% around its default value while readjusting only the overall normalization. These calculations show that the low-to-intermediate p_T spectral shapes remain robust, with the main variations absorbed into the core-corona weights; this supports that the thermal component is not an artifact of parameter tuning. A full comparison against numerical viscous hydrodynamics lies outside the scope of the present analytical study, whose purpose is to supply a transparent benchmark; we will nevertheless add a brief discussion of this limitation and its implications for the model assumptions. revision: partial
-
Referee: The abstract states that the model describes the data for most observables while noting deviations at high p_T for J/ψ. Because the Hubble flow parameter, freeze-out temperature, and core-corona separation weights are free parameters typically tuned to spectra, the reported agreement may largely reflect the fit rather than independent predictive power. The paper should report quantitative fit metrics (e.g., χ² per degree of freedom) and show results with the corona component turned off to demonstrate that the analytical hydro component supplies the correct shape on its own.
Authors: We accept that quantitative goodness-of-fit measures and an explicit separation of components would improve clarity. In the revised version we will report χ² per degree of freedom for the p_T spectra fits to both ALICE and CMS data. We have also prepared additional plots (to be included as a new figure or in the main text) that isolate the pure thermal (core) contribution. These show that the hydrodynamic component alone reproduces the low-p_T rise and the characteristic flow-induced hardening, while the corona term is required only for the high-p_T tail; this explicitly demonstrates the independent role of the analytical hydro part. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper presents an analytical hydrodynamic model (cylindrical symmetry, boost-invariant expansion, Hubble-like transverse flow) combined with Cooper-Frye freeze-out and a core-corona split to compute quarkonium pT spectra. The provided text states that the model 'describes the measurements' over a broad pT range and reproduces ratios, but contains no quoted equations or sections showing that any central result (e.g., spectra shapes or ratios) is obtained by fitting a parameter to the very same data and then relabeling the output as a prediction. No self-citation chains are invoked as load-bearing uniqueness theorems, and the derivation chain rests on explicit physical assumptions rather than reducing to its inputs by construction. This is the normal case of a model with adjustable parameters being compared to data.
Axiom & Free-Parameter Ledger
free parameters (3)
- Hubble-like transverse flow parameter
- Freeze-out temperature
- Core-corona separation parameters
axioms (2)
- domain assumption Cylindrical symmetry with boost-invariant longitudinal expansion
- domain assumption Cooper-Frye formalism on a constant-temperature freeze-out hypersurface
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The medium evolution is described assuming cylindrical symmetry with boost-invariant longitudinal expansion and Hubble-like transverse flow. Quarkonium spectra are evaluated using the Cooper-Frye formalism on a constant-temperature freeze-out hypersurface, supplemented by a core-corona approach
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel contradicts?
contradictsCONTRADICTS: the theorem conflicts with this paper passage, or marks a claim that would need revision before publication.
The model parameters are determined via a χ² minimization procedure... T≈160 MeV and ⟨v_T⟩≈0.60... T≈224 MeV and ⟨v_T⟩≈0.56... with three parameters total
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]
- [2]
-
[3]
H. Satz, J. Phys. G32, R25 (2006)
work page 2006
-
[4]
L. Kluberg and H. Satz, Relativistic Heavy Ion Physics V ol.23(Springer-Verlag, Berlin, Heidelberg, 2010)
work page 2010
-
[5]
R. Rapp, D. Blaschke, and P. Crochet, Prog. Part. Nucl. Phys.65, 209 (2010)
work page 2010
-
[6]
J. Zhao, K. Zhou, S. Chen, and P. Zhuang, Prog. Part. Nucl. Phys.114, 103801 (2020)
work page 2020
- [7]
-
[8]
Rothkopf, Heavy quarkonium in extreme conditions, Phys
A. Rothkopf, Heavy quarkonium in extreme conditions, Phys. Rep.858, 1 (2020)
work page 2020
- [9]
-
[10]
P. P. Bhaduri, P. Hegde, H. Satz, and P. Tribedy, Lecture Notes in Physics V ol.785(Springer, Berlin, 2010), p. 179
work page 2010
-
[11]
M. C. Abreuet al.(NA50 Collaboration), Phys. Lett. B 477, 28 (2000)
work page 2000
-
[12]
Arnaldiet al.(NA60 Collaboration), Phys
R. Arnaldiet al.(NA60 Collaboration), Phys. Rev. Lett 99, 132032 (2007)
work page 2007
-
[13]
Adareet al.(PHENIX Collaboration), Phys
A. Adareet al.(PHENIX Collaboration), Phys. Rev. Lett.98, 232301 (2007)
work page 2007
-
[14]
Adareet al.(PHENIX Collaboration), Phys
A. Adareet al.(PHENIX Collaboration), Phys. Rev. C 84, 054912 (2011)
work page 2011
-
[15]
Adamet al.(STAR Collaboration), Phys
J. Adamet al.(STAR Collaboration), Phys. Lett. B797, 134917 (2019)
work page 2019
-
[16]
B. B. Abelevet al.(ALICE Collaboration), Phys. Lett. B734, 314 (2014)
work page 2014
-
[17]
Adamet al.(ALICE Collaboration), J
J. Adamet al.(ALICE Collaboration), J. High Energy Phys.05(2016) 179
work page 2016
-
[18]
Adamet al.(ALICE Collaboration), Phys
J. Adamet al.(ALICE Collaboration), Phys. Lett. B 766, 212 (2017)
work page 2017
-
[19]
Acharyaet al.(ALICE Collaboration), J
S. Acharyaet al.(ALICE Collaboration), J. High En- ergy Phys.02, 41 (2020)
work page 2020
-
[20]
Acharyaet al.(ALICE Collaboration), Phys
S. Acharyaet al.(ALICE Collaboration), Phys. Lett. B 849, 138451 (2024)
work page 2024
-
[21]
Acharyaet al.(ALICE Collaboration), Phys
S. Acharyaet al.(ALICE Collaboration), Phys. Rev. Lett.132, 042301 (2024)
work page 2024
-
[22]
Chatrchyanet al.(CMS Collaboration), J
S. Chatrchyanet al.(CMS Collaboration), J. High En- ergy Phys.05, 063 (2012)
work page 2012
-
[23]
A. M. Sirunyanet al.(CMS Collaboration), Phys. Rev. Lett.118, 162301 (2017)
work page 2017
-
[24]
A. M. Sirunyanet al.(CMS Collaboration), Eur. Phys. J. C78, 509 (2018)
work page 2018
-
[25]
Aaboudet al.(ATLAS Collaboration), Eur
M. Aaboudet al.(ATLAS Collaboration), Eur. Phys. J. C78, 762 (2018)
work page 2018
-
[26]
Aaijet al.(LHCb Collaboration), Phys
R. Aaijet al.(LHCb Collaboration), Phys. Rev. C105, L032201 (2022)
work page 2022
-
[27]
Acharyaet al.(ALICE Collaboration), Phys
S. Acharyaet al.(ALICE Collaboration), Phys. Lett. B 822, 136579 (2021)
work page 2021
-
[28]
Tumasyanet al.(CMS Collaboration), Phys
A. Tumasyanet al.(CMS Collaboration), Phys. Lett. B 790, 270–293 (2019)
work page 2019
-
[29]
Tumasyanet al.(CMS Collaboration), Phys
A. Tumasyanet al.(CMS Collaboration), Phys. Rev. Lett.120, 142301 (2018)
work page 2018
-
[30]
Tumasyanet al.(CMS Collaboration), Phys
A. Tumasyanet al.(CMS Collaboration), Phys. Rev. Lett.133, 022302 (2024)
work page 2024
-
[31]
Aadet al.(ATLAS Collaboration), Phys
G. Aadet al.(ATLAS Collaboration), Phys. Rev. C107, 054912 (2023)
work page 2023
-
[32]
X. Du, M. He, and R. Rapp, Phys. Rev. C96, 054901 (2017). 7
work page 2017
- [33]
-
[34]
K. A. Bugaev, M. Gazdzicki, and M. I. Gorenstein, Phys. Lett. B544, 127 (2002)
work page 2002
-
[35]
A. Andronic, P. Braun-Munzinger, M. K. Köhler, K. Redlich, and J. Stachel, Phys. Lett. B797, 134836 (2019)
work page 2019
-
[36]
Khachatryanet al.(CMS Collaboration), Phys
V . Khachatryanet al.(CMS Collaboration), Phys. Lett. B770, 357 (2017)
work page 2017
-
[37]
K. Reygers, A. Schmah, A. Berdnikova, and X. Sun, Phys. Rev. C101, 064905 (2020)
work page 2020
-
[38]
Acharyaet al., (ALICE Collaboration), Phys
S. Acharyaet al., (ALICE Collaboration), Phys. Rev. Lett.123, 192301 (2019)
work page 2019
-
[39]
A. M. Sirunyanet al.(CMS Collaboration), Phys. Lett. B819, 136385 (2021)
work page 2021
-
[40]
Acharyaet al., (ALICE Collaboration), Phys
S. Acharyaet al., (ALICE Collaboration), Phys. Rev. Lett.119, 242301 (2017)
work page 2017
-
[41]
Acharyaet al., (ALICE Collaboration), J
S. Acharyaet al., (ALICE Collaboration), J. High En- ergy Phys.02, 12 (2019)
work page 2019
-
[42]
Acharyaet al., (ALICE Collaboration), J
S. Acharyaet al., (ALICE Collaboration), J. High En- ergy Phys.10, 141 (2020)
work page 2020
-
[43]
M. He, B. Wu, and R. Rapp, Phys. Rev. Lett.128, 162301 (2022)
work page 2022
-
[44]
P. P. Bhaduri, N. Borghini, A. Jaiswal, and M. Strick- land, Phys. Rev. C100, 051901(R) (2019)
work page 2019
-
[45]
P. P. Bhaduri, M. Alqahtani, N. Borghini, A. Jaiswal, and M. Strickland, Eur. Phys. J. C81, 585 (2021)
work page 2021
-
[46]
A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, Nucl. Phys. A789, 334 (2007)
work page 2007
-
[47]
A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, Phys. Lett. B652, 259 (2007)
work page 2007
-
[48]
A. Andronic, P. Braun-Munzinger, M. K. Köhler, and J. Stachel, Nucl. Phys. A982, 759 (2019)
work page 2019
-
[49]
A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, Nature (London)561, 321 (2018)
work page 2018
- [50]
-
[51]
S. Chatterjee, R. M. Godbole, and S. Gupta, Phys. Lett. B727, 554 (2013)
work page 2013
- [52]
- [53]
-
[54]
S. Lin, J. Liao, Nucl. Phys. A 837, 195-209 (2010)
work page 2010
-
[55]
M. S. Ali, D. Biswas, A. Jaiswal, S. K. Singh, Eur. Phys. J. C85, 30 (2025)
work page 2025
- [56]
-
[57]
Acharyaet al.(ALICE Collaboration), J
S. Acharyaet al.(ALICE Collaboration), J. High En- ergy Phys.10, 084 (2019)
work page 2019
-
[58]
Acharyaet al.(ALICE Collaboration), Eur
S. Acharyaet al.(ALICE Collaboration), Eur. Phys. J. C83, 61 (2023)
work page 2023
- [59]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.