Pith. sign in

REVIEW 4 major objections 4 minor 1 cited by

Charged particle multiplicity distributions derived from the Principle of Maximal Entropy

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read Maximum entropy alone yields measured particle multiplicity shapes

desk verdict The paper's central derivation doesn't hold: on the discrete support of multiplicity, max entropy gives a geometric, not exponential, distribution, and the fit cuts the low-n region where they differ. read the letter →

arxiv 2505.23491 v1 pith:ATOE5LJX submitted 2025-05-29 hep-ph

Sándor Lökös This is my paper
classification hep-ph
keywords chargedparticlemultiplicitymaximumentropyShannonnegativebinomialdistributionGammaPoisson-GammamixturerapiditywindowsLHCbdata
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

The paper sets out to show that the measured charged-particle multiplicity distributions in high-energy collisions can be derived from the principle of maximum entropy alone, with no dynamical input beyond normalization and a fixed mean. In narrow rapidity windows, where the mean is approximately constant, the maximum-entropy distribution is the exponential $p(n)=\lambda e^{-\lambda n}$, and the paper finds that this fits LHCb data with an acceptable $\chi^2$ while the extracted Shannon entropy stays roughly constant across rapidity. For wider rapidity windows, the paper argues that convoluting $k$ such exponentials yields a Gamma distribution, and the Poisson-Gamma mixture turns into the negative binomial distribution, offering a derivation of the well-known NBD shape that is more general than the clan model. A sympathetic reader would care because the result would mean that two ubiquitous empirical features of multiplicity data are consequences of information-theoretic indifference, not of specific production mechanisms.

What carries the argument

The machinery is the variational entropy functional $F[p,\alpha,\beta] = -\int_0^\infty p(n)\ln p(n)\,dn + \alpha\left(\int_0^\infty p(n)\,dn - 1\right) + \beta\left(\int_0^\infty n p(n)\,dn - \mu\right)$, whose maximization produces $p(n)=\lambda e^{-\lambda n}$ with $\lambda=1/\mu$; the normalization-only version of the same functional gives the uniform $p(n)=1/N$ used for the initial partonic state. The second ingredient is the Poisson-Gamma mixture identity: if $Y\sim\Gamma(k,\lambda)$ and $X\sim\mathrm{Poisson}(\mu)$ with $\mu=Y$, then $P(X=x\mid\mu)=\mathrm{NBD}(k,\bar n=k/\lambda)$, which converts the $k$-fold convolution of exponentials into the observed negative binomial shape.

What would settle it

Measure a charged-particle multiplicity distribution in a narrow rapidity window with high statistics and test whether $\ln P(n)$ is strictly linear in $n$ when the low-multiplicity points ($n<3$) that the current fit excludes are included; a significant curvature or a low-$n$ excess would falsify the exponential prediction, and a fitted mean that drifts across sub-windows would falsify the constant-mean assumption on which it rests.

Watch

Extended reading notes

Core claim

The central discovery claimed is that maximum Shannon entropy, subject only to normalization and a fixed mean $\mu$, reproduces the charged-particle multiplicity data: in a narrow rapidity window the unique maximizer is the exponential distribution $p(n)=\lambda e^{-\lambda n}$ with $\lambda=1/\mu$, and this form describes the measured LHCb multiplicities in 0.5-unit rapidity windows, with the fitted entropy consistent with a constant. The same principle is then extended to wider windows: if the exponential is the building block of each narrow slice, then the wider-window distribution is a $k$-fold convolution of exponentials, i.e., a Gamma distribution $\Gamma(k,\lambda)$; mixing that Gamma with a Poisson yields the negative binomial distribution, with parameter relations $k_\Gamma=k_{\mathrm{NBD}}$ and $\lambda_\Gamma=k_{\mathrm{NBD}}/\langle n\rangle$. Thus the paper claims that the famous NBD shape of charged multiplicities is not an input but an output of the maximum-entropy requirement.

Load-bearing premise

The load-bearing premise is that the multiplicity count $n$ can be treated as a continuous variable on $[0,\infty)$ and that normalization plus a fixed mean, with the mean value imported from data, are the only constraints; if a continuous distribution with only those two constraints is not the right description of the measurement, the exponential form and the NBD explanation lose their stated foundation.

Editorial extensions

If this is right

  • The exponential tail observed in narrow rapidity windows becomes a generic information-theoretic prediction, requiring no model of particle production.
  • The same maximum-entropy argument yields the negative binomial distribution in wider windows, making the NBD shape a consequence of repeated exponential convolutions rather than of clan dynamics.
  • The NBD parameter $k$ should grow linearly with the width of the rapidity window and the entropy should saturate, both of which the paper notes are consistent with existing data.
  • The extracted Shannon entropy in narrow windows should be constant across rapidity, offering a direct final-state counterpart to the initial-state entanglement entropy conjectured in earlier work.

