{"id":"4b3b7ffe-7e27-4d29-b41f-0efd606e0808","arxiv_id":"2507.10665","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using a line-of-sight transport formalism, the paper finds a 3.1 sigma tension between the KM3NeT ultra-high-energy event and IceCube non-observation under a diffuse power-law flux, with stronger tension for steady point sources.","lead":"A transport-equation framework connects the neutrino flux at Earth's surface to muon rates at neutrino telescopes, including muons from up-scatters and tau decays. Applied to the KM3NeT 220 PeV event and IceCube's 12-year null observation, it finds a 3.1 sigma tension for a diffuse power-law source, larger for steady point sources and smaller for transients.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported 3.1σ diffuse tension rests on a statistical conversion error: using the stated χ² calibration with B≈21 gives about 2σ, not 3.1σ.","rationale":"The paper develops a transparent first-principles transport-equation formalism and presents the numerical ingredients in detail, which is a genuine strength. The reader's conditional verdict focuses on simplified detector modeling and the 9.5-year vs 12-year exposure mismatch. Those are real but secondary: they affect the expected ratio Δ at the tens-of-percent level and can be clarified with modest effort. The most load-bearing issue is instead the statistical calibration that turns the Bayes factor B≈21 into the quoted 3.1σ. The formula in the Table 3 footnote uses Q(1; 2 ln B) = B⁻², whereas the survival function of a χ²₂ distribution at 2 ln B is Q(1; ln B) = B⁻¹. This is a factor-of-two error in the exponent of the p-value, changing the significance from about 3.1σ to about 2.0σ. Since the abstract and Table 1 present 3.1σ as the central result, this is an internal inconsistency that must be settled before the paper's main claim can be accepted. The verdict should remain conditional, but the condition should now include recomputing the significance with the correct chi-square calibration or with a direct frequentist approach. If the recalculation confirms roughly 2σ, the paper's claims of tension should be softened accordingly.","tokens_in":45392,"tokens_out":23069,"duration_ms":269735,"concrete_test":"Recompute the diffuse power-law significance from the Bayes factor B ≈ 21 using the correct χ²₂ survival function, p = exp(−ln B) = 1/B, instead of exp(−2 ln B). If the quoted significance drops from 3.1σ to roughly 2.0σ, the headline tension is overstated and the text and tables must be corrected. As a cross-check, compute a frequentist p-value directly from the Poisson likelihoods: under the IceCube-fit diffuse flux, use the ratio Δ ≈ 70 and the observed zero IceCube events and one KM3NeT event to evaluate the joint probability, and compare with Ref. [20].","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative claim is the 3.1σ tension quoted in the Abstract and Table 1. In Section 5.3 and the Table 3 footnote, the Bayes factor is converted to a significance via Δσ = √2 erfc⁻¹(Q(#dof/2; 2 ln B)). For #dof = 2, Q(1; 2 ln B) = exp(−2 ln B) = B⁻². But the survival function of a chi-square with two degrees of freedom at x = 2 ln B is Q(1; ln B) = exp(−ln B) = B⁻¹. The paper therefore uses B⁻² as the p-value instead of B⁻¹. With the reported B ≈ 21, the stated formula gives 3.1σ, while the correct χ²₂ calibration gives p = 1/21 ≈ 0.048, corresponding to about 2.0σ. This directly affects the headline claim of a '3.1σ tension' and the conclusion of a 'stark contrast' between KM3NeT and IceCube. Moreover, the conversion of a Bayes factor to a chi-square is itself a heuristic; the significance should be verified with a direct frequentist calculation or a posterior probability statement. This is an internal, checkable error in the statistical analysis, not a modeling or data-handling ambiguity.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a first-principles transport-equation formalism, in the line-of-sight approximation, to compute the muon event rate at neutrino telescopes from the neutrino flux at the Earth's surface. It includes both muon production from νμ charged-current scattering and from ντ→τ→μ decay chains, and it derives master formulas, Eqs. (4.7) and (4.13), that depend on the surface neutrino flux and Standard Model/environmental inputs. The formalism is applied to the KM3-230213A event in comparison with the IceCube non-observation under several flux hypotheses: diffuse power-law, energy-localized diffuse, point source, transient point source, τ-neutrino point source, and a BSM two-state scenario. The headline result is a claimed 3.1σ tension between KM3NeT and IceCube for a diffuse power-law flux, with a Bayes factor of about 21, stronger tension for point sources, and reduced tension for transient sources.","tokens_in":45661,"tokens_out":7051,"duration_ms":81371,"significance":"If the quantitative claims survive scrutiny, the paper provides a transparent and physically motivated alternative to full detector Monte Carlo for UHE neutrino telescopes. The explicit treatment of tau propagation, Earth density profiles, and the derivation of effective-area-like quantities from first principles are genuine strengths, and the master formulas are stated in a form that can be reused and checked. The paper is also honest about its main approximations, notably the unit-efficiency spherical detector and the neglect of tau energy loss. However, the central quantitative significance is weakened by a statistical conversion error and by an exposure inconsistency between the fitted and predicted IceCube data, so the headline conclusion cannot be accepted as it stands.","major_comments":[{"comment":"The conversion of the Bayes factor into a Gaussian significance is incorrect. The quoted formula Δσ = √2 erfc⁻¹(Q(#dof/2; 2 ln B)) with #dof = 2 evaluates Q(1; 2 ln B) = exp(−2 ln B) = B⁻². The survival probability of a chi-square distribution with two degrees of freedom at χ² = 2 ln B is instead Q(1; ln B) = exp(−ln B) = B⁻¹. For the reported B ≈ 21, the correct two-sided significance is about 2.0σ, not 3.1σ. Because the abstract and Table 1 present the 3.1σ value as the central result, the significance must be recalibrated, and the related statements ('stark contrast', and the σ values for point and transient sources) revised accordingly.","section":"§5.3 and Table 3 footnote"},{"comment":"The IceCube likelihood in Section 5.3 is fit to the 9.5-year dataset of Ref. [8], using 'the last 31 bins in the right panel of Fig. 1 of [8]', while Section 5.1 and Table 2 state a 12-year IceCube exposure at the time of KM3-230213A, and Eq. (5.5) uses the combination [T A]_IC. If the flux normalization fitted to 9.5 years is used together with the 12-year exposure to predict the UHE-bin counts without an explicit rescaling, the expected numbers entering the Bayes factor are biased by a factor of 12/9.5 ≈ 1.26. The paper should specify which exposure is used for each step and make the rescaling explicit and consistent.","section":"§5.1, §5.3, Table 2"},{"comment":"The single high-energy bin [10 PeV, 400 PeV] is introduced specifically to contain the KM3-230213A event after that event was observed. The reported tension is therefore conditional on a bin choice made with knowledge of the data. If the bin has a physics motivation independent of the event, that motivation should be stated; otherwise the significance should be accompanied by a trials-factor correction or be described as exploratory. This issue compounds the calibration problem identified above.","section":"§5.2"}],"minor_comments":[{"comment":"The contour labels 'egg' and 'vegan bacon' are informal and unexplained; they should be replaced with standard posterior-density labels or a legend description.","section":"Fig. 13 caption"},{"comment":"The caption defines Δσ but not the meanings of B and ̃B in a self-contained way; the definitions in the text should be repeated briefly in the caption.","section":"Table 3 caption"},{"comment":"The statement that more detailed efficiency choices 'do not change the result of Section 5' would be easier to trust if a quantitative robustness bound or a short scan were shown.","section":"§4.1"},{"comment":"The approximation bτ = 0 is stated, but Section 5.5 later explains that tau energy loss can affect the conclusions; a quantitative estimate of the size of the bτ correction on Eq. (4.13) would be useful.","section":"§4.4"},{"comment":"There is a typo in the heading 'statystical treatment', which should read 'statistical treatment'.","section":"§5 heading"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the 3.1σ conversion is correct and important: the paper's headline significance is overstated by about 1σ. The exposure mismatch between the 9.5-year fit and the 12-year prediction is also present and should be corrected. These are fixable within the scope of the manuscript, and the underlying transport-equation formalism and Bayes-factor comparison are a useful contribution, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper is worth reading for the transport machinery, but the central claim needs a fix. The line-of-sight Boltzmann solution and the analytic event-rate formulas (Eqs. 4.7–4.14) are genuinely useful, and the day-averaged point-source transparency treatment is new relative to the existing literature. The application to KM3-230213A is clear, and the BSM section is properly labeled as speculation.\n\nThe problem is statistics. The abstract and Table 1 report a 3.1σ tension from a Bayes factor of about 21. The conversion formula in the Table 3 footnote is wrong as written: they define Δσ = √2 erfc⁻¹(Q(#dof/2; 2 ln B)). For #dof=2 that is Q(1; 2 ln B) = exp(−2 ln B) = B⁻², giving p ≈ 1/441 ≈ 0.0023, which corresponds to about 3.1σ. But the asymptotic chi-squared survival with 2 dof evaluated at the statistic 2 ln B is Q(1; ln B) = exp(−ln B) = 1/B ≈ 0.048, which is about 2.0σ. So the headline significance is overstated by roughly 1σ. This is an internal, checkable arithmetic error, not a modeling ambiguity.\n\nBeyond that, the detector modeling is simplified: unit-efficiency sphere, a single post hoc UHE bin [10, 400] PeV, and the IceCube fit uses 9.5-year data while predictions use 12-year exposure without explicit rescaling. Those are minor-to-moderate issues; they could shift the event ratio by O(1) factors but would likely not erase a ~2σ tension. The transport formalism itself is internally consistent, the inputs (PREM, MadGraph cross-sections, PDF fits) are standard, and the paper is transparent about its approximations. The point-source result — stronger tension than the diffuse case, contrary to Ref. [20] — is interesting and worth checking.\n\nWho should read this: anyone doing UHE neutrino phenomenology or planning next-generation detectors; the analytic formulas could be reused. If I referee, I would ask for a corrected significance, a direct frequentist check or posterior probability statement, and a clarification of the exposure bookkeeping. I would not desk-reject; the formalism deserves a serious review. My verdict would be conditional on those fixes.","headline":"Solid transport formalism, but the headline 3.1σ tension is inflated by a Bayes-factor-to-significance conversion error; corrected calibration gives about 2σ.","tokens_in":46213,"tokens_out":4390,"would_cite":true,"duration_ms":46476,"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 first-principles line-of-sight transport calculation shows the KM3NeT 220 PeV muon is in 3.1 sigma tension with IceCube under a diffuse power-law neutrino flux, with only a transient source capable of lowering the tension to 1.6 sigma.","keywords":["ultra-high-energy neutrinos","neutrino telescopes","transport equation","line-of-sight approximation","KM3-230213A","IceCube","KM3NeT","Bayes factors"],"falsifier":"Recompute the expected IceCube count in the 10 to 400 PeV bin using the detector's simulated effective area and the full 12-year livetime rather than a geometric disk; if the predicted ratio to KM3NeT drops from about 70 toward about 10, the 3.1 sigma tension falls below 3 sigma. More directly, a single IceCube muon above 10 PeV from the KM3-230213A sky region would falsify the paper's central claim.","tokens_in":45154,"feed_emoji":"🌊","tokens_out":9517,"duration_ms":104485,"temperature":0.7,"pith_summary":"The paper derives the expected muon rate at a neutrino telescope from the transport equation in the line-of-sight approximation, so that the only free parameters are the normalization, spectral index, and directionality of the neutrino flux at Earth's surface. It applies this machinery to the KM3-230213A event, a 120+110-60 PeV muon track seen by KM3NeT at elevation 0.54 degrees, and compares the predicted rates with IceCube's non-observation of similar events. The central quantitative result is a 3.1 sigma tension between the two experiments under a diffuse power-law flux, corresponding to a Bayes factor of about 21 and an expected event ratio of about 70 in the 10 to 400 PeV bin. A steady point source sharpens the tension to 3.8 sigma, while a transient point source would relax it to 1.6 sigma. The formalism also clarifies how detector volume, Earth density, and lepton energy loss enter the event rate, and it provides a roadmap for including beyond-Standard-Model muon sources.","feed_headline":"One 220 PeV neutrino clashes with IceCube at 3.1 sigma","feed_subtitle":"Diffuse power-law sources fail a 3.1-sigma test; transient sources drop the tension to 1.6 sigma.","key_machinery":"The load-bearing object is the muon phase-space distribution $f_\\mu$ computed by the line-of-sight solution of the Boltzmann transport equation: one integrates along the straight neutrino trajectory from the crust to the detector, using the Earth's density profile, the deep-inelastic-scattering neutrino cross section, and the charged-lepton energy-loss coefficients as inputs. The master outputs are Eq. (4.7) for muon-neutrino events and Eq. (4.13) for tau-neutrino events, each an integral over the line of sight of a transport function times the surface flux. These formulas expose the effective detector volume as $A_{\\rm disk}(R_{\\rm det} + 1/[b_\\mu(2+\\gamma-\\lambda)])$, so the inverse radiative stopping length $1/b_\\mu$, not the geometric volume, sets the collecting power for through-going muons.","core_discovery":"The central claim is that the differential muon event rate at a neutrino telescope can be written, without Monte Carlo event generation, as an integral along the observer's line of sight through Earth of transport functions acting on the surface neutrino flux, with muons produced either by neutrino up-scattering near the detector or by tau decays at a distance. For a diffuse power-law flux the resulting master formulas, Eq. (4.7) for $\\nu_\\mu$ and Eq. (4.13) for $\\nu_\\tau$, imply an IceCube-to-KM3NeT expected-event ratio near 70 in the ultra-high-energy bin, an effective-volume enhancement $1/b_\\mu$ beyond the geometric detector volume, and a 3.1 $\\sigma$ tension between the KM3NeT event and the IceCube null. The same calculation gives a stronger tension for steady point sources and for energy-localized diffuse sources, and a milder 1.6 $\\sigma$ tension if the source is transient.","pith_inferences":["Editorially, the same line-of-sight integrals can be applied to other large neutrino telescopes and to future high-energy upgrades; the framework's practical payoff is that nuisance parameters such as the cross-section slope, Earth density, and energy-loss coefficients enter as calculable theory inputs rather than as Monte Carlo systematics.","The 3.1 sigma number is tied to treating IceCube's exposure as a single zero-event bin; a revised fit using a longer livetime or a different high-energy binning could move the significance, so the exact sigma should be read as a model-dependent estimate.","A discriminating test the paper leaves implicit is that if more ultra-high-energy events accumulate, the energy distribution distinguishes transient from steady sources: a transient source populates a narrow energy window, whereas a power-law diffuse source predicts a particular falling spectrum.","If the tension persists, the paper's BSM discussion points to a concrete signature: a long-lived particle with decay length between the two detectors' chord lengths would preferentially produce muons at KM3NeT, visible as a high-energy track without a strong gamma-ray counterpart."],"forward_implications":["For a diffuse power-law flux, the expected ratio of ultra-high-energy events between IceCube and KM3NeT is about 70 and only mildly dependent on the spectral index; with IceCube's fit as reference this translates into a Bayes factor of 21 and a 3.1 sigma tension.","A steady point source at the KM3-230213A position increases the ratio to about 140 and the tension to 3.8 sigma, so point sources do not relieve the discrepancy.","An energy-localized diffuse source behaves like a diffuse power law, with ratio about 67 and tension 2.4 sigma, while a transient point source lowers the ratio to about 11 and the tension to 1.6 sigma, the only considered Standard Model scenario that substantially reduces it.","Tau-neutrino point sources do not favor KM3NeT once tau energy loss is accounted for; even a pure $\\nu_\\tau$ source keeps the event ratio above about 70 at the relevant energies.","Because the effective volume grows as $1/b_\\mu$, differences in muon energy loss in ice versus water materially change the expected-event ratio; ignoring this enhancement would push the ratio closer to 100."],"supporting_citations":[{"why":"Reports the KM3-230213A event with its reconstructed energy and direction, the datum being tested against IceCube.","marker":"[15]"},{"why":"Supplies the 9.5-year IceCube muon-neutrino flux measurement used to fit the diffuse power-law reference.","marker":"[8]"},{"why":"Provides updated IceCube event directions and the exposure information used for the ultra-high-energy null bin.","marker":"[9]"},{"why":"Is the earlier tension analysis whose 3.1 sigma result this paper reproduces and refines with its own statistical criterion.","marker":"[20]"},{"why":"Underlies the treatment of neutrino attenuation and tau propagation along chords through Earth used in the line-of-sight transport functions.","marker":"[23]"},{"why":"Defines the reference Earth density profile that sets the nucleon number density along each line of sight.","marker":"[24]"},{"why":"Constrains the small-x PDF exponent that controls the high-energy neutrino cross-section uncertainty band.","marker":"[27]"},{"why":"Gives the range of small-x power-law behaviors used to estimate the cross-section extrapolation uncertainty.","marker":"[30]"},{"why":"Supplies the charged-lepton energy-loss coefficients for bremsstrahlung, pair production, and deep inelastic scattering used in $b_\\mu$ and $b_\\tau$.","marker":"[40]"}],"fun_headline_variants":["3.1 sigma tension between KM3NeT event and IceCube null","Transient sources may explain neutrino tension at 1.6 sigma","Muon rate from neutrino flux via transport equation","KM3NeT 220 PeV event clashes with IceCube at 3.1 sigma","Line-of-sight transport links muon rate to neutrino flux"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on modeling IceCube's ultra-high-energy exposure as a single Poisson bin from 10 to 400 PeV with zero observed events, and on treating both detectors as unit-efficiency spheres with geometric transverse area.","fun_headline_variants_meta":{"raw":{"variants":["3.1 sigma tension between KM3NeT event and IceCube null","Transient sources may explain neutrino tension at 1.6 sigma","Muon rate from neutrino flux via transport equation","KM3NeT 220 PeV event clashes with IceCube at 3.1 sigma","Line-of-sight transport links muon rate to neutrino flux"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000898,"raw_usage":{"total_tokens":3913,"prompt_tokens":1038,"completion_tokens":2875,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":2781}},"tokens_in":654,"tokens_out":2875,"duration_ms":20974,"temperature":1.0,"reasoning_tokens":2781,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:29:35.791982+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the expected IceCube count in the 10 to 400 PeV bin using the detector's simulated effective area and the full 12-year livetime rather than a geometric disk; if the predicted ratio to KM3NeT drops from about 70 toward about 10, the 3.1 sigma tension falls below 3 sigma. More directly, a single IceCube muon above 10 PeV from the KM3-230213A sky region would falsify the paper's central claim.","supporting_citations":[{"cited_title":"Aiello et al.,Observation of an ultra-high-energy cosmic neutrino with KM3NeT, Nature 638 (2025) 376–382","cited_arxiv_id":null,"evidence_quote":"Reports the KM3-230213A event with its reconstructed energy and direction, the datum being tested against IceCube."},{"cited_title":"Observing EeV neutrinos through the Earth: GZK and the anomalous ANITA events","cited_arxiv_id":"1909.10487","evidence_quote":"Underlies the treatment of neutrino attenuation and tau propagation along chords through Earth used in the line-of-sight transport functions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the reference Earth density profile that sets the nucleon number density along each line of sight."},{"cited_title":"Koehne, K","cited_arxiv_id":null,"evidence_quote":"Supplies the charged-lepton energy-loss coefficients for bremsstrahlung, pair production, and deep inelastic scattering used in $b_\\mu$ and $b_\\tau$."}],"review_version":1}