---
title: Decaying Dark Matter Cosmological Constraints
url: https://www.emergentmind.com/topics/decaying-dark-matter-cosmological-constraints
type: topic
---

# Decaying Dark Matter Cosmological Constraints

Decaying dark matter (DDM) cosmologies postulate that a fraction of the dark matter sector is unstable and decays into lighter species—typically assumed to be invisible, relativistic “dark radiation.” Such models are motivated by particle physics frameworks and are also explored as possible solutions to small-scale structure anomalies and tensions in ΛCDM, such as the $S_8$ and $H_0$ discrepancies. The decay process alters the expansion history, the growth and distribution of cosmic structures, CMB anisotropies, and leaves potentially observable signatures across a wide array of cosmic data sets. Cosmological probes at both linear and nonlinear scales, together with laboratory and astrophysical limits, allow for stringent constraints on DDM lifetimes, branching fractions, and daughter particle properties.

## 1. Theoretical Framework and Parameterization

DDM models fundamentally split the dark sector into stable (SDM) and decaying (DCDM) components. The primary parameters are:

- **Fractional abundance**: $f \equiv \Omega_{\mathrm{dcdm}}/\Omega_{\mathrm{dm}}$ is the initial fraction of (cold) DM that is unstable.
- **Decay rate (lifetime)**: $\Gamma$ in Gyr$^{-1}$, $\tau = \Gamma^{-1}$.
- **Nature of decay products**: Most cosmological analyses consider invisible dark radiation (massless or very light, noninteracting), though models with decays to massive daughters also exist.
- **Energy splitting**: For two-body decays (parent mass $m_0$ to massive daughter $m_2$ and massless $\phi$), the fraction of rest-mass energy transferred to relativistic products is $\epsilon = \frac{1}{2}\big[1-(m_2/m_0)^2\big]$, with associated kick velocity $v_k \approx \epsilon c$.

The homogeneous evolution equations (in conformal or cosmic time) for decaying and stable dark matter, and the dark radiation component, are
\[
\begin{aligned}
\frac{d\rho_\mathrm{dcdm}}{dt} &= -3H\rho_\mathrm{dcdm} - \Gamma \rho_\mathrm{dcdm}, \\
\frac{d\rho_\mathrm{dr}}{dt} &= -4H\rho_\mathrm{dr} + \epsilon \Gamma \rho_\mathrm{dcdm},
\end{aligned}
\]
with similar equations for massive daughters in nontrivial decay kinematics. The impact on perturbations is captured in a modified Boltzmann hierarchy, with additional source terms determined by the decay kinematics and branching.

## 2. Observational Signatures and Sensitive Probes

The specific cosmological consequences of DDM depend on the decay rate, the fraction of decaying DM, and the products' properties. Major observable effects include:

- **CMB anisotropies**: DDM modifies the integrated Sachs-Wolfe effect (late-time potential decay), alters the lensing kernel, and shifts the acoustic peak structure by changing the matter–radiation content. The effect is most visible at low multipoles ($\ell \lesssim 30$) and, for sufficiently early decays, can shift the position of the peaks. The derived constraints primarily come from Planck temperature, polarization, and lensing power spectra [1606.02073, 2205.05636, 2006.03678, 2307.00418, 2102.12498].

- **Structure growth and matter power spectrum**: The removal of clustered matter and injection of free-streaming (relativistic or semi-relativistic) daughters suppresses the amplitude of matter fluctuations, leading to a reduction in $\sigma_8$ and a suppression of $P(k)$ at scales $k$ above the daughter free-streaming horizon. The scale and depth of suppression encode both $\tau$ and $v_k$ (or $\epsilon$) [1402.2972, 1503.05682, 1011.4970, 1001.3870, 2102.12498].

- **Lyman-$\alpha$ forest**: Decay-generated free-streaming suppresses small-scale power at high redshift ($z\sim2-4$), providing some of the strongest constraints, especially for kick velocities $V_k \sim 10-100$ km/s and lifetimes less than the Hubble time [1309.7354].

- **Halo abundance, internal structure, and substructure**: N-body simulations show that sufficient DDM suppresses the abundance of massive halos and subhalos, alters the concentration–mass relation, and can lead to observable features in the cluster mass function and galaxy counts. This directly impacts cluster abundance studies, galaxy clustering, and Milky Way satellite populations [1011.4970, 1001.3870, 1503.05682, 1906.09112].