Reading between the lines

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

  • The derivation implies that any collision system with a fixed charged multiplicity mean in a small phase-space window should show the same exponential shape, so the exponential tail is not a distinctive signature of the underlying production mechanism but a baseline expected from entropy maximization.
  • If the constant-mean assumption were relaxed to allow a mean that varies within the window, the resulting maximum-entropy distribution would be a generalized exponential whose shape could be used to measure the rapidity dependence of the mean directly from the data.
  • The same machinery could be applied to joint multiplicity distributions in multiple rapidity windows by imposing cross-window correlations as additional constraints, potentially exposing short- and long-range correlation parameters without a dynamical model.
  • One testable extension is to apply the Gamma-Poisson mixture with the parameter relations of Eq. (7) to existing wide-window data sets and check whether the predicted $k$--$\langle n\rangle$ correlation matches the measured NBD parameter correlation.
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

4 major / 4 minor

Summary. The paper proposes that charged-particle multiplicity distributions measured in high-energy collisions follow from the Principle of Maximum Entropy (POME) with no a priori physical assumptions. It treats multiplicity n as a continuous variable on [0,∞), first maximizing Shannon entropy under normalization alone and claiming to recover a uniform distribution p(n)=1/N for the initial partonic state, then adding a mean constraint and obtaining the exponential distribution p(n)=λe^{-λn}. The exponential is fitted to LHCb multiplicity data in narrow rapidity windows, and the extracted Shannon entropy S=1-ln(λ) is reported to be approximately constant. The paper then extends the argument by convolving exponential distributions k times to obtain a Gamma distribution and, via a Poisson-Gamma mixture, a negative binomial distribution (NBD), which it presents as a natural explanation of the observed NBD shape.

Significance. If the central claim were sound, the paper would offer a strikingly parsimonious derivation of the exponential and NBD forms of multiplicity distributions, connecting them to information-theoretic maximal entropy without dynamical input. The manuscript does contain a correct statement of the standard max-entropy result for a continuous distribution on [0,∞) with fixed mean, and the comparison to LHCb data, though limited, is a useful sanity check. However, the central claim as stated is not supported: the normalization-only maximization on an infinite domain is ill-posed, the discrete nature of charged multiplicities is ignored without any error estimate, and the apparent 'constant entropy' is a monotone function of the fitted mean, so it is not an independent prediction. These issues directly undermine the advertised assumption-free derivation, so the paper does not, in its present form, establish its main conclusion.

major comments (4)
  1. [Eqs. (2)-(3), Sec. 2.1] The normalization-only maximization on [0,∞) is not well defined. A constant density p(n)=e^{α-1} on an infinite interval cannot be normalized to unity: the integral ∫_0^∞ dn diverges, so no finite Lagrange multiplier enforces the constraint. Consequently, the claimed derivation of p(n)=1/N has no valid variational formulation. This is load-bearing because it is used to recover the uniform partonic distribution in Eq. (1) and to identify N with the number of states.
  2. [Eqs. (4)-(5), Sec. 2.2 and Fig. 1] The exponential form p(n)=λe^{-λn} is obtained by maximizing continuous entropy on [0,∞) with a fixed mean, but charged multiplicities are integers. On the actual support n∈{0,1,2,...}, the same max-entropy problem with normalization and mean μ yields the geometric distribution p(n)=(1/(1+μ))(μ/(1+μ))^n, not the exponential. The paper provides no estimate of the error introduced by the continuum approximation, and the fit in Fig. 1 excludes n<3, precisely the region where discrete and continuous forms differ most. This unsupported discretization is an a priori assumption that contradicts the claim of an assumption-free derivation.
  3. [Sec. 4, Eq. (7)] The identification of the number of convolutions k with the number of partonic states N in Eq. (1) is asserted without derivation. Convolving exponential distributions k times gives a Gamma distribution, but k is a free parameter (the NBD shape parameter) that is not fixed by POME. The statement that 'the origin of the NBD shape lies in the convolution of maximal entropy distributions repeated k=N times' therefore rests on an additional assumption, not on the maximization principle itself.
  4. [Sec. 3 and Eq. (5)] The reported constant Shannon entropy is not an independent prediction. From Eq. (5), S=1-ln(λ)=1+ln(μ), so the extracted entropy is a monotone function of the fitted mean μ. If μ is approximately constant across rapidity windows, S is constant by construction; the agreement with H1 therefore provides no independent confirmation of the max-entropy picture beyond the constancy of the fitted mean.
