---
title: Cascade Decaying Dark Matter
url: https://www.emergentmind.com/topics/cascade-decaying-dark-matter
type: topic
---

# Cascade Decaying Dark Matter

Searching arXiv for recent and relevant papers on cascade decaying dark matter and related frameworks.
Cascade decaying dark matter denotes a class of dark-sector frameworks in which the observable consequences of dark matter are controlled by a sequence of decays, conversions, or induced breakups rather than a single direct transition to Standard Model states. In the literature, this includes multi-step dark-sector ladders ending in \(f\bar f\) final states, relic-abundance mechanisms driven by decays of nearly degenerate partners, continuous decays in a gapped continuum, and medium-induced disassembly of composite dark matter while traversing the Earth [1503.01773] [1607.03110] [2105.07035] [2512.16043].

## 1. Terminological scope and conceptual core

A common core runs through otherwise distinct constructions. Dark matter is not treated as an isolated stable particle with one dominant annihilation or decay channel into the Standard Model. Instead, the relevant dynamics are redistributed across intermediate states, so that either the relic abundance, the late-time signal, or the detector response is shaped by a cascade. In the simplest dark-sector ladder, a decay chain takes the form
\[
\chi \to \phi_n\phi_n \to 2\times \phi_{n-1}\phi_{n-1}\to \cdots \to 2^{n-1}\times\phi_1\phi_1 \to 2^n\times f\bar f,
\]
with \(\phi_i\) denoting lighter dark states and \(f\bar f\) the Standard Model endpoint [1503.01773].

The expression also appears in neutrino astronomy in a second, detector-level sense. PeV-scale decaying dark matter can yield neutrino fluxes that are dominated by IceCube cascade events, because \(\nu_e\) charged-current interactions, most \(\nu_\tau\) charged-current interactions, and all neutral-current interactions produce shower-like topologies rather than tracks [1308.1105]. A third usage arises in medium-induced cascades: loosely bound composite dark matter can dissociate constituent by constituent while crossing the Earth, producing a scattering cascade whose signatures are non-collinear multiscatters, timing-separated events, and possible coincidences across underground laboratories [2512.16043].

These usages are not identical, but they are structurally related. In all cases, the salient observable is set by energy flow through intermediate states or intermediate interactions, not by a one-step process.

## 2. Multi-step dark-sector cascades and spectral kinematics

The most developed kinematic framework is the multi-step dark-sector cascade studied for gamma-ray phenomenology. For decays, the first step is \(\chi\to\phi_n\phi_n\), and the photon spectrum is obtained from the decay differential rate according to
\[
\frac{1}{N_\gamma}\frac{dN_\gamma}{dE_\gamma}
=
\frac{1}{\Gamma}\frac{d\Gamma}{dE_\gamma}.
\]
For the same internal cascade, the decay spectrum has the same shape as the annihilation spectrum, except that the initial dark matter particle is twice as heavy [1503.01773].

The cascade is parameterized by
\[
\epsilon_i=\frac{2m_i}{m_{i+1}},\qquad \epsilon_f=\frac{2m_f}{m_1}.
\]
A hierarchical step has \(\epsilon_i\ll 1\), so the daughters are highly boosted; a degenerate step has \(\epsilon_i\to 1\), so the daughters are nearly at rest in the parent frame. For an \(n\)-step chain ending in a fermion of mass \(m_f\),
\[
m_\chi = 2^n\frac{m_f}{\epsilon_f\,\epsilon_1\epsilon_2\cdots\epsilon_n},
\]
which implies the self-consistency bound
\[
m_\chi \ge 2^n\frac{m_f}{\epsilon_f}.
\]
This fixes the maximum number of hierarchical steps allowed for a given \(m_\chi\) and final state [1503.01773].

In the hierarchical limit, the central recursion is the boost convolution
\[
\frac{dN_\gamma}{dx_i}
=
2\int_{x_i}^{1}\frac{dx_{i-1}}{x_{i-1}}\frac{dN_\gamma}{dx_{i-1}}
+\mathcal{O}(\epsilon_i^2),
\]
which recursively maps the \((i-1)\)-step spectrum to the \(i\)-step spectrum. Each hierarchical step doubles multiplicity, lowers the characteristic photon energy by roughly a factor of \(2\), and broadens the spectrum. By contrast, a degenerate step does not broaden the spectrum; it only doubles multiplicity and halves the energy scale. This distinction is central to model building, because hierarchical steps control spectral broadening while degenerate steps raise the mass scale without changing the spectral shape [1503.01773].

