The recent IceCube publication claims the observation of cosmic neutrinos with energies down to $\sim 10$ TeV, reinforcing the growing evidence that the neutrino flux in the 10-100 TeV range is unexpectedly large. Any conceivable source of these neutrinos must also produce a $\gamma$-ray flux which degrades in energy en route to the Earth and contributes to the extragalactic $\gamma$-ray background measured by the Fermi satellite. In a quantitative multimessenger analysis, featuring minimalistic assumptions, we find a $\geq 3\sigma$ tension in the data, reaching $\sim 5\sigma$ for cosmic neutrinos extended down to $\sim 1$ TeV, interpreted as evidence for a population of hidden cosmic-ray accelerators.

Heavy sterile neutrinos with masses ${\mathcal O}(100)$ MeV mixing with active neutrinos can be produced in the core of a collapsing supernova (SN). In order to avoid an excessive energy loss, shortening the observed duration of the SN 1987A neutrino burst, we show that the active-sterile neutrino mixing angle should satisfy $\sin^2 \theta \lesssim 5 \times 10^{-7}$. For a mixing with tau flavour, this bound is much stronger than the ones from laboratory searches. Moreover, we show that in the viable parameter space the decay of such "heavy" sterile neutrinos in the SN envelope would lead to a very energetic flux of daughter active neutrinos; if not too far below current limits, this would be detectable in large underground neutrino observatories, like Super-Kamiokande, as a (slightly time-delayed) high-energy bump in the spectrum of a forthcoming Galactic SN event.

This white paper summarizes the activities of the Brazilian community concerning dark matter physics and highlights the importance of financial support to Brazilian groups that are deeply involved in experimental endeavours. The flagships of the Brazilian dark matter program are the Cherenkov Telescope Array, DARKSIDE, SBN and LHC experiments, but we emphasize that smaller experiments such as DAMIC and CONNIE constitute important probes to dark sectors as well and should receive special attention. Small experimental projects showing the potential to probe new regions of parameter space of dark matter models are encouraged. On the theoretical and phenomenological side, some groups are devoted to astrophysical aspects such as the dark matter density profile while others explore the signature of dark matter models at colliders, direct and indirect detection experiments. In summary, the Brazilian dark matter community that was born not long ago has grown tremendously in the past years and now plays an important role in the hunt for a dark matter particle.

The ANITA experiment has registered two anomalous events that can be interpreted as $\nu_\tau$ or $\bar{\nu}_\tau$ with a very high energy of $\mathcal{O}(0.6)$~EeV emerging from deep inside the Earth. At such high energies, the Earth is opaque to neutrinos so the emergence of these neutrinos at such large zenith angles is a mystery. In our paper, we present a model that explains the two anomalous events through a $L_e -L_\tau$ gauge interaction involving two new Weyl fermions charged under the new gauge symmetry. We find that, as a bonus of the model, the lighter Weyl fermion can be a dark matter component. We discuss how the ANITA observation can be reconciled with the IceCube and Auger upper bounds. We also demonstrate how this model can be tested in future by collider experiments.

In view of the IceCube's 6-year high-energy starting events (HESE) sample, we revisit the possibility that the updated data may be better explained by a combination of neutrino fluxes from dark matter decay and an isotropic astrophysical power-law than purely by the latter. We find that the combined two-component flux qualitatively improves the fit to the observed data over a purely astrophysical one, and discuss how these updated fits compare against a similar analysis done with the 4-year HESE data. We also update fits involving dark matter decay via multiple channels, without any contribution from the astrophysical flux. We find that a DM-only explanation is not excluded by neutrino data alone. Finally, we also consider the possibility of a signal from dark matter annihilations and perform analogous analyses to the case of decays, commenting on its implications.

Different types of core-collapse supernovae (SNe) have been considered as candidate sources of high-energy cosmic neutrinos. Stripped-envelope SNe, including energetic events like hypernovae and super-luminous SNe, are of particular interest. They may harbor relativistic jets, which are capable of explaining the diversity among gamma-ray bursts (GRBs), low-luminosity GRBs, ultra-long GRBs, and broadline Type Ib/c SNe. Using the six-year IceCube data on high-energy starting events (HESEs), we perform an unbinned maximum likelihood analysis to search for spatial and temporal coincidences with 222 samples of SNe Ib/c. We find that the present data are consistent with the background only hypothesis, by which we place new upper constraints on the isotropic-equivalent energy of cosmic rays, ${\mathcal E}_{\rm cr}\lesssim{10}^{52}~{\rm erg}$, in the limit that all SNe are accompanied by on-axis jets. Our results demonstrate that not only upgoing muon neutrinos but also HESE data enable us to constrain the potential contribution of these SNe to the diffuse neutrino flux observed in IceCube. We also discuss implications for the next-generation neutrino detectors such as IceCube-Gen2.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

