Hints of Dark Matter Spikes in Low-mass X-ray Binaries: a critical assessment

Through $N$-body simulations, this study critically assesses claims that dark matter spikes explain anomalous period decays in three low-mass X-ray binaries, ruling out previously proposed shallow profiles and establishing that significantly steeper density slopes ($\gamma \gtrsim 2.15$) are required to account for the observed orbital evolution.

The Gist

Imagine you are watching a dance between two partners: a massive, invisible giant (a Black Hole) and a tiny, glowing star. They are locked in a tight embrace, spinning around each other so fast that they complete a full circle in just a few hours.

According to the rules of physics we know, this dance should slow down very gradually over millions of years, like a spinning top losing energy. However, astronomers have noticed something strange with three specific pairs of these dancers. They are spiraling inward and speeding up their dance 100 times faster than anyone expected. It's as if the spinning top suddenly started spinning out of control.

For a long time, scientists tried to explain this with standard "brakes" like magnetic fields or gas friction, but nothing fit the data perfectly.

The New Theory: The "Dark Matter Crowd"

Recently, some researchers proposed a wild idea: What if these black holes are surrounded by a super-dense crowd of Dark Matter?

Think of Dark Matter as an invisible, ghostly gas that fills the universe. Usually, it's very thin, like fog in a vast field. But near a black hole, this gas might pile up into a sharp, steep mountain peak called a "spike."

If the black hole is moving through this dense spike, the invisible particles would act like thick molasses or a crowded dance floor. As the star orbits, it would constantly bump into these invisible particles. This creates a drag force (called Dynamical Friction) that sucks energy out of the orbit, causing the star to spiral inward much faster.

The Problem: The "Feedback Loop"

The authors of this paper, Francesca Scarcella and Bradley Kavanagh, decided to test this idea. But they didn't just do math on a piece of paper; they built a massive, virtual universe using a supercomputer to simulate these dances.

Here is the twist they discovered: The dancers change the crowd.

In previous studies, scientists assumed the Dark Matter crowd was static and unchanging. But in reality, as the star and black hole spin through the Dark Matter, they act like a blender in a smoothie. Their movement kicks the Dark Matter particles out of the way, clearing a hole in the "molasses."

• The Analogy: Imagine running through a dense forest. At first, the trees (Dark Matter) are thick, and you get slowed down a lot. But as you run, you knock the trees over and clear a path. Suddenly, the forest isn't dense anymore; it's a cleared path. The drag force disappears because you've destroyed the very thing that was slowing you down.

The authors found that if the Dark Matter spike isn't incredibly steep and dense to begin with, the binary system would clear it out so quickly that the "drag" would vanish before we could even observe the fast orbital decay.

The Results: How Dense Must the Crowd Be?

The team ran thousands of simulations to see what kind of Dark Matter spike could survive long enough to explain the observations.

1. Shallow Spikes Fail: If the Dark Matter density rises slowly (a "shallow" spike), the system clears it out in a few years. The drag disappears, and the orbit doesn't decay fast enough. This rules out the theories that suggested a gentle slope.
2. Steep Spikes are Required: To explain the fast orbital decay, the Dark Matter spike must be extremely steep (like a sheer cliff face rather than a gentle hill). The density must be so high that even after the system clears out some of it, there's still enough "molasses" left to keep the drag strong.

They calculated that the "steepness" of this Dark Matter mountain needs to be much higher than previously thought.

The Big "But" (The Catch)

Even if the math works, there is a huge problem with the story: How did this crowd get there?

• Standard Black Holes: The black holes in these systems are likely "normal" black holes formed from dying stars. The standard story of how they form involves violent explosions and chaotic movements. It is very hard to imagine how a delicate, dense cloud of Dark Matter could survive that violence. It's like trying to keep a house of cards intact while a tornado is tearing through the room.
• Primordial Black Holes: The only scenario where this makes sense is if these black holes are "Primordial"—meaning they were born in the very first moments of the Big Bang, long before stars existed. These ancient black holes are expected to naturally come with these dense Dark Matter spikes.

The Conclusion

The paper is a critical reality check. It says:

"The idea that Dark Matter is causing these fast orbital decays is not impossible, but it requires the Dark Matter to be much denser and steeper than we thought. Furthermore, it implies these black holes might be ancient, primordial relics from the dawn of time, rather than normal stellar black holes."

If future observations confirm these extreme densities, it would be a massive discovery, potentially proving that some black holes are primordial and that Dark Matter behaves in very specific, extreme ways near them. But for now, the "Dark Matter explanation" is still a hypothesis that needs to survive the test of time and more simulations.

Technical Summary

Here is a detailed technical summary of the paper "Hints of Dark Matter Spikes in Low-mass X-ray Binaries: a critical assessment" by Francesca Scarcella and Bradley J. Kavanagh.

1. Problem Statement

Three Galactic black-hole low-mass X-ray binaries (LMXBs)—XTE J1118+480 (System A), A0620–00 (System B), and Nova Muscae 1991 (System C)—exhibit orbital period decay rates ($\dot{P}$) that are 100 times faster than predicted by standard mechanisms (gravitational wave emission, magnetic braking, and mass loss).

Recent studies (e.g., Chan & Lee) proposed that dynamical friction (DF) caused by dense Dark Matter (DM) spikes surrounding the black holes could explain these anomalies. However, these previous claims relied on static DM density profiles and neglected the feedback effect: the energy extracted from the binary by dynamical friction is injected into the DM spike, depleting the density near the orbit and reducing the friction force. The authors aim to critically assess the DM spike hypothesis by performing $N$-body simulations that self-consistently account for this feedback.

2. Methodology

The authors utilized NbodyIMRI, a dedicated $N$-body code designed to study binaries embedded in DM spikes.

