Pith. sign in

REVIEW 3 major objections 6 minor 24 references

Excited cluster states: A new source for proton number fluctuations in the high baryon density regime

T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Excited nuclear cluster decays shift proton number fluctuation ratios by up to 100 percent at low collision energies.

desk verdict A useful new baseline for CBM-era proton cumulants, but the headline corrections rest on an unverified at-most-one-proton assumption that the cited decay tables likely contradict. read the letter →

arxiv 2412.04994 v1 pith:GIZ7GSAL submitted 2024-12-06 nucl-th hep-phnucl-ex

Boris Tomasik , Marcus Bleicher This is my paper · ORCID
classification nucl-thhep-phnucl-ex PACS 25.75.-q24.60.-k
keywords protonnumberfluctuationsQCDphasediagramrelativisticheavy-ioncollisionsexcitednuclearclustershadronresonancegashigher-ordercumulantsCBM
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper identifies a previously neglected source of event-by-event proton number fluctuations in heavy-ion collisions: the decay products of excited nuclear cluster states. Using a statistical (hadron-resonance-gas) model, the authors include A=5–7 excited clusters as proton-feeding resonances and compute the first six cumulant ratios of the net-proton distribution for collision energies $\sqrt{s_{NN}}$ from 2 to 5 GeV. They find that the clusters change every cumulant ratio, from about 1% for the variance-to-mean ratio to roughly 100% for the sixth-to-second ratio at the lowest energy. The correction is substantial below $\sqrt{s_{NN}} \approx 3.5$ GeV and falls away above about 5 GeV, moving each ratio closer to unity. The authors' point is that planned CBM measurements of high-order proton cumulants can be compared with lattice QCD and effective-model predictions only if this feeddown is included in the baseline.

What carries the argument

The binomial feeddown cumulant machinery: for each resonance species $R$, the proton-number cumulant-generating function is built from the probability $p_R$ that a decay produces one proton, with direct protons treated as $p_R = 1$. Derivatives of this function give the first six proton cumulants as sums over all resonance species—including the excited clusters—of terms in $p_R$ and the species' grand-canonical number cumulants. The excited clusters enter through the same binomial formulas, and their decay probabilities and masses are taken from the imported light-nucleus tables. Because the cumulant ratios are volume-independent, the calculation needs only the chemical freeze-out temperature $T$ and baryochemical potential $\mu_B$, which are fixed by a parametrization of $\sqrt{s_{NN}}$; this machinery is what converts the decay tables into the reported corrections.

What would settle it

Check the A=5–7 decay tables used in the paper for a state with a two-proton branch—6Be decaying to an alpha plus two protons would suffice—and if such a branch exists, recompute the sixth-to-second cumulant ratio with the true multinomial decay probabilities; the result would show whether the reported 100% correction survives.

Watch

Extended reading notes

Core claim

The paper claims that the standard hadron-resonance-gas baseline for proton cumulant ratios has been missing an entire class of sources: the decay products of excited nuclear cluster states with A=5–7. Treating these clusters exactly like resonances that decay into at most one proton, and using the binomial cumulant formulas derived from the cumulant-generating function, the authors calculate the first six cumulant ratios of the (net-)proton number at chemical freeze-out. Including the clusters raises every ratio toward unity; at the lowest inspected energy the sixth-to-second cumulant ratio nearly doubles, while the variance-to-mean ratio moves by only about one percent. The effect is visible for $\sqrt{s_{NN}}$ below about 3.5 GeV and becomes negligible above about 5 GeV. The authors frame the result as a necessary improvement to baseline calculations for CBM.

Load-bearing premise

The calculation assumes that every excited nuclear cluster included in the decay tables produces at most one proton; the paper does not verify this channel-by-channel, and if any included A=5–7 state has a two-proton decay branch the binomial cumulant formulas and the reported correction sizes would need to be recomputed.

Editorial extensions

