The discovery of high-energy astrophysical neutrinos by IceCube has opened a new window to the Universe. However, the origin of these neutrinos is still a mystery, and some of them could be a result of dark matter interactions such as decay. Next generation gigaton water-Cherenkov neutrino telescope, KM3NeT, is expected to offer significantly improved energy resolution in the cascade channel, and advantageous viewing condition to the Galactic Center; both important for searches of dark matter decay signals. We study the sensitivity of KM3NeT on dark matter decays by performing a mock likelihood analysis for both cascade and track type events, taking into account both angular and energy information. We find that, combining both channels, KM3NeT is expected to produce world leading limits on dark matter decay lifetime in the PeV mass range, and could test some of the dark matter hints in the current IceCube data.

The hypothesis of two different components in the high-energy neutrino flux observed with IceCube has been proposed to solve the tension among different data-sets and to account for an excess of neutrino events at 100 TeV. In addition to a standard astrophysical power-law component, the second component might be explained by a different class of astrophysical sources, or more intriguingly, might originate from decaying or annihilating dark matter. These two scenarios can be distinguished thanks to the different expected angular distributions of neutrino events. Neutrino signals from dark matter are indeed expected to have some correlation with the extended galactic dark matter halo. In this paper, we perform angular power spectrum analyses of simulated neutrino sky maps to investigate the two-component hypothesis with a contribution from dark matter. We provide current constraints and expected sensitivity to dark matter parameters for future neutrino telescopes such as IceCube-Gen2 and KM3NeT. The latter is found to be more sensitive than IceCube-Gen2 to look for a dark matter signal at low energies towards the galactic center. Finally, we show that after 10 years of data-taking, they will firmly probe the current best-fit scenario for decaying dark matter by exploiting the angular information only.

Astrophysical neutrino events have been measured in the last couple of years, which show an isotropic distribution, and the current discussion is their astrophysical origin. We use both isotropic and anisotropic components of the diffuse neutrino data to constrain the contribution of a broad number of extra-galactic source populations to the observed neutrino sky. We simulate up-going muon neutrino events by applying statistical distributions for the flux of extragalactic sources, and by Monte Carlo method we exploit the simulation for current and future IceCube, IceCube-Gen2 and KM3NeT exposures. We aim at constraining source populations by studying their angular patterns, for which we assess the angular power spectrum. We leave the characteristic number of sources ($N_{\star}$) as a free parameter, which is roughly the number of neutrino sources over which the measured intensity is divided. With existing two-year IceCube data, we can already constrain very rare, bright sources with $N_{\star}\lesssim$100. This can be improved to $N_{\star}\lesssim 10^4$-$10^5$ with IceCube-Gen2 and KM3NeT with ten-year exposure, constraining the contribution of BL Lacs ($N_{\star}=6\times10^{2}$). On the other hand, we can constrain weak sources with large number densities, like starburst galaxies ($N_{\star} = 10^{7}$), if we measure an anisotropic neutrino sky with future observations.

We report the detection of the pulsed signal of the radio-quiet magnetar-like pulsar PSR J1846-0258 in the high-energy \gr-ray data of the Fermi Large Area Telescope (Fermi LAT). We produced phase-coherent timing models exploiting RXTE PCA and Swift XRT monitoring data for the post- (magnetar-like) outburst period from 2007 August 28 to 2016 September 4, with independent verification using INTEGRAL ISGRI and Fermi GBM data. Phase-folding barycentric arrival times of selected Fermi LAT events from PSR J1846-0258, resulted in a 4.2 sigma detection (30--100 MeV) of a broad pulse consistent in shape and aligned in phase with the profiles that we measured with Swift XRT (2.5--10 keV), INTEGRAL ISGRI (20--150 keV) and Fermi GBM (20--300 keV). The pulsed flux (30--100 MeV) is (3.91 +/- 0.97)E-9 photons/(cm^2 s MeV). Declining significances of the INTEGRAL ISGRI 20--150 keV pulse profiles suggest fading of the pulsed hard X-ray emission during the post-outburst epochs. We revisited with greatly improved statistics the timing and spectral characteristics of PSR B1509-58 as measured with the Fermi LAT. The broad-band pulsed emission spectra (from 2 keV up to GeV energies) of PSR J1846-0258 and PSR B1509-58 can be accurately described with similarly curved shapes, with maximum luminosities at 3.5 +/- 1.1 MeV (PSR J1846-0258) and 2.23 +/- 0.11 MeV (PSR B1509-58). We discuss possible explanations for observational differences between Fermi LAT detected pulsars that reach maximum luminosities at GeV energies, like the second magnetar-like pulsar PSR J1119-6127, and pulsars with maximum luminosities at MeV energies, which might be due to geometric differences rather than exotic physics in high-B fields.