Detecting dark matter substructure with lensed quasars in optical bands

Basic Information

• Paper ID: 2510.09148
• Title: Detecting dark matter substructure with lensed quasars in optical bands
• Authors: Jianxiang Liu, Kai Liao, Yan Gong
• Classification: astro-ph.CO astro-ph.GA
• Publication Date: Preprint, October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.09148v1

Abstract

The flux ratios of multiple images in strong gravitational lensing systems provide a powerful tool for detecting dark matter substructure. Optical flux ratios of lensed quasars are typically affected by stellar microlensing effects, and therefore studies of dark matter substructure have often relied on sufficiently extended emission regions to avoid microlensing. To expand the available wavelength range for studying dark matter substructure through flux ratios and reduce dependence on specific instruments, this paper directly addresses the challenges posed by microlensing and proposes a method for detecting dark matter substructure using optical flux ratios of lensed quasars. The authors selected 100 strong lensing systems (including 90 double-image and 10 quadruple-image systems) as a representative sample and employed the Kolmogorov-Smirnov (KS) test as the statistical method. By introducing different types of dark matter substructure into these strong lensing systems, the authors demonstrate that quadruple-image systems alone provide the strongest constraints on dark matter, with tens to hundreds of independent flux ratio measurements from quadruple-image systems being sufficient to study dark matter substructure properties and constrain dark matter parameters.

Research Background and Motivation

Core Issues

The fundamental properties of dark matter remain unknown. Although the cold dark matter (CDM) model is highly successful on large scales, it faces several challenges on small scales, such as the missing satellites problem, the cusp-core problem, and the too-big-to-fail problem. Alternative models such as warm dark matter (WDM) and fuzzy dark matter (FDM) have been proposed to address these discrepancies.

Research Significance

Strong gravitational lensing, as a pure gravitational probe, is highly sensitive to the matter distribution between the observer and the source, providing a powerful and complementary method for constraining dark matter properties. Flux ratio anomalies are particularly sensitive to small-scale mass perturbations because they are related to the second derivatives of the lensing potential.

Limitations of Existing Methods

Due to the effects of stellar microlensing, previous studies have typically avoided using optical flux ratios to detect dark matter substructure, instead focusing on flux ratios measured from sufficiently large emission regions, such as narrow-line regions, radio emission, and warm dust regions. This limits the available wavelength range and increases dependence on specific instruments.

Research Motivation

With the implementation of large-scale survey projects such as LSST, Euclid, and CSST, a large number of strongly lensed quasars are expected to be discovered. LSST is expected to detect approximately 3,500 strongly lensed quasars, providing valuable opportunities for statistical analysis of dark matter using optical flux ratio data.

Core Contributions

1. Novel Methodology: Developed an optical flux ratio analysis method that simultaneously accounts for dark matter substructure and stellar microlensing effects
2. Statistical Analysis Framework: Established a statistical method for detecting dark matter substructure using the Kolmogorov-Smirnov test
3. Sample Optimization Strategy: Demonstrated that quadruple-image systems provide stronger constraints compared to double-image systems and all systems combined
4. Quantification of Observational Requirements: Determined the number of flux ratio observations needed to distinguish between different dark matter models
5. Multi-wavelength Advantages: Demonstrated that multi-wavelength flux ratios can significantly reduce the required number of quadruple-image systems

Methodology Details

Task Definition

Input: Optical flux ratio measurements of strongly lensed quasar systems Output: Detection of dark matter substructure existence and parameter constraints Constraints: Must account for the combined effects of stellar microlensing, observational errors, and model uncertainties

Model Architecture

1. Strong Lensing System Modeling

An elliptical NFW profile describes the dark matter component, an elliptical Hernquist profile represents the baryonic component, plus external shear. The lensing equation is solved using the lenstronomy software package.

2. Dark Matter Substructure Modeling

CDM Subhalos:

