2025-11-22T11:49:15.930592

Dark Matter Subhalos and Higher Order Catastrophes in Gravitational Wave Lensing

Vujeva, Ezquiaga, Gilman et al.
Gravitational lensing is an invaluable probe of the nature of dark matter, and the structures it forms. Lensed gravitational waves in particular allow for unparalleled sensitivity to small scale structures within the lenses, due to the precise time resolution in combination with the continuous monitoring of the entire sky. In this work, we show two distinct ways of using strongly lensed gravitational waves to identify the presence of dark matter subhalos: \emph{i)} through higher order caustics generating high relative magnification ($μ_r > 2$), short time delay image pairs that break the caustic universality relations of single dark matter halos, which occur for $\sim 1-10$ percent of strongly lensed events in our cold dark matter models, and \emph{ii)} through the presence of more than three highly magnified images, which occur for $\sim 0.01-1$ percent of the same simulated events. We find that these results are highly sensitive to the concentrations of subhalos in our simulations, and more mildly to their number densities. The presence of low-mass subhalos increases the probability of observing wave-optics lensing in lensed gravitational waves, which is studied by solving the diffraction integral with the stationary phase approximation, as well as numerically. We also report distinct quantitative and qualitative differences in the distributions of relative magnifications and time delays for subhalo populations with increased number densities or concentrations. With the upcoming detection of strongly lensed events by ground- and space- based detectors, comparisons against these simulated distributions will provide insight into the nature of dark matter.
academic

Dark Matter Subhalos and Higher Order Catastrophes in Gravitational Wave Lensing

Basic Information

  • Paper ID: 2510.14953
  • Title: Dark Matter Subhalos and Higher Order Catastrophes in Gravitational Wave Lensing
  • Authors: Luka Vujeva, Jose María Ezquiaga, Daniel Gilman, Srashti Goyal, Miguel Zumalacárregui
  • Classification: astro-ph.CO, astro-ph.GA, gr-qc
  • Publication Date: October 17, 2025
  • Paper Link: https://arxiv.org/abs/2510.14953

Abstract

This paper investigates the effects of gravitational lensing on gravitational waves, particularly how dark matter subhalos influence gravitational wave observations through strong gravitational lensing phenomena. The study demonstrates that the presence of dark matter subhalos can be identified through two mechanisms: 1) image pairs with high relative magnification (μr > 2) and short time delays produced by higher-order caustics, which violate the caustic universality relations of single dark matter halos; 2) the production of more than three highly magnified images. The findings reveal that these results are highly sensitive to subhalo concentration, with relatively lower sensitivity to number density.

Research Background and Motivation

Problem Definition

Traditional gravitational lensing studies primarily focus on single, smooth dark matter halos (such as SIS, NFW, or SIE profiles), yet realistic galaxy-scale lenses contain numerous dark matter subhalos. The presence of these subhalos significantly alters the observational characteristics of gravitational wave lensing, but systematic studies have been lacking.

Research Significance

  1. Dark Matter Detection: Gravitational wave lensing provides a unique tool for detecting small-scale dark matter structures, capable of detecting dark matter subhalos with masses as low as 10⁷M☉
  2. Precise Temporal Resolution: Gravitational wave detectors achieve millisecond-level timing precision, far exceeding the day-level time delay precision of electromagnetic wave lensing
  3. Dark Matter Model Discrimination: Different dark matter models (cold, warm, fuzzy dark matter, etc.) exhibit significant differences in small-scale structure formation

Limitations of Existing Methods

  • Most studies assume smooth, single-lens configurations
  • Lack of systematic modeling of subhalo effects
  • Electromagnetic observations are limited by spatial resolution, making it difficult to distinguish higher-order caustic effects

Core Contributions

  1. First Systematic Study: Comprehensive investigation of dark matter subhalo effects on strong gravitational wave lensing
  2. Identification Criteria: Proposes two explicit subhalo identification standards: short time-delay image pairs with μr > 2 and more than three highly magnified images
  3. Higher-Order Caustic Theory: Applies swallow-tail and butterfly caustic theory to gravitational wave lensing analysis
  4. Wave Optical Effects: Investigates subhalo impacts on gravitational wave diffraction and interference phenomena
  5. Parameter Sensitivity Analysis: Quantifies the dependence of detection probability on subhalo concentration and number density

