{"id":"d6b4c44b-9fb8-414e-b172-a68bff34ca84","arxiv_id":"2411.16356","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"IceCube can test massless sterile neutrinos lighter than active neutrinos in the inverted-ordering regime, with sensitivity to sin^2 2theta24 ~ 10^-2 at |Delta m^2_41| ~ 0.2 eV^2, a region opened by the Farzan-Hannestad mechanism.","lead":"This paper works out how IceCube could detect massless sterile neutrinos that are lighter than ordinary neutrinos, a situation usually forbidden by cosmology but allowed in a model that eases the neutrino mass bound. It presents a new sensitivity curve for this inverted ordering case and says IceCube data already start to probe the model's allowed region.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central quantitative claim—IceCube 10.7-year data already test sIO at sin^2 2θ24 ~ 10^-2—rests on an unpublished scaling of the public 1-year sNO analysis to 10.7 years, without uncertainty quantification. If that scaling is wrong, the reach claim is unsupported.","rationale":"The reader's CONDITIONAL verdict is appropriate. The core physical idea—massless sterile neutrinos in the sIO configuration with Ns ~ 50–70 can evade the cosmological bound and be probed through the MSW resonance in the neutrino channel—is coherent and builds on the established Farzan–Hannestad mechanism. The paper is transparent that the 10.7-year data are not public and that the sensitivity curve is obtained by rescaling a 1-year public analysis. However, the paper's headline numerical result depends entirely on that scaling, including the assumption that the newer event selection, backgrounds, and systematics are equivalent to the older ones, and the final curve cites a private communication. The sIO signal specifically enhances the neutrino rather than antineutrino channel, so the 1-year sNO-based validation does not fully certify the sIO reach. This is a single identifiable point of failure: if the 10.7-year selection or systematics differ, the quoted reach could shift by a factor of order two or more. This does not warrant REJECT, because the qualitative argument is solid and the limitations are acknowledged; it does warrant keeping the verdict at CONDITIONAL pending release or independent validation of the 10.7-year sensitivity.","tokens_in":16587,"tokens_out":5334,"duration_ms":53866,"concrete_test":"Have the author release the event-level simulation and 10.7-year livetime scaling (or the unpublished fit behind Ref. [107]) and recompute the sIO sensitivity with an independent code using the public 1-year IceCube data release, including a 5–10% systematic nuisance on the neutrino/antineutrino ratio and on the overall normalization. If the 90% CL boundary at |Δm^2_41| = 0.2 eV^2 moves above sin^2 2θ24 = 2 × 10^-2, the quantitative claim that IceCube already tests θs ≳ 10^-2 for Ns = 50–70 is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. III B builds the 10.7-year sensitivity by (1) reproducing the public 1-year IceCube exclusion curve from [74], (2) checking that simulated event rates fall inside the official 1-year simulation band, which is for the sterile-normal-ordering (sNO) search, and (3) scaling the rates by livetime to 10.7 years. The 10.7-year data selection, background model, and systematics from [77, 78] are not public, and the final numerical result cites Ref. [107] as a private communication. The sIO signal differs from sNO in the channel where the MSW resonance occurs (neutrino versus antineutrino), so the 1-year sNO calibration does not validate the neutrino/antineutrino event ratio or the selection response to the sIO signal. No nuisance parameters or systematic-error band are attached to the blue sensitivity curves. The qualitative MSW argument for an sIO advantage is sound; what is load-bearing is the quantitative statement that IceCube sensitivity reaches sin^2 2θ24 >~ 10^-2 at |Δm^2_41| ~ 0.2 eV^2. Without release of the 10.7-year event-level information or an independent reproduction, the specific 'already started testing' claim is not verifiable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the Farzan–Hannestad mechanism, in which massless sterile neutrinos interacting through a new U(1)_Z' gauge boson replace a fraction of the cosmic neutrino background, thereby alleviating the cosmological bound on the sum of neutrino masses. It combines the model's allowed parameter region with a non-unitarity bound and a two-generation mapping to estimate the sensitivity of IceCube atmospheric neutrino observations to a massless sterile neutrino in the sterile-inverted-ordering (sIO) case. The central numerical result is that the 10.7-year IceCube data should be sensitive to sin^2 2theta_24 >~ 10^{-2} at |Delta m^2_41| ~ 0.2 eV^2, which would test the scenario with N_s ~ 50-70 massless sterile neutrinos and flavour-universal mixing theta_s >~ 6-7 x 10^{-3}.","tokens_in":16956,"tokens_out":5396,"duration_ms":51728,"significance":"If the quantitative sensitivity claim survives scrutiny, the paper identifies a genuinely new and testable parameter region: sterile neutrinos lighter than the active ones, which are usually excluded by cosmology, become allowed in the Farzan–Hannestad framework. The qualitative MSW argument that the sIO signal is enhanced in the neutrino channel and therefore statistically more efficient than the standard sNO search is sound and clearly presented. The paper also makes its approximation scheme explicit in Eq. (18) and gives a transparent interpretation of the two-generation parameters in terms of the model's N_s and theta_s. The main weakness is that the specific reach numbers rest on a scaling of the public one-year IceCube analysis to 10.7 years, with no public 10.7-year data, no uncertainty bands, and a citation to a private communication. The conceptual framework is valuable, but the headline 'already started testing' claim is not independently verifiable from the manuscript as written.","major_comments":[{"comment":"The 10.7-year sensitivity curves, which carry the paper's central quantitative claim, are produced by scaling the public one-year IceCube analysis to 10.7 years and by citing Ref. [107], a private communication. The text states that the latest IceCube data used in Refs. [77,78] are not public, and no systematic-uncertainty band is attached to the blue curves. Because the statements 'IceCube shows a good sensitivity' and 'IceCube has already started testing the model' depend directly on these curves, the analysis needs either a public, reproducible derivation of the 10.7-year sensitivity or an explicit downgrade of the claim to a projected estimate pending public data.","section":"Sec. III B, Fig. 1"},{"comment":"The calibration against the official IceCube simulation bands is performed for the sNO case, for which the MSW resonance occurs in the antineutrino channel. The sIO signal instead has the resonance in the neutrino channel, so the event rates and the selection response to neutrino versus antineutrino events are not validated by the one-year sNO simulation bands. Consequently, the sIO sensitivity curves carry an additional source of uncertainty that is not covered by the stated calibration procedure, and this uncertainty directly affects the quoted reach at |Delta m^2_41| ~ 0.2 eV^2.","section":"Sec. III B, calibration step"},{"comment":"The mapping sin^2 2theta_24 = (Theta_s Theta_s^dagger)_mu mu and Delta m^2_41 = -mbar^2 is explicitly derived in the high-energy, small-mixing, degenerate-active-mass limit. This is a reasonable approximation for the IceCube analysis, but the paper does not assess how the reach changes when the approximation degrades, for example at the lower-energy end of the sample or for non-degenerate active masses. Since the sensitivity curves are used to draw the quantitative conclusion, a validity check or an estimate of the resulting systematic shift should be shown for the quoted parameter region.","section":"Eqs. (18)-(19)"}],"minor_comments":[{"comment":"The phrase 'the latest 10.7-year atmospheric neutrino data [107]' is misleading because Ref. [107] is a private communication and not a public data release; the text should distinguish between IceCube's public data, the author's simulation, and the private communication.","section":"Sec. III B"},{"comment":"The reference entry 'M. Maltoni (2024)' should be labelled as a private communication, following standard practice, so that its non-public status is immediately clear.","section":"Reference [107]"},{"comment":"The grey-shaded cosmological regions are said to correspond to different numbers N_s of massless sterile generations, but the curve labels in the figure do not indicate which N_s values are used; adding a legend or explicit text would improve readability.","section":"Fig. 1, right panel"},{"comment":"The symbol mbar^2 is used for the squared average active-neutrino mass, but mbar as an averaged mass is introduced only later in the text; a short definition before Eq. (18) would avoid confusion with the mass-square difference.","section":"Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the conceptual MSW argument is appealing. My main concern is the reproducibility of the central quantitative claim: it relies on a non-public scaling and a private communication, with no uncertainty quantification. I would suggest asking the author to either make the 10.7-year sensitivity derivation publicly available (for example, as a code or a detailed technical appendix) or to soften the 'already started testing' wording to reflect that this is an estimated projected reach. The referee report should maintain a constructive tone, since the qualitative idea is worth publishing once the quantitative basis is clarified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Briefly: this paper opens a genuinely new window — sterile-inverted ordering (sIO) in IceCube atmospheric data — and its qualitative argument is solid, but the quantitative reach claim rests on a private rescaling of the public 1-year analysis. Treat the 10.7-year blue curves as a plausible projection, not a measured limit.\n\nThe new pieces are real. The mapping in Eq. (19) — sin²θ24 = (ΘsΘs†)_µµ and Δm²41 = −m̄² — cleanly connects the Ns-massless-sterile model to the 2-generation oscillation fit. The sIO sensitivity curve for IceCube is, as far as I know, the first in the literature. The MSW point is physically right: in sIO the resonance is in the neutrino channel, and the neutrino cross section is larger, so the sensitivity gain over sNO is expected. I checked that the author does not oversell this: Eq. (18) is explicitly labeled approximate, valid in the high-energy small-mixing limit.\n\nThe soft spot is exactly where the stress-test puts it. The 10.7-year sensitivity is built by reproducing the 1-year exclusion curve, verifying the simulated event rates against the official sNO simulation band, and then scaling by livetime. That is a reasonable calibration for sNO, but it does not validate the neutrino/antineutrino event ratio or the selection response for an sIO signal, since the 1-year official band was for the sNO topology. No systematic uncertainty band is attached to the blue curves. The bottom-line statement that 'IceCube has already started testing the model' is thus an interpretation of a projected sensitivity, not an observation. The paper is transparent about this — it cites a private fit by Maltoni [107] and says the data are not public — but the specific numerical reach depends entirely on that private input.\n\nThe recap of the cosmological bound alleviation is not new but is well done, and the KATRIN/KamLAND-Zen side constraints on the open parameter region are a useful addition. The flavour-universal assumption is a limitation the author flags and defers to future work; fair enough.\n\nWho should read it: neutrino phenomenologists and IceCube analysts interested in the sIO channel. The paper deserves a serious referee. The new parameter region and the mapping are worth publishing; the sensitivity estimate will need refinement as official 10.7-year data become public. If I were handling it, I would send to peer review and ask for (i) a clear statement that the curves are projected sensitivity, not an observed bound, and (ii) maximum detail on the rescaling procedure, ideally with a cross-check against a public 10.7-year dataset once available. I would not desk-reject on the private-data point alone.\n\nOverall: worth a serious referee.","headline":"A genuinely new sIO sensitivity estimate for IceCube whose qualitative MSW argument is sound, but the headline reach rests on an unpublished rescaling of the 1-year public analysis.","tokens_in":17453,"tokens_out":3491,"would_cite":true,"duration_ms":32605,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that a scenario with dozens of massless sterile neutrinos can relax the cosmological bound on neutrino masses, and that IceCube's atmospheric-neutrino data are already sensitive enough to test the newly opened parameter…","keywords":["sterile neutrinos","cosmological neutrino mass bound","Farzan-Hannestad mechanism","IceCube","atmospheric neutrinos","inverted mass ordering","MSW resonance","neutrino mass sum"],"falsifier":"Release of the IceCube collaboration’s official 10.7-year atmospheric-neutrino analysis in the sterile-inverted-ordering channel would settle the central quantitative claim: if the official bound at $|\\Delta m^2_{41}| \\sim 0.2$ eV$^2$ is weaker than $\\sin^2 2\\theta_{24} \\sim 10^{-2}$, the rescaled sensitivity curve is unsupported. Alternatively, a confirmed signal in that region would require a cosmologically consistent mechanism such as the Farzan-Hannestad scenario to be viable.","tokens_in":10,"feed_emoji":"🧊","tokens_out":8212,"duration_ms":129452,"temperature":0.7,"pith_summary":"The paper revisits a 2016 proposal by Farzan and Hannestad in which a fraction of the cosmic neutrino background is replaced by massless fermions, loosening the cosmological upper bound on the sum of neutrino masses. A concrete realization introduces dozens of massless sterile neutrinos interacting through a new light gauge boson; this model remains viable with an active-sterile mixing as large as $\\theta_s \\sim 0.1$. Because the bound is relaxed, the usually forbidden case where sterile neutrinos are lighter than the active ones — the 'sterile inverted ordering' — becomes open, and the paper computes, for the first time, IceCube's sensitivity to that case. The numerical result is that 10.7 years of IceCube data can already probe mixing $\\sin^2 2\\theta_{24} \\gtrsim 10^{-2}$ at mass-squared difference $\\sim 0.2$ eV$^2$, which corresponds to testing the model with 50–70 massless sterile species and average active mass 0.3–0.5 eV. If correct, the scenario gives a concrete, lab-testable target for a possible future conflict between cosmological and laboratory neutrino-mass measurements.","feed_headline":"Massless sterile neutrinos enter IceCube's reach","feed_subtitle":"A scenario that eases the cosmic mass bound predicts a new signal that 10.7 years of IceCube data can already probe.","key_machinery":"The central object is the 2-generation MSW-affected disappearance probability for atmospheric muon neutrinos, $P_{\\mu\\mu} \\simeq 1 - 4(\\Theta_s \\Theta_s^\\dagger)_{\\mu\\mu} \\left[\\frac{\\bar{m}^2}{\\bar{m}^2 + a_{\\text{NC}}}\\right]^2 \\sin^2\\frac{(\\bar{m}^2 + a_{\\text{NC}})L}{4E}$, which links the model's effective active-sterile mixing $(\\Theta_s \\Theta_s^\\dagger)_{\\mu\\mu}$ to the fit parameters $\\sin^2 2\\theta_{24}$ and $\\Delta m^2_{41} = -\\bar{m}^2$. The identity that carries the argument is the mapping between the $N_s$-flavoured model and the 2-generation fit: $\\sin^2 2\\theta_{24} = (\\Theta_s \\Theta_s^\\dagger)_{\\mu\\mu}$ with flavour-universal mixing $(\\Theta_s \\Theta_s^\\dagger)_{\\alpha\\beta} = N_s \\theta_s^2$, which lets the author interpret IceCube's $\\theta_{24}$ sensitivity as a test of the model's $\\theta_s$ for a given $N_s$. The MSW matter potential $a_{\\text{NC}} = -\\sqrt{2} G_F n_n E$ is what selects the neutrino channel for enhancement in the sIO case, which is the reason the sensitivity is better than in the ordinary normal-ordering search.","core_discovery":"The paper's central claim is that the inverted-ordering sterile neutrino (massless sterile states with $m_{\\text{active}}^2 > m_{\\text{sterile}}^2 = 0$) — normally ruled out by cosmology because it would imply $\\sum m_i \\sim 1$–3 eV — becomes viable once the Farzan-Hannestad mechanism is invoked, and that IceCube's atmospheric-neutrino observation can test it. Using a 2-generation fit with parameters $\\sin^2 2\\theta_{24}$ and $\\Delta m^2_{41}$, the author calibrates a simulation against the public 1-year IceCube analysis and scales it to 10.7 years, finding sensitivity to $\\sin^2 2\\theta_{24} \\gtrsim 10^{-2}$ at $|\\Delta m^2_{41}| \\sim 0.2$ eV$^2$ in the sterile-inverted-ordering case. With the flavour-universal identification $\\sin^2 \\theta_{24} = (\\Theta_s \\Theta_s^\\dagger)_{\\mu\\mu} = N_s \\theta_s^2$ and $\\Delta m^2_{41} = -\\bar{m}^2$, this translates to sensitivity to $\\theta_s \\gtrsim 6$–$7 \\times 10^{-3}$ for $N_s = 50$–70 and $\\bar{m} \\sim 0.3$–0.5 eV, meaning IceCube has already started testing the model. The better sensitivity in the sIO case comes from the MSW resonance occurring in the neutrino (rather than antineutrino) channel, where the detection cross-section is larger.","pith_inferences":["If the official IceCube 10.7-year analysis is released, the sIO sensitivity curve can be checked directly; this would also calibrate the author's rescaling procedure for future searches for lighter sterile states.","The flavour-universal assumption is a strong simplification; the same phenomenological approach could map flavour-dependent mixings onto the 2-generation fit, potentially opening or closing specific flavour channels.","A detection of an sIO signal would not by itself identify the Farzan-Hannestad mechanism; corroborating evidence from cosmological probes, such as a modified effective number of relativistic species, would be needed to distinguish it from other light-sterile scenarios."],"forward_implications":["If the sensitivity estimate is right, current IceCube data already constrain or exclude the flavour-universal model for $N_s \\sim 50$–70 with $\\theta_s \\gtrsim 10^{-2}$ and $\\bar{m} \\sim 0.3$–0.5 eV.","A null result in this sIO search would push the model toward either smaller mixing or a larger number of sterile species, with $N_s \\gtrsim 90$–150 needed if the tighter DESI-based cosmological bound is adopted.","A positive signal would motivate a reconsideration of the cosmological neutrino-mass bound within extended cosmologies such as the Farzan-Hannestad scenario.","The sIO search region overlaps with the KATRIN beta-decay bound, which already excludes the opened parameter region for $N_s \\gtrsim 70$, so combining lab and cosmological probes can cross-check the model.","The same MSW-in-the-neutrino-channel logic applies to future atmospheric-neutrino searches, potentially extending the reach beyond IceCube's current configuration."],"supporting_citations":[{"why":"Proposes the Farzan-Hannestad mechanism and the formula (Eq. 1) that relaxes the cosmological bound on the sum of neutrino masses.","marker":"[52]"},{"why":"Realizes the mechanism with massless sterile neutrinos and a U(1)′ gauge boson, and provides the allowed parameter region including the Addendum’s relaxed upper bound on the active-sterile mixing.","marker":"[51]"},{"why":"The public 1-year IceCube analysis used to calibrate the author’s simulation setup before rescaling to 10.7 years.","marker":"[74]"},{"why":"Private communication providing the 10.7-year IceCube atmospheric-neutrino sensitivity curve used for the main numerical result.","marker":"[107]"},{"why":"Identifies that the MSW resonance in the inverted-ordering case occurs in the neutrino channel, giving a detection cross-section advantage for the sIO search.","marker":"[65]"},{"why":"DESI BAO data yielding the tight cosmological bound that the scenario must alleviate, and used to estimate the required number of sterile species.","marker":"[22]"},{"why":"Snowmass summary of non-unitarity bounds that constrains the effective mixing $N_s \\theta_s^2$ and sets the reference constraint on $\\theta_s$.","marker":"[103]"},{"why":"KATRIN beta-decay bound that excludes part of the opened parameter region for $N_s \\gtrsim 70$ and sets the upper limit on the average active neutrino mass.","marker":"[1]"}],"fun_headline_variants":["IceCube can probe massless sterile neutrinos in inverted ordering","First test of inverted-ordering sterile neutrinos with IceCube","Massless sterile neutrino scenario enters IceCube's reach","Cosmic bound eased, sterile neutrinos testable at IceCube","IceCube sensitivity to massless sterile neutrinos revealed"],"cache_read_input_tokens":19584,"weakest_assumption_plain":"The quantitative sensitivity estimate assumes that scaling the public 1-year IceCube event rates to 10.7 years reproduces the official analysis, with the same data selection, backgrounds, and systematics, even though the 10.7-year data are not public and the curve rests on a private communication.","fun_headline_variants_meta":{"raw":{"variants":["IceCube can probe massless sterile neutrinos in inverted ordering","First test of inverted-ordering sterile neutrinos with IceCube","Massless sterile neutrino scenario enters IceCube's reach","Cosmic bound eased, sterile neutrinos testable at IceCube","IceCube sensitivity to massless sterile neutrinos revealed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000485,"raw_usage":{"total_tokens":2466,"prompt_tokens":1093,"completion_tokens":1373,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":1290}},"tokens_in":709,"tokens_out":1373,"duration_ms":9911,"temperature":1.0,"reasoning_tokens":1290,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:13:13.750048+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Release of the IceCube collaboration’s official 10.7-year atmospheric-neutrino analysis in the sterile-inverted-ordering channel would settle the central quantitative claim: if the official bound at $|\\Delta m^2_{41}| \\sim 0.2$ eV$^2$ is weaker than $\\sin^2 2\\theta_{24} \\sim 10^{-2}$, the rescaled sensitivity curve is unsupported. Alternatively, a confirmed signal in that region would require a cosmologically consistent mechanism such as the Farzan-Hannestad scenario to be viable.","supporting_citations":[{"cited_title":"Sterile Neutrino Oscillations: The Global Picture","cited_arxiv_id":"1303.3011","evidence_quote":"The public 1-year IceCube analysis used to calibrate the author’s simulation setup before rescaling to 10.7 years."}],"review_version":1}