• Generated using SASHIMI-C
• Employ a tidally truncated NFW (TNFW) profile: $\rho(r, r_s, r_t) = \frac{\rho_s}{(r/r_s)(1 + r/r_s)^2} \cdot \frac{r_t^2}{r^2 + r_t^2}$
• Mass range considered: $10^7 \leq M_{\text{subhalo}}/M_{\odot} \leq 10^{10}$

FDM Fluctuations:

• de Broglie wavelength: $\lambda_{dB} = 150 \left(\frac{10^{-22} \text{eV}}{m_\psi}\right) \left(\frac{M_h}{10^{12} M_{\odot}}\right)^{-1/3} \text{pc}$
• Surface mass density variance: $\sigma_\Sigma^2(\boldsymbol{\xi}) = \lambda_{dB}\sqrt{\pi} \int_{-\infty}^{\infty} \rho_{DM}^2(z, \boldsymbol{\xi})dz$

3. Microlensing Simulation

Magnification maps are generated using MULES, accounting for a Salpeter initial mass function with mean stellar mass $\langle M_* \rangle = 0.2 M_{\odot}$. The accretion disk is modeled as a Gaussian distribution, with the half-light radius based on a thin disk model: $R_s = 9.7 \times 10^{15} \text{cm} \left(\frac{\lambda_{\text{rest}}}{\mu\text{m}}\right)^{4/3} \left(\frac{M_{BH}}{10^9 M_{\odot}}\right)^{2/3} \left(\frac{L}{\eta L_E}\right)^{1/3}$

Technical Innovations

1. Comprehensive Perturbation Modeling: First systematic unified modeling of dark matter substructure, stellar microlensing, and observational errors
2. Statistical Testing Method: Employs KS test to quantify the distinguishability of different dark matter models
3. Sample Selection Strategy: Discovers the superiority of quadruple-image systems, providing guidance for future observational strategies

Experimental Setup

Dataset

• Sample Source: LSST strongly lensed quasar simulation catalog
• Sample Composition: 100 strong lensing systems (90 double-image + 10 quadruple-image)
• Parameter Distribution: Median Einstein radius of 0.56 arcseconds, lens and source redshifts covering cosmologically relevant ranges

Evaluation Metrics

• Primary Metric: KS test p-value ($p_{ks}$)
• Significance Level: 95% confidence ($p_{ks} < 0.05$)
• Flux Ratio Definition: Flux ratio of minimum image to saddle-point image $\mu_{\text{minimum}}/\mu_{\text{saddle}}$

Comparison Methods

1. Null Hypothesis: Considering only microlensing and observational errors
2. CDM Model: Different subhalo mass fractions and profiles
3. FDM Model: Different boson masses ($m_\psi = 10^{-22}$ eV and $5 \times 10^{-22}$ eV)

Implementation Details

• Error Estimation: $\sigma_{SL} \approx 0.1$ mag, $\sigma_{obs} \approx 0.1$ mag
• Total Error: $\sigma_{error} = \sqrt{2(\sigma_{SL}^2 + \sigma_{obs}^2)} \approx 0.2$ mag
• Accretion Disk Size: 3 and 6 light-days half-light radius
• Statistical Repetitions: 100 repetitions per configuration

Experimental Results

Main Results

1. Detectability of Dark Matter Substructure

Using All Strong Lensing Systems:

• CDM subhalos: Require approximately 1100 (3 light-days) and 600 (6 light-days) flux ratio observations
• FDM fluctuations ($m_\psi = 10^{-22}$ eV): Require approximately 700 (3 light-days) and 300 (6 light-days) observations

Using Only Quadruple-Image Systems:

• CDM subhalos: Require only 90 (3 light-days) and 50 (6 light-days) flux ratio observations
• FDM fluctuations: Required number of observations significantly reduced

2. Dark Matter Parameter Discrimination Ability

CDM Parameter Constraints:

• Point mass vs TNFW: Require approximately 100 (3 light-days) and 50 (6 light-days) observations
• Double substructure mass fraction: Require approximately 1100 (3 light-days) and 500 (6 light-days) observations

