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.

Warm dark matter cosmologies have been widely studied as an alternative to the cold dark matter paradigm, the characteristic feature being a suppression of structure formation on small cosmological scales. A very similar situation occurs if standard cold dark matter particles are kept in local thermal equilibrium with a, possibly dark, relativistic species until the universe has cooled down to keV temperatures. We perform a systematic phenomenological study of this possibility, and classify all minimal models containing dark matter and an arbitrary radiation component that allow such a late kinetic decoupling. We recover explicit cases recently discussed in the literature and identify new classes of examples that are very interesting from a model-building point of view. In some of these models dark matter is inevitably self-interacting, which is remarkable in view of recent observational support for this possibility. Hence, dark matter models featuring late kinetic decoupling have the potential not only to alleviate the missing satellites problem but also to address other problems of the cosmological concordance model on small scales, in particular the cusp-core and too-big-too-fail problems, in some cases without invoking any additional input.

The annihilation of non-relativistic dark matter particles at tree level can be strongly enhanced by the radiation of an additional gauge boson. This is particularly true for the helicity-suppressed annihilation of Majorana particles, like neutralinos, into fermion pairs. Surprisingly, and despite the potentially large effect due to the strong coupling, this has so far been studied in much less detail for the internal bremsstrahlung of gluons than for photons or electroweak gauge bosons. Here, we aim at bridging that gap by presenting a general analysis of neutralino annihilation into quark anti-quark pairs and a gluon, allowing e.g. for arbitrary neutralino compositions and keeping the leading quark mass dependence at all stages in the calculation. We find in some cases largely enhanced annihilation rates, especially for scenarios with squarks being close to degenerate in mass with the lightest neutralino, but also notable distortions in the associated antiproton and gamma-ray spectra. Both effects significantly impact limits from indirect searches for dark matter and are thus important to be taken into account in, e.g., global scans. For extensive scans, on the other hand, full calculations of QCD corrections are numerically typically too expensive to perform for each point in parameter space. We present here for the first time an efficient, numerically fast implementation of QCD corrections, extendable in a straight-forward way to non-supersymmetric models, which avoids computationally demanding full one-loop calculations or event generator runs and yet fully captures the leading effects relevant for indirect dark matter searches. In this context, we also present updated constraints on dark matter annihilation from cosmic-ray antiproton data. Finally, we comment on the impact of our results on relic density calculations.

We use WMAP 9-year and other CMB data to constrain cosmological models where the primordial perturbations have both an adiabatic and a (possibly correlated) neutrino density (NDI), neutrino velocity (NVI), or cold dark matter density (CDI) isocurvature component. For NDI and CDI we use both a phenomenological approach, where primordial perturbations are parametrized in terms of amplitudes at two scales, and a slow-roll two-field inflation approach, where slow-roll parameters are used as primary parameters. For NVI we use only the phenomenological approach, since it is difficult to imagine a connection with inflation. We find that in the NDI and NVI cases larger isocurvature fractions are allowed than in the corresponding models with CDI. For uncorrelated perturbations, the upper limit to the primordial NDI (NVI) fraction is 24% (20%) at k = 0.002 Mpc^{-1} and 28% (16%) at k = 0.01 Mpc^{-1}. For maximally correlated (anticorrelated) perturbations, the upper limit to the NDI fraction is 3.0% (0.9%). The nonadiabatic contribution to the CMB temperature variance can be as large as 10% (-13%) for the NDI (NVI) modes. Bayesian model comparison favors pure adiabatic initial mode over the mixed primordial adiabatic and NDI, NVI, or CDI perturbations. At best, the betting odds for a mixed model (uncorrelated NDI) are 1:3.4 compared to the pure adiabatic model. For the phenomenological generally correlated mixed models the odds are about 1:100, whereas the slow-roll approach leads to 1:13 (NDI) and 1:51 (CDI).