Methodology Details

Task Definition

Investigate how dark matter subhalos affect the observational characteristics of strong gravitational wave lensing, particularly the distribution relationships between relative magnification (μr) and time delay (ΔT).

Model Architecture

Composite Lens Model

  1. Main Halo Model: Uses elliptical singular isothermal sphere (eSIS) to describe galaxy-scale main halo
    • Lens potential: ψ(x⃗) = √(x₁² + x₂²/q²)
    • Einstein radius: ~9.4 kpc (lens redshift zL = 0.5)
    • Ellipticity: q = 0.9
  2. Subhalo Distribution: Generates subhalo properties using pyHalo code
    • Mass function: d²N/dlogmdA = Σsub(m/m₀)^(-α)
    • NFW density profile: ρ(r) = ρs/(r/rs)(1 + r/rs)²
    • Concentration-mass relation: Adopts Diemer-Joyce relation

Wave Optical Treatment

  1. Diffraction Integral: F(ω) = (ω/2πi)∫d²x⃗e^(iωT(x⃗,y⃗))
  2. Geometric Optics Approximation: Uses stationary phase approximation in high-frequency limit
  3. Time-Domain Magnification Factor: Computed via GLoW code for I(τ) and G(τ)

Technical Innovations

Higher-Order Caustic Analysis

  1. Swallow-Tail Caustic (A₄): Generated by perturbed fold caustics
    • Generating function: V = (1/5)x⁵ + (1/3)ax³ + (1/2)bx² + cx
    • Produces up to 4 real images
  2. Butterfly Caustic (A₅): Generated by perturbed cusp caustics
    • Generating function: V = (1/6)x⁶ + (1/4)ax⁴ + (1/3)bx³ + (1/2)cx² + dx
    • Produces up to 5 real images

Breakdown of Universality Relations

Traditional single-lens caustics obey strict universality relations:

  • Fold caustic: μr = 1 (constant)
  • Cusp caustic: μr ≤ 2 (upper limit)

The presence of subhalos disrupts these relations, making μr > 2 possible.

Experimental Setup

Simulation Parameters

  1. Subhalo Mass Range: 10⁷ - 10⁹M☉
  2. Concentration Models:
    • Baseline: Diemer-Joyce relation (cDJ)
    • High concentration: c = 3×cDJ, 6×cDJ
  3. Number Density Parameters: Σsub = 0.025, 0.050, 0.075

Selection Criteria

  • Magnification threshold: |μ| > 5 (only subhalos near critical curves considered)
  • Source plane sampling: High-resolution regions with minimum magnification μsr > 25

Computational Tools

  • pyHalo: Subhalo property generation
  • GLoW: Gravitational wave lensing calculations (including wave optical effects)

Experimental Results

Main Statistical Results

ModelP(μr > 2)P(Nim > 5)
Baseline Model9.0×10⁻³2.8×10⁻⁴
c = 3×cDJ8.0×10⁻²6.4×10⁻³
c = 6×cDJ8.7×10⁻²1.6×10⁻²
Σsub = 0.0501.0×10⁻²1.7×10⁻⁴
Σsub = 0.0754.0×10⁻²3.1×10⁻⁴

Key Findings

1. Breakdown of Universality Relations

  • Approximately 1% of strong lensing events in the baseline model exhibit short time-delay image pairs with μr > 2
  • This proportion rises to 8-9% in high-concentration models
  • Maximum observed relative magnification reaches μr ~ 10

2. High-Image-Number Events

  • Events with more than 5 images have probability ~10⁻⁴ in the baseline model
  • Probability increases to ~10⁻² in high-concentration models
  • These events are primarily produced by nested caustics

3. Parameter Sensitivity

  • Concentration Sensitivity: Results are highly sensitive to subhalo concentration, with significant enhancement at c = 3×cDJ
  • Number Density Impact: Increasing subhalo number density also enhances detection probability, but less effectively than increasing concentration
  • Mass Threshold: Primary contributions come from subhalos with mass > 10⁸M☉

