---
title: Spiked Dark Matter Profile
url: https://www.emergentmind.com/topics/spiked-dark-matter-profile
type: topic
---

# Spiked Dark Matter Profile

Searching arXiv for recent and foundational papers on dark matter spikes around black holes.
A spiked dark matter profile is a centrally concentrated dark-matter overdensity that develops around a black hole, most often a supermassive black hole (SMBH), when the hole reshapes the phase-space distribution of the surrounding halo. In the canonical adiabatic-growth picture, an initial inner halo cusp is transformed into a much steeper inner density law inside a characteristic spike radius, but subsequent work has shown that this idealized description is not universal: relativistic dynamics, stellar heating, annihilation saturation, self-interactions, primordial-black-hole formation channels, and the microscopic nature of dark matter can all alter the spike or even turn enhancement into depletion [1701.00067] [2412.17919].

## 1. Canonical definition and analytic form

The standard starting point is an initial cuspy halo, often written as
\[
\rho_{\rm cusp}(r)=\rho_0\Bigl(\frac{r}{r_0}\Bigr)^{-\gamma_c},
\]
with inner slope \(\gamma_c\sim1.0\!-\!1.5\). If an SMBH grows adiabatically at the halo center, the density inside the spike radius \(R_{\rm sp}\) steepens to
\[
\rho_{\rm sp}(r)\propto r^{-\gamma_{sp}},\qquad
\gamma_{sp}=\frac{9-2\gamma_c}{4-\gamma_c},
\]
which gives \(\gamma_{sp}\approx2.3\!-\!2.4\) for \(\gamma_c=1.0\!-\!1.5\). For an NFW seed with \(\gamma_c=1\), this reduces to the familiar \(\gamma_{sp}=7/3\) result [1701.00067].

A widely used phenomenological representation is a single power law inside the spike radius,
\[
\rho(r)=\rho_{\rm sp}\Bigl(\frac{r}{r_{\rm sp}}\Bigr)^{-\gamma_{\rm sp}},
\]
valid for \(r<r_{\rm sp}\), with \(\rho_{\rm sp}\) defined at \(r=r_{\rm sp}\). In some gravitational-wave applications, the normalization is traded for \(\Sigma\equiv \rho_{\rm sp}r_{\rm sp}\), since \(\Sigma\) sets the overall normalization of the inner potential [2508.03803].

The canonical adiabatic result is not unique. For an initially finite-central-density or cored halo, the classic Gondolo–Silk-type construction gives a shallower spike, \(\rho\sim r^{-3/2}\), rather than \(r^{-7/3}\). Fully relativistic treatments of constant-density initial configurations recover precisely this \(r^{-3/2}\) behavior just outside the innermost bound region [2412.17919] [2503.09104].

The definition of “spike” is therefore structural rather than tied to a single slope: it denotes the black-hole-induced inner rearrangement of the halo, usually expressed as a broken profile that matches onto the outer halo at \(R_{\rm sp}\), but whose interior behavior depends on both dynamics and microphysics.

## 2. Characteristic scales, inner cutoffs, and saturation structure

For the Galactic Center, representative parameters used in spike studies are \(M_{\rm BH}\simeq4\times10^6\,M_\odot\), Schwarzschild radius \(r_{\rm Sch}\sim4\times10^{-7}\,\mathrm{pc}\), stellar-velocity dispersion \(\sigma\sim93\;\mathrm{km/s}\), and influence radius \(r_h\sim2\;\mathrm{pc}\). In idealized adiabatic models this implies
\[
R_{\rm sp}(0)\simeq0.2\,r_h\approx0.4\;\mathrm{pc},
\]
while a related Galactic-center parametrization gives \(R_h\simeq1.7\) pc and \(R_{\rm sp}\simeq0.2R_h\approx0.34\) pc for \(M_{\rm BH}=4.3\times10^6\,M_\odot\) and \(v_0\simeq105\) km/s [1701.00067] [2410.16379].