After the first four years of data taking, the IceCube neutrino telescope has observed 54 high-energy starting events (HESE) with deposited energies between 20 TeV and 2 PeV. The background from atmospheric muons and neutrinos is expected to be of about 20 events, all below 100 TeV, thus pointing towards the astrophysical origin of about 8 events per year in that data set. However, their precise origin remains unknown. Here, we perform a detailed analysis of this event sample (considering simultaneously the energy, hemisphere and topology of the events) by assuming two contributions for the signal events: an isotropic power-law flux and a flux from decaying heavy dark matter. We fit the mass and lifetime of the dark matter and the normalization and spectral index of an isotropic power-law flux, for various decay channels of dark matter. We find that a significant contribution from dark matter decay is always slightly favored, either to explain the excess below 100 TeV, as in the case of decays to quarks or, as in the case of neutrino channels, to explain the three multi-PeV events. Also, we consider the possibility to interpret all the data by dark matter decays only, considering various combinations of two decay channels. We show that the decaying dark matter scenario provides a better fit to HESE data than the isotropic power-law flux.

The installation of forward detectors in CMS and ATLAS turn the LHC to an effective photon-photon collider. The elastic scattering of the beam-protons via the emission of photons, which can be identified by tagging the intact protons in the forward detectors, provides a powerful diagnostic of the central production of new particles through photon-photon annihilation. In this letter we study the central production of dark matter particles and the potential of LHC to constrain the cross section of this process. By virtue of the crossing symmetry, this limit can immediately be used to constrain the production of monochromatic gamma-rays in dark matter annihilation, a smoking gun signal under investigation in indirect dark matter searches. We show that with the integrated luminosity $\mathcal{L}=30~{\rm fb}^{-1}$ in LHC at center-of-mass energy $\sqrt{s}=$ 13 TeV, for dark matter masses $\sim (50-600)$ GeV, a model-independent constraint on the cross section of dark matter annihilation to monochromatic gamma-rays at the same order of magnitude as the current Fermi-LAT and the future limits from CTA, can be obtained.

The very high energy Galactic $\gamma$-ray sky is partially opaque in the ($0.1-10$) PeV energy range. In the light of the recently detected high energy neutrino flux by IceCube, a comparable very high energy $\gamma$-ray flux is expected in any scenario with a sizable Galactic contribution to the neutrino flux. Here we elaborate on the peculiar energy and anisotropy features imposed upon these very high energy $\gamma$-rays by the absorption on the cosmic microwave background photons and Galactic interstellar light. As a notable application of our considerations, we study the prospects of probing the PeV-scale decaying DM scenario, proposed as a possible source of IceCube neutrinos, by extensive air shower (EAS) cosmic ray experiments. In particular, we show that anisotropy measurements at EAS experiments are already sensitive to $\tau_{\rm DM}\sim {\cal O}(10^{27})$~s and future measurements, using better gamma/hadron separation, can improve the limit significantly.

In the light of the new IceCube data on the (yet unidentified) astrophysical neutrino flux in the PeV and sub-PeV range, we present an update on the status of decaying dark matter interpretation of the events. In particular, we develop further the angular distribution analysis and discuss the perspectives for diagnostics. By performing various statistical tests (maximum likelihood, Kolmogorov-Smirnov and Anderson-Darling tests) we conclude that currently the data show a mild preference (below the two sigma level) for the angular distribution expected from dark matter decay vs. the isotropic distribution foreseen for a conventional astrophysical flux of extragalactic origin. Also, we briefly develop some general considerations on heavy dark matter model building and on the compatibility of the expected energy spectrum of decay products with the IceCube data, as well as with existing bounds from gamma-rays. Alternatively, assuming that the IceCube data originate from conventional astrophysical sources, we derive bounds on both decaying and annihilating dark matter for various final states. The lower limits on heavy dark matter lifetime improve by up to an order of magnitude with respect to existing constraints, definitively making these events---even if astrophysical in origin---an important tool for astroparticle physics studies.