Wave Optical Effects

Analysis of time-domain magnification factor I(τ) reveals:

  • Images near subhalos exhibit unique interference and diffraction characteristics
  • Both regular and singular parts show structures different from single-lens cases
  • These effects are more pronounced at low frequencies

Related Work

Electromagnetic Wave Lensing Studies

  • Quasar Flux Anomalies: Mao & Schneider (1998) and others studied substructure effects on quasar images
  • Image Position Deviations: Observed image positions deviate from smooth model predictions
  • Over-magnification Phenomena: Total magnification of sources near caustics exceeds expectations

Gravitational Wave Lensing Theory

  • Large-Scale Structure: Primarily focuses on galaxy and galaxy cluster-scale lensing effects
  • Compact Objects: Studies interference effects of stars and black holes
  • Wave Optics: Takahashi & Nakamura (2003) and others established foundational theoretical frameworks

Dark Matter Subhalo Simulations

  • N-body Simulations: Aquarius, Symphony, and other projects provide subhalo statistical properties
  • Observational Constraints: Strong lensing observations constrain dark matter model parameters

Conclusions and Discussion

Main Conclusions

  1. Detection Standards: Proposes two explicit subhalo detection criteria, providing theoretical foundation for future observations
  2. Parameter Constraints: Gravitational wave lensing exhibits high sensitivity to subhalo concentration, enabling constraints on dark matter models
  3. Observational Prospects: With the construction of next-generation detectors such as Einstein Telescope, Cosmic Explorer, and LISA, detection of strongly lensed gravitational waves becomes feasible

Limitations

  1. Computational Constraints: Only considers subhalos near critical curves, neglecting line-of-sight subhalos
  2. High-Frequency Calculations: Difficulty computing complete frequency response in LIGO/LISA bands due to highly oscillatory integrals
  3. Simplified Assumptions: Neglects tidal evolution effects and models all subhalos as NFW profiles

Future Directions

  1. Complete Modeling: Include line-of-sight subhalos and tidal effects
  2. High-Frequency Algorithms: Develop more efficient computational methods for high-frequency diffraction integrals
  3. Multi-Messenger Observations: Joint analysis combining electromagnetic and gravitational wave observations
  4. Alternative Dark Matter Models: Investigate effects of fuzzy dark matter, self-interacting dark matter, and other models

In-Depth Evaluation

Strengths

  1. High Innovation: First systematic study of subhalo effects on gravitational wave strong lensing, filling an important theoretical gap
  2. Rigorous Methodology: Combines precise wave optical calculations with higher-order caustic theory, with solid theoretical foundation
  3. High Practical Value: Proposed detection criteria can be directly applied to future gravitational wave observations
  4. Comprehensive Parameter Study: Systematically analyzes effects of different subhalo parameters on observational results

Weaknesses

  1. Limited Computational Scope: Only considers partial subhalo distributions due to computational complexity
  2. Incomplete Frequency Coverage: Wave optical effects at high frequencies require further refinement
  3. Simplified Models: Some physical processes (such as tidal effects) are treated with simplifications

Impact

  1. Theoretical Contribution: Opens new research directions in gravitational wave astronomy
  2. Observational Guidance: Provides important reference for scientific objectives of next-generation gravitational wave detectors
  3. Dark Matter Research: Offers new methods for detecting small-scale dark matter structures
  4. Interdisciplinary Value: Connects gravitational wave physics, cosmology, and dark matter theory

Applicable Scenarios

  1. Next-Generation Gravitational Wave Detectors: Einstein Telescope, Cosmic Explorer, LISA, etc.
  2. Strong Lensing Event Analysis: Provides theoretical framework for identifying and analyzing strongly lensed gravitational wave events
  3. Dark Matter Model Testing: Constrains parameters of different dark matter models through observational data
  4. Multi-Messenger Astronomy: Joint analysis combining electromagnetic observations

References

This paper cites 93 relevant references, covering important works in gravitational lensing theory, gravitational wave physics, dark matter simulations, and observational cosmology. The main theoretical foundations include Schneider et al.'s "Gravitational Lenses," Arnold's "Catastrophe Theory," and recent gravitational wave lensing research results.