Pith. sign in

REVIEW 2 major objections 5 minor 1 cited by

The paper claims that the standard centrality bin width correction, applied to proton-number cumulants, overcorrects and suppresses genuine fluctuation signals whenever baryon resonance decays correlate proton number with charged multiplici

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 →

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

T0 review · deepseek-v4-flash

2026-08-03 22:08 UTC pith:VZ2ZRY5O

load-bearing objection Useful analytic caution about CBWC overcorrection, but the quantitative plots are for a full-multiplicity class, not the finite centrality bins STAR uses. the 2 major comments →

arxiv 2511.11869 v2 pith:VZ2ZRY5O submitted 2025-11-14 nucl-th

To bin or not to bin: does binning in multiplicity reliably suppress unwanted volume fluctuations?

Bengt Friman , Volker Koch This is my paper
classification nucl-th
keywords centrality bin width correctionvolume fluctuationsnet-proton cumulantsbaryon resonancesbivariate Poisson distributionheavy-ion collisionsfactorial cumulants
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Centrality bin width correction (CBWC) is the standard procedure for removing volume fluctuations from measured cumulants of (net-)proton number in heavy-ion collisions. The paper shows within an analytically solvable model that CBWC does not just remove volume fluctuations: it also subtracts a spurious contribution whenever the proton number and the charged-particle multiplicity are correlated by resonance decays. When the resonance correlation strength λ exceeds roughly the mean number of primordial protons, the corrected cumulants fall below the true fixed-volume cumulants—for the variance, κ₂^CBWC = λ + ν + ν²/(k+λ+μ) − λ²/(λ+μ) < λ + ν = κ₂^true—and the deficit grows for higher-order cumulants. This matters because these cumulants are the observables used to search for critical-point structure in the QCD phase diagram; a corrected signal that is biased low could mimic or mask the physics of interest.

Core claim

The central result is that CBWC, which computes corrected cumulants as weighted averages of per-bin cumulants, preserves the moments but not the cumulants of the particle-number distribution. In the model—a bivariate Poisson distribution for protons and charged pions (primordial means ν and μ, plus resonance decays generating λ correlated pairs) folded over a Gamma distribution of the reduced volume with shape parameter k—the corrected cumulants acquire extra terms proportional to λ²/(λ+μ) and ν²/(k+λ+μ). For λ ≳ ν the negative λ term dominates, so κₙ^CBWC is smaller than the true fixed-volume cumulant κₙ^true = λ + ν. For net protons, the procedure additionally introduces a negative covaria

What carries the argument

The centrality bin width correction (CBWC): κₙ^CBWC = Σᵢ pᵢ κₙ,ᵢ, averaging per-multiplicity-bin cumulants weighted by the bin's event fraction. The model couples a bivariate Poisson distribution (primordial protons/pions with means ν, μ plus resonance-decay correlations with strength λ) to a Gamma distribution of the reduced volume (shape k). The generating-function analysis yields closed-form cumulants, and the cancellation of volume fluctuations is controlled by the balance between the terms ν²/(k+λ+μ) and λ²/(λ+μ); when the second exceeds the first, CBWC overcorrects.

Load-bearing premise

The analytical CBWC cumulants are computed for a centrality class that sums over the entire multiplicity distribution—the bin weights are the full-distribution probabilities—so the quantitative predictions do not yet account for the finite multiplicity intervals used in real analyses, and truncating to those intervals could alter the cancellation terms that drive the overcorrection.

What would settle it