minor comments (4)
  1. [Sec. 1] There is a typo: 'he essential idea' should be 'The essential idea'. Also, the phrase 'by its nature' is used twice in one sentence in Sec. 2.
  2. [Fig. 1] The axis labels and data-range notation in the left panel are difficult to read due to formatting of the probability axis; the plot would benefit from standard logarithmic axis labels and an explicit statement of the fitting range.
  3. [Sec. 4] The sentence 'The Gamma distribution is known to be the scaling function of the negative binomial distribution' is unclear; it would be more precise to state that the Gamma distribution is the continuous mixing distribution in a Poisson-Gamma mixture.
  4. [Sec. 4, Eq. (7)] The relation k_Γ = k_NBD and λ_Γ = k_NBD/⟨n⟩ is stated without derivation; a brief derivation would help the reader follow the parameter mapping.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the central POME derivation is a standard maximum-entropy calculation whose parameter is fitted from data, not predicted from the result.

full rationale

The paper's core step (Eqs. 4-5) maximizes Shannon differential entropy on [0,∞) subject to normalization and a fixed mean μ. The exponential solution p(n)=λe^{-λn} with λ=1/μ is the standard variational result; μ is an input constraint imported from data, and the functional form is the output. This is not circular: the mean is not claimed to be predicted by POME. The NBD extension (Eq. 6) is a standard Gamma-Poisson mixture, and the number of convolutions k is acknowledged to be a free parameter not fixed by POME; the paper even notes that the integer-k prediction is not confirmed by data. The extracted Shannon entropy S=1-lnλ=1+lnμ is a deterministic function of the fitted mean, so the observed constant entropy is a restatement of the fitted μ values rather than an independent prediction; however, the paper presents it as an extraction, not as a prediction, and it is not load-bearing for the main derivation. No load-bearing self-citations or imported uniqueness theorems appear. Concerns about the continuous support, the non-normalizability of Eq. (3), and the exclusion of n<3 in the fit are correctness or modeling issues, not circularity.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The central derivation rests on constraints whose values come from the data, on a continuous approximation to a discrete observable, and on an unproven identification between the number of convolutions k and the partonic state count N. The ledger is therefore not empty: mu and k carry the predictive weight, and the normalization-only uniform result in Eq. (3) is not a well-defined maximum on an infinite domain.

free parameters (3)
  • mu (mean charged multiplicity per rapidity bin) = example mu = 2.10 +/- 0.08 for 2 < eta < 2.5
    Eq. (5) sets lambda=1/mu; mu is imported from the measured distribution for each bin, so the exponential scale is fitted, not predicted.
  • k (number of convolutions / NBD shape parameter) = not fitted in this paper; in standard NBD fits k is a free parameter
    Sec. 4 identifies k with the partonic state count N and with the NBD parameter, but does not derive its value; the paper admits the integer-k consequence is not confirmed by data.
  • N (number of equally probable partonic states) = taken from xg(x,Q)+xSigma(x,Q) in Eq. (1)
    N is inherited from previous entanglement papers, not derived by POME; it is reused as the number of convolutions in Sec. 4.
assumptions (6)
  • domain assumption Only constraints needed are normalization and a fixed mean; higher moments carry no information
    Used to get the exponential in Eq. (5); not justified from data or physics, and the fit itself excludes low-n points.
  • domain assumption Multiplicity n treated as continuous on [0, infinity)
    Eqs. (2)-(5) use integrals; discrete measured multiplicities are later reintroduced via Poisson mixture in Eq. (6).
  • ad hoc to paper Normalization-only maximization yields a well-defined uniform distribution p=1/N
    Eq. (3) states p=1/N, but the integral in Eq. (2) runs to infinity, where a constant density is not normalizable.
  • ad hoc to paper The number of convolutions k equals the number of partonic states N from Eq. (1)
    Sec. 4 identifies the Gamma shape parameter with N without derivation; paper admits the integer-k prediction is not confirmed.
  • domain assumption Mean is constant within each 0.5-unit-wide rapidity window
    Sec. 2.2 and Fig. 1 rely on this; the approximation breaks for wider windows, as discussed in Sec. 4.
  • standard math Poisson-Gamma mixture maps continuous Gamma to discrete NBD
    Eq. (6) uses a standard result; this step is mathematically sound.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Charged particle multiplicity distributions derived from the Principle of Maximal Entropy." pith.science (2026). https://pith.science/paper/ATOE5LJX

@misc{pith2026250523491,
  author       = {Pith},
  title        = {Pith review of: Charged particle multiplicity distributions derived from the Principle of Maximal Entropy},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ATOE5LJX}},
  note         = {Machine review of arXiv:2505.23491}
}
read the original abstract