Practical spike profiles are almost always piecewise. One standard three-zone form consists of an inner annihilation plateau, an intermediate spike, and an outer cusp:
\[
\rho(r)=
\begin{cases}
\rho_{\rm core}, & 10\,r_{\rm Sch}<r\le r_{\rm core},\\
\rho_0\,(r/R_{\rm sp})^{-\gamma_{\rm sp}}, & r_{\rm core}<r\le R_{\rm sp},\\
\rho_0\,(r/R_{\rm sp})^{-\gamma_c}, & r>R_{\rm sp}.
\end{cases}
\]
The core radius is set by equating the annihilation time to the spike age,
\[
\frac{\rho(r_{\rm core})}{m_\chi}\langle\sigma v\rangle\sim (\tau\,t_{\rm heat})^{-1}.
\]
Closely related constructions define an annihilation-saturation density
\[
\rho_{\rm sat}\simeq \frac{m_\chi}{\langle\sigma v\rangle\,t_{\rm BH}},
\]
and a saturation radius \(R_{\rm sat}\) via \(\rho_{\rm sat}=\rho_{\rm spike}(R_{\rm sat})\) [1701.00067] [2410.16379].

In relativistic and orbit-based treatments, the innermost spike does not necessarily terminate in a flat plateau. Several Galactic-center analyses use a softened inner law
\[
\rho(r)\propto r^{-1/2}
\]
for \(4R_S\le r<R_{\rm sat}\), combined with a relativistic suppression factor such as \((1-4R_S/r)^3\) or \((1-2R_S/r)^{3/2}\), and set \(\rho=0\) below \(4R_S\) [2602.23348] [2410.16379].

These inner prescriptions matter because all indirect-detection observables scale with \(\rho^2\). In gamma-ray analyses the line-of-sight integral is
\[
J(\theta)=\int_{\rm l.o.s.}\rho^2(r)\,ds,
\]
and for a Majorana WIMP with \(\langle\sigma v\rangle=c_0+c_1(v/c)^2\), the velocity-suppressed term is enhanced near the black hole through the local virial relation \((v/c)^2\simeq r_{\rm Sch}/(2r)\) [1701.00067].

## 3. Formation channels and noncanonical spike classes

The literature now contains several distinct spike classes whose inner slopes and even qualitative behavior differ from the canonical adiabatic cusp.

| Scenario | Inner behavior | Defining feature |
|---|---|---|
| Adiabatic SMBH growth in cuspy halo | \(\rho\propto r^{-\gamma_{\rm sp}}\), \(\gamma_{\rm sp}=(9-2\gamma)/(4-\gamma)\) | Classical Gondolo–Silk construction |
| Adiabatic growth in cored/finite-density halo | \(\rho\propto r^{-3/2}\) | Finite central density |
| PBH mini-spike | \(\rho\propto r^{-9/4}\) between \(r_{\rm kd}\) and \(r_{\rm sp}\) | Early-universe turnaround/infall |
| Relativistic Bondi scalar spike | \(\rho_0\propto r^{-1.20},\,r^{-1.08},\,r^{-1.00}\) piecewise | Self-interacting dark scalar accretion |

For primordial black holes (PBHs), the mini-spike arises from early-universe infall rather than galactic adiabatic compression. In the cold-DM limit, the profile takes the form
\[
\rho(r)=
\begin{cases}
\alpha\,\rho_{\rm kd}, & 2r_s<r<r_{\rm kd},\\
\alpha\,\rho_{\rm sp}(r/r_{\rm sp})^{-9/4}, & r_{\rm kd}<r<r_{\rm sp},\\
\rho_{\rm bg}(r), & r>r_{\rm sp},
\end{cases}
\]
with \(\gamma_{\rm sp}=9/4\approx2.25\) and \(\alpha\simeq1.526\). This \(r^{-9/4}\) law is traced to adiabatic compression of an initially uniform medium and to self-similar secondary infall after matter–radiation equality [2406.07624].

When the orbital distribution and annihilation are treated self-consistently in PBH spikes, the central annihilation-modified region is not a flat plateau. For an initial exponential angular-momentum distribution of nearly radial orbits, orbit-averaged depletion produces a weak central cusp,
\[
\rho(r)\propto r^{-\gamma},\qquad \gamma\approx0.7,
\]
rather than \(\gamma=0\). The associated change in annihilation luminosity is an \(\mathcal{O}(1)\) correction rather than a qualitative breakdown of the spike picture [2407.10225].

