The Depletion of Collisionless Dark Matter Spikes

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Cosmic Expectation vs. Reality

Imagine you are looking for a hidden treasure (Dark Matter) around a giant, invisible whirlpool in space (a Supermassive Black Hole).

For decades, scientists believed that as a black hole slowly grew, it would act like a vacuum cleaner, sucking in the surrounding dark matter and creating a super-dense "spike" right next to it. Think of this like a snowdrift piling up high against a house wall during a blizzard. If this "snowdrift" existed, it would be so dense that when smaller objects (like stars or smaller black holes) spiraled into the big black hole, they would drag against this thick dark matter, creating a unique "fingerprint" in the gravitational waves they emit.

The new discovery in this paper is: That snowdrift probably doesn't exist. Or at least, it gets melted away long before we can see it. The authors show that the environment around these black holes is much more chaotic and "messy" than previously thought, and it actively destroys these dense spikes.

The Two "Thieves" Stealing the Dark Matter

The paper identifies two main mechanisms that strip away the dark matter, acting like two different thieves stealing from a bank vault.

Thief #1: The Heavyweights (Mass Segregation)

The Old Idea: Scientists used to think the stars around a black hole were all roughly the same size, like a crowd of people all wearing the same size shoes. In this calm crowd, the dark matter would slowly settle into a dense pile.

The New Reality: The authors realized the crowd isn't uniform. It's a mix of tiny mice (small stars) and heavy elephants (massive black holes and giant stars).

The Analogy: Imagine a crowded dance floor where heavy dancers (massive stars) naturally sink to the center because they are heavy, while the light dancers (small stars) get pushed to the edges. This is called Mass Segregation.
The Result: Because the heavy stars crowd the center, they bump into each other and the dark matter much more frequently. It's like a mosh pit. These constant, chaotic bumps heat up the dark matter, kicking it out of the center.
The Outcome: Instead of a steep, dense spike, the dark matter settles into a much flatter, lower-density cloud (called a Bahcall-Wolf profile) in less than a billion years. The "snowdrift" has been flattened into a thin layer of dust.

Thief #2: The Slingshot (EMRIs)

The Scenario: Deep inside the center, where the dark matter is supposed to be thickest, there are smaller black holes (stellar-mass black holes) spiraling into the giant one. These are called EMRIs (Extreme Mass-Ratio Inspirals).

The Mechanism:

The Analogy: Imagine the dark matter particles are like ping-pong balls floating in a room. The spiraling black hole is a giant, fast-moving bowling ball. As the bowling ball swings through the room, it doesn't just push the ping-pong balls; it grabs them with gravity and flings them out of the room like a slingshot.
The Result: Every time a small black hole spirals in, it ejects a chunk of the dark matter. Since the dark matter is "collisionless" (the ping-pong balls don't bump into each other to reform the pile), once they are kicked out, they are gone forever.
The Outcome: Over billions of years, the repeated "slingshot" effect of these spiraling black holes empties the inner core of the dark matter spike completely.

Why This Matters for LISA (The Space Telescope)

The European Space Agency is building a space telescope called LISA (Laser Interferometer Space Antenna) to listen to gravitational waves. Scientists were very excited because they thought: "If we find a black hole with a dense dark matter spike, the gravitational waves will sound different, and we can prove what dark matter is!"

The Bad News:
This paper shows that for most black holes, the dark matter spike has already been destroyed by the two "thieves" mentioned above.

If the spike is gone, the gravitational waves will look just like they would in a vacuum (empty space).
The "fingerprint" of dark matter will be too faint to detect.

The Good News (Sort of):
This actually helps us understand the universe better. It tells us that we shouldn't waste time looking for these spikes in the places we thought were most likely. Instead, we need to look at:

Very young black holes (at very high redshifts, before the "thieves" had time to work).
Very specific types of dark matter that might stick together (self-interacting) and refill the pile.

Summary in a Nutshell

The Dream: Black holes should be surrounded by super-dense spikes of dark matter that leave a clear signal in gravitational waves.
The Reality Check:
- Heavy stars mix up the crowd, kicking dark matter out of the center (like a mosh pit).
- Spiraling black holes act as slingshots, flinging the remaining dark matter into deep space.
The Conclusion: The "perfect" dark matter spikes are likely rare or non-existent in the nearby universe. This makes it much harder for the LISA telescope to detect dark matter using gravitational waves, narrowing down the search to very specific, rare scenarios.

Essentially, the universe is too messy and chaotic to let these perfect dark matter piles survive long enough for us to easily find them.

1. Problem Statement

The existence of dense "spikes" of collisionless dark matter (DM) surrounding massive black holes (MBHs) has long been a theoretical prediction, most notably in the Gondolo–Silk model. In this paradigm, the adiabatic growth of an MBH within a pre-existing DM cusp steepens the density profile to $\rho \propto r^{-7/3}$ . Such spikes are considered prime candidates for detection via:

Gravitational Waves (GWs): Inducing measurable dephasing in Extreme Mass-Ratio Inspirals (EMRIs) detectable by LISA.
Indirect Detection: Enhancing annihilation signals (gamma-rays/neutrinos) for self-annihilating DM.

However, previous studies often relied on idealized assumptions: single-mass stellar populations and simplified relaxation models. This paper challenges the robustness of these spikes by addressing two critical, realistic factors often neglected:

Multi-mass stellar dynamics: Realistic nuclear star clusters (NSCs) contain a spectrum of stellar masses, leading to mass segregation, which accelerates relaxation.
Strong encounters in the loss cone: The innermost regions are dominated by discrete, strong interactions between DM particles and compact objects (stellar-mass black holes, sBHs) spiraling into the MBH (EMRIs), rather than the weak, cumulative encounters described by standard Fokker–Planck (FP) theory.

The authors ask: Do realistic stellar dynamics and EMRI interactions deplete DM spikes before they can be observed by LISA?

2. Methodology

The authors employ a two-regime approach to model the evolution of the DM spike within the MBH's sphere of influence (SoI):

A. The Stellar Sphere (Outer Region): Orbit-Averaged Fokker–Planck Simulations

Code: They use PHASEFLOW (part of the AGAMA package), a 1D orbit-averaged, isotropic Fokker–Planck solver.
Stellar Models: They compare two scenarios:
1. Single-effective-mass model: A traditional approach using a single "effective" stellar mass.
2. Multi-mass model: A realistic distribution of stellar masses (0.1 to 100 $M_\odot$ ) following the Kroupa initial mass function, including stellar-mass black holes.
Physics: The simulations track the co-evolution of stars and DM, calculating two-body relaxation times and dynamical friction. They specifically look for the transition to the steady-state Bahcall–Wolf (BW) profile ( $\rho \propto r^{-3/2}$ ).
Initial Conditions: The DM spike starts with a Gondolo–Silk profile ( $\rho \propto r^{-7/3}$ ). The stellar component starts with a Hernquist profile.

B. The Loss Cone (Inner Region): Post-Newtonian 3-Body Simulations

Context: Inside the loss cone (radii $\lesssim 10^2 R_{BH}$ ), the FP approximation breaks down. Dynamics are governed by strong, discrete encounters.
Model: They simulate three-body interactions (MBH + sBH/EMRI + DM particle) using the multistar code, which includes 3.5 Post-Newtonian (PN) corrections.
Mechanism: They calculate the probability of DM particles being ejected via gravitational slingshots during close encounters with EMRIs.
Scaling: They establish that the ejection fraction scales as $q^2$ (where $q = M_{EMRI}/M_{BH}$ ).
Integration: They integrate these ejection rates over the lifetime of the MBH (up to $z=3$ , or 2.14 Gyr) using EMRI distributions derived from their FP simulations.

3. Key Contributions

Acceleration of Relaxation via Mass Segregation: The paper demonstrates that multi-mass stellar populations relax significantly faster than single-mass models. Mass segregation drives heavy stars inward, increasing central density and scattering rates, which in turn heats the DM component.
Irreversible Depletion Mechanism: The authors identify that EMRIs eject DM particles to low binding energies. Because collisionless DM cannot efficiently repopulate these depleted phase-space regions (due to the extremely long DM-DM relaxation time), the depletion is irreversible.
Quantification of Depletion Rates: They provide a quantitative framework linking EMRI rates, MBH masses, and the remaining fraction of DM density, showing that even conservative EMRI rates lead to massive depletion.

4. Key Results

A. Stellar Sphere Evolution (Outer Spike)

Time-scales: In multi-mass models, the relaxation time is shorter by a factor of 10–100 compared to single-mass models.
Profile Evolution: While single-mass models maintain the steep $\rho \propto r^{-7/3}$ spike for several Gyr, multi-mass models drive the DM density to the much flatter Bahcall–Wolf profile ( $\rho \propto r^{-3/2}$ ) within $\lesssim 1$ Gyr.
Density Reduction: By $t \sim 1$ Gyr, the central DM density in multi-mass models drops by 3 to 5 orders of magnitude compared to the initial Gondolo–Silk prediction.

B. Loss Cone Evolution (Inner Spike)

Slingshot Ejection: Repeated EMRI passages eject DM particles. The depletion is exponential with respect to the total accumulated mass of EMRIs.
Critical Mass Threshold: Systems with $M_{BH} \lesssim 1.7 \times 10^6 M_\odot$ relax within the 2.14 Gyr window, leading to a specific EMRI distribution that maximizes depletion.
Depletion Factors:
- For typical EMRI rates ( $\sim 100-300$ Gyr $^{-1}$ ), the DM density within the LISA band ( $r < r_{p, LISA}$ ) is reduced by several orders of magnitude for $M_{BH} \lesssim 10^6 M_\odot$ .
- Even for low rates ( $\sim 3-10$ Gyr $^{-1}$ ), MBHs with $M_{BH} \lesssim 10^5 M_\odot$ experience near-complete depletion.
- A single IMRI event around an Intermediate Mass Black Hole (IMBH) of $10^3-10^4 M_\odot$ could clear the spike entirely.

C. Impact on GW Dephasing

Dephasing Suppression: Existing predictions for DM-induced GW dephasing range from $10^2$ to $10^6$ radians.
Detectability: A density reduction factor of $f_{rem} \lesssim 10^{-6}$ (which the authors find likely) suppresses the accumulated phase shift to $< 1$ radian.
Conclusion: This falls below the detection threshold for LISA, rendering the DM signal indistinguishable from vacuum evolution for most viable MBH masses.

5. Significance and Implications

Redefining the Search Space: The paper drastically narrows the parameter space where collisionless DM spikes can be detected. The "sweet spot" for LISA (MBHs of $10^5-10^6 M_\odot$ ) is likely devoid of detectable spikes due to EMRI-driven depletion.
Re-evaluating Indirect Detection: The depletion of central densities implies that non-detection of gamma-ray signals from the Galactic Center may not rule out DM annihilation, but rather confirm that spikes do not survive in realistic environments.
Modeling Necessity: Future studies of DM signatures must account for realistic stellar mass functions and the cumulative effect of EMRIs. Ignoring these leads to a severe overestimation of DM density and detectability.
Exceptions: Detectable spikes might only survive at very high redshifts ( $z \gg 3$ ), in systems without nuclear clusters, or if DM is self-interacting (SIDM), where collisions could replenish the depleted phase space.

In summary, Sharpe et al. argue that the canonical Gondolo–Silk spike is a transient feature that is rapidly eroded by the very stellar dynamics and compact object populations that define realistic galactic nuclei, effectively closing the window for detecting collisionless DM spikes via LISA for the foreseeable future.