FDM Parameter Constraints:

• $m_\psi = 10^{-22}$ eV vs $5 \times 10^{-22}$ eV: Require approximately 400 (3 light-days) and 180 (6 light-days) observations

Key Findings

1. Quadruple-Image System Advantage: Quadruple-image systems alone provide stronger constraints than systems including double-image systems
2. Accretion Disk Size Effect: Larger accretion disks are less affected by microlensing and require fewer observations
3. Perturbation Distribution Differences: CDM with increased subhalo mass fraction primarily thickens the PDF tail, while FDM with enhanced fluctuations changes the overall distribution shape

Experimental Validation

Verification through inclusion of CDM subhalos in the mass range $10^5-10^{10} M_{\odot}$ confirms that neglecting low-mass subhalos does not affect the main conclusions.

Related Work

Main Research Directions

1. Traditional Methods: Rely on extended regions such as narrow-line regions and radio emission to avoid microlensing
2. Substructure Detection: Utilize image positions, flux ratios, and time delay anomalies
3. Statistical Analysis: KS test applications in primordial black hole research provide inspiration

Innovations of This Work

Compared to existing work, this paper is the first to systematically address the microlensing problem in the optical band, expanding the wavelength range for dark matter substructure research and reducing dependence on specific instruments.

Conclusions and Discussion

Main Conclusions

1. Optical flux ratios can be used to detect dark matter substructure despite microlensing effects
2. Quadruple-image systems provide significantly stronger constraints compared to double-image systems
3. Tens to hundreds of flux ratio measurements from quadruple-image systems are sufficient to constrain dark matter parameters
4. Multi-wavelength observations can reduce the required number of quadruple-image systems by more than a factor of 6

Limitations

1. Simplified Assumptions: Neglects line-of-sight structure and FDM subhalo contributions
2. Sample Limitations: Based on simulation catalogs; actual observations may involve additional complexities
3. Error Estimation: May underestimate modeling errors for double-image systems
4. Time-Varying Effects: Does not account for intrinsic quasar variability

Future Directions

1. Complete analysis including line-of-sight structure and FDM subhalos
2. Optimize multi-wavelength strategies utilizing LSST six-band observations
3. Combine time delay measurements to reduce intrinsic variability effects
4. Extend to larger samples and more precise error modeling

In-Depth Evaluation

Strengths

1. Methodological Innovation: First systematic solution to the microlensing problem in the optical band, opening new research pathways
2. Statistical Rigor: Employs KS test to provide quantitative detectability analysis
3. Practical Value: Provides observational strategy guidance for upcoming large-scale surveys
4. Technical Comprehensiveness: Integrates multiple physical effects with complete modeling

Weaknesses

1. Model Simplification: Neglects some important physical processes to highlight the core methodology
2. Sample Representativeness: The sample of 100 systems is relatively small and may contain selection bias
3. Parameter Exploration: Considers only limited dark matter models and parameter space
4. Verification Limitations: Based on simulated data, lacks validation with actual observations

Impact

1. Disciplinary Contribution: Provides new observational tools and analysis methods for dark matter research
2. Practical Value: Directly serves the scientific objectives of major projects like LSST and Euclid
3. Method Generalization: Statistical framework extensible to other lensing phenomena studies
4. Policy Impact: May influence future survey project observational strategy formulation

Applicable Scenarios

1. Large-Scale Surveys: Particularly suitable for projects like LSST, Euclid, and CSST
2. Dark Matter Constraints: Provides new discrimination tools for CDM vs alternative models
3. Lensing Cosmology: Expands the application range of strong lensing in cosmology
4. Instrument Planning: Provides scientific requirement inputs for future space telescope design

References

The paper cites 151 relevant references, covering important works in dark matter theory, strong gravitational lensing observations, microlensing simulations, and statistical analysis across multiple fields, reflecting the comprehensive and cutting-edge nature of the research.