Dark matter (DM) as a thermal relic of the primordial plasma is increasingly pressured by direct and indirect searches, while the same production mechanism in a decoupled sector is much less constrained. We extend the standard treatment of the freeze-out process to such scenarios and perform precision calculations of the annihilation cross section required to match the observed DM abundance. We demonstrate that the difference to the canonical value is generally sizeable, and can reach orders of magnitude. Our results directly impact the interpretation of DM searches in hidden sector scenarios.

Axion-like particles with masses in the keV-GeV range have a profound impact on the cosmological evolution of our Universe, in particular on the abundance of light elements produced during Big Bang Nucleosynthesis. The resulting limits are complementary to searches in the laboratory and provide valuable additional information regarding the validity of a given point in parameter space. A potential drawback is that altering the cosmological history may potentially weaken or even fully invalidate these bounds. The main objective of this article is therefore to evaluate the robustness of cosmological constraints on axion-like particles in the keV-GeV region, allowing for various additional effects which may weaken the bounds of the standard scenario. Employing the latest determinations of the primordial abundances as well as information from the cosmic microwave background we find that while bounds can indeed be weakened, very relevant robust constraints remain.

Traditional direct searches for dark matter, looking for nuclear recoils in deep underground detectors, are challenged by an almost complete loss of sensitivity for light dark matter particles. Consequently, there is a significant effort in the community to devise new methods and experiments to overcome these difficulties, constantly pushing the limits of the lowest dark matter mass that can be probed this way. From a model-building perspective, the scattering of sub-GeV dark matter on nucleons essentially must proceed via new light mediator particles, given that collider searches place extremely stringent bounds on contact-type interactions. Here we present an updated compilation of relevant limits for the case of a scalar mediator, including a new estimate of the near-future sensitivity of the NA62 experiment as well as a detailed evaluation of the model-specific limits from Big Bang nucleosynthesis. We also derive updated and more general limits on DM particles upscattered by cosmic rays, applicable to arbitrary energy- and momentum dependences of the scattering cross section. Finally we stress that dark matter self-interactions, when evaluated beyond the common s-wave approximation, place stringent limits independently of the dark matter production mechanism. These are, for the relevant parameter space, generically comparable to those that apply in the commonly studied freeze-out case. We conclude that the combination of existing (or expected) constraints from accelerators and astrophysics, combined with cosmological requirements, puts robust limits on the maximally possible nuclear scattering rate. In most regions of parameter space these are at least competitive with the best projected limits from currently planned direct detection experiments.

We present the first $N$-body simulations that adapt the equations of smoothed particle hydrodynamics to capture the effect of dark matter self-interactions which are too frequent to be resolved explicitly. The relevant energy transfer equations are derived, the appropriate thermal conductivity is determined and the effects of different smoothing kernels are studied. We apply our framework to simulate the formation of isothermal cores in isolated dark matter haloes and determine the core growth rate as a function of the self-scattering cross section. Our approach may be combined with explicit simulations of rare scatterings in order to simulate accurately the effects of arbitrary dark matter self-interactions in future cosmological simulations.

Thermal dark matter at the MeV scale faces stringent bounds from a variety of cosmological probes. Here we perform a detailed evaluation of BBN bounds on the annihilation cross section of dark matter with a mass $1\,\text{MeV} \lesssim m_\chi \lesssim 1\,\text{GeV}$. For $p-wave suppressed annihilations, constraints from BBN turn out to be significantly stronger than the ones from CMB observations, and are competitive with the strongest bounds from other indirect searches. We furthermore update the lower bound from BBN on the mass of thermal dark matter using improved determinations of primordial abundances. While being of similar strength as the corresponding bound from CMB, it is significantly more robust to changes in the particle physics model.

The idea of dark matter in the form of primordial black holes has seen a recent revival triggered by the LIGO detection of gravitational waves from binary black hole mergers. In this context, it has been argued that a large initial clustering of primordial black holes can help alleviate the strong constraints on this scenario. In this work, we show that on the contrary, with large initial clustering the problem is exacerbated and constraints on primordial black hole dark matter become overwhelmingly strong.

