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

Basic Information

• Paper ID: 2510.11635
• Title: Hints of Dark Matter Spikes in Low-mass X-ray Binaries: a critical assessment
• Authors: Francesca Scarcella, Bradley J. Kavanagh (Instituto de Física de Cantabria, UC-CSIC)
• Classification: hep-ph astro-ph.HE gr-qc
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.11635

Abstract

Three black hole low-mass X-ray binary systems (LMXBs) in the Milky Way exhibit orbital period decay rates that cannot be easily explained by standard mechanisms. Recent studies have claimed that the anomalous period decay in two of these systems can be explained by dynamical friction arising from extremely high dark matter densities surrounding the black holes. This paper provides a critical assessment of these claims through N-body simulations of binary systems embedded in dense dark matter "spikes." The study simulates the previously studied systems XTE J1118+480 and A0620-00, as well as Nova Muscae 1991, which is examined in this context for the first time. The simulations reveal that feedback effects from the dark matter distribution play a crucial role and exclude previously claimed shallow dark matter spikes. The research establishes lower bounds on the steepness parameter γ of the dark matter density profile required to explain the period decay in these LMXBs, requiring γ ≳ 2.15-2.20 for XTE J1118+480 and A0620-00, and γ ≳ 2.3 for Nova Muscae 1991.

Research Background and Motivation

Problem Statement

1. Anomalous orbital decay phenomenon: Three black hole low-mass X-ray binary systems in the Milky Way exhibit anomalously rapid orbital period decay, more than 100 times faster than predicted by standard mechanisms
2. Insufficiency of standard mechanisms: Standard angular momentum loss mechanisms—gravitational radiation, system mass loss, and magnetic braking—cannot explain the observed rapid decay rates
3. Controversy over dark matter explanations: Some studies propose that dark matter spikes around black holes might explain these anomalies through dynamical friction, but lack adequate consideration of feedback effects

Research Significance

1. Dark matter detection: The environment around black holes provides a unique laboratory for detecting dark matter effects, with dark matter densities potentially 20 orders of magnitude higher than in galactic halos
2. Theoretical verification: Testing the consistency between dark matter spike theory predictions and observational data
3. Multi-messenger astrophysics: Providing theoretical foundations for dark matter signals in future gravitational wave detections

Limitations of Existing Methods

1. Neglect of feedback effects: Previous studies based on the Chandrasekhar dynamical friction formula did not account for feedback from binary motion on the dark matter distribution
2. Semi-analytical approximations: Existing HaloFeedback formulations may overestimate feedback effects
3. Lack of comprehensive simulations: Absence of unified N-body simulation studies for all three systems

Core Contributions

1. First comprehensive N-body simulation study: Systematic N-body simulations of three LMXB systems using the NbodyIMRI code
2. Precise modeling of feedback effects: Accurately captures the dissipation of dark matter spikes and the impact of feedback effects on dynamical friction
3. Stringent constraints: Establishes more rigorous lower bounds on the steepness parameter γ of the dark matter density profile
4. First study of Nova Muscae 1991: Examines this system in a dark matter context for the first time
5. Exclusion of shallow spike models: Rules out previously proposed shallow dark matter spikes through energy balance analysis and numerical simulations

Methodology

Task Definition

Investigate the orbital evolution of binary systems embedded in dark matter spikes, using N-body simulations to determine the minimum dark matter density profile steepness required to explain the observed anomalous period decay.

Model Architecture

Dark Matter Spike Model

Dark matter density distribution follows a power-law profile:

ρ_DM(r) = {
  ρ_sp(r) = ρ_bkg(r/r_sp)^(-γ), r ≤ r_sp
  ρ_bkg, r ≥ r_sp
}

where γ is the density profile steepness parameter and r_sp is the spike radius.

N-body Simulation Framework

Simulations performed using the NbodyIMRI code:

1. Particle count: 20,000 dark matter particles per simulation (25,000 for γ=2.45)
2. Time step: ~10,000 steps per orbit, employing hierarchical time-stepping for efficiency
3. Softening: Black hole softening radius r_soft = 3×10³r_isco, with physical radius softening for the star
4. Simulation duration: 6,000 orbital periods (approximately several years of physical time)

Dynamical Friction Calculation

Chandrasekhar dynamical friction formula:

F_DF = -4πG²M²_⋆ρξ(v)ln Λ/v²

Corresponding period decay rate:

Ṗ = -6G_N(M_⋆/M_BH)P²ρ(r_orb)ξ(v)ln Λ

Technical Innovations

1. Complete feedback modeling: N-body simulations naturally incorporate acceleration and displacement of dark matter particles, accurately capturing spike dissipation
2. Multi-system comparative study: Simultaneous investigation of three systems with different mass ratios, validating result generality
3. Energy balance analysis: Estimates disruption timescales by comparing spike binding energy with energy injection rates
4. Statistical error handling: Robust estimation methods using 68% confidence intervals to handle outliers

Experimental Setup

Research Objects

Three black hole low-mass X-ray binary systems:

Parameter Space

• Dark matter spike steepness: γ ∈ 1.9, 2.45, step size 0.05
• Per γ value: 16 independent realizations to quantify statistical errors
• Background density: Based on Milky Way NFW profile, ρ_bkg ~ 0.3-0.4 GeV cm⁻³

