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.

Observational constraints on gamma rays produced by the annihilation of weakly interacting massive particles around primordial black holes (PBHs) imply that these two classes of Dark Matter candidates cannot coexist. We show here that the successful detection of one or more PBHs by radio searches (with the Square Kilometer Array) and gravitational waves searches (with LIGO/Virgo and the upcoming Einstein Telescope) would set extraordinarily stringent constraints on virtually all weak-scale extensions of the Standard Model with stable relics, including those predicting a WIMP abundance much smaller than that of Dark Matter. Upcoming PBHs searches have in particular the potential to rule out almost the entire parameter space of popular theories such as the minimal supersymmetric standard model and scalar singlet Dark Matter.

With several proposed MeV gamma-ray telescopes on the horizon, it is of paramount importance to perform accurate calculations of gamma-ray spectra expected from sub-GeV dark matter annihilation and decay. We present hazma, a python package for reliably computing these spectra, determining the resulting constraints from existing gamma-ray data, and prospects for upcoming telescopes. For high-level analyses, hazma comes with several built-in dark matter models where the interactions between dark matter and hadrons have been determined in detail using chiral perturbation theory. Additionally, hazma provides tools for computing spectra from individual final states with arbitrary numbers of light leptons and mesons, and for analyzing custom dark matter models. hazma can also produce electron and positron spectra from dark matter annihilation, enabling precise derivation of constraints from the cosmic microwave background.

Multiple space-borne cosmic ray detectors have detected line-like features in the electron and positron spectra. Most recently, the DAMPE collaboration reported the existence of such a feature at 1.4 TeV, sparking interest in a potential dark matter origin. Such quasi-monochromatic features, virtually free of any astrophysical background, could be explained by the annihilation of dark matter particles in a nearby dark matter clump. Here, we explore the consistency of producing such spectral features with dark matter annihilation from the standpoint of dark matter substructure statistics, constraints from anisotropy, and constraints from gamma-ray emission. We demonstrate that if indeed a high-energy, line-like feature in the electron-positron spectrum originates from dark matter annihilation in a nearby clump, a significant or even dominant fraction of the dark matter in the Solar System likely stems from the clump, with dramatic consequences for direct dark matter searches.

We demonstrate that the tentative detection of a few antihelium events with the Alpha Magnetic Spectrometer (AMS) on board the International Space Station can, in principle, be ascribed to the annihilation or decay of Galactic dark matter, when accounting for uncertainties in the coalescence process leading to the formation of antinuclei. We show that the predicted antiproton rate, assuming the antihelium events came from dark matter, is marginally consistent with AMS data, as is the antideuteron rate with current available constraints. We argue that a dark matter origin can be tested with better constraints on the coalescence process, better control of misidentified events, and with future antideuteron data.

We investigate the relation between the annihilation of dark matter (DM) particles into lepton pairs and into 2-body final states including one or two photons. We parametrize the DM interactions with leptons in terms of contact interactions, and calculate the loop-level annihilation into monochromatic gamma rays, specifically computing the ratio of the DM annihilation cross sections into two gamma rays versus lepton pairs. While the loop-level processes are generically suppressed in comparison with the tree-level annihilation into leptons, we find that some choices for the mediator spin and coupling structure lead to large branching fractions into gamma-ray lines. This result has implications for a dark matter contribution to the AMS-02 positron excess. We also explore the possibility of mediators which are charged under a dark symmetry and find that, for these loop-level processes, an effective field theory description is accurate for DM masses up to about half the mediator mass.

Cosmic-ray anti-nuclei provide a promising discovery channel for the indirect detection of particle dark matter. Hadron showers produced by the pair-annihilation or decay of Galactic dark matter generate anti-nucleons which can in turn form light anti-nuclei. Previous studies have only focused on the spectrum and flux of low energy antideuterons which, although very rarely, are occasionally also produced by cosmic-ray spallation. Heavier elements ($A\geq3$) have instead entirely negligible astrophysical background and a primary yield from dark matter which could be detectable by future experiments. Using a Monte Carlo event generator and an event-by-event phase space analysis, we compute, for the first time, the production spectrum of \antihe and \antiT for dark matter annihilating or decaying to $b\bar{b}$ and ${W^+}{W^-}$ final states. We then employ a semi-analytic model of interstellar and heliospheric propagation to calculate the \antihe flux as well as to provide tools to relate the anti-helium spectrum corresponding to an arbitrary antideuteron spectrum. Finally, we discuss prospects for current and future experiments, including GAPS and AMS-02.