Meta-stable dark sector particles decaying into electrons or photons may non-trivially change the Hubble rate, lead to entropy injection into the thermal bath of Standard Model particles and may also photodisintegrate light nuclei formed in the early universe. We study generic constraints from Big Bang Nucleosynthesis on such a setup, with a particular emphasis on MeV-scale particles which are neither fully relativistic nor non-relativistic during all times relevant for Big Bang Nucleosynthesis. We apply our results to a simple model of self-interacting dark matter with a light scalar mediator. This setup turns out to be severely constrained by these considerations in combination with direct dark matter searches and will be fully tested with the next generation of low-threshold direct detection experiments.

Light vector mediators can naturally induce velocity-dependent dark matter self-interactions while at the same time allowing for the correct dark matter relic abundance via thermal freeze-out. If these mediators subsequently decay into Standard Model states such as electrons or photons however, this is robustly excluded by constraints from the Cosmic Microwave Background. We study to what extent this conclusion can be circumvented if the vector mediator is stable and hence contributes to the dark matter density while annihilating into lighter degrees of freedom. We find viable parts of parameter space which lead to the desired self-interaction cross section of dark matter to address the small-scale problems of the collisionless cold dark matter paradigm while being compatible with bounds from the Cosmic Microwave Background and Big Bang Nucleosynthesis observations.

Dark matter in the cosmological concordance model is parameterised by a single number, describing the covariantly conserved energy density of a non-relativistic fluid. Here we test this assumption in a model-independent and conservative way by considering the possibility that, at any point during the cosmological evolution, dark matter may be converted into a non-interacting form of radiation. This scenario encompasses, but is more general than, the cases where dark matter decays or annihilates into these states. We show that observations of the cosmic microwave background allow to strongly constrain this scenario for any conversion time after big bang nucleosynthesis. We discuss in detail, both from a Bayesian and frequentist point of view, in which sense adding large-scale structure observations may even provide a certain preference for a conversion of dark matter to radiation at late times. Finally we apply our general results to a specific particle physics realisation of such a scenario, featuring late kinetic decoupling and Sommerfeld-enhanced dark matter annihilation. We identify a small part of parameter space that both mitigates the tension between cosmic microwave and large-scale structure data and allows for velocity-dependent dark matter self-interactions strong enough to address the small-scale problems of structure formation.

The relative sensitivity of different direct detection experiments depends sensitively on the astrophysical distribution and particle physics nature of dark matter, prohibiting a model-independent comparison. The situation changes fundamentally if two experiments employ the same target material. We show that in this case one can compare measurements of an annual modulation and exclusion bounds on the total rate while making no assumptions on astrophysics and no (or only very general) assumptions on particle physics. In particular, we show that the dark matter interpretation of the DAMA/LIBRA signal can be conclusively tested with COSINUS, a future experiment employing the same target material. We find that if COSINUS excludes a dark matter scattering rate of about $0.01\,\text{kg}^{-1}\,\text{days}^{-1}$ with an energy threshold of $1.8\,$keV and resolution of $0.2\,$keV, it will rule out all explanations of DAMA/LIBRA in terms of dark matter scattering off sodium and/or iodine.

We study constraints from Big Bang Nucleosynthesis on inert particles in a dark sector which contribute to the Hubble rate and therefore change the predictions of the primordial nuclear abundances. We pay special attention to the case of MeV-scale particles decaying into dark radiation, which are neither fully relativistic nor non-relativistic during all temperatures relevant to Big Bang Nucleosynthesis. As an application we discuss the implications of our general results for models of self-interacting dark matter with light mediators.