Fermionic dark matter introduces a different departure from the standard power-law template. In the finite-temperature RAR model, dilute Boltzmannian fermions reproduce the classical \(\rho\sim r^{-3/2}\) spike, but semi-degenerate fermions with a dense core plus halo do not. The numerical profile has \(\alpha(r)\equiv-d\ln\rho/d\ln r\approx1.5\) immediately outside the capture region, then decreases to values \(\lesssim1\) over \(0.01\,r_{\rm sp}\lesssim r\lesssim r_{\rm sp}\). In this class, the SMBH does not always enhance the density; for sufficiently large black-hole mass and sufficiently heavy fermions it can instead deplete it [2412.17919].

A further noncanonical class appears in self-interacting dark scalar models treated through relativistic Bondi accretion. There the density cannot be fit by a single power law, but admits a three-segment form:
\[
\rho_0(r)\propto
\begin{cases}
r^{-1.20}, & r\sim r_s,\\
r^{-1.08}, & r_s<r<r_B,\\
r^{-1.00}, & r\sim r_B,
\end{cases}
\]
which is much shallower than the \(r^{-7/4}\) profile associated with Coulomb-like self-interactions [2112.05160].

These alternatives make clear that “spike” is not synonymous with “\(r^{-7/3}\)”. Inference based on a universal broken power law is therefore model dependent.

## 4. Relaxation, heating, depletion, and the possibility of density reduction

The main astrophysical mechanism that degrades a canonical spike is gravitational scattering by stars in the nuclear cluster. A widely used phenomenological model writes
\[
\rho(r,t)\approx \rho(r,0)\,e^{-\tau/2},\qquad \tau\equiv t/t_{\rm heat},\qquad t_{\rm heat}\sim10^9\;\mathrm{yr},
\]
so that the spike radius evolves as
\[
R_{\rm sp}(t)=R_{\rm sp}(0)\exp\!\Bigl[-\frac{\tau}{2(\gamma_{\rm sp}-\gamma_c)}\Bigr].
\]
For \(\tau\sim10\), the depletion factor is \(\kappa=e^{-\tau/2}\sim10^{-2.2}\), suppressing the normalization and shrinking the spike radius [1701.00067].

More broadly, violent processes such as galactic mergers, dynamical heating by stars, or self-interactions can flatten the spike to \(\gamma_{\rm sp}\sim0.5\). This places the classic adiabatic prediction \(\gamma_{\rm sp}\simeq7/3\) at one end of a much wider physical range [2508.03803].

Annihilation itself also regulates the innermost density. In Galactic-center orbit analyses, steep Gondolo–Silk spikes are eroded into a weak \(r^{-1/2}\) cusp inside the saturation radius, with
\[
\rho_{\max}\simeq \frac{m_{\rm DM}}{\langle\sigma v\rangle\,t_{\rm BH}}.
\]
For S2-based fits, the surviving NFW-seed spike is compatible with the stellar-orbit constraints only if
\[
\langle\sigma v\rangle \gtrsim 7.7\times10^{-27}\,\frac{\rm cm^3}{\rm s}\,
\Bigl(\frac{m_{\rm DM}}{100\,\rm GeV}\Bigr)
\]
at \(95\%\) confidence, so that annihilation has already reduced the density enough to evade the dynamical bounds [2303.09284].

The possibility of depletion rather than enhancement is particularly explicit in fermionic models. For semi-degenerate RAR halos, the spike peak density first rises and then falls as \(M_{\rm BH}\) increases; there is a threshold \(M_{\rm BH}^*(m_f)\) above which \(\rho_{\max}<\rho_0\). Moreover, there is a critical fermion mass \(m_f^*\simeq300\) keV such that for \(m_f>m_f^*\) no adiabatic spike exceeds the original central density [2412.17919].

A common misconception is therefore that a central black hole necessarily creates a larger annihilation signal. The contemporary literature does not support that as a generic statement: the sign and magnitude of the effect depend on heating history, annihilation saturation, initial halo structure, and dark-matter phase-space properties.

## 5. Observational diagnostics and empirical constraints

