REVIEW 3 major objections 5 minor 44 references
Probing vorticity and fluctuations in a rotating hadron resonance gas at LHC energy
T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read A rotating hadron resonance gas predicts an enhanced proton-to-pion ratio, and matching it to LHC data bounds the fireball vorticity at freeze-out between about 0.03 and 0.09 GeV.
desk verdict Useful rotating-HRG phenomenology with a genuinely new susceptibility study, but the headline vorticity extraction has a load-bearing internal inconsistency: the peripheral upper limit violates the paper's own causality bound R*omega <= 1. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the rotating ideal hadron resonance gas: a grand-canonical gas of all known hadrons up to 2.5 GeV inside a cylinder of radius $R$, with pressure written as a sum over orbital angular momentum modes $l$ and a squared Bessel function $J_\nu^2(k_r r)$ enforcing cylindrical symmetry. The load-bearing term is the rotation-induced effective chemical potential $\mu_{\rm eff}=(l+s_i)\omega$ in the dispersion relation $E_{l,i}=\sqrt{k_r^2+k_z^2+m_i^2}-(l+s_i)\omega$, which preferentially populates high-spin, high-mass states. The causality constraint $R\omega\leq1$ caps the allowed vorticity and, through the first zeros $\zeta_{l,1}$ of the Bessel functions, introduces an infrared cutoff $\Lambda^{\rm IR}_l=\zeta_{l,1}\omega$ on transverse momenta. This machinery converts a given $\omega$ into definite hadron yields, yield ratios, and susceptibility ratios.
What would settle it
Run a hadronic afterburner with baryon annihilation but no rotation on a hydrodynamic profile of Pb+Pb collisions at 5.02 TeV: if it reproduces the ALICE $p/\pi$ centrality trend, the reported $\omega$ values cannot be cleanly identified from this observable. Alternatively, measure $K^+/\pi^+$ in the same centrality bins; the model predicts a decrease with $\omega$, so a flat or increasing $K^+/\pi^+$ at the $p/\pi$-inferred $\omega$ values would refute the rotation interpretation.
Extended reading notes
Core claim
At fixed chemical freeze-out conditions ($T=158$ MeV, $\mu_B=0$), the rotating hadron resonance gas shifts each single-particle energy by $\mu_{\rm eff}=(l+s_i)\omega$, where $l$ is the orbital angular momentum quantum number and $s_i$ the spin. That shift lowers the effective energy of high-angular-momentum and high-spin states, producing a mass–spin ordering in the yield ratios: $p/\pi^+$, $\Delta^{++}/\pi^+$, $\Omega/\pi^+$, and $\rho^+/\pi^+$ increase with $\omega$, whereas $K^+/\pi^+$ decreases. Taking the centrality trend of the ALICE $p/\pi$ ratio as a proxy for growing vorticity from central to peripheral collisions, and assuming the $0$–$5\%$ bin has negligible rotation, the paper obtains an upper bound on freeze-out vorticity of about $0.028$ GeV ($\sim10^{22}$ s$^{-1}$) for central and $0.088$ GeV ($\sim10^{23}$ s$^{-1}$) for peripheral Pb+Pb collisions. It further reports that rotation enhances quadratic and mixed conserved-charge susceptibilities, leaves $\chi_3/\chi_2$ nearly unchanged, and generates an $\omega$-dependent split between $-3\chi^{BS}_{11}/\chi^S_2$ and $3\chi^{QS}_{11}/\chi^S_2$.
Load-bearing premise
The extraction rests on the assumption that the $0$–$5\%$ most central collisions have negligible vorticity and that the remaining centrality variation of the measured $p/\pi$ ratio is caused by rotation rather than by baryon annihilation or other centrality-dependent effects.
Editorial extensions
If this is right
- Thermal-model fits that omit rotation will systematically misread the proton yield: part of the $p/\pi$ enhancement now attributed to freeze-out parameters would be reassigned to vorticity, shifting the extracted $T$ and $\mu_B$.
- The same matching procedure applied to $\Omega/\pi$, $\Delta^{++}/\pi$, and $\rho^+/\pi$ should show a mass–spin ordering with centrality, a signature that distinguishes rotation from baryon annihilation.
- The predicted suppression of $K^+/\pi^+$ with $\omega$ provides a consistency check on the vorticity values extracted from $p/\pi$.
- Conserved-charge fluctuation observables such as $\chi^B_2$ and $\chi^{BS}_{11}$ increase sharply with rotation while $\chi_3/\chi_2$ stays flat, giving independent data-comparable probes of freeze-out vorticity.
- The reported upper bounds, $\sim0.028$ GeV central and $\sim0.088$ GeV peripheral, set the scale at which rotation must be considered in freeze-out physics at LHC energies.
Reading between the lines
- A decisive cross-check would be to repeat the extraction at several collision energies: since deposited angular momentum grows with $\sqrt{s_{NN}}$ and impact parameter, the inferred $\omega$ should grow accordingly if the interpretation is right, and should vanish for $p$+$p$ collisions where no vorticity is expected.
- The model's near-flat prediction for $K^+/\pi^+$ across centrality already constrains how large $\omega$ can be at freeze-out; a combined fit to both $p/\pi$ and $K^+/\pi$ would tighten the upper bounds.
- If the same vorticity drives both hadron abundances and global polarization, the centrality dependence of $\Lambda$ and $\bar\Lambda$ polarization from the STAR measurement could be compared point-by-point with the $\omega$ values extracted here.
- The extracted central value $\omega\sim0.028$ GeV sits close to the causality bound $R\omega\leq1$ for $R\approx6$ fm, so central collisions may already probe the model's physical ceiling; smaller systems or larger radii would clarify whether the bound or the data set the limit.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper extends the ideal hadron resonance gas model to a rigidly rotating system by modifying the single-particle energy with a rotational term -(l+s)ω (Eq. 2) and imposing a causality bound Rω≤1 (Section II). At fixed T=158 MeV and μ_B=0, the authors compute particle densities, hadron-to-pion ratios, and conserved-charge susceptibilities as functions of ω, finding strong enhancement for high-spin and high-mass hadrons and a moderate suppression of K+/π+. They then use the ALICE centrality dependence of (p+pbar)/(π+ + π-) in Pb+Pb at 5.02 TeV, normalized to the 0-5% bin, to map each centrality to a value of ω; this yields approximately 0.028 GeV for central and 0.088 GeV for peripheral collisions. The final part presents the ω dependence of χ2, χ3, χ4, mixed susceptibilities, and ratios such as χ3/χ2 and χ4/χ2, and argues these are rotation-sensitive observables.
Significance. If the central extraction were sound, the paper would provide a quantitative estimate of freeze-out vorticity from measured hadron yields, a useful input to the rotating-QCD program, and the fluctuation section offers falsifiable predictions for future measurements. The authors are transparent about the main assumptions: negligible vorticity in the 0-5% bin, rotation as one possible contributor to the p/π centrality trend, and survival of ω to chemical freeze-out. The model curves are internally consistent for ω within the stated causality bound, and the qualitative mass/spin ordering of the rotational enhancement is clear. However, the headline peripheral value lies outside the model's own validity domain, and the extraction is a mapping through a precomputed model curve rather than an independent test; these issues need to be fixed before the quantitative conclusion can be accepted.
major comments (3)
- [Sec. III B / Sec. II] The reported peripheral vorticity ω≈0.088 GeV violates the causality constraint Rω≤1 stated in Section II. The text assigns peripheral Pb+Pb collisions a maximum radius R=12.5 GeV^-1, giving Rω≈1.1>1. The peripheral point therefore lies in the region where the rotational term in Eq. (2) can exceed the free-particle energy component and where the IR-cutoff replacement in Eq. (4) is not a controlled approximation. The extracted peripheral value should be capped at 1/R=0.080 GeV, or the geometry used in the calculation should be changed so that every reported point satisfies the bound; either way, the abstract, the summary, and Fig. 5(b) must be revised.
- [Sec. III B] The extraction procedure computes p/π(ω) from the model in Fig. 4 and then inverts the measured ALICE p/π centrality ratios through that same curve. Consequently, the statement that rotation may explain the observed centrality dependence is a restatement of the inversion assumption rather than a model test. The extraction also assumes that the 0-5% bin has negligible vorticity and that the centrality variation of p/π is caused by rotation rather than by baryon annihilation or other final-state effects. The authors acknowledge the annihilation ambiguity, but the 'upper limit' wording in the abstract and Section IV should be qualified throughout, and the degeneracy should be quantified or at least discussed more concretely, for example by estimating the expected annihilation contribution.
- [Sec. II, Eq. (4) / Fig. 4] The model curves used for the extraction are computed at a fixed radial distance r=3 GeV^-1, while the centrality geometry is described by R=30 GeV^-1 for central and R=12.5 GeV^-1 for peripheral collisions. The paper does not clarify whether r or R is the radius that enters the causality bound and the IR cutoff in Eq. (4), nor why a single fixed r is appropriate for all centralities. If r=3 GeV^-1 is the actual size of the system in the calculation, the stated peripheral geometry is not being used; if R is the relevant radius, the extraction must be redone with the per-centrality boundary condition. This ambiguity directly affects the extracted ω values and should be resolved.
minor comments (5)
- [Abstract and Sec. IV] The values 0.028 and 0.088 GeV are described as 'upper limits', but under the assumptions listed in Sec. III B they are better characterized as conditional estimates; this wording should be adjusted for accuracy.
- [Eqs. (1), (3), and (4)] The definition of the IR cutoff Λ_IR_l = ζ_{l,1}ω and the modified integration in Eq. (4) should specify the lower limit of the k_r integral explicitly and explain how the l=0 mode is treated, since the current notation is not sufficient for reproducing the numerical results.
- [Fig. 5 and Sec. III B] The experimental quantity is (p+pbar)/(π+ + π-) from Ref. [37], but the text frequently writes p/π; please use consistent notation and state explicitly whether the calculated p/π(ω) curve used for the extraction includes feed-down contributions in the same way as the experimental ratio.
- [Fig. 5 caption] The uncertainty propagation is described only as 'uncertainties associated with the p/π ratio are propagated'; please specify whether the shaded band includes uncorrelated uncertainties only and how the correlation between the numerator and denominator in the ratio was handled.
- [References and text] There are several typographical issues: 'V orticity' in Ref. [5], 'Barun-Mumnzinger' in Ref. [31] should be 'Braun-Munzinger', 'Rev. Part. Phys. C38' in Ref. [35] should be corrected, and 'asses' in Sec. III C should be 'assess'.
Circularity Check
The p/pi-to-vorticity extraction is a calibration, and the claim that rotation 'may explain' the p/pi centrality dependence is a restatement of that calibration; the fluctuation and K/pi results retain independent content.
-
fitted input called prediction
[Section III B, Figure 5]
"For each centrality class, the p/π ratios, normalized to the corresponding values of the 0−5% centrality interval, is mapped with the one given by nonzero vorticity values shown in Fig 4. This correspondence enables the extraction of the vorticity strength at a given centrality."
The p/pi ratio is first computed as a function of omega (Fig. 4). The measured centrality-dependent p/pi ratios are then inverted through this same function to read off omega. Therefore the later statement, in Section III B, that rotation 'may significantly influence the hadron yields and potentially explain the observed centrality dependence of the p/π yield ratio' is not a test of the model: the data were used to set omega, so the agreement is enforced by construction. The only independent inputs are the assumed zero-vorticity 0-5% baseline and the model's functional shape; the extracted omega values inherit the data and provide no additional confirmation of the rotation mechanism for p/pi itself. The separate K/pi and susceptibility calculations are not forced by this step.
full rationale
The paper's central quantitative result is an extraction, not a prediction: omega is obtained by equating the model's p/pi ratio to the ALICE centrality data. The abstract and summary honestly call this an 'estimate' or 'upper limit', and the text explicitly assumes 0-5% centrality has negligible vorticity and lists baryon annihilation as a competing mechanism. Thus the circularity is limited to the rhetorical step of saying the fitted curve 'explains' the data. No self-citation is load-bearing: the rotating-HRG formulas are also attributed to non-overlapping references, and the susceptibility and K/pi sections are self-contained model calculations. A separate internal-consistency concern exists: if the peripheral radius R=12.5 GeV^-1 were used, the extracted omega=0.088 GeV would exceed the stated causality bound R*omega<=1 (giving 1.1>1); however, the paper fixes r=3 GeV^-1 for the calculations, making R*omega=0.264<1. That ambiguity is a correctness risk, not a circularity. Overall, the paper is mostly a calibration study with one overreaching explanatory claim, so a moderate score is appropriate.
Assumptions & free parameters
free parameters (2)
- Vorticity omega per centrality =
0.028 to 0.088 GeV
- Radial distance r in the Bessel integrand =
3.0 GeV^-1
assumptions (7)
- domain assumption Ideal gas of non-interacting hadrons and resonances from the PDG list up to 2.5 GeV.
- domain assumption Rotation modifies the single-particle energy as E = sqrt(k_r^2+k_z^2+m_i^2) - (l+s_i)omega.
- domain assumption Causality bound R*omega <= 1 and the Bessel infrared cutoff Lambda_IR = zeta_{l,1} omega.
- domain assumption Chemical freeze-out at T=158 MeV with mu_B=mu_Q=mu_S=0 for Pb-Pb at 5.02 TeV.
- ad hoc to paper The 0-5% central collision bin has negligible vorticity.
- ad hoc to paper The measured p/pi centrality dependence is attributed to rotation rather than to baryon annihilation or other centrality effects.
- ad hoc to paper Vorticity survives undiminished until chemical freeze-out.
Cite this review
Pith. "Pith review of Probing vorticity and fluctuations in a rotating hadron resonance gas at LHC energy." pith.science (2026). https://pith.science/paper/MPFQEP67
@misc{pith2026260810561,
author = {Pith},
title = {Pith review of: Probing vorticity and fluctuations in a rotating hadron resonance gas at LHC energy},
year = {2026},
howpublished = {\url{https://pith.science/paper/MPFQEP67}},
note = {Machine review of arXiv:2608.10561}
}
abstract
A large vorticity produced in non-central ultra-relativistic heavy-ion collisions induces an effective chemical potential in both partonic and hadronic matter, thereby influencing the quark-hadron transition and its associated properties. In this work, we investigate the influence of rotation on hadron yields within the framework of Hadron Resonance Gas (HRG) model. Our results show that vorticity significantly modifies hadron yields and their ratios. Most notably, rotation enhances the $p/\pi^+$ ratio while suppressing the $K^+/\pi^+$ ratio, suggesting that these observables may serve as sensitive probes of the rotational properties of the medium created in heavy-ion collisions. To quantify the effect of rotation, the calculated dependence of the $p/\pi^+$ ratio on vorticity is compared with the centrality dependence of the $p/\pi$ ratio measured by the ALICE collaboration in Pb+Pb collisions at $\sqrt{s_{_{NN}}}$ = 5.02 TeV. From this comparison, we estimate the maximum vorticity produced at freeze-out for different collision centralities. In case of peripheral collisions, the freeze-out vorticity ($\omega$) is found to reach an upper bound value of approximately 0.088 GeV, whereas for central collisions it is about 0.028 GeV. Furthermore, we investigate the effect of rotation on fluctuations of conserved quantities and their correlations. This study provides a quantitative framework for assessing the role of rotation in the thermodynamics of hadronic matter and its phenomenological consequences for heavy-ion collisions.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
- [1]
-
[2]
D. H. Rischke, Prog. Part. Nucl. Phys.52, 197 (2004)
work page 2004
- [3]
- [4]
-
[5]
V orticity and spin polarization in heavy ion collisions: Transport models,
X.-G. Huang, J. Liao, Q. Wang, and X.-L. Xia, “V orticity and spin polarization in heavy ion collisions: Transport models,” in Strongly Interacting Matter under Rotation, edited by F. Becat- tini, J. Liao, and M. Lisa (Springer International Publishing, Cham, 2021) pp. 281–308
work page 2021
-
[6]
W. T. Deng and X. G. Huang, Phys. Rev. C93, 064907 (2016)
work page 2016
-
[7]
Z. T. Liang and X. N. Wang, Phys. Rev. Lett.94, 102301 (2016)
work page 2016
- [8]
Show all 44 references
-
[9]
Adamczyket al.(STAR), Nature548, 62 (2017)
L. Adamczyket al.(STAR), Nature548, 62 (2017)
2017
-
[10]
Fujimoto, K
Y . Fujimoto, K. Fukushima, and Y . Hidaka, Phys. Lett. B816, 136184 (2021)
2021
-
[11]
Huang, Rept
X.-G. Huang, Rept. Prog. Phys.79, 076302 (2016)
2016
-
[12]
G. S. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. D. Katz, S. Krieg, A. Schafer, and K. K. Szabo, JHEP02, 044 (2012)
2012
-
[13]
B. Fu, S. Y . F. Liu, L. Pang, H. Song, and Y . Yin, Phys. Rev. Lett.127, 142301 (2021)
2021
-
[14]
Fukushima, D
K. Fukushima, D. E. Kharzeev, and H. J. Warringa, Phys. Rev. D78, 074033 (2008)
2008
-
[15]
Mukherjee, D
G. Mukherjee, D. Dutta, and D. K. Mishra, Eur. Phys. J. C84, 258 (2024)
2024
-
[16]
V ovchenko, Phys
V . V ovchenko, Phys. Rev. C110, 034914 (2024)
2024
-
[17]
Padhan, K
N. Padhan, K. K. Pradhan, A. Chatterjee, and R. Sahoo, Phys. Lett. B879, 140604 (2026)
2026
-
[18]
Mukherjee, D
G. Mukherjee, D. Dutta, and D. K. Mishra, Phys. Lett. B846, 138228 (2023)
2023
-
[19]
Ebihara, K
S. Ebihara, K. Fukushima, and K. Mameda, Phys. Lett. B764, 94 (2017)
2017
-
[20]
D. T. Son and P. Surowka, Phys. Rev. Lett.103, 191601 (2009)
2009
-
[21]
Yamamoto and Y
A. Yamamoto and Y . Hirono, Phys. Rev. Lett.111, 081601 (2013)
2013
-
[22]
H.-L. Chen, K. Fukushima, X.-G. Huang, and K. Mameda, Phys. Rev. D93, 104052 (2016)
2016
-
[23]
M. N. Chernodub and S. Gongyo, J. High. Energy. Phys.1701, 136 (2017)
2017
-
[24]
X. Wang, M. Wei, Z. Li, and M. Huang, Phys. Rev. D99, 016018 (2019)
2019
-
[25]
Cleymans, H
J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, Phys. Rev. C73, 034905 (2006)
2006
-
[26]
Bazavovet al., Phys
A. Bazavovet al., Phys. Rev. Lett.109, 192302 (2012)
2012
-
[27]
Marczenko, M
M. Marczenko, M. Szyma ´nski, P. M. Lo, B. Karmakar, P. Huovinen, C. Sasaki, and K. Redlich, Phys. Rev. C110, 065203 (2024)
2024
-
[28]
Andronic, P
A. Andronic, P. Braun-Munzinger, J. Stachel, and M. Winn, Phys. Lett. B718, 80 (2012)
2012
-
[29]
Karsch and K
F. Karsch and K. Redlich, Phys. Lett. B695, 136 (2011)
2011
-
[30]
P. Garg, D. K. Mishra, P. K. Netrakanti, B. Mohanty, A. K. Mohanty, and B. K. Singh, Phys. Lett. B726, 691 (2013)
2013
-
[31]
Barun-Mumnzinger, J
P. Barun-Mumnzinger, J. Stachel, J. P. Wessels, and N. Xu, Phys. Lett. B344, 43 (1995)
1995
-
[32]
Gupta, D
S. Gupta, D. Mallick, D. K. Mishra, B. Mohanty, and N. Xu, Phys. Lett. B829, 137021 (2022)
2022
-
[33]
M. N. Chernodub, Phys. Rev. D103, 054027 (2021)
2021
-
[34]
Adamczyket al.(STAR), Phys
L. Adamczyket al.(STAR), Phys. Rev. C96, 044904 (2017)
2017
-
[35]
K. A. Oliveet al., Rev. Part. Phys. C38, 090001 (2014)
2014
-
[36]
Andronic, P
A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, Nature561, 321 (2018)
2018
-
[37]
Acharyaet al.(ALICE), Phys
S. Acharyaet al.(ALICE), Phys. Rev. C101, 044907 (2020)
2020
-
[38]
Ejiri, F
S. Ejiri, F. Karsch, and K. Redlich, Phys. Lett. B633, 275 (2006)
2006
-
[39]
Ding, J.-B
H.-T. Ding, J.-B. Gu, A. Kumar, and S.-T. Li, Phys. Rev. D 112, 094508 (2025)
2025
-
[40]
V . Koch, A. Majumder, and J. Randrup, Phys. Rev. Lett.95, 182301 (2005)
2005
-
[41]
Jeon and V
S. Jeon and V . Koch, Phys. Rev. Lett.85, 2076 (2000)
2000
-
[42]
Majumder and B
A. Majumder and B. Muller, Phys. Rev. C74, 054901 (2006)
2006
-
[43]
Bazavovet al.(HotQCD), Phys
A. Bazavovet al.(HotQCD), Phys. Rev. D86, 034509 (2012)
2012
-
[44]
Astrakhantsev, V
N. Astrakhantsev, V . V . Braguta, A. Y . Kotov, and A. A. Roenko, Phys. Rev. D109, 094511 (2024)
2024
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.