Zee-Burst: A New Probe of Neutrino Non-Standard Interactions at IceCube

TL;DR

This paper introduces a novel method to detect neutrino non-standard interactions using ultra-high energy neutrino data at IceCube, focusing on the Zee model which predicts observable resonance features from light charged scalars.

Contribution

It proposes using the Zee model's charged scalars to identify NSI effects in UHE neutrino spectra at IceCube and IceCube-Gen2, a new approach in neutrino physics.

Findings

Potential to observe Glashow-like resonance features in UHE neutrino spectrum.

Probing a significant portion of NSI parameter space with current and future detectors.

Light charged scalars in the Zee model can induce detectable NSI effects.

Abstract

We propose a new way to probe non-standard interactions (NSI) of neutrinos with matter using the ultra-high energy (UHE) neutrino data at current and future neutrino telescopes. We consider the Zee model of radiative neutrino mass generation as a prototype, which allows two charged scalars -- one $SU(2)_L$-doublet and one singlet, both being leptophilic, to be as light as 100 GeV, thereby inducing potentially observable NSI with electrons. We show that these light charged Zee-scalars could give rise to a Glashow-like resonance feature in the UHE neutrino event spectrum at the IceCube neutrino observatory and its high-energy upgrade IceCube-Gen2, which can probe a sizable fraction of the allowed NSI parameter space.

Equations

$$h^{+}$$

$$H^{+}$$

$$\sin 2\phi =\frac{-\sqrt{2}v\mu }{m_{H^{+}}^2-m_{h^{+}}^2}\text{leavevmode},$$

$$\displaystyle-{\cal L}_{Y}\ \supset\$$

$$+Y_{\alpha \beta }{\tilde{H}}_2^i{L}_{\alpha }^j{\ell }_{\beta }^c{ϵ}_{ij}+H.c.,$$

$$M_{\nu }=\kappa (f{M}_{\ell }Y+Y^T{M}_{\ell }{f}^T),$$

$$\kappa =\frac{1}{16{\pi }^2}\sin 2\phi \log \left(\frac{m_{h^{+}}^2}{m_{H^{+}}^2}\right).$$

$${\mathcal{L}}_Y\supset Y_{\alpha \beta }(h^{-}\sin \phi +H^{-}\cos \phi ){\nu }_{\alpha }{\ell }_{\beta }^c+H.c.$$

$$N_i=T\int d\Omega {\int }_{E_i^{\min }}^{E_i^{\max }}dE\sum _{\alpha }{\Phi }_{{\nu }_{\alpha }}(E)A_{{\nu }_{\alpha }}(E,\Omega ).$$

$${\sigma }_{\nu (\bar{\nu })N}^{\mathrm{CC}}\approx 3{\sigma }_{\nu (\bar{\nu })N}^{\mathrm{NC}}\simeq 2.7\times 10^{-36}\mathrm{c}{\mathrm{m}}^2{\left(\frac{E_{\nu }}{\mathrm{GeV}}\right)}^{0.4}.$$

$$\displaystyle\sigma_{\rm Glashow}(s)\ =\$$

$$\times \frac{s/m_W^2}{(s-m_W^2{)}^2+(m_W{\Gamma }_W{)}^2},$$

$$\displaystyle\sigma_{\rm Zee}(s)\ =\$$

$$\times \frac{s/m_X^2}{(s-m_X^2{)}^2+(m_X{\Gamma }_X{)}^2},$$

$${\epsilon }_{\alpha \beta }=\frac{Y_{\alpha e}{Y}_{\beta e}^{\star }}{4\sqrt{2}{G}_F}\left(\frac{{\sin }^2\phi }{m_{h^{+}}^2}+\frac{{\cos }^2\phi }{m_{H^{+}}^2}\right),$$

Full text

††preprint: OSU-HEP-19-05