If this is right

  • At $\sqrt{s_{NN}}$ below about 3.5 GeV, hadron-resonance-gas baselines that omit excited cluster feeddown understate the sixth-to-second-order proton cumulant ratio by up to a factor of two.
  • All the cumulant ratios examined move closer to unity when excited clusters are included, with the higher-order ratios affected most strongly.
  • Above about 5 GeV the correction is negligible, so high-energy fluctuation baselines do not need to be changed.
  • The scaling variable $r_{ij}$ collapses the corrected ratios onto two energy-dependent curves, meaning the correction can be summarized as a single factor applied to the distance of each ratio from unity—up to one sixth at 2.4 GeV.
  • Direct comparisons of CBM proton-fluctuation data with lattice QCD or effective-model results at FAIR energies require this feeddown to be part of the model baseline.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The reported size of the correction is sensitive to the assumption that every excited cluster emits at most one proton; if the imported decay tables contain a two-proton channel such as 6Be decaying to an alpha plus two protons, the binomial formulas would need to be replaced by multinomial ones and the 100% sixth-order effect could change.
  • Because all the $r_{ij}$ curves collapse onto two universal curves, one could parametrize the cluster correction for any cumulant ratio as a function of (1 minus the ratio) and apply it to existing hadron-resonance-gas calculations without recomputing the full thermal model.
  • The same feeddown logic should carry over to other conserved-charge cumulants, so excited cluster decays may also affect net-baryon, electric-charge, or strangeness fluctuation ratios at low beam energies—an extension the paper does not pursue.
  • A direct data-level test: if high-order proton cumulant ratios in central Au+Au at the lowest CBM energies rise toward unity in a way standard hadron-resonance-gas models do not reproduce, that rise would match the cluster-feeddown mechanism described here.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper estimates the contribution of decays of excited nuclear cluster states to event-by-event proton number cumulant ratios in the statistical model, for sqrt(s_NN) between 2 and 5 GeV. The formalism treats each cluster as a resonance that produces at most one proton, uses the binomial cumulant formulas (12a)-(12f) from [23], and takes the cluster input from [19]. It reports corrections from about 1% for the second-to-first cumulant ratio to about 100% for the sixth-to-second cumulant ratio at the lowest energies, with the effect becoming negligible above 3.5 GeV. The authors argue that such corrections provide an improved baseline for CBM comparisons with lattice QCD and effective models.

Significance. If correct, this is a useful quantitative baseline improvement for proton cumulant measurements in the high-baryon-density regime. The calculation is non-circular: the freeze-out parameters come from an independent fit to hadron multiplicities [24], the cluster input comes from external nuclear data evaluations via [19], and no fluctuation observable is fitted. The cumulant recursion is transparent, the volume cancels in the reported ratios, and the final ratios constitute a falsifiable prediction for CBM. The main risk is technical rather than conceptual: the binomial treatment rests on an unverified at-most-one-proton assumption, and the p_R values and state selection are not documented. The significance is moderate: this is an incremental but potentially important baseline correction, not a new phenomenon or a new formalism.

major comments (3)
  1. [Sec. 2, Eqs. (11)-(12)] The cumulant-generating function (11) and the explicit formulas (12a)-(12f) rely on the statement that 'all included excited clusters produce at most one proton.' The manuscript does not verify this against its adopted input: the A=5-7 evaluation [21] lists the 6Be ground state with a dominant alpha+p+p decay channel, and the text does not state that such states are excluded from Tables I/II of [19]. If any included state has a two-proton decay channel, p_R is not the probability of producing exactly one proton, the binomial factorial-cumulant structure of (12) is invalid for that species, and the reported corrections—including the roughly 100% change in kappa_H sigma^4 at low energy—must be recomputed with a compound or multinomial decay distribution. This is load-bearing because the central quantitative claim follows directly from these formulas.
  2. [Sec. 3, cluster input and p_R] The paper gives no explicit list of the excited cluster states included, no branching ratios, and no prescription for converting the decay tables in [19]/[21] into the p_R values entering Eqs. (12). Without this information the calculation is not reproducible, and one cannot assess how strongly the headline 100% correction depends on the inclusion or exclusion of specific states. The authors should provide the state list with adopted p_R values and state explicitly how multi-body decay channels are treated.
  3. [Sec. 2, Eq. (14)] Eq. (14) assumes resonances and clusters have vanishing widths, yet the excited cluster states taken from [21] often have substantial widths. The finite-width effect on the mean and higher cumulants in Eqs. (15)-(16) is not estimated. Since the paper highlights corrections at the 1% to 100% level, a quantitative argument that finite widths change the reported ratios by less than, say, the claimed effect size is needed before the numbers can be considered robust.