Indirect searches remain the most developed probe because annihilation rates scale as \(\rho^2\). For the Galactic Center, a benchmark Fermi-LAT comparison integrates the flux over \(1\!-\!100\) GeV and a small angle \(\sim0.1^\circ\), using the point source Sgr A\(^*\) with observed flux \(\Phi_{\rm Fermi}\simeq2.18\times10^{-8}\,\mathrm{ph/cm^2/s}\). In such analyses, the relative enhancement over a pure NFW profile can exceed \(10^2\!-\!10^5\), especially when the \(c_1v^2\) term dominates. Yet the allowed dark-matter parameter space depends strongly on the spike assumptions: for a depleted spike with \(\gamma_c\le1.3\), even a thermal relic with \(c_0\simeq3\times10^{-26}\,\mathrm{cm^3/s}\) produces flux several orders of magnitude below \(\Phi_{\rm Fermi}\), whereas idealized adiabatic spikes exclude \(m_\chi\lesssim(25,80,240)\) GeV for \(\gamma_c=1.0,1.1,1.2\) [1701.00067].

More recent line-search recasts sharpen this conclusion for thermal WIMPs. Photon-line constraints from Fermi-LAT and MAGIC exclude the entire Gondolo–Silk range \(\gamma_{\rm sp}\in[2.25,2.5]\) for \(m_\chi\in[10\,\mathrm{GeV},10^5\,\mathrm{GeV}]\), even when the \(\gamma\gamma\) branching fraction is only \(1\%\). Neutrino-line searches with IceCube give a complementary bound, with the strongest limit \(\gamma_{\rm sp}\lesssim2.35\) around \(m_\chi\sim2\!-\!4\) TeV [2602.23348].

Radio and microwave emission provide a second \(\rho^2\)-weighted channel. For a spike with \(\gamma_{\rm sp}=7/3\), the small-angle surface-brightness scaling is \(I_\nu(\theta)\propto\theta^{-11/3}\), much steeper than the NFW expectation \(I_\nu(\theta)\propto\theta^{-1}\). In the specific Galactic-center synchrotron calculations with direct \(e^+e^-\) annihilation, a \(10\) GeV WIMP plus spike yields \(I_\nu(0.1^\circ)\simeq10^6\) Jy sr\(^{-1}\) at \(30\) GHz, compared with \(10^2\) Jy sr\(^{-1}\) for NFW. Planck and future SKA-like observations were proposed as probes of \(R_{\rm sp}\), \(r_{\rm sat}\), and \(\gamma_{\rm sp}\) [1311.0139].

The \(511\) keV bulge line has also been modeled with a Galactic-center spike. Using a four-zone piecewise profile with either a Gondolo–Silk spike or a stellar-heated spike softened to \(\gamma_{\rm sp}=1.5\), dark matter with mass up to approximately \(20\) MeV can reproduce most of the observed bulge intensity while remaining compatible with disk-emission and in-flight-annihilation constraints, provided the disk component is dominated by an astrophysical low-energy positron source [2410.16379].

Dynamical probes constrain the spike independently of annihilation physics. Fits to Keck and VLT S-star data rule out an initial slope \(\gamma\gtrsim0.92\) for the generalized NFW spike at \(95\%\) confidence. For \(\gamma=1\), the analyses find \(R_{\rm sp}<15.7\) pc at \(95\%\), excluding the corresponding Gondolo–Silk prediction \(R_{\rm sp}^{\rm GS}\simeq18.6\) pc, and also obtain \(\gamma_{\rm sp}<2.32\) at \(95\%\), excluding \(\gamma_{\rm sp}^{\rm GS}=7/3\) [2303.09284].

Gravitational-wave observables are emerging as a particularly clean spike diagnostic. For binaries orbiting inside an SMBH spike, the dark-matter-induced center-of-mass acceleration produces a secular Doppler modulation of the waveform. Fisher plus Bayesian forecasts for LISA and DECIGO indicate that \(\gamma_{\rm sp}\) can be measured to a few-percent precision when \(\gamma_{\rm sp}\gtrsim1.8\), improving to \(\sim1\!-\!2\%\) for \(\gamma_{\rm sp}\sim7/3\), while even \(\gamma_{\rm sp}\sim1.5\) can be constrained at the \(\sim10\%\) level. In that setup, dynamical friction and tidal effects are negligible compared with the conservative dark-matter potential [2508.03803].

A distinct extragalactic dynamical route is reverberation mapping. In a sample of fourteen AGN, five objects show \(1\!-\!2\sigma\) evidence for an enclosed mass that grows with radius across multiple emission lines. A joint fit gives a preferred universal dark-matter profile
\[
\rho_{\rm DM}(r)=\rho_0\,r^{-\gamma},\qquad \gamma_{\rm global}\simeq1.6\pm0.2,
\]
consistent with a mildly relaxed spike [2506.10122].