Introduction.– The observation of ultra-high energy (UHE) neutrinos at the IceCube neutrino observatory Aartsen:2013bka ; Aartsen:2013jdh ; Aartsen:2014gkd ; Aartsen:2015zva ; Aartsen:2017mau ; Aartsen:2019kpk has commenced a new era in neutrino astrophysics. Understanding all aspects of these UHE neutrino events, including their sources, energy flux, flavor composition, propagation, and detection, is of paramount importance to both astrophysics and particle physics communities Anchordoqui:2013dnh ; Ahlers:2018mkf . A simple, single-component unbroken power-law flux $\Phi(E_{\nu})=\Phi_{0}(E_{\nu}/100\leavevmode\nobreak\ {\rm TeV})^{-\gamma}$ gives a reasonably good fit to the high-energy starting event (HESE) component of the IceCube data, with the latest best-fit values of $\Phi_{0}=(6.45^{+1.46}_{-0.46})\times 10^{-18}\ {\rm GeV}^{-1}{\rm cm}^{-2}{\rm s}^{-1}{\rm sr}^{-1}$ and $\gamma=2.89^{+0.20}_{-0.19}$ at $1\sigma$ significance Schneider:2019ayi . Any anomalous features in the observed event spectrum could potentially be used as a probe of fundamental physics. One such anomalous feature could be in the form of a new resonance. The purpose of this Letter is to show that such a new resonance can arise naturally in the popular Zee model of radiative neutrino masses Zee:1980ai ; Babu:2019mfe , which contains two charged scalars. We refer to this Zee-scalar resonance as the ‘Zee-burst’.

Within the SM, the only resonance IceCube is sensitive to is the Glashow resonance Glashow:1960zz , where electron anti-neutrinos hitting the target electrons in ice could produce an on-shell $W$-boson: $\bar{\nu}_{e}e^{-}\to W^{-}\to{\rm anything}$. The energy of the incoming neutrino required to make this resonance happen is fixed at $E_{\nu}=m_{W}^{2}/2m_{e}=6.3$ PeV. One candidate Glashow event was identified in a partially-contained PeV event (PEPE) search with deposited energy of $5.9\pm 0.18$ PeV Aartsen:2017mau ; Taboada:2018 , but has not been included in the event spectrum yet Aartsen:2019kpk . The non-observation of Glashow events might be still consistent with the SM expectations within the error bars, given the uncertainty in the source type ($pp$ versus $p\gamma$), as well as ($\nu_{e}$, $\nu_{\mu}$, $\nu_{\tau}$) flavor composition (1:2:0 vs 0:1:0) Bhattacharya:2011qu ; Barger:2012mz ; Biehl:2016psj ; Sahu:2016qet ; Sui:2018bbh . On the other hand, the possibility of observing a $Z$-boson resonance ($Z$-burst) at IceCube due to UHE anti-neutrinos interacting with non-relativistic relic neutrinos Weiler:1982qy is bleak, as the required incoming neutrino energy in this case turns out to be $E_{\nu}=m_{Z}^{2}/2m_{\nu}\gtrsim 10^{23}$ eV, well beyond the Greisen–Zatsepin–Kuzmin cut-off energy of $\sim 5\times 10^{19}$ eV for the UHE cosmic rays Greisen:1966jv ; Zatsepin:1966jv –the most likely progenitors of the UHE neutrinos (for related discussion, see Ref. Fodor:2001qy ). An interesting alternative is the existence of secret neutrino interactions with a light (MeV-scale) $Z^{\prime}$ Araki:2014ona ; Araki:2015mya ; Kamada:2015era ; DiFranzo:2015qea or light neutrinophilic neutral scalar Ioka:2014kca ; Ng:2014pca ; Ibe:2014pja , in which case the resonance could again fall in the multi-TeV to PeV range which will be accessible at IceCube. Heavy (TeV-scale) resonances induced by neutrino-nucleon interactions mediated by exotic charged particles, such as leptoquarks Barger:2013pla ; Dey:2017ede ; Becirevic:2018uab ; Dorsner:2019vgp , or squarks in $R$-parity violating supersymmetry Carena:1998gd ; Dev:2016uxj ; Collins:2018jpg ; Dev:2019ekc . have also been discussed. In this Letter, we propose the possibility of light charged scalar resonances at IceCube, which are intimately related to neutrino mass generation Zee:1980ai , as well as observable non-standard interactions (NSI) Wolfenstein:1977ue (for a recent update, see Ref. Dev:2019anc ).

