{"id":"fb9e2aad-3a79-457f-ab08-8c88eafa9cde","arxiv_id":"2412.16348","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Measuring the entropy of final-state hadron multiplicities in H1 data gives a Pomeron intercept of 0.322 ± 0.007, consistent with the value from inclusive DIS cross-section scaling.","lead":"The authors used particle-multiplicity measurements from the H1 experiment at HERA to compute an entropy and extract the Pomeron intercept, a parameter controlling how gluon densities grow at small momentum fraction. The value agrees with the one obtained from inclusive deep-inelastic cross-section scaling, suggesting a new, simpler way to find geometric scaling in collider data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equating the multiplicity entropy S_mult with the partonic scaling entropy S_parton is assumed, not derived; a hadronization/fragmentation test is needed before λ_entropy can be called direct.","rationale":"The reader’s weakest_assumption identifies LHPD and ⟨z⟩ independence as the fragile step; I agree that this is the key conceptual link. My concern is slightly broader: the equality S_mult = S_parton + constant is not even a direct consequence of the single-particle momentum-rescaling argument, because S_mult is an entropy over multiplicity counts rather than over the partonic transverse-momentum distribution. Additional numerical issues (N=0 exclusion, NBD deviations at high N, and unpropagated fit uncertainties) reinforce the need for a hadronization-level test. The paper’s cross-check with λ_σ = 0.329 ± 0.025 is useful supporting evidence, and the entropy-derived value is plausibly consistent with model-independent estimates, but it does not by itself validate the hadronization step. A generator/toy-shower test with a known input λ would settle whether the multiplicity entropy slope actually recovers the input intercept. This does not change the reader’s conditional verdict: the article should be accepted pending such a test or an explicit quantitative estimate of the hadronization correction. My agreement with the reader is partial because the crucial assumption is the identity of the two entropy functionals, with ⟨z⟩-independence being one part of that assumption rather than the whole of it.","tokens_in":11789,"tokens_out":7366,"duration_ms":72428,"concrete_test":"Generate events with a known scaling gluon distribution, e.g. Q_s(x) = (x0/x)^{λ/2} with input λ = 0.33, using PYTHIA 8 (or a dipole shower plus a fragmentation model) in the H1 Q²/y bins of Table I; construct charged-hadron P(N), compute S_mult = −Σ P(N)ln P(N), fit S_mult versus log(1/x) per Q² bin, and compare the recovered slope to the input λ. If the recovered slope differs by more than the quoted ±0.007, the S_mult = S_parton + constant identification fails; if it recovers λ, the concern is resolved. As part of the same test, check whether ⟨z⟩ or the correction 2ln(1/⟨z⟩) varies with x within a fixed Q² bin.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result Eq. (15) rests on identifying S_mult from the NBD-fitted multiplicity distribution with S_parton = −∫d²kT P log P of Eq. (8). The paper’s only bridge is local parton-hadron duality plus the statement S_hadron = S_parton + 2log(1/⟨z⟩), with ⟨z⟩ assumed independent of x. But S_mult is a functional of the count distribution P(N), while S_parton is a functional of the single-parton transverse-momentum distribution; the rescaling argument for single-particle spectra does not, by itself, prove that the Shannon entropy of P(N) differs from S_parton only by an x-independent constant. Even granting LHPD, the computed S_mult excludes the N=0 bin and renormalizes via α (Table I), and the paper concedes that one-component NBD deviates for N≳40, both of which alter the entropy and can bias the slope versus log(1/x). Thus the claim that the entropy slope 'directly provides λ' is an assumption about hadronization, not an established consequence of geometric scaling.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes that the Boltzmann-Gibbs entropy of final-state charged-hadron multiplicities in DIS, S_mult = -∑ P(N) log P(N), inherits the geometric scaling of the partonic transverse-momentum distribution, leading to S_mult = C + λ log(1/x). Using single-component NBD fits to H1 multiplicity data in four Q² bins, the authors extract λ_entropy = 0.322 ± 0.007, compare it with λ_σ = 0.329 ± 0.025 from a dipole-model fit to inclusive γ*p cross sections, and argue that scaling entropy is a model-independent way to determine the Pomeron intercept and detect geometric scaling.","tokens_in":11920,"tokens_out":4797,"duration_ms":41989,"significance":"If the identification of S_mult with the partonic scaling entropy can be justified, the result is valuable: it would provide a new observable, accessible from multiplicity measurements alone, that exhibits geometric scaling and yields a Pomeron intercept consistent with inclusive DIS scaling. The paper is careful to compare with existing model-dependent extractions and to place the result in the context of saturation-model predictions. The use of public H1 data and a transparent NBD parametrization makes the analysis easy to follow, although no code or covariance matrices are provided.","major_comments":[{"comment":"The quoted λ_entropy = 0.322 ± 0.007 is presented without an explicit error-propagation formula from the NBD fits in Table I. The slope of S_mult versus log(1/x) is fitted from four y-bins per Q² bin, but no covariance matrix or weighted-average prescription is given, and the low χ²/dof values in Table I (e.g., 0.018) suggest that the NBD parameter uncertainties do not account for bin-to-bin correlations. Please provide the full propagation from the NBD parameters to S_mult and to the final λ, including systematic uncertainties from the α renormalization and from the choice of χ² definition in Eq. (14).","section":"Section III, Eq. (15) and Fig. 3"},{"comment":"The identification of the multiplicity entropy S_mult with the partonic entropy S_parton is the load-bearing assumption, but it is only justified by local parton-hadron duality and by the statement S_hadron = S_parton + 2 log(1/⟨z⟩) with ⟨z⟩ independent of x. As the paper itself notes, the additive constant 'can have a dependence on Q²'; if ⟨z⟩ depends on x or Q², the slope of S_mult versus log(1/x) is contaminated and the extracted λ is biased. A quantitative hadronization-robustness test is required — for example, repeating the extraction with a modified ⟨z⟩(x) ansatz, with an alternative P(N) parametrization, or with N=0 included — before the claim that Eq. (15) 'directly provides λ' can be accepted.","section":"Section II, Eqs. (8)–(10)"},{"comment":"The fits exclude the N=0 bin and renormalize the NBD through α, and the paper acknowledges that one-component NBD deviates for N≳40 at small x. These cuts alter the value of S_mult differently across y-bins, so they can bias the slope λ even if each individual fit is acceptable. Please quantify the sensitivity of λ_entropy to (i) the N=0 exclusion, (ii) the high-N truncation, and (iii) the use of a two-component NBD, and state the resulting systematic error on Eq. (15).","section":"Section III, Fig. 1 and Table I"}],"minor_comments":[{"comment":"The phrase 'high-order corrections like log n(1/x)' should be written as log^n(1/x) or explained in words, since the current notation is ambiguous.","section":"Section III, paragraph after Fig. 2"},{"comment":"The sentence 'The fit for P (N ) versus N is perform for N >1' should read 'is performed for N > 1'.","section":"Fig. 1 caption"},{"comment":"The phrase 'Its important to note' should be 'It is important to note'.","section":"Section III, discussion after Fig. 6"},{"comment":"It would help to list the number of data points and the number of degrees of freedom for each NBD fit, since χ²/dof values below 0.05 are otherwise difficult to interpret.","section":"Table I"},{"comment":"The three panels use λ=0.29, 0.33, and 0.37, but the text says the minimum is λ_σ = 0.329 ± 0.025; please state explicitly whether λ=0.33 is the central value used in the scaling line and describe how the ratio data/theory was computed, including the normalization and the data selection.","section":"Section III, Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The core comparison is interesting and potentially publishable after revisions. The main risk is not the cross-check with the inclusive cross-section but the hadronization bridge between S_mult and S_parton, which the paper itself flags as dependent on ⟨z⟩ and on the validity of one-component NBD. I would ask for a quantitative sensitivity analysis and a proper uncertainty propagation before endorsing the quoted λ_entropy. No concerns about citation behavior or scope beyond the report."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The paper does something new: it takes the Shannon entropy of NBD-fitted H1 multiplicity distributions, plots it against log(1/x), and reads off the slope as the Pomeron intercept λ, getting 0.322 ± 0.007. It then compares that with λ from inclusive cross-section scaling in the same kinematic region, 0.329 ± 0.025, and finds agreement. That cross-check is the strongest part of the paper. If the entropy slope works, it gives you a cheap observable for detecting geometric scaling in multiplicity data alone, useful for LHC pp and pA contexts where inclusive cross-section scaling is harder to measure.\n\nThe soft underbelly is the bridge between hadron multiplicity entropy and partonic entropy. The paper says S_hadron = S_parton + 2 log(1/⟨z⟩), and assumes ⟨z⟩ doesn't depend on x. That's an assumption, not a derivation, and the stress-test note is right: S_mult is a functional of the counting distribution P(N), while S_parton is the Shannon entropy of a single-parton kT distribution. Unless you have a stronger hadronization argument than LHPD plus a constant shift, the slope can be contaminated by any x-dependence in ⟨z⟩ or in the shape of P(N). The paper is upfront that this is conditional on hadronization, so I don't call it circular, but 'model-independent' oversells it.\n\nThe statistics are a bit hand-wavy. The reported ±0.007 has no error budget. The χ2/dof values are suspiciously low (0.018 in one bin), which usually means the used uncertainties are too big or the fit isn't a meaningful chi-square. They also drop N=0 and renormalize with α, but don't show how that affects the entropy slope. Given ⟨N⟩ around 3-6 in several bins, P(0) isn't tiny; if it varies with x, the slope is biased. These are fixable with a sensitivity check and a proper propagation of NBD parameter uncertainties. The one-component NBD issue for N≳40 is acknowledged and probably minor here because the tail is small, but they still extrapolate the fit to compute entropy.\n\nBottom line: this is a worthwhile contribution to saturation phenomenology, not a finished measurement. The idea is sensible, the comparison is informative, and the paper is honest about its caveats. A serious referee should ask for: (1) a hadronization robustness test (e.g., try a different ⟨z⟩ behavior or Monte Carlo fragmentation), (2) uncertainty propagation, (3) N=0 sensitivity. If those come back clean, the result becomes much more convincing. Yes, send it to peer review.","headline":"A promising but methodologically loose extraction of the Pomeron intercept from multiplicity entropy; the agreement with cross-section scaling is real, but the error bar and hadronization assumptions need work before I'd trust the central value.","tokens_in":12549,"tokens_out":3881,"would_cite":true,"duration_ms":34734,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The slope of the scaling entropy from hadron multiplicities directly measures the Pomeron intercept, with $\\lambda = 0.322 \\pm 0.007$ from H1 data.","keywords":["geometric scaling","scaling entropy","Pomeron intercept","parton saturation","negative binomial distribution","hadron multiplicity","deep inelastic scattering","H1 data"],"falsifier":"Recompute the entropy $S_{\\rm mult} = -\\sum_N P(N)\\log P(N)$ directly from the binned H1 multiplicity data without the single-component NBD extrapolation to high multiplicities; if the slope of $S_{\\rm mult}$ versus $\\log(1/x)$ differs from $0.322 \\pm 0.007$ outside the quoted uncertainty, the high-multiplicity tail of the fit is carrying the result.","tokens_in":11479,"feed_emoji":"📈","tokens_out":8260,"duration_ms":68652,"temperature":0.7,"pith_summary":"The paper claims that the scaling entropy computed from final-state hadron multiplicity distributions in deep inelastic scattering is, under the scaling hypothesis, equal in slope to the partonic entropy, and that this slope directly gives the Pomeron intercept $\\lambda$. Analyzing H1 charged-hadron multiplicities with a negative binomial fit, the authors extract $\\lambda_{\\rm entropy} = 0.322 \\pm 0.007$, which agrees with $\\lambda_{\\sigma} = 0.329 \\pm 0.025$ obtained from the geometric scaling of the inclusive DIS cross section. If true, this makes the Pomeron intercept determinable from multiplicity data alone, without modeling the gluon distribution, and provides a model-independent probe of geometric scaling and parton saturation.","feed_headline":"Hadron multiplicity entropy gives Pomeron intercept 0.322","feed_subtitle":"The entropy slope in H1 data matches cross-section scaling, offering a model-independent probe of saturation.","key_machinery":"The machinery is the identity connecting partonic scaling entropy to hadron multiplicity entropy. Assuming the scaling form $P(x,k_T^2) \\sim x^{\\lambda} f(k_T^2/x^{\\lambda})$, the Tsallis/Boltzmann-Gibbs entropy of the partonic transverse momentum distribution becomes $S_{\\rm parton} = C + \\lambda \\log(1/x)$, so the slope in $\\log(1/x)$ is exactly $\\lambda$. Experimentally, the entropy is computed from the negative binomial distribution (NBD) fit to the measured charged-hadron multiplicity: $S_{\\rm mult} = -\\sum_N P(N) \\log P(N)$. Local parton-hadron duality is invoked to assert that hadronization adds only a constant shift, $S_{\\rm hadron} = S_{\\rm parton} + 2\\log(1/\\langle z\\rangle)$, leaving the slope unchanged. The slope extracted from the linear fit of $S_{\\rm mult}$ against $\\log(1/x)$ is the reported $\\lambda_{\\rm entropy}$, and the same $\\lambda$ parametrizes the saturation scale $Q_s^2(x) \\sim x^{-\\lambda}$.","core_discovery":"The central discovery is a direct link between the slope of the Boltzmann-Gibbs entropy of final-state hadron multiplicities and the Pomeron intercept. For a probability distribution that satisfies geometric scaling, $P(x,k_T^2) \\sim x^{\\lambda} f(k_T^2/x^{\\lambda})$, the partonic entropy grows as $S = C + \\lambda \\log(1/x)$; the paper shows that the entropy from the measured multiplicity $P(N)$ follows the same linear behavior in the H1 kinematic range, and fits its slope to obtain $\\lambda_{\\rm entropy} = 0.322 \\pm 0.007$. This value coincides with the $\\lambda$ obtained from the scaling of the inclusive $\\gamma^* p$ cross section ($\\lambda_{\\sigma} = 0.329 \\pm 0.025$), indicating that the two observables are equivalent probes of the same small-$x$ dynamics. The paper presents this as evidence that scaling entropy is a more efficient and more model-independent way to detect geometric scaling and determine the Pomeron intercept than the traditional cross-section analysis.","pith_inferences":["One extension the authors do not pursue is applying the same entropy-slope method to diffractive or exclusive data at a future electron-ion collider, where the scaling hypothesis has different kinematic reach; a deviation in $\\lambda$ would reveal where geometric scaling breaks down.","The assumption that $\\langle z\\rangle$ is independent of $x$ could be tested directly with hadronization-aware Monte Carlo generators that include intrinsic transverse momentum or shower effects; if the extracted $\\lambda$ shifts when $\\langle z\\rangle$ varies, the method would require a hadronization correction.","A stronger universal claim follows if the result holds: $\\lambda \\approx 0.32$ is the same in deep inelastic scattering and hadroproduction, so the entropy observable may serve as a cleaner comparator between HERA and LHC than the cross-section scaling used previously."],"forward_implications":["If the method is correct, the Pomeron intercept $\\lambda$ can be measured from multiplicity distributions alone, bypassing the need to model unintegrated gluon distributions.","The same entropy-slope analysis can be applied to $pp$ and $pA$ collisions at LHC energies, where geometric scaling is harder to isolate in cross sections.","The agreement between $\\lambda_{\\rm entropy}$ and $\\lambda_{\\sigma}$ suggests a universal small-$x$ scaling, providing a target that saturation models must reproduce.","Because the entropy slope is extracted from an integrated quantity, it is insensitive to high-$k_T$ details; differential $p_T$ spectra would be needed to expose those, a limitation the paper itself notes.","The method can be used to detect geometric scaling in data sets where the scaling variable is not known a priori, by checking for linearity of $S_{\\rm mult}$ versus $\\log(1/x)$."],"supporting_citations":[{"why":"Provides the power-law gluon distribution and scaling variable used in the inclusive cross-section fit, and the MPM model with $\\lambda\\approx 0.33$.","marker":"[8]"},{"why":"The H1 data set on charged-hadron multiplicities from which $P(N)$ and the experimental entropy are derived.","marker":"[38]"},{"why":"Gives the local parton-hadron duality statement that lets the partonic entropy slope be read off hadron multiplicities.","marker":"[53]"},{"why":"Introduces the scaling-entropy formalism connecting entropy to the saturation scale in hadronic collisions, the basis of Eq. (9).","marker":"[11]"},{"why":"A model-independent quality-factor analysis of DIS scaling with $\\lambda\\approx 0.32$, used as a consistency check.","marker":"[6]"},{"why":"Another model-independent scaling analysis with a compatible $\\lambda$, supporting the universality of the entropy result.","marker":"[7]"},{"why":"Supplies the NLO BK prediction $\\lambda\\approx 0.3$, a benchmark for the extracted Pomeron intercept.","marker":"[70]"}],"fun_headline_variants":["Entropy slope yields precise pomeron intercept","Hadron entropy matches scaling for pomeron intercept","Scaling entropy: a precise probe of pomeron intercept","Pomeron intercept pinned by multiplicity entropy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that transforming partons into hadrons adds only a constant to the entropy, independent of the momentum fraction $x$; if the average hadron momentum fraction varies with $x$ or $Q^2$, the slope measured from multiplicity data is contaminated and the extracted $\\lambda$ is biased.","fun_headline_variants_meta":{"raw":{"variants":["Entropy slope yields precise pomeron intercept","Hadron entropy matches scaling for pomeron intercept","Scaling entropy: a precise probe of pomeron intercept","Pomeron intercept pinned by multiplicity entropy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000378,"raw_usage":{"total_tokens":1972,"prompt_tokens":868,"completion_tokens":1104,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":1042}},"tokens_in":484,"tokens_out":1104,"duration_ms":8251,"temperature":1.0,"reasoning_tokens":1042,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:40:09.938481+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the entropy $S_{\\rm mult} = -\\sum_N P(N)\\log P(N)$ directly from the binned H1 multiplicity data without the single-component NBD extrapolation to high multiplicities; if the slope of $S_{\\rm mult}$ versus $\\log(1/x)$ differs from $0.322 \\pm 0.007$ outside the quoted uncertainty, the high-multiplicity tail of the fit is carrying the result.","supporting_citations":[{"cited_title":"Systematics of geometric scaling","cited_arxiv_id":"hep-ph/0610435","evidence_quote":"Provides the power-law gluon distribution and scaling variable used in the inclusive cross-section fit, and the MPM model with $\\lambda\\approx 0.33$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"A model-independent quality-factor analysis of DIS scaling with $\\lambda\\approx 0.32$, used as a consistency check."},{"cited_title":"Quantitative Study of Geometrical Scaling in Deep Inelastic Scattering at HERA","cited_arxiv_id":"1211.5305","evidence_quote":"Another model-independent scaling analysis with a compatible $\\lambda$, supporting the universality of the entropy result."}],"review_version":1}