Pith. Sign in

REVIEW

Multiplicity Distribution in the Color Glass Condensate Formalism

High Energy Physics - Phenomenology · Submitted 22 Oct 2026 · Published version · 78 references

Pith's verdict: Competent parameter reduction via kT factorization, but the cutoff determination is circular—worth engaging but not a breakthrough.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 23, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.30558 (version 1)
Pith IDA6IIKPHY, Pith's permanent record of this paper and its review
CitationH. R. Martins-Fontes. "Multiplicity Distribution in the Color Glass Condensate Formalism." arXiv:2610.30558 (2026). https://pith.science/paper/A6IIKPHY
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

A four-parameter double negative binomial distribution, constrained by kT factorization, describes charged-particle multiplicity in pp collisions at LHC energies.

The reading

The paper tries to establish that the charged-particle multiplicity distribution in high-energy proton-proton collisions can be described by a double negative binomial distribution (DNBD) with only four free parameters, instead of the usual six, by linking the parameters to the kT factorization formalism of QCD. It introduces a cutoff Λ=1.4 GeV that cleanly separates the kT factorization integral into a soft and a semihard component, and these components are directly identified with the products of DNBD parameters (λα⟨n⟩_s and λ(1-α)⟨n⟩_sh). This connection is the paper's central claim: if it holds, then the multiplicity distribution across a wide range of LHC energies (0.9–13 TeV) and pseudorapidity windows can be fitted with just four parameters, yielding χ² per degree of freedom often below unity. Moreover, the same framework gives a consistent description of pseudorapidity distributions and allows an investigation of KNO scaling, which the paper finds holds only for the most central pseudorapidity window.

What carries the argument

The key machinery is the decomposition of the kT factorization formula (equation 3.1) into two integrals over transverse momentum p_⊥, one from p_⊥,min to Λ (soft component) and one from Λ to p_⊥,max (semihard component), as shown in equation (3.14). These integrals, after proper normalization, are identified with the average multiplicities λα⟨n⟩_s and λ(1-α)⟨n⟩_sh appearing in the double negative binomial distribution (equation 4.7). This identification provides two constraints that reduce the six free parameters (α, λ, ⟨n⟩_s, ⟨n⟩_sh, k_s, k_sh) to four (α, λ, k_s, k_sh). The cutoff Λ is determined by fitting the products λα⟨n⟩_s and λ(1-α)⟨n⟩_sh obtained from a prior six-parameter DNBD fit to the data, yielding Λ=1.4 GeV. The unintegrated gluon distribution used in the kT factorization is taken from the GBW model, with parameters from the literature.

Extended reading notes

Core claim

The core discovery is that the double negative binomial distribution, traditionally requiring six free parameters to fit the shoulder-like multiplicity distributions at LHC, can be reduced to four parameters by imposing constraints derived from the kT factorization formalism. Specifically, two relations (equation 4.10) tie the products λα⟨n⟩_s and λ(1-α)⟨n⟩_sh to the integrals of the kT factorization expression for the soft and semihard components, which are separated by a fixed cutoff Λ. Using the Golec-Biernat–Wüsthoff (GBW) model for the unintegrated gluon distribution, the paper finds that Λ=1.4 GeV produces excellent fits to CMS and ALICE data across energies 0.9–13 TeV and pseudorapidity windows from |η|<0.5 to |η|<3.4. The parameter α (fraction of soft events) decreases with energy in a consistent way across all windows, and the average multiplicities of soft and semihard components are related by ⟨n⟩_sh ≈ 3⟨n⟩_s. KNO scaling is shown to hold only for the narrowest window |η|<0.5, and the violation is linked to the energy dependence of the combined shape parameter k_total.

Load-bearing premise

The cutoff Λ=1.4 GeV is determined by fitting the products λα⟨n⟩_s and λ(1-α)⟨n⟩_sh from a prior six-parameter DNBD fit to the same data, which assumes that the six-parameter fit gives the true values, making the parameter reduction partially circular.

Editorial extensions

