{"id":"5a8f23b2-e064-4e9a-a53b-1abc0adf0ce4","arxiv_id":"1908.09668","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The low-muon-number tail slope in simulated air showers is presented as a direct observable for the high-energy neutral-pion spectrum of the first proton-air interaction.","lead":"Using computer simulations, this paper shows that the low-muon-number tail of cosmic-ray air showers is set by the first proton-air collision. The slope of that tail could let observatories measure neutral-pion production at energies beyond the LHC, offering a new handle on the muon puzzle.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The abstract's 'direct link' between the muon-tail slope Λ_mu and the neutral-pion spectrum is not established: the Sec. 3.2 simulation changes total E_had and π0 shape simultaneously, and the π0 mapping is admitted to be model-dependent.","rationale":"The reader's weakest assumption was that the calibration between Λ_mu and Λ_had is universal and usable on data, and that the pion-level mapping retains model dependence. My concern sharpens this: even if the Λ_mu-to-Λ_had calibration were perfectly universal, the further step to the neutral-pion spectrum is confounded because the only demonstrated sensitivity is to the total hadronic/electromagnetic energy split. The Sec. 3.2 simulation changes the π0 spectrum and the energy split simultaneously, so it cannot distinguish a genuine spectral-shape sensitivity from a trivial energy-balance effect. This is a real soft spot in the central 'direct link' claim, and it is located in the paper's own admission that the π0-level function is model-dependent and preliminary. However, the paper is an ICRC proceedings contribution that explicitly frames the π0 connection as an initial study, and the qualitative program—using the low-N_mu tail to access first-interaction energy flow—is plausible and supported by the Sec. 3.1 simulations. The correct verdict remains CONDITIONAL: the abstract should be tempered to say the slope is sensitive to the first-interaction hadronic/electromagnetic energy split and, with additional model-controlled work, may constrain the high-energy π0 spectrum. My concern does not move the verdict away from CONDITIONAL, so UNCHANGED is appropriate. The proposed twin-spectrum test would settle whether the claimed spectral link is real or an artifact of the energy-split confound.","tokens_in":6251,"tokens_out":4580,"duration_ms":51362,"concrete_test":"Using Sibyll 2.3c, generate a reference set of proton showers at E = 10^19 eV and θ = 67 degrees. Construct two modified π0 production scenarios: (A) suppress high-x_L π0 production as in Fig. 5, and (B) reweight the π0 spectrum so that the per-shower distribution of E_had/E0 (or α1) is statistically identical to the reference while the high-x_L shape is changed. Compare the reconstructed Λ_mu among reference, A, and B. If B shows no significant change in Λ_mu, the muon-tail slope is sensitive only to the energy split, not to the π0 spectrum shape. Additionally, repeat scenario A with EPOS-LHC and QGSJET-II-04 to quantify the model dependence admitted in Sec. 3.2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the ground-level muon tail slope Λ_mu carries information about the shape of the high-energy neutral-pion spectrum, not merely about the total energy split E_had/E0 in the first interaction. Section 3.1 establishes, at best, a model-independent connection between Λ_mu and the slope of the E_had/E0 distribution, Λ_had. Section 3.2 then argues that because E_em = 1 - E_had is fed mainly by π0, Λ_had is connected to the π0 spectrum. But E_had is a single number per shower: its complement is the total electromagnetic energy, not the spectral shape of the pions. The only simulation shown in Sec. 3.2 modifies the inclusive π0 cross-section by suppressing high-x_L production; this necessarily changes the total electromagnetic energy and hence E_had/E0 and α1 (Eq. 2.1). The observed change in the N_mu tail (Fig. 5) is therefore fully compatible with Λ_mu tracking only the hadronic/electromagnetic energy split, with no sensitivity to the shape of the π0 spectrum beyond its contribution to that split. Moreover, the paper itself states that the derived monotonic function relating Λ_mu to the π0 spectrum 'has some dependence on the details of the hadronic interaction models,' while the abstract claims a 'direct link' without that caveat. Thus the load-bearing step from a measured slope to the first-interaction π0 energy spectrum is not supported by the presented evidence; the supported claim is sensitivity to the first-interaction energy flow into the hadronic sector.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper claims that the slope of the low-number tail of the muon-number distribution, Λ_mu, measured at the ground in ultra-high-energy cosmic-ray air showers is a direct probe of properties of the first hadronic interaction, specifically the fraction of primary energy transferred to the hadronic component and, more strongly, the high-energy tail of the neutral-pion energy spectrum. The study uses CONEX simulations of proton-induced showers at 10^19 eV and 67° zenith angle with a 1 GeV muon threshold. Section 3.1 establishes a simulation-based calibration between Λ_mu and the slope Λ_had of the first-interaction hadronic-energy fraction distribution, and shows that Λ_mu can be extracted in a mixed-composition scenario with 20% muon-number smearing. Section 3.2 reports that suppressing high-x_L neutral-pion production in SIBYLL 2.3c changes the N_mu tail, and the authors argue that this connects Λ_mu to the pion spectrum. The paper also discusses experimental precision and connections to LHC forward measurements.","tokens_in":6582,"tokens_out":3281,"duration_ms":35347,"significance":"If the claimed connection holds, this observable would provide a new way to access first-interaction physics at center-of-mass energies around 100 TeV, beyond current accelerators, using existing or planned cosmic-ray observatories such as the Pierre Auger Observatory. The paper has notable strengths: the simulation setup is clearly stated, statistical bands are shown, the mixed-composition and detector-smearing analysis is a useful feasibility check, and the comparison across three post-LHC hadronic models is informative. The central limitation, however, is that the evidence supports sensitivity to the first-interaction hadronic/electromagnetic energy split, not yet a model-independent sensitivity to the neutral-pion spectral shape. The abstract's phrase 'direct link' overstates what the presented simulations establish.","major_comments":[{"comment":"The simulation that suppresses high-x_L neutral-pion production changes not only the spectral shape of the π0 distribution but also the total electromagnetic energy fraction E_em = 1 - E_had, and hence α1 in Eq. (2.1). The observed change in the N_mu tail is therefore fully compatible with Λ_mu tracking only the hadronic/electromagnetic energy split, with no demonstrated sensitivity to the shape of the π0 spectrum beyond its contribution to that split. To support the central claim, the authors should show a control test where the π0 spectral shape is varied while E_had/E0 (or α1) is held fixed, or provide a quantitative decomposition showing that Λ_mu carries information about the shape independently of the first moment.","section":"Sec. 3.2, Fig. 5"},{"comment":"The calibration between Λ_mu and Λ_had is obtained by selecting simulated showers from a large ensemble and fitting the response, rather than from a physical derivation or an out-of-sample test. The claim that the relation is independent of hadronic models is supported only by three post-LHC models that share common assumptions about the leading-particle and energy-flow behavior. The figure shows model-dependent lines, so the 'independently of the hadronic interaction models' assertion requires a stronger demonstration, for example a model with a deliberately different first-interaction energy-flow distribution or a cross-validation on independent simulation sets.","section":"Sec. 3.1, Fig. 3 (left)"},{"comment":"The abstract claims that the slope of the low-N_mu tail is 'a direct link to the high energy spectrum of neutral pions,' but the text explicitly states that the derived function relating Λ_mu to the π0 spectrum 'has some dependence on the details of the hadronic interaction models.' This is an internal inconsistency in the level of claim. Either the abstract should be softened to reflect the model-dependent calibration, or the paper should provide an estimate of the systematic uncertainty and demonstrate that the model dependence does not affect the qualitative conclusion.","section":"Abstract vs. Sec. 3.2, final paragraph"}],"minor_comments":[{"comment":"There is an encoding artifact in the text: 'protons.¢aMoreover' should read 'protons. Moreover'. Also, 'the slope of theNµ' is missing a space.","section":"Sec. 4"},{"comment":"The caption reads 'Conversion between Λµ and Λα,' but the text and axes use Λ_had; please make the notation consistent.","section":"Fig. 3 (left), caption"},{"comment":"The symbols m and m_tot are not defined precisely enough; it is unclear whether they refer to charged multiplicity, total hadronic multiplicity, or something else. Please clarify in the text.","section":"Eq. (2.2)"},{"comment":"The paper uses the words 'prove' and 'proved' for statements obtained from a finite set of simulations; 'demonstrate' or 'show' would be more accurate and appropriate for a proceedings contribution.","section":"Introduction and Conclusions"}],"recommendation":"major_revision","confidential_remarks":"This is an ICRC proceedings paper, so the length and level of proof are constrained. The main scientific concern is the gap between the abstract's 'direct link' to the neutral-pion spectrum and the actual simulation evidence, which conflates spectral shape with the total electromagnetic energy fraction. This is fixable within the scope of the paper by adding a control simulation or by revising the claim in the abstract and conclusions. I would not reject, but the current overstatement is load-bearing for the paper's stated central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nShort version: this is a competent proceedings paper whose abstract overreaches. The new pieces are the explicit interpretation connecting the muon-number tail slope to the neutral-pion spectrum and a feasibility study with mixed composition and smearing. What the paper actually shows cleanly is that Lambda_mu tracks the first-interaction hadronic energy fraction E_had/E0 across three post-LHC models, and that this tail can be extracted even with 20% smearing and an admixture of heavier primaries. That part is useful and reasonably solid.\n\nThe soft spot is the step from Lambda_had to the pi0 spectrum. In Sec. 3.2 the authors modify SIBYLL by suppressing high-x_L pi0 production. That necessarily changes the total electromagnetic energy, hence E_had/E0 and alpha1. So the observed shift in the N_mu tail is fully compatible with Lambda_mu responding only to the hadronic/electromagnetic energy split, not to the shape of the pion spectrum. The paper itself concedes in Sec. 3.2 that the mapping function 'has some dependence on the details of the hadronic interaction models,' which undercuts the abstract's 'direct link.' For a single number E_had, its complement gives only the total electromagnetic energy, not the spectral shape. So the load-bearing claim is not established by the simulation as presented.\n\nMinor: the tail fit range is unspecified, and no code or data are provided; that is typical for a proceedings preview but still limits reproducibility. The suggested cross-check with LHC forward experiments is well placed and worth pursuing.\n\nWho is this for? People working on muon production in UHECR showers and the muon puzzle. They can cite the calibration and feasibility plots, but should not repeat the 'direct link' claim without qualification. If this were submitted to a journal, I would send it to a referee with a request to temper the abstract and add a test that changes the pi0 spectrum shape while keeping E_had/E0 fixed. As an ICRC proceedings contribution it is acceptable, but the abstract should not be the takeaway.\n\nRecommendation: engage with the paper and flag the overclaim. It deserves referee time as a full paper only if the required discriminating test is added.","headline":"Competent proceedings paper with a useful calibration result, but the abstract's 'direct link' to the neutral-pion spectrum overstates what the simulation actually shows.","tokens_in":7145,"tokens_out":2418,"would_cite":false,"duration_ms":25944,"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 low muon-number tail at the ground is a direct probe of the high-energy neutral-pion spectrum of the first proton-air collision.","keywords":["ultra-high energy cosmic rays","extensive air showers","muon number fluctuations","neutral pion production","first hadronic interaction","hadronic interaction models","muon problem","forward LHC physics"],"falsifier":"If, on real data with good statistics, the low-$N_\\mu$ tail slope changed with atmospheric depth or with the amount of later shower development while the first-interaction energy flow were held fixed, the claimed direct link would fail; a direct test would be to compare the $\\Lambda_{\\rm had}$ inferred from $\\Lambda_\\mu$ at approximately $10^{17}$ eV with the same quantity measured by forward LHC experiments at 13 TeV.","tokens_in":6020,"feed_emoji":"🌌","tokens_out":9229,"duration_ms":88593,"temperature":0.7,"pith_summary":"Ultra-high-energy cosmic-ray air showers hide their first collision: by the time a proton hits an air nucleus at center-of-mass energy around 100 TeV, the details of that interaction are buried under tens of generations of secondary particles. This paper argues that one fragment of the ground signal escapes that burial. Using simulations with three post-LHC hadronic models, it shows that the exponential low-value tail of the shower-to-shower distribution of muon number, $N_\\mu$, is set by the energy flowing into the hadronic sector in the very first interaction, not by the later cascade. The paper then connects that first-interaction energy split to the inclusive production of high-energy neutral pions, and shows that measuring the tail slope $\\Lambda_\\mu$ at a ground observatory would constrain the neutral-pion spectrum at $\\sqrt{s}\\sim100$ TeV, beyond accelerator reach. The practical payoff is a ground-based window into multiparticle production at energies no collider currently accesses.","feed_headline":"Muon tail slope exposes first-collision pion spectrum","feed_subtitle":"Simulated showers show the low-muon tail tracks the first hadronic collision, beyond accelerator reach","key_machinery":"The paper's central object is the exponential low tail of the ground muon-number distribution for proton-induced showers, quantified by its slope $\\Lambda_\\mu$. That slope is interpreted through a weighted first-interaction energy variable $\\alpha_1 = \\sum_i (E_{{\\rm had},i}/E_0)^\\beta$, with $\\beta = \\log(m)/\\log(m_{\\rm tot})$ from the Heitler-Matthews cascade model; this variable encodes how the first interaction splits energy between the hadronic and electromagnetic sectors. A calibration curve between $\\Lambda_\\mu$ and $\\Lambda_{\\rm had}$, obtained by reweighting simulated showers, is what lets a ground measurement stand in for the first-interaction energy flow. The neutral-pion spectrum enters through $E_{\\rm had} = 1 - E_{\\rm em}$, since $\\pi^0 \\to \\gamma\\gamma$ feeds the electromagnetic component and the highest-energy pions carry the information about fast leading particles.","core_discovery":"The central claim is that the shape of the shower-to-shower muon-number distribution, specifically its exponential tail at low $N_\\mu$, is controlled by the first hadronic interaction of the primary proton. The paper defines $\\Lambda_\\mu$ as the slope of this tail and shows, by reweighting simulated showers from a large ensemble, that $\\Lambda_\\mu$ is tied to $\\Lambda_{\\rm had}$, the slope of the distribution of the hadronic energy fraction $E_{\\rm had}/E_0$ deposited by that first interaction. Because the electromagnetic sector is fed almost entirely by neutral pions, fluctuations of $E_{\\rm had}/E_0$ are the same as fluctuations of the electromagnetic energy fraction, so $\\Lambda_{\\rm had}$ in turn reflects the inclusive high-energy tail of the neutral-pion spectrum. Simulation with QGSJET-II.04, EPOS-LHC, and SIBYLL 2.3c shows that modifying the neutral-pion inclusive cross-section at large $x_L$ moves $\\Lambda_\\mu$, and that a mixed composition of 25% protons, 50% helium, and 25% nitrogen still permits a clean measurement of the proton tail when the muon number is smeared by 20%.","pith_inferences":["A natural extension the authors do not spell out: the same tail analysis could be applied to the electromagnetic component, such as shower-maximum depth or ground electrons, to separate first-interaction energy flow from subsequent shower physics and cross-check the muon-based calibration.","If the model-dependent step between $\\Lambda_{\\rm had}$ and the neutral-pion spectrum cannot be sharpened, the practical reach of the method may be limited to detecting qualitative deviations from standard-model expectations, such as an anomalous leading-particle energy fraction, rather than a precise spectrum.","The technique effectively converts the atmosphere into a forward detector for $\\sqrt{s}\\sim100$ TeV proton-air collisions; a dedicated measurement of $\\Lambda_\\mu$ as a function of primary energy could map the onset of any new high-energy behavior.","An obvious testable extension is to apply the tail-slope method to existing high-statistics data sets and compare the derived $\\Lambda_{\\rm had}$ at overlapping energies with forward LHC measurements; any mismatch would pinpoint either the calibration or new physics."],"forward_implications":["A ground array that collects on the order of 3000 proton-like showers could distinguish among the tested hadronic interaction models using the low-$N_\\mu$ tail slope alone.","The same measurement at cosmic-ray energies around $10^{17}$ eV could be confronted with LHC forward measurements at $\\sqrt{s}=13$ TeV, offering a direct accelerator-to-cosmic-ray cross-check.","If the calibration holds, $\\Lambda_\\mu$ becomes a measurement of the fluctuation of the hadronic energy fraction in the first interaction at $\\sqrt{s}\\sim100$ TeV, where no collider data exist.","A suppression of neutral-pion production at large $x_L$ shifts the muon tail, so the tail slope provides a test of high-rapidity multiparticle production and of possible violations of longitudinal scaling at ultra-high energies.","The measurement remains feasible even under a pessimistic mixed composition with abundant helium, provided enough showers are recorded."],"supporting_citations":[{"why":"Earliest work by the same authors establishing that shower-to-shower muon-number fluctuations are driven by the first interaction's hadronic energy fraction; this paper extends that result to the tail slope.","marker":"[8]"},{"why":"Heitler-Matthews cascade model used to derive the power-law index $\\beta$ in the weighted first-interaction energy variable $\\alpha_1$.","marker":"[9]"},{"why":"QGSJET-II-04 simulations provide one of the three hadronic models used to test the universality of the $\\Lambda_\\mu$ to $\\Lambda_{\\rm had}$ calibration.","marker":"[5]"},{"why":"EPOS-LHC simulations provide the second model and are used for the mixed-composition precision study.","marker":"[6]"},{"why":"SIBYLL 2.3c simulations provide the third model and supply the nominal neutral-pion spectrum that is modified to test the tail's sensitivity.","marker":"[7]"},{"why":"The proton-air cross-section measurement from the $X_{\\max}$ tail supplies the template for extracting a proton tail slope in the presence of mixed composition.","marker":"[10]"},{"why":"Composition measurements in the relevant energy range justify the pessimistic 25/50/25 proton/helium/nitrogen scenario used to test feasibility.","marker":"[11]"}],"fun_headline_variants":["Muon tail metering yields first-collision pion spectrum","Low-muon tail reveals ultra-high-energy pion yield","Shower muons fingerprint proton-air collision pions","Tail slope ties ground muons to first-interaction pions","Muon statistics trace first-hit pion spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The link between the measured low-muon tail slope and the first-interaction hadronic-energy slope is universal across hadronic models and remains valid in real mixed-composition data; the paper demonstrates it only by reweighting simulated showers, not by analytic derivation.","fun_headline_variants_meta":{"raw":{"variants":["Muon tail metering yields first-collision pion spectrum","Low-muon tail reveals ultra-high-energy pion yield","Shower muons fingerprint proton-air collision pions","Tail slope ties ground muons to first-interaction pions","Muon statistics trace first-hit pion spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000614,"raw_usage":{"total_tokens":2858,"prompt_tokens":951,"completion_tokens":1907,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":1829}},"tokens_in":567,"tokens_out":1907,"duration_ms":12129,"temperature":1.0,"reasoning_tokens":1829,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:05:24.109429+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If, on real data with good statistics, the low-$N_\\mu$ tail slope changed with atmospheric depth or with the amount of later shower development while the first-interaction energy flow were held fixed, the claimed direct link would fail; a direct test would be to compare the $\\Lambda_{\\rm had}$ inferred from $\\Lambda_\\mu$ at approximately $10^{17}$ eV with the same quantity measured by forward LHC experiments at 13 TeV.","supporting_citations":[{"cited_title":"Probing the energy spectrum of hadrons in proton air interactions at ultrahigh energies through the fluctuations of the muon content of extensive air showers","cited_arxiv_id":"1803.05699","evidence_quote":"Earliest work by the same authors establishing that shower-to-shower muon-number fluctuations are driven by the first interaction's hadronic energy fraction; this paper extends that result to the tail slope."},{"cited_title":"Matthews, A Heitler model of extensive air showers, Astropart","cited_arxiv_id":null,"evidence_quote":"Heitler-Matthews cascade model used to derive the power-law index $\\beta$ in the weighted first-interaction energy variable $\\alpha_1$."}],"review_version":1}