Evaluation Metrics

1. Period decay rate: Extract Ṗ_γ from simulated orbital evolution
2. Density dissipation: Monitor dark matter density changes over time
3. Confidence intervals: Determine parameter constraints using 68% confidence level

Experimental Results

Main Results

Dark Matter Spike Dissipation

• Rapid initial phase: Density dissipates by 2-3 orders of magnitude within the first few hundred orbits
• Quasi-steady state: Relatively stable configuration achieved after several thousand orbits
• Mass ratio dependence: Systems with larger mass ratios (Nova Muscae 1991) exhibit stronger dissipation

Period Decay Rate Constraints

Lower bounds on γ for the three systems:

• XTE J1118+480: γ ≳ 2.19
• A0620-00: γ ≳ 2.15
• Nova Muscae 1991: γ ≳ 2.32

Density Requirements

Accounting for feedback effects, required initial dark matter densities are approximately 3 orders of magnitude higher than naive Chandrasekhar analysis predicts.

Long-term Evolution Tests

Extended simulations to 15,000 orbits show:

• Ṗ values continue to decline slowly
• Best-fit γ values may increase by ~0.05
• Systems have not yet reached complete steady state

Disruption Timescale Analysis

Spike disruption timescales estimated through energy balance:

• γ < 2: T_disr < 20 years (observational timescale)
• γ = 2.5: T_disr ~ 10⁵-10⁷ years

Related Work

Dark Matter Spike Theory

1. Adiabatic growth: r⁻⁷/³ profile proposed by Gondolo & Silk (1999)
2. Primordial black holes: Expected steeper γ = 9/4 profiles
3. Realistic modeling: Improved descriptions from Bertone et al. (2025)

Binary Orbital Evolution

1. Standard mechanisms: Gravitational radiation, magnetic braking, mass loss
2. Circumbinary disks: Resonant torque mechanisms
3. Enhanced magnetic braking: CARB prescription

Numerical Methods

1. HaloFeedback: Semi-analytical feedback prescription
2. N-body simulations: Recent IMRI/EMRI studies
3. Gravitational wave detection: Future multi-messenger observation prospects

Conclusions and Discussion

Main Conclusions

1. Feedback effects are crucial: Dark matter distribution feedback significantly affects dynamical friction efficiency
2. Shallow spikes excluded: Shallow spikes with γ < 2.15-2.3 cannot sustain observed decay rates
3. Consistency check: The γ values required by all three systems are remarkably consistent, approaching theoretical expectations
4. Density constraints: Extremely high dark matter densities (~10¹⁵-10²⁰ GeV cm⁻³) required at orbital radii

Limitations

1. Timescale constraints: Simulation duration (several years) far shorter than system age (Gyr)
2. Formation mechanism: Formation mechanism of dark matter spikes around stellar-mass black holes remains unclear
3. Long-term evolution: Extrapolation to longer timescales needed to determine final steady state
4. Uncertainties: Systematic errors in observational parameters not fully considered

Future Directions

1. Extended simulations: Extend to longer timescales to study ultimate evolution
2. Formation mechanisms: Investigate possible formation pathways for dark matter spikes in LMXBs
3. Primordial black holes: Detailed investigation of feasibility in PBH-spike scenarios
4. Multi-messenger observations: Joint analysis combining gravitational wave and electromagnetic observations

In-Depth Evaluation

Strengths

1. Rigorous methodology: First systematic study of this problem using complete N-body simulations
2. Advanced techniques: NbodyIMRI code specifically designed for such problems
3. Reliable results: Multi-system consistency and statistical treatment enhance credibility
4. Comprehensive analysis: Multi-level analysis from energy balance to numerical simulations

Weaknesses

1. Timescale challenges: Enormous discrepancy between simulation timescale and actual system age
2. Theoretical difficulties: Theoretical challenges in dark matter spike formation mechanisms
3. Parameter dependence: Result dependence on softening parameters and resolution requires further verification
4. Alternative explanations: Discussion of purely astrophysical explanations relatively brief

Impact

1. Methodological contribution: Establishes new standards for binary evolution studies in dark matter environments
2. Theoretical constraints: Provides important observational constraints for dark matter spike theory
3. Future observations: Provides theoretical foundation for dark matter signals in gravitational wave detection
4. Interdisciplinary significance: Connects particle physics, astrophysics, and gravitational wave astronomy

Applicable Scenarios

1. Similar binary systems: Methods generalizable to other anomalously evolving tight binaries
2. Gravitational wave sources: Provides theoretical guidance for dark matter searches with future detectors like LISA
3. Dark matter research: Opens new pathways for studying dark matter properties in extreme environments
4. Numerical methods: NbodyIMRI code applicable to broader range of related problems

References

The paper contains 84 references spanning dark matter theory, binary evolution, numerical methods, and observational astrophysics, reflecting the interdisciplinary nature and deep theoretical foundation of this research.

Overall Assessment: This is a high-quality research paper of significant importance in dark matter astrophysics. Through rigorous numerical simulations, the authors provide a critical assessment of a controversial dark matter detection claim, not only correcting deficiencies in previous research but also providing important theoretical constraints and methodological contributions for future development in this field. Despite challenges such as timescale extrapolation, the scientific value and technical rigor of the research are both high.