As a prototypical example, we take the Zee model Zee:1980ai – one of the most popular radiative neutrino mass models, which contains an $SU(2)_{L}$-singlet charged scalar $\eta^{\pm}$ and an $SU(2)_{L}$-doublet scalar $H_{2}$, in addition to the SM-like Higgs doublet $H_{1}$. The original version of the Zee model Zee:1980ai is fully consistent with neutrino oscillation data Herrero-Garcia:2017xdu (for explicit neutrino mass fits, see Ref. Babu:2019mfe ), although the Wolfenstein version of the model Wolfenstein:1980sy which assumes a $Z_{2}$ symmetry, thus making the diagonal entries of the neutrino mass matrix vanishing, is excluded by oscillation data Koide:2001xy ; He:2003ih . Furthermore, it was pointed out in Ref. Babu:2019mfe that both the singlet and doublet charged scalar components can be as light as $\sim 100$ GeV, while satisfying all existing theoretical and experimental constraints in both charged and neutral scalar sectors. More interestingly, such light charged scalars can lead to sizable diagonal NSI of neutrinos with electrons, with the maximum allowed values of the NSI parameters $(\varepsilon_{ee},\varepsilon_{\mu\mu},\varepsilon_{\tau\tau})=(8\%,3.8\%,43\%)$. We show here that the possibility of having a resonance feature with these light charged Zee-scalars (‘Zee burst’) provides a new probe of NSI at high-energy IceCube, complementary to the low-energy neutrino oscillation and scattering experiments.

Light charged scalars in the Zee model.– In the Higgs basis Davidson:2005cw , only the neutral component of $H_{1}$ gets a vacuum expectation value $\langle H_{1}^{0}\rangle=v\simeq 246.2$ GeV, while $H_{2}$ is parametrized as $H_{2}=(H_{2}^{+},(H_{2}^{0}+iA^{0})/\sqrt{2})$. The charged scalars $\{H_{2}^{+},\eta^{+}\}$ mix in the physical basis to give rise to the physical charged scalar mass eigenstates

[TABLE]

with the mixing angle $\varphi$ given by

[TABLE]

where $\mu$ is the dimensionful coefficient of the cubic term $\mu H_{1}^{i}H_{2}^{j}\epsilon_{ij}\eta^{-}$ in the scalar potential, with $\{i,j\}$ being the $SU(2)_{L}$ indices and $\epsilon_{ij}$ being the $SU(2)_{L}$ antisymmetric tensor.

The leptonic Yukawa couplings are given by the Lagrangian

[TABLE]

where $\{\alpha,\beta\}$ are flavor indices, $\ell^{c}$ denotes the left-handed antilepton fields, and $\widetilde{H}_{a}=i\tau_{2}H^{\star}_{a}$ ($a=1,2)$ with $\tau_{2}$ being the second Pauli matrix. The neutrino mass is generated at one-loop level and is given by

[TABLE]

where $M_{\ell}=\widetilde{Y}v/\sqrt{2}$ is the charged lepton mass matrix and $\kappa$ is a one-loop factor given by

[TABLE]

According to Eq. (4), the product of the Yukawa couplings $f$ and $Y$ is constrained by the neutrino oscillation data, which allows for only one of these couplings to be of order one. We will adopt the choice $Y\sim{\cal O}(1)$ and $f\ll 1$, which maximizes the neutrino NSI in the model Babu:2019mfe .

For the IceCube phenomenology, we are specifically interested in the light charged scalar scenario. This is confronted with several theoretical and experimental constraints, such as charge breaking minima, electroweak precision tests, charged lepton flavor violation (cLFV), collider constraints from LEP and LHC, lepton universality tests and monophoton constraints. It was shown Babu:2019mfe that both $h^{+}$ and $H^{+}$ charged scalars can be as light as 100 GeV, while satisfying all these constraints. The main constraints for light charged scalars come from direct searches at LEP, which are applicable as long as $Y_{\alpha e}\neq 0$ for any flavor $\alpha$. More stringent limits from lepton universality tests in $W$ decays LEP:2003aa will apply if $Y_{ee}\neq 0$, restricting the charged scalars masses to above 130 GeV Babu:2019mfe . In what follows, we will consider the scenario where $Y_{\tau e}\neq 0$ and $Y_{\alpha\tau}\neq 0$ for $\alpha=e$ or $\mu$, which satisfies all constraints for $m_{h^{+}}=100$ GeV, and at the same time, allows for the largest NSI effect.

Signature at IceCube.– Expanding the last term in Eq. (3), we get

[TABLE]

For $\beta=e$, this will induce neutrino-electron interactions mediated by the charged scalars $h^{-}$ and $H^{-}$. For $E_{\nu}=m_{h^{-}(H^{-})}^{2}/2m_{e}$, this will lead to an $h^{-}(H^{-})$-resonance (Zee-burst) at IceCube. There is no interference with the SM Glashow process (even for $\alpha=e$), because the Zee burst involves only right-handed electrons. Thus, depending on the mass spectrum of $h^{-}$ and $H^{-}$, we would expect either one or two additional resonance peaks in the IceCube energy spectrum. We will consider two benchmark scenarios: (i) $m_{h^{-}}\approx m_{H^{-}}$, so that the two peaks are indistinguishable, i.e. contribute to the same energy bin, and (ii) $\Delta m_{h}\equiv m_{H^{-}}-m_{h^{-}}=30$ GeV, so that the two peaks are distinguishable (i.e. their dominant contributions fall in different energy bins).