Self-interacting dark matter may have striking astrophysical signatures, such as observable offsets between galaxies and dark matter in merging galaxy clusters. Numerical N-body simulations used to predict such observables typically treat the galaxies as collisionless test particles, a questionable assumption given that each galaxy is embedded in its own dark matter halo. To enable a more accurate treatment we develop an effective description of small dark matter haloes taking into account the two major effects due to dark matter self-scatterings: deceleration and evaporation. We point out that self-scatterings can have a sizeable impact on the trajectories of galaxies, diminishing the separation between galaxies and dark matter in merging clusters. This effect depends sensitively on the underlying particle physics, in particular the angular dependence of the self-scattering cross section, and cannot be predicted from the momentum transfer cross section alone.

Dark matter particles interacting via the exchange of very light spin-0 mediators can have large self-interaction rates and obtain their relic abundance from thermal freeze-out. At the same time, these models face strong bounds from direct and indirect probes of dark matter as well as a number of constraints on the properties of the mediator. We investigate whether these constraints can be consistent with having observable effects from dark matter self-interactions in astrophysical systems. For the case of a mediator with purely scalar couplings we point out the highly relevant impact of low-threshold direct detection experiments like CRESST-II, which essentially rule out the simplest realization of this model. These constraints can be significantly relaxed if the mediator has CP-violating couplings, but then the model faces strong constraints from CMB measurements, which can only be avoided in special regions of parameter space.

A reliable comparison of different dark matter (DM) searches requires models that satisfy certain consistency requirements like gauge invariance and perturbative unitarity. As a well-motivated example, we study two-mediator DM (2MDM). The model is based on a spontaneously broken $U(1)'$ gauge symmetry and contains a Majorana DM particle as well as two $s$-channel mediators, one vector (the $Z'$) and one scalar (the dark Higgs). We perform a global scan over the parameters of the model assuming that the DM relic density is obtained by thermal freeze-out in the early Universe and imposing a large set of constraints: direct and indirect DM searches, monojet, dijet and dilepton searches at colliders, Higgs observables, electroweak precision tests and perturbative unitarity. We conclude that thermal DM is only allowed either close to an $s$-channel resonance or if at least one mediator is lighter than the DM particle. In these cases a thermal DM abundance can be obtained although DM couplings to the Standard Model are tiny. Interestingly, we find that vector-mediated DM-nucleon scattering leads to relevant constraints despite the velocity-suppressed cross section, and that indirect detection can be important if DM annihilations into both mediators are kinematically allowed.

We show that simplified models used to describe the interactions of dark matter with Standard Model particles do not in general respect gauge invariance and that perturbative unitarity may be violated in large regions of the parameter space. The modifications necessary to cure these inconsistencies may imply a much richer phenomenology and lead to stringent constraints on the model. We illustrate these observations by considering the simplified model of a fermionic dark matter particle and a vector mediator. Imposing gauge invariance then leads to strong constraints from dilepton resonance searches and electroweak precision tests. Furthermore, the new states required to restore perturbative unitarity can mix with Standard Model states and mediate interactions between the dark and the visible sector, leading to new experimental signatures such as invisible Higgs decays. The resulting constraints are typically stronger than the 'classic' constraints on DM simplified models such as monojet searches and make it difficult to avoid thermal overproduction of dark matter.

Self-interactions of dark matter particles can potentially lead to an observable separation between the dark matter halo and the stars of a galaxy moving through a region of large dark matter density. Such a separation has recently been observed in a galaxy falling into the core of the galaxy cluster Abell 3827. We estimated the DM self-interaction cross section needed to reproduce the observed effects and find that the sensitivity of Abell 3827 has been significantly overestimated in a previous study. Our corrected estimate is $\tilde{\sigma}/m_\text{DM} \sim 3\:\text{cm}^2\:\text{g}^{-1}$ when self-interactions result in an effective drag force and $\sigma/m_\text{DM} \sim 1.5\:\text{cm}^2\:\text{g}^{-1}$ for the case of contact interactions, in some tension with previous upper bounds.