Repeat the model calculation with the multiplicity sums restricted to a finite centrality window [M1, M2] and the bin weights renormalized within that window; if the corrected cumulants no longer fall below the true fixed-volume cumulants for large λ, the overcorrection is an artifact of integrating over all multiplicities.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • CBWC-corrected cumulants underestimate the true fixed-volume cumulants whenever the resonance correlation strength λ is comparable to or larger than the primordial proton number ν; the model puts the crossover near λ ≈ ν.
  • The suppression is stronger for higher-order cumulants, so κ₃ and κ₄—the observables most sensitive to critical behavior—are the most vulnerable to being artificially reduced.
  • For net-proton cumulants, CBWC induces a negative proton–antiproton covariance and, in the strong-correlation limit, can flip the sign of κ₃ or drive κ₄ negative relative to the true positive values.
  • The same analysis shows CBWC converges to the true cumulants when correlations are absent (λ=0) or when primordial pions dominate the multiplicity (μ ≫ λ, ν), so the failure is specific to resonance-dominated correlations.
  • Because the authors verify numerically that bin widths around ΔM ≈ 10 converge to the unit-width analytic results, the qualitative overcorrection conclusion is robust against realistic binning.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The paper's analytic cumulants integrate the multiplicity over all values, so the quantitative ratios and the location of the crossover (λ ≈ ν) apply to a centrality class spanning the entire multiplicity range, not to the narrow 0–5%/5–10% windows used in experiments; redoing the sum over a truncated interval with renormalized pⱼ could shift the curves.
  • An experimental test is within reach: measure the ratio of CBWC-corrected to uncorrected cumulants for proton or net-proton numbers as a function of collision energy; if the thermal-model estimate of resonance-to-primordial proton production follows the paper's scenario, the ratio should dip below one at the energies where λ ≈ ν.
  • The factorial-cumulant relation F₂^CBWC = ν²/(k+λ+μ) − λ²/(λ+μ) suggests a targeted observable: a negative factorial cumulant of the proton distribution at large λ would be a signature of overcorrection, distinguishing it from genuine critical fluctuations which typically enter with positive factorial cumulants.
  • A natural extension, given the model's lack of critical dynamics, is to add a genuine critical contribution and compute how much of it CBWC would wash out; the generating-function machinery developed here makes that contamination calculable.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. This paper studies the Centrality Bin Width Correction (CBWC) in a solvable stochastic model. The volume is Gamma-distributed; at fixed volume, proton number and pion multiplicity follow a bivariate Poisson distribution, with an extra correlated component modeling Delta-resonance decays into a proton and a pion. The authors derive analytic expressions for the 'true' fixed-volume cumulants, the volume-fluctuating cumulants, and the CBWC-corrected cumulants for protons and net protons. They find that CBWC removes volume fluctuations when the resonance component is absent, but for sufficiently strong resonance correlations it overcorrects and yields cumulants below the true values. The paper argues this implies CBWC can suppress the physics of interest in heavy-ion experiments.

Significance. The question is timely and consequential because CBWC is the standard correction used by the STAR collaboration. The model is transparent, with explicit generating-function derivations in appendices and several nontrivial limits correctly reproduced (e.g., the Skellam limit for net protons at mu_B=0, the N proportional to M limit, and the k->infinity limit). The paper does not fit any parameter to the result, and the algebra can be checked by the reader. If the central claim were established for the finite centrality classes used in experiments, this would be an important caveat. At present the quantitative evidence is computed for the full multiplicity distribution, not for a finite [M1,M2] class, which weakens the conclusion's applicability.

major comments (2)
  1. [§IV and Appendix C, Eqs. (33), (44)-(45)] The CBWC cumulants are evaluated with the bin weight p_j = P(M_j) from the full multiplicity distribution and sums over all M_j. This is not the CBWC of a finite centrality class [M1,M2] defined in §II, where p_j must be renormalized to the class. Because kappa_{n,j} in Eq. (43) is linear in M_j, the class CBWC differs from Eq. (45) by terms proportional to <M>_class - (lambda+mu). For the top-left Fig. 2 parameters (nu=21, mu=157, k=83, lambda/nu=1.4), a 0-5% most-central interval gives <M>_class ~ 237 versus <M> = lambda+mu ~ 186, yielding kappa_CBWC,2/kappa_true ~ 1.16 instead of the ~0.94 shown in Fig. 2. The asymptotic lambda->infinity suppression survives, but the quantitative ratios in Figs. 2-3 do not apply to STAR-like centrality bins.
  2. [§VI, Eq. (67), and Conclusion] The statement that overcompensation 'starts already at moderate correlation strengths' is based on the full-range calculation. The estimate in the previous comment shows that for a 0-5% bin the correction has the opposite sign at lambda/nu=1.4, so the threshold for overcorrection depends on the centrality interval. The authors should either compute truncated CBWC cumulants or restrict the claim to the full minimum-bias sample; as written, the central conclusion is not established for the experimental context that motivates the paper.
