pith. sign in

arxiv: 2606.21519 · v1 · pith:BLDBLXSXnew · submitted 2026-06-19 · 🌌 astro-ph.CO · astro-ph.HE· gr-qc

Wave-optics imprints of dark matter subhalos on strongly lensed gravitational waves. II. Saddle images and detectability

Shin'ichiro Ando This is my paper

Pith reviewed 2026-06-26 13:32 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.HEgr-qc
keywords dark matter subhalosstrongly lensed gravitational waveswave opticsLISAsaddle imagescold dark mattermatched filteramplification factor
0
0 comments X

The pith

Subhalos thread galaxy lenses and imprint percent-level frequency-dependent distortions on both minimum and saddle images of strongly lensed gravitational waves.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends wave-optics calculations of subhalo imprints from the minimum image to the saddle image. It demonstrates that a population of cold dark matter subhalos produces percent-level amplitude and phase modulations in both image parities, with the average magnification increased for the minimum and decreased for the saddle. A matched-filter search that removes lens-parameter uncertainties detects the combined signal above 5 sigma in 62 percent of Monte Carlo realizations when the source sits close to the caustic at small impact parameter. Folding the per-event significances with optimistic strong-lensing rates gives an expected 10 to 20 substructure detections across the LISA mission, opening a new window on subhalos in the 10^4 to 10^7 solar-mass range.

Core claim

Across an ensemble of cold dark matter subhalo realizations, subhalos induce percent-level amplitude and phase modulations in both image parities, while the mean (de)magnification splits by parity: the minimum is net magnified and the saddle net demagnified. Demodulating the macro-image interference recovers the per-image modulations, and a matched-filter analysis that projects out the lens parameters yields a combined detection above 5 sigma in 62 percent of realizations for fiducial massive-black-hole-binary sources of total mass about 10^6 solar masses at redshift 1.5, provided the source lies close to the lens caustic at small impact parameter y_src less than or equal to 0.1. Folding the

What carries the argument

Time-domain evaluation of the amplification factor at saddle points, which handles open equal-arrival-time contours as a small difference of large terms, followed by matched-filter projection that removes lens-parameter uncertainties.

If this is right

  • Mean magnification splits by image parity, with minima net magnified and saddles net demagnified.
  • Demodulation of macro-image interference isolates the per-image subhalo modulations.
  • The approach probes substructure masses between 10^4 and 10^7 solar masses that are inaccessible to electromagnetic observations.
  • Strongly lensed gravitational waves become a complementary probe of dark-matter substructure.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the fraction of sources near caustics is lower than assumed, the total number of detections would scale down proportionally.
  • Parity-dependent mean shifts could help separate subhalo signals from other lens-model uncertainties in future analyses.
  • The same time-domain saddle method could be applied to other wave-optics problems involving open contours.

Load-bearing premise

The source must lie close to the lens caustic with impact parameter y_src less than or equal to 0.1 and the strong-lensing rate forecasts must be realized at the optimistic level used.

What would settle it

A survey of strongly lensed gravitational-wave events with sources satisfying y_src less than or equal to 0.1 that shows far fewer than 62 percent yielding combined 5-sigma detections would falsify the projected detection rate.

Figures

Figures reproduced from arXiv: 2606.21519 by Shin'ichiro Ando.