To estimate the modification to the event spectrum, we compute the number of events in a given energy bin $i$ as

[TABLE]

Here $T$ is the exposure time for which we use $T_{0}=2653$ days, corresponding to 7.5 years of live data taking at IceCube Aartsen:2019kpk ; $\Omega$ is the solid angle of coverage and we integrate over the whole sky; $E$ is the electromagnetic-equivalent deposited energy which is an approximately linear function of the incoming neutrino energy Palladino:2018evm ; the limits of the energy integration $E_{i}^{\rm min}$ and $E_{i}^{\rm max}$ give the size of the $i$th deposited energy bin over which the expected number of events is being calculated; $\Phi_{\nu_{\alpha}}(E)$ is the differential astrophysical neutrino+anti-neutrino flux for flavor $\alpha$, for which we use a simple, single-component unbroken power-law, isotropic flux $\Phi(E_{\nu})=\Phi_{0}(E_{\nu}/E_{0})^{-\gamma}$ with the IceCube best-fit values of $\Phi_{0}=6.45\times 10^{-18}\ {\rm GeV}^{-1}{\rm cm}^{-2}{\rm s}^{-1}{\rm sr}^{-1}$ and $\gamma=2.89$ Schneider:2019ayi ; and $A_{\nu_{\alpha}}$ is the effective area per energy per solid angle for the neutrino flavor $\nu_{\alpha}$, which includes the effective neutrino-matter cross section, number density of target nucleons/electrons and acceptance rates for the shower and track events. In presence of new interactions as in Eq. (6), only the neutrino-electron cross section gets modified, which in turn affects the effective area. For the SM interactions only, we use the publicly available flavor-dependent effective area integrated over solid angle from Ref. Aartsen:2017mau (for 2078 days of IceCube data), along with a $67\%$ increase in the acceptance (for 2653 days of data) Aartsen:2019epb . In presence of non-SM interactions as in Eq. (6), we rescale the effective area accordingly by taking the ratio of the cross sections, assuming that the acceptance remains the same.

In the SM, neutrinos interact with nucleons via charged- and neutral-current processes. In the energy range of interest, the corresponding deep inelastic scattering cross sections can be approximated by Gandhi:1995tf

[TABLE]

In addition, there are subdominant antineutrino-electron interactions, except in the energy range of 4.6–7.6 PeV, when the $\bar{\nu}_{e}$–$e^{-}$ interaction becomes important due to the Glashow resonance Glashow:1960zz . In the vicinity of the resonance, the dominant piece of the cross section can be expressed by a Breit–Wigner distribution as Barger:2014iua :

[TABLE]

where $s=2m_{e}E_{\nu}$ and $\Gamma_{W}$ is the total width of the $W$ boson with ${\rm BR}(W^{-}\to\bar{\nu}_{e}e^{-})=10.7\%$ and ${\rm BR}(W^{-}\to{\rm had})=67.4\%$ Tanabashi:2018oca . At resonance, Eq. (Zee-Burst: A New Probe of Neutrino Non-Standard Interactions at IceCube) gives $\sigma_{\rm Glashow}(E_{\nu}=6.3\>{\rm PeV})=3.4\times 10^{-31}\rm\>cm^{2}$, about 240 times larger than $\sigma_{\nu(\bar{\nu})N}^{\rm CC}(E_{\nu}=6.3\>{\rm PeV})\approx 1.4\times 10^{-33}\>\rm cm^{2}$. However, due to the narrowness of the resonance and the $E_{\nu}^{-\gamma}$ nature of the astrophysical neutrino flux, the ratio of the reconstructed events between the resonance-induced $\bar{\nu}_{e}$-$e$ and non-resonant $\nu(\bar{\nu})$-$N$ interactions is not so pronounced in the event spectrum, as shown by the red-shaded histograms in Fig. 1. For instance, for $E_{\nu}>4$ PeV, $N_{\rm Res}/N_{\rm non-Res}\sim 2.05$ giving a total of about 0.3 events in the Glashow bin for the IceCube best-fit flux. Also shown in Fig. 1 (gray shaded) are the total expected atmospheric background (from atmospheric muons and neutrinos, as well as the charm contribution) and the 7.5 year IceCube data Schneider:2019ayi . The vertical line at 60 TeV denotes the low-energy cutoff for the HESE analysis, i.e. the bins below this energy are not considered in the fitting process.