minor comments (5)
  1. [Eqs. (17)-(18)] The notation delta M[i] is used both for a bin and as a summation range; clarifying that the sum over i partitions the interval would help.
  2. [Eqs. (43) and (45)] The summation index j is reused for the Stirling-index sum and for the bin index; please use different symbols to avoid confusion.
  3. [Fig. 1 caption] The curves labeled 'Full (*0.07)' and 'True (*0.07)' are rescaled by a factor 0.07, but the caption does not explain why the scaling is applied; the reader may mistake it for part of the data.
  4. [Eq. (67)] The arrow '-\rightarrow' is non-standard; use a standard limit arrow such as '\to'.
  5. [Appendix C, Eq. (C5)] The text writes 'kk' where it means k^k; please typeset as k^k throughout (see also Eq. (44)).

Circularity Check

0 steps flagged

No significant circularity: the CBWC overcorrection is a derived property of the stated Poisson+Gamma model, not an input.

full rationale

The derivation is self-contained. Starting from explicit model assumptions — a Gamma-distributed reduced volume Q(x) (Eq. 29), bivariate/tri-variate Poisson sources (Eqs. 30 and 53), and volume-scaled means — the paper computes the fixed-volume ('true') cumulants, the volume-fluctuated cumulants, and the bin-conditioned cumulants via generating functions (Eqs. 32–45 and Appendices C–D). The central inequalities, e.g. κ_CBWC_2 = λ+ν+ν²/(k+λ+μ)−λ²/(λ+μ) < λ+ν for large λ (Eq. 67), and the ratios in Figs. 2–3 follow algebraically from Eqs. (47)–(49) and (63); no parameter is fitted to the CBWC/true ratio or to the conclusion. The self-references [10] and [14] are used for standard cumulant identities and the volume-fluctuation decomposition, but the needed identities are re-derived in the text (Eqs. 39–49 and Appendix D), so these citations are not load-bearing. The finite-centrality-interval issue noted by the skeptic — Eq. (45) sums over the full multiplicity distribution while Sec. II defines a class restricted to [M1,M2] — is a modeling/quantitative-validity caveat, not a circularity, because no term in the derivation is defined in terms of the conclusion. The paper also flags its own limitations (unit-width bins, semi-realistic model, no direct comparison with data) and reports a numerical convergence check for bin width, which supports that the result is an honest model calculation rather than a renamed or fitted input.

Axiom & Free-Parameter Ledger

5 free parameters · 6 axioms · 0 invented entities

The paper is a parameter study: the model parameters are inputs chosen to illustrate different collision energies, and none is fitted to the paper's target (the 'true' cumulants). The only structurally load-bearing choice that is neither derived nor flagged is summing over all multiplicities rather than over a finite centrality interval.

free parameters (5)
  • Gamma shape parameter k = k=83 in Figs. 2-3
    Sets the volume fluctuation strength (relative variance 1/k); chosen by hand, with no derivation or citation for the value.
  • primordial proton Poisson mean nu = 21 (7.7 GeV), 8.5 (39 GeV), 14 (LHC)
    Input mean for the proton Poisson source, taken from thermal model scenarios (Ref. [17]) to choose realistic baryon densities.
  • primordial antiproton Poisson mean nu_bar = 2.3 (39 GeV), 14 (LHC)
    Input mean for the antiproton Poisson source in the net-proton model.
  • pion multiplicity Poisson mean mu = 157, 360, 1150 in figures
    Input mean for the charged-pion multiplicity, chosen from thermal model particle yields.
  • resonance correlation means lambda, lambda_bar = varied in Figs. 2-3; e.g., lambda=18, lambda_bar=2.5
    Controls the correlation between multiplicity and (anti)proton number; the x-axis of the figures scans these values to probe the overcorrection threshold.
