{"id":"11531b3c-6c48-4486-9939-315db5ce9d58","arxiv_id":"2502.08798","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A systematic library of modified hadronic simulations shows that matching Auger muon and Xmax data is possible only with a combination of cross-section, multiplicity, and elasticity changes that conflicts with other measurements.","lead":"This paper builds a library of 750,000 simulated cosmic-ray air showers with deliberately tweaked hadronic interaction parameters. It shows that current models can be adjusted to match some Pierre Auger data, but a fully consistent description remains out of reach, which matters for designing next-generation observatories.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Section 3's pure-superposition ansatz for iron (each nucleon at E/A) is the load-bearing assumption for the proton-iron merit factors and for the claim that only all three modified parameters fit the Auger data; its impact is untested.","rationale":"The reader's weakest-assumption analysis points to the pure-superposition treatment of nucleus-air interactions, and I concur. This is the point where the paper's quantitative claims for iron primaries and for the mixed-composition Auger comparison are least secure. The statement in Section 3 that 'the effects of modifications for iron primaries are always smaller than for proton primaries' is a direct consequence of the 1/A energy assignment, and the merit-factor correlations in Figs. 2 and 3 inherit this consequence. If the prescription were replaced by a Glauber-based or collective treatment, the iron-shower response to a given nucleon-level modification could change nonlinearly. Because the Auger [8] data are not proton-only, the claimed necessity of changing all three parameters could depend on this choice. The authors are fully transparent: they flag the assumption in the text and in the Conclusions, so this is not a hidden flaw. The proton-only conclusions (the p-air cross-section tension [10] and the Xmax-fluctuation tension [11]) do not depend on this assumption and support the 'consistent description remains elusive' message. The same transparency, however, means the paper should not be read as establishing the iron-mediated claims beyond the chosen ansatz; the proposed Glauber/alternative-energy test would settle the sensitivity. This does not change the reader's CONDITIONAL verdict: the paper is a useful resource, but its quantitative statements need either public data on the fit or an explicit sensitivity analysis to the superposition assumption.","tokens_in":6575,"tokens_out":11632,"duration_ms":117044,"concrete_test":"Generate a small subset of iron showers (e.g., the 25-point grid in sigma/N/eta that brackets the claimed Auger-agreement region) with two alternative prescriptions: (a) Glauber-Gribov effective energy per nucleon instead of E/A, and (b) full proton-level modification applied to every nucleon as an upper bound. Recompute the fractional change in N_mu for iron, the Fe/p merit factors of Fig. 2 (at 49 degrees and 1000 m), and the qualitative fit region for Auger [8]. If the parameter region requiring all three changes is reproduced under (a) or (b), the superposition assumption is not load-bearing; if the region moves by more than a grid step or a two-parameter solution appears, the iron-mediated conclusions must be restricted to the superposition ansatz.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 implements nucleus-air interactions as a pure superposition of nucleon-air interactions (explicitly stated at the start of Section 3 and in the Conclusions), with each nucleon taken at 1/A of the primary energy. This is not a neutral detail: the paper uses it to conclude that 'the effects of modifications for iron primaries are always smaller than for proton primaries', and that conclusion drives the proton-iron merit-factor correlations in Figs. 2 and 3. The Auger [8] data used for the headline fit are composition-mixed, so the finding that only a triple change (decreased sigma, increased N and eta) describes the data is partly determined by how the modified nucleon interactions propagate to iron showers. If Glauber shadowing, leading-nucleon energy sharing, or collective effects make the iron response larger or nonlinear in the modification parameters, the acceptable region in parameter space could shift and the 'only with changes to all three parameters' statement could fail. The authors explicitly flag the assumption, so this is not a hidden error; but the paper does not quantify how much the conclusions depend on it. The proton-only tensions with [10] and [11] are unaffected, so the 'consistent description remains elusive' message survives; the risk is confined to the iron-mediated claims.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This conference paper describes a large CORSIKA simulation library in which the high-energy hadronic interaction parameters (cross-section, multiplicity, and elasticity) of Sibyll 2.3d are modified in an energy-dependent, accelerator-consistent way. For proton and iron primaries at a single primary energy of 10^18.7 eV and several zenith angles, the authors study how ground-level observables, Xmax, and proton/iron separation merit factors respond to the modifications. Their main physics claim is that the Pierre Auger Observatory data of Ref. [8] can be reproduced by a modification of Sibyll 2.3d only if all three parameters are changed together (decreased cross-section, increased multiplicity, increased elasticity), and that this required combination is already in tension with other Auger measurements. The remaining paper is devoted to illustrating how the library can quantify modeling uncertainties for future observatory designs, focusing on distance-dependent effects and proton/iron discrimination.","tokens_in":6901,"tokens_out":3605,"duration_ms":39776,"significance":"If the main claim is substantiated, the paper provides a useful demonstration that the UHECR muon and Xmax tensions can be simultaneously accommodated by a retuned hadronic model, but only at the price of new tensions, reinforcing the conclusion that a fully consistent description remains elusive. The simulation library itself is a potentially valuable community resource: it systematically spans 75 modified interaction scenarios with a fully three-dimensional treatment, and the authors are explicit that their parameter choices respect accelerator constraints and avoid abrupt energy dependence. The paper also honestly flags its main modeling assumption about nuclear projectiles. However, the paper does not actually show the comparison to the Auger data that supports its central claim, reports no statistical uncertainties, and presents only a single primary energy and a single baseline model. These omissions make the significance conditional on the availability of the underlying comparison, which would need to be documented for the claims to be independently assessable.","major_comments":[{"comment":"The central claim that \"the hadronic model Sibyll 2.3d can indeed be modified to describe this particular Auger data well within our set parameter space, but only with changes to all three parameters\" is not supported by any displayed comparison. No figure or table shows the Auger data of Ref. [8] together with the modified simulations, no goodness-of-fit statistic is given, and the reader is not told which observables, zenith angles, or composition assumptions enter the comparison. This statement drives the abstract and the conclusions, so it must be documented quantitatively, for example with a plot of the allowed/excluded parameter combinations in the (σ, N, η) space.","section":"§1"},{"comment":"The iron/proton merit-factor results and the statement that \"the effects of modifications for iron primaries are always smaller than for proton primaries\" rest on the pure-superposition assumption, with each nucleon treated at 1/A of the primary energy. The paper explicitly acknowledges this assumption, but it does not test or bound its impact. If nuclear effects such as Glauber shadowing, leading-nucleon energy sharing, or collective behavior make the iron response nonlinear in the modification parameters, the correlations in Figs. 2 and 3 could change, and the conclusion that only a triple modification fits the Auger data would need re-evaluation. The proton-only conclusions would survive, but the composition-related claims require either a test against an alternative nuclear treatment or a stronger caveat stating which conclusions are robust.","section":"§3"},{"comment":"All quantitative statements are made at a single primary energy (10^18.7 eV) with 1000 showers per setting, and no statistical uncertainties are reported on mean profiles, merit factors, or relative changes. Since the text itself notes that shower-to-shower fluctuations are generally larger than differences between modifications, the reader cannot judge whether the reported trends, such as the 10% muon increase or the tight correlations in Fig. 2, are significant. Error bars or confidence intervals should be added at least for the quantities that support the main claims.","section":"§1 and Figs. 1–3"},{"comment":"The parameter space is not defined in this paper: the threshold energy, the f19 values for cross-section, multiplicity, and elasticity, and the exact 75 combinations are deferred to previous works. Without a table or appendix listing these values, the phrases \"well within our set parameter space\" and \"reasonable limits\" cannot be checked, and the claim that the modifications stay within accelerator constraints is not self-contained in this manuscript.","section":"§1"}],"minor_comments":[{"comment":"Typographical errors should be corrected, including \"the none of current hadronic interaction models\" and \"for for D_X > 200 g/cm2.\"","section":"§1"},{"comment":"The figures are dense: in Fig. 2 the color coding for zenith angles and parameter values is difficult to separate, and the baseline (unmodified) point is not marked. Please mark the unmodified model explicitly and consider enlarging the legend or using a separate panel for the baseline.","section":"Figs. 2 and 3"},{"comment":"The sentence stating that the EM energy density and charged-particle number density are conserved \"within a few percent\" at fixed D_X is only supported by visual inspection of Fig. 1; the relevant numerical values should be given, including the range of D_X over which this statement holds.","section":"§2"},{"comment":"Equation (1) defines the merit factor, but the text should clarify whether the quoted merit factors use the mean and standard deviation over simulated showers only, or also include systematic spreads due to the different modifications themselves.","section":"§3"},{"comment":"The conclusion would benefit from explicitly separating the results that are robust (proton-only muon and Xmax tensions) from those that depend on the superposition assumption for iron, since the text currently mixes the two levels of robustness.","section":"§4"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings-style paper whose main novelty is the systematic library and the qualitative insight that a simultaneous description of the Auger data requires a triple modification. The missing comparison to the Auger data is the key gap: a reader cannot verify the central claim from the manuscript alone. The superposition caveat is honestly stated, but its consequences for the iron-related conclusions deserve a sensitivity test or at least a much more prominent caveat. I would not reject the paper, because the library construction and the proton-only tension discussion are useful and likely correct; however, the presented evidence for the headline claim is currently insufficient."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jan, here's my read. The genuinely useful part of this paper is the 3D simulation library: 75 modified hadronic settings, proton and iron, five zenith angles, 750k CORSIKA showers, with the modifications tied to accelerator constraints. The authors show that 3D geometry matters for detector response—the sensitivity to cross-section versus multiplicity changes with core distance and zenith—and that the usual three-model spread (Sibyll/EPOS/QGSJET) may not bracket the truth. That is a real contribution for anyone estimating model uncertainty in future observatory design.\n\nThe paper is also honest. It flags the superposition assumption for nuclear primaries, and it does not pretend the modifications are unique. The main physical claim—that a retuned Sibyll 2.3d can describe the Auger hybrid data [8] only by decreasing cross-section and increasing both multiplicity and elasticity, and that this retuning conflicts with the proton-air cross-section [10] and Xmax fluctuations [11]—is plausible and consistent with the cited prior papers.\n\nThe soft spots are real but proportionate. First, this proceedings text does not show the actual comparison to the Auger data or any statistical uncertainties; the 'only with changes to all three' statement is asserted rather than demonstrated here. Second, the iron-shower conclusions rest on the superposition ansatz, each nucleon at E/A. The stress-test note is correct that this is load-bearing for the merit-factor correlations and for the composition-mixed fit; the authors flag it but do not quantify how the conclusions shift under alternatives like Glauber shadowing or energy sharing. Third, the library covers a single primary energy (10^18.7 eV) and one baseline model, so extrapolating to other energies is an assumption.\n\nNone of this sinks the paper. The 'consistent description remains elusive' message survives even if the iron-mediated claims are wrong, because the proton-only tensions with [10] and [11] are independent of the superposition ansatz. For a proceedings contribution, the scope is right and the claims are modest. What is missing is that the library itself is not public and the parameters are only summarized by reference to earlier papers; that limits reproducibility.\n\nWho should read it: people designing the next-generation UHECR observatory and anyone who wants a practical alternative to the three-model uncertainty bracket. I'd send it to peer review—it deserves referee time—but as a journal article it would need the actual data comparison and error bars, not as a proceedings note. My recommendation: engage with the library, cite it for the merit-factor analysis, and treat the superposition caveat as an explicit limitation rather than a hidden flaw.","headline":"Useful 3D simulation library with honest caveats; the iron-mediated conclusions rest on a flagged but untested superposition assumption.","tokens_in":7417,"tokens_out":4219,"would_cite":true,"duration_ms":39668,"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":"A retuned hadronic model can match Auger's muon and Xmax data, but only with simultaneous changes to cross-section, multiplicity, and elasticity, and those changes strain against other measurements.","keywords":["ultra-high-energy cosmic rays","air showers","hadronic interaction models","muon content","depth of shower maximum","Pierre Auger Observatory","CORSIKA simulations","Sibyll 2.3d"],"falsifier":"A direct test is to compare the best-fit modified model's predictions for the energy dependence of the muon signal and for the slope and fluctuations of the $X_{\\rm max}$ distribution against Auger's full energy scan: if no single set of the three parameters reproduces all of these simultaneously, the claimed simultaneous description fails. A future measurement of the proton-air cross-section at around $10^{18.7}$ eV that rules out the fitted value would also falsify the central finding.","tokens_in":6299,"feed_emoji":"🔭","tokens_out":6306,"duration_ms":61781,"temperature":0.7,"pith_summary":"The paper asks whether the discrepancies between current hadronic interaction models and ultra-high-energy cosmic-ray observations can be removed by modest, energy-dependent retuning of ordinary hadronic physics. Using a library of 750,000 simulated showers built from 75 combinations of changes to cross-section, multiplicity, and elasticity applied to the Sibyll 2.3d model, the authors find that the Pierre Auger Observatory's hybrid data can indeed be reproduced, but only by decreasing the cross-section and increasing both multiplicity and elasticity. That specific combination already conflicts with Auger's own proton-air cross-section measurement and with measured fluctuations in shower maximum, so a fully consistent description of all observations remains out of reach. The library itself is the lasting product: it offers a realistic spread of hadronic-model freedom for quantifying systematic uncertainties and for testing the design of future observatories.","feed_headline":"A single retuned hadron model fits Auger's muon and Xmax data","feed_subtitle":"The fit requires lower cross-section plus higher multiplicity and elasticity, straining against other measurements.","key_machinery":"The central object is a parameterized modification scheme applied inside the CORSIKA simulation code to the hadronic model Sibyll 2.3d. Modifications to elasticity, multiplicity, and cross-section take effect above an energy threshold and grow logarithmically with energy, reaching a chosen value at $10^{19}$ eV; 75 combinations of these settings were simulated for proton and iron primaries at $10^{18.7}$ eV across five zenith angles, giving 750,000 showers. For nuclear projectiles, the scheme assumes superposition: an iron nucleus is treated as independent nucleons, each modified at $1/56$ of the primary energy, which automatically makes modification effects smaller for iron than for proton primaries. This machinery lets the authors compare modified predictions directly with Auger observations and drives their distance-dependent, detector-geometry-dependent conclusions about muon counts and proton-iron separation power.","core_discovery":"The authors discover that, within their predefined parameter space, the hadronic model Sibyll 2.3d can be modified to describe Auger's simultaneous measurements of the muon content at ground and the depth of shower maximum, but only if all three interaction parameters are changed together: cross-section down, multiplicity up, and elasticity up. The increase in elasticity is the hardest to reconcile with other data: Auger's proton-air cross-section measurement constrains the cross-section-elasticity plane such that a lower cross-section should be accompanied by a lower elasticity, and a higher elasticity enlarges the shower-to-shower fluctuations of $X_{\\rm max}$ for protons beyond what Auger's profile measurements appear to allow. The paper therefore concludes that no combination in the library gives a fully consistent description of all UHECR observations, while the library itself provides a realistic representation of the freedom in modeling hadronic interactions at ultra-high energies.","pith_inferences":["One extension the authors do not pursue is to use the library to build a model-independent mapping from measured ground-signal ratios to primary mass; the strong correlations they report at 1000 m suggest such a mapping could be more robust to hadronic uncertainty than absolute muon counts.","The finding that muon counts can shift by up to 10% at fixed $DX = X - X_{\\rm max}$ while the electromagnetic energy density stays nearly fixed implies that future arrays with different detector spacing may infer different compositions from the same cosmic-ray flux, providing a cross-check of the library's assumptions.","If the superposition assumption for nuclei is relaxed in later work, the proton-only results would survive but the iron-based conclusions and merit factors would need to be recomputed; a full treatment of nuclear effects is the natural next test.","Because the parameter space explored here is deliberately mechanism-agnostic, the best-fit region can serve as a target for more specific physical models that aim to explain the same Auger data through collective effects or modified hadronization."],"forward_implications":["If the combined modification is correct, the long-standing muon deficit and the $X_{\\rm max}$ discrepancy can be produced by ordinary hadronic physics retuned within accelerator bounds, without invoking dramatically new physics at ultra-high energies.","The three standard models (Sibyll 2.3d, EPOS-LHC, QGSJET-II-04) may all lie on the same side of the data, so using their spread as the modeling uncertainty could underestimate the true range of possibilities.","For a future observatory, proton-iron separation power depends on the detector type and on its distance from the shower axis; quantities such as muon-to-electromagnetic ratios behave differently at 500 m and 1000 m, so design choices can be tested quantitatively with this library.","A consistent description of UHECR showers remains elusive: the modifications that fit the hybrid Auger data are already in tension with the proton-air cross-section and $X_{\\rm max}$ fluctuation measurements.","The strong correlations between the iron-proton merit factor and the change in muon number for protons mean that composition measurements and hadronic-model uncertainties cannot be treated independently in future analyses."],"supporting_citations":[{"why":"Provides the Auger hybrid data on muon content and $X_{\\rm max}$ that the modified models must reproduce.","marker":"[8]"},{"why":"Defines the Sibyll 2.3d baseline hadronic model that the paper modifies.","marker":"[9]"},{"why":"Supplies the macroscopic method of applying energy-dependent modifications to elasticity, multiplicity, and cross-section.","marker":"[5]"},{"why":"CORSIKA is the Monte Carlo code used to generate the three-dimensional air-shower library.","marker":"[6]"},{"why":"The Auger proton-air cross-section measurement constrains the cross-section-elasticity plane and conflicts with the required parameter combination.","marker":"[10]"},{"why":"The Auger measurement of $X_{\\rm max}$ fluctuations provides the constraint that increased elasticity may violate.","marker":"[11]"},{"why":"Earlier systematic exploration by the same authors establishes the parameter space and its effects on air-shower observables.","marker":"[1]"}],"fun_headline_variants":["Hadron tweaks fit Auger but clash with other UHECR data","No consistent hadron model for all UHECR observations yet","Auger muon and Xmax force cross-section down, multiplicity up","Modified hadron interactions: fit Auger, but elasticity strains"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that an iron nucleus hitting air behaves exactly as a superposition of independent modified nucleon-air interactions, with each nucleon carrying one fifty-sixth of the primary energy; if collective nuclear effects break this picture, the iron-related conclusions and merit-factor results would change, while the proton-only conclusions would remain intact.","fun_headline_variants_meta":{"raw":{"variants":["Hadron tweaks fit Auger but clash with other UHECR data","No consistent hadron model for all UHECR observations yet","Auger muon and Xmax force cross-section down, multiplicity up","Modified hadron interactions: fit Auger, but elasticity strains"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000866,"raw_usage":{"total_tokens":3749,"prompt_tokens":934,"completion_tokens":2815,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":2749}},"tokens_in":550,"tokens_out":2815,"duration_ms":21817,"temperature":1.0,"reasoning_tokens":2749,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T23:38:12.853292+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test is to compare the best-fit modified model's predictions for the energy dependence of the muon signal and for the slope and fluctuations of the $X_{\\rm max}$ distribution against Auger's full energy scan: if no single set of the three parameters reproduces all of these simultaneously, the claimed simultaneous description fails. A future measurement of the proton-air cross-section at around $10^{18.7}$ eV that rules out the fitted value would also falsify the central finding.","supporting_citations":[{"cited_title":"Blazek, J","cited_arxiv_id":null,"evidence_quote":"Earlier systematic exploration by the same authors establishes the parameter space and its effects on air-shower observables."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"CORSIKA is the Monte Carlo code used to generate the three-dimensional air-shower library."}],"review_version":1}