
measurements might not produce as large deviations, but show a trend being significantly
lower compared to CMB data. To account for this, and leaving aside the possibility of
a statistical fluctuation, there either needs to be some unaccounted systematic error (see
e.g. [11]) or an alternative to ΛCDM featuring a suppression of the matter power spectrum
in the k∼0.1−1h/Mpc regime.
One model to achieve such a suppression is the Decaying Cold Dark Matter (DCDM)
model. It is based on the hypothesis that dark matter can decay on cosmological time-
scales into secondary dark sector particles. The decay products are assumed to be effectively
stable on cosmological scales and, like the dark matter itself, sufficiently weakly coupled to
visible matter to escape (in-)direct detection. However, the kinetic energy released in the
decay process counteracts the growth of structures and leads to a suppression of the power
spectrum. The model has been mainly investigated in two variants: a decay into massless
secondaries that act as dark radiation (DR), or into a massless and a massive daughter.
Depending on the mass splitting between mother and massive daughter particles, the latter
acts as warm dark matter (WDM) being gradually produced in the decay process in the
late universe. For both variants, the evolution in the early universe is identical to ΛCDM,
thereby preserving its success in explaining the CMB and LSS on very large scales. Both
models were studied regarding the S8and also the Hubble tension [12], taking CMB, BAO
and recently also galaxy clustering data into account. The decay into massless secondaries
was investigated e.g. in [13–23], allowing also the possibility that only a fraction fof the
dark matter decays. It was found that the decay into purely DR will most likely not be
able to solve cosmological tensions and requires minimum lifetimes of around ∼200 Gyrs for
f= 1. The latest works [22,23] for this model confirm this even more while also providing
tight constraints. For lifetimes shorter than the age of the universe, [22] finds f < 2.16%
and for f→1 a lower bound of τ > 250 Gyrs [22,23]. As a further variant, also the decay
of warm dark matter mother particles into massless dark radiation has been considered [24],
but similar to the decay of CDM into massless daughters, this setup was found to neither
solve the H0nor S8tensions [25].
The model with a decay of CDM into WDM and DR, on which we mainly focus in this
work, is apart from the lifetime τdescribed by the mass splitting parameter
=1
21−m2
M2,(1.1)
involving the mass of the mother (M) and the massive daughter (m) particle. The results
regarding the Hubble tension are similar compared to the massless case, implying that it is
probably not to be resolved with DCDM [26–28] (see also [29,30] for earlier work). The
situation for the S8tension is however not so clear. While [28] suggests that this tension can
also not be addressed, [31], which includes an improved treatment of perturbations, finds that
it actually can be lessened for τ∼55 Gyrs and ∼0.7% based on Planck CMB, BAO, RSD
and SN Ia data. In two follow-up works a newly developed code for much faster computation
of the DCDM power spectra is used. This allows for a more in depth analysis, like in [32],
where a mild preference for DCDM is found depending on the priors for S8. The latest
work [22] also includes full-shape information from BOSS DR12 galaxy clustering, and finds
that DCDM can ease the S8tension even though it is not performing significantly better than
ΛCDM when disregarding KiDS data. However, when including KiDS, DCDM is preferred,
and the best-fit model occurs for a lifetime of τ∼120 Gyrs and ∼1.2%.
– 2 –