Search

Traditional computations of the dark matter (DM) relic abundance, for models where attractive self-interactions are mediated by light force-carriers and bound states exist, rely on the solution of a coupled system of classical on-shell Boltzmann equations. This idealized description misses important thermal effects caused by the tight coupling among force-carriers and other charged relativistic species potentially present during the chemical decoupling process. We develop for the first time a comprehensive ab initio derivation for the description of DM long-range interactions in the presence of a hot and dense plasma background directly from nonequilibrium quantum field theory. Our results clarify a few conceptional aspects of the derivation and show that under certain conditions the finite temperature effects can lead to sizable modifications of DM Sommerfeld-enhanced annihilation and bound-state decay rates. In particular, the scattering and bound states get strongly mixed in the thermal plasma environment, representing a characteristic difference from a pure vacuum theory computation. The main result of this work is a novel differential equation for the DM number density, written down in a form which is manifestly independent under the choice of what one would interpret as a bound or a scattering state at finite temperature. The collision term, unifying the description of annihilation and bound-state decay, turns out to have in general a nonquadratic dependence on the DM number density. This generalizes the form of the conventional Lee-Weinberg equation which is typically adopted to describe the freeze-out process. We prove that our number density equation is consistent with previous literature results under certain limits. In the limit of vanishing finite temperature corrections our central equation is fully compatible with the classical on-shell Boltzmann equation treatment. So far, finite temperature corrected annihilation rates for long-range force systems have been estimated from a method relying on linear response theory. We prove consistency between the latter and our method in the linear regime close to chemical equilibrium. [ABSTRACT FROM AUTHOR]

We explore the phenomenology of having a second epoch of dark matter annihilation into dark radiation long after the standard thermal freeze-out. Such a hidden reannihilation process could affect visible sectors only gravitationally. As a concrete realization we consider self-interacting dark matter (SIDM) with a light force mediator coupled to dark radiation. We demonstrate how resonantly Sommerfeld enhanced cross sections emerge to induce the reannihilation epoch. The effect is a temporally local modification of the Hubble expansion rate, and we show that the cosmic microwave background (CMB) measurements--as well as other observations--have a high sensitivity to observe this phenomenon. Special attention is given to the model region where late kinetic decoupling and strong self-interactions can alleviate several small-scale problems in the cold dark matter paradigm at the same time. Interestingly, we find that reannihilation might here also simultaneously lower the tension between CMB and low-redshift astronomical observations of H0 and σ8. Moreover, we identify reannihilation as a clear signature to discriminate between the phenomenologically otherwise almost identical vector and scalar mediator realizations of SIDM. [ABSTRACT FROM AUTHOR]

Calculating the abundance of thermally produced dark matter particles has become a standard procedure, with sophisticated methods guaranteeing a precision that matches the percent-level accuracy in the observational determination of the dark matter density. Here, we point out that one of the main assumptions in the commonly adopted formalism, namely local thermal equilibrium during the freeze-out of annihilating dark matter particles, does not have to be satisfied in general. We present two methods for how to deal with such situations, in which the kinetic decoupling of dark matter happens so early that it interferes with the chemical decoupling process: (i) an approximate treatment in terms of a coupled system of differential equations for the leading momentum moments of the dark matter distribution, and (ii) a full numerical solution of the Boltzmann equation in phase space. For illustration, we apply these methods to the case of scalar singlet dark matter. We explicitly show that even in this simple model the prediction for the dark matter abundance can be affected by up to 1 order of magnitude compared to the traditional treatment. [ABSTRACT FROM AUTHOR]