axioms (6)
  • domain assumption Volume fluctuations follow a Gamma distribution with unit mean and variance 1/k (Eq. 29).
    Used throughout Section IV and Appendices A, C, D; an analytic convenience that is not derived from data.
  • domain assumption At fixed volume, the joint distribution of proton number and multiplicity (and antiprotons) is a bi-/tri-variate Poisson distribution (Eqs. 30, 53).
    Central model assumption; a Poisson number of resonances decaying into proton-plus-pion creates the shared component lambda that produces correlations.
  • domain assumption All Poisson means scale linearly with reduced volume x (Eq. 31 and Section V).
    Volume-proportional particle production; the authors note this is the leading term for large systems.
  • domain assumption The multiplicity M contains only pions, not protons or antiprotons.
    Matches the STAR analysis where charged-particle multiplicity excludes protons.
  • ad hoc to paper The sums defining CBWC run over the full multiplicity distribution with no [M1,M2] cutoff (Eqs. 33, 44-45).
    Section II defines a finite centrality interval, but the analytical implementation sums over all M and normalizes p_j over the whole distribution; this restriction is never implemented.
  • standard math Stirling numbers of the first and second kind and the Faa di Bruno formula are used to differentiate composite generating functions (Appendices C-D).
    Standard combinatorics used for cumulant expansion; not in dispute.

pith-pipeline@v1.3.0-alltime-deepseek · 18097 in / 17356 out tokens · 157599 ms · 2026-08-03T22:08:26.540478+00:00 · methodology

0 comments
Cite this review

Pith. "Pith review of To bin or not to bin: does binning in multiplicity reliably suppress unwanted volume fluctuations?." pith.science (2026). https://pith.science/paper/VZ2ZRY5O

@misc{pith2026251111869,
  author       = {Pith},
  title        = {Pith review of: To bin or not to bin: does binning in multiplicity reliably suppress unwanted volume fluctuations?},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VZ2ZRY5O}},
  note         = {Machine review of arXiv:2511.11869}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

In this study, we examine the effect of the so-called Centrality Bin Width Correction (CBWC) on the measurement of (net-)proton number cumulants in nucleus-nucleus collisions. We present an analytically tractable model, which includes correlations between multiplicity and proton number similar to those generated by the decay of baryon resonances. Within this model, we analyze the circumstances under which the CBWC method correctly removes the undesired effects of volume or impact parameter fluctuations. Additionally, we explore situations where the method fails and produces misleading results.

Figures

Figures reproduced from arXiv: 2511.11869 by Bengt Friman, Volker Koch.