This directly modifies indirect-detection expectations. Hard 0-step spectra such as \(e^+e^-\) can become broad GeV-scale gamma-ray spectra after several steps, whereas soft 0-step spectra such as \(b\bar b\) require fewer steps. In the Galactic Center excess analysis, hierarchical cascades preferred \(m_\chi\) in the range \(20\)–\(150\ \mathrm{GeV}\) for all final states, while degenerate steps admitted much higher masses; the paper also emphasized that the observed morphology disfavors decaying, rather than annihilating, explanations because decays trace \(\rho\) rather than \(\rho^2\) [1503.01773].

## 3. Decay-controlled relic abundance and cosmological histories

Several frameworks use decays not merely as a late-time signal, but as the mechanism setting the dark matter abundance.

| Framework | Defining mechanism | Characteristic consequence |
|---|---|---|
| Co-decay | Stable \(A\), unstable \(B\), \(AA\leftrightarrow BB\), out-of-equilibrium \(B\) decay | Relic abundance set by decay-driven depletion |
| Vev flip-flop | \(\chi\) decay allowed only while \(\langle S_3\rangle\neq 0\) between weak-scale transitions | Temporary instability reduces an initially excessive abundance |
| Top-down injection | Defect-sourced injection \(Q(t)=Q_*(t/t_*)^{p-4}\) | Continuous non-thermal dark matter production |

In co-decaying dark matter, the dark sector contains a stable state \(A\) and a nearly degenerate unstable state \(B\), with rapid \(AA\leftrightarrow BB\) keeping them in chemical equilibrium while \(B\) decays out of equilibrium to the Standard Model. The density is depleted exponentially through the decay of \(B\), rather than by ordinary Boltzmann suppression, and the required annihilation cross section is therefore boosted relative to the standard WIMP case [1607.03110]. In the cosmological extension of this framework, the non-relativistic dark sector naturally drives an early matter-dominated phase whose duration is
\[
N_{\rm RH}-N_{\rm MD}\simeq \frac{1}{3}\log\!\left(\frac{\xi^4 m^4}{M_{\rm pl}^2\Gamma_B^2}\right)+1,
\]
and whose small-scale consequence is enhanced perturbation growth and microhalo formation. The characteristic microhalo mass is set by the reheating horizon,
\[
M_{\rm RH}\simeq 10^3\,M_\oplus\left(\frac{3~{\rm MeV}}{T_{\rm RH}}\right)^3
\left(\frac{10.75}{g_\ast(T_{\rm RH})}\right)^{1/2},
\]
which can generate substantial boost factors for indirect detection and point-source-like gamma-ray signals [1711.04773].

A distinct implementation is the “vev flip-flop” scenario. Here a fermionic singlet \(\chi\) freezes out while still relativistic and is therefore initially overabundant. As the universe cools, a scalar mediator \(S_3\) develops a vacuum expectation value, breaks the stabilizing symmetry, mixes \(\chi\) with electroweak-triplet fermions, and opens decay channels. When the Higgs subsequently acquires its own vacuum expectation value, portal terms restore the dark symmetry and make \(\chi\) stable again. The relic density is therefore set by the integrated decay rate during a finite temperature window bounded by two weak-scale phase transitions [1608.07578].

A third possibility is non-thermal top-down production from decaying topological defects. In that case, the dark matter source term is time dependent, with injection power
\[
Q(t)=Q_*\left(\frac{t}{t_*}\right)^{p-4},
\]
leading to a Boltzmann equation
\[
\frac{dY}{dx}=-\frac{A}{x^{n+2}}\left(Y^2-Y_{\rm eq}^2\right)+\frac{B}{x^{4-2p}}.
\]
For all \(p<3/2\), the paper derives a closed-form asymptotic yield and shows that topological defects can be the principal source of dark matter even when standard freeze-out underproduces the relic density, potentially producing large annihilation boost factors [1311.1637].

