Probing the Distribution and Nature of Dark Matter Around Supermassive Black Holes from EMRI and IMRI Gravitational Waves

Abstract: The distribution of dark matter in the immediate vicinity of supermassive black holes remains poorly understood despite its importance for galaxy evolution and precision tests of gravity. Future space-based gravitational-wave observatories offer a unique opportunity to probe this relativistic regime through the inspiral of compact objects into supermassive black holes. Building upon our previously constructed exact Einstein-cloud solutions within General Relativity, we develop a fully relativistic framework to investigate the gravitational-wave signatures of collisionless dark-matter halos surrounding supermassive black holes. The framework provides a unified treatment of orbital dynamics, adiabatic inspiral, accumulated gravitational-wave cycles, waveform phase evolution, signal-to-noise ratio, and waveform mismatch for extreme- and intermediate-mass-ratio inspirals (EMRIs/IMRIs). As a representative application, we specialize the formalism to Model I. We show that relativistic dark-matter halos can produce measurable modifications to the accumulated gravitational-wave cycles, waveform phase, signal-to-noise ratio, and waveform mismatch, leading to consistent conclusions regarding detectability. By separating conservative modifications of the spacetime geometry from dissipative effects due to relativistic dynamical friction, we find that the observable signatures are dominated by the former, while the latter remains negligible for the halo models considered. These results demonstrate that future gravitational-wave observations by LISA and similar missions may provide a powerful probe of the relativistic distribution and physical nature of dark matter around supermassive black holes.

• The paper develops a relativistic Einstein-cloud framework showing that dark matter modifies EMRI and IMRI orbital frequencies, inspiral phases, and gravitational-wave signals around supermassive black holes.
• For a 10^6-solar-mass black hole, densities of 0.3–10 GeV cm⁻³ produce measurable dephasing or waveform mismatch within roughly 1.5–3.6 years, while stronger densities reach a one-radian phase shift before SNR 8.
• The paper finds conservative changes to spacetime geometry dominate the observable signal by many orders of magnitude over dynamical friction, emphasizing the need for environmental effects in future LISA waveform models.

Motivation and context

The distribution of dark matter (DM) in the strong-gravity regime immediately surrounding supermassive black holes (SMBHs) is among the least constrained aspects of galactic astrophysics. Electromagnetic probes—rotation curves, lensing, stellar orbits near Sagittarius A$^\ast$, and horizon-scale imaging—constrain DM only indirectly or on larger scales. This paper develops a fully relativistic framework for assessing whether extreme- and intermediate-mass-ratio inspirals (EMRIs/IMRIs), observed by LISA-class detectors, can directly probe the relativistic DM environment of SMBHs. The work builds on the authors' previously constructed exact Einstein-cloud solutions (Shen et al., 2023), which are static, spherically symmetric spacetimes sourced by collisionless massive particles on circular timelike geodesics.

The Einstein-cloud halo model

The background is a static, spherically symmetric metric characterized by a lapse function $f(r)$ and Misner–Sharp mass $m(r)$. The stress-energy tensor is anisotropic, $T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$, with vanishing radial pressure arising naturally from the circular-geodesic microscopic dynamics rather than from an ad hoc equation of state. Einstein's equations yield three independent relations among $f$, $m$, $\rho$, and $P$, so one additional input closes the system: the density profile

$$\rho(r)=\rho_0\,\frac{\left(1-\frac{4M}{r}\right)^n H(r-4M)}{\left(\frac{r}{a}\right)^\gamma \left[1+\left(\frac{r}{a}\right)^\alpha\right]^{(\beta-\gamma)/\alpha}},$$

a relativistic generalization of the double-power-law family that unifies cored ($\gamma=0$) and cuspy ($f(r)$0) halos. The Heaviside factor enforces the relativistic inner boundary $f(r)$1, motivated by the fully relativistic phase-space analyses of Sadeghian et al. (Sadeghian et al., 2013) and Speeney et al. (Speeney et al., 2022), which showed that no stationary equilibrium density exists inside $f(r)$2 for collisionless particles in Schwarzschild spacetime. This contrasts with phenomenological profiles such as NFW or Hernquist, which predict nonvanishing density arbitrarily close to the horizon when extrapolated inward. As a representative application, the paper adopts Model I, corresponding to $f(r)$3—a cuspy model with an exact analytical mass function $f(r)$4 outside $f(r)$5.

Inspiral formalism

The compact object of mass $f(r)$6 is treated as a test particle on an adiabatic sequence of quasicircular geodesics. Circular-orbit energy, angular momentum, and frequency are derived exactly in the halo geometry; notably, $f(r)$7 generalizes Kepler's third law, and all observables depend on the halo through this modified orbital frequency. The ISCO follows from $f(r)$8, reducing to $f(r)$9 in vacuum. Dissipation includes both quadrupolar gravitational-wave emission and relativistic dynamical friction (DF), using a Chandrasekhar-type prescription augmented by the relativistic correction factor $m(r)$0 (Speeney et al., 2022). Four diagnostics are computed: accumulated GW cycles, waveform phase, accumulated SNR, and waveform mismatch via the noise-weighted inner product with the distinguishability threshold $m(r)$1. Crucially, each diagnostic is decomposed into a conservative geometric contribution ($m(r)$2, $m(r)$3) and a dissipative DF contribution ($m(r)$4, $m(r)$5).

Numerical results for Model I

The benchmark system has $m(r)$6, $m(r)$7, halo scale radius $m(r)$8 kpc, and densities $m(r)$9–$T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$0. Since $T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$1 throughout the inspiral region $T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$2, the small-$T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$3 expansions of Appendix A are highly accurate.

The central quantitative findings are:

Diagnostic	Key result
Cycle shift	Always negative; grows monotonically with $T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$4
One-cycle dephasing time	$T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$5 yr ($T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$6), $T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$7 yr ($T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$8), $T^\mu{}_\nu = \mathrm{diag}[-\rho, 0, P, P]$9 yr ($f$0)
Phase threshold ($f$1 rad)	Reached before SNR $f$2 for $f$3; at $f$4 yr for $f$5
Mismatch distinguishability	$f$6 yr for $f$7; $f$8 yr for $f$9; not reached within 4 yr for $m$0

All four diagnostics yield mutually consistent conclusions regarding detectability. The strongest and most consequential claim of the paper is the decomposition result: across cycles, phase, and mismatch, the observable signatures are overwhelmingly dominated by the conservative modification of the spacetime geometry, while relativistic dynamical friction remains many orders of magnitude smaller and negligible throughout the inspiral. This implies that the principal observable effect of an Einstein-cloud environment is encoded in the accumulated GW phase through the modified orbital dynamics, not in environmental dissipation—a conclusion that contrasts with much of the literature emphasizing DF-driven dephasing in dark-matter spikes.

Limitations and open questions

Several caveats qualify these results. First, the analysis assumes quasicircular inspirals; eccentric inspirals, which are generic for EMRIs, are left unexplored. Second, the leading-order quadrupole flux is used rather than self-force or black-hole perturbation-theory fluxes, though the authors note the balance-law framework would accommodate improved fluxes. Third, the local Chandrasekhar-type DF prescription neglects the global response of the halo and any halo modification induced by the inspiral—an approximation the authors acknowledge explicitly, although its impact is muted given the dominance of the conservative effect. Fourth, the mismatch criterion is fixed-parameter; a full Bayesian parameter-estimation analysis allowing intrinsic parameters to vary remains outstanding. Fifth, results are specialized to Model I and nonrotating SMBHs; generalization to Kerr backgrounds and other Einstein-cloud solutions is deferred. Finally, the detectability conclusions depend on assumed halo densities and the Coulomb logarithm ($m$1 in the timing estimates).

Conclusion

This paper provides a unified, fully relativistic treatment of how exact collisionless DM halos modify EMRI/IMRI gravitational-wave signals, demonstrating consistent detectability across four independent diagnostics within nominal LISA mission durations for halo densities $m$2 around $m$3–$m$4 primaries. Its most robust finding—that conservative spacetime deformation, rather than dynamical friction, dominates the observable signature—reframes how environmental effects should be modeled in EMRI waveform templates and establishes gravitational waves as a direct probe of the strong-field DM distribution complementary to electromagnetic observations.