Figure 1
Figure 1. Figure 1: FIG. 1. The probability distribution for [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Ratio of CBWC cumulants to their true values, [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Ratio of CBWC cumulants to their true values, [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Exploring the QCD phase diagram through correlations and fluctuations

    nucl-th 2025-12 accept novelty 3.0

    A synthesis of the QCD critical point search: theoretical estimates cluster at T_C≈100–120 MeV, μ_B≈550–650 MeV, and BES-II cumulant data match non-critical baselines above 10 GeV while 7.7–9 GeV deviations remain une...

Reference graph

Works this paper leans on

21 extracted references · 16 linked inside Pith · cited by 1 Pith paper

  1. [1]

    Bleicher, S

    M. Bleicher, S. Jeon, and V. Koch, Phys. Rev.C62, 061902 (2000), arXiv:hep-ph/0006201 [hep-ph]

    Pith/arXiv arXiv 2000

  2. [2]

    Bzdak, V

    A. Bzdak, V. Koch, and V. Skokov, Phys. Rev. C87, 014901 (2013), arXiv:1203.4529 [hep-ph]

    Pith/arXiv arXiv 2013

  3. [3]

    Vovchenko, O

    V. Vovchenko, O. Savchuk, R. V. Poberezhnyuk, M. I. Gorenstein, and V. Koch, Phys. Lett. B811, 135868 (2020), arXiv:2003.13905 [hep-ph]

    Pith/arXiv arXiv 2020

  4. [4]

    Vovchenko, R

    V. Vovchenko, R. V. Poberezhnyuk, and V. Koch, JHEP10, 089 (2020), arXiv:2007.03850 [hep-ph]

    Pith/arXiv arXiv 2020

  5. [5]

    Braun-Munzinger, B

    P. Braun-Munzinger, B. Friman, K. Redlich, A. Rustamov, and J. Stachel, Nucl. Phys. A1008, 122141 (2021), arXiv:2007.02463 [nucl-th]

    Pith/arXiv arXiv 2021

  6. [6]

    Braun-Munzinger, K

    P. Braun-Munzinger, K. Redlich, A. Rustamov, and J. Stachel, JHEP08, 113 (2024), arXiv:2312.15534 [nucl-th]

    Pith/arXiv arXiv 2024

  7. [7]

    Kitazawa and M

    M. Kitazawa and M. Asakawa, Phys. Rev.C85, 021901 (2012), arXiv:1107.2755 [nucl-th]

    Pith/arXiv arXiv 2012

  8. [8]

    Kitazawa and M

    M. Kitazawa and M. Asakawa, Phys.Rev.C86, 024904 (2012), arXiv:1205.3292 [nucl-th]

    Pith/arXiv arXiv 2012

  9. [9]

    Event by event fluctuations,

    S. Jeon and V. Koch, “Event by event fluctuations,” inQuark-gluon plasma. Vol. 3, edited by X. W. R. Hwa (World Scientific, 2004) pp. 430–490, arXiv:hep-ph/0304012

    Pith/arXiv arXiv 2004

  10. [10]

    Skokov, B

    V. Skokov, B. Friman, and K. Redlich, Phys. Rev.C88, 034911 (2013), arXiv:1205.4756 [hep-ph]

    Pith/arXiv arXiv 2013

  11. [11]

    Rustamov, J

    A. Rustamov, J. Stroth, and R. Holzmann, Nucl. Phys. A1034, 122641 (2023), arXiv:2211.14849 [nucl-th]

    Pith/arXiv arXiv 2023

  12. [12]

    Holzmann, V

    R. Holzmann, V. Koch, A. Rustamov, and J. Stroth, Nucl. Phys. A1050, 122924 (2024), arXiv:2403.03598 [nucl-th]

    Pith/arXiv arXiv 2024

  13. [13]

    X. Luo, J. Xu, B. Mohanty, and N. Xu, J. Phys. G40, 105104 (2013), arXiv:1302.2332 [nucl-ex]

    Pith/arXiv arXiv 2013

  14. [14]

    Friman and K

    B. Friman and K. Redlich, (2022), arXiv:2205.07332 [nucl-th]

    Pith/arXiv arXiv 2022

  15. [15]

    Riordan, Bull

    J. Riordan, Bull. Amer. Math. Soc.52, 664 (1946)

    1946

  16. [16]

    Comtet,Advanced Combinatorics(D

    L. Comtet,Advanced Combinatorics(D. Riedel, Dordrecht, 1974)

    1974

  17. [17]

    Vovchenko and H

    V. Vovchenko and H. Stoecker, Comput. Phys. Commun.244, 295 (2019), arXiv:1901.05249 [nucl-th]

    Pith/arXiv arXiv 2019

  18. [18]

    J. T. Campbell, Proceedings of the Edinburgh Mathematical Society4, 18–26 (1934)

    1934

  19. [19]

    Kawamura, Kodai Math

    K. Kawamura, Kodai Math. Sem. Rep.25, 246 (1973)

    1973

  20. [20]

    J. G. Skellam, Journal of the Royal Statistical Society109, 296 (1946)

    1946

  21. [21]

    Schumann, (2019), arXiv:1903.03899 [math.CA]

    A. Schumann, (2019), arXiv:1903.03899 [math.CA]

    Pith/arXiv arXiv 2019