Taken together, these constructions show that cascade-decay ideas are not restricted to late decays of a pre-existing relic; they can reorganize the entire thermal history.

## 4. Indirect signatures: gamma rays, neutrinos, morphology, and detector cascades

In gamma-ray phenomenology, the decay flux is obtained by replacing the annihilation kernel with the decay kernel,
\[
\Phi_{\rm dec}(E_\gamma,l,b)
=
\frac{1}{4\pi m_\chi\tau_\chi}\frac{dN_\gamma}{dE_\gamma}\,J_{\rm dec}(l,b),
\qquad
J_{\rm dec}(l,b)=\int_0^\infty ds\,\rho(r),
\]
so the spectral machinery of dark-sector cascades carries over directly, but the morphology changes from \(\rho^2\) to \(\rho\). That distinction was explicitly noted as a reason the observed Galactic Center excess disfavors simple decaying explanations even though the cascade spectra themselves remain viable in other contexts [1503.01773].

For high-energy neutrinos, PeV-scale decaying dark matter was proposed as an explanation of the early IceCube events. A benchmark model with
\[
m_{\rm DM}=3.2~{\rm PeV},\qquad \tau_{\rm DM}=2\times10^{27}~{\rm s},
\]
and branching fractions \(b_{\rm H}=0.12\) to \(\nu_e\bar\nu_e\) and \(1-b_{\rm H}=0.88\) to \(q\bar q\) produces a hard cutoff at \(E_\nu=m_{\rm DM}/2\simeq1.6\ {\rm PeV}\), a bump near PeV energies, a dip between about \(0.3\) and \(1\ {\rm PeV}\), and a populated tail at \(30\)–\(400\ {\rm TeV}\). After oscillations, the flavor ratio at Earth is approximately \(1:1:1\), so the signal is cascade dominated in IceCube topology as well as shaped by QCD and electroweak particle cascades in the decay products [1308.1105].

Subsequent dedicated IceCube searches used six years of track data and two years of cascade data and found no significant excess attributable to decaying dark matter. The resulting lower limits exclude lifetimes shorter than \(10^{28}\ {\rm s}\) at \(90\%\) CL for dark matter masses above \(10\ {\rm TeV}\), with hard channels most strongly constrained [1804.03848]. KM3NeT was then projected to improve these tests in the PeV regime because of its Galactic Center visibility and improved cascade angular resolution; combining tracks and cascades, it was expected to produce world-leading limits on the decay lifetime and to test some of the dark matter interpretations of IceCube data [2007.03692].

A recurring source of confusion is that “cascade” in this literature can refer simultaneously to the dark-sector decay chain, the electroweak or hadronic shower generated by the final state, and the detector topology of the observed neutrino event. All three usages occur in the decaying-dark-matter neutrino literature.

## 5. Electromagnetic cascades before recombination and precision constraints

Cascade decays that inject energetic photons or \(e^\pm\) before recombination are constrained by CMB spectral distortions, but only after the full electromagnetic cascade is followed. For a decaying species \(X\) with fractional abundance \(f_X\), lifetime \(\tau_X\), and standard dark matter density \(\rho_{\rm DM}\), the injected power is
\[
\frac{dE}{dt}=f_X\,\rho_{\rm DM}\,\frac{1}{\tau_X}e^{-t/\tau_X}.
\]
The channels explicitly studied were \(X\to\gamma\gamma\) and \(X\to e^-e^+\), with the subsequent electromagnetic cascade evolved exactly rather than collapsed into instantaneous heating [1903.04503].

The distinction is important because relativistic electromagnetic cascades generate non-thermal relativistic spectral distortions, or \(ntr\)-type distortions, whose shape and amplitude differ from the standard \(y\)-, \(i\)-, and \(\mu\)-type templates. A significant fraction of the injected energy can be shifted into the high-frequency Wien tail, where COBE/FIRAS is less sensitive. As a result, the usual approximation that all injected energy becomes \(yim\)-type heating can overstate the constraint substantially. For decays at \(z_X\approx 5000\), the full treatment weakens the bound on \(f_X\) by a factor of about \(4.7\) for \(X\to e^-e^+\) at \(m_X\approx 4\ {\rm MeV}\), and by a factor of about \(4.1\) for \(X\to\gamma\gamma\) at \(m_X\approx 3\ {\rm MeV}\); more generally, the paper states that the relaxation can be as large as a factor of \(5\) [1903.04503].

