{"id":"6d5f7ec6-5780-426f-9ca0-b0cbe6a01dfa","arxiv_id":"astro-ph/0403630","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"A full 3D Monte Carlo of atmospheric neutrino production agrees with the 1D calculation above 5 GeV except for azimuthal asymmetries caused by geomagnetic bending of muons in the atmosphere.","lead":"The Bartol atmospheric-neutrino group rebuilt its one-dimensional cascade code as a full three-dimensional Monte Carlo that tracks particles across the whole globe and bends them in the geomagnetic field. Running both versions with identical physics input, they show the cheap 1D calculation is adequate above about 5 GeV, with residual azimuthal differences from muon bending that persist to higher energy.","discovery_kind":"extension","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The >5 GeV 'essentially identical' claim is a null result stated qualitatively: zenith agreement is shown visually in one bin (6.3–10 GeV) but never given a residual bound, and the surviving muon-bending effect is assumed purely azimuthal although Kamioka's field is inclined.","rationale":"The reader's weakest_assumption — the 500 km flat-detector quadratic-fit correction (Eq. 2, Fig. 15) reported without a propagated uncertainty — is a real methodological soft spot, but it is mislocated relative to the headline claim. §3.4 states explicitly that the correction C is largest at 300–500 MeV (below 0.5% for r = 500 km) and 'at higher energies, C becomes negligible.' The central claim lives entirely above 5 GeV, where the detector-size concern does not bite. So I agree the correction lacks an error bar but disagree that it is the load-bearing threat to abstract:s4. The genuinely load-bearing gap is that the paper's central assertion is a null result expressed qualitatively. A methods paper whose stated purpose is to tell users when 1D suffices owes a number: the residual 3D−1D zenith difference above 5 GeV, with its statistical error. Fig. 3 shows overlap in a single 6.3–10 GeV bin, and the 4.0–6.3 GeV bin straddles the quoted threshold, so '5 GeV' is a stated boundary rather than a measured one. The companion issue is the azimuthal-only framing: the paper demonstrates (via NM, Fig. 8) that bending causes the azimuthal difference, but never shows that bending produces no zenith residual at high energy, even though Kamioka's inclined field makes a vertical bending component expected. Neither point overturns the result — the same-physics-input comparison is a genuine strength, the pseudo-1D/original-1D cross-check validates the weighting, and the toy-generator test against Lipari [9] (Fig. 1) is independent support. The acknowledged limitations (single hadronic model, primary spectrum discrepant with AMS/BESS below 50 GeV, no released code) are already flagged by the authors and by the reader and are not new findings. My concern is a request for one quantitative number that the existing runs can supply, so I would leave the verdict at ACCEPT/UNCHANGED rather than move it to CONDITIONAL.","tokens_in":19377,"tokens_out":4633,"duration_ms":76291,"concrete_test":"Re-analyze the existing 3D, pseudo-1D and NM runs (no new simulation needed): for every zenith bin with Eν > 5 GeV, compute the fractional difference (3D−1D)/1D and (3D−NM)/NM with propagated Monte-Carlo errors, and report whether each is consistent with zero at the 1% level. If any zenith bin above 5 GeV shows a >2% difference outside errors, the 'identical above 5 GeV' claim weakens and the threshold must be raised; if the 3D−NM zenith residuals vanish above 5 GeV while the azimuthal ones persist, the 'purely azimuthal bending' attribution is confirmed. Either outcome converts the qualitative null result into a falsifiable bound.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim (abstract:s4) is a null result — 1D and 3D agree above 5 GeV except in azimuth — and null results need a quantified sensitivity bound. The paper supplies figures, not numbers. §3 states 'At neutrino energies exceeding 5 GeV the two zenith angle distributions become identical' and the Conclusions repeat 'completely negligible for Eν > 5 GeV', but no fractional residual (3D−1D)/1D with its Monte-Carlo error is ever quoted for any zenith bin above 5 GeV. The only high-energy zenith evidence is Fig. 3 right panel, bin 6 (6.3–10 GeV), where the curves overlap within error bars; bin 5 (4.0–6.3 GeV) straddles the stated threshold, so the '5 GeV' boundary is asserted rather than demonstrated. The second, sharper issue: the paper attributes the entire high-energy residual to muon bending and frames it as azimuthal ('small effects in the azimuthal distributions'). But the geomagnetic field at Kamioka is inclined, so a muon bending perpendicular to B has a vertical component as well as a horizontal one; a purely azimuthal residual is an assumption, not a shown result. The NM control (Fig. 8) establishes that bending drives the azimuthal difference, but no figure decomposes the high-energy 3D−1D difference into zenith versus azimuth channels. Because the paper's practical value is licensing 1D codes for high-energy zenith analyses, the size of any zenith residual above 5 GeV is exactly what a user needs and exactly what is left unquantified. This is a quantification gap, not a refutation: the figures are consistent with agreement at the few-percent level, and the same-physics-input design correctly isolates the geometric/transport difference so hadronic-model uncertainty largely cancels in the 3D/1D ratio.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The paper describes the extension of the Bartol one-dimensional atmospheric-neutrino Monte Carlo (Agrawal et al., Barr et al.) to a full three-dimensional treatment in which primaries are injected over the whole globe, secondaries are tracked with transverse momentum, and charged particles are bent in the geomagnetic field during propagation in the atmosphere. The authors document the technique in detail: backtracking with Liouville's theorem for the primary cutoff, an adjustable-step field integrator using IGRF parameterizations, a flat circular sampling detector of radius 500 km with 'binlet' weighting to control the 1/cos(theta_z) divergence near the horizon, and five operating modes (3D, NM, pseudo-1D, 1D, original-1D) that allow the transition from 1D to 3D to be made in steps within one code. Results are given for Kamioka, Soudan and SNO. The central claim is that 1D and 3D zenith distributions are essentially identical above 5 GeV, with residual differences confined to azimuthal distributions and attributed to muon bending in the geomagnetic field.","tokens_in":19690,"tokens_out":5081,"duration_ms":91958,"significance":"This is a careful and useful methods paper. Its value lies less in any single number than in the systematic decomposition it provides: the pseudo-1D mode reproduces the original 1D code, which validates the weighting machinery that the 3D calculation depends on, and the NM mode isolates the role of in-atmosphere bending from the purely geometrical transverse-momentum effect. The demonstration that the low-energy horizontal enhancement is geometrical while the East-West asymmetry is geomagnetic, and that the latter persists to 10 GeV and differs between neutrinos and antineutrinos, is a genuinely useful result for oscillation analyses and for planning detectors with charge identification. The paper is also candid about its own approximations, including the outdated primary spectrum and the deferred comparison with other calculations. The practical conclusion, that 1D codes remain adequate for high-energy zenith analyses but not for azimuthal observables, is well motivated by the figures, though it would be more useful if quantified.","major_comments":[{"comment":"The headline claim of abstract:s4 is a null result, and null results need a stated sensitivity bound. §3 asserts that 'At neutrino energies exceeding 5 GeV the two zenith angle distributions become identical' and the Conclusions repeat 'completely negligible for Eν > 5 GeV', but no fractional residual (3D-1D)/1D with its Monte-Carlo error is quoted for any zenith bin above 5 GeV. The only high-energy zenith evidence is Fig. 3, right panel, bin 6 (6.3-10 GeV), where the histograms overlap within error bars; bin 5 (4.0-6.3 GeV) straddles the stated threshold, so the 5 GeV boundary is asserted rather than demonstrated. Since the paper's practical function is to license 1D codes for high-energy zenith analyses, the size of the residual above threshold is exactly what a user needs. I recommend a short table giving the 3D/1D ratio and its Monte-Carlo error per zenith bin for at least the 4-6.3, 6.3-10 and >10 GeV ranges, so that the 5 GeV boundary is derived from the numbers rather than from visual overlap.","section":"§3, Fig. 3 and Conclusions","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract attributes the surviving high-energy difference to 'small effects in the azimuthal distributions', but no figure decomposes the high-energy 3D-1D difference into zenith and azimuth channels. The NM control in Fig. 8 establishes that bending drives the azimuthal difference, which is convincing; it does not establish that the bending-induced residual is purely azimuthal. The geomagnetic field at Kamioka is inclined, so a muon bending perpendicular to B has a vertical component as well as a horizontal one, and a zenith residual from bending is plausible a priori. The claim that the residual is azimuthal is therefore an assumption rather than a shown result. A decomposition of the >5 GeV 3D-1D difference into its zenith and azimuth components, for example by plotting the zenith distribution of the 3D/NM ratio above 5 GeV alongside the azimuthal one, would settle this and would sharpen the guidance the paper gives to users.","section":"§3.2 and abstract:s4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The 3D flux is extracted from a flat circular sampling detector of radius 500 km, and the assertion that this averaging introduces less than 0.5% bias rests on a quadratic fit phi(x) = A0 + A1 x + A2 x^2 to the flux variation with latitude, evaluated through Eq. (2). Fig. 13 shows the flux itself varying by a factor of 3 across that region at 300-500 MeV, so the quoted correction is a small residual of a large gradient and depends on the quadratic form being adequate in every zenith bin. The fit uncertainty is not propagated into the quoted 0.5%, and the statement that at higher energies 'C becomes negligible' is not backed by a figure. This matters because the 'identical within statistical differences of about a per cent' statement for angle-averaged fluxes above 1 GeV in §3 is at the same level as the unquantified detector-averaging bias. I ask the authors to attach an uncertainty to C from the fit residuals, and to state explicitly whether the 0.5% figure is a bound or a central value.","section":"§3.4, Eq. (2), Fig. 15","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The expression '<pT>/Eπ ~ 300 MeV/Eπ ~ 0.1 Eν GeV radians' is dimensionally garbled; the final term appears to be intended as 0.1/Eν(GeV) radians. Please correct the units and check the numerical coefficient against the stated 300 MeV.","section":"Eq. (1)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The sentence '80 binlets are used in the range -1 < cos(theta_z) < 1, so the maximum variation in weights is slightly less than 1:80' is unclear. I take it to mean that the smallest bin-centre |cos(theta_z)| is 1/80, so the largest weight factor is 80; if so, please say that directly, since as written it reads as a statement about the ratio between adjacent bins, which is not bounded by 80 near the horizon.","section":"§2.2, binlet weighting","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The author of Ref. [20] is spelled 'Plyaskin' in the reference list but 'Playskin' in the text. Please make the spelling consistent.","section":"§2.1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The paper states that the surface of the Earth is treated as sea level everywhere but that detectors are assigned their actual altitude, and then lists SNO at an altitude of -1.7 km. Please clarify what 'altitude' means for an underground detector and confirm that the SNO value has no effect on the flux other than through the negligible baseline, so that a reader does not mistake it for a production-altitude input.","section":"§2.2, Earth model","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The paper surveys other 3D calculations in §2.1 but presents no numerical comparison with any of them, deferring this to a later publication. Given that the title advertises a 3D calculation of atmospheric neutrinos rather than only a 1D/3D comparison, at least one number, for example the angle-averaged flux against Honda et al. [16] or Tserkovnyak et al. [17] at a common site and epoch, would let a reader judge whether the absolute normalization is credible and not only whether it differs from the authors' own 1D result.","section":"§2.1 and §3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The only validation of the full 3D chain against an external benchmark is Fig. 1, which compares to analytic formulae from Ref. [9], a coauthor's paper. This is a legitimate and useful closure test, but it is a self-consistency check within the collaboration rather than independent validation. The text should say so plainly, and the deferred comparison with other groups becomes correspondingly more important.","section":"Fig. 1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The ~3% sub-GeV excess in the 3D angle-averaged flux is smaller than the spread among primary-spectrum measurements that the authors themselves discuss in §2.2. It would help to state this explicitly so that readers do not over-interpret the 3% figure as a correction that must be applied regardless of spectral uncertainties.","section":"§3 and Conclusions","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The 1D curves in these figures are drawn without error bars, while the 3D points carry them. Since several conclusions rest on the two being 'identical', please state the statistical uncertainty on the 1D samples, or note in the captions that it is negligible compared with the 3D errors.","section":"Figures 3, 6, 7, 9","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"There are typographical errors in the author affiliations ('Dipartmento di Fisica', 'Universit´ a di Roma') and the date line reads 'October 30, 2018' on a 2004 manuscript. These should be corrected in the final version.","section":"Title page","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision"},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"This is the Barr/Gaisser/Lipari/Robbins/Stanev 3D atmospheric neutrino paper, and its value is exactly what the title says: a controlled comparison, same physics input, of 3D versus 1D across five operating modes of one code. The result you'd cite is that above 5 GeV the zenith distributions agree while flavour-dependent East-West asymmetries from muon bending survive to 10 GeV. That is a useful, practical statement, and it is delivered with an unusually explicit approximation budget for 2004: detector size, binlet width, atmosphere model, altitude, all examined and bounded. The internal cross-checks are the strongest part. Pseudo-1D reproduces the original 1D code, which validates the weighting machinery, and the toy generator reproduces Lipari's analytic formulae in Fig. 1. That is real evidence, not decoration. The NM mode cleanly isolates bending as the driver of the azimuthal residual, and the attribution to muon bending rather than the primary cutoff is well argued in §3.2 and Fig. 9. Credit where due: this is honest, systematic work, and the authors flag their own input weaknesses (one hadronic model, a primary spectrum they note is discrepant with AMS/BESS below 50 GeV). On the stress-test concern: it half-lands. The >5 GeV claim is a null result and the paper does not quote a fractional (3D-1D)/1D residual with its Monte-Carlo error for any zenith bin above threshold; Fig. 3 bin 6 shows overlap within error bars, and bin 5 straddles the stated 5 GeV line, so the boundary is asserted more than demonstrated. That is a quantification gap, not a refutation, and the figures are consistent with few-percent agreement. The sharper point, that an inclined field at Kamioka could give a zenith component to the bending residual, is not addressed by any decomposition figure, but Fig. 3's convergence does bound it visually. So: a real but modest gap. The detector-size correction is the other soft spot: a quadratic-fit residual of a factor-of-3 gradient, reported below 0.5% with no propagated fit uncertainty. Minor, but worth a referee question. No code or tabulated fluxes are released, which limits reproducibility. Who should read it: anyone running or interpreting atmospheric neutrino flux calculations, and oscillation analysts deciding whether 1D is adequate. It deserves a serious referee, and I'd cite it.","headline":"A careful, well-cross-checked methods paper that establishes when 1D atmospheric neutrino calculations suffice; the main claim is sound but stated more sharply than the numbers behind it.","tokens_in":20531,"tokens_out":999,"duration_ms":17512},"pith_extraction":null,"created_at":"2026-11-13T08:47:39.244548+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"pipeline_version":"pith-pipeline@v1.5.0-archive-selfhost","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":null,"supporting_citations":[],"lane_accuracy":{"status":"not yet measured","lane":{"model":"selfhosted/deepseek-v4.1-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}