Now in presence of light charged scalars, we expect a new resonance for $\bar{\nu}_{\alpha}e^{-}\to X^{-}\to\rm anything$ (where $X^{-}=h^{-},H^{-}$ for the Zee model) with a cross section similar to Eq (Zee-Burst: A New Probe of Neutrino Non-Standard Interactions at IceCube):

[TABLE]

where $\Gamma_{X}=\sum_{\alpha\beta}|Y_{\alpha\beta}|^{2}{\rm sin^{2}\varphi}\>m_{X}/16\pi$ is the total decay width of $X$. The factor of 1/3, compared to Eq. (Zee-Burst: A New Probe of Neutrino Non-Standard Interactions at IceCube), is due to the difference in the degrees of polarization between scalar and vector bosons.

In Fig. 1, we consider a benchmark case with $m_{h^{-}}\approx m_{H^{-}}=100$ GeV, so that the two new resonances due to $h^{-}$ and $H^{-}$ coincide, and thus, maximize the effect in the bin containing the resonance energy $E_{\nu}=m_{h^{-}}^{2}/2m_{e}$, as shown by the light, medium and dark blue-shaded histograms corresponding to three illustrative values of $Y_{\tau e}=1,0.5,0.25$ respectively. The excess events due to this new resonance mostly populate the energy bins between 7.6–12.9 PeV, distinguishable from those dominated by the Glashow resonance bin (4.6–7.6 PeV), and the effect is more pronounced for larger Yukawa couplings, as expected from Eq. (10). Here we have taken the maximal mixing $\varphi=\pi/4$ and ${\rm BR}(h^{-}\to\bar{\nu}_{\tau}e)=60\%$, ${\rm BR}(h^{-}\to\bar{\nu}_{\beta}\tau)=40\%$ (with $\beta=e$ or $\mu$) for a fixed $Y_{\tau e}$ given above and accordingly chosen $Y_{\beta\tau}$, while all other Yukawa couplings $Y_{\alpha\beta}$ are taken to be much smaller than one to satisfy the cLFV constraints Babu:2019mfe . Note that as we increase the mass difference $\Delta m_{h}\equiv m_{H^{-}}-m_{h^{-}}$, the two peaks start populating different bins, but because of the falling power-law flux, the effect is more pronounced in the smallest resonance energy bin. Also note that we cannot make $\Delta m_{h}$ exactly zero, otherwise the neutrino mass vanishes [cf. Eq. (5)].

From Fig. 1, it is clear that for a given charged scalar mass $m_{h^{-}}$, the Yukawa coupling $Y_{\tau e}$ cannot be made arbitrarily large without spoiling the best-fit to the observed IceCube HESE data. We can use this fact to derive new IceCube constraints in the $m_{h^{-}}-Y_{\tau e}$ plane, as shown in Fig. 2 by the thick black contours. The curve labeled ‘IC 1$T_{0}$’ represents the parameter set which would give rise to one event when summed over the last three bins considered by IceCube best-fit ($4.6<E_{\nu}/{\rm PeV}<10$) with the current exposure $T_{0}=2653$ days Schneider:2019ayi , and the other curves are with increased exposures of $2T_{0}$, $4T_{0}$, $10T_{0}$ and $50T_{0}$ respectively, keeping the other parameters in Eq. (7) the same. The left panel is for $m_{h+}\approx m_{H^{+}}$ and the right panel is for $m_{H^{+}}-m_{h^{+}}=30$ GeV. This explains the appearance of one ‘dip’ in the left panel (corresponding to one resonance for $h^{-}$ and $H^{-}$ combined) and two ‘dips’ in the right panel (corresponding to two distinct resonances for $h^{-}$ and $H^{-}$).

Probing NSI.– The same Yukawa interactions in Eq. (6) lead to neutrino NSI with electrons, given by Babu:2019mfe

[TABLE]