minor comments (6)
  1. [Sec. 3, Eq. (19)] For i=1 the definition r_i1 is undefined, since the denominator becomes 1 - <N>_c/<N>_c = 0; the text should specify that Eq. (19) applies for i>=2 or otherwise define the i=1 case.
  2. [Fig. 3] The caption states that the curves collapse to 'two curves' but does not identify which curve corresponds to ratios with <N>_c in the denominator and which to ratios with <(Delta N)^2>_c in the denominator; add a legend or explanatory note.
  3. [Sec. 2, text above Eq. (11)] The symbol 'NR' appears as 'the number of protons resulting from decays of NR resonances'; this should be typeset as N_R (resonance multiplicity) to avoid confusion with the resonance species label R.
  4. [Sec. 2, Eqs. (6)-(7)] The terms 'hyperskewness' and 'hyperkurtosis' are not standard in heavy-ion literature; please define them or cite standard usage, as some readers will confuse the fifth and sixth standardized cumulants with skewness and kurtosis of hyperdistributions.
  5. [Sec. 3, text after Fig. 3] The phrase 'cuts up to one sixth from the distance of the cumulant ratio from unity' is unclear; please rephrase to state precisely what the scaling variable r_ij measures and what the 'one sixth' value implies quantitatively.
  6. [References] Reference [19] is incomplete: the volume/page is missing ('Phys. Lett. B (2020) 135746'); please provide the full citation.

Circularity Check

0 steps flagged · score 1.0 of 10

No load-bearing circularity: the reported cumulant corrections follow from published decay tables and externally fitted freeze-out parameters, with no parameter tuned to the reported fluctuation ratios.

full rationale

The derivation chain is self-contained. The input quantities are (i) the freeze-out T and baryochemical potential parametrization adopted from Vovchenko et al. [24], which was fitted to hadron multiplicities rather than to proton cumulants; (ii) the list of stable and excited clusters from Tables I and II of Vovchenko et al. [19]; and (iii) the standard binomial compound cumulant-generating function, Eq. (11), whose resulting formulas (12a)-(12f) are displayed in the paper rather than imported as a black box. No parameter appearing in the reported ratios is fitted to those ratios, so the 1% to 100% corrections are a direct model consequence, not a re-fit in disguise. The self-citation to [23] for the formalism is not load-bearing because the relevant equations are reproduced in the present text. One non-circular caveat: Section 2's assertion that all included excited clusters produce at most one proton is not checked against the decay tables of [21]; if a state such as 6Be can decay to alpha + p + p, the binomial p_R structure would be an input assumption rather than a validated fact. That is a correctness and validity risk, not a circularity.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

No new particles, forces, or conserved quantities are introduced. The model inputs are existing nuclear states and a standard HRG framework. The free parameters are the externally fitted freeze-out curve and the branching probabilities imported from prior compilations.

free parameters (2)
  • chemical freeze-out parameters a, b, c, d, e = a=0.157 GeV, b=0.087 GeV^-1, c=0.092 GeV^-3, d=1.477 GeV, e=0.343 GeV^-1
    Taken from Vovchenko et al. [24]; fitted to hadron multiplicities, not to proton cumulants. They set T(sqrt(s_NN)) and mu_B(sqrt(s_NN)) used in all cumulant evaluations.
  • decay-to-proton probabilities p_R for each excited cluster = not tabulated here; taken as branching ratios in Tables I/II of [19]
    The central cumulant formulas depend on p_R; their values are imported from [19] and not listed here, so numerical reproduction requires the external tables.
assumptions (4)
  • domain assumption Grand-canonical hadron resonance gas with zero-width resonances describes the freeze-out state at sqrt(s_NN)=2-5 GeV.
    Invoked in Section 2 (Eq. 14) and stated in the conclusion; the model is the basis of all cumulants.
  • domain assumption Chemical equilibrium with T and mu_B from the parametrization of [24].
    Eqs. (17)-(18) and Table 1 set the freeze-out conditions; the fit was made to hadron multiplicities, not to fluctuations.
  • ad hoc to paper Every included excited cluster decays into at most one proton.
    Section 2 after Eq. (11); justifies binomial cumulant formulas. Not verified against the decay tables.
  • domain assumption Proton and antiproton numbers are independent so net-proton cumulants add as in Eq. (1).
    Uses independent GCE species; no cross-correlations are considered.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Excited cluster states: A new source for proton number fluctuations in the high baryon density regime." pith.science (2026). https://pith.science/paper/GIZ7GSAL

@misc{pith2026241204994,
  author       = {Pith},
  title        = {Pith review of: Excited cluster states: A new source for proton number fluctuations in the high baryon density regime},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GIZ7GSAL}},
  note         = {Machine review of arXiv:2412.04994}
}
abstract

We calculate the contribution of the decay products of excited nuclear cluster states to the event-by-event fluctuations of protons in the energy range from $\sqrt{s_{NN}}=2-5$~GeV within the statistical model. We find that the inclusion of the excited nuclear clusters yields corrections to all cumulant ratios, ranging from 1\% for ratio of second to first-order cumulant to 100\% for the sixth to second order cumulant towards the lowest inspected energy. As expected the contribution of excited cluster states is most important at low energies $\sqrt{s_{NN}}<3.5$~GeV and becomes negligible at higher collision energies. Especially in light of the expected ultra-high precision data from CBM at FAIR, this new contribution is important to allow for a quantitative comparison with (potentially later available) lattice QCD or effective model results.

