Several studies have pointed out an excess in the AMS-02 antiproton spectrum at rigidities of 10-20 GV. Its spectral properties were found to be consistent with a dark-matter particle of mass 50-100 GeV which annihilates hadronically at roughly the thermal rate. In this work, we reinvestigate the antiproton excess including all relevant sources of systematic errors. Most importantly, we perform a fully realistic derivation of the correlations in the AMS-02 systematic error which could potentially "fake" a dark-matter signal. The dominant systematics in the relevant rigidity range originate from uncertainties in the cross sections for absorption of cosmic rays within the detector material. For the first time ever, we calculate their correlations within the full Glauber-Gribov theory of inelastic scattering. The AMS-02 correlations enter our spectral search for dark matter in the form of covariance matrices which we make publicly available for the cosmic-ray community. We find that the global significance of the antiproton excess is reduced to below 1 $\sigma$ once all systematics, including the derived AMS-02 error correlations, are taken into account. No significant preference for a dark-matter signal in the AMS-02 antiproton data is found in the mass range 10-10000 GeV.

We perform a global fit within the pseudo-Nambu-Goldstone dark matter (DM) model emerging from an additional complex scalar singlet with a softly broken global U(1) symmetry. Leading to a momentum-suppressed DM-nucleon cross section at tree level, the model provides a natural explanation for the null results from direct detection experiments. Our global fit combines constraints from perturbative unitarity, DM relic abundance, Higgs invisible decay, electroweak precision observables and latest Higgs searches at colliders. The results are presented in both frequentist and Bayesian statistical frameworks. Furthermore, post-processing our samples, we include the likelihood from gamma-ray observations of $\mathit{Fermi}$-LAT dwarf spheroidal galaxies and compute the one-loop DM-nucleon cross section. We find two favoured regions characterised by their dominant annihilation channel: the Higgs funnel and annihilation into Higgs pairs. Both are compatible with current $\mathit{Fermi}$-LAT observations, and furthermore, can fit the slight excess observed in four dwarfs in a mass range between about 30-300 GeV. While the former region is hard to probe experimentally, the latter can partly be tested by current observations of cosmic-ray antiprotons as well as future gamma-ray observations.

Global fits of primary and secondary cosmic-ray (CR) fluxes measured by AMS-02 have great potential to study CR propagation models and search for exotic sources of antimatter such as annihilating dark matter (DM). Previous studies of AMS-02 antiprotons revealed a possible hint for a DM signal which, however, could be affected by systematic uncertainties. To test the robustness of such a DM signal, in this work we systematically study two important sources of uncertainties: the antiproton production cross sections needed to calculate the source spectra of secondary antiprotons and the potential correlations in the experimental data, so far not provided by the AMS-02 Collaboration. To investigate the impact of cross-section uncertainties we perform global fits of CR spectra including a covariance matrix determined from nuclear cross-section measurements. As an alternative approach, we perform a joint fit to both the CR and cross-section data. The two methods agree and show that cross-section uncertainties have a small effect on the CR fits and on the significance of a potential DM signal, which we find to be at the level of $3\sigma$. Correlations in the data can have a much larger impact. To illustrate this effect, we determine possible benchmark models for the correlations in a data-driven method. The inclusion of correlations strongly improves the constraints on the propagation model and, furthermore, enhances the significance of the DM signal up to above $5\sigma$. Our analysis demonstrates the importance of providing the covariance of the experimental data, which is needed to fully exploit their potential.

We consider dark matter (DM) with very weak couplings to the standard model (SM), such that its self-annihilation cross section is much smaller than the canonical one, $\langle\sigma v\rangle_{\chi\chi} \ll 10^{-26}\mathrm{cm}^3/\mathrm{s}$. In this case DM self-annihilation is negligible for the dynamics of freeze-out and DM dilution is solely driven by efficient annihilation of heavier accompanying dark sector particles provided that DM maintains chemical equilibrium with the dark sector. This chemical equilibrium is established by conversion processes which require much smaller couplings to be efficient than annihilation. The chemical decoupling of DM from the SM can either be initiated by the freeze-out of annihilation, resembling a co-annihilation scenario, or of conversion processes, leading to the scenario of conversion-driven freeze-out. We focus on the latter and discuss its distinct phenomenology.

We present MadDM v.3.0, a numerical tool to compute particle dark matter observables in generic new physics models. The new version features a comprehensive and automated framework for dark matter searches at the interface of collider physics, astrophysics and cosmology and is deployed as a plugin of the MadGraph5_aMC@NLO platform, inheriting most of its features. With respect to the previous version, MadDM v.3.0 can now provide predictions for indirect dark matter signatures in astrophysical environments, such as the annihilation cross section at present time and the energy spectra of prompt photons, cosmic rays and neutrinos resulting from dark matter annihilation. MadDM indirect detection features support both $2\to2$ and $2 \to n$ dark matter annihilation processes. In addition, the ability to compare theoretical predictions with experimental constraints is extended by including the Fermi-LAT likelihood for gamma-ray constraints from dwarf spheroidal galaxies and by providing an interface with the nested sampling algorithm PyMultinNest to perform high dimensional parameter scans efficiently. We validate the code for a wide set of dark matter models by comparing the results from MadDM v.3.0 to existing tools and results in the literature.

Non-thermalized dark matter is a cosmologically valid alternative to the paradigm of weakly interacting massive particles. For dark matter belonging to a $Z_2$-odd sector that contains in addition a thermalized mediator particle, dark matter production proceeds in general via both the freeze-in and superWIMP mechanism. We highlight their interplay and emphasize the connection to long-lived particles at colliders. For the explicit example of a colored t-channel mediator model we map out the entire accessible parameter space, cornered by bounds from the LHC, big bang nucleosynthesis and Lyman-alpha forest observations, respectively. We discuss prospects for the HL- and HE-LHC.

We present a comprehensive analysis of top-philic Majorana dark matter that interacts via a colored t-channel mediator. Despite the simplicity of the model -- introducing three parameters only -- it provides an extremely rich phenomenology allowing us to accommodate the relic density for a large range of coupling strengths spanning over six orders of magnitude. This model features all `exceptional' mechanisms for dark matter freeze-out, including the recently discovered conversion-driven freeze-out mode, with interesting signatures of long-lived colored particles at colliders. We constrain the cosmologically allowed parameter space with current experimental limits from direct, indirect and collider searches, with special emphasis on light dark matter below the top mass. In particular, we explore the interplay between limits from Xenon1T, Fermi-LAT and AMS-02 as well as limits from stop, monojet and Higgs invisible decay searches at the LHC. We find that several blind spots for light dark matter evade current constraints. The region in parameter space where the relic density is set by the mechanism of conversion-driven freeze-out can be conclusively tested by R-hadron searches at the LHC with 300\,fb$^{-1}$.

Cosmic-ray observations provide a powerful probe of dark matter annihilation in the Galaxy. In this paper we derive constraints on heavy dark matter from the recent precise AMS-02 antiproton data. We consider all possible annihilation channels into pairs of standard model particles. Furthermore, we interpret our results in the context of minimal dark matter, including higgsino, wino and quintuplet dark matter. We compare the cosmic-ray antiproton limits to limits from $\gamma$-ray observations of dwarf spheroidal galaxies and to limits from $\gamma$-ray and $\gamma$-line observations towards the Galactic center. While the latter limits are highly dependent on the dark matter density distribution and only exclude a thermal wino for cuspy profiles, the cosmic-ray limits are more robust, strongly disfavoring the thermal wino dark matter scenario even for a conservative estimate of systematic uncertainties.

Observations of cosmic rays are a sensitive probe of dark matter annihilation in our Galaxy. In this article we present an analysis of the AMS-02 antiproton data, reducing cosmic-ray propagation uncertainties by fitting at the same time dark matter and propagation parameters. The result exhibits a possible hint for dark matter pointing to an annihilation cross section close to the thermal value. We investigate the compatibility of this signal with a dark matter interpretation of the Galactic center excess seen in the Fermi-LAT gamma-ray data and discuss implications for dark matter models.

Chemical equilibrium is a commonly made assumption in the freeze-out calculation of coannihilating dark matter. We explore the possible failure of this assumption and find a new conversion-driven freeze-out mechanism. Considering a representative simplified model inspired by supersymmetry with a neutralino- and sbottom-like particle we find regions in parameter space with very small couplings accommodating the measured relic density. In this region freeze-out takes place out of chemical equilibrium and dark matter self-annihilation is thoroughly inefficient. The relic density is governed primarily by the size of the conversion terms in the Boltzmann equations. Due to the small dark matter coupling the parameter region is immune to direct detection but predicts an interesting signature of disappearing tracks or displaced vertices at the LHC. Unlike freeze-in or superWIMP scenarios, conversion-driven freeze-out is not sensitive to the initial conditions at the end of reheating.

A possible hint of dark matter annihilation has been found in Cuoco, Korsmeier and Kr\"amer (2017) from an analysis of recent cosmic-ray antiproton data from AMS-02 and taking into account cosmic-ray propagation uncertainties by fitting at the same time dark matter and propagation parameters. Here, we extend this analysis to a wider class of annihilation channels. We find consistent hints of a dark matter signal with an annihilation cross-section close to the thermal value and with masses in range between 40 and 130 GeV depending on the annihilation channel. Furthermore, we investigate in how far the possible signal is compatible with the Galactic center gamma-ray excess and recent observation of dwarf satellite galaxies by performing a joint global fit including uncertainties in the dark matter density profile. As an example, we interpret our results in the framework of the Higgs portal model.

We perform a global fit within the inert doublet model taking into account experimental observables from colliders, direct and indirect dark matter searches and theoretical constraints. In particular, we consider recent results from searches for dark matter annihilation-induced gamma-rays in dwarf spheroidal galaxies and relax the assumption that the inert doublet model should account for the entire dark matter in the Universe. We moreover study in how far the model is compatible with a possible dark matter explanation of the so-called Galactic center excess. We find two distinct parameter space regions that are consistent with existing constraints and can simultaneously explain the excess: One with dark matter masses near the Higgs resonance and one around 72 GeV where dark matter annihilates predominantly into pairs of virtual electroweak gauge bosons via the four-vertex arising from the inert doublet's kinetic term. We briefly discuss future prospects to probe these scenarios.

We show that solar $\gamma$-ray observations can provide a complementary probe of Dark Matter in scenarios where the interactions with the Standard Model proceed via long-lived mediators. For illustration we consider a simplified model which provides solar $\gamma$-ray fluxes observable with the next generation $\gamma$-ray telescopes, while complying with the existing experimental constraints. Our results suggest that solar $\gamma$-ray fluxes can be orders of magnitude larger than the ones from the Galactic center, while being subject to low backgrounds.

We show that observations of solar $\gamma$-rays offer a novel probe of dark matter in scenarios where interactions with the visible sector proceed via a long-lived mediator. As a proof of principle, we demonstrate that there exists a class of models which yield solar $\gamma$-ray fluxes observable with the next generation of $\gamma$-ray telescopes, while being allowed by a variety of current experimental constraints. The parameter space allowed by big bang nucleosynthesis and beam dump experiments naturally leads to mediator lifetimes sufficient to produce observable solar $\gamma$-ray signals. The model allows for solar $\gamma$-ray fluxes up to orders of magnitude larger compared to dwarf spheroidal galaxies, without reaching equilibrium between dark matter annihilation and capture rate. Our results suggest that solar $\gamma$-ray observations are complementary, and in some cases superior, to existing and future dark matter detection efforts.

We analyse the excess in the $\gamma$-ray emission from the center of our galaxy observed by Fermi-LAT in terms of dark matter annihilation within the scalar Higgs portal model. In particular, we include the astrophysical uncertainties from the dark matter distribution and allow for unspecified additional dark matter components. We demonstrate through a detailed numerical fit that the strength and shape of the $\gamma$-ray spectrum can indeed be described by the model in various regions of dark matter masses and couplings. Constraints from invisible Higgs decays, direct dark matter searches, indirect searches in dwarf galaxies and for $\gamma$-ray lines, and constraints from the dark matter relic density reduce the parameter space to dark matter masses near the Higgs resonance. We find two viable regions: one where the Higgs-dark matter coupling is of ${\cal O}(10^{-2})$, and an additional dark matter component beyond the scalar WIMP of our model is preferred, and one region where the Higgs-dark matter coupling may be significantly smaller, but where the scalar WIMP constitutes a significant fraction or even all of dark matter. Both viable regions are hard to probe in future direct detection and collider experiments.

We consider a simplified model for Majorana fermion dark matter and explore constraints from direct, indirect and LHC collider searches. The dark matter is assumed to couple to the Standard Model through a vector mediator with axial-vector interactions. We provide detailed analyses of LHC mono-jet searches and IceCube limits on dark matter annihilation in the Sun. In particular, we develop a method for calculating limits on simplified WIMP dark matter models from public IceCube data, which are only available for a limited number of Standard Model final states. We demonstrate that LHC and IceCube searches for Majorana dark matter are complementary and derive new limits on the dark matter and mediator masses, including in addition constraints from LHC di-jet searches, direct detection and the dark matter relic density.

Despite its importance for modeling the homogeneous hot early universe very little is experimentally known about the magnitude of the reheating temperature, leaving an uncertainty of remarkable 18 orders of magnitude. In this paper we consider a general class of polynomial inflaton potentials up to fourth order. Employing a Monte Carlo scan and imposing theoretical and experimental constraints we derive a robust lower limit on the energy scale at the end of inflation, $V_\text{end}^{1/4} > 3 \times 10^{15}$ GeV for sizable tensor modes, $r \geq 10^{-3}$. If the reheating phase is perturbative and matter dominated, this translates into a lower bound on the reheating temperature, yielding $T_\text{rh} > 3 \times 10^8 \; (7 \times 10^2)$ GeV for gravitational inflaton decay through a generic dimension five (six) operator.

The presence of high reheating temperatures in the thermal history of the universe challenges supersymmetric scenarios owing to the gravitino problem. We revise a general R-parity conserving gravitino dark matter scenario with a stau as the next-to-lightest superparticle being particularly constrained by searches for heavy stable charged particles at the LHC. Imposing a variety of experimental and theoretical constraints we show that points with T_R > 10^9 GeV survive only in a very particular corner of parameter space.