REVIEW 3 major objections 3 minor 44 references
Determining the Duration of the Hadronic Stage at RHIC-BES Energies via Resonance Suppression Using a Full Set of Rate Equations
T0 review · 3 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The hadronic stage in heavy-ion collisions lasts 2-4 times longer than STAR's estimate because resonance regeneration makes the K*/K ratio fall linearly rather than exponentially.
desk verdict Regeneration makes K*/K decay roughly linear and implies a 2-4x longer hadronic stage, but Eq. (12)'s universal slope is asserted, not demonstrated, and the advertised collision-loss term is absent from the solved network. 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 machinery is a network of 23 coupled rate equations for hadron species in partial chemical equilibrium, where each resonance's multiplicity evolves under detailed balance between its decay $h^* \to h + X$ and its regeneration $h + X \to h^*$. Equilibrium multiplicities are fixed at chemical freeze-out using a standard parametrization of $T_{\rm ch}$ and $\mu_{B,\rm ch}$, and the fireball expansion enters through a parametrized volume $V(t)$ with conserved entropy, baryon number, and strangeness. The load-bearing output is the nearly energy-independent slope of the $K^*/K$ ratio versus time, which converts a measured ratio difference into a hadron phase lifetime through Eq. (12).
What would settle it
Measure the $K^*/K$ ratio in a centrality or system where the hadron phase duration is independently known (e.g., from two-particle HBT radii) and check whether Eq. (12) with the universal slope reproduces that duration; a systematic mismatch beyond the stated uncertainties would show that the slope is not universal.
Extended reading notes
Core claim
The central finding is that the time evolution of the $K^*/K$ ratio after chemical freeze-out is controlled by a counterbalance between resonance decay and regeneration, and this balance makes the ratio decrease approximately linearly at a slope that is essentially independent of collision energy. Because regeneration compensates for decays, the ratio does not follow the simple exponential decay law used in earlier estimates. The difference between the $K^*/K$ ratio at chemical and kinetic freeze-out, divided by the universal slope of about $0.02~\mathrm{fm}^{-1}$, yields the lifetime of the hadronic stage, which the authors find to be 2-4 times larger than the STAR collaboration's previous extraction.
Load-bearing premise
The slope of about $-0.02~\mathrm{fm}^{-1}$ is extracted for central Au+Au collisions and then used as a universal constant for all centralities and energies, with no demonstration that it is centrality-independent and no estimate of its uncertainty.
Editorial extensions
If this is right
- The hadronic phase at RHIC-BES energies is roughly 2-4 times longer than the STAR collaboration's estimate, resolving the tension between resonance suppression and other lifetime probes.
- Analyses that treat resonance suppression as a pure exponential decay will systematically underestimate hadron phase durations, so the linear formula provides a corrected tool for future measurements.
- The energy independence of the slope means a single constant can be used across the beam-energy scan, simplifying extraction of lifetimes at NA61, STAR, and FAIR.
- The agreement with independent 4-8 fm/c estimates suggests that the rate-equation network captures the dominant physics of the late hadronic stage, at least for the $K^*$ channel.
Reading between the lines
- The same detailed-balance mechanism should apply to other resonance ratios such as $\rho/\pi$ or $\Delta/N$; if those ratios also evolve linearly, the universal-slope method could be extended to multi-resonance consistency checks.
- One could test the universality claim directly by comparing central and peripheral collisions: if the slope changes with system size or centrality, Eq. (12) would need a centrality-dependent slope.
- Because the slope is anchored to the assumed volume expansion law and the chemical freeze-out time from UrQMD, a direct measurement of the time dependence of $K^*/K$ (for instance through femtoscopy) could validate the linear evolution independently.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the time evolution of the K*/K ratio during the hadronic phase in Au+Au collisions at RHIC-BES energies using a network of 23 coupled rate equations with partial chemical equilibrium and detailed balance. The authors report that the K*/K ratio decreases approximately linearly rather than exponentially, because regeneration counteracts decay, and they extract an essentially energy-independent slope d(K*/K)/dt ≈ -0.02 fm^-1 (Section IV, Fig. 3). They propose an improved relation, Eq. (12), between the K*/K suppression and the hadronic-phase lifetime, and apply it to STAR data to obtain lifetimes that are 2-4 times larger than the STAR estimates based on the exponential formula in Eq. (1).
Significance. If the central claims hold, the paper offers a more realistic method for extracting the hadronic-stage duration from resonance suppression and identifies a systematic bias in previous exponential-decay estimates. The explicit rate-equation framework with detailed balance, the reproduction of centrality and energy dependence of the STAR K*/K data in Fig. 1, and the clear linear behavior shown in Fig. 3 are genuine strengths. The proposed relation Eq. (12) is simple and directly usable by STAR, NA61, and future FAIR experiments. However, the quantitative validity of the result depends on the completeness of the loss terms and on the universality of the slope, both of which need additional support before the factor-2-4 statement can be considered robust.
major comments (3)
- [Section II, Eqs. (2) and (4)] The loss rate due to collisions, L_coll, announced in Eq. (2) is never implemented in the equations actually solved. The evolution equation Eq. (4) contains only the decay/regeneration channel h* ⇌ X + h, and no collisional dissociation term is defined or discussed in Section II. In the dense hadronic medium at RHIC-BES energies, channels such as K* + N → K + N or K* + π → K_1 → ... would add an extra loss mechanism that makes d(K*/K)/dt more negative than -0.02 fm^-1. Since Eq. (12) divides the measured suppression by this slope, omitting collisional losses could make the inferred Δt_hadronic systematically too large, and the claimed factor of 2-4 relative to STAR would rest on an unverified completeness assumption. The authors should either implement L_coll in the network or provide a quantitative estimate of its magnitude and demonstrate that it is negligible for the slope.
- [Section IV, Eq. (12), Figs. 3 and 4] The slope s = d(K*/K)/dt ≈ -0.02 fm^-1 is extracted from central Au+Au collisions in Fig. 3, but Eq. (12) is then applied to all centralities and energies in Fig. 4 without demonstrating that the slope is centrality independent. The slope also depends on the volume parametrization Vch in Eq. (9), the chemical freeze-out time tch taken from UrQMD, and the transverse expansion time t_perp in Eq. (11); none of these dependencies is quantified. The paper provides no uncertainty on s and no sensitivity analysis, although any change in s linearly rescales all extracted lifetimes. The statement in Section II that t_perp has only minor influence is not supported by a calculation. Please show the slope as a function of centrality, provide an uncertainty band for s, and propagate that band into Fig. 4.
- [Section III, Fig. 1] The benchmarking against STAR data in Fig. 1 does not independently validate the slope used in Eq. (12). The kinetic freeze-out point in Fig. 3 is read off when the model's K*/K matches the same STAR data that enter the numerator of Eq. (12), and the slope is taken from the same model run. This consistency is necessary but not sufficient to certify the slope as a universal constant. An independent cross-check, such as comparing the model's K*/K evolution against a transport-code calculation with controlled resonance processes, would substantially strengthen the claim that the slope is a reliable input for lifetime extraction.
minor comments (3)
- [Fig. 1 caption] The caption states that filled symbols with bands denote the rate-equation calculations, but the meaning of the band is not defined in the text; please specify whether it reflects numerical uncertainty, centrality-bin width, or a model variation.
- [Fig. 4] The vertical axis extends to negative values and some STAR points appear at negative hadronic lifetimes, but the text does not explain how a negative duration is to be interpreted; please address this for clarity.
- [Section II] The paper refers to a network of 23 rate equations and to 'the specific rate equations described in [24]', but the included species and the relevant decay/regeneration channels are not itemized; a table listing the 23 species and their principal channels would improve reproducibility.
Circularity Check
No significant circularity: the slope in Eq. (12) is a model output, not a fit to the target lifetimes.
full rationale
This paper's central relation (Eq. 12) converts the measured K*/K suppression into a hadronic-phase duration using a slope d(K*/K)/dt ≈ -0.02 fm^-1 that is read off the rate-equation time evolution in Fig. 3. This is a model output, not a fitted parameter: the rate equations are stated in Section II with decay widths from PDG [37], initial conditions from the Cleymans parametrization [38], the volume from Pan-Pratt [40], and tch from UrQMD [41]; no parameter is adjusted to reproduce the STAR K*/K data. The STAR data are used only to identify the kinetic freeze-out time (circles in Fig. 3) and as the numerator in Eq. 12. The linear time evolution and the value of the slope are consequences of the regeneration term in Eq. (4), independent of the final lifetime values. The self-citations to [24] (rate-equation framework), [41] (UrQMD tch), and [10,42,43] (comparison estimates) provide modeling inputs or external benchmarks rather than a uniqueness theorem or a forced conclusion. The omission of Lcoll from Eq. (2) is a completeness/sensitivity concern, not a circular reduction. Therefore no load-bearing argument reduces to its own inputs.
Assumptions & free parameters
free parameters (4)
- Vch(x) coefficients =
2351.7, -10009.5, 13878.4, exponent -1.36
- t_ch (chemical freeze-out time) =
Not quoted, taken from UrQMD [41]
- t_perp (transverse expansion time) =
6.5 fm/c
- Linear slope s = d(K*/K)/dt =
0.02 fm^-1
assumptions (6)
- domain assumption Partial chemical equilibrium: total yields of stable hadrons are fixed by non-equilibrium chemical potentials, while resonances evolve by decays and regeneration.
- standard math Detailed balance: the cross section for h + X -> h* and the decay width of h* -> h + X are related by Eq. (6), with equilibrium multiplicities.
- domain assumption Decay widths and cross sections are taken from Particle Data Tables and are thermally averaged.
- domain assumption Volume expansion follows V(t) = V_ch (t/t_ch) (t_perp^2 + t^2)/(t_perp^2 + t_ch^2) from Pan and Pratt [40].
- domain assumption Chemical freeze-out conditions are given by the Cleymans parametrizations, Eqs. (7) and (8), and the centrality scaling of Vch with charged particle multiplicity.
- ad hoc to paper The slope d(K*/K)/dt is essentially independent of collision energy and centrality, so a single value can be used in Eq. (12).
Cite this review
Pith. "Pith review of Determining the Duration of the Hadronic Stage at RHIC-BES Energies via Resonance Suppression Using a Full Set of Rate Equations." pith.science (2026). https://pith.science/paper/UG4NACT7
@misc{pith2026250103893,
author = {Pith},
title = {Pith review of: Determining the Duration of the Hadronic Stage at RHIC-BES Energies via Resonance Suppression Using a Full Set of Rate Equations},
year = {2026},
howpublished = {\url{https://pith.science/paper/UG4NACT7}},
note = {Machine review of arXiv:2501.03893}
}
abstract
We present realistic estimates for the duration of the hadronic stage in central Au+Au reactions in the RHIC-BES energy regime. To this aim, we employ a full set of coupled rate equations to describe the time evolution of the system from chemical to kinetic freeze-out. Combined with the recently measured data by the STAR collaboration on $K^*/K$ ratios, we show that the previous estimates substantially underestimated the duration of this stage due to the omission of the regeneration of hadron resonances. We provide an improved relation between the $K^*/K$ ratio at chemical and kinetic freeze-out and the life time of the hadronic phase. The calculated improved life times are now in line with estimates from other methods and are relevant for the NA61 and STAR collaborations and for upcoming experiments at the FAIR facility.
Figures
Reference graph
Works this paper leans on
-
[1]
A. Sorensen, K. Agarwal, K. W. Brown, Z. Chajecki, P. Danielewicz, C. Drischler, S. Gandolfi, J. W. Holt, M. Kaminski and C. M. Ko, et al.Prog. Part. Nucl. Phys. 134, 104080 (2024) doi:10.1016/j.ppnp.2023.104080 [arXiv:2301.13253 [nucl-th]]
arXiv 2024
-
[2]
M. Bleicher and J. Aichelin, Phys. Lett. B 530, 81-87 (2002) doi:10.1016/S0370-2693(02)01334-5 [arXiv:hep- ph/0201123 [hep-ph]]
-
[3]
A. Andronic, P. Braun-Munzinger and J. Stachel, Nucl. Phys. A 772, 167-199 (2006) doi:10.1016/j.nuclphysa.2006.03.012 [arXiv:nucl- th/0511071 [nucl-th]]
-
[4]
R. Hanbury Brown and R. Q. Twiss, Nature 178, 1046- 1048 (1956) doi:10.1038/1781046a0
-
[5]
M. A. Lisa, S. Pratt, R. Soltz and U. Wiede- mann, Ann. Rev. Nucl. Part. Sci. 55, 357-402 (2005) doi:10.1146/annurev.nucl.55.090704.151533 [arXiv:nucl- ex/0505014 [nucl-ex]]
arXiv 2005
-
[6]
G. Torrieri and J. Rafelski, Phys. Lett. B 509, 239-245 (2001) doi:10.1016/S0370-2693(01)00492-0 [arXiv:hep- ph/0103149 [hep-ph]]
-
[7]
J. Rafelski, J. Letessier and G. Torrieri, Phys. Rev. C 64, 054907 (2001) [erratum: Phys. Rev. C 65, 069902 (2002)] doi:10.1103/PhysRevC.64.054907 [arXiv:nucl-th/0104042 [nucl-th]]
arXiv 2001
-
[8]
C. Le Roux, F. S. Navarra and L. M. Abreu, Phys. Lett. B 817, 136284 (2021) doi:10.1016/j.physletb.2021.136284 [arXiv:2101.07302 [hep-ph]]
arXiv 2021
Show all 44 references
-
[9]
A. G. Knospe, C. Markert, K. Werner, J. Steinheimer and M. Bleicher, Phys. Rev. C 104, no.5, 054907 (2021) doi:10.1103/PhysRevC.104.054907 [arXiv:2102.06797 [nucl-th]]
2021 arXiv
-
[10]
A. G. Knospe, C. Markert, K. Werner, J. Steinheimer and M. Bleicher, Phys. Rev. C 93, no.1, 014911 (2016) doi:10.1103/PhysRevC.93.014911 [arXiv:1509.07895 [nucl-th]]
2016 arXiv
-
[11]
Ilner, J
A. Ilner, J. Blair, D. Cabrera, C. Markert and E. Bratkovskaya, Phys. Rev. C 99, no.2, 024914 (2019) doi:10.1103/PhysRevC.99.024914 [arXiv:1707.00060 [hep-ph]]
2019 arXiv
- [12]
-
[13]
Chabane, L
A. Chabane, L. Engel, T. Reichert, J. Steinheimer and M. Bleicher, [arXiv:2409.08639 [hep-ph]]
-
[14]
Adams et al
J. Adams et al. [STAR], Phys. Rev. C 71, 064902 (2005) doi:10.1103/PhysRevC.71.064902 [arXiv:nucl- ex/0412019 [nucl-ex]]
2005
-
[15]
Markert, J
C. Markert, J. Phys. G 31, S169-S178 (2005) doi:10.1088/0954-3899/31/4/021 [arXiv:nucl-ex/0503013 [nucl-ex]]
2005 arXiv
-
[16]
B. I. Abelev et al. [STAR], Phys. Rev. C 78, 044906 (2008) doi:10.1103/PhysRevC.78.044906 [arXiv:0801.0450 [nucl-ex]]
2008 arXiv
-
[17]
Agakishiev et al
G. Agakishiev et al. [HADES], Eur. Phys. J. A 49, 34 (2013) doi:10.1140/epja/i2013-13034-7
2013 doi
-
[18]
A. G. Knospe [ALICE], J. Phys. Conf. Ser. 612, no.1, 012064 (2015) doi:10.1088/1742-6596/612/1/012064 [arXiv:1501.03797 [nucl-ex]]
2015 arXiv
-
[19]
Acharya et al
S. Acharya et al. [ALICE], Phys. Rev. C 99, 024905 (2019) doi:10.1103/PhysRevC.99.024905 [arXiv:1805.04361 [nucl-ex]]
2019 arXiv
-
[20]
Acharya et al
S. Acharya et al. [ALICE], Eur. Phys. J. C 83, no.5, 351 (2023) doi:10.1140/epjc/s10052-023-11475-1 [arXiv:2205.13998 [nucl-ex]]
2023 arXiv
-
[21]
Acharya et al
S. Acharya et al. [ALICE], Phys. Rev. C 109, no.4, 044902 (2024) doi:10.1103/PhysRevC.109.044902 [arXiv:2308.16119 [nucl-ex]]
2024 arXiv
-
[22]
Koz lowski [NA61/SHINE], [arXiv:2409.20229 [nucl- ex]]
B. Koz lowski [NA61/SHINE], [arXiv:2409.20229 [nucl- ex]]
-
[23]
Motornenko, V
A. Motornenko, V. Vovchenko, C. Greiner and H. Stoecker, Phys. Rev. C 102, no.2, 024909 (2020) doi:10.1103/PhysRevC.102.024909 [arXiv:1908.11730 [hep-ph]]
2020 arXiv
-
[24]
Neidig, K
T. Neidig, K. Gallmeister, C. Greiner, M. Bleicher and V. Vovchenko, Phys. Lett. B 827, 136891 (2022) doi:10.1016/j.physletb.2022.136891 [arXiv:2108.13151 [hep-ph]]
2022
-
[25]
Adler et al
C. Adler et al. [STAR], Phys. Rev. C 66, 061901 (2002) doi:10.1103/PhysRevC.66.061901 [arXiv:nucl- ex/0205015 [nucl-ex]]
2002
-
[26]
Abdallah et al
M. Abdallah et al. [STAR], Phys. Rev. C 107, no.3, 034907 (2023) doi:10.1103/PhysRevC.107.034907 [arXiv:2210.02909 [nucl-ex]]
2023
-
[27]
Torrieri and J
G. Torrieri and J. Rafelski, J. Phys. G 28, 1911- 1920 (2002) doi:10.1088/0954-3899/28/7/346 [arXiv:hep- ph/0112195 [hep-ph]]
2002
-
[28]
Bebie, P
H. Bebie, P. Gerber, J. L. Goity and H. Leutwyler, Nucl. Phys. B 378, 95-128 (1992) doi:10.1016/0550- 3213(92)90005-V
1992 doi
-
[29]
Xu and R
X. Xu and R. Rapp, Eur. Phys. J. A 55, no.5, 68 (2019) doi:10.1140/epja/i2019-12757-7 [arXiv:1809.04024 [nucl- th]]
2019 arXiv
-
[30]
Vovchenko, K
V. Vovchenko, K. Gallmeister, J. Schaffner-Bielich and C. Greiner, Phys. Lett. B 800, 135131 (2020) doi:10.1016/j.physletb.2019.135131 [arXiv:1903.10024 [hep-ph]]
2020
-
[31]
Tomasik, P
B. Tomasik, P. Hillmann and M. Bleicher, Phys. Rev. C 104, no.4, 044907 (2021) doi:10.1103/PhysRevC.104.044907 [arXiv:2107.03830 [nucl-th]]
2021 arXiv
-
[32]
P. Koch, B. Muller and J. Rafelski, Phys. Rept. 142, 167-262 (1986) doi:10.1016/0370-1573(86)90096-7
1986 doi
-
[33]
H. W. Barz, B. L. Friman, J. Knoll and H. Schulz, Nucl. Phys. A 484, 661-684 (1988) [erratum: Nucl. Phys. A 492, 663-663 (1989)] doi:10.1016/0375-9474(88)90315-6
1988 doi
-
[34]
Noronha-Hostler, C
J. Noronha-Hostler, C. Greiner and I. A. Shovkovy, Phys. Rev. Lett. 100, 252301 (2008) 6 doi:10.1103/PhysRevLett.100.252301 [arXiv:0711.0930 [nucl-th]]
2008 arXiv
-
[35]
Noronha-Hostler, M
J. Noronha-Hostler, M. Beitel, C. Greiner and I. Shovkovy, Phys. Rev. C 81, 054909 (2010) doi:10.1103/PhysRevC.81.054909 [arXiv:0909.2908 [nucl-th]]
2010 arXiv
-
[36]
E. W. Kolb, Front. Phys. 69, 1-547 (1990) Taylor and Francis, 2019, ISBN 978-0-429-49286-0, 978-0-201-62674- 2 doi:10.1201/9780429492860
1990 doi
-
[37]
P. A. Zyla et al. [Particle Data Group], PTEP 2020, no.8, 083C01 (2020) doi:10.1093/ptep/ptaa104
2020 doi
-
[38]
Cleymans, H
J. Cleymans, H. Oeschler, K. Redlich and S. Wheaton, J. Phys. G 32, S165-S170 (2006) doi:10.1088/0954- 3899/32/12/S21 [arXiv:hep-ph/0607164 [hep-ph]]
2006 arXiv
-
[39]
Andronic, P
A. Andronic, P. Braun-Munzinger and J. Stachel, Acta Phys. Polon. B 40, 1005-1012 (2009) [arXiv:0901.2909 [nucl-th]]
2009 arXiv
-
[40]
Pan and S
Y. Pan and S. Pratt, Phys. Rev. C 89, no.4, 044911 (2014) doi:10.1103/PhysRevC.89.044911
2014 doi
-
[41]
Reichert, G
T. Reichert, G. Inghirami and M. Bleicher, Eur. Phys. J. A 56, no.10, 267 (2020) doi:10.1140/epja/s10050-020- 00273-y [arXiv:2007.06440 [nucl-th]]
2020 arXiv
-
[42]
Q. Li, M. Bleicher and H. Stocker, J. Phys. G 34, 2037-2044 (2007) doi:10.1088/0954-3899/34/9/012 [arXiv:0706.2091 [nucl-th]]
2007 arXiv
-
[43]
Graef, M
G. Graef, M. Bleicher and Q. Li, Phys. Rev. C 85, 044901 (2012) doi:10.1103/PhysRevC.85.044901 [arXiv:1203.4071 [nucl-th]]
2012 arXiv
-
[44]
Lisa, Acta Phys
M. Lisa, Acta Phys. Polon. B 47, 1847 (2016) doi:10.5506/APhysPolB.47.1847 [arXiv:1607.06188 [nucl- th]]
2016 arXiv
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.