- **Redshift-space distortions and kSZ effect**: Direct measurement of growth (e.g., $f\sigma_8$), and velocity fields (kSZ) are also suppressed in DDM, with future surveys projected to reach sub-percent level sensitivity to DDM signatures [1908.02668].

- **Astrophysical gamma rays and cosmic-ray constraints**: In scenarios with decays to SM particles, especially photons, electrons, or quarks, one obtains stringent bounds on $\tau$ (up to $\sim 10^{28}$ s) from the isotropic gamma-ray background, X-ray and cosmic-ray measurements [1811.05988, 1205.6783, 1610.06933].

## 3. Constraints Across Decay Regimes

The constraints are typically segmented by the regime of the decay lifetime:

| Regime                | Main Bound  | Principal Data Sets                            | Physical Consequence             |
|-----------------------|-------------|-----------------------------------------------|----------------------------------|
| Short-lived (pre-recomb.) | $f\lesssim2-5\%$ (lifetime $\tau<10^6$ yr) | CMB, $H_0$, BAO                    | Early energy injection, alters acoustic scales, shifts $H_0$ up by $\sim1$ km/s/Mpc at best |
| Intermediate ($10^3$ yr $< \tau < t_0$) | $f\lesssim3-4\%$ ($\tau \sim t_0$) | Planck TTTEEE (+lensing+BAO), LSS, clusters | Enhanced ISW, late-time $P(k)$ suppression, $S_8$ and $H_0$ remain largely unaffected |
| Long-lived ($\tau \gg t_0$) | $f\Gamma \lesssim (5.9-15.9)\times10^{-3}$ Gyr$^{-1}$ ($\tau>150-250$ Gyr) | Planck TTTEEE+lensing+BAO           | Effectively indistinguishable from CDM; tightest bound on instability fraction |

These values are robust to moderate extensions (varying $M_\nu$, curvature, tensor sector) [2205.05636, 1606.02073, 2203.07440]. If only a fraction $f$ of the DM is unstable, for $\tau \gtrsim t_0$ the bound is effectively on the product $f\Gamma$.

## 4. Statistical Methodologies and Analysis Pipelines

Modern analyses employ Boltzmann codes (CLASS/CAMB), often with custom implementations of the DDM Boltzmann hierarchy to accommodate the energy/momentum flow between species with arbitrary masses [1606.02073, 1402.2972, 1510.06026, 1106.1984]. Parameter estimation is carried out via Markov Chain Monte Carlo techniques (e.g., MontePython, CosmoMC), varying standard ΛCDM and DDM-specific parameters under flat or log priors. The principal observational datasets include:

- CMB: Planck (2015 and 2018) TT, TE, EE, lensing
- Large-scale structure: BAO (BOSS, WiggleZ, eBOSS), $f\sigma_8$, full-shape galaxy power spectra (EFTofLSS)
- Supernovae: Pantheon, JLA
- Lyman-$\alpha$: SDSS, VHS, BOSS 1D flux power
- Weak lensing: KiDS, CFHTLenS, DES
- SZ cluster counts: Planck
- Cosmic shear: KiDS450
- Astrophysical indirect detection: Fermi IGRB, AMS-02, extragalactic X-ray background

Nuisance parameters are marginalized according to the implementations of each survey's likelihood [2205.05636, 1606.02073, 2203.07440, 2102.12498].

## 5. Physical Implications and Tensions

### $S_8$ and $H_0$ Tensions

Decaying DM generally suppresses $\sigma_8$ and increases $H_0$ compared to ΛCDM, aligning directionally with the observed tensions. However, multiple analyses consistently show that the improvement in tension is modest and statistically insignificant with current data [1606.02073, 1908.02668, 2307.00418, 2102.12498]:
- $H_0$ can increase by $\sim1$ km/s/Mpc—but is still $2.5-3\sigma$ below local measurements.
- $S_8$ can be shifted down by up to 3\%, but not enough to fully reconcile lensing and CMB-inferred constraints.

The preferred parameter spaces capable of further alleviating these tensions are essentially excluded by CMB and LSS observables due to excessive ISW effects and $P(k)$ suppression (see [2006.03678, 2307.00418]).

