{"id":"be727363-17d0-4aab-9981-e3112b1cdc4f","arxiv_id":"1909.00375","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Wounded quark sources with shadowing in a Glauber Monte Carlo model reproduce the centrality dependence of multiplicity fluctuations measured by NA49 in Pb+Pb collisions for 80 to 200 projectile participants.","lead":"This paper compares computer models of particle production in lead-lead collisions at CERN's SPS accelerator to explain why the event-to-event spread in particle counts changes with collision geometry. The authors find that a model with wounded quarks and shadowing matches the measured spread, while the standard wounded-nucleon model does not.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim depends on a per-quark source distribution from a pp fit with chi2/Ndof approx 6; if a better-fitting source distribution changes the 80<N_proj^p<200 scaled variance, the shadowing conclusion is not robust.","rationale":"I read the paper as an honest, exploratory GMC study. Its strongest claim is modest and explicitly qualified in Sec. VI, where the authors state that the shadowed WQM \"partially reproduces\" the centrality dependence and that peripheral discrepancies grow. The shadowing parameter lambda is a phenomenological parameter, which is legitimate; the fact that lambda is fixed using the average multiplicity rather than the variance is genuine independent support. The weakest link is that the model's only fixed input for particle emission, P_H, comes from a fit the authors themselves call \"rather poor\" (chi2/Ndof ~ 6), while the abstract claims it \"describe[s] quite well.\" Because the scaled variance in the shadowed WQM is dominated by source-number fluctuations, it is sensitive to the exact shape of P_H, which is not uniquely determined by a poor pp fit. The paper also provides no error bars or chi2 for the quoted 80 < N_proj^p < 200 variance agreement, so the central claim is not quantitatively pinned down. The reader's identified weakest assumption concerns independent superposition more broadly; my concern sharpens one concrete way that assumption can fail, namely a mis-specified P_H fixed from inadequate pp data. The proposed re-fit test would settle whether the shadowing conclusion is robust or an artifact of the poor pp fit. This does not move the verdict: the claim is already CONDITIONAL, and the concern supports that conditionality rather than overturning it.","tokens_in":10346,"tokens_out":18271,"duration_ms":191553,"concrete_test":"Re-fit the NA49 pp multiplicity distribution (Fig. 3) with an alternative single-source distribution, for example a two-component NB mixture, that achieves chi2/Ndof ~ 1; keep the wounded-quark Glauber part of GLISSANDO unchanged; fix lambda by matching the Pb+Pb average multiplicity as in Sec. V; and recompute the centrality dependence of the scaled variance. If the 80 < N_proj^p < 200 variance curves shift by more than the NA49 published errors from Ref. [16], the central claim depends on the poorly constrained pp fit; if the curves remain inside the data, the concern does not land.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing input is the per-source multiplicity distribution P_H fixed by the WQM fit to pp data in Sec. V. The authors report the fit as \"rather poor\" with chi2/Ndof ~ 6, in tension with the abstract's \"describe quite well.\" This P_H is then used unchanged in Pb+Pb, while the shadowing strength lambda is adjusted only to reproduce the average multiplicity. The scaled-variance curve in Fig. 6b is therefore a prediction only if this P_H is the correct source distribution. The fitted NB is nearly Poissonian (omega_s ~ 1.04), so the centrality rise of omega is generated essentially by the source-number fluctuation term, which is sensitive to the shape and parameters of P_H. A chi2/Ndof of about 6 means that alternative source distributions compatible with the pp data could change this term materially. In addition, no uncertainties or chi2 are given for the 80 < N_proj^p < 200 agreement in the right panel of Fig. 6, so the claim that the model \"reproduces\" the centrality dependence is asserted visually rather than established quantitatively.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the Glauber Monte Carlo approach with wounded nucleon and wounded quark sources to describe the centrality dependence of the scaled variance of charged-particle multiplicity measured by NA49 in Pb+Pb collisions at sqrt(s_NN)=17.3 GeV. The authors first show that the Wounded Nucleon Model, with source parameters fitted to the pp multiplicity distribution, reproduces the average multiplicity but fails for the scaled variance. They then introduce a Wounded Quark Model in which each wounded quark emits particles according to a negative binomial distribution fitted to pp data, and implement a shadowing suppression factor S(n,lambda)=exp(-n*lambda) (Eq. 7) to reduce the overproduction of average multiplicity in Pb+Pb. With lambda=0.95, the model is reported to reproduce (in Sec. VI, 'partially reproduces') the centrality dependence of the scaled variance for 80 < N_proj^p < 200, while deviating for more peripheral collisions. The paper concludes that subnucleonic sources with shadowing provide a better description of the NA49 fluctuation data than the standard wounded nucleon picture.","tokens_in":10636,"tokens_out":3104,"duration_ms":32187,"significance":"If the central claim holds, the paper gives a simple, physically motivated explanation for a long-standing puzzle: the non-trivial centrality dependence of multiplicity fluctuations observed by NA49, which standard event generators and the wounded nucleon model fail to describe. The strength of the paper is its transparent use of the compound-distribution formalism, Eq. (2), and the explicit demonstration that fluctuations in the number of sources, rather than in the source strength, drive the rise of the scaled variance; the pedagogical appendix makes this point clearly. The use of pp and Pb+Pb data at the same experimental acceptance is a further advantage. However, because the WQM fit to pp data is poor and the shadowing parameter is adjusted to the Pb+Pb average multiplicity, the prediction is not as clean as the abstract suggests. The central claim is visually supported but lacks quantitative goodness-of-fit assessment, so the paper in its current form is a suggestive phenomenological study rather than a definitive demonstration.","major_comments":[{"comment":"The central claim that the shadowed WQM 'reproduces' the centrality dependence of the scaled variance for 80 < N_proj^p < 200 is based on visual inspection of the right panel of Fig. 6. No chi-squared value, confidence band, or other quantitative measure is provided for this range, despite the fact that the model curve and data both have nontrivial structure. The authors should add a quantitative goodness-of-fit measure (e.g., chi^2/N_dof) for the quoted centrality range, or explicitly state that the agreement is only qualitative.","section":"Sec. V, Fig. 6b"},{"comment":"The WQM fit to the pp multiplicity distribution is reported as 'rather poor' with chi^2/N_dof ~ 6. The per-source negative binomial parameters (langle N_NB rangle = 0.53, k = 14) obtained from this poor fit are then used unchanged as the source distribution P_H in Pb+Pb. Since the scaled variance of Eq. (2) depends directly on the source-level variance omega_s, and since the fitted omega_s is close to 1.04 (nearly Poissonian), the centrality rise of omega is generated almost entirely by the source-number fluctuation term, which is sensitive to the shape of P_H. Alternative source distributions compatible with the pp data could therefore change the predicted scaled variance materially. The authors should assess this sensitivity, for example by repeating the Pb+Pb calculation with other P_H parameter sets that still describe the pp data within uncertainties.","section":"Sec. V, Fig. 3 and Eqs. (3)-(6)"},{"comment":"The shadowing parameter lambda = 0.95 is chosen to bring the model's average multiplicity into agreement with the Pb+Pb data (Fig. 6a). This means the model is not independent of the Pb+Pb data it describes; the scaled-variance prediction is conditional on a parameter tuned to the same dataset. The paper should show how the predicted scaled variance in Fig. 6b varies with lambda, or at least provide the range of lambda allowed by the average-multiplicity fit, so that the reader can judge whether the agreement in the range 80 < N_proj^p < 200 is robust or accidental.","section":"Sec. V, Eq. (7) and Fig. 6a"},{"comment":"The abstract states that the model 'reproduces the centrality dependence of scaled variance', while Sec. VI says the model 'partially reproduces' it and that 'the discrepancy between data and WQM predictions grows when going towards peripheral collisions'. Since the abstract is the central claim, it should carry the same qualification as the conclusions. The reader should not have to reach Sec. VI to discover that the claim is limited to 80 < N_proj^p < 200 and that the model fails for peripheral collisions.","section":"Abstract and Sec. VI"}],"minor_comments":[{"comment":"There is a typo in the first bullet: 'Wounded nuleons' should be 'Wounded nucleons'.","section":"Sec. VI"},{"comment":"In the paragraph describing the NA49 measurement, 'multiplicity distributions of charged articles' should read 'charged particles'.","section":"Sec. II"},{"comment":"The text says the WNM scaled variance shows a 'small monotonic increase' with decreasing participant number, but Fig. 2b appears to show a very weak increase; please clarify whether this is a monotonic trend within statistical fluctuations.","section":"Sec. IV, Fig. 2"},{"comment":"The vertical axis label 'QW' is not defined in the caption or text before the figure; define it as the average number of wounded quarks per nucleon in the caption.","section":"Sec. V, Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is a plausible phenomenological study, but the central claim is weakened by the poor WQM fit to pp (chi^2/N_dof ~ 6) and by the absence of a quantitative comparison in the key range. The authors should be encouraged to add a sensitivity study and to align the abstract with the qualified conclusion in Sec. VI. If those issues are addressed, the paper could be suitable for publication in a specialized journal; in its current form, the overclaim in the abstract should not pass unchallenged."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper applies an existing wounded-quark model with shadowing to the NA49 scaled-variance data and shows, honestly, that it does better than wounded nucleons for mid-central Pb+Pb. But the load-bearing pp fit is poor (chi2/Ndof ~ 6) and lambda is tuned to the average multiplicity, so the scaled-variance agreement is less robust than the abstract suggests.\n\nWhat's new: the demonstration that subnucleonic sources with shadowing can produce the observed centrality dependence of omega for 80 < N_proj^p < 200. This is a real result, even if the machinery is from Bozek et al. and Chatterjee et al. The paper is also appropriately honest: the bullets in Sec. VI say 'partially reproduces,' the peripheral mismatch is stated, and the appendix makes the compound-distribution point cleanly.\n\nSoft spots. First, the WQM fit to pp is called 'rather poor' in Sec. V but 'quite well' in the abstract. That is a factual inconsistency. Second, the per-source NB parameters come from that poor fit, and the scaled variance prediction is sensitive to the shape of P_H through the source-number-fluctuation term. With chi2/Ndof ~ 6, there is room for alternative source distributions that still fit pp but change the Pb+Pb omega. The stress-test note gets this right. Third, lambda=0.95 is adjusted to reproduce the average multiplicity, so the model is not parameter-free; the scaled variance is not directly fitted, which is good, but the test is weaker than an independent prediction. Fourth, the agreement in Fig. 6b for 80<N_proj^p<200 is asserted visually; no uncertainties or chi2 are given. That is a minor issue given the data, but it should be quantified.\n\nThe central claim as qualified is plausible. I would not call it fully established, but it is a fair phenomenological point that the source-number fluctuations from quark degrees of freedom help. The paper deserves a serious referee; it is not a desk reject. The referee should ask the authors to fix the fit-quality language, provide uncertainties or at least a band for the scaled variance from the pp fit parameters, and discuss how robust the shadowing conclusion is to alternative P_H shapes.\n\nMy verdict: conditional accept, with requested clarifications.","headline":"Modest, honest application of wounded-quark+shadowing to NA49 scaled variance; the claim holds qualitatively but the poor pp fit and fitted lambda leave real slack.","tokens_in":11175,"tokens_out":1990,"would_cite":false,"duration_ms":18328,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["02.50.Ey","05.10.Ln","12.40.Ee"],"model":"deepseek-v4-flash","headline":"Shadowed quark sources reproduce Pb+Pb multiplicity fluctuations where wounded-nucleon models fail.","keywords":["multiplicity fluctuations","scaled variance","Fano factor","Glauber Monte Carlo","wounded quark model","wounded nucleon model","shadowing","NA49"],"falsifier":"Measure the scaled variance in Pb+Pb collisions at fixed numbers of both projectile and target participants, binning tightly in both, at the same acceptance. If the variance is largely reduced once target participants are fixed, the source-number fluctuation term $\\omega_k$ is not the origin of the data; if it remains high, the quark-source fluctuation mechanism is supported.","tokens_in":10119,"feed_emoji":"⚛️","tokens_out":8332,"duration_ms":74073,"temperature":0.7,"pith_summary":"This paper asks whether the non-trivial centrality dependence of charged-particle multiplicity fluctuations measured in fixed-target Pb+Pb collisions at $\\sqrt{s_{NN}}=17.3$ GeV can be explained by the number and character of particle-emitting sources. It argues that the Wounded Nucleon Model, which fits the average multiplicity, fails for the scaled variance (Fano factor) of the multiplicity distribution. Treating the sources as wounded quarks, with the per-source distribution fixed by proton-proton data, overproduces the average multiplicity; adding shadowing, in which a quark source behind other sources emits less, brings both the average multiplicity and the scaled variance into agreement for central to mid-peripheral collisions. The claim is that subnucleonic sources with shadowing, rather than nucleon sources alone, are what generate the observed multiplicity fluctuations.","feed_headline":"Shadowed quark sources explain Pb+Pb multiplicity fluctuations","feed_subtitle":"Wounded-nucleon models miss the centrality trend; quark sources with shadowing match the measured Fano factor.","key_machinery":"The load-bearing object is the compound (superposition) multiplicity distribution: $N=\\sum_{i=1}^{N_p} n_i$, built from a Glauber Monte Carlo distribution of the number of sources $N_p$ and a per-source Negative Binomial distribution $P_H$ fixed by a fit to proton-proton data. The scaled variance splits as $\\omega=\\omega_s+\\langle N_s\\rangle\\omega_k$ (Eq. 2), so fluctuations in the number of sources multiply the per-source contribution. For the Wounded Quark Model the relevant sources are wounded quarks, meaning quarks inside nucleons that have collided at least once; shadowing modifies each source's emission by the factor $\\exp(-n\\lambda)$, where $n$ counts quark sources ahead of it in the same nucleus and $\\lambda=0.95$ is the fitted suppression parameter. This machinery converts a fixed participant number into a fluctuating source count and thereby generates the extra variance the data require.","core_discovery":"On the paper's own terms, the central discovery is that the centrality dependence of the scaled variance $\\omega$ of the charged-particle multiplicity distribution in Pb+Pb collisions at fixed projectile participant number is controlled by the number and fluctuation of subnucleonic sources, not by nucleon woundedness. The paper constructs compound distributions $P(N)$ from a Negative Binomial source distribution $P_H$ and a Glauber-determined distribution of sources $P_S$, so that $\\omega = \\omega_s + \\langle N_s\\rangle \\omega_k$. It finds that the Wounded Nucleon Model gives $\\omega$ nearly independent of centrality, contradicting data, while the Wounded Quark Model's additional quark-source fluctuations raise $\\omega$ at all centralities. The specific new result is that only after applying an exponential shadowing suppression $S(n,\\lambda)=\\exp(-n\\lambda)$ to quark sources, with $\\lambda=0.95$, does the model reproduce the measured centrality trend for $80 < N_{\\rm proj}^p < 200$; the shortfall at more peripheral collisions is left as an open trigger-process.","pith_inferences":["Inference: Because the formalism predicts the full shape of $P(N)$ rather than just its scaled variance, higher normalized moments such as skewness and kurtosis of the same data provide a sharper, currently unused test of the Negative Binomial per-source plus shadowing assumptions.","Inference: The shadowing parameter $\\lambda$ should be extractable from other observables, such as centrality-dependent mean transverse momentum or two-particle correlations; if different observables require different $\\lambda$, the mechanism would need revision.","Inference: The centrality window $80<N_{\\rm proj}^p<200$ where the model works suggests a natural experiment: measure multiplicity fluctuations in smaller collision systems at the same energy, where quark-counting effects and the proposed trigger process should move the predicted breakpoint."],"forward_implications":["The failure of the Wounded Nucleon Model for the scaled variance is read as evidence that particle production sources are subnucleonic already at SPS energies, not as a signal from a new final-state mechanism.","Fixing the per-source distribution from pp data, the model predicts both the average multiplicity and the scaled variance of Pb+Pb collisions in the range $80<N_{\\rm proj}^p<200$ once $\\lambda$ is set, with no additional centrality-dependent parameter needed there.","Beyond $N_{\\rm proj}^p<80$ the model systematically undershoots the data, so the paper points to an additional mechanism that increases the number of effective sources (sea quarks or gluons) only in the most peripheral range.","The extra variance from fluctuations in the number of quark sources is the reason $\\omega$ in the Wounded Quark Model exceeds the Wounded Nucleon Model value, as demonstrated by the paper's compound-distribution appendix."],"supporting_citations":[{"why":"This reference defines the wounded-nucleon model used as the baseline, which the paper shows fails for the scaled variance.","marker":"[10]"},{"why":"This reference supplies the pp and Pb+Pb multiplicity distributions and centrality-dependent scaled variances that the model must reproduce.","marker":"[16]"},{"why":"This reference provides the wounded-quark Glauber framework, including quark-quark collision geometry and the average number of wounded quarks per nucleon.","marker":"[15]"},{"why":"This reference supplies the exponential source-shadowing suppression factor $S(n,\\lambda)=e^{-n\\lambda}$ used to lower the average multiplicity.","marker":"[36]"},{"why":"This reference establishes wounded-quark scaling for multiplicity data, motivating quark sources over nucleon sources.","marker":"[7]"},{"why":"This reference describes the Glauber Monte Carlo generator used to simulate wounded nucleon and wounded quark configurations.","marker":"[34]"}],"fun_headline_variants":["Shadowed quarks fix Pb+Pb multiplicity variance","Quark sources with shadowing match Pb+Pb data","Why wounded nucleons fail: quark shadowing works","Nucleon model misses, quark shadowing nails Pb+Pb","Quark shadowing beats nucleons for Pb+Pb fluctuations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim rests on the assumption that each wounded quark emits particles independently according to the same probability law measured in proton-proton collisions, so that only the number of quark sources and their shadowing change with centrality; if sources are correlated, or if emission depends on local density, the predicted scaled variance can fail.","fun_headline_variants_meta":{"raw":{"variants":["Shadowed quarks fix Pb+Pb multiplicity variance","Quark sources with shadowing match Pb+Pb data","Why wounded nucleons fail: quark shadowing works","Nucleon model misses, quark shadowing nails Pb+Pb","Quark shadowing beats nucleons for Pb+Pb fluctuations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000573,"raw_usage":{"total_tokens":2732,"prompt_tokens":993,"completion_tokens":1739,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":1657}},"tokens_in":609,"tokens_out":1739,"duration_ms":65405,"temperature":1.0,"reasoning_tokens":1657,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:54:22.840147+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the scaled variance in Pb+Pb collisions at fixed numbers of both projectile and target participants, binning tightly in both, at the same acceptance. If the variance is largely reduced once target participants are fixed, the source-number fluctuation term $\\omega_k$ is not the origin of the data; if it remains high, the quark-source fluctuation mechanism is supported.","supporting_citations":[{"cited_title":"Bialas, W","cited_arxiv_id":null,"evidence_quote":"This reference defines the wounded-nucleon model used as the baseline, which the paper shows fails for the scaled variance."},{"cited_title":"Centrality Dependence of Charged Particle Multiplicity at Mid-Rapidity in Au+Au Collisions at sqrt(s_NN) = 130 GeV","cited_arxiv_id":"nucl-ex/0105011","evidence_quote":"This reference supplies the pp and Pb+Pb multiplicity distributions and centrality-dependent scaled variances that the model must reproduce."},{"cited_title":"GLISSANDO 3: GLauber Initial-State Simulation AND mOre..., ver. 3","cited_arxiv_id":"1901.04484","evidence_quote":"This reference supplies the exponential source-shadowing suppression factor $S(n,\\lambda)=e^{-n\\lambda}$ used to lower the average multiplicity."},{"cited_title":"Kharzeev, E","cited_arxiv_id":null,"evidence_quote":"This reference establishes wounded-quark scaling for multiplicity data, motivating quark sources over nucleon sources."}],"review_version":1}