Recent theoretical results renewed the interest in charged particle multiplicity distributions. The Shannon entropy of such distributions is conjectured to be related to the entanglement or von Neumann entropy of partonic quantum system. In this paper, we show that the measured charged particle multiplicities can be derived from the principle of maximum entropy (POME or MAXENT) without any a priori physical assumption. The approach provides a natural explanation for the well-known negative binomial shape of the measured distributions.

Figures

Figures reproduced from arXiv: 2505.23491 by the authors.

Figure 1
Figure 1. An example fit of the multiplicity distributions (left) measured by LHCb [17] in [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗

Discussion (0). Sign in to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Deep inelastic scattering as a probe of entanglement: the complete QCD dipole cascade

    hep-ph 2026-07 conditional novelty 6.0 of 10

    The Shannon entropy of dipole multiplicities from the full Levin–Lublinsky equation in DIS reproduces the H1 hadron entropy, growing linearly with ln(1/x) and described by S = ln(2/3⟨n⟩) + 0.85.

Reference graph

Works this paper leans on

21 extracted references · 5 canonical work pages · cited by 1 Pith paper

  1. [1]

    D. E. Kharzeevet al., Phys. Rev. D95, 114008 (2017) [arXiv:1702.03489]

    arXiv 2017

  2. [2]

    O. K. Bakeret al., Phys. Rev. D98, 054007 (2018) [arXiv:1712.04558]

    arXiv 2018

  3. [3]

    Tuet al., Phys

    Z. Tuet al., Phys. Rev. Lett.124, 062001 (2020) [arXiv:1904.11974]

    arXiv 2020

  4. [4]

    Hentschinski et al., Eur

    M. Hentschinski et al., Eur. Phys. J. C82, 111 (2022)[Erratum: Eur.Phys.J.C 83, 1147 (2023)] [arXiv:2110.06156]

    arXiv 2022

  5. [5]

    D. E. Kharzeevet al., Phys. Rev. D104, L031503 (2021) [arXiv:2102.09773]

    arXiv 2021

  6. [6]

    Zhang et al., Phys

    K. Zhang et al., Phys. Rev. D105, 014002 (2022) [arXiv:2110.04881]

    arXiv 2022

  7. [7]

    Hentschinski et al., Eur

    M. Hentschinski et al., Eur. Phys. J. C82, 1147 (2022) [arXiv:2207.09430]

    arXiv 2022

  8. [8]

    Hentschinski et al., Phys

    M. Hentschinski et al., Phys. Rev. Lett.131, 241901 (2023) [arXiv:2305.03069]

    arXiv 2023

Show all 21 references
  1. [9]

    Khachatryan et al., JHEP 01, 079 (2011) [arXiv:1011.5531]

    V. Khachatryan et al., JHEP 01, 079 (2011) [arXiv:1011.5531]

    2011 arXiv

  2. [10]

    Aad et al., Eur

    G. Aad et al., Eur. Phys. J. C76, 403 (2016) [arXiv:1603.02439]

    2016 arXiv

  3. [11]

    Acharyaet al., Eur

    S. Acharyaet al., Eur. Phys. J. C77, 852 (2017) [arXiv:1708.01435]

    2017 arXiv

  4. [12]

    Andreev et al., Eur

    V. Andreev et al., Eur. Phys. J. C81, 212 (2021) [arXiv:2011.01812]

    2021 arXiv

  5. [13]

    E. T. Jaynes, Phys. Rev.106, 620 (1957)

    1957

  6. [14]

    E. T. Jaynes, Phys. Rev.108, 171 (1957)

    1957

  7. [15]

    E. D. Malazaet al., Phys. Lett. B266, 169 (1991). 5

    1991

  8. [16]

    F. S. Navarra, O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, Physica A344, 568 (2004) [arXiv:hep-ph/0312136]

    2004 arXiv

  9. [17]

    Aaij et al., Eur

    R. Aaij et al., Eur. Phys. J. C74, 2888 (2014) [arXiv:1402.4430]

    2014 arXiv

  10. [18]

    R. E. Ansorgeet al., Z. Phys. C43, 357 (1989)

    1989

  11. [19]

    M. R. Adamset al., Z. Phys. C61, 179 (1994)

    1994

  12. [20]

    Giovanniniet al., Nucl

    A. Giovanniniet al., Nucl. Phys. B Proc. Suppl.25, 115 (1992)

    1992

  13. [21]

    Simak et al., Physics Letters B206, 159 (1988)

    V. Simak et al., Physics Letters B206, 159 (1988). 6

    1988

Pith tools

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