REVIEW 3 major objections 5 minor 35 references
Temperature-fluctuation cumulants in interacting hadron gases match lattice QCD below the chiral crossover, split along freeze-out at low energies, and trace the nuclear liquid-gas Widom line through a c2 minimum and a c3 sign change.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-04 01:27 UTC pith:EUMIYBRO
load-bearing objection First HRG calculation of temperature-fluctuation cumulants; plausible and honest, but the nuclear liquid-gas fingerprint rests on an unvalidated high-density VDW extrapolation. the 3 major comments →
Interaction fingerprints in temperature-fluctuation cumulants from the QCD crossover to nuclear liquid-gas criticality
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that temperature cumulants of a hadron resonance gas are nearly insensitive to the interaction scheme at high collision energies and at zero baryon density, but become sharply model-dependent in the low-energy freeze-out region. The van der Waals model, which reproduces the nuclear ground state, predicts that c2 has a minimum along the nuclear liquid-gas Widom line and that c3 changes sign across it, with the minimum deepening and moving to lower temperatures as the chemical potential approaches the liquid-gas critical point near (Tc, mu_c) ≈ (19.7, 908) MeV. The paper identifies the baryon sector as the carrier of this interaction dependence through a meson-baryon secto
What carries the argument
The machinery is the set of temperature-fluctuation cumulants c_n defined through dimensionless pressure derivatives χ_n = T^(n−4) ∂^n p/∂T^n at fixed chemical potentials, with c2 = 1/χ2, c3 = −χ3/χ2^3, and c4 = 3χ3^2/χ2^5 − χ4/χ2^4. These formulas convert the equation of state into observables whose sensitivity grows with each derivative, so interaction effects that are small in the pressure show up clearly in the cumulants. The van der Waals HRG with parameters a = 329 MeV fm^3 and b = 3.42 fm^3, fixed by the nuclear ground state, supplies a nuclear liquid-gas critical point and a Widom line along which the variance minimum and the c3 sign change develop. A sector decomposition into meson
Load-bearing premise
The central liquid-gas fingerprint rests on the assumption that the van der Waals hadron resonance gas with parameters fixed by the nuclear ground state (a = 329 MeV fm^3, b = 3.42 fm^3) remains quantitatively valid for baryonic matter at the densities reached along low-energy freeze-out (mu_B up to about 910 MeV); the paper does not validate this equation of state against finite-density lattice QCD or other independent high-density data, and the freeze-out parametrization is
What would settle it
A finite-density lattice QCD calculation of the baryon pressure at chemical potentials mu_B = 850–910 MeV and temperatures T = 20–50 MeV would settle whether the van der Waals equation of state is approximately correct; if the lattice pressure deviates markedly from the VDW prediction, the c2 minimum and c3 sign change are model artifacts. Alternatively, measurements of mean-pT fluctuation cumulants in fixed-target Au+Au collisions at sqrt(s_NN) = 2–3 GeV that fail to show the predicted c2 minimum and positive-c3 band as the trajectory approaches the Widom line would falsify the fingerprint.
If this is right
- At zero baryon density, c2 is fixed directly by the lattice specific heat, and the interacting HRG models match lattice QCD on the hadronic side of the chiral crossover (up to T ≈ 155–165 MeV); above that temperature the hadronic description ceases to apply.
- Along the chemical freeze-out curve the three models are close and nearly energy independent for sqrt(s_NN) ≥ 20 GeV, but separate strongly below about 15 GeV; the separation is carried almost entirely by the baryon sector and reflects the competition between repulsive and attractive interactions.
- In the van der Waals model, c2 develops a minimum that deepens as the chemical potential approaches the nuclear liquid-gas critical point, and c3 changes sign in a narrow band that follows the high-temperature side of the Widom line; this pattern is a hadronic-baseline signature, distinct from expected chiral-critical non-monotonic structures.
- The low-energy continuation of the chemical freeze-out curve runs parallel to and only 5–9 MeV below the Widom line for sqrt(s_NN) between 2 and 3 GeV, so low-energy fixed-target collisions probe a region where temperature cumulants may be sensitive to nuclear liquid-gas thermodynamics.
- Even the ideal HRG produces a mild non-monotonicity along the freeze-out curve (a shallow c2 minimum near sqrt(s_NN) ≈ 11 GeV) purely from the trajectory, and the extrema of c2, c3, and c4 fall inside the fixed-target window of 3.0–7.7 GeV where non-monotonic pT correlations have been reported; the paper positions these cumulants as a quantitative hadronic baseline for such measurements.
Where Pith is reading between the lines
- Editorial inference: If the predicted c2 minimum and c3 sign change survive embedding in a realistic dynamical evolution, event-by-event mean-pT cumulants could become a direct experimental probe of the nuclear liquid-gas Widom line at low collision energies, complementing conserved-charge cumulants.
- Editorial inference: The sector decomposition suggests a testable separation: because meson cumulants are model-independent and vary slowly with energy while baryon cumulants grow and interact, the energy dependence of the third or fourth cumulant of mean-pT fluctuations could in principle constrain the strength of baryonic repulsion and attraction separately.
- Editorial inference: The paper does not validate the van der Waals equation of state against finite-density lattice QCD at mu_B ≈ 850–910 MeV; a lattice computation of the baryon pressure in that region would either substantiate the Widom-line fingerprint or show it to be a model artifact.
- Editorial inference: The same cumulant construction could be applied to equations of state that include strangeness or that target the chiral critical region; the paper's finding that temperature derivatives amplify interaction effects suggests that higher-order cumulants of mean-pT distributions may be more discriminating than the variance alone in beam-energy scans.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript defines scaled temperature-fluctuation cumulants c2 through c6 as nonlinear functions of temperature derivatives of the pressure (Eqs. (2)-(6)) and evaluates them in three hadronic models: the ideal hadron resonance gas, a lattice-constrained baryon-only excluded-volume HRG, and a baryonic van der Waals HRG. It compares these models with lattice QCD at zero baryon chemical potential, follows the cumulants along the Andronic chemical freeze-out curve, and analyzes the nuclear liquid-gas region. The central results are that the models agree with lattice thermodynamics up to the chiral crossover, separate below about 15 GeV along freeze-out, and, in the VDW-HRG, c2 develops a minimum and c3 changes sign near the nuclear liquid-gas Widom line. The paper frames these as interaction fingerprints that may be relevant to STAR fixed-target measurements of mean transverse-momentum fluctuations.
Significance. This is the first HRG-based calculation of temperature-fluctuation cumulants, and the derivative observables are shown to be considerably more sensitive to interaction physics than the pressure itself. The sector decomposition in Appendix A.1 and the direct c2 comparison to published lattice specific-heat data are genuine strengths. The central liquid-gas prediction is concrete and falsifiable with fixed-target data. However, the significance is conditional: the zero-density comparison is solid, while the headline liquid-gas fingerprint rests on a high-density VDW extrapolation that is not independently validated, and the low-energy freeze-out trajectory is a continuation without propagated uncertainty.
major comments (3)
- [Sec. 4.4, Eq. (9), Fig. 4] The c2 minimum and c3 sign change are computed with VDW-HRG parameters a=329 MeV fm^3 and b=3.42 fm^3 fixed solely by the nuclear ground state. The region mu_B ~ 700-910 MeV and T ~ 20-100 MeV is not confronted with independent finite-density constraints, and Sec. 5 itself calls for trajectories constrained by finite-density lattice QCD. Because these features are the paper's headline result, the revision should include a sensitivity analysis of a and b over the range allowed by nuclear ground-state properties and, where possible, a comparison with an alternative high-density equation of state. Without such tests, the abstract's claim that temperature cumulants connect to thermodynamic structures 'in the QCD phase diagram' is stronger than the evidence presented.
- [Sec. 4.3/4.4, Eq. (10), Fig. 4(c)] The chemical freeze-out parametrization of Eq. (10) is continued below sqrt(s_NN)=3 GeV without propagating any uncertainty. The continuation lies only 5-9 MeV below the Widom line, so the statement that fixed-target energies probe nuclear liquid-gas thermodynamics depends on this extrapolation. Please specify the fitted range of Eq. (10) and either propagate the thermal-fit uncertainties or compare with an alternative low-energy freeze-out curve obtained from HADES/STAR thermal fits.
- [Sec. 4.2, Fig. 2, Appendix A.2] The lattice bands for c3 and c4 are obtained by differentiating the HotQCD and WB analytic parametrizations; they do not include a full propagation of the published lattice uncertainties and correlations. The paper acknowledges this, but the abstract's first comparison claims agreement with lattice QCD for all three models up to the crossover. To support that claim quantitatively, the c3/c4 comparison should either use direct lattice determinations of the higher temperature derivatives or include a proper covariance propagation; otherwise the agreement claim should be explicitly restricted to c2.
minor comments (5)
- [Eqs. (2)-(6)] The power notation in the rendered equations is difficult to parse. Use explicit parentheses, e.g., (chi_2)^5, so that expressions like chi_2^5 are unambiguous.
- [Sec. 2] Please clarify whether chi_2 is per unit volume and what volume is implicit in c2=1/chi_2, so that the relation <(Delta T)^2>=T^2/C_V is unambiguous.
- [Table 1] State the fitted domain of the Andronic parametrization. The 2.4 GeV row is described as a continuation, but if the original fit includes this energy, the terminology should be adjusted.
- [Appendix A.2] The grey curves labeled 'mu_B=0 LQCD-EoS fits mapped by T_fo' need a clearer explanation: are they the mu_B=0 lattice curves evaluated at the freeze-out temperature, and if so, why is this a meaningful high-energy reference?
- [Fig. 4(c)] The dashed c3=0 contour is visually difficult to distinguish from the color-map features. Consider increasing the line width or using a contrasting color.
Circularity Check
No significant circularity: the cumulant fingerprints are emergent outputs of independently calibrated EoS models.
full rationale
The derivation chain is self-contained and does not recycle its inputs as predictions. The cumulant definitions in Eqs. (2)-(4) are taken from Ref. [25], but the paper evaluates them from the pressures of three HRG models; no temperature-cumulant data are used to fix any parameter. The EV repulsion b is constrained by lattice-QCD conserved-charge fluctuations (Refs. [28,29]) and the VDW parameters a,b by the nuclear ground state (Refs. [28,30], Eq. (9)); these are independent of c2,c3,c4. The lattice comparisons in Fig. 2 use published HotQCD specific-heat points and analytic EoS parametrizations as external benchmarks, not as fit targets. The c2 minimum, c3 sign change, and low-energy freeze-out separation in Figs. 3-4 are emergent consequences of the calibrated VDW/EV equations of state: the c3=0 contour is a separate derivative calculation and is displaced from the c2 minimum by the explicit T powers in Eq. (2), so it is not a trivial derivative of the minimum. The only self-citation (Ref. [4], a review with author overlap) is contextual and non-load-bearing. The identification of the c2 valley as the Widom line in Fig. 4(c) is a labeling choice for a computed minimum locus, not an input used to derive that minimum. The acknowledged extrapolation of Eq. (10) below sqrt(s_NN)=3 GeV is a validity limitation, not a circular step.
Axiom & Free-Parameter Ledger
free parameters (4)
- EVHRG excluded-volume parameter b =
b = 1.0 fm^3 (central; interval 0.4-1.0 fm^3 considered)
- VDW-HRG attraction a =
a = 329 MeV fm^3
- VDW-HRG repulsion b =
b = 3.42 fm^3
- Chemical freeze-out curve coefficients =
T_fo: 158.4, 2.60, 0.45; mu_B: 1307.5, 0.288 (Eq. 10)
axioms (6)
- domain assumption Temperature-fluctuation cumulants follow from entropy derivatives of w=T s-p, giving Eqs. (4)-(6).
- domain assumption The hadron resonance gas with PDG2020 states up to 2.6 GeV (433 states, charm excluded) approximates the hadronic partition function.
- domain assumption Baryon-only EV/VDW sector structure: mesons stay ideal, baryon and antibaryon sectors are separate, and no meson-baryon or baryon-antibaryon interactions are included.
- domain assumption The van der Waals EoS with a=329 MeV fm^3 and b=3.42 fm^3 reproduces the nuclear ground state and yields a liquid-gas critical point near (19.7, 908) MeV.
- domain assumption Analytic HotQCD/WB parametrizations of the lattice EoS can be differentiated to give c3, c4, c5, c6 at mu_B=0.
- ad hoc to paper Chemical freeze-out follows Eq. (10) with strangeness neutrality and n_Q=0.4 n_B, including below sqrt(s_NN)=3 GeV.
read the original abstract
Temperature fluctuations of the matter produced in relativistic heavy-ion collisions are related to event-by-event fluctuations of the mean transverse momentum and, through the specific heat, to the QCD equation of state. We calculate the temperature-fluctuation cumulants $c_2$, $c_3$, and $c_4$ in the ideal hadron resonance gas (HRG), in an excluded-volume HRG consistent with lattice QCD constraints on baryon repulsion, and in a van der Waals HRG that reproduces the nuclear ground state. At zero baryon density all three models agree with lattice QCD thermodynamics up to the chiral crossover and separate above it. Along the chemical freeze-out curve the models remain close for $\sqrt{s_{NN}} \gtrsim 20$ GeV and separate strongly at lower energies, where baryonic interactions dominate the thermal response. In the van der Waals model the variance $c_2$ develops a minimum along the Widom line of the nuclear liquid-gas transition and $c_3$ changes sign across it. Temperature cumulants thus connect mean-transverse-momentum fluctuations to thermodynamic structures in two regions of the QCD phase diagram.
Figures
Reference graph
Works this paper leans on
-
[1]
M. A. Stephanov, K. Rajagopal, E. V . Shuryak, Signatures of the tricritical point in QCD, Phys. Rev. Lett. 81 (1998) 4816.arXiv:hep-ph/9806219
Pith/arXiv arXiv 1998
-
[2]
M. A. Stephanov, K. Rajagopal, E. V . Shuryak, Event- by-event fluctuations in heavy ion collisions and the QCD critical point, Phys. Rev. D 60 (1999) 114028.arXiv: hep-ph/9903292
Pith/arXiv arXiv 1999
-
[3]
A. Bzdak, S. Esumi, V . Koch, J. Liao, M. Stephanov, N. Xu, Mapping the phases of quantum chromodynam- ics with beam energy scan, Phys. Rept. 853 (2020) 1. arXiv:1906.00936
Pith/arXiv arXiv 2020
-
[4]
A. Pandav, D. Mallick, B. Mohanty, Search for the QCD critical point in high energy nuclear collisions, Prog. Part. Nucl. Phys. 125 (2022) 103960.arXiv:2203.07817
Pith/arXiv arXiv 2022
-
[5]
W.-j. Fu, J. M. Pawlowski, F. Rennecke, QCD phase struc- ture at finite temperature and density, Phys. Rev. D 101 (2020) 054032.arXiv:1909.02991
Pith/arXiv arXiv 2020
-
[6]
Adam, et al., Nonmonotonic energy dependence of net- proton number fluctuations, Phys
J. Adam, et al., Nonmonotonic energy dependence of net- proton number fluctuations, Phys. Rev. Lett. 126 (2021) 092301.arXiv:2001.02852
arXiv 2021
-
[7]
M. S. Abdallah, et al., Measurements of proton high order cumulants in √sNN =3 GeV Au+Au collisions and im- plications for the QCD critical point, Phys. Rev. Lett. 128 (2022) 202303.arXiv:2112.00240
arXiv 2022
-
[8]
B. Aboona, et al., Beam energy dependence of fifth- and sixth-order net-proton number fluctuations in Au+Au col- lisions at RHIC, Phys. Rev. Lett. 130 (2023) 082301. arXiv:2207.09837. 6
arXiv 2023
-
[9]
B. E. Aboona, et al., Precision measurement of net- proton-number fluctuations in Au+Au collisions at RHIC, Phys. Rev. Lett. 135 (2025) 142301.arXiv:2504. 00817
2025
-
[10]
P. Braun-Munzinger, B. Friman, K. Redlich, A. Rustamov, J. Stachel, Relativistic nuclear collisions: Establishing a non-critical baseline for fluctuation measurements, Nucl. Phys. A 1008 (2021) 122141.arXiv:2007.02463
Pith/arXiv arXiv 2021
-
[11]
L. D. Landau, E. M. Lifshitz, Statistical Physics, Part 1, Pergamon Press, Oxford, 1980
1980
-
[12]
Stodolsky, Temperature fluctuations in multiparticle production, Phys
L. Stodolsky, Temperature fluctuations in multiparticle production, Phys. Rev. Lett. 75 (1995) 1044
1995
-
[13]
Gavin, Traces of thermalization from transverse mo- mentum fluctuations in nuclear collisions, Phys
S. Gavin, Traces of thermalization from transverse mo- mentum fluctuations in nuclear collisions, Phys. Rev. Lett. 92 (2004) 162301.arXiv:nucl-th/0308067
Pith/arXiv arXiv 2004
-
[14]
R. Korus, S. Mrowczynski, M. Rybczynski, Z. Wlodar- czyk, Transverse momentum fluctuations due to tempera- ture variation in high-energy nuclear collisions, Phys. Rev. C 64 (2001) 054908.arXiv:nucl-th/0106041
Pith/arXiv arXiv 2001
-
[15]
S. Acharya, et al., Skewness and kurtosis of mean trans- verse momentum fluctuations at the LHC energies, Phys. Lett. B 850 (2024) 138541.arXiv:2308.16217
Pith/arXiv arXiv 2024
-
[16]
G. Aad, et al., Disentangling sources of momentum fluc- tuations in Xe+Xe and Pb+Pb collisions with the ATLAS detector, Phys. Rev. Lett. 133 (2024) 252301.arXiv: 2407.06413
Pith/arXiv arXiv 2024
-
[17]
Adams, et al., Incident energy dependence of pt corre- lations at RHIC, Phys
J. Adams, et al., Incident energy dependence of pt corre- lations at RHIC, Phys. Rev. C 72 (2005) 044902.arXiv: nucl-ex/0504031
Pith/arXiv arXiv 2005
-
[18]
S. S. Adler, et al., Measurement of nonrandom event- by-event fluctuations of average transverse momentum in√sNN =200 GeV Au+Au and p+p collisions, Phys. Rev. Lett. 93 (2004) 092301.arXiv:nucl-ex/0310005
Pith/arXiv arXiv 2004
-
[19]
T. Anticic, et al., Energy dependence of transverse mo- mentum fluctuations in Pb+Pb collisions at the CERN SPS, Phys. Rev. C 79 (2009) 044904.arXiv:0810.5580
Pith/arXiv arXiv 2009
-
[20]
D. Adamova, et al., Event-by-event fluctuations of the mean transverse momentum in 40, 80 and 158 a GeV/c Pb-Au collisions, Nucl. Phys. A 727 (2003) 97.arXiv: nucl-ex/0305002
Pith/arXiv arXiv 2003
-
[21]
S. Basu, S. Chatterjee, R. Chatterjee, T. K. Nayak, B. K. Nandi, Specific heat of matter formed in relativistic nu- clear collisions, Phys. Rev. C 94 (2016) 044901.arXiv: 1601.05631
Pith/arXiv arXiv 2016
-
[22]
B. E. Aboonaet al.[STAR Collaboration], Nonmono- tonicity of transverse momentum correlations in Au+Au collisions at RHIC, Phys. Rev. Lett. (in press) (2026). arXiv:2604.06434,doi:10.1103/2xsn-rgx3
Pith/arXiv arXiv 2026
-
[23]
F. G. Gardim, G. Giacalone, J.-Y . Ollitrault, The mean transverse momentum of ultracentral heavy-ion collisions: A new probe of hydrodynamics, Phys. Lett. B 809 (2020) 135749.arXiv:1909.11609
Pith/arXiv arXiv 2020
-
[24]
L. Zhang, J. Chen, C. Zhang, Energy dependence of trans- verse momentum fluctuations in Au+Au collisions from a multiphase transport model, Phys. Rev. C 111 (2025) 024911.arXiv:2501.08209
Pith/arXiv arXiv 2025
-
[25]
J. Chen, W.-j. Fu, S. Yin, C. Zhang, High-order fluctuations of temperature in hot QCD matter, arXiv:2504.06886 [hep-ph] (2025).arXiv:2504.06886
Pith/arXiv arXiv 2025
-
[26]
H. Liu, P. Wu, H.-M. Liu, P.-C. Chu, Fluctuation of temperature in the Polyakov-loop extended Nambu–Jona- Lasinio model, Universe 12 (2026) 37.arXiv:2512. 23934
2026
-
[27]
D. H. Rischke, M. I. Gorenstein, H. Stoecker, W. Greiner, Excluded volume effect for the nuclear matter equation of state, Z. Phys. C 51 (1991) 485
1991
-
[28]
V . V ovchenko, M. I. Gorenstein, H. Stoecker, van der waals interactions in hadron resonance gas: From nu- clear matter to lattice QCD, Phys. Rev. Lett. 118 (2017) 182301.arXiv:1609.03975
Pith/arXiv arXiv 2017
-
[29]
J. M. Karthein, V . Koch, C. Ratti, V . V ovchenko, Con- straining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges, Phys. Rev. D 104 (2021) 094009.arXiv:2107. 00588,doi:10.1103/PhysRevD.104.094009
-
[30]
V . V ovchenko, D. V . Anchishkin, M. I. Gorenstein, van der waals equation of state with Fermi statistics for nu- clear matter, J. Phys. A 48 (2015) 305001.arXiv: 1501.03785
Pith/arXiv arXiv 2015
-
[31]
S. Samanta, B. Mohanty, Criticality in a hadron resonance gas model with the van der waals interaction, Phys. Rev. C 97 (2018) 015201.arXiv:1709.04446
Pith/arXiv arXiv 2018
-
[32]
Bazavov, et al., Equation of state in (2+1)-flavor QCD, Phys
A. Bazavov, et al., Equation of state in (2+1)-flavor QCD, Phys. Rev. D 90 (2014) 094503.arXiv:1407.6387
Pith/arXiv arXiv 2014
-
[33]
S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, K. K. Szabo, Full result for the QCD equation of state with 2+1 flavors, Phys. Lett. B 730 (2014) 99.arXiv: 1309.5258
Pith/arXiv arXiv 2014
-
[34]
S. Borsanyi, Z. Fodor, J. N. Guenther, P. Parotto, A. Pasz- tor, C. Ratti, V . V ovchenko, C. H. Wong, Lattice QCD constraints on the critical point from an improved preci- sion equation of state, arXiv:2502.10267 [hep-lat] (2025). arXiv:2502.10267
Pith/arXiv arXiv 2025
-
[35]
A. Andronic, P. Braun-Munzinger, K. Redlich, J. Stachel, Decoding the phase structure of QCD via particle pro- duction at high energy, Nature 561 (2018) 321.arXiv: 1710.09425. 7 Appendix A. Particle-sector diagnostics and higher orders Appendix A.1. Meson and baryon sectors To identify the origin of the model hierarchy, we evaluate the cumulants from sect...
Pith/arXiv arXiv 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.