## 6. Relativistic, numerical, and geometric refinements

Relativistic phase-space treatments do not simply reproduce the Newtonian broken-power-law profile. In Schwarzschild geometry, a numerical fit to the spike density over \(r\ge4GM\) is
\[
\rho_{\rm sp}(r)=\kappa\,(r/GM)^{-\omega}(1-4GM/r)^{\eta},
\]
with best-fit parameters that are nearly universal across several halo and black-hole configurations: \(\omega=2.00\pm0.05\) and \(\eta=2.1\pm0.1\). This corresponds to a quasi-isothermal \(r^{-2}\) cusp over \(r\sim5\!-\!50\,GM\) together with a relativistic cutoff at \(r=4GM\) [2403.18529].

Spin modifies the spike further. In the exact Kerr geometry, black-hole rotation increases the dark-matter density close to the hole despite angular-momentum transfer to the halo. For \(\gamma_c=1\), the nonrotating relativistic value is \(\gamma_{\rm sp}=7/3\), while an empirical fit gives
\[
\gamma_{\rm sp}(a)\simeq \gamma_{\rm sp}(0)\Bigl[1+\delta(\gamma_c)\frac{a}{m}\Bigr],
\]
with \(\delta(\gamma_c)\approx0.10(2-\gamma_c)\); at \(a/m=0.8\), one finds \(\gamma_{\rm sp}\simeq2.36\) and a normalization increase of about \(25\%\) at the influence radius for \(\gamma_c=1\) [1707.06302].

The most significant numerical revision to the classical picture comes from fully simulated cold-dark-matter spikes in Hernquist halos. Using 218 SWIFT N-body runs, an empirical one-parameter “spike + Hernquist” profile was proposed:
\[
\rho_{\rm final}(r;\mu)=\rho_H(r)\Bigl[\beta(\mu)+(r/r_{\rm sp}(\mu))^{\gamma_{\rm sp}+1}\Bigr],
\]
where \(\mu=M_{\rm BH}/M_{\rm tot}\), \(\beta(\mu)=1-1.10\,\mu^{0.884}\), and
\[
r_{\rm sp}/a=0.789\,[\mu/(\mu+0.156)]^{1.25}.
\]
This differs sharply from the classic \(r_{\rm sp}\propto \mu^{1/2}\) scaling: the simulations give \(r_{\rm sp}\propto \mu^{1.25}\) at low \(\mu\) and saturation at \(r_{\rm sp}\to0.79a\) for \(\mu\gg0.1\), together with \(10\!-\!20\%\) outer-halo depletion for \(\mu\lesssim0.2\) [2411.12007].

Halo modeling can dominate observational inference. In M87, the assumption of an NFW host plus \(\gamma_{\rm sp}=7/3\) spike yields extremely strong annihilation limits, excluding thermal \(s\)-wave dark matter up to \(\mathcal{O}(100)\) TeV. But adopting a cored Burkert halo changes the normalization at the spike radius by about five orders of magnitude; in that case the smooth halo overwhelmingly dominates the annihilation signal, with \(\bar Q\ll \bar J\), and the spike contribution becomes negligible [1505.00785] [2411.18751].

At the level of spacetime backreaction, the spike is not accurately represented by rest-mass density alone. In an Einstein-cluster treatment of a Hernquist-seeded Milky-Way spike, the kinetic term contributes about \(50\%\) of the peak energy density, the stress tensor is mildly anisotropic, all standard energy conditions are satisfied, and the resulting metric deviations from Schwarzschild reach fractional levels \(10^{-5}\!-\!10^{-6}\), approximately \(2.5\times\) larger than when only the mass density is used as the source [2511.12570].

The present status is therefore mixed. Classical power-law spikes remain a useful organizing approximation, but current work shows that their slope, normalization, radial extent, and observational impact are all highly contingent on the underlying halo model, the black-hole growth history, the relativistic orbital structure, and the nature of the dark matter itself. This suggests that robust use of spiked profiles in indirect detection, stellar dynamics, or gravitational-wave inference requires system-specific modeling rather than a universal Gondolo–Silk template.

Source: https://www.emergentmind.com/topics/spiked-dark-matter-profile