This result has direct relevance for cascade-decaying dark matter more broadly. If a multi-step dark-sector chain ends in electromagnetic injection before recombination, the detailed internal cascade matters only through the final redshift distribution and energy spectrum of SM photons and \(e^\pm\). Once those are specified, the observable CMB distortion must be computed from the full electromagnetic cascade, not inferred from a pure-heating ansatz.

## 6. Extended variants: continuum spectra, instanton-mediated decay, and medium-induced disassembly

The idea of cascade-decaying dark matter extends beyond finite ladders of particles. In gapped continuum dark matter, the dark sector consists of a continuum of states above a mass gap \(\mu_0\), with spectral density
\[
\langle 0|\Phi(p)\Phi^\dagger(-p)|0\rangle
=
\int \frac{d\mu^2}{2\pi}\frac{i\,\rho(\mu^2)}{p^2-\mu^2+i\epsilon},
\qquad
\rho(\mu^2)=0\ \text{for}\ \mu^2<\mu_0^2.
\]
In the weakly interacting continuum model, dark matter interacts through a \(Z\)-portal, reproduces the observed relic density, avoids direct detection because continuum kinematics render low-energy scattering intrinsically inelastic, and exhibits continuous decays throughout cosmological history as well as cascade decays of continuum states produced at colliders [2105.07035]. Here the cascade is not a finite chain but a continuous flow in mass space from heavier continuum modes toward the gap.

A different route to metastable decays uses non-perturbative dark dynamics. In the dark-instanton model, a global symmetry protects a TeV-scale dark matter candidate \(\psi\) at the perturbative level, but an \(SU(2)_D\) instanton induces an operator violating that symmetry and generating leptophilic decays. The dominant channel is \(\psi\to \ell^+\ell^-\nu\), and the lifetime can naturally be tuned to \(\sim 10^{26}\ {\rm s}\) by choosing the dark gauge coupling \(g_D\) near \(1.15\) for \(m_\psi=3.5\ {\rm TeV}\) and \(v_D=4\ {\rm TeV}\) [1008.0642]. The same construction can be generalized into a genuine dark-sector cascade if some of the heavier \(\chi^{(i)}\) fermions that appear in the instanton vertex are taken lighter than \(\psi\), so that the instanton first produces on-shell dark daughters which then decay to the Standard Model.

The most radical extension is medium-induced. For loosely bound composite dark matter, each constituent–nucleus scatter in the Earth can impart enough energy to dissociate a constituent, so the Earth acts as a cascade medium. The mean free path is
\[
\lambda_{\rm mfp}(\vec x)=\left(\sum_A n_A(\vec x)\sigma_{Ad}\right)^{-1},
\]
the distance to the next scatter is sampled as
\[
L=-\ln(1-\zeta)\,\lambda_{\rm mfp},
\]
and the resulting random walk produces a spread radius
\[
R_{Sp}^{(N_s)}\approx
\left(\sum_A n_A\sigma_{Ad}\right)^{-1}
\sum_{n=1}^{N_s}\sin\!\left(\sqrt{n}\,\frac{\bar m_A}{m_d}\right).
\]
Underground signatures then include non-collinear multiple scatters in a single detector, parameter-dependent timing separations of multiscatter events, and coincident signals in different laboratories [2512.16043]. In this setting, “decay” is not an intrinsic instability but an interaction-triggered disassembly of the dark state.

The literature therefore uses cascade decaying dark matter as a genuinely broad organizing concept. It encompasses finite dark-sector ladders, decay-driven relic-abundance mechanisms, continuum spectra with perpetual intra-sector decays, and composite states whose cascades are activated by ordinary matter. What unifies these otherwise disparate constructions is the replacement of one-step dark-matter phenomenology by hierarchical energy transport through dark-sector intermediates.

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