Figure 1
Figure 1. Figure 1: FIG. 1. Iso-arrival-time (iso-Fermat) topology for three realizations, one per row, selected to span the range of subhalo [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. WO signal of the [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Ensemble distribution of the net (de)magnification dc [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The same as Fig [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Physical origin of the parity asymmetry. (a) dc versus the local tidal strength [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. How a dark matter substructure imprint enters a strongly lensed MBHB signal, for a strongly imprinted realization [the [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Distribution of the subhalo detection significance [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Detectability versus source impact parameter [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
read the original abstract

Wave-optics interference in strongly lensed gravitational waves is a new interferometric probe of dark matter substructure: a subhalo population threading a galaxy-scale lens imprints frequency-dependent distortions on the amplification factor of each macro image. In a companion paper (arXiv:2603.04267), we computed these imprints for the magnified minimum image. Here, we extend the calculation to the saddle-point image and we assess the detectability of the combined signal with the Laser Interferometer Space Antenna (LISA). Evaluating the amplification factor at a saddle is numerically delicate, because the equal-arrival-time contours are open and the subhalo signal is a small difference of large terms; we present a time-domain method that resolves it. Across a Monte Carlo ensemble of cold dark matter subhalo realizations, subhalos induce percent-level amplitude and phase modulations in both image parities, while the mean (de)magnification splits by parity: the minimum is net magnified and the saddle net demagnified. Demodulating the macro-image interference recovers the per-image modulations, and a matched-filter analysis that projects out the lens parameters yields a combined detection above $5\sigma$ in $62\%$ of realizations for fiducial massive-black-hole-binary sources of total mass $\sim10^{6}\,M_\odot$ at redshift $1.5$, provided the source lies close to the lens caustic at small impact parameter $y_{\rm src}\lesssim0.1$. Folding these naive per-event significances through optimistic strong-lensing rate forecasts yields $10$-$20$ substructure detections over the LISA mission. Strongly lensed gravitational waves are thus a sensitive, complementary probe of substructure at $10^{4}$-$10^{7}\,M_\odot$ scales inaccessible to electromagnetic observations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript extends wave-optics calculations of dark matter subhalo imprints on strongly lensed gravitational waves from minimum to saddle images. It introduces a time-domain method to evaluate the amplification factor at saddles despite open equal-arrival-time contours. Monte Carlo ensembles of CDM subhalos show percent-level amplitude and phase modulations in both parities, with mean (de)magnification splitting by parity. After demodulating macro-image interference, a matched-filter analysis projecting out lens parameters yields combined >5σ detections in 62% of realizations for fiducial ~10^6 M_⊙ MBHB sources at z=1.5 with y_src ≲0.1; folding through optimistic strong-lensing rates projects 10-20 substructure detections over the LISA mission.

Significance. If the numerical method holds and the geometric/rate assumptions are justified, the work establishes strongly lensed GWs as a complementary probe of subhalos at 10^4-10^7 M_⊙ scales. The parity-dependent effects, time-domain saddle method, and matched-filter projection are technically novel contributions that could enable new constraints inaccessible to electromagnetic observations.

major comments (2)
  1. [Abstract] Abstract (final paragraph): the 62% fraction and 10-20 detection forecast are stated only for sources satisfying y_src ≲0.1 and are obtained by folding per-event significances through optimistic strong-lensing rate forecasts; the manuscript must quantify the fraction of sources meeting the impact-parameter threshold or include sensitivity tests on the rate model, as both directly scale the headline yield.
  2. [Time-domain method] Description of the time-domain method: the method is introduced to resolve numerical delicacy arising from open contours and small differences of large terms at saddles, yet no explicit convergence tests, error propagation, or stability metrics across the Monte Carlo ensemble are referenced; this leaves the robustness of the reported percent-level modulations and the 62% detection rate unverified.
minor comments (1)
  1. The abstract refers to 'fiducial' source parameters (total mass ~10^6 M_⊙, z=1.5); clarify whether these are held fixed or varied when reporting the 62% fraction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the technical contributions. We address each major comment below and will revise the manuscript to strengthen the presentation of the results.

read point-by-point responses

  1. Referee: [Abstract] Abstract (final paragraph): the 62% fraction and 10-20 detection forecast are stated only for sources satisfying y_src ≲0.1 and are obtained by folding per-event significances through optimistic strong-lensing rate forecasts; the manuscript must quantify the fraction of sources meeting the impact-parameter threshold or include sensitivity tests on the rate model, as both directly scale the headline yield.

    Authors: The abstract already conditions the quoted numbers on y_src ≲0.1. To address the scaling concern, we will add a short paragraph in the discussion section that (i) notes the dependence of the y_src threshold fraction on the assumed lens and source populations and (ii) presents a brief sensitivity test varying the strong-lensing rate model by factors of a few. These additions will make the conditional nature of the forecast explicit without altering the core results. revision: yes

  2. Referee: [Time-domain method] Description of the time-domain method: the method is introduced to resolve numerical delicacy arising from open contours and small differences of large terms at saddles, yet no explicit convergence tests, error propagation, or stability metrics across the Monte Carlo ensemble are referenced; this leaves the robustness of the reported percent-level modulations and the 62% detection rate unverified.

    Authors: Convergence tests (varying contour cutoff radius, time sampling, and integration tolerances) were performed internally and confirmed that the reported percent-level modulations are stable to ≲0.1% across the ensemble. These checks were not documented in the submitted manuscript. We will add a concise methods subsection (or short appendix) summarizing the convergence criteria, error estimates, and stability metrics to allow readers to assess robustness directly. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results from independent Monte Carlo and matched-filter analysis

full rationale

The paper derives its claims via forward Monte Carlo realizations of CDM subhalos, a new time-domain method for evaluating the saddle amplification factor, demodulation of macro-image interference, and a matched-filter projection that removes lens parameters to obtain per-realization detection significances. These steps are self-contained and do not reduce any output quantity to a fitted parameter or prior result by construction. The companion-paper citation supplies only the minimum-image baseline and is not invoked to justify the saddle extension or the 62% detection fraction; both the geometric condition y_src ≲ 0.1 and the lensing-rate forecast are external inputs, not internal circular steps. No self-definitional, fitted-input, or uniqueness-theorem patterns appear in the supplied text.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete; typical CDM subhalo mass-function parameters and the strong-lensing optical depth are implicitly taken from prior literature without independent verification here.

free parameters (1)
  • y_src impact parameter threshold
    The 62% detection fraction is conditioned on y_src ≲ 0.1; this geometric cut is chosen rather than derived.
axioms (1)
  • domain assumption Cold dark matter subhalo population statistics are correctly described by the Monte Carlo realizations used
    The ensemble is drawn from a standard CDM model; the abstract does not test alternative dark matter models.

pith-pipeline@v0.9.1-grok · 5876 in / 1465 out tokens · 16600 ms · 2026-06-26T13:32:18.395403+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

90 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    The iso-arrival contour does not close. As Eq. (11) shows, the open-arc template is not a finite closed- loop value but diverges, growing logarithmically with the patch radiusR c at which the contour is cut

  2. [2]

    Toward the saddle delay (s→0), the templateI quad diverges as−2 √µln|s|, andδIis a percent-level residual sitting on top of this large, singular background

    The subhalo signal rides on a divergent back- ground. Toward the saddle delay (s→0), the templateI quad diverges as−2 √µln|s|, andδIis a percent-level residual sitting on top of this large, singular background. Recovering it requires sub- tracting the macro background accurately where it is largest, rather than reading the signal off directly as one can a...

  3. [3]

    An NFW subhalo falls off too slowly for its imprint to be compact, leaving a non-cancelling tail inδI that cannot be traced to arbitrarily large delay

    The perturbing potential reaches to infinity. An NFW subhalo falls off too slowly for its imprint to be compact, leaving a non-cancelling tail inδI that cannot be traced to arbitrarily large delay. We therefore adopt the truncated-NFW profile of Eq. (5): its finite total mass gives a single clean far- field logarithmψ→mlnrthat can be subtracted analytical...

  4. [4]

    (11) is the co-area inte- gral (9) of the bare local quadratic

    The patch-bounded quadratic-saddle template The templateI quad of Eq. (11) is the co-area inte- gral (9) of the bare local quadratic. In the Hessian eigen- frame,ϕ−ϕ sad =λ +x2 +/2− |λ −|x2 −/2, the co-area in- tegrand reduces to a constant in the contour’s natural parameterζ, dℓ |∇ϕ| = √µ dζ,(A1) thesamefor both image types; only the range ofζand its tri...

  5. [5]

    Tracing the open saddle contour The open saddle contour is traced as four half-arcs (two hyperbola branches, two vertex halves). Each half-arc is integrated as an ordinary differential equation along the contour, in the co-area parameterσwithdI/dσ=R, by an adaptive, error-controlled Runge–Kutta integrator (GSL’s RK8(7), absolute and relative tolerance 10 ...

  6. [6]

    To first order the co-area integral [Eq

    Far-field tail: co-area derivation and Fourier transform Write the full arrival-time surface as the quadratic sad- dle plus the subhalo potential,ϕ=ϕ quad +ψ. To first order the co-area integral [Eq. (9)] is perturbed by δI(s) =− d ds I ϕquad=s ψ dℓ |∇ϕquad| ≡ − d ds ⟨ψ⟩s.(A4) For the open saddle the weightdℓ/|∇ϕ quad|is the s-independent constantg ∞ = 1/...

  7. [7]

    Numerical fringes from the Filon transform The residualδIis transformed in two pieces: a Filon quadrature over the window|s| ≤s hi and the analytic tail eTbeyond it (App. A 3). Handled naively, each piece imprints a fringe on the envelope that is common to all re- alizations (hence not physical subhalo signal), and we re- move both at negligible cost. The...

    2000

  8. [8]

    Matched-filter identities and the per-image noise floor We collect the matched-filter identities underlying Sec. V. With the one-sided power spectral densityS n defined by⟨˜n(f)˜n∗(f ′)⟩=S n(f)δ(f−f ′)/2 and the re- ality condition ˜n(−f) = ˜n ∗(f), the likelihood-natural two-sided inner product reduces, for real signals, to the one-sided form of Eq. (20)...

  9. [9]

    Validating the image separation: the demodulation test The claim of Sec. V B that the two macro images sep- arate losslessly is verified by a direct round-trip onF tot that is deterministic and source-model-independent: no waveform, sensitivity curve, or noise enters. For each of the 1000 realizations we take the per-image factors√µj e−iπnj /2[1+η j(f)], ...

  10. [10]

    Subhalo sampling Heavy (m >10 9 M⊙) and light (102–109 M⊙) subhalos are drawn from the semi-analyticSASHIMImodel [52– 54], with the three-dimensional number densitynsub(r)∝ (r2 +R 2 s )−3/2. Light subhalos near an image are sam- pled in a projected disc of radiusR near with the expected count fixed by the local cylinder weight wcyl =w sub πR2 sample Iz Vh...

  11. [11]

    Truncated-NFW lens and mass match Each subhalo is then= 1 truncated-NFW profile of Eq. (5), whose enclosed mass within radiusc r s, in units of 4πρsr3 s , is I1(c) = c2 (c2 + 1)2 (c2 −1) lnc+πc−(c 2 + 1) .(C2) SASHIMIfixes the truncation assuming an abrupt three-dimensional cut,I ∞(c) = ln(1 +c)−c/(1 +c). To preserve the bound mass exactly under the profi...

  12. [12]

    The macro fieldψ macro collects the host, galaxy, and heavy subhalos together with the light subhalos that the thresholds of Sec

    External-field decomposition The potential is split asψ=ψ macro +δψ. The macro fieldψ macro collects the host, galaxy, and heavy subhalos together with the light subhalos that the thresholds of Sec. II B classify as GO, whileδψretains only the WO light subhalos. Because the GO light subhalos perturb the arrival-time surface, we re-solve for the stationary...

  13. [13]

    Ando, Wave-Optics Imprints of Dark Matter Subha- los on Strongly Lensed Gravitational Waves, (2026), arXiv:2603.04267 [astro-ph.CO]

    S. Ando, Wave-Optics Imprints of Dark Matter Subha- los on Strongly Lensed Gravitational Waves, (2026), arXiv:2603.04267 [astro-ph.CO]

    arXiv 2026

  14. [14]

    J. S. Bullock and M. Boylan-Kolchin, Small-Scale Chal- lenges to the ΛCDM Paradigm, Ann. Rev. Astron. As- trophys.55, 343 (2017), arXiv:1707.04256 [astro-ph.CO]

    arXiv 2017

  15. [15]

    Zavala and C

    J. Zavala and C. S. Frenk, Dark matter haloes and sub- haloes, Galaxies7, 81 (2019), arXiv:1907.11775 [astro- ph.CO]

    arXiv 2019

  16. [16]

    Springel, J

    V. Springel, J. Wang, M. Vogelsberger, A. Ludlow, A. Jenkins, A. Helmi, J. F. Navarro, C. S. Frenk, and S. D. M. White, The Aquarius Project: the subhalos of galactic halos, Mon. Not. Roy. Astron. Soc.391, 1685 (2008), arXiv:0809.0898 [astro-ph]

    Pith/arXiv arXiv 2008

  17. [17]

    M. R. Lovell, C. S. Frenk, V. R. Eke, A. Jenkins, L. Gao, and T. Theuns, The properties of warm dark matter haloes, Mon. Not. Roy. Astron. Soc.439, 300 (2014), arXiv:1308.1399 [astro-ph.CO]

    Pith/arXiv arXiv 2014

  18. [18]

    Dekker, S

    A. Dekker, S. Ando, C. A. Correa, and K. C. Y. Ng, Warm dark matter constraints using Milky Way satellite observations and subhalo evolution modeling, Phys. Rev. D106, 123026 (2022), arXiv:2111.13137 [astro-ph.CO]

    arXiv 2022

  19. [19]

    Tulin and H.-B

    S. Tulin and H.-B. Yu, Dark Matter Self-interactions and Small Scale Structure, Phys. Rept.730, 1 (2018), arXiv:1705.02358 [hep-ph]

    Pith/arXiv arXiv 2018

  20. [20]

    S. Ando, S. Horigome, E. O. Nadler, D. Yang, and H.-B. Yu, SASHIMI-SIDM: semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales, JCAP02, 053, arXiv:2403.16633 [astro-ph.CO]

    arXiv

  21. [21]

    S. Ando, K. Hayashi, S. Horigome, M. Ibe, and S. Shirai, Stringent Constraints on Self-Interacting Dark Matter Using Milky-Way Satellite Galaxies Kinematics, (2025), arXiv:2503.13650 [astro-ph.CO]

    arXiv 2025

  22. [22]

    L. Hui, J. P. Ostriker, S. Tremaine, and E. Witten, Ul- tralight scalars as cosmological dark matter, Phys. Rev. D95, 043541 (2017), arXiv:1610.08297 [astro-ph.CO]

    Pith/arXiv arXiv 2017

  23. [23]

    Banik, J

    N. Banik, J. Bovy, G. Bertone, D. Erkal, and T. J. L. de Boer, Evidence of a population of dark subhaloes fromGaiaand Pan-STARRS observations of the GD- 1 stream, Mon. Not. Roy. Astron. Soc.502, 2364 (2021), arXiv:1911.02662 [astro-ph.GA]

    arXiv 2021

  24. [24]

    Vegettiet al., Strong Gravitational Lensing as a Probe of Dark Matter, Space Sci

    S. Vegettiet al., Strong Gravitational Lensing as a Probe of Dark Matter, Space Sci. Rev.220, 58 (2024), arXiv:2306.11781 [astro-ph.CO]

    arXiv 2024

  25. [25]

    J. D. Simon, The Faintest Dwarf Galaxies, Ann. Rev. As- tron. Astrophys.57, 375 (2019), arXiv:1901.05465 [astro- ph.GA]

    arXiv 2019

  26. [26]

    E. O. Nadler, S. Birrer, D. Gilman, R. H. Wechsler, X. Du, A. Benson, A. M. Nierenberg, and T. Treu, Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies, Astrophys. J.917, 7 (2021), arXiv:2101.07810 [astro- ph.CO]

    arXiv 2021

  27. [27]

    Takahashi and T

    R. Takahashi and T. Nakamura, Wave effects in grav- itational lensing of gravitational waves from chirping binaries, Astrophys. J.595, 1039 (2003), arXiv:astro- ph/0305055

    arXiv 2003

  28. [28]

    Jung and C

    S. Jung and C. S. Shin, Gravitational-Wave Fringes at LIGO: Detecting Compact Dark Matter by Gravi- tational Lensing, Phys. Rev. Lett.122, 041103 (2019), arXiv:1712.01396 [astro-ph.CO]

    Pith/arXiv arXiv 2019

  29. [29]

    K.-H. Lai, O. A. Hannuksela, A. Herrera-Mart´ ın, J. M. Diego, T. Broadhurst, and T. G. F. Li, Discovering intermediate-mass black hole lenses through gravita- tional wave lensing, Phys. Rev. D98, 083005 (2018), arXiv:1801.07840 [gr-qc]

    Pith/arXiv arXiv 2018

  30. [30]

    Tambalo, M

    G. Tambalo, M. Zumalac´ arregui, L. Dai, and M. H.-Y. Cheung, Gravitational wave lensing as a probe of halo properties and dark matter, Phys. Rev. D108, 103529 (2023), arXiv:2212.11960 [astro-ph.CO]. 19

    arXiv 2023

  31. [31]

    Lin, J.-d

    X.-y. Lin, J.-d. Zhang, L. Dai, S.-J. Huang, and J. Mei, Detecting strong gravitational lensing of gravitational waves with TianQin, Phys. Rev. D108, 064020 (2023), arXiv:2304.04800 [gr-qc]

    arXiv 2023

  32. [32]

    Savastano, G

    S. Savastano, G. Tambalo, H. Villarrubia-Rojo, and M. Zumalacarregui, Weakly lensed gravitational waves: Probing cosmic structures with wave-optics features, Phys. Rev. D108, 103532 (2023), arXiv:2306.05282 [gr- qc]

    arXiv 2023

  33. [33]

    Zumalac´ arregui, Lens Stochastic Diffraction: A Sig- nature of Compact Structures in Gravitational-Wave Data, (2024), arXiv:2404.17405 [gr-qc]

    M. Zumalac´ arregui, Lens Stochastic Diffraction: A Sig- nature of Compact Structures in Gravitational-Wave Data, (2024), arXiv:2404.17405 [gr-qc]

    arXiv 2024

  34. [34]

    Zumalac´ arregui and X

    M. Zumalac´ arregui and X. Shan, Effective description of lensed gravitational waves diffracted by stellar fields, (2026), arXiv:2606.17765 [astro-ph.HE]

    Pith/arXiv arXiv 2026

  35. [35]

    Urrutia and V

    J. Urrutia and V. Vaskonen, Dark timbre of grav- itational waves, Phys. Rev. D111, 123047 (2025), arXiv:2402.16849 [gr-qc]

    arXiv 2025

  36. [36]

    Villarrubia-Rojo, S

    H. Villarrubia-Rojo, S. Savastano, M. Zumalac´ arregui, L. Choi, S. Goyal, L. Dai, and G. Tambalo, Gravitational lensing of waves: Novel methods for wave-optics phenom- ena, Phys. Rev. D111, 103539 (2025), arXiv:2409.04606 [gr-qc]

    arXiv 2025

  37. [37]

    Chen and Y

    Z. Chen and Y. Lu, Gravitational Lensing of Gravita- tional Waves from Astrophysical Sources: Theory, De- tection, and Applications, Res. Astron. Astrophys.26, 062001 (2026), arXiv:2605.06321 [astro-ph.HE]

    Pith/arXiv arXiv 2026

  38. [38]

    Fairbairn, J

    M. Fairbairn, J. Urrutia, and V. Vaskonen, Microlensing of gravitational waves by dark matter structures, JCAP 07, 007, arXiv:2210.13436 [astro-ph.CO]

    arXiv

  39. [39]

    Guo and Y

    X. Guo and Y. Lu, Probing the nature of dark matter via gravitational waves lensed by small dark matter ha- los, Phys. Rev. D106, 023018 (2022), arXiv:2207.00325 [astro-ph.CO]

    arXiv 2022

  40. [40]

    Brando, S

    G. Brando, S. Goyal, S. Savastano, H. Villarrubia-Rojo, and M. Zumalac´ arregui, Signatures of dark and baryonic structures on weakly lensed gravitational waves, Phys. Rev. D111, 024068 (2025), arXiv:2407.04052 [gr-qc]

    arXiv 2025

  41. [41]

    C ¸ alı¸ skan, L

    M. C ¸ alı¸ skan, L. Ji, R. Cotesta, E. Berti, M. Kamionkowski, and S. Marsat, Observability of lensing of gravitational waves from massive black hole binaries with LISA, Phys. Rev. D107, 043029 (2023), arXiv:2206.02803 [astro-ph.CO]

    arXiv 2023

  42. [42]

    C ¸ alı¸ skan, N

    M. C ¸ alı¸ skan, N. Anil Kumar, L. Ji, J. M. Ezquiaga, R. Cotesta, E. Berti, and M. Kamionkowski, Probing wave-optics effects and low-mass dark matter halos with lensing of gravitational waves from massive black holes, Phys. Rev. D108, 123543 (2023), arXiv:2307.06990 [astro-ph.CO]

    arXiv 2023

  43. [43]

    Oguri and R

    M. Oguri and R. Takahashi, Amplitude and phase fluc- tuations of gravitational waves magnified by strong grav- itational lensing, Phys. Rev. D106, 043532 (2022), arXiv:2204.00814 [astro-ph.CO]

    arXiv 2022

  44. [44]

    Sereno, A

    M. Sereno, A. Sesana, A. Bleuler, P. Jetzer, M. Volon- teri, and M. C. Begelman, Strong lensing of gravitational waves as seen by LISA, Phys. Rev. Lett.105, 251101 (2010), arXiv:1011.5238 [astro-ph.CO]

    Pith/arXiv arXiv 2010

  45. [45]

    Oguri, Effect of gravitational lensing on the distribu- tion of gravitational waves from distant binary black hole mergers, Mon

    M. Oguri, Effect of gravitational lensing on the distribu- tion of gravitational waves from distant binary black hole mergers, Mon. Not. Roy. Astron. Soc.480, 3842 (2018), arXiv:1807.02584 [astro-ph.CO]

    Pith/arXiv arXiv 2018

  46. [46]

    Guti´ errez and M

    J. Guti´ errez and M. Lagos, Strong-lensing rates of mas- sive black hole binaries in LISA, Phys. Rev. D112, 123512 (2025), arXiv:2510.02061 [astro-ph.CO]

    arXiv 2025

  47. [47]

    Hu and Y.-L

    W.-R. Hu and Y.-L. Wu, The Taiji Program in Space for gravitational wave physics and the nature of gravity, Natl. Sci. Rev.4, 685 (2017)

    2017

  48. [48]

    Z. Luo, Y. Wang, Y. Wu, W. Hu, and G. Jin, The Taiji program: A concise overview, PTEP2021, 05A108 (2021)

    2021

  49. [49]

    Luoet al.(TianQin), TianQin: a space-borne gravi- tational wave detector, Class

    J. Luoet al.(TianQin), TianQin: a space-borne gravi- tational wave detector, Class. Quant. Grav.33, 035010 (2016), arXiv:1512.02076 [astro-ph.IM]

    Pith/arXiv arXiv 2016

  50. [50]

    J. M. Ezquiaga, R. K. L. Lo, and L. Vujeva, Diffraction around caustics in gravitational wave lensing, Phys. Rev. D112, 043544 (2025), arXiv:2503.22648 [gr-qc]

    arXiv 2025

  51. [51]

    Vujeva, J

    L. Vujeva, J. M. Ezquiaga, D. Gilman, S. Goyal, and M. Zumalac´ arregui, Dark Matter Subhalos and Higher Order Catastrophes in Gravitational Wave Lensing, (2025), arXiv:2510.14953 [astro-ph.CO]

    arXiv 2025

  52. [52]

    Mao and P

    S.-d. Mao and P. Schneider, Evidence for substructure in lens galaxies?, Mon. Not. Roy. Astron. Soc.295, 587 (1998), arXiv:astro-ph/9707187

    Pith/arXiv arXiv 1998

  53. [53]

    P. L. Schechter and J. Wambsganss, Quasar microlens- ing at high magnification and the role of dark matter: Enhanced fluctuations and suppressed saddlepoints, As- trophys. J.580, 685 (2002), arXiv:astro-ph/0204425

    Pith/arXiv arXiv 2002

  54. [54]

    Dalal and C

    N. Dalal and C. S. Kochanek, Direct detection of CDM substructure, Astrophys. J.572, 25 (2002), arXiv:astro- ph/0111456

    arXiv 2002

  55. [55]

    C. S. Kochanek and N. Dalal, Tests for substructure in gravitational lenses, Astrophys. J.610, 69 (2004), arXiv:astro-ph/0302036

    Pith/arXiv arXiv 2004

  56. [56]

    J. M. Diego, O. A. Hannuksela, P. L. Kelly, T. Broad- hurst, K. Kim, T. G. F. Li, G. F. Smoot, and G. Pagano, Observational signatures of microlensing in gravitational waves at LIGO/Virgo frequencies, Astron. Astrophys. 627, A130 (2019), arXiv:1903.04513 [astro-ph.CO]

    Pith/arXiv arXiv 2019

  57. [57]

    Mishra, A

    A. Mishra, A. K. Meena, A. More, S. Bose, and J. S. Bagla, Gravitational lensing of gravitational waves: effect of microlens population in lensing galaxies, Mon. Not. Roy. Astron. Soc.508, 4869 (2021), arXiv:2102.03946 [astro-ph.CO]

    arXiv 2021

  58. [58]

    S. M. C. Yeung, M. H. Y. Cheung, E. Seo, J. A. J. Gais, O. A. Hannuksela, and T. G. F. Li, Detectability of mi- crolensed gravitational waves, Mon. Not. Roy. Astron. Soc.526, 2230 (2023), arXiv:2112.07635 [gr-qc]

    arXiv 2023

  59. [59]

    X. Shan, G. Li, X. Chen, W. Zheng, and W. Zhao, Wave effect of gravitational waves intersected with a microlens field: A new algorithm and supplementary study, Sci. China Phys. Mech. Astron.66, 239511 (2023), arXiv:2208.13566 [astro-ph.CO]

    arXiv 2023

  60. [60]

    Aghanimet al.(Planck), Planck 2018 results

    N. Aghanimet al.(Planck), Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]

    Pith/arXiv arXiv 2018

  61. [61]

    J. F. Navarro, C. S. Frenk, and S. D. M. White, The Structure of cold dark matter halos, Astrophys. J.462, 563 (1996), arXiv:astro-ph/9508025

    Pith/arXiv arXiv 1996

  62. [62]

    J. F. Navarro, C. S. Frenk, and S. D. M. White, A Uni- versal density profile from hierarchical clustering, Astro- phys. J.490, 493 (1997), arXiv:astro-ph/9611107

    Pith/arXiv arXiv 1997

  63. [63]

    C. A. Correa, J. S. B. Wyithe, J. Schaye, and A. R. Duffy, The accretion history of dark matter haloes – III. A phys- ical model for the concentration–mass relation, Mon. Not. Roy. Astron. Soc.452, 1217 (2015), arXiv:1502.00391 [astro-ph.CO]. 20

    Pith/arXiv arXiv 2015

  64. [64]

    Hiroshima, S

    N. Hiroshima, S. Ando, and T. Ishiyama, Modeling evolu- tion of dark matter substructure and annihilation boost, Phys. Rev. D97, 123002 (2018), arXiv:1803.07691 [astro- ph.CO]

    Pith/arXiv arXiv 2018

  65. [65]

    S. Ando, T. Ishiyama, and N. Hiroshima, Halo Substruc- ture Boosts to the Signatures of Dark Matter Annihi- lation, Galaxies7, 68 (2019), arXiv:1903.11427 [astro- ph.CO]

    arXiv 2019

  66. [66]

    S. Ando, A. Geringer-Sameth, N. Hiroshima, S. Hoof, R. Trotta, and M. G. Walker, Structure formation mod- els weaken limits on WIMP dark matter from dwarf spheroidal galaxies, Phys. Rev. D102, 061302 (2020), arXiv:2002.11956 [astro-ph.CO]

    arXiv 2020

  67. [67]

    E. A. Baltz, P. Marshall, and M. Oguri, Analytic models of plausible gravitational lens potentials, JCAP01, 015, arXiv:0705.0682 [astro-ph]

    Pith/arXiv arXiv

  68. [68]

    Tambalo, M

    G. Tambalo, M. Zumalac´ arregui, L. Dai, and M. H.-Y. Cheung, Lensing of gravitational waves: Efficient wave- optics methods and validation with symmetric lenses, Phys. Rev. D108, 043527 (2023), arXiv:2210.05658 [gr- qc]

    arXiv 2023

  69. [69]

    T. T. Nakamura and S. Deguchi, Wave Optics in Grav- itational Lensing, Prog. Theor. Phys. Suppl.133, 137 (1999)

    1999

  70. [70]

    Carrillo Gonzalez, V

    M. Carrillo Gonzalez, V. De Luca, A. Garoffolo, J. Parra- Martinez, and M. Trodden, Scattering perspective on gravitational lensing, Phys. Rev. D113, 024024 (2026), arXiv:2511.15797 [hep-th]

    arXiv 2026

  71. [71]

    Ephremidze, M

    N. Ephremidze, M. Kamionkowski, and C. Dvorkin, Fast Fourier Transform evaluation of the Fresnel integral for gravitational-wave lensing, (2026), arXiv:2603.12333 [astro-ph.CO]

    arXiv 2026

  72. [72]

    Ulmer and J

    A. Ulmer and J. Goodman, Femtolensing: Beyond the semiclassical approximation, Astrophys. J.442, 67 (1995), arXiv:astro-ph/9406042

    Pith/arXiv arXiv 1995

  73. [73]

    Dai and T

    L. Dai and T. Venumadhav, On the waveforms of gravitationally lensed gravitational waves, (2017), arXiv:1702.04724 [gr-qc]

    Pith/arXiv arXiv 2017

  74. [74]

    Iserles and S

    A. Iserles and S. P. Nørsett, Efficient quadrature of highly oscillatory integrals using derivatives, Proc. Roy. Soc. A 461, 1383 (2005)

    2005

  75. [75]

    L. S. Finn, Detection, measurement and gravitational radiation, Phys. Rev. D46, 5236 (1992), arXiv:gr- qc/9209010

    arXiv 1992

  76. [76]

    Cutler and E

    C. Cutler and E. E. Flanagan, Gravitational waves from merging compact binaries: How accurately can one extract the binary’s parameters from the inspiral wave form?, Phys. Rev. D49, 2658 (1994), arXiv:gr- qc/9402014

    arXiv 1994

  77. [77]

    J. M. Ezquiaga, D. E. Holz, W. Hu, M. Lagos, and R. M. Wald, Phase effects from strong gravitational lensing of gravitational waves, Phys. Rev. D103, 064047 (2021), arXiv:2008.12814 [gr-qc]

    arXiv 2021

  78. [78]

    S. Husa, S. Khan, M. Hannam, M. P¨ urrer, F. Ohme, X. Jim´ enez Forteza, and A. Boh´ e, Frequency-domain gravitational waves from nonprecessing black-hole bina- ries. I. New numerical waveforms and anatomy of the sig- nal, Phys. Rev. D93, 044006 (2016), arXiv:1508.07250 [gr-qc]

    Pith/arXiv arXiv 2016

  79. [79]

    S. Khan, S. Husa, M. Hannam, F. Ohme, M. P¨ urrer, X. Jim´ enez Forteza, and A. Boh´ e, Frequency-domain gravitational waves from nonprecessing black-hole bi- naries. II. A phenomenological model for the ad- vanced detector era, Phys. Rev. D93, 044007 (2016), arXiv:1508.07253 [gr-qc]

    Pith/arXiv arXiv 2016

  80. [80]

    Robson, N

    T. Robson, N. J. Cornish, and C. Liu, The construction and use of LISA sensitivity curves, Class. Quant. Grav. 36, 105011 (2019), arXiv:1803.01944 [astro-ph.HE]

    Pith/arXiv arXiv 2019

Showing first 80 references.