where $G_{F}$ is the Fermi coupling constant. In Fig. 2, we show the predictions for $\varepsilon_{\tau\tau}$ by thin black dotted contours. Here again we have taken the maximal mixing case with $\varphi=\pi/4$ to get the largest possible NSI. The shaded regions are all excluded: blue shaded by direct LEP searches lepsusy ; Abbiendi:2013hk and lepton universality (LU) tests in tau decays Tanabashi:2018oca ; green shaded by LEP dilepton searches LEP:2003aa ; Abbiendi:2003dh ; purple shaded (dashed) by LEP monophoton searches off (on) $Z$-pole Acciarri:1998vf ; Achard:2003tx ; red shaded by BOREXINO Agarwalla:2019smc , orange shaded by global fit to neutrino oscillation plus COHERENT data Esteban:2018ppq , and brown shaded by IceCube atmospheric neutrino data Esmaili:2013fva ; Day:2016shw . For more details on these exclusion regions, see Ref. Babu:2019mfe . Note that the atmospheric neutrino data only constrains $|\varepsilon_{\tau\tau}-\varepsilon_{\mu\mu}|<9.3\%$ Esmaili:2013fva ; Day:2016shw , which in the Zee model is equivalent to a bound on $\varepsilon_{\tau\tau}$ itself, because both $\varepsilon_{\tau\tau}$ and $\varepsilon_{\mu\mu}$ cannot be large simultaneously due to stringent cLFV constraints. One can do similar analysis for other $\varepsilon_{\alpha\beta}$, which are however restricted to be less than a few % Babu:2019mfe , and hence, are not so promising for IceCube.

We should comment here that the LEP dilepton constraints LEP:2003aa shown in Fig. 2 (green shaded region) are equally applicable to the extra neutral CP-even and odd scalars ($H,A$) present in the Zee model, since they could modify the $e^{+}e^{-}\to\ell_{\alpha}^{+}\ell_{\alpha}^{-}$ cross section via $t$-channel mediation through the Yukawa couplings $Y_{\alpha e}$. Moreover, these neutral scalars are required to be quasi-degenerate with the doublet charged scalar $H^{+}$ in order to satisfy the electroweak $T$-parameter constraint Babu:2019mfe .

From Fig. 2, we see that the existing constraints on NSI are stronger than the current sensitivity of high-energy IceCube data. However, the (non)observation of a resonance-like feature in the future IceCube HESE data could provide a complementary probe of the allowed NSI parameter space, which can even supersede the future DUNE sensitivities (shown by the upper and lower blue solid lines for 300 and 850 kt.MW.yr exposures, respectively dev_pondd ). We note here that an exposure of $10T_{0}$ does not necessarily require 75 years of IceCube running, as a number of factors could improve the conservative projected IceCube limits shown here in a non-linear fashion. For instance, the future data in all the bins may not scale proportionately to the current data and may turn out to be in better agreement with the SM prediction, thus restricting even further any room for new physics contribution. Similarly, the energy-dependent acceptance rate might improve in the future (as it did by 67% from two to seven years of data Aartsen:2019epb ), thereby increasing the effective area, and hence, the ‘effective’ exposure time defined here at a rate faster than linear. Finally, the proposed IceCube-Gen2 with 10 km3 detector volume Aartsen:2014njl could increase the total effective exposure by about an order of magnitude. At the very least, combining IceCube data with the future KM3NeT data Adrian-Martinez:2016fdl could increase the effective exposure by a factor of two.

Before concluding, we remark that for heavier charged scalars, the resonance energy will be shifted to higher values at which IceCube will become less sensitive, given an isotropic power-law spectrum. However, if there exists powerful transient sources of UHE neutrinos, then IceCube, as well as current and next-generation radio-Cherenkov neutrino detectors, such as ARA Allison:2011wk , ARIANNA Barwick:2014pca , ANITA Gorham:2019guw , GNO Avva:2016ggs and RNO Aguilar:2019jay , could be sensitive to electrophilic charged scalars up to a TeV or so (corresponding to the resonance energy of EeV), as might occur e.g. in left-right symmetric model Boyarkin:2017rte . The possibility of a larger flux at higher energies, together with better energy resolution of the IceCube detectors, might help distinguishing the degenerate versus non-degenerate charged-scalar mass spectrum by exploiting the ‘dip’ features in Fig. 2.

Conclusion.– We have proposed a new way to probe light charged scalars using a Glashow-like resonance feature in the ultra-high energy neutrino data at IceCube and its future extensions. The same interactions that lead to the new signature at IceCube also give rise to observable non-standard interactions of neutrinos with matter, so that the UHE neutrinos provide a complementary probe of NSI. Taking the popular Zee model of radiative neutrino mass as a prototypical example, we have provided an explicit realization of this idea.