Figures

Figures reproduced from arXiv: 2412.04994 by the authors.

Figure 1
Figure 1. Dependence of volume-independent ratios of the first four cumulants of the net-proton number distribution on the collision energy. The relation to the [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Dependence of volume-independent ratios that include the fifth and sixth cumulants of the net-proton number distribution on the collision energy. The [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. The scaling variable ri j, defined in Eqs. (19) and (20) as depending on collision energy, for various ratios of the cumulants. 4. Conclusion We have calculated the contribution of the decay products of excited nuclear cluster states to the event-by-event fluctua￾tions of protons in the energy range from √ sNN = 2 − 5 GeV. To this aim we assumed validity of the statistical model, with parameters adjusted to describe… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 8 canonical work pages

  1. [19]

    V ovchenko, B

    V . V ovchenko, B. Dönigus, B. Kardan, M. Lorenz, H. Stoecker, Feed- down contributions from unstable nuclei in relativistic heavy-ion col- lisions, Phys. Lett. B (2020) 135746. arXiv:2004.04411, doi:10. 1016/j.physletb.2020.135746

    arXiv 2020

  2. [21]

    D. R. Tilley, C. M. Cheves, J. L. Godwin, G. M. Hale, H. M. Hofmann, J. H. Kelley, C. G. Sheu, H. R. Weller, Energy levels of light nuclei A=5, A =6, A =7, Nucl. Phys. A 708 (2002) 3–163. doi:10.1016/ S0375-9474(02)00597-3

    work page 2002

  3. [23]

    Tomasik, P

    B. Tomasik, P. Hillmann, M. Bleicher, Proton number fluctuations in par- tial chemical equilibrium, Phys. Rev. C 104 (4) (2021) 044907. arXiv: 2107.03830, doi:10.1103/PhysRevC.104.044907

    arXiv 2021

  4. [24]

    V ovchenko, V

    V . V ovchenko, V . V . Begun, M. I. Gorenstein, Hadron multiplicities and chemical freeze-out conditions in proton-proton and nucleus-nucleus collisions, Phys. Rev. C 93 (6) (2016) 064906. arXiv:1512.08025, doi:10.1103/PhysRevC.93.064906. 5

    arXiv 2016

  5. [1]

    B. P. Abbott, et al., GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral, Phys. Rev. Lett. 119 (16) (2017) 161101. arXiv:1710.05832, doi:10.1103/PhysRevLett.119.161101

    arXiv 2017

  6. [2]

    Sorensen, et al., Dense nuclear matter equation of state from heavy- ion collisions, Prog

    A. Sorensen, et al., Dense nuclear matter equation of state from heavy- ion collisions, Prog. Part. Nucl. Phys. 134 (2024) 104080. arXiv:2301. 13253, doi:10.1016/j.ppnp.2023.104080

    arXiv 2024

  7. [3]

    Odyniec, Beam Energy Scan Program at RHIC (BES I and BES II) – Probing QCD Phase Diagram with Heavy-Ion Collisions, PoS CORFU2018 (2019) 151

    G. Odyniec, Beam Energy Scan Program at RHIC (BES I and BES II) – Probing QCD Phase Diagram with Heavy-Ion Collisions, PoS CORFU2018 (2019) 151. doi:10.22323/1.347.0151

    work page doi:10.22323/1.347.0151 2019

  8. [4]

    Stephanov, QCD critical point: recent developments, in: 11th In- ternational Worksop on QCD - Theory and Experiment, 2024

    M. Stephanov, QCD critical point: recent developments, in: 11th In- ternational Worksop on QCD - Theory and Experiment, 2024. arXiv: 2410.02861

    arXiv 2024

Show all 24 references
  1. [5]

    M. A. Stephanov, K. Rajagopal, E. V . Shuryak, Signatures of the tricritical point in QCD, Phys. Rev. Lett. 81 (1998) 4816–4819. arXiv:hep-ph/ 9806219, doi:10.1103/PhysRevLett.81.4816

    1998 doi

  2. [6]

    R. A. Lacey, Indications for a Critical End Point in the Phase Diagram for Hot and Dense Nuclear Matter, Phys. Rev. Lett. 114 (14) (2015) 142301. arXiv:1411.7931, doi:10.1103/PhysRevLett.114.142301

    2015 arXiv

  3. [7]

    Adhikary, Search for the QCD critical point by NA61 /SHINE at the CERN SPS, in: 10th International Conference on New Frontiers in Physics, 2023

    H. Adhikary, Search for the QCD critical point by NA61 /SHINE at the CERN SPS, in: 10th International Conference on New Frontiers in Physics, 2023. arXiv:2308.04254

    2023 arXiv

  4. [8]

    Herold, M

    C. Herold, M. Nahrgang, I. Mishustin, M. Bleicher, Formation of droplets with high baryon density at the QCD phase transition in expanding matter, Nucl. Phys. A 925 (2014) 14–24.arXiv:1304.5372, doi:10.1016/j. nuclphysa.2014.01.010

    2014 arXiv

  5. [9]

    M. A. Stephanov, Non-Gaussian fluctuations near the QCD critical point, Phys. Rev. Lett. 102 (2009) 032301. arXiv:0809.3450, doi:10. 1103/PhysRevLett.102.032301

    2009 arXiv

  6. [10]

    Bazavov, et al., Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model, Phys

    A. Bazavov, et al., Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model, Phys. Rev. D 86 (2012) 034509. arXiv:1203.0784, doi:10.1103/PhysRevD.86.034509

    2012 arXiv

  7. [11]

    Nahrgang, T

    M. Nahrgang, T. Schuster, M. Mitrovski, R. Stock, M. Bleicher, Net- baryon-, net-proton-, and net-charge kurtosis in heavy-ion collisions within a relativistic transport approach, Eur. Phys. J. C 72 (2012) 2143. arXiv:0903.2911, doi:10.1140/epjc/s10052-012-2143-6

    2012 arXiv

  8. [12]

    Braun-Munzinger, K

    P. Braun-Munzinger, K. Redlich, A. Rustamov, J. Stachel, The imprint of conservation laws on correlated particle production, JHEP 08 (2024) 113. arXiv:2312.15534, doi:10.1007/JHEP08(2024)113

    2024 arXiv

  9. [13]

    V . A. Kuznietsov, M. I. Gorenstein, V . Koch, V . V ovchenko, Coordinate versus momentum cuts and e ffects of collective flow on critical fluc- tuations, Phys. Rev. C 110 (1) (2024) 015206. arXiv:2404.00476, doi:10.1103/PhysRevC.110.015206

    2024 arXiv

  10. [14]

    Nahrgang, S

    M. Nahrgang, S. Leupold, C. Herold, M. Bleicher, Nonequilibrium chiral fluid dynamics including dissipation and noise, Phys. Rev. C 84 (2011) 024912. arXiv:1105.0622, doi:10.1103/PhysRevC.84.024912

    2011 arXiv

  11. [15]

    Skokov, B

    V . Skokov, B. Friman, K. Redlich, V olume Fluctuations and Higher Order Cumulants of the Net Baryon Number, Phys. Rev. C 88 (2013) 034911. arXiv:1205.4756, doi:10.1103/PhysRevC.88.034911

    2013 arXiv

  12. [16]

    Holzmann, V

    R. Holzmann, V . Koch, A. Rustamov, J. Stroth, Controlling volume fluc- tuations for studies of critical phenomena in nuclear collisions, Nucl. Phys. A 1050 (2024) 122924. arXiv:2403.03598, doi:10.1016/j. nuclphysa.2024.122924

    2024 arXiv

  13. [17]

    Sombun, J

    S. Sombun, J. Steinheimer, C. Herold, A. Limphirat, Y . Yan, M. Bleicher, Higher order net-proton number cumulants dependence on the centrality definition and other spurious e ffects, J. Phys. G 45 (2) (2018) 025101. arXiv:1709.00879, doi:10.1088/1361-6471/aa9c6c

    2018 arXiv

  14. [18]

    Fecková, J

    Z. Fecková, J. Steinheimer, B. Tomášik, M. Bleicher, Net-proton num- ber kurtosis and skewness in nuclear collisions: Influence of deuteron formation, Phys. Rev. C 92 (6) (2015) 064908. arXiv:1510.05519, doi:10.1103/PhysRevC.92.064908

    2015 arXiv

  15. [20]

    D. R. Tilley, H. R. Weller, G. M. Hale, Energy levels of light nuclei A = 4, Nucl. Phys. A 541 (1992) 1–104. doi:10.1016/0375-9474(92) 90635-W

    1992 doi

  16. [22]

    Nahrgang, M

    M. Nahrgang, M. Bluhm, P. Alba, R. Bellwied, C. Ratti, Impact of reso- nance regeneration and decay on the net-proton fluctuations in a hadron resonance gas, Eur. Phys. J. C 75 (12) (2015) 573. arXiv:1402.1238, doi:10.1140/epjc/s10052-015-3775-0

    2015 arXiv

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.