Since upcoming telescopes will observe thousands of strong lensing systems, creating fully-automated analysis pipelines for these images becomes increasingly important. In this work, we make a step towards that direction by developing the first end-to-end differentiable strong lensing pipeline. Our approach leverages and combines three important computer science developments: (a) convolutional neural networks, (b) efficient gradient-based sampling techniques, and (c) deep probabilistic programming languages. The latter automatize parameter inference and enable the combination of generative deep neural networks and physics components in a single model. In the current work, we demonstrate that it is possible to combine a convolutional neural network trained on galaxy images as a source model with a fully-differentiable and exact implementation of gravitational lensing physics in a single probabilistic model. This does away with hyperparameter tuning for the source model, enables the simultaneous optimization of nearly one hundred source and lens parameters with gradient-based methods, and allows the use of efficient gradient-based posterior sampling techniques. These features make this automated inference pipeline potentially suitable for processing a large amount of data. By analyzing mock lensing systems with different signal-to-noise ratios, we show that lensing parameters are reconstructed with percent-level accuracy. More generally, we consider this work as one of the first steps in establishing differentiable probabilistic programming techniques in the particle astrophysics community, which have the potential to significantly accelerate and improve many complex data analysis tasks.

We propose a multi-messenger probe of QCD axion Dark Matter based on observations of black hole-neutron star binary inspirals. It is suggested that a dense Dark Matter spike may grow around intermediate mass black holes ($10^{3}-10^{5} \mathrm{\,M_{\odot}}$). The presence of such a spike produces two unique effects: a distinct phase shift in the gravitational wave strain during the inspiral and an enhancement of the radio emission due to the resonant axion-photon conversion occurring in the neutron star magnetosphere throughout the inspiral and merger. Remarkably, the observation of the gravitational wave signal can be used to infer the Dark Matter density and, consequently, to predict the radio emission. We study the projected reach of the LISA interferometer and next-generation radio telescopes such as the Square Kilometre Array. Given a sufficiently nearby system, such observations will potentially allow for the detection of QCD axion Dark Matter in the mass range $10^{-7}\,\mathrm{eV}$ to $10^{-5}\,\mathrm{eV}$.

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.

The origin of neutrino flux observed in IceCube is still mainly unknown. Typically two flux components are assumed, namely: atmospheric neutrinos and an unknown astrophysical term. In principle the latter could also contain a top-down contribution coming for example from decaying dark matter. In this case one should also expect prompt and secondary gammas as well. This leads to the possibility of a multimessenger analysis based on the simultaneous comparison of the Dark Matter hypothesis both with neutrino and high energy gamma rays data. In this paper, we analyze, for different decaying Dark Matter channels, the 7.5 years IceCube HESE data, and compare the results with previous exclusion limits coming from Fermi data. Finally, we test whether the Dark Matter hypothesis could be further scrutinised by using forthcoming high energy gamma rays experiments.

The difficult task of observing Dark Matter subhaloes is of paramount importance since it would constrain Dark Matter particle properties (cold or warm relic) and confirm once again the longstanding $\Lambda$CDM model. In the near future the new generation of ground and space surveys will observe thousands of strong gravitational lensing systems providing a unique probe of Dark Matter substructures. Here, we describe a new strong lensing analysis pipeline that combines deep Convolutional Neural Networks with physical models and exploits traditional sampling techniques such as Hamiltonian Monte Carlo. Using simulated strong gravitational lensing systems, we discuss first results and characterize the accuracy of the reconstruction of the main lensing parameters.

We show how 21 cm cosmology can test relic neutrino radiative decays into sterile neutrinos. Using recent EDGES results, we derive constraints on the lifetime of the decaying neutrinos. If the EDGES anomaly will be confirmed, then there are two solutions, one for much longer and one for much shorter lifetimes than the age of the universe, showing how relic neutrino radiative decays can explain the anomaly in a simple way. We also show how to combine EDGES results with those from radio background observations, showing that potentially the ARCADE 2 excess can be also reproduced together with the EDGES anomaly within the proposed non-standard cosmological scenario. Our calculation of the specific intensity at the redshifts probed by EDGES can be also applied to the case of decaying dark matter and it also corrects a flawed expression appeared in previous literature.

We perform the first statistical combined analysis of the diffuse neutrino flux observed by ANTARES (nine-year) and IceCube (six-year) by assuming a single astrophysical power-law flux. The combined analysis reduces by a few percent the best-fit values for the flux normalization and the spectral index. Both data samples show an excess in the same energy range (40--200 TeV), suggesting the presence of a second component. We perform a goodness-of-fit test to scrutinize the null assumption of a single power-law, scanning different values for the spectral index. The addition of the ANTARES data reduces the $p$-value by a factor 2$\div$3. In particular, a single power-law component in the neutrino flux with the spectral index deduced by the six-year up-going muon neutrinos of IceCube is disfavored with a $p$-value smaller than $10^{-2}$.

The assumption of a single astrophysical power-law flux to explain the IceCube 6-year HESE extraterrestrial events yields a large spectral index that is in tension with gamma-ray observations and the 6-year up-going muon neutrinos data. Adopting a spectral index belonging to the range $\left[2.0,2.2\right]$, which is compatible with the one deduced by the analysis performed on the 6-year up-going muon neutrinos data and with $p$-$p$ astrophysical sources, the latest IceCube data show an up to $2.6\,\sigma$ excess in the number of events in the energy range 40--200 TeV. We interpret such an excess as a decaying Dark Matter signal and we perform a likelihood-ratio statistical test to compare the two-component scenario with respect to the single-component one.

The recent study on the the 6-year up-going muon neutrinos by the IceCube Collaboration support the hypothesis of a two-component scenario explaining the diffuse TeV-PeV neutrino flux. Once a hard astrophysical power-law is considered, an excess in the IceCube data is shown in the energy range 10-100 TeV (low-energy excess}). By means of a statistical analysis on the neutrino energy spectrum and on the angular distribution of neutrino arrival directions, we characterize a two-component neutrino flux where decaying/annihilating Dark Matter particles provide a contribution to the IceCube observations.

The high energy events observed at the IceCube Neutrino Observatory have triggered many investigations interpreting the highly energetic neutrinos detected as decay products of heavy unstable Dark Matter particles. However, while very detailed treatments of the IceCube phenomenology exist, only a few references focus on the (non-trivial) Dark Matter production part -- and all of those rely on relatively complicated new models which are not always testable directly. We instead investigate two of the most minimal scenarios possible, where the operator responsible for the IceCube events is directly involved in Dark Matter production. We show that the simplest (four-dimensional) operator is not powerful enough to accommodate all constraints. A more non-minimal setting (at mass dimension six), however, can do both fitting all the data and also allowing for a comparatively small parameter space only, parts of which can be in reach of future observations. We conclude that minimalistic approaches can be enough to explain all data required, while complicated new physics seems not to be required by IceCube.

The recent study on the the 6-year up-going muon neutrinos by the IceCube Collaboration and the multi-messenger analyses support the hypothesis of a two-component scenario explaining the diffuse TeV-PeV neutrino flux. Depending on the steepness of the astrophysical power-law, an excess in the IceCube data is shown in the energy range 10-100 TeV (low-energy excess) or at PeV (high-energy excess). In both cases, we characterize a two-component neutrino flux where decaying Dark Matter particles provide a contribution to the IceCube observations.

The 2-years MESE IceCube events show a slightly excess in the energy range 10-100 TeV with a maximum local statistical significance of 2.3$\sigma$, once a hard astrophysical power-law is assumed. A spectral index smaller than 2.2 is indeed suggested by multi-messenger studies related to $p$-$p$ sources and by the recent IceCube analysis regarding 6-years up-going muon neutrinos. In the present paper, we propose a two-components scenario where the extraterrestrial neutrinos are explained in terms of an astrophysical power-law and a Dark Matter signal. We consider both decaying and annihilating Dark Matter candidates with different final states (quarks and leptons) and different halo density profiles. We perform a likelihood-ratio analysis that provides a statistical significance up to 3.9$\sigma$ for a Dark Matter interpretation of the IceCube low energy excess.

We analyze the scenario where the IceCube high energy neutrino events are explained in terms of an extraterrestrial flux due to two different components: a contribution coming from know astrophysical sources for energies up to few hundreds TeV and a top-down contribution originated by the decay of heavy dark matter particles with a mass of few PeV. Contrary to previous approaches, we consider a leptophilic three-body decay that dominates at PeV energies due to the absence of quarks in the final state. We find that the theoretical predictions of such a scenario are in a slightly better agreement with the IceCube data if the astrophysical component has a cut-off at about 100 TeV. This interpretation of IceCube data can be easily tested in the near future since the decaying dark matter scenario predicts a sharp cut-off at PeV energy scale and the observation of an anisotropy towards Galactic Center of our Galaxy in contrast with the isotropic astrophysical flux.

In this Letter we focus our attention on the IceCube events in the energy range between 60 and 100 TeV, which show an order 2-sigma excess with respect to a power-law with spectral index 2. We analyze the possible origin of such an excess by comparing the distribution of the arrival directions of IceCube events with the angular distributions of simply distributed astrophysical galactic/extragalactic sources, as well as with the expected flux coming from DM interactions (decay and annihilation) for different DM profiles. The statistical analysis performed seems to disfavor the correlation with the galactic plane, whereas rules out the DM annihilation scenario only in case of small clumpiness effect. The small statistics till now collected does not allow to scrutinize the cases of astrophysical isotropic distribution and DM decay scenarios. For this reason we perform a forecast analysis in order to stress the role of future Neutrino Telescopes.