We investigate the possibility that cosmic-ray electron cooling through dark matter-electron scatterings contributes to the low radiative efficiency observed in radio-loud galaxies such as M87. Light dark matter can scatter efficiently off electrons in M87, lowering the observed bolometric luminosity compared to astrophysical expectations. This consideration allows us to probe previously unexplored regions of the parameter space of dark matter-electron interactions. We further model the cosmic-ray electron distribution by numerically solving a diffusion equation along the jet and find that efficient dark matter-electron interactions can induce a flattening of the spectral index at different distances from the central supermassive black hole, in better alignment with radio observations from M87.
Dark Matter-Electron Interactions Alter the Luminosity and Spectral Index of M87
- Paper ID: 2510.12877
- Title: Dark Matter-Electron Interactions Alter the Luminosity and Spectral Index of M87
- Authors: Abdelaziz Hussein (MIT), Gonzalo Herrera (MIT, Harvard, Virginia Tech)
- Classification: hep-ph astro-ph.CO astro-ph.HE
- Publishing Institution: FERMILAB-PUB-25-0729-T
- Paper Link: https://arxiv.org/abs/2510.12877v1
This paper investigates whether cosmic ray electron cooling induced by dark matter-electron scattering can explain the low radiative efficiency phenomenon observed in radio-active galaxies (such as M87). Light dark matter can effectively scatter with electrons in M87, reducing the observed thermal luminosity relative to astrophysical predictions. This consideration enables exploration of previously unexplored regions in the dark matter-electron interaction parameter space. The study further models the cosmic ray electron distribution by numerically solving the diffusion equation along the jet, finding that effective dark matter-electron interactions can cause spectral index flattening at different distances from the central supermassive black hole, achieving better consistency with M87's radio observations.
- Silent Black Hole Problem: The observed thermal luminosity is many orders of magnitude lower than the Eddington luminosity, an open question in galactic astrophysics
- Limitations of Traditional Models: Existing accretion models (such as ADAF) and GRMHD simulations face difficulties in explaining observed fluxes and spectral index distributions
- Model Simplification Assumptions: GRMHD simulations typically assume thermal Maxwell-Juttner distributions, neglecting non-thermal effects such as magnetic reconnection and turbulent dissipation
- Explore the potential impact of dark matter particles on black hole observations
- Utilize the dark matter spike structure near supermassive black holes
- Develop new dark matter detection methods, particularly for sub-GeV mass regions
- Address velocity suppression issues in direct detection experiments (relativistic electrons vs. non-relativistic dark matter)
- Propose a New Dark Matter Detection Mechanism: Utilize cosmic ray electron cooling induced by dark matter-electron scattering in M87 to constrain dark matter parameters
- Establish Luminosity-Cooling Timescale Relationship: Establish quantitative constraint relationships by correlating observed low thermal luminosity with dark matter cooling effects
- Develop Spectral Index Modeling Method: Numerically solve the diffusion equation to predict the impact of dark matter interactions on jet spectral indices
- Obtain Competitive Constraints: Surpass existing direct detection experimental constraints in certain parameter regions
- Explore Unconstrained Parameter Space: Particularly for velocity-suppressed regions of Majorana dark matter
Employ an extended Standard Model containing an additional U(1)' gauge symmetry:
L=−41FμνFμν−41Fμν′F′μν−2ϵFμνF′μν+2mA′2Aμ′A′μ+Aμ′gχχˉγμχ
where A′ is the dark photon, ϵ is the kinetic mixing parameter, and gχ is the dark matter coupling strength.
The differential cross section for dark matter-electron scattering is:
dTχdσ=4π(Ee2−me2)(mA′2+2mχTχ)2ϵ2e2gχ2×[2Ee2mχ−(me2+mχ(2Ee+mχ))Tχ+mχTχ2]
Synchrotron Cooling Timescale:
τsyn=[3Ee4cσT8πB2γeβe2]−1
Dark Matter-Electron Scattering Cooling Timescale:
τχe=[−Ee1(dtdEe)χe]−1
where:
(dtdEe)χe=−mχρDM∫0TχmaxdTχTχdTχdσχe
Based on the principle that dark matter cooling cannot reduce luminosity below observed values:
log(LeddLbol)=−τχeτsync≈−9
Parameters are considered excluded when τsync/τχe≥10.
Solve the diffusion equation:
dτdN(γ,τ)=Q(γ,τ)−3badiN(γ,τ)−∂γ∂b(γ,τ)N(γ,τ)
where the total energy loss rate is:
b=bSM+bDM=−badiγ−bsyncγ2+me1(dtdE)χe
Employ a dark matter spike formed by adiabatic growth of the supermassive black hole:
ρDM(r)=ρRgγ(r)(rRsp)γsp
For M87: Rsp≈4×105Rs, γsp=7/3
Consider two representative values:
- Conservative Estimate: ⟨BVLBI⟩=7G (based on average VLBI measurements)
- Aggressive Estimate: BMAGIC=0.003G (based on MAGIC telescope SED fitting)
- Thermal Luminosity: Lbol∼1042 erg/s
- Eddington Luminosity: Ledd≈8.45×1047 erg/s
- Spectral Index Observations: From KaVA and VLBA observations at 22 GHz and 43 GHz
Dark Photon Parameter Space Constraints:
- For mχ=15 eV: Constrain ϵ≳10−12 in the mA′∼10−6 eV region
- For mχ=15 keV: Stricter constraints covering larger parameter space
DM-Electron Scattering Cross Section Constraints:
- Dirac Dark Matter: Surpass direct detection constraints when mχ≤6 MeV
- Majorana Dark Matter: More significant constraints due to velocity suppression
- Using BMAGIC, constraints are 2-4 orders of magnitude stronger than future direct detection experiments
Dark matter interactions lead to spectral index flattening:
- mχ=100 MeV, mA′=150 GeV: No significant effect
- mχ=10 MeV, mA′=30 GeV: Predictions consistent with observations within 1σ range
- Lighter dark matter masses lead to over-flattening, inconsistent with observations
- BBN Constraints: mχ≳10 MeV (from effective relativistic degrees of freedom measurements)
- Direct Detection: DAMIC-M (mχ≲20 MeV), PandaX-4T (20-200 MeV), CRESST-III (>200 MeV)
- Thermal Relic Abundance: Covers portions of thermally-produced dark matter parameter space
- Gondolo & Silk (1999): First proposed the concept of dark matter spikes near supermassive black holes
- Subsequent studies considered effects of self-annihilation, stellar heating, and mergers on spike structure
- Recent studies (Herrera & Murase 2024, etc.): Utilize cosmic ray cooling in AGN to constrain dark matter
- This work extends to luminosity constraints and spectral index modeling
- Recent breakthroughs by DAMIC-M and PandaX-4T in sub-GeV dark matter detection
- This method circumvents velocity suppression issues in laboratory experiments
- New Constraint Method: M87's low luminosity can serve as a new tool for constraining dark matter-electron interactions
- Parameter Space Exploration: Provides stronger constraints than direct detection in certain regions
- Spectral Index Diagnostics: First proposal to utilize spectral index observations at different distances to constrain dark matter interactions
- Model Verification: Certain dark matter parameter configurations better explain M87's radio observations
- Magnetic Field Uncertainty: Constraint strength highly depends on magnetic field strength assumptions (varying by 3-4 orders of magnitude)
- Dark Matter Density Distribution: Spike structure may be disrupted by astrophysical processes
- Standard Model Process Understanding: Spectral index constraints depend on accurate modeling of electron acceleration mechanisms
- Single Source Limitation: Based solely on M87; extension to larger AGN samples needed
- AGN Sample Extension: Apply framework to more radio-active galaxies
- Multi-Wavelength Observations: Combine simultaneous observations at different wavelengths to distinguish dark matter effects from conventional astrophysical processes
- Density Distribution Determination: Develop particle-independent methods to determine dark matter density distribution
- Theoretical Improvements: Improve understanding of standard model processes through particle-in-cell simulations
- High Originality: First to connect AGN luminosity suppression with dark matter interactions, opening new detection pathways
- Complete Methodology: Forms comprehensive analytical framework from theoretical modeling to numerical computation
- Practical Value: Provides currently strongest constraints in certain parameter regions, particularly for Majorana dark matter
- Observation-Oriented: Closely integrates multi-frequency observational data from M87
- Technical Rigor: Considers relativistic effects, magnetic field distribution, dark matter density distribution and other critical factors
- Large Systematic Uncertainties: Uncertainties in magnetic field and dark matter density lead to wide constraint ranges
- Model Dependence: Results strongly depend on specific dark matter models and astrophysical assumptions
- Verification Difficulty: Hard to independently verify dark matter cooling hypothesis; model degeneracies exist
- Limited Applicability: Method primarily applicable to AGN with strong jets, limiting universality
- Interdisciplinary Bridging: Connects particle physics and astrophysics, promoting cross-disciplinary research
- Methodological Contribution: Provides new perspectives for dark matter indirect detection
- Observational Incentive: May stimulate more precise observations of AGN jets
- Theoretical Development: Advances theoretical understanding of AGN physics and dark matter interactions
- Sub-GeV Dark Matter Research: Particularly suitable for exploring light dark matter where direct detection is challenging
- Majorana Dark Matter: Circumvents velocity suppression, providing effective constraint means
- Radio Astronomy: Provides new physical interpretation framework for radio observations
- Multi-Messenger Astrophysics: Studies dark matter properties combining different observational methods
This paper cites 79 relevant references spanning multiple fields including dark matter theory, AGN physics, and direct detection experiments. Key references include Gondolo & Silk (1999) on pioneering dark matter spike work and recent results from DAMIC-M and PandaX-4T experiments.