If this is right

  • The four-parameter DNBD can be applied to predict multiplicity distributions at future collider energies or in wider pseudorapidity windows, based on the observed systematic energy dependence of α and the k parameters.
  • The paper's methodology provides a quantitative way to separate soft and semihard contributions in particle production, which may be used to study the transition between perturbative and non-perturbative QCD regimes.
  • The finding that KNO scaling holds only for the central pseudorapidity window suggests that the violation at larger windows is driven by the energy dependence of k_total, not by α, offering a new perspective on scaling violations.
  • The weak dependence of fit quality on the precise value of Λ (tested at 1.0, 1.4, 2.0 GeV) indicates that the reduction to four parameters is robust, and Λ can be treated as a fixed scale.
  • The constant ratio ⟨n⟩_sh ≈ 3⟨n⟩_s across energies and windows provides a simple phenomenological rule that future models of multiparticle production should reproduce.

Reading between the lines

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

  • The reduction to four parameters, while impressive, relies on the prior six-parameter fit to determine Λ, which introduces a mild circularity. An independent determination of Λ (e.g., from the saturation scale or from the point where the running coupling becomes perturbative) would strengthen the claim.
  • The decreasing trend of α with energy across all pseudorapidity windows suggests that semihard processes become increasingly dominant at higher energies, consistent with expectations from parton saturation. This trend could be extrapolated to predict α at future higher-energy colliders.
  • The methodology could be extended to proton-nucleus or nucleus-nucleus collisions using appropriate unintegrated gluon distributions (e.g., from the rcBK equation), potentially testing whether the same Λ and parameter reduction holds in different collision systems.
  • The observation that KNO scaling is satisfied only in the central window, and that k_total is the controlling parameter, may indicate that the semi-hard component has a more complex energy dependence at large pseudorapidity, possibly related to the onset of non-linear effects.
Original abstract

The multiplicity distribution (MD) of charged particles in high-energy proton-proton collisions remains an intricate aspect of quantum chromodynamics (QCD). Traditional models, such as the Negative Binomial Distribution (NBD), were fitted and successfully described MDs at lower energies, but fail at higher energies, such as those probed by the Large Hadron Collider (LHC), where a ``shoulder-like'' structure emerges. To address this, we employ the Double Negative Binomial Distribution (DNBD), a weighted combination of two NBDs, representing contributions from soft and semi-hard processes. This model typically involves six free parameters, but by incorporating the $k_T$ factorization formalism and introducing an energy scale $\Lambda$, we effectively reduce the number of free parameters from six to four while maintaining a strong fit to high-energy MD data. The unintegrated gluon distribution used in this work is derived from the Golec-Biernat-W\"usthoff (GBW) model. Additionally, we investigate the KNO scaling and use the $k_T$ factorization formalism to describe the pseudo-rapidity distribution. Our approach demonstrates excellent agreement with recent experimental data, yielding low $\chi^2$ values and highlighting the potential for improvements using more realistic gluon distribution model. This analysis provides valuable insights into the interplay of soft and semi-hard processes in high-energy collisions.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-23: 1 request, 6,393 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.88 · human participants 0.01 · controlled comparison 0.14 · randomized assignment 0.01 · states quantitative result 0.83 · compares with prior results 0.86 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The multiplicity distribution (MD) of charged particles in high-energy proton-proton collisions remains an intricate aspect of quantum chromodynamics (QCD). Traditional models, such as the Negative Binomial Distribution (NBD), were fitted and successfully described MDs at lower energies, but fail at higher energies, such as those probed by the Large Hadron Collider (LHC), where a ``shoulder-like'' structure emerges. To address this, we employ the Double Negative Binomial Distribution (DNBD), a weighted combination of two NBDs, representing contributions from soft and semi-hard processes. This model typically involves six free parameters, but by incorporating the $k_T$ factorization formalism and introducing an energy scale $\Lambda$, we effectively reduce the number of free parameters from six to four while maintaining a strong fit to high-energy MD data. The unintegrated gluon distribution used in this work is derived from the Golec-Biernat-W\"usthoff (GBW) model. Additionally, we investigate the KNO scaling and use the $k_T$ factorization formalism to describe the pseudo-rapidity distribution. Our approach demonstrates excellent agreement with recent experimental data, yielding low $\chi^2$ values and highlighting the potential for improvements using more realistic gluon distribution model. This analysis provides valuable insights into the interplay of soft and semi-hard processes in high-energy collisions.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The multiplicity distribution (MD) of charged particles in high-energy proton-proton collisions remains an intricate aspect of quantum chromodynamics (QCD). background 1.00 0.04 0.02 asserted 1.00
s1 Traditional models, such as the Negative Binomial Distribution (NBD), were fitted and successfully described MDs at lower energies, but fail at higher energies, such as those probed by the Large Hadron Collider (LHC), where a ``shoulder-like'' structure emerges. background 0.60 0.07 0.04 asserted 1.00
s2 To address this, we employ the Double Negative Binomial Distribution (DNBD), a weighted combination of two NBDs, representing contributions from soft and semi-hard processes. method 0.84 0.59 0.06 asserted 1.00
s3 This model typically involves six free parameters, but by incorporating the $k_T$ factorization formalism and introducing an energy scale $\Lambda$, we effectively reduce the number of free parameters from six to four while maintaining a strong fit to high-energy MD data. method 0.67 0.92 0.22 asserted 1.00
s4 The unintegrated gluon distribution used in this work is derived from the Golec-Biernat-W\"usthoff (GBW) model. method 1.00 0.63 0.03 asserted 1.00
s5 Additionally, we investigate the KNO scaling and use the $k_T$ factorization formalism to describe the pseudo-rapidity distribution. method 0.65 0.91 0.05 asserted 0.97
s6 Our approach demonstrates excellent agreement with recent experimental data, yielding low $\chi^2$ values and highlighting the potential for improvements using more realistic gluon distribution model. result 0.99 0.95 0.67 asserted 0.99
s7 This analysis provides valuable insights into the interplay of soft and semi-hard processes in high-energy collisions. interpretation 1.00 0.56 0.16 asserted 0.89

What kind of paper the classifier reports this is: theoretical result 0.65 · evidence from analysis of existing data 0.98.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. 7: The DGLAP equations in the ln 1/x − Q2 plane [18] 1.5 Conclusions In this chapter, we introduced the Deep Inelastic Scattering formalism and the corre￾sponding kinematics. The proton structure function F2 was seen to depend on the quark and anti-quark distributions. It was shown that the parton model predicts Bjorken scaling, but the inclusion of QCD corrections leads to its violation. These corrections lead to the DGLAP evolution equations, which describe the evolution of the distribution functions with Q2 . The DGLAP predicts a strong growth of the gluon distribution xg(x, Q2 )…
Figure 2. 4: Evolution of a dipole with rapidity 2.6 Phenomenological models An analytical solution to the Balitsky-Kovchegov (BK) equation is not known for the entire phase space it describes. The few solutions that exist in the literature are valid only in the asymptotic regions of the evolution. Consequently, several phenomenological mod￾els have been developed to model the scattering amplitude, ensuring that they incorporate the expected characteristics of N, such as: • Color transparency: In the limit r → 0 (r ≪ 1/Qs), the scattering amplitude
Figure 3. 1: Scheme of kT factorization in pp collisions, as implied by equation (3.1). From (3.1), we obtain the multiplicity per unit of rapidity by performing the integral over p⊥ dN dy = 1 S Z p⊥max p⊥min d 2 p⊥E dσ d 3p , (3.4) where S denotes a typical interaction area. As proposed in [54], the dominant contribution
Figure 4. 1: Negative binomial distribution for several values of
Figure 5. 4: Behavior of α with increasing energy for various pseudo-rapidity windows (−ηc < η < +ηc). (a) correspond to the case Λ = 1 GeV and (b) to the case Λ = 2 GeV. The lines are just to guide the eyes. The data points in blue refer to INEL>0 events. All the others refer to NSD events. 5.3 KNO scaling In order to verify whether KNO scaling is observed in the LHC data, we plotted the multiplicity distributions computed with (4.7) and presented in terms of the KNO variable z = n/⟨n⟩. The resulting curves are shown in

Figures reproduced from arXiv:2610.30558 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.