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.

The ANITA balloon experiment has observed two EeV-energy, upgoing events originating from well below the horizon. This is puzzling, because (i) no Standard Model (SM) particle is expected to survive passage through Earth at such energies and incident angles, and (ii) no such events were reported by IceCube. In this paper, we address both these issues by invoking a beyond SM interpretation of the EeV events as due to the decay of a long-lived bino in the R-parity violating (RPV) supersymmetry. In particular, a TeV-scale slepton/squark can be resonantly produced through the interaction of the EeV neutrino with electrons/nucleons inside Earth, that decays to a light, long-lived bino, which survives the propagation through Earth matter before decaying back to neutrinos, leptons and/or quarks, thus producing upgoing air showers in the atmosphere near ANITA. We find that the ANITA events can be explained with a GeV-scale bino and ${\cal O}(0.1)$ RPV couplings, which are consistent with all existing high and low-energy constraints. We also find that an isotropic neutrino flux is inadequate for a beyond the SM explanation of this kind, and an anisotropic flux must be invoked. Finally, we also address the apparent tension of these observations with IceCube. Various aspects of our interpretation are testable in the near future at different frontiers, such as by the LHC, Belle II, ANITA-IV and IceCube.

We perform a combined likelihood analysis for the IceCube 6-year high-energy starting events (HESE) and 8-year throughgoing muon events above 10 TeV using a two-component neutrino flux model. The two-component flux can be motivated either from purely astrophysical sources or due to a beyond Standard Model contribution, such as decaying heavy dark matter. We find that the astrophysical plus dark matter interpretation is mildly preferred by the current data over the purely astrophysical explanation. As for the astrophysical neutrinos, we consider two different source flavor compositions corresponding to the standard pion decay and muon-damped pion decay sources. We find that the latter is slightly preferred over the former as the high-energy component, while the low-energy component does not show any such preference. We also take into account the multi-messenger gamma-ray constraints and find that our two-component fit is compatible with these constraints, whereas the single-component power-law bestfit to the HESE data is ruled out. The astrophysical plus dark matter interpretation of the two-component flux is found to be mildly preferred by the current data and the gamma-ray constraints over the purely astrophysical explanation.

The DAMPE $e^+ e^-$ excess at around 1.4 TeV could be explained in the type-II seesaw model with a scalar dark mater $D$ which is stabilized by a discrete $Z_2$ symmetry. The simplest scenario is the annihilation $DD \to H^{++} H^{--}$ followed by the subsequent decay $H^{\pm\pm} \to e^\pm e^\pm$, with both the DM and triplet scalars roughly 3 TeV with a small mass splitting. In addition to the Drell-Yan process at future 100 TeV hadron colliders, the doubly-charged components could also be produced at lepton colliders like ILC and CLIC in the off-shell mode, and mediate lepton flavor violating processes $e^+ e^- \to \ell_i^\pm \ell_j^\mp$ (with $i \neq j$). A wide range of parameter space of the type-II seesaw could be probed, which are well below the current stringent lepton flavor constraints.