{"id":"f6004343-9d0c-4f26-b490-a991df52653c","arxiv_id":"2610.30558","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"The double negative binomial distribution, constrained by kT factorization, describes LHC multiplicity distributions with four free parameters instead of six.","lead":"This paper uses the Color Glass Condensate formalism to describe how many charged particles are produced in high-energy proton-proton collisions. By linking a statistical model to QCD theory, the authors reduce the number of free parameters and still match LHC data accurately.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"Determination of the cutoff Lambda via a global fit to the unconstrained six-parameter DNBD is circular; the claimed parameter reduction from six to four is partially self-referential because Lambda is fitted to the same data, making the reduction less independent than stated.","rationale":"The reader's weakest assumption is the circular determination of Lambda, which I endorse as the most load-bearing concern. The central argument—that kT factorization provides a physical basis for reducing the DNBD's free parameters—requires Lambda to be fixed independently. The paper's own procedure (fitting to the six-parameter results) undermines this physical interpretation, even though the tests with alternative Lambda values show that the fits remain good. However, those tests do not break the circularity because they lack an external anchor. If Lambda can be fixed from, e.g., the saturation scale without loss of fit quality, the claim becomes robust; if not, the reduction is partially circular and the predictive advantage is overstated. The authors acknowledge the limitation, so the verdict remains CONDITIONAL (unchanged). The concrete test I propose would resolve whether the circularity is essential or merely a convenience.","tokens_in":54504,"tokens_out":4214,"duration_ms":44091,"concrete_test":"Fix Lambda from the saturation scale Q_s at a typical small-x value (e.g., using the GBW parametrization Q_s^2 = Q0^2 (x0/x)^lambda with Q0=0.6 GeV, x0=0.01, lambda=0.205, giving Q_s ~ 1.2–1.5 GeV at LHC energies). Then refit all multiplicity distributions with the four remaining free parameters (lambda, alpha, k_s, k_sh) and compare the chi-squared/dof and the energy dependence of alpha to the values in Tables 5.2 and 5.3. If the fits are of comparable quality, the circularity is not essential; if they degrade significantly, the model's success depends on the circular calibration.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that kT factorization reduces the DNBD from six to four free parameters hinges on the cutoff Lambda = 1.4 GeV, which is determined (Section 5.2) by fitting Equation (4.10) to the products lambda*alpha*<n>_s and lambda*(1-alpha)*<n>_sh obtained from a full six-parameter DNBD fit to the same data. The authors state: 'In this work, we assume that the \"true\" values of lambda, alpha, <n>_s and <n>_sh are those obtained through a full fitting of the DNBD to the multiplicity distribution data.' This makes the parameter reduction circular: Lambda is calibrated to reproduce the six-parameter fit, so the four-parameter model inherits the fit quality rather than deriving it independently from kT factorization. The authors acknowledge this ('The global fit should be regarded only as an initial step to determine it') and test other values (1 and 2 GeV), but those alternatives also lack independent physical justification. The tests show similar chi-squared values, but parameters like alpha change significantly (Figures 5.4a,b), indicating the model is not uniquely constrained. Without an independent determination of Lambda—e.g., from the saturation scale—the claim of a predictive parameter reduction is weakened.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"This MSc dissertation proposes that the double negative binomial distribution (DNBD) for charged-particle multiplicity in pp collisions, when combined with the kT-factorization formalism and the GBW unintegrated gluon distribution, can be reduced from six free parameters to four. The reduction is achieved by identifying the soft and semi-hard components of the kT-factorized expression (Eq. 3.14) with the corresponding DNBD terms via the relations in Eq. (4.10). The separation scale Λ is first determined (Λ=1.4 GeV) by fitting those relations to the parameter products extracted from a full six-parameter DNBD fit. The remaining four parameters (λ, α, k_s, k_sh) are then fitted to LHC multiplicity data over √s = 0.9–13 TeV and several pseudorapidity windows. The fits have χ²/dof values often below unity. The paper also examines KNO scaling and the energy dependence of the fitted parameters.","tokens_in":54906,"tokens_out":2125,"duration_ms":23391,"significance":"If the parameter reduction were truly independent, the work would provide a compact QCD-motivated description of a key soft-QCD observable, connecting the CGC framework to the empirically successful DNBD. The fits are of good quality and the decreasing trend of α with energy is a plausible physical result. However, the significance is tempered by the circularity in the determination of Λ and by the fact that the reduction is demonstrated only within the GBW model. The paper is a valuable phenomenological study but does not yet constitute a sharp test of the CGC formalism.","major_comments":[{"comment":"The determination of Λ is circular: Λ is chosen so that the kT-factorization soft and semi-hard components (Eq. 3.14) reproduce the products λα⟨n⟩_s and λ(1−α)⟨n⟩_sh obtained from a full six-parameter DNBD fit to the same multiplicity data (sentence b344: 'we assume that the “true” values of λ,α,⟨n⟩_s and ⟨n⟩_sh are those obtained through a full fitting of the DNBD'). The four-parameter model therefore inherits, rather than independently generates, the fit quality of the six-parameter model. The authors acknowledge the issue (sentence b347: 'The global fit should be regarded only as an initial step') and test Λ=1 and 2 GeV, but those alternatives also lack independent physical justification. Without an independent cross-check, the claimed reduction from six to four parameters is partly self-referential and the predictive power is overstated.","section":"§5.2, Eq. (4.10) and surrounding text","authors_acknowledge":"The global fit should be regarded only as an initial step to determine it.","authors_acknowledge_p":0.66},{"comment":"Most χ²/dof values are below unity, which the authors attribute to the inclusion of both statistical and systematic errors (sentence b363: 'most of the values of χ²/dof are below unity, which can be attributed to the large errors of the data points'). While conservative errors are welcome, χ²/dof significantly below unity can also indicate overfitting or an inflated error budget. The paper should discuss whether the fit parameters are overconstrained, especially given that the four-parameter model is anchored to the six-parameter fit, and should provide a stability test (e.g., leave-one-out cross-validation or splitting the dataset).","section":"§5.2, Tables 5.2 and 5.3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The parameters K/S used in the multiplicity analysis are taken from fits to pseudorapidity distributions (Table 5.1). However, the multiplicity data and the pseudorapidity data do not always share the same event class, pseudorapidity window, or experimental collaboration. The authors note (sentence b368) that for some windows the mismatch leads to larger χ². This inconsistency makes the analysis less controlled; the authors should provide a sensitivity study or restrict the analysis to matching datasets.","section":"§5.1 and §5.2","authors_acknowledge":"In these cases, the normalization parameter K/Sis fitted to pseudo-rapidity distribution data covering−2.4< η <2.4, which may not provide the best description for multiplicity distribution data in wider windows.","authors_acknowledge_p":0.91}],"minor_comments":[{"comment":"The caption says 'The data points in blue refer to INEL>0 events. All the others refer to NSD events.' The blue points (INEL>0) are at η_c = 0.8 for 5.02 and 13 TeV. In Fig. 5.2, these appear with different vertical offsets, but the reader cannot easily separate the two event classes in the main plot. Adding distinct marker shapes or a separate panel would improve clarity.","section":"Figure 5.2 (a) and (b)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract states that the parameter reduction is 'effectively' achieved, but the text later acknowledges the circularity. The abstract should briefly note that Λ is determined from the data, to avoid misleading a casual reader.","section":"Abstract and §5.4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Equation (3.6) is an approximate form of the full kT factorization, justified (sentence b265) 'Keeping only these kinematic regions yields the following approximate expression.' The paper uses this approximation throughout, but the impact of the approximation on the extracted Λ and on the fit quality is not quantified. A brief estimate of the uncertainty introduced by the approximation would strengthen the analysis.","section":"Eq. (3.6) and §3.2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Tables A.1 and A.3 show that for Λ=1 GeV the ⟨n⟩_sh exceeds ⟨n⟩_s, while for Λ=2 GeV the opposite occurs. The authors remark (sentence b385) that the origin is not understood and may be related to the fitting algorithm. This should be explored further; an uncontrolled dependence on the fitting algorithm suggests that the decomposition is not robust.","section":"Table 5.2 and Table A.1, A.3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The KNO scaling plots show the model curves for each energy, but the experimental data points are not shown. Overlaying a subset of the data in one or two windows would demonstrate how well the model captures the energy dependence.","section":"§5.3, Fig. 5.5","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision","confidential_remarks":"The paper is a dissertation chapter rather than a standalone research article; it reads like a competent MSc work. The novelty is incremental — the connection between kT factorization and DNBD parameters is interesting, but the circularity in the Λ determination undermines the central claim of parameter reduction. The authors should either (a) determine Λ from an independent observable (e.g., the saturation scale Q_s(x) at the relevant x) and show that the four-parameter model still fits, or (b) reframe the claim to state that the model provides a good fit when Λ is calibrated on the DNBD products, and then discuss the resulting energy trends as a posteriori QCD-motivated parametrizations rather than a genuine reduction of free parameters. Without such a revision, the paper would not pass a strong phenomenological journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"The paper connects kT factorization with the double negative binomial distribution (DNBD), reducing the free parameters from six to four via a cutoff Λ = 1.4 GeV. The fits to LHC multiplicity distributions are genuinely good, with χ²/dof often below unity, and the energy dependence of the soft fraction α is consistent across pseudorapidity windows. That is a solid piece of phenomenology worth a look.\n\nThe main soft spot is the determination of Λ. As the authors explicitly state, Λ is obtained by fitting the constrained model to the products λα⟨n⟩_s and λ(1-α)⟨n⟩_sh derived from the unconstrained six-parameter DNBD fit. That is circular: the reduction in parameters is calibrated to reproduce the same data, so the claim of independent predictive power is weakened. The authors acknowledge this (Section 5.2) and test Λ = 1 and 2 GeV, finding similar fit quality. Those tests help, but they do not remove the circularity—they just show the fit is not hugely sensitive to Λ. Without an independent determination, e.g., from the saturation scale, the parameter reduction is partly self-referential.\n\nOn the positive side, the paper is well-organized and the analysis is thorough. The KNO scaling discussion is sensible, and the observation that k_total is constant only in the central window is a useful detail. The authors are honest about limitations, including the use of an approximate form of kT factorization and the mismatch between pseudorapidity windows in the K/S fits.\n\nWho should read this? Phenomenologists working on soft/hard separation in QCD and multiplicity distributions. It is a competent master's thesis study that could become a decent journal paper after tightening the Λ determination or at least making the circularity more transparent. I would send it to peer review—it has enough substance and honest discussion to merit referee time. I would not cite it in my own work until the cutoff is better justified, but I would recommend it to a student as a case study in applying kT factorization to multiplicity data.","headline":"Competent parameter reduction via kT factorization, but the cutoff determination is circular—worth engaging but not a breakthrough.","tokens_in":55361,"tokens_out":1566,"would_cite":false,"duration_ms":20618,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek/deepseek-v4-flash","headline":"A four-parameter double negative binomial distribution, constrained by kT factorization, describes charged-particle multiplicity in pp collisions at LHC energies.","keywords":["multiplicity distribution","double negative binomial distribution","kT factorization","Color Glass Condensate","GBW model","soft and semihard components","KNO scaling","proton-proton collisions"],"falsifier":"If new pp collision data at higher energies (e.g., 14 TeV or beyond) or in wider pseudorapidity windows cannot be fitted by the four-parameter DNBD with Λ=1.4 GeV and the same GBW-based components, then the claimed reduction is not universal. Conversely, an independent experimental determination of the soft/semihard separation scale (e.g., from the transverse momentum spectrum or from event shapes) that disagrees significantly with Λ=1.4 GeV would challenge the interpretation.","tokens_in":55980,"feed_emoji":"⚛️","tokens_out":3876,"duration_ms":39733,"temperature":0.7,"pith_summary":"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.","feed_headline":"kT factorization cuts DNBD parameters from six to four for LHC pp multiplicity","feed_subtitle":"A Λ=1.4 GeV cutoff separates soft and semihard components, yielding χ²/dof often below unity across energies 0.9–13 TeV and seven…","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduces the kT factorization formalism for gluon production in hadronic collisions, which is the foundation of the paper's two-component decomposition.","marker":"[24]"},{"why":"Provides the Golec-Biernat–Wüsthoff (GBW) model for the unintegrated gluon distribution, which the paper uses to compute the soft and semihard components.","marker":"[44]"},{"why":"Develops the double negative binomial distribution (DNBD) as a superposition of soft and semihard components, which the paper constrains with kT factorization.","marker":"[67]"},{"why":"Provides the ALICE collaboration data on charged-particle multiplicity distributions at √s=0.9–8 TeV, against which the paper fits its four-parameter DNBD.","marker":"[69]"},{"why":"Provides the CMS collaboration data on charged-particle multiplicity distributions at √s=0.9, 2.36, and 7 TeV, used in the paper's fits.","marker":"[75]"},{"why":"Supplies the running coupling prescription and the GBW gluon density formula with (1-x)^4 factors, which the paper adopts in the kT factorization.","marker":"[53]"},{"why":"Establishes KNO scaling, which the paper tests using its fitted DNBD parameters to determine energy independence.","marker":"[58]"}],"fun_headline_variants":["kT factorization cuts DNBD parameters from six to four at LHC","Shoulder-like pp multiplicity fit by kT factorized DNBD at Λ=1.4 GeV","Four-parameter DNBD describes LHC pp data across 0.9–13 TeV","KNO scaling holds only in narrow |η|<0.5 for LHC multiplicity"],"cache_read_input_tokens":0,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["kT factorization cuts DNBD parameters from six to four at LHC","Shoulder-like pp multiplicity fit by kT factorized DNBD at Λ=1.4 GeV","Four-parameter DNBD describes LHC pp data across 0.9–13 TeV","KNO scaling holds only in narrow |η|<0.5 for LHC multiplicity"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":9e-05,"raw_usage":{"cost":8.9964e-05,"is_byok":false,"cost_details":{"upstream_inference_cost":8.9964e-05,"upstream_inference_prompt_cost":5.367852e-05,"upstream_inference_completions_cost":3.628548e-05},"total_tokens":1437,"prompt_tokens":1074,"completion_tokens":363,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":0,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":270}},"tokens_in":1074,"tokens_out":363,"cost_source":"openrouter_usage","duration_ms":4725,"image_pages":[],"temperature":1.0,"billed_attempts":1,"reasoning_tokens":270,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-23T14:46:37.403626+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"falsifier":"If new pp collision data at higher energies (e.g., 14 TeV or beyond) or in wider pseudorapidity windows cannot be fitted by the four-parameter DNBD with Λ=1.4 GeV and the same GBW-based components, then the claimed reduction is not universal. Conversely, an independent experimental determination of the soft/semihard separation scale (e.g., from the transverse momentum spectrum or from event shapes) that disagrees significantly with Λ=1.4 GeV would challenge the interpretation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the kT factorization formalism for gluon production in hadronic collisions, which is the foundation of the paper's two-component decomposition."},{"cited_title":"Golec-Biernat and M","cited_arxiv_id":null,"evidence_quote":"Provides the Golec-Biernat–Wüsthoff (GBW) model for the unintegrated gluon distribution, which the paper uses to compute the soft and semihard components."},{"cited_title":"Possible scenarios for soft and semi-hard components structure in central hadron-hadron collisions in the TeV region","cited_arxiv_id":"hep-ph/9810446","evidence_quote":"Develops the double negative binomial distribution (DNBD) as a superposition of soft and semihard components, which the paper constrains with kT factorization."},{"cited_title":"Charged-particle multiplicities in proton-proton collisions at $\\sqrt{s}$ = 0.9 to 8 TeV","cited_arxiv_id":"1509.07541","evidence_quote":"Provides the ALICE collaboration data on charged-particle multiplicity distributions at √s=0.9–8 TeV, against which the paper fits its four-parameter DNBD."},{"cited_title":"Charged particle multiplicities in pp interactions at sqrt(s) = 0.9, 2.36, and 7 TeV","cited_arxiv_id":"1011.5531","evidence_quote":"Provides the CMS collaboration data on charged-particle multiplicity distributions at √s=0.9, 2.36, and 7 TeV, used in the paper's fits."},{"cited_title":"Gluon saturation in $pA$ collisions at the LHC: KLN model predictions for hadron multiplicities","cited_arxiv_id":"1111.3031","evidence_quote":"Supplies the running coupling prescription and the GBW gluon density formula with (1-x)^4 factors, which the paper adopts in the kT factorization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes KNO scaling, which the paper tests using its fitted DNBD parameters to determine energy independence."}],"review_version":1}