### Small-Scale Structure

For DDM models designed to address the overabundance of small galaxies or too-dense halo cores, Lyman-$\alpha$ forest and halo mass function measurements robustly exclude lifetimes $<40$ Gyr for kick velocities $V_k>20-100$ km/s, consistent with the need to avoid substantial suppression of small-scale power [1011.4970, 1503.05682, 1309.7354, 1402.2972].

### Astrophysical Constraints

For models with decays to SM particles, the most stringent bounds arise from extragalactic gamma-ray backgrounds and antiproton fluxes: $\tau > (1-5)\times10^{28}$ s for a variety of decay channels and $m_{DM}\sim$ GeV–EeV [1811.05988, 1205.6783, 1610.06933], which are far stronger than cosmological limits for such scenarios.

## 6. Open Questions, Future Prospects, and Model Extensions

- **Nonlinear structure**: Present limits rely predominantly on linear cosmological observables or simplified halo models. Improved N-body simulations (including baryons and velocity effects of daughter particles) are actively pursued to robustly characterize DDM on nonlinear and small scales [1503.05682, 1906.09112].
- **Complex decay channels**: Phenomenology with massive daughters, partial/fractional decay, and time-dependent rates is being systematically explored. Some mixed decay models may permit additional parameter space, but generally only for highly fine-tuned or subdominant fractions [1402.2972, 1106.1984, 1510.06026, 2203.07440].
- **Next-generation data**: Upcoming observatories (Euclid, DESI, LSST, CMB-S4) will measure growth, lensing, velocity fields, and the CMB lensing power to unprecedented precision, shrinking the allowed parameter volume for DDM or potentially uncovering signatures in the event of deviations from ΛCDM [1606.02073, 2205.05636, 1908.02668, 2102.12498].
- **Systematics vs new physics**: There is an active effort to disentangle unmodeled systematics in low-$z$ observations (for instance, in cluster mass calibrations or lensing systematics) from true new-physics signatures [1606.02073].
- **PBH mergers**: DDM constraints are particularly restrictive for scenarios where PBHs constitute all DM, as merging events efficiently generate "dark radiation" exceeding permitted fractions [1606.02073].

## 7. Representative Numerical Constraints

Selected robust 95% C.L. constraints derived from recent literature:

| Dataset Combination                        | Regime                | Bound                                                          | Reference            |
|---------------------------------------------|-----------------------|----------------------------------------------------------------|----------------------|
| Planck18 TT,TE,EE + lensing + BAO          | $f_{\rm dcdm}=1$      | $\tau > 246$ Gyr                                               | [2205.05636]         |
| As above, fraction $f$ decaying            | $f\ll1$, $\tau \gg t_0$ | $\tau/f > 246$ Gyr                                             | [2205.05636]         |
| Planck15 + BAO + WiggleZ                   | long-lived ($\xi$)    | $f\Gamma < 5.9\times10^{-3}$ Gyr$^{-1}$, $\tau>170$ Gyr        | [1606.02073]         |
| Planck15 + BAO + RSD                       | short-lived           | $f<2.73\%$                                                     | [1908.02668]         |
| SDSS Lyman-$\alpha$                        | $\tau < 10$ Gyr       | $V_k<30-70$ km/s                                               | [1309.7354]          |
| Lyman-$\alpha$, large $V_k$                | $V_k\gg 100$ km/s     | $\tau>40$ Gyr                                                  | [1309.7354]          |
| Planck15+BAO+SZ+weak lensing               | full DDM              | $\tau>175$ Gyr                                                 | [1906.09112]         |
| Fermi IGRB (γ-ray, SM decays)              | GeV–EeV DM            | $\tau>10^{28}$ s                                               | [1811.05988, 1205.6783]  |

## References

- [1606.02073]  
- [2205.05636]
- [2006.03678]
- [1402.2972]
- [1503.05682]
- [1011.4970]
- [1001.3870]
- [2203.07440]
- [2102.12498]
- [1908.02668]
- [1309.7354]
- [1106.1984]
- [1811.05988]
- [1205.6783]
- [1906.09112]
- [2307.00418]
- [1610.06933]

Source: https://www.emergentmind.com/topics/decaying-dark-matter-cosmological-constraints