REVIEW 3 major objections 4 minor 18 references
The Equation of State and Multiparticle Production
T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The net-baryon distribution seen at LHC midrapidity is fixed before thermalization, at a freezeout temperature of at least about 300 MeV, so neither the free-quark model nor the hadron resonance gas fully explains it.
desk verdict A sharp speculative claim built on a real finite-volume calculation; the transport premise needs quantitative support before the central inference holds. 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 key object is the newly introduced net-baryon number freezeout temperature $T_{Bf}$, the temperature of the fireball at the moment $\tau\sim2$ fm/c when the net-baryon number in the midrapidity domain $|y|\lesssim4$ becomes conserved; before that moment the fireball is described grand canonically and after it canonically. The formal machinery is the fugacity expansion $Z_{GC}(\theta)=\sum_B Z_C(B)e^{B\theta}$ and the cumulant generating function $K(t)=\ln[Z_{GC}(t)/Z_{GC}(0)] = (\hat p(t)-\hat p(0))\nu$, which convert the equation of state into predictions for $\kappa_2,\kappa_4,\ldots$. Two distributions are compared: the free-quark-model mass function, whose large-$B$ tail behaves as $\exp[-(3N_C/4)\sqrt[3]{3\pi^2B^4/(\nu N_F)}]$, and the HRG Skellam-type distribution $P(B)\propto e^{-(N_b+N_{\bar b})}(N_b/N_{\bar b})^{B/2} I_B(2\sqrt{N_bN_{\bar b}})$. The freezeout temperature carries the argument because once $B$ is fixed, the hadronic stage cannot rearrange the distribution, so the observed cumulants report the early high-temperature ensemble.
What would settle it
A transport simulation of Pb+Pb at $\sqrt{s_{NN}} = 5.02$ TeV with a realistic baryon diffusion coefficient that shows the net-baryon number inside $|y| < 4$ changing appreciably after 2 fm/c would falsify the early-freezeout premise; alternatively, a measurement of net-baryon cumulant ratios whose dependence on rapidity-window width matches the HRG-with-conservation prediction would show the distribution is not frozen before thermalization.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that the probability distribution of the net-baryon number at midrapidity at LHC energies is established at an early, high-temperature stage and then frozen: the net-baryon freezeout temperature is $T_{Bf} \gtrsim 300$ MeV, well above the pseudocritical temperature of the chiral crossover. The argument has three legs. The measured baryon chemical potential is below about 1 MeV, so baryons from the incoming nuclei cannot be entering midrapidity at late times. The net-baryon rapidity distribution at lower energies has two stopping peaks plus a plateau, and by $\sqrt{s_{NN}}\approx 5$ TeV the domain $|y|\lesssim 4$ should have zero average net-baryon number, with the number fixed at $\tau\approx 2$ fm/c. In that early high-temperature regime the equation of state is close to the free-quark model, giving cumulant ratios $\chi_4/\chi_2 = 1/(2\pi^2 N_c^2)$ and $\chi_n/\chi_2=0$ for $n>2$, whereas the measured ratios match the Skellam-like HRG values $\chi_{2n}/\chi_2 = 1$. The paper infers that neither model is complete: the free-quark model suppresses large fluctuations too strongly, and the hadron resonance gas would require rapid baryon transport into midrapidity that the small chemical potential rules out. The observed distribution must therefore be formed before thermalization, with independent production of quarks and antiquarks, for which the sea quarks of the colliding nuclei are a candidate source.
Load-bearing premise
The argument stands on the assumption that after about 2 fm/c the net-baryon number in the rapidity window $|y| \lesssim 4$ is fixed, because baryons from the incoming nuclei do not diffuse into it in appreciable numbers — a premise inferred from the measured near-zero chemical potential and the shape of the rapidity distributions, not computed from a transport model.
Editorial extensions
If this is right
- If the central claim is correct, the measured net-baryon cumulants at $\sqrt{s_{NN}}\approx 5$ TeV report the equation of state at $T\gtrsim 300$ MeV, so lattice-QCD comparisons at those temperatures, not at $T\approx155$ MeV, are the relevant check.
- The hadron resonance gas model, despite matching the data, would describe the wrong dynamics: it requires baryon exchange with the rest of phase space that the measured chemical potential near zero forbids.
- The distribution in $B$ must be produced by early-time, pre-thermal production of independent quark and antiquark excitations, making sea quarks of the colliding nuclei a concrete candidate mechanism.
- Cumulant ratios at LHC energies should differ from those at beam energies below 200 GeV, where baryon stopping populates midrapidity and late-time HRG logic applies.
- Finite-volume and transverse-momentum acceptance corrections to the free-quark model are sizable but still leave a dramatic gap from the observed cumulant ratios.
Reading between the lines
- Inference: A transport calculation with a realistic baryon diffusion coefficient at 5.02 TeV would directly test the freezeout-time premise; if diffusion moves net baryons into $|y|<4$ after 2 fm/c, the hadron resonance gas with baryon-number conservation corrections could survive with $T_f\approx155$ MeV.
- Inference: The sea-quark hypothesis predicts that the net-baryon distribution should depend on the nuclear content of the incoming ions, so comparing Pb+Pb with p+Pb or with isobar collisions at the same energy could discriminate it from late-time Skellam noise.
- Inference: If $T_{Bf}$ is genuinely about 300 MeV, the same early-frozen fluctuations should be insensitive to the width of the rapidity acceptance once $|y|\lesssim4$ is covered, whereas HRG-with-conservation predicts a characteristic acceptance dependence; this is a direct experimental discriminator.
- Inference: The argument implies that the traditional chemical freezeout at about 155 MeV describes abundances of hadron species but not the net-baryon cumulants; separating the two freezeouts could change how the LHC data are compared to lattice QCD.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript argues that the net-baryon number distribution measured by ALICE in Pb-Pb collisions at sqrt(s_NN) ~ 5 TeV cannot be explained by either the free-quark model (FQM) or the hadron resonance gas (HRG) model at freezeout. It introduces the concept of a net-baryon-number freezeout temperature T_Bf, estimated to be at least about 300 MeV, at which the net-baryon number in the midrapidity domain becomes fixed. The paper performs an explicit finite-volume computation of the FQM partition function in Section 3 and shows that finite-volume and transverse-momentum cutoffs do not remove the disagreement between FQM and the ALICE data. The central positive inference, that the observed distribution is formed before thermalization, relies on the assumption that baryon transport into and within the experimental acceptance is negligible after tau ~ 2 fm/c. The manuscript itself repeatedly labels the arguments qualitative and speculative, and it proposes sea-quark contributions as a possible but undeveloped explanation.
Significance. If the central inference were correct, ALICE net-baryon cumulants would provide a direct probe of the early, high-temperature phase of the collision, and neither the standard HRG nor the FQM would suffice to describe the data. This would be an interesting and potentially important proposal. The paper also contains a useful concrete negative result: the explicit finite-volume computation in Section 3 shows that FQM predictions remain far from the ALICE cumulant ratios even after finite-volume and acceptance effects are included, supporting the claim that FQM alone cannot explain the data. However, the positive claim about early formation depends on an unquantified transport premise, so the paper is best read as a speculative proposal rather than an established conclusion.
major comments (3)
- [Section 4, second paragraph] The inference that the observed fluctuations are related to T_Bf requires that the net-baryon number in the actual ALICE acceptance be frozen after tau ~ 2 fm/c. The manuscript establishes only that the average net-baryon number vanishes and is conserved in the larger domain |y| < y_c ~ 4. Conservation in that domain does not imply that B in a narrower midrapidity acceptance is constant: baryon-number fluctuations can migrate across the acceptance boundary inside the domain, changing the measured B while preserving the total in D. No calculation of the baryon-number diffusion current or of the relaxation of acceptance-level cumulants is provided. This missing step is load-bearing for the central claim that the ALICE cumulants probe the early high-temperature phase.
- [Section 4, third paragraph] The argument that final-stage baryon transport into midrapidity is excluded by the measured chemical potential mu < 1 MeV [16] is not sufficient. The chemical potential constrains the first moment of the net-baryon distribution, namely the mean net-baryon density, but it does not constrain the variance or higher cumulants in the acceptance. A zero-mean diffusion process could generate HRG-like cumulant ratios while maintaining mu ~ 0, for instance through a baryon-number conservation correction of the Bzdak-Koch-Skokov type [11] with a small acceptance-mixing parameter. Therefore the paper's rejection of the HRG explanation at T_f ~ T_pc is not established by the cited argument.
- [Section 5, Conclusions] The claim that the FQM distribution 'cannot be rapidly rearranged' to the HRG distribution 'in view of the net-baryon number conservation' is too strong. Conservation of the total net-baryon number in the full domain does not freeze the probability distribution in a subvolume; cumulants can evolve through diffusion even while the total is conserved. The conclusion that the observed distribution should be attributed to an early stage therefore depends on an unproven invariance of the acceptance-level cumulants under the later evolution. Without an estimate of the relaxation time for acceptance-level baryon cumulants, the central conclusion remains conditional.
minor comments (4)
- [Throughout] There are numerous typographical errors, including 'intiguing' in the Introduction, 'paramtric' in Section 2, 'probabulity' in Section 2, 'freezout' in several places, and 'nonvanising' and 'midtapidity' in Section 4. These should be corrected in a revision.
- [Section 3, paragraph on the cube division] The assumption that the fireball can be divided into about 40 independent cubes, each in contact with a thermostat at T ~ 300 MeV, is introduced without discussion. Since the finite-volume results in Table 1 are a central quantitative part of the negative FQM claim, the sensitivity of the cumulant ratios to the number of cubes and to the independence assumption should be stated.
- [Eq. (13)] Since odd-order cumulants vanish by C-parity, the statement 'chi_n/chi_2 = 0 at n > 2' should be phrased as 'chi_{2n}/chi_2 = 0 for n > 2' to avoid ambiguity about odd n.
- [Section 2, Eq. (9)] The normalization constant C(nu) in the FQM probability mass function is stated but not given in closed form; for reproducibility of the finite-volume computation it would help to state how C(nu) is determined in the direct computation in Section 3.
Circularity Check
No significant circularity: the finite-volume FQM computation and HRG comparison rest on external data and standard model parameters; the central conclusion depends on an external CGC transport premise, not on a fit or self-citation chain.
full rationale
The paper's quantitative FQM prediction (Table 1) is obtained by direct computation of the zero-triality partition function for massless quarks with standard N_F=2, N_C=3 and spin 1/2, with no parameter fitted to the ALICE data. The ALICE cumulant ratios are an external experimental benchmark, and the HRG/Skellam predictions are standard literature results. The central inference (T_Bf >~ 300 MeV) follows from the external CGC baryon-stopping picture [15] plus the measured mu < 1 MeV bound [16]; whether that transport premise is quantitatively correct is a physics-risk issue, not a circular reduction. The paper explicitly labels the sea-quark mechanism as speculative and says the LHC B-distribution is 'yet to be explained.' Self-citations [4] and [7] supply an FQM probability formula and lattice EoS inputs, but the paper's own finite-volume computation does not depend on them, and nothing is fitted and then renamed as a prediction. Therefore no step reduces by construction to its own input.
Assumptions & free parameters
free parameters (4)
- ν ~ 30 (dimensionless volume per cube) =
~30
- τ ~ 2 fm/c (net-baryon number freezeout time) =
2 fm/c
- R = 7 fm (fireball radius) =
7 fm
- yc ~ 4 (rapidity window boundary) =
4
assumptions (5)
- standard math The grand canonical partition function determines the canonical net-baryon probability through the fugacity expansion Z_GC(θ) = Σ_B Z_C(B) e^{Bθ}.
- domain assumption At T > T_RW ~ 200 MeV the QCD equation of state follows the free quark model, as supported by lattice QCD.
- ad hoc to paper At LHC energies the net-baryon number in the rapidity window |y| less than about 4 is zero on average and conserved after τ ~ 2 fm/c.
- ad hoc to paper The measured tiny baryon chemical potential (μ < 1 MeV) implies negligible net-baryon transport into midrapidity at the final stage.
- standard math C-parity symmetry makes the net-baryon probability even, P_B = P_{-B}.
invented entities (1)
-
Net-baryon number freezeout temperature T_Bf
Cite this review
Pith. "Pith review of The Equation of State and Multiparticle Production." pith.science (2026). https://pith.science/paper/ODQJ6LSH
@misc{pith2026241118190,
author = {Pith},
title = {Pith review of: The Equation of State and Multiparticle Production},
year = {2026},
howpublished = {\url{https://pith.science/paper/ODQJ6LSH}},
note = {Machine review of arXiv:2411.18190}
}
read the original abstract
We discuss the distribution of fireballs produced in heavy-ion collisions in the net-baryon number and argue that neither the Free-Quark Model (FQM) nor the Hadron Resonance Gas (HRG) model can provide a comprehensive explanation of the distribution observed at the LHC. The concept of net-baryon number freezeout temperature is suggested and the role of sea quarks as a possible source of net-baryon number fluctuations is emphasized.
Reference graph
Works this paper leans on
-
[16]
Acharya S. et al. [ALICE Collaboration] Measurements of chemical po- tentials in Pb-Pb collisions at √sNN = 5 . 02 TeV. – – 2023. – – 11. – – arXiv:2311.13332. 10
arXiv 2023
-
[11]
Baryon number conservation and the cu- mulants of the net proton distribution // Phys
Bzdak A., Koch V., Skokov V. Baryon number conservation and the cu- mulants of the net proton distribution // Phys. Rev. C. – – 2013. – – V. 87, no. 1. – – P. 014901. – – arXiv:1203.4529 [hep-ph]
arXiv 2013
-
[1]
Behera N.K. et al. [ALICE Collaboration] Higher moment fluctuations of identified particle distributions from ALICE // Nucl. Phys. A. – – 2019 . – – V. 982. – – P. 851–854. – – arXiv:1807.06780
work page Pith review arXiv 2019
-
[2]
Gauge Theories With Imaginary Chemical Potential and the Phases of QCD // Nucl
Roberge A., Weiss N. Gauge Theories With Imaginary Chemical Potential and the Phases of QCD // Nucl. Phys. B. – – 1986. – – V. 275. – – P. 734–745
work page 1986
-
[3]
Dyons and Roberge - Weiss transition in lattice QCD
Bornyakov V.G., Boyda D.L., Goy V.A., Ilgenfritz E.M., Mart emyanov B.V., Molochkov A.V., Nakamura A., Nikolaev A.A., Zakharov V.I. Dyons and Roberge - Weiss transition in lattice QCD // EPJ Web Conf. – – 2017. – – V. 137. – – P. 03002. – – arXiv:1611.07789
work page Pith review arXiv 2017
-
[4]
Numerical Study of the Roberge-Weiss Transition
Bornyakov V.G., Gerasimeniuk N.V., Goy V.A., Korneev A.A., Molochkov A.V., Nakamura A., Rogalyov R.N. Numerical study of the Roberge-Weiss transition // Phys. Rev. D. – – 2023. – – V. 107, no. 1. – – P. 014508 . – – arXiv:2203.06159
work page Pith review arXiv 2023
-
[5]
Triality in QCD at zero and finite tem- perature: A New direction // Nucl
Faber M., Borisenko O., Zinovev G. Triality in QCD at zero and finite tem- perature: A New direction // Nucl. Phys. B. – – 1995. – – V. 444. – – P . 563–
work page 1995
-
[6]
[HotQCD Collaboration] The QCD Equation of State to O(µ 6 B) from Lattice QCD // Phys
Bazavov A., et al. [HotQCD Collaboration] The QCD Equation of State to O(µ 6 B) from Lattice QCD // Phys. Rev. D. – – 2017. – – V. 95, no. 5. – – P. 054504. – – arXiv:1701.04325
arXiv 2017
Show all 18 references
-
[7]
Net-Baryon Probability Distributions from Lattice Simulations // Particles
Rogalyov R., Goy V. Net-Baryon Probability Distributions from Lattice Simulations // Particles. – – 2023. – – V. 6, no. 3. – – P. 771–783
2023
-
[8]
Net- proton probability distribution in heavy ion collisions // Phys
Braun-Munzinger P., Friman B., Karsch F., Redlich K., Skoko v V. Net- proton probability distribution in heavy ion collisions // Phys. Rev. C. – –
-
[9]
Net-charge probability distributions in heavy ion collisions at chem- ical freeze-out // Nucl
Braun-Munzinger P., Friman B., Karsch F., Redlich K., Skoko v V. Net-charge probability distributions in heavy ion collisions at chem- ical freeze-out // Nucl. Phys. A. – – 2012. – – V. 880. – – P. 48–64. – – arXiv:1111.5063 [hep-ph]
2012 arXiv
-
[10]
Probing freeze-out conditions in heavy ion collisions with moments of charge fluctuations // Phys
Karsch F., Redlich K. Probing freeze-out conditions in heavy ion collisions with moments of charge fluctuations // Phys. Lett. B. – – 2011. – – V. 695. – – P. 136–142. – – arXiv:1007.2581 [hep-ph]
2011 arXiv
-
[12]
Probing initial baryon stopping and equa- tion of state with rapidity-dependent directed flow of identified par ticles // Phys
Du L., Shen C., Jeon S., Gale C. Probing initial baryon stopping and equa- tion of state with rapidity-dependent directed flow of identified par ticles // Phys. Rev. C. – – 2023. – – V. 108, no. 4. – – P. L041901. – – arXiv:2211.16408
2023 arXiv
-
[13]
Rapidity scan approach for net-baryon cumulants with a statistical thermal model
Li J., Du L., Shi S. Rapidity scan approach for net-baryon cumulants with a statistical thermal model. – – 2023. – – 11. – – arXiv:2311.11374. 9
2023 arXiv
-
[14]
Videbaek F. et al. [BRAHMS Collaboration] Overview and Recent Results from BRAHMS // Nucl. Phys. A. – – 2009. – – V. 830. – – P. 43C–50C. – – arXiv:0907.4742 [nucl-ex]
2009 arXiv
-
[15]
Spacetime picture of baryon stopping in the color-glass condensate // Phys
McLerran L.D., Schlichting S., Sen S. Spacetime picture of baryon stopping in the color-glass condensate // Phys. Rev. D. – – 2019. – – V. 99, n o. 7. – – P. 074009. – – arXiv:1811.04089
2019 arXiv
-
[576]
– – arXiv:hep-ph/9504264
-
[2011]
– – V. 84. – – P. 064911. – – arXiv:1107.4267 [hep-ph]
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.