• Simulation Setup:

  • Systems: Simulations were run for the three specific LMXBs using observed parameters (BH mass, stellar mass, orbital period, and radius) from Table I.
  • DM Profile: The initial DM density was modeled as a power-law spike $\rho_{DM} \propto r^{-\gamma}$ extending from a truncation radius $r_{sp}$ (defined by the BH's gravitational influence) out to a cutoff radius.
  • Parameters: The spike slope $\gamma$ was varied between 1.9 and 2.45. For each $\gamma$, 16 different realizations (random initializations) were generated.
  • Physics: The code tracks the Newtonian evolution of the binary and DM pseudo-particles. It includes gravitational forces between the binary components and between the binary and DM particles. Crucially, it neglects DM-DM pairwise interactions (justified by the low total DM mass relative to the binary) and excludes gravitational wave emission to isolate the DF effect.
  • Softening: To prevent numerical singularities, the BH and star potentials were softened. The BH softening radius was set to $3 \times 10^3 r_{isco}$, and the star's softening radius matched its physical radius.
  • Duration: Simulations ran for 6,000 orbits (approx. a few years). A subset of simulations for System A was extended to 15,000 orbits to check for long-term stability.

• Analysis:

  • The orbital period decay rate $\dot{P}$ was extracted from the simulation data after the initial transient phase (where the spike is rapidly depleted) settled into a quasi-stationary state.
  • Confidence intervals (68% CL) were established by analyzing the spread of results across the 16 realizations, filtering out outliers caused by rare close encounters.

3. Key Contributions

• Inclusion of Feedback: This is the first study to use full $N$-body simulations to quantify the depletion of DM spikes due to the binary's motion in the context of LMXBs. Previous analytic estimates (e.g., the HaloFeedback formalism) were found to potentially overestimate the resilience of the spike.
• Re-evaluation of Shallow Spikes: The study demonstrates that "shallow" spikes (low $\gamma$) cannot sustain the observed period decays because they are disrupted on timescales shorter than the observation window ($\sim 20$ years).
• New Constraints on System C: The paper provides the first DM spike analysis for Nova Muscae 1991, a system with a higher mass ratio ($q \approx 0.08$) compared to the others, leading to stronger feedback effects.
• Distinction from PBH Scenarios: The authors discuss the implications for Primordial Black Holes (PBHs), noting that while standard astrophysical BH formation scenarios struggle to host such spikes, PBHs naturally possess steep spikes ($\gamma \approx 2.25$) that remain viable candidates.

4. Key Results

A. Spike Depletion and Feedback

• Rapid Initial Depletion: Simulations show a rapid depletion of DM density around the orbit during the first few hundred orbits ($\sim 1$ year).
• Quasi-Stationary State: After the initial phase, the density stabilizes at a level roughly two orders of magnitude lower than the initial profile.
• Mass Ratio Dependence: The depletion is most severe in System C (highest mass ratio), where density drops by nearly three orders of magnitude. This confirms that feedback is more pronounced for larger mass ratios.

B. Constraints on Spike Slope ($\gamma$)

The simulations revealed that the period decay rate $|\dot{P}|$ increases with the steepness of the spike ($\gamma$). By comparing simulated $\dot{P}$ values with observational constraints, the authors established lower limits on $\gamma$ required to explain the anomalies:

• System A (XTE J1118+480): $\gamma \gtrsim 2.19$
• System B (A0620–00): $\gamma \gtrsim 2.15$
• System C (Nova Muscae 1991): $\gamma \gtrsim 2.32$

Crucial Finding: The previously proposed shallow spikes with $\gamma \approx 1.71 - 1.85$ (suggested by Chan & Lee) are ruled out. These profiles are depleted too quickly to sustain the observed decay rates over the 20-year observation period.

C. Required DM Densities

To match the observed $\dot{P}$ after accounting for feedback, the initial DM density at the orbital radius must be significantly higher than naive estimates:

• The required densities are roughly three orders of magnitude larger than those inferred from the static Chandrasekhar formula (Eq. 2) which ignores feedback.
• For System C, despite its much larger observed decay rate, the required DM density is comparable to System A due to the stronger depletion effects.

D. Long-term Evolution

Extended simulations (15,000 orbits) showed that $\dot{P}$ continues to decrease slowly as the spike is further depleted. This implies that the lower limits on $\gamma$ derived from the 6,000-orbit runs are conservative; the true required $\gamma$ might be even steeper if the system has been evolving for millions of years.

5. Significance and Conclusion

• Critical Assessment: The paper does not confirm the existence of DM spikes but establishes that only very steep spikes ($\gamma > 2.15$) could theoretically explain the data.
• Formation Mechanism Implications:
  • Astrophysical BHs: Standard formation scenarios (adiabatic growth in a halo) typically predict spikes that are disrupted by the binary's formation or evolution. The required steepness and density challenge these standard models.
  • Primordial Black Holes (PBHs): PBHs are expected to form with steep spikes ($\gamma = 2.25$) and high densities ($\sim 10^{-6}$ g/cm$^3$). The results are consistent with a PBH origin for these systems, though long-term depletion over Gyr timescales requires further study.
• Alternative Explanations: The authors note that purely astrophysical mechanisms (e.g., circumbinary disks or convection-boosted magnetic braking) also struggle to consistently explain all three systems simultaneously, leaving the DM spike hypothesis (specifically for PBHs) as a viable, though unproven, alternative.

Final Conclusion: The DM spike hypothesis cannot be ruled out, but it requires steep density profiles ($\gamma \gtrsim 2.2$) and implies either a Primordial Black Hole origin or a non-standard formation history for these LMXBs. Future work must address the long-term evolution of these spikes over Gyr timescales to definitively confirm or exclude this interpretation.