The recent high-statistics high-energy atmospheric neutrino data collected by IceCube open a new window to probe new physics scenarios that are suppressed in lower energy neutrino experiments. In this paper we analyze the IceCube atmospheric neutrino data to constrain the Violation of Equivalence Principle (VEP) in the framework of three neutrinos with non-universal gravitational couplings. In this scenario the effect of VEP on neutrino oscillation probabilities can be parametrized by two parameters $\Delta \gamma_{21}\equiv \gamma_2-\gamma_1$ and $\Delta\gamma_{31}\equiv \gamma_3-\gamma_1$, where $\gamma_i$'s denote the coupling of neutrino mass eigenstates to gravitational field. By analyzing the latest muon-tracks data sets of IceCube-40 and IceCube-79, besides providing the 2D allowed regions in $(\phi\Delta\gamma_{21},\phi\Delta\gamma_{31})$ plane, we obtain the upper limits $|\phi\Delta\gamma_{21}| < 9.1\times 10^{-27}$ (at 90\% C.L.) which improves the previous limit by $\sim4$ orders of magnitude and $|\phi\Delta\gamma_{31}| \lesssim 6\times 10^{-27}$ (at 90\% C.L.) which improves the current limit by $\sim1$ order of magnitude. Also we discuss in detail and analytically the effect of VEP on neutrino oscillation probabilities.

We study the impact of the existence of an eV-mass scale sterile neutrino---with parameters in the ballpark of what is required to fit the laboratory anomalies---on the early time profile of the electron neutrino and antineutrino fluxes associated to a core-collapse supernova (SN). In particular, we focus on the universal feature of neutronization burst expected in the first tens of ms of the signal: provided that a detector with sufficient sensitivity is available, it is well-known that in the 3 neutrino framework the detection of the neutronization burst in neutrino channel would signal inverted mass hierarchy. This conclusion is dramatically altered in the presence of a sterile neutrino: we study here both analytically and numerically the region in parameter space where this characteristic signal disappears, mimicking normal hierarchy expectations. Conversely, the detection of a peak consistent with expectations for inverted mass hierarchy would exclude the existence of a sterile state over a much wider parameter space than what required by laboratory anomalies fits, or even probed by detectors coming on-line in the near future. Additionally, we show the peculiar alteration in the energy-time double differential flux, with a delayed peak appearing for kinematical reasons, which might offer a remarkable signature in case of favorable parameters and for a high statistics detection of a Galactic SN. We also comment on additional potentially interesting effects in the electron antineutrino channel, if more than one angle in the active-sterile sector is non-vanishing. As an ancillary result that we derived in the technical resolution of the equations, in an appendix we report the Cayley-Hamilton formalism for the evolution of a four neutrino system in matter, generalizing existing results in the literature.

We study oscillations of the high energy atmospheric neutrinos in the Earth into sterile neutrinos with the eV-scale mass. The MSW resonance and parametric enhancement of the $\bar{\nu}_\mu\to\bar{\nu}_s$ oscillations lead to distortion of the zenith angle distribution of the muon-track events which can be observed by IceCube. Due to matter effect, the IceCube signal depends not only on the mixing element $U_{\mu 4}$ relevant for LSND and MiniBooNE but also on $U_{\tau 4}$ and the CP-violating phase $\delta_{24}$. We show that the case with $U_{\tau 4} = \delta_{24} = 0$ leads to the weakest IceCube signal and therefore should be used to bound $U_{\mu 4}$. We compute the zenith angle distributions of the $\nu_\mu-$events for different energy intervals in the range (0.1 - 10) TeV and find that inclusion of the energy information (binning in energy) improves the sensitivity to $\nu_s$ drastically. We estimate that with already collected (during 3 - 4 years) IceCube statistics the bound $|U_{\mu 4}|^2 < 0.01$ ($99\%$ C.L.) can be established and the mixing required by LSND and MiniBooNE can be excluded at $(4 - 6) \sigma$ confidence level.

Recent observations by IceCube, notably two PeV cascades accompanied by events at energies ~ (30-400) TeV, are clearly in excess over atmospheric background fluxes and beg for an astroparticle physics explanation. Although some models of astrophysical accelerators can account for the observations within current statistics, intriguing features in the energy and possibly angular distributions of the events make worth exploring alternatives. Here, we entertain the possibility of interpreting the data with a few PeV mass scale decaying Dark Matter, with lifetime of the order of 10^27 s. We discuss generic signatures of this scenario, including its unique energy spectrum distortion with respect to the benchmark $E_\nu^{-2}$ expectation for astrophysical sources, as well as peculiar anisotropies. A direct comparison with the data show a good match with the above-mentioned features. We further discuss possible future checks of this scenario.

Medium baseline reactor experiments (Double Chooz, Daya Bay and RENO) provide a unique opportunity to test the presence of light sterile neutrinos. We analyze the data of these experiments in the search of sterile neutrinos and also test the robustness of theta_13 determination in the presence of sterile neutrinos. We show that existence of a light sterile neutrino state improves the fit to these data moderately. We also show that the measured value of theta_{13} by these experiments is reliable even in the presence of sterile neutrinos, and the reliability owes significantly to the Daya Bay and RENO data. From the combined analysis of the data of these experiments we constrain the mixing of a sterile neutrino with mass squared difference in the range of (10^{-3}-10^{-1}) eV^2 to sin^2 2 theta_{14} <0.1 at 95 % C.L..

The atmospheric neutrino data collected by the IceCube experiment and its low-energy extension DeepCore provide a unique opportunity to probe the neutrino sector of the Standard Model. In the low energy range the experiment have observed neutrino oscillations, and the high energy data are especially sensitive to signatures of new physics in the neutrino sector. In this context, we previously demonstrated the unmatched potential of the experiment to reveal the existence of light sterile neutrinos. The studies are routinely performed in the simplest 3+1 model concentrating on disappearance of muon neutrinos of TeV energy as a result of their mixing with a sterile neutrino. We here extend this analysis to include cascade events that are secondary electromagnetic and hadronic showers produced by neutrinos of all flavors. We find that it is possible to probe the complete parameter space of 3+1 model, including the poorly constrained mixing of the sterile neutrino to tau neutrinos. We show that $\nu_\tau-\nu_s$ mixing results into a unique signature in the data that will allow IceCube to obtain constraints well below the current upper limits.

The source of Ultra High Energy Cosmic Rays (UHECR) is still an unresolved mystery. Up until recently, sources of Gamma Ray Bursts (GRBs) had been considered as a suitable source for UHECR. Within the fireball model, the UHECR produced at GRBs should be accompanied with a neutrino flux detectable at the neutrino telescope such as IceCube. Recently, IceCube has set an upper bound on the neutrino flux accompanied by GRBs about 3.7 times below the prediction. We investigate whether this deficit can be explained by the oscillation of the active neutrinos to sterile neutrinos en route from the source to the detectors within the pseudo-Dirac scenario. We then discuss the implication of this scenario for diffuse supernova relic neutrinos.

Two of the most fundamental properties of the dark matter particle, the mass and the lifetime, are only weakly constrained by the astronomical and cosmological evidence of dark matter. We derive in this paper lower limits on the lifetime of dark matter particles with masses in the range 10 TeV-10^15 TeV from the non-observation of ultrahigh energy neutrinos in the AMANDA, IceCube, Auger and ANITA experiments. For dark matter particles which produce neutrinos in a two body or a three body decay, we find that the dark matter lifetime must be longer than O(10^26-10^28) s for masses between 10 TeV and the Grand Unification scale. Finally, we also calculate, for concrete particle physics scenarios, the limits on the strength of the interactions that induce the dark matter decay.

The recent global fit of short baseline neutrino oscillation data favors the presence of one (or more) sterile neutrino state which leads to new mass splitting \Delta m^2 ~1 eV^2. We consider the effect of this new states on the evolution of neutrinos from the dark matter annihilation inside the Sun. We show that neutrinos with energy E_\nu > 100 GeV undergo resonant active-sterile oscillation which depletes the flux of neutrinos arriving at the Earth. As an example of this effect, we present the oscillation probabilities for the case of monochromatic neutrinos from the direct annihilation of dark matter particles to neutrinos and the depletion due to the presence of sterile neutrinos.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.