{"id":"d5f5ea3e-c66d-4757-897c-52d793b7ca5a","arxiv_id":"1908.02779","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Light charged scalars in the Zee model would produce a Glashow-like resonance peak in IceCube's PeV neutrino spectrum, providing a complementary probe of neutrino non-standard interactions.","lead":"This paper proposes that lightweight charged particles from the Zee model of neutrino mass could create a resonant burst of events in the IceCube neutrino telescope, alongside the known Glashow resonance. This offers a new way to test non-standard neutrino interactions with matter using ultra-high-energy cosmic neutrinos.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Projected IceCube sensitivity uses a one-event threshold that ignores the SM background; the claim of probing a sizable NSI fraction is not statistically substantiated.","rationale":"I read the paper in good faith. The Zee-burst mechanism is physically plausible: a 100 GeV charged scalar produced in bar-nu_tau e^- collisions at ~10 PeV gives a resonant enhancement, and the same couplings generate NSI. The resonance cross section and event rate estimates are broadly consistent with the Glashow resonance formalism. The paper's most important weakness is not the physics mechanism but the statistical interpretation of the projected sensitivity. The 'one event' contours in Fig. 2 ignore the SM background, which is not negligible at these energies. At T0 the SM gives ~0.2–0.3 events in the relevant bin, and at 10 T0 it gives ~2–3 events, so a one-signal-event criterion is both uncalibrated and, at large exposure, below the mean background. A proper Poisson treatment could substantially weaken the projected reach. The reader's weakest_assumption concerned the validity of the 100 GeV parameter point from Ref. [11]; that is a valid external dependency, but it is not the single most load-bearing concern because the paper's quantitative claim fails (or at least is unsubstantiated) if the statistical analysis is inadequate, even assuming the parameter point is valid. Since the reader already issued CONDITIONAL and listed the simplified sensitivity analysis as a weakness, my concern does not change the verdict but sharpens the condition: the projections need a full statistical treatment before the 'sizable fraction' claim can be accepted.","tokens_in":11576,"tokens_out":26854,"duration_ms":289389,"concrete_test":"Recompute the IceCube sensitivity curves in Fig. 2 using a binned Poisson likelihood or Feldman–Cousins upper limits that include the SM astrophysical background in the resonance bins, using the same flux normalization, effective areas, and exposure scaling. Derive the 90% CL exclusion/discovery curves in the (m_h+, |Yτe| sinφ) plane for T0, 2T0, 4T0, 10T0, and 50T0. If the required |Yτe| at a given mass increases by more than ~30% relative to the one-event contours, the paper's claim that a sizable fraction of the allowed NSI parameter space is probed is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that IceCube-Gen2 can probe a sizable fraction of the allowed NSI parameter space—rests on the sensitivity curves in Fig. 2, which are defined as the parameter set yielding one expected signal event in the combined resonance bins. This threshold does not account for the SM astrophysical background in those same bins. For the best-fit flux, the SM already contributes roughly 0.2–0.3 events in the 7.6–12.9 PeV bin at the current exposure T0; at 10 T0 the background grows to ~2–3 events. A one-event signal threshold is not a valid discovery or exclusion statistic: it ignores Poisson fluctuations of the background, and at large exposures it lies below the mean background itself. Therefore the projected sensitivity could shift significantly when a proper binned Poisson likelihood or Feldman–Cousins treatment is applied. A related but separate concern is the reliance on Ref. [11] for the viability of the 100 GeV scalar parameter point; if that point is invalid the signature disappears. However, the statistical issue is the more load-bearing internal weakness because it directly controls whether the headline claim of probing a 'sizable fraction' of NSI parameter space is supported by the presented analysis.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter proposes a new IceCube signature, the 'Zee-burst,' arising from light charged scalars in the Zee model of radiative neutrino mass. The scalars h^+ and H^+, if as light as 100 GeV, can produce Glashow-like resonances in neutrino-electron scattering at energies E_nu = m_h^2/2m_e in the few-to-tens PeV range. The authors compute the expected event rates for benchmark parameters (m_h+ ~ m_H+ = 100 GeV, maximal mixing phi = pi/4, Y_tau e nonzero), compare the resulting spectra with 7.5 years of IceCube HESE data, and derive projected sensitivity curves in the (m_h+, |Y_tau e| sin phi) plane for exposures T0, 2T0, 4T0, 10T0, 50T0. They then map these sensitivities onto the NSI parameter epsilon_tau tau and conclude that IceCube and IceCube-Gen2 can probe a sizable fraction of the currently allowed NSI parameter space. The central claim is that the same Yukawa couplings responsible for neutrino mass and NSI produce a directly observable resonance feature at neutrino telescopes.","tokens_in":11787,"tokens_out":6587,"duration_ms":75025,"significance":"If the Zee-burst is real, it would be a genuinely new probe of the charged scalars that generate neutrino mass, and it would connect low-energy NSI bounds to a high-energy resonance search in a parameter-free way within the model. The paper is forward-looking and uses a forward calculation from the Lagrangian; it also uses publicly available effective areas and standard cross-section formulas, which is a strength. The mapping from the resonance rate to epsilon_tau tau via Eq. (11) is explicit and makes the NSI connection concrete. However, the quantitative significance of the projected reach is weakened by the simplified sensitivity definition: the curves in Fig. 2 are based on a one-event threshold that ignores the SM background and Poisson fluctuations, and the effective-area rescaling procedure is not fully specified. The model-parameter benchmark is also inherited from the companion paper Ref. [11] without reproduction. These issues affect the headline claim but are addressable in a revision.","major_comments":[{"comment":"The projected sensitivity curves labeled 'IC 1T0' through 'IC 50T0' are defined as the parameter set giving one expected signal event summed over the resonance bins, with no treatment of the SM background or of Poisson fluctuations. For the IceCube best-fit flux, the SM expectation in the relevant PeV bins is not negligible at the larger exposures: at 10 T0 the mean background in the same bins is already of order a few events, so a one-event threshold does not correspond to a valid discovery or exclusion statistic. Because the central claim that IceCube and IceCube-Gen2 can probe a sizable fraction of the allowed NSI parameter space rests directly on these curves, the analysis should be redone using a binned Poisson likelihood or a Feldman-Cousins prescription that includes the SM background and its fluctuations.","section":"Fig. 2 and 'Signature at IceCube'"},{"comment":"The text says that in the presence of the new interactions the effective area is rescaled 'by taking the ratio of the cross sections,' but the publicly available effective area from Ref. [5] includes both neutrino-nucleon and neutrino-electron interactions. Since the Zee scalars modify only the neutrino-electron piece, the paper must specify how the electron contribution is separated from the nucleon contribution before rescaling. Without this decomposition the event rates in Fig. 1 and the sensitivity curves in Fig. 2 are not reproducible, and the size of the Zee-burst excess could be either over- or underestimated depending on the assumed ratio.","section":"Eq. (7) and effective-area rescaling"},{"comment":"The benchmark point with m_h+ = 100 GeV, Y_tau e nonzero and Y_alpha tau nonzero (alpha = e or mu) is asserted to satisfy all constraints solely on the authority of Ref. [11]; the relevant LEP, lepton-universality, and cLFV limits are not reproduced or summarized in this manuscript. Since the entire Zee-burst signature disappears if this parameter point is not viable, the paper should either display the relevant exclusion regions or explicitly state that the sensitivity projection inherits all assumptions of Ref. [11] and is conditional on those constraints.","section":"Section 'Light charged scalars in the Zee model'"},{"comment":"The event-rate calculation uses a single fixed power-law flux with the IceCube best-fit values of Phi_0 and gamma, and the projected sensitivities do not propagate the uncertainties in these parameters or any systematic uncertainties in the effective area. Because the signal rate is directly proportional to the flux normalization and the background depends on the spectral index, the 'sizable fraction' claim should be accompanied by a range of projected reaches obtained by varying the flux parameters within their 1-sigma uncertainties.","section":"Eq. (7) and flux assumptions"}],"minor_comments":[{"comment":"The sentence 'squarks in R-parity violating supersymmetry [34–37]. have also been discussed' has a misplaced period and should be corrected.","section":"Introduction, paragraph 2"},{"comment":"The width formula Gamma_X = sum |Y_alpha beta|^2 sin^2 phi m_X / 16 pi is written for a generic X, but Eq. (6) shows that h^- couples with sin phi and H^- couples with cos phi. The formula should distinguish the two mass eigenstates; it is only correct at the phi = pi/4 benchmark used in the figures.","section":"Eq. (10) and width formula"},{"comment":"The text states that 'we cannot make Delta m_h exactly zero, otherwise the neutrino mass vanishes,' yet the benchmark for Fig. 1 is labeled m_h+ ~ m_H+ = 100 GeV. The size of the mass splitting used in the 'degenerate' benchmark should be quantified so the reader can assess the two-peak separation.","section":"Benchmark mass degeneracy"},{"comment":"The notation 'N_Res/N_non-Res' in the text is not defined, and the caption of Fig. 2 refers to 'thick black curves' while the figure contains multiple curve styles; the captions should define all curve labels and the meaning of the shaded regions.","section":"Fig. 1 and Fig. 2 captions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well-motivated and the proposal is interesting, but the quantitative claims in the abstract and Fig. 2 are not yet backed by a statistically valid sensitivity analysis. The one-event threshold is the most serious internal issue; a proper binned Poisson treatment could shift the projected reach substantially. I also recommend the editors ask the authors to clarify the effective-area rescaling procedure and to state explicitly the dependence on the companion paper for the benchmark parameter space."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look. The paper proposes a Glashow-like resonance from 100 GeV charged scalars in the Zee model — the Zee-burst — and connects it to NSI with electrons. As far as I can tell, that specific signature has not been proposed before; the previous literature covers Glashow, Z-burst, leptoquarks, RPV squarks, but not the Zee scalars. The calculation is straightforward: Breit-Wigner cross section, rescale IceCube effective area by cross-section ratio, compute event rates. They use public data and clearly state their flux assumptions. Credit where due: the idea is clean and the link between the resonance and neutrino mass generation is nontrivial and appealing.\n\nThe soft spots are real, but they are concentrated in the sensitivity estimate. The 'IC 1T0' curves in Fig. 2 are defined by one expected signal event in the combined resonance bins. That threshold ignores the SM astrophysical background in those bins. For the best-fit flux, the background alone is already ~0.2–0.3 events in the 7.6–12.9 PeV bin at T0, and grows to ~2–3 events at 10 T0. A one-event signal threshold is not a valid discovery or exclusion criterion, and at high exposure it sits below the mean background. So the headline claim that IceCube-Gen2 can probe a sizable NSI fraction is not yet statistically substantiated. This is fixable with a proper binned Poisson likelihood or Feldman-Cousins treatment, but it has to be done.\n\nThe other dependency is the 100 GeV parameter point taken from Ref. [11]. The paper asserts that Y_tau e != 0 and Y_alpha tau != 0 satisfy all constraints, but the proof is in the companion paper. That is a normal division of labor, but it makes the signature contingent on a result the reader has to go elsewhere to check.\n\nOverall: the physics idea is interesting and likely correct in broad strokes; the quantitative projection needs work. This is a paper for IceCube phenomenologists and anyone working on NSI in radiative mass models. It deserves a serious referee, not a desk reject, and the referee should push for a background-aware statistical analysis and an explicit statement of the Ref. [11] parameter-space validity.","headline":"A genuinely new IceCube resonance proposal from Zee-model scalars, with a real but fixable statistical weakness in the projected sensitivity.","tokens_in":12355,"tokens_out":2143,"would_cite":true,"duration_ms":23475,"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":"Light charged scalars in the Zee model of neutrino mass could create a Glashow-like resonance—dubbed the Zee-burst—in IceCube's ultra-high-energy spectrum, giving a new handle on neutrino non-standard interactions.","keywords":["Zee model","Zee-burst","neutrino non-standard interactions","IceCube","Glashow resonance","radiative neutrino mass","ultra-high-energy neutrinos","tau neutrino NSI"],"falsifier":"Search IceCube's high-energy starting events in the deposited-energy bins between 7.6 and 12.9 PeV: if the counted events with the best-fit single-power-law flux match the Standard Model prediction and no excess appears at the level expected for $m_{h^\\pm}=100$ GeV and $\\lvert Y_{\\tau e}\\rvert\\sin\\phi=0.4$, the benchmark Zee-burst scenario is excluded.","tokens_in":11370,"feed_emoji":"🧊","tokens_out":12118,"duration_ms":113745,"temperature":0.7,"pith_summary":"Ultra-high-energy neutrinos at IceCube could reveal a new resonance, the Zee-burst, if the charged scalars of the Zee radiative neutrino-mass model are light enough. The paper shows that an antineutrino scattering on a target electron can produce one of these scalars on-shell at neutrino energy $E_\\nu = m_X^2/(2m_e)$, placing a 100 GeV scalar's resonance at about 10 PeV. Because the same Yukawa couplings that generate neutrino mass at one loop also produce neutrino–electron non-standard interactions, the presence or absence of this resonance constrains the non-standard-interaction parameter space, especially $\\varepsilon_{\\tau\\tau}$, which can be as large as 43% in the model. The authors compute event spectra using IceCube's best-fit astrophysical flux and show that current 7.5-year data already start to bound the relevant couplings, while future exposures and IceCube-Gen2 could probe a sizable fraction of the allowed non-standard-interaction region.","feed_headline":"A 'Zee-burst' could add a 10 PeV resonance to IceCube's spectrum","feed_subtitle":"The scalars behind neutrino mass could add a visible bump near 10 PeV—a new window on neutrino NSI.","key_machinery":"The load-bearing object is the Zee-burst resonance: an $s$-channel process $\\bar\\nu_\\alpha e^- \\to X^- \\to$ anything, with $X^- = h^-$ or $H^-$, whose Breit–Wigner cross section is $\\sigma_{\\rm Zee}(s) = 8\\pi \\Gamma_X^2 \\, \\mathrm{BR}(X^-\\to \\bar\\nu_\\alpha e^-)\\, \\mathrm{BR}(X^-\\to \\mathrm{all}) \\, \\frac{s/m_X^2}{(s-m_X^2)^2+(m_X\\Gamma_X)^2}$. It is the scalar analogue of the Glashow resonance; the factor $8\\pi$ instead of $24\\pi$ reflects the single polarization degree of a spin-0 propagator rather than the three of the $W$. The resonance position $E_\\nu=m_X^2/(2m_e)$ converts a 100 GeV scalar mass into a bump near 9.8 PeV. The same Yukawa matrix $Y$ enters the one-loop neutrino mass formula $M_\\nu = \\kappa(f M_\\ell Y + Y^T M_\\ell f^T)$ and the non-standard-interaction parameter $\\varepsilon_{\\alpha\\beta} = \\frac{Y_{\\alpha e}Y^*_{\\beta e}}{4\\sqrt{2}G_F}\\left(\\frac{\\sin^2\\phi}{m_h^2}+\\frac{\\cos^2\\phi}{m_H^2}\\right)$, which is the link between IceCube event rates and the low-energy non-standard-interaction parameter space.","core_discovery":"The paper's central claim is that the Zee-burst, the charged-scalar analogue of the Glashow resonance, is a viable and clean high-energy probe of neutrino non-standard interactions. In the Zee model the physical charged scalars $h^\\pm$ and $H^\\pm$ couple to neutrinos and charged leptons; for an electron target the process $\\bar\\nu_\\alpha e^- \\to X^- \\to$ anything becomes resonant at $E_\\nu = m_X^2/(2m_e)$. With $m_{h^\\pm}\\simeq m_{H^\\pm}\\simeq 100$ GeV the resonance sits at $E_\\nu\\simeq 9.8$ PeV and populates the 7.6–12.9 PeV deposited-energy bins of IceCube's high-energy starting events, adjacent to the Glashow bin. Because the Zee coupling involves right-handed electrons, there is no interference with the $W$-mediated Glashow process, so the excess is a separate bump. Adopting the viable benchmark with $Y_{\\tau e}\\neq 0$, $Y_{\\alpha\\tau}\\neq0$ for $\\alpha=e,\\mu$, $Y_{ee}=0$, maximal mixing $\\phi=\\pi/4$, and branching ratios fixed by the companion analysis, the authors derive sensitivity contours in the $m_{h^\\pm}$–$\\lvert Y_{\\tau e}\\rvert\\sin\\phi$ plane and translate them into projected reach on $\\varepsilon_{\\tau\\tau}$. They conclude that observation or non-observation of the Zee-burst in IceCube and IceCube-Gen2 can probe a sizable fraction of the allowed non-standard-interaction parameter space and may supersede DUNE's projected sensitivity for $\\varepsilon_{\\tau\\tau}$.","pith_inferences":["The same resonance logic transfers to any model with a light leptophilic charged scalar: IceCube's spectrum becomes a generic map from non-standard-interaction parameters to a PeV bump, so the search is not limited to the Zee model.","If the Zee-burst is observed, the peak energy fixes the scalar mass, and the relative weights of the $h^-$ and $H^-$ contributions constrain the mixing angle $\\phi$ and the mass splitting, information that low-energy oscillation and scattering experiments cannot easily resolve.","A null search at IceCube-Gen2 combined with future oscillation and coherent-scattering limits would corner the light-scalar Zee model, because all three observables are governed by the same Yukawa matrix through the non-standard-interaction relation.","The projected reach assumes an isotropic single-power-law flux; a harder or anisotropic flux, or a dedicated point-source stacking search, could improve the sensitivity beyond the stated contours."],"forward_implications":["If the benchmark Zee-burst is real, IceCube should see an excess in the 7.6–12.9 PeV bins whose size grows with $\\lvert Y_{\\tau e}\\rvert\\sin\\phi$; the current 7.5-year data already rule out the largest couplings.","A non-observation at future exposures would exclude a substantial slice of the allowed non-standard-interaction parameter space, particularly $\\varepsilon_{\\tau\\tau}$ at the tens-of-percent level for scalar masses near 100 GeV.","The two benchmark mass patterns give distinct signatures: a degenerate $h^-/H^-$ pair produces one combined peak, while a 30 GeV splitting produces two dips, so the shape of any excess can diagnose the scalar mass spectrum.","With IceCube-Gen2's roughly tenfold exposure, or a combination with KM3NeT, the Zee-burst search becomes competitive with, and for $\\varepsilon_{\\tau\\tau}$ potentially stronger than, DUNE's projected sensitivity.","Heavier scalars shift the resonance beyond IceCube's reach for an isotropic flux, but transient astrophysical sources could make TeV-scale scalars accessible through EeV neutrinos in radio-Cherenkov detectors."],"supporting_citations":[{"why":"Introduces the Zee model with two charged scalars and one-loop radiative neutrino mass that the paper adopts as its prototype.","marker":"[10]"},{"why":"Supplies the viable 100 GeV charged-scalar benchmark satisfying LEP, lepton-universality and charged-lepton-flavor-violation constraints, plus the non-standard-interaction predictions used here.","marker":"[11]"},{"why":"Defines the Glashow resonance that the Zee-burst generalizes and that sets the Standard Model resonance baseline.","marker":"[12]"},{"why":"Provides the IceCube best-fit single-power-law astrophysical flux and 7.5-year high-energy starting event data used for all event-rate calculations.","marker":"[9]"},{"why":"Provides the flavor-dependent IceCube effective areas that the paper rescales to compute reconstructed event spectra.","marker":"[5]"},{"why":"Gives the Standard Model neutrino-nucleon deep-inelastic cross sections that form the non-resonant background.","marker":"[48]"},{"why":"Provides the Breit–Wigner form of the Glashow resonance cross section that the Zee-burst cross-section formula adapts.","marker":"[49]"},{"why":"Defines IceCube-Gen2, whose larger exposure underlies the projected future sensitivity of the Zee-burst search.","marker":"[61]"}],"fun_headline_variants":["Zee-burst: new probe of neutrino NSI at IceCube","Zee scalars could add a 10 PeV bump to IceCube's spectrum","Charged scalars may leave a Glashow-like peak in IceCube","Zee-burst could reveal neutrino NSI with a 10 PeV peak"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The signature stands on the existence of a viable Zee-model parameter point with $m_{h^\\pm}=100$ GeV, $Y_{ee}=0$, $Y_{\\tau e}\\neq 0$, and $Y_{\\alpha\\tau}\\neq0$ that satisfies all LEP, lepton-universality, charged-lepton-flavor-violation, and scalar-sector constraints; the paper takes this point from its companion paper [11].","fun_headline_variants_meta":{"raw":{"variants":["Zee-burst: new probe of neutrino NSI at IceCube","Zee scalars could add a 10 PeV bump to IceCube's spectrum","Charged scalars may leave a Glashow-like peak in IceCube","Zee-burst could reveal neutrino NSI with a 10 PeV peak"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000971,"raw_usage":{"total_tokens":4183,"prompt_tokens":1052,"completion_tokens":3131,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":668,"completion_tokens_details":{"reasoning_tokens":3045}},"tokens_in":668,"tokens_out":3131,"duration_ms":21534,"temperature":1.0,"reasoning_tokens":3045,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:34:28.240377+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Search IceCube's high-energy starting events in the deposited-energy bins between 7.6 and 12.9 PeV: if the counted events with the best-fit single-power-law flux match the Standard Model prediction and no excess appears at the level expected for $m_{h^\\pm}=100$ GeV and $\\lvert Y_{\\tau e}\\rvert\\sin\\phi=0.4$, the benchmark Zee-burst scenario is excluded.","supporting_citations":[{"cited_title":"A Theory of Lepton Number Violation, Neutrino Majorana Mass, and Oscillation,","cited_arxiv_id":null,"evidence_quote":"Introduces the Zee model with two charged scalars and one-loop radiative neutrino mass that the paper adopts as its prototype."},{"cited_title":"Resonant Scattering of Antineutrinos,","cited_arxiv_id":null,"evidence_quote":"Defines the Glashow resonance that the Zee-burst generalizes and that sets the Standard Model resonance baseline."},{"cited_title":"Characterization of the Astrophysical Diffuse Neutrino Flux with IceCube High-Energy Starting Events","cited_arxiv_id":"1907.11266","evidence_quote":"Provides the IceCube best-fit single-power-law astrophysical flux and 7.5-year high-energy starting event data used for all event-rate calculations."}],"review_version":1}