When a dark matter halo moves through a background of dark matter particles, self-interactions can lead to both deceleration and evaporation of the halo and thus shift its centroid relative to the collisionless stars and galaxies. We study the magnitude and time evolution of this shift for two classes of dark matter self-interactions, viz. frequent self-interactions with small momentum transfer (e.g. due to long-range interactions) and rare self-interactions with large momentum transfer (e.g. contact interactions), and find important differences between the two cases. We find that neither effect can be strong enough to completely separate the dark matter halo from the galaxies, if we impose conservative bounds on the self-interaction cross-section. The majority of both populations remain bound to the same gravitational potential and the peaks of their distributions are therefore always coincident. Consequently any apparent separation is mainly due to particles which are leaving the gravitational potential, so will be largest shortly after the collision but not observable in evolved systems. Nevertheless the fraction of collisions with large momentum transfer is an important characteristic of self-interactions, which can potentially be extracted from observational data and provide an important clue as to the nature of dark matter.

The CDMS-II collaboration has reported 3 events in a Si detector, which are consistent with being nuclear recoils due to scattering of Galactic dark matter particles with a mass of about 8.6 GeV and a cross-section on neutrons of about 2 x 10^-41 cm^2. While a previous result from the XENON10 experiment has supposedly ruled out such particles as dark matter, we find by reanalysing the XENON10 data that this is not the case. Some tension remains however with the upper limit placed by the XENON100 experiment, independently of astrophysical uncertainties concerning the Galactic dark matter distribution. We explore possible ways of ameliorating this tension by altering the properties of dark matter interactions. Nevertheless, even with standard couplings, light dark matter is consistent with both CDMS and XENON10/100.

Improved limits as well as tentative claims for dark matter annihilation into gamma-ray lines have been presented recently. We study the direct detection cross section induced from dark matter annihilation into two photons in a model-independent fashion, assuming no additional couplings between dark matter and nuclei. We find a striking non-standard recoil spectrum due to different destructively interfering contributions to the dark matter nucleus scattering cross section. While in the case of s-wave annihilation the current sensitivity of direct detection experiments is insufficient to compete with indirect detection searches, for p-wave annihilation the constraints from direct searches are comparable. This will allow to test dark matter scenarios with p-wave annihilation that predict a large di-photon annihilation cross section in the next generation of experiments.

We study the impact of the assumed velocity distribution of galactic dark matter particles on the interpretation of results from nuclear recoil detectors. By converting experimental data to variables that make the astrophysical unknowns explicit, different experiments can be compared without implicit assumptions concerning the dark matter halo. We extend this framework to include the annual modulation signal, as well as multiple target elements. Recent results from DAMA, CoGeNT and CRESST-II can be brought into agreement if the velocity distribution is very anisotropic and thus allows a large modulation fraction. However constraints from CDMS and XENON cannot be evaded by appealing to such astrophysical uncertainties alone.

DAMA observes an annual modulation in their event rate, as might be expected from dark matter scatterings, while CoGeNT has reported evidence for a similar modulation. The simplest interpretation of these findings in terms of dark matter-nucleus scatterings is excluded by other direct detection experiments. We consider the robustness of these exclusions with respect to assumptions regarding the scattering and find that isospin-violating inelastic dark matter helps alleviate this tension and allows marginal compatibility between experiments. Isospin-violation can significantly weaken the XENON constraints, while inelasticity enhances the annual modulation fraction of the signal, bringing the CoGeNT and CDMS results into better agreement.

A ~5 GeV `dark baryon' with a cosmic asymmetry similar to that of baryons is a natural candidate for the dark matter. We study the possibility of generating such a state through dynamical electroweak symmetry breaking, and show that it can share the relic baryon asymmetry via sphaleron interactions, even though it has no electroweak interactions. The scattering cross-section on nucleons, estimated in analogy to QCD, is within reach of underground direct detection experiments.

We perform an extensive study of the DAMA annual modulation data in the context of inelastic dark matter. We find that inelastic dark matter with mass m > 15 GeV is excluded at the 95% confidence level by the combination of DAMA spectral information and results from other direct detection experiments. However, at smaller m, inelastic dark matter constitutes a possible solution to the DAMA puzzle.