{"id":"1842f66b-f359-42b8-949a-07f16f487c4b","arxiv_id":"2411.19022","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"An analytic effective theory of the sterile-active MSW resonance in the (3+1) model reveals that only three channels, P(νe→νe), P(ν̄e→ν̄e), and P(ν̄μ→ν̄τ), carry unsuppressed on-peak resonance effects.","lead":"For neutrinos traveling through dense matter, a proposed extra sterile neutrino would create a resonance near 10 TeV. This paper builds a simplified theory of that resonance and maps out, channel by channel, which flavor transitions show the strongest effect and which events a neutrino telescope should count.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The anti-neutrino local-resonance approximation is the load-bearing premise: the unsuppressed P(ν̄μ→ν̄τ) central to the three-path cascade claim is derived from sequential diagonalization of overlapping resonances and is not numerically checked, so the central hierarchy remains conditional.","rationale":"Reading the paper in good faith, its aim is a qualitative global picture of the SA resonance, and the effective-theory construction is transparent and instructive. The S-matrix texture-zero structure and the resulting channel-by-channel suppression hierarchy are genuinely useful organizing results. The SA limit itself is well motivated because Δm²31/a is tiny, and the numerical spot-checks that are shown give roughly 1% agreement in the tested channels. However, the strongest claim includes a specific anti-neutrino appearance channel, P(ν̄μ→ν̄τ), as one of the three dominant cascade paths. That formula is obtained only under the local-resonance approximation, which the authors explicitly flag as potentially inaccurate in a wide fraction of parameter space. Crucially, this channel is absent from the numerical validation plots, and no full (3+1) comparison is provided for it. This is not a dispute with external consensus; it is an internal validation gap at the point where the central claim is most novel. The proposed check, a full numerical comparison of P(ν̄μ→ν̄τ) and P(ν̄μ→ν̄μ) in the relevant and overlapping parameter regions, would settle whether the concern lands. Because the reader's conditional verdict already identifies exactly this weakest assumption and recommends the same kind of test, my stress-test does not change the verdict: it confirms that conditional is the right call until the check is run.","tokens_in":46936,"tokens_out":4622,"duration_ms":61120,"concrete_test":"Using the same full (3+1) numerical solver as in Sec. 4.5, compute P(ν̄μ→ν̄τ) and P(ν̄μ→ν̄μ) over a baseline/energy grid spanning the SA region, for example L = 100-10000 km, E = 1-30 TeV, Δm²41 = 1 eV², ρ = 5.5 g/cm³, and θ14 = θ24 = θ34 = 10°, and compare with Eqs. (9.1) and (9.2). Require agreement within the 1-2% accuracy claimed for the effective theory. In addition, run the same comparison in the overlapping-resonance regime, e.g., scan Δm²41 over 0.5-10 eV² and θ24, θ34 over 5-15°, to see whether deviations grow where the 2-4 and 3-4 resonances overlap. If P(ν̄μ→ν̄τ) deviates by more than about 10% on peak, the unsuppressed-appearance claim and the three-path hierarchy do not survive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that cascade events come dominantly from three paths rests on the anti-neutrino-channel result that P(ν̄μ→ν̄τ) is unsuppressed on the SA resonance (Eq. 9.2), alongside P(ν̄μ→ν̄μ) (Eq. 9.1). These formulas are obtained by diagonalizing the 2-4 and 3-4 level crossings one by one (Sec. 8.1), the 'local resonance approximation.' The authors themselves warn that this approximation 'may not hold in a good accuracy' because the two resonances can overlap in a substantial fraction of parameter space (Secs. 8.1, 11.2). Unlike the νSM case, the relevant crossings are not widely separated, and no exact or non-sequential treatment is provided. The numerical validation in Sec. 4.5 and Fig. 2 tests only νμ→νμ, νμ→νe, νe→νe and the corresponding anti-neutrino disappearance/ν̄μ→ν̄e channels; it does not include P(ν̄μ→ν̄τ), the very appearance channel that is the novel basis for the cascade recommendation. Since the claimed hierarchy (Tables 1 and 2) and the three-component fit (Sec. 9.3) inherit this approximation, an unchecked failure there would change the phenomenological conclusion even if the effective theory itself is accurate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the MSW resonance induced by an eV-scale sterile neutrino in the (3+1) model, the \"sterile-active (SA) resonance,\" at E ~ 1-10 TeV in Earth matter. Starting from the (3+1) Hamiltonian in matter, the authors take the SA limit in which the large matter potential freezes the νSM oscillations; the resulting effective theory depends only on the sterile mixing angles θ14, θ24, θ34 and Δm²41. They then develop a perturbative treatment of this effective theory, using texture zeros of the S matrix and smallness of sj4 to identify a flavor/event-type hierarchy. The principal phenomenological claims are: (i) in the neutrino channel, the resonance effect is concentrated in P(νe→νe) with no sj4 suppression on resonance; (ii) in the anti-neutrino channel, P(ν̄μ→ν̄μ) and P(ν̄μ→ν̄τ) are unsuppressed on resonance; (iii) cascade events in neutrino telescopes should therefore be sought through three paths (νe→νe, ν̄e→ν̄e, ν̄μ→ν̄τ), motivating a three-component fit. A numerical comparison with the (3+1) code is presented for six channels at a benchmark point.","tokens_in":47265,"tokens_out":4769,"duration_ms":82091,"significance":"If correct, the paper provides a useful global map of the SA resonance and makes a concrete, falsifiable recommendation for IceCube/KM3NeT searches: neutrino-channel tracks are not promising, while νe disappearance and ν̄μ→ν̄τ appearance are the key cascade signatures. Strengths include the explicit derivation of the SA limit with an estimate of sub-asymptotic corrections, the unified treatment of neutrino and anti-neutrino channels, and the self-contained analytic S-matrix computation in the appendices. The effective theory itself is shown to agree with the (3+1) numerics at about the 1% level for the tested benchmark. However, the central phenomenological recommendation rests on the anti-neutrino appearance probability P(ν̄μ→ν̄τ), which is derived under the local resonance approximation and is not included in the numerical validation. That gap, together with the absence of a universal expansion parameter, makes the hierarchy in Tables 1 and 2 conditional.","major_comments":[{"comment":"The unsuppressed on-peak formula for P(ν̄μ→ν̄τ) in Eq. (9.2) is a central element of the \"three origins of the cascade events\" recommendation in Sec. 9.3, but it is derived under the local resonance approximation in which the 2-4 and 3-4 level crossings are diagonalized sequentially. The authors state explicitly in Secs. 8.1 and 11.2 that this approximation \"may not hold in a good accuracy\" because the two resonances can overlap over a substantial fraction of parameter space. The numerical validation in Sec. 4.5 and Fig. 2 covers P(νμ→νμ), P(νμ→νe), P(νe→νe) and the corresponding anti-neutrino channels, but it does not include P(ν̄μ→ν̄τ) or any anti-neutrino appearance channel. Since a failure of the sequential diagonalization would directly affect the claimed cascade hierarchy in Tables 1 and 2, I ask the authors to add a numerical comparison of the analytic P(ν̄μ→ν̄τ) against the (3+1) code over the resonance region, and to quantify the parameter-space region where the local resonance approximation is reliable.","section":"Sec. 8.1, Eq. (9.2), Sec. 11.2"},{"comment":"The numerical test of the effective theory is performed for a single parameter point (Δm²41 = 1 eV², θ14 = θ24 = θ34 = 10°, ρ = 5.5 g/cm³) and compares the effective-theory Hamiltonian with the full (3+1) code. The analytic perturbative expressions developed in Secs. 6-9, including the key probabilities in Eqs. (7.1), (7.2), (9.1), and (9.2), are not directly benchmarked against either the effective-theory numerics or the (3+1) code. The paper's qualitative conclusions are based on the perturbative hierarchy, so I request either a direct comparison of the analytic probability formulas with the numerical effective theory, or an explicit statement that the quoted 1% accuracy applies only to the effective theory itself and not to the subsequent perturbative approximations.","section":"Sec. 4.5 and Fig. 2"},{"comment":"As the authors acknowledge in Sec. 6.2, the perturbative framework lacks a universal expansion parameter, and the smallness of the correction terms is not guaranteed by the Hamiltonian decomposition alone. The channel hierarchy in Tables 1 and 2 relies on the perturbative order counting in powers of sj4, with sj4 ~ 0.1 as a rough benchmark. Because the next-to-leading order terms are not parametrically suppressed by a small expansion parameter, the paper would benefit from an explicit estimate or numerical demonstration that the retained first- and second-order terms dominate the omitted higher orders for the benchmark values used in the phenomenological discussion.","section":"Sec. 6.2 and Sec. 11.2"}],"minor_comments":[{"comment":"The phrase \"the cascade events dominantly comes from\" should be \"the cascade events dominantly come from\".","section":"Abstract and Sec. 9.3"},{"comment":"The notation \"ec14\" for the matter-affected cosine can be misread as a product \"e c14\"; consider introducing an explicit symbol such as c̃14 or a parenthetical definition at first use.","section":"Sec. 2.1 and throughout"},{"comment":"The labels \"NMO IMO\" appear below the two panels, but the individual panels are not clearly identified; labeling each panel directly would improve readability.","section":"Fig. 1"},{"comment":"The definition of F in Eq. (4.13) refers to P_eff-th, but the caption of Fig. 2 and the text in Sec. 11 describe the comparison slightly differently (effective theory versus analytic formulas); please make the object being compared unambiguous.","section":"Sec. 4.5 and Fig. 2"},{"comment":"Reference [54] cites a conference talk by name only; consider replacing it with a citable proceeding or archival entry, or remove it if not essential.","section":"References"},{"comment":"The sentence \"this is the situation we never met in the νSM\" would read more naturally as \"this is a situation never encountered in the νSM\".","section":"Sec. 8.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of JHEP, and the authors' own warnings in Secs. 8.1 and 11.2 about the local resonance approximation are the key issue. I would not recommend rejection: the effective theory construction is valuable and partially validated. The central concern is that the anti-neutrino appearance prediction P(ν̄μ→ν̄τ), which underpins the cascade recommendation, is neither numerically checked nor protected by a controlled expansion. Closing that gap with a dedicated numerical comparison should be feasible within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a real analytic advance, and it is honest about its biggest weakness. The effective theory construction is clean, and the texture-zero S-matrix plus the s_j4 suppression tables give a genuinely new bird-eye view of the SA resonance. But the headliner phenomenological claim—that cascades come dominantly from three paths, with P(ν̄μ→ν̄τ) unsuppressed—rides on a local-resonance approximation that the authors themselves say may not hold when the 2-4 and 3-4 crossings overlap. They never test P(ν̄μ→ν̄τ) numerically, and the only numerical check is at a single parameter point. So the hierarchy is plausible, not proven.\n\nWhat I like: simultaneous treatment of all channels in one framework; the S-matrix row/column zero structure is a good organizing principle; the suppression tables are useful and clearly explained; the T-invariance discussion is careful. The paper also does not oversell: it flags the lack of a universal expansion parameter, the small-angle assumption, and the local-resonance approximation explicitly. The derivation of the effective theory via the SA limit is transparent, and the comparison to the (3+1) code, while limited, shows 1% level agreement for the channels tested.\n\nSoft spots: (1) the anti-neutrino appearance channel central to the three-path claim is never checked against a full (3+1) numerical solution; (2) the sequential diagonalization is an approximation with no quantitative estimate of its failure in the overlapping regime; (3) the perturbative expansion has no small parameter except the assumed small sterile angles, so the suppression tables are heuristic, not a controlled series. These are all acknowledged, but the first two are load-bearing for the phenomenological punchline.\n\nWho should read it: people constructing search strategies for IceCube/KM3Net, and anyone working on analytic neutrino oscillation theory. The framework will save them time even if the specific numerical results need re-checking.\n\nRecommendation: yes, take it to peer review. The novelty is real, the analysis is careful, and the main gap—a numerical test of P(ν̄μ→ν̄τ) and the overlapping-resonance regime—is something a referee can reasonably require. I would not desk-reject it, but I would send it back for that check.","headline":"A genuinely new analytic framework for the SA resonance, honest about its main limitation: the cascade recommendation rests on an untested sequential-diagonalization approximation.","tokens_in":47830,"tokens_out":2485,"would_cite":true,"duration_ms":24800,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.60.Pq"],"model":"deepseek-v4-flash","headline":"The paper establishes that near the sterile-active resonance the (3+1) model reduces to an effective theory with only the sterile mixing angles and $\\Delta m^2_{41}$, with the resonance effect concentrated in three cascade paths.","keywords":["sterile-active resonance","(3+1) model","MSW resonance","effective theory","neutrino telescopes","cascade events","flavor oscillation probabilities","matter effects"],"falsifier":"Numerically diagonalize the exact (3+1) Hamiltonian in the anti-neutrino channel over the parameter region where the 2-4 and 3-4 resonances overlap, and compare the exact $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\mu)$ and $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\tau)$ with the analytic local-resonance expressions; if the discrepancy exceeds the few-percent level or the on-peak unsuppression disappears, the central hierarchy claim fails.","tokens_in":46719,"feed_emoji":"🔭","tokens_out":7801,"duration_ms":65493,"temperature":0.7,"pith_summary":"The paper claims that in the (3+1) model with an eV-scale sterile neutrino, the sterile-active resonance at energies around 1–10 TeV in the Earth can be captured by an effective theory that depends only on the sterile mixing angles and $\\Delta m^2_{41}$, because the large matter potential freezes the standard three-flavor oscillations. Working in that effective theory, the authors derive a perturbative expansion whose “texture zeros” in the flavor-space $S$ matrix organize the resonance effect into a clear hierarchy: at on-peak, the resonance is unsuppressed in $P(\\nu_e\\to\\nu_e)$, $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\mu)$, and $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\tau)$, while muon-neutrino track channels in the neutrino channel are suppressed by the small sterile mixing angles. The paper concludes that cascade events at neutrino telescopes would dominantly come from three paths with different resonance signatures, and proposes a three-component fit to separate them. If this picture is right, it tells the experimental program precisely which channels carry the sterile-resonance signal and which are nearly blind to it.","feed_headline":"Sterile neutrino resonance points to three cascade paths","feed_subtitle":"At TeV energies the effect is strong in electron survival and anti-muon to anti-tau, weak in muon tracks.","key_machinery":"The machinery is a limiting procedure that freezes the $\\nu$SM oscillations and a perturbative diagonalization of the resulting effective theory. In the SA region the matter potential makes the matter-affected angles $\\tilde\\theta_{13}$ and $\\tilde\\theta_{12}$ collapse to $\\pi/2$ (neutrino channel) or $0$ (anti-neutrino channel), turning the active 3-flavor rotations into discrete permutations and leaving an effective Hamiltonian built only from the sterile mixing angles and $\\Delta m^2_{41}$. The argument then runs on the “texture zeros” of the flavor-basis $S$ matrix — entries that vanish at zeroth or first order and dictate how many powers of $\\sin\\theta_{j4}$ multiply each probability — together with a smallness assumption $\\sin\\theta_{j4}\\sim 0.1$. In the anti-neutrino channel the two resonances (2-4 then 3-4) are diagonalized one by one, the “local resonance approximation,” which is the step that carries the unsuppressed on-peak results.","core_discovery":"On the paper's own terms, the central discovery is a global qualitative map of where the sterile-active resonance acts. After taking the SA limit, the effective Hamiltonian reduces to a $\\mathrm{diag}(0,0,0,\\Delta m^2_{41})$ mass structure plus the matter potentials, with the state space reshuffled: in the neutrino channel the 1-4 rotation acts on the physical 3-4 crossing, while in the anti-neutrino channel no reshuffling occurs. Diagonalizing with the matter-affected resonance angles $\\tilde\\theta_{14}$ (neutrino channel) and then $\\tilde\\theta_{24}$, $\\tilde\\theta_{34}$ (anti-neutrino channel) yields explicit probabilities in which the channel-dependent suppression by $\\sin\\theta_{j4}$ is controlled by the positions of zeros in the $S$ matrix. The resulting hierarchy is that the unsuppressed resonance effects live in $P(\\nu_e\\to\\nu_e)$, $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\mu)$, and $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\tau)$, with $\\bar\\nu_e$ decoupling at leading order, so the cascade events originate from three paths. A byproduct is that all computed probabilities respect $T$-invariance: the matter-affected phases equal their vacuum counterparts, so the phases cancel.","pith_inferences":["A concrete testable extension would be to run a neutrino-telescope cascade simulation at 1–10 TeV with a three-component template fit; if the enhancement, depletion, and flat components do not separate cleanly, the proposed fit would need a more elaborate treatment of neutrino/anti-neutrino event separation.","The paper leaves open whether the local resonance approximation survives an exact treatment of the overlapping 2-4 and 3-4 crossings; comparing the analytic $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\mu)$ to an exact numerical diagonalization in the overlap region would quantify the risk to the unsuppressed on-peak claim.","If the $T$-invariance found here is not an artifact of the sequential diagonalization, then any observed asymmetry between $\\bar\\nu_\\mu\\to\\bar\\nu_\\tau$ and $\\bar\\nu_\\tau\\to\\bar\\nu_\\mu$ in this energy band would be evidence for physics beyond the (3+1) model or for the breakdown of the approximation.","The same machinery can be adapted to inverted mass ordering by swapping the order of the 2-4 and 3-4 rotations; the qualitative prediction is that the three-path hierarchy persists with the roles of the two crossings exchanged."],"forward_implications":["In the neutrino channel, track events from $P(\\nu_\\mu\\to\\nu_\\mu)$ are a weak probe of the SA resonance, while the unsuppressed $P(\\nu_e\\to\\nu_e)$ depletion in the subdominant $\\nu_e$ flux is the neutrino-side signal.","In the anti-neutrino channel, both the track channel $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\mu)$ and the cascade appearance channel $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\tau)$ carry unsuppressed resonance enhancement and are the golden observables.","Cascade events have three origins with distinct resonance behaviours — enhancement in $\\bar\\nu_\\mu\\to\\bar\\nu_\\tau$, depletion in $\\nu_e\\to\\nu_e$, and no effect in $\\bar\\nu_e\\to\\bar\\nu_e$ — so a three-component fit can separate the contributions.","The computed probabilities are $T$-invariant in the SA region, with all CP-like phases cancelling, so this symmetry can be used as a consistency check on the framework.","The effective theory reproduces the full (3+1) numerical oscillation code to about 1% accuracy over the SA resonance region, so the analytic hierarchy is quantitatively usable."],"supporting_citations":[{"why":"Established the sterile-active resonance in a neutrino telescope as a search target, fixing the physical scenario the paper analyzes.","marker":"[24]"},{"why":"Supplies the matter-perturbation method used to embed and then freeze the standard three-flavor oscillations in the SA region.","marker":"[48]"},{"why":"Pointed out the leading-order decoupling of the $\\bar\\nu_e$ row in the anti-neutrino channel, which the paper's $S$-matrix structure reproduces.","marker":"[35]"},{"why":"Showed that $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\tau)$ is enhanced when both $\\theta_{24}$ and $\\theta_{34}$ are nonzero, motivating the cascade analysis.","marker":"[39]"},{"why":"Provided the track/cascade event classification and cascade-appearance strategy that the paper's phenomenology adopts.","marker":"[43]"},{"why":"Recent IceCube search whose 95% CL closed contour motivates the parameter region and the comparison of resonance channels.","marker":"[32]"}],"fun_headline_variants":["Sterile resonance maps three cascade paths","Three paths dominate sterile-active resonance","Global picture: sterile resonance's three paths","Resonance effects pinpointed to three cascade paths"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the anti-neutrino 2-4 and 3-4 level crossings can be diagonalized sequentially, one at a time, despite overlapping over a wide fraction of parameter space; the unsuppressed on-peak probabilities in the anti-neutrino channel are derived under that local resonance approximation.","fun_headline_variants_meta":{"raw":{"variants":["Sterile resonance maps three cascade paths","Three paths dominate sterile-active resonance","Global picture: sterile resonance's three paths","Resonance effects pinpointed to three cascade paths"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000307,"raw_usage":{"total_tokens":1815,"prompt_tokens":1059,"completion_tokens":756,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":675,"completion_tokens_details":{"reasoning_tokens":702}},"tokens_in":675,"tokens_out":756,"duration_ms":7051,"temperature":1.0,"reasoning_tokens":702,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:36:44.380696+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically diagonalize the exact (3+1) Hamiltonian in the anti-neutrino channel over the parameter region where the 2-4 and 3-4 resonances overlap, and compare the exact $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\mu)$ and $P(\\bar\\nu_\\mu\\to\\bar\\nu_\\tau)$ with the analytic local-resonance expressions; if the discrepancy exceeds the few-percent level or the on-peak unsuppression disappears, the central hierarchy claim fails.","supporting_citations":[],"review_version":1}