pith. sign in

arxiv: 2510.11790 · v3 · submitted 2025-10-13 · 🌀 gr-qc · astro-ph.CO· astro-ph.GA· astro-ph.HE

False Alarm Rates in Detecting Gravitational Wave Lensing from Astrophysical Coincidences: Insights with Model-Independent Technique GLANCE

Pith reviewed 2026-05-18 07:15 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.COastro-ph.GAastro-ph.HE
keywords gravitational wave lensingfalse alarm ratebinary black hole mergersmodel-independent detectionastrophysical coincidencesLIGO sensitivitytime delay analysis
0
0 comments X

The pith

Simulations show only 0.01 percent of binary black hole pairs falsely appear as lensed gravitational waves at long delays.

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

The paper tests how often pairs of unlensed gravitational wave events from merging binary black holes get misidentified as lensed pairs because their sky locations overlap by chance. It applies the GLANCE pipeline to a large simulated population of such events and measures the resulting false alarm rate at different signal-to-noise thresholds and time delays. A reader would care because true lensing events are rare and valuable for studying intervening matter, so any technique must keep false positives low enough that detections can be trusted. The work maps how this false positive rate changes with source mass, time delay, and magnification, showing the conditions under which GLANCE can separate real lensing from ordinary astrophysical coincidences using current LIGO data.

Core claim

Applying the model-independent GLANCE technique to a simulated population of merging binary black holes, the authors find that approximately 0.01 percent of event pairs can be falsely classified as lensed when using a signal-to-noise ratio threshold of 1.5, particularly for time delays of around 1000 days or more. They map the false alarm rate distribution across the parameter space of source masses, delay times, and lensing magnification factors, demonstrating the conditions under which GLANCE can confidently identify true lensed pairs with existing LIGO sensitivity.

What carries the argument

The GLANCE pipeline, which performs model-independent consistency checks on pairs of events to test for gravitational lensing signatures while quantifying contamination from unlensed astrophysical coincidences.

If this is right

  • GLANCE can be used to set reliable detection thresholds for lensed gravitational wave pairs with current LIGO sensitivity.
  • The false alarm rate remains low enough across wide ranges of source masses and magnifications to support searches for time delays of 1000 days or longer.
  • The same simulation framework directly informs the expected contamination levels for next-generation detectors that will record many more events.
  • Parameter-dependent false alarm maps allow observers to adjust search criteria based on the properties of candidate pairs.

Where Pith is reading between the lines

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

  • If the low false alarm rate holds when GLANCE is applied to real events, it would reduce the statistical burden of confirming individual lensing candidates through follow-up observations.
  • Extending the same simulation approach to other source classes such as neutron-star mergers could reveal whether false alarm rates change with different mass and distance distributions.
  • Longer observing runs will naturally produce more pairs at large time delays, so the reported scaling of false alarm rate with delay time sets a clear target for how much additional data can be tolerated before contamination becomes noticeable.

Load-bearing premise

The simulated population of merging binary black holes accurately represents the true astrophysical distribution of sources, noise properties, and sky-localization errors that GLANCE will encounter in real LIGO data.

What would settle it

Running GLANCE on actual LIGO data and finding a false positive fraction that is much higher or lower than the simulated 0.01 percent rate at comparable time delays and thresholds would indicate that the simulation does not capture the relevant astrophysical or detector effects.

Figures

Figures reproduced from arXiv: 2510.11790 by Aniruddha Chakraborty, Suvodip Mukherjee.

Figure 1
Figure 1. Figure 1: FIG. 1. In this figure, we present the outline of the workflow for GW lensing false alarm rate (FAR) estimation with [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. In this figure, we show a schematic diagram of GW lensing by a massive object in the GO-lensing regime. All distances [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. In this figure, we show that given a pair of ‘similar’ GW sources, there can be two way to explain. First guess is that [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. In this figure, we have categorized the possible scenarios that can account for the false positives when detecting lensing [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. In this figure, we have demonstrated the sky-localization error for two events (blue and red regions) and the region [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. In this figure we show the chirp-mass distribution and delay time distribution of the sky-overlapping pairs. The chirp [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. In this figure, we present the lensing false positive probability (FPP in percentage) as a function of lensing SNR [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. In this figure, in the time-delay vs chirp-mass plane, we show the regions which are more susceptible to false lensing [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. In this figure, we show the time-delay vs lensing threshold SNR (which is proportional to the square root of the [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. The figure shows the primary mass distribution chosen for BBH to be a power law (with negative exponent) + [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
read the original abstract

Lensing of gravitational waves (GWs) due to intervening massive astrophysical systems between the source and the observer is an inevitable consequence of the general theory of relativity, which can produce multiple GW events with overlapping sky localization error. However, the confirmed detection of such a unique astrophysical phenomenon is challenging due to several sources of contamination, ranging from detector noise to astrophysical uncertainties. Robust model-independent search techniques that can mitigate noise contamination have been developed in the past. In this study, we explore the astrophysical uncertainty associated with incorrectly classifying a pair of unlensed GW events as a lensed pair and the associated false alarm rate (FAR) depending on the GW source properties. To understand the effect of unlensed astrophysical GW sources in producing false lensing detections, we perform a model-independent test using the pipeline GLANCE on a simulated population of merging binary black holes (BBHs). We find that $\sim$ 0.01% of the pair of events can be falsely classified as lensed with a lensing threshold signal-to-noise ratio of 1.5, appearing at a time delay between the pair of events of $\sim$ 1000 days or more. We show the FAR distribution for the parameter space of the GW source masses, delay time, and lensing magnification parameter over which the model-independent technique GLANCE can confidently detect lensed GW pair with the current LIGO detector sensitivity. In the future, this technique will be useful in understanding the lensing FAR for next-generation GW detectors, which can observe more GW sources.

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 / 2 minor

Summary. The manuscript applies the model-independent GLANCE pipeline to a simulated population of unlensed binary black hole mergers to quantify the false alarm rate arising from astrophysical coincidences. It reports that ∼0.01% of event pairs are falsely classified as lensed at a lensing SNR threshold of 1.5, particularly for time delays ≳1000 days, and presents the FAR distribution over source masses, delay times, and magnification parameters at current LIGO sensitivity.

Significance. If the simulations are representative, the work supplies a concrete estimate of contamination in searches for gravitationally lensed GW events, which is relevant for interpreting candidate pairs in existing LIGO data and for planning analyses with next-generation detectors. The model-independent character of GLANCE is a methodological strength that could help isolate astrophysical false positives from instrumental effects.

major comments (2)
  1. [Simulation description] Simulation description: the paper states that a simulated population of merging BBHs was fed into GLANCE but supplies no details on the component-mass distribution, redshift range, merger-rate density, or the modeling of sky-localization uncertainties and detector noise. Because the reported ∼0.01% false-alarm fraction is obtained by counting triggers in the tails of these distributions, the absence of these specifications directly affects the reliability of the central quantitative claim.
  2. [Results] Results: the ∼0.01% figure is presented without error bars, bootstrap uncertainties, or the total number of simulated pairs, making it impossible to judge the statistical precision of the false-alarm rate or its sensitivity to simulation choices.
minor comments (2)
  1. [Abstract] The abstract refers to showing the FAR distribution over parameter space but does not indicate whether this is conveyed in a figure, table, or text summary.
  2. [Methods] Clarify the precise definition of the 'lensing threshold signal-to-noise ratio of 1.5' and its relation to the GLANCE detection statistic.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments, which have helped us improve the clarity and completeness of our manuscript. We address each major comment below and outline the revisions we will make.

read point-by-point responses
  1. Referee: Simulation description: the paper states that a simulated population of merging BBHs was fed into GLANCE but supplies no details on the component-mass distribution, redshift range, merger-rate density, or the modeling of sky-localization uncertainties and detector noise. Because the reported ∼0.01% false-alarm fraction is obtained by counting triggers in the tails of these distributions, the absence of these specifications directly affects the reliability of the central quantitative claim.

    Authors: We agree with the referee that detailed specifications of the simulation are necessary to evaluate the robustness of the reported false alarm rate. In the revised version of the manuscript, we will add a comprehensive description of the simulated population, including the component-mass distribution, redshift range, merger-rate density, and the modeling of sky-localization uncertainties and detector noise. This will provide the necessary context for interpreting the ∼0.01% false-alarm fraction. revision: yes

  2. Referee: Results: the ∼0.01% figure is presented without error bars, bootstrap uncertainties, or the total number of simulated pairs, making it impossible to judge the statistical precision of the false-alarm rate or its sensitivity to simulation choices.

    Authors: We acknowledge that including the total number of simulated pairs and statistical uncertainties is important for assessing the precision of our results. We will revise the manuscript to report the total number of event pairs considered in the simulation and provide error estimates on the false alarm rate, for example using Poisson statistics or bootstrap methods. Additionally, we will discuss the sensitivity of the results to variations in simulation parameters. revision: yes

Circularity Check

0 steps flagged

No significant circularity; FAR obtained from independent Monte Carlo simulation

full rationale

The central result (~0.01% false lensing classification rate) is produced by feeding an independently generated simulated population of unlensed BBH events into the GLANCE pipeline and directly counting triggers. This is a forward Monte Carlo measurement, not a fit to the target data nor a quantity defined in terms of itself. No step in the reported derivation reduces by construction to the inputs via self-definition, fitted-parameter renaming, or load-bearing self-citation chains. The simulation fidelity assumption affects correctness but does not create circularity in the derivation chain.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the fidelity of the simulated BBH population and the assumption that GLANCE's response to these synthetic events captures all relevant real-world systematics.

free parameters (2)
  • lensing SNR threshold
    Fixed at 1.5 to define a positive lensing classification in the false-alarm test.
  • minimum time delay
    Set near 1000 days to identify the regime where false positives appear.
axioms (1)
  • domain assumption The simulated merging binary black hole population statistically matches the true astrophysical distribution of sources, sky localizations, and noise properties.
    Invoked when generating the unlensed event catalog used to measure false positives.

pith-pipeline@v0.9.0 · 5836 in / 1460 out tokens · 57635 ms · 2026-05-18T07:15:54.607784+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

71 extracted references · 71 canonical work pages · 10 internal anchors

  1. [1]

    The noise can also produce artefacts that appear like substructures on the GW signal, thus posing as microlensing modulations

    Additionally, noise in the data can be similar to a GW signal, thus posing as strongly lensed pair. The noise can also produce artefacts that appear like substructures on the GW signal, thus posing as microlensing modulations. A review article on different effects that can cause hin- drance in the detection of a lensed GW (or a lensed pair), can be found ...

  2. [2]

    The waveform model for the GW used here is IMRPhenomPv2 [41, 42]

    Generation of GW events : We use GWSIM package [33], to get a set of BBH merger events given cosmology model and BBH population pa- rameters, given a 3 yrs observation period with 0.75 duty factors in a combination of three detectors: H1, L1 and V1. The waveform model for the GW used here is IMRPhenomPv2 [41, 42]. The 64s long data chunks are sampled at 2...

  3. [3]

    Selection of GW events : We use publicly- available O4 noise sensitivity curves of the H1, L1 and V1 detectors [43] to generate time-domain noise. We generate Gaussian noise in the time- domain and injected the signals of each event and calculate the matched-filter SNR 14 for all these events and select the events with matched-filter SNR ≥ 8 in at least t...

  4. [4]

    We estimated a set of three parameters RA, Dec and coalescence time to estimate their sky-localization error with nested sampler dynesty [44]

    Estimation of sky-localization error : We used BILBY for the selected events, to estimate their sky-localization errors. We estimated a set of three parameters RA, Dec and coalescence time to estimate their sky-localization error with nested sampler dynesty [44]. Since the characteristics (i.e. frequency evolution along with amplitude and phase evolution)...

  5. [5]

    sky-localization error of a pair of events

    Calculation of sky-overlap: We used a 300 x 150 uniform grid in the RA-Dec plane to find whether any of the grid points is lying inside the 90% C.I. sky-localization error of a pair of events. For the case, it has at least one grid point, we call these two events as a sky-overlapping event pair. We search for all possible event pairs to find the sky- over...

  6. [6]

    The RA-Dec points were also required to align the signals in different detectors

    Performing cross-correlation on sky- overlapping event pairs : We start with the sky-overlapping event pairs and reconstruct their polarizations from the overlapping sky-grid points and on a grid on the polarization angle containing 10 equispaced points between 0 and π. The RA-Dec points were also required to align the signals in different detectors. Data...

  7. [7]

    Abdikamalov, G

    E. Abdikamalov, G. Pagliaroli, and D. Radice, Gravita- tional waves from core-collapse supernovae, in Handbook of Gravitational Wave Astronomy (Springer Singapore,

  8. [8]

    Vartanyan, A

    D. Vartanyan, A. Burrows, T. Wang, M. S. B. Coleman, and C. J. White, Gravitational-wave signature of core- collapse supernovae, Phys. Rev. D 107, 103015 (2023), arXiv:2302.07092 [astro-ph.HE]

  9. [9]

    K. S. Thorne, Gravitational waves, in 1994 DPF Summer Study on High-energy Physics: Particle and Nuclear As- trophysics and Cosmology in the Next Millenium (Snow- mass 94) (1995) pp. 0160–184, arXiv:gr-qc/9506086

  10. [10]

    Belczynski, V

    K. Belczynski, V. Kalogera, and T. Bulik, A Compre- hensive study of binary compact objects as gravitational wave sources: Evolutionary channels, rates, and physical properties, Astrophys. J. 572, 407 (2001), arXiv:astro- ph/0111452

  11. [11]

    Riles, Searches for continuous-wave gravitational ra- diation, Living Rev

    K. Riles, Searches for continuous-wave gravitational ra- diation, Living Rev. Rel. 26, 3 (2023), arXiv:2206.06447 [astro-ph.HE]

  12. [12]

    Pagliaro, M

    G. Pagliaro, M. A. Papa, J. Ming, J. Lian, D. Tsuna, C. Maraston, and D. Thomas, Continuous Gravitational Waves from Galactic Neutron Stars: Demography, De- tectability, and Prospects, Astrophys. J.952, 123 (2023), arXiv:2303.04714 [gr-qc]

  13. [13]

    B. P. Abbott et al. (LIGO Scientific, Virgo), Observation of Gravitational Waves from a Binary Black Hole Merger, Phys. Rev. Lett. 116, 061102 (2016), arXiv:1602.03837 [gr-qc]

  14. [14]

    Abbott et al

    R. Abbott et al. (KAGRA, VIRGO, LIGO Scientific), Open Data from the Third Observing Run of LIGO, Virgo, KAGRA, and GEO, Astrophys. J. Suppl. 267, 29 (2023), arXiv:2302.03676 [gr-qc]

  15. [15]

    GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run

    R. Abbott et al. (KAGRA, VIRGO, LIGO Scien- tific), GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run, Phys. Rev. X 13, 041039 (2023), arXiv:2111.03606 [gr-qc]

  16. [16]

    Abbottet al.(LIGO Scientific, KAGRA, VIRGO), As- trophys

    R. Abbott et al. (LIGO Scientific, KAGRA, VIRGO), Observation of Gravitational Waves from Two Neutron Star–Black Hole Coalescences, Astrophys. J. Lett. 915, L5 (2021), arXiv:2106.15163 [astro-ph.HE]

  17. [17]

    A. G. Abac et al. (LIGO Scientific, Virgo, KAGRA), Ob- servation of Gravitational Waves from the Coalescence of a 2.5–4.5 M ⊙ Compact Object and a Neutron Star, Astrophys. J. Lett. 970, L34 (2024), arXiv:2404.04248 [astro-ph.HE]

  18. [18]

    B. P. Abbott et al. (LIGO Scientific, Virgo), GW170817: Observation of Gravitational Waves from a Binary Neu- tron Star Inspiral, Phys. Rev. Lett. 119, 161101 (2017), arXiv:1710.05832 [gr-qc]. 14

  19. [19]

    Y. Wang, A. Stebbins, and E. L. Turner, Gravitational lensing of gravitational waves from merging neutron star binaries, Phys. Rev. Lett. 77, 2875 (1996)

  20. [20]

    T. T. Nakamura, Gravitational lensing of gravitational waves from inspiraling binaries by a point mass lens, Phys. Rev. Lett. 80, 1138 (1998)

  21. [21]

    Bartelmann, Gravitational lensing, Classical and Quantum Gravity 27, 233001 (2010)

    M. Bartelmann, Gravitational lensing, Classical and Quantum Gravity 27, 233001 (2010)

  22. [22]

    L. Dai, T. Venumadhav, and K. Sigurdson, Effect of lens- ing magnification on the apparent distribution of black hole mergers, Phys. Rev. D 95, 044011 (2017)

  23. [23]

    Schneider, J

    P. Schneider, J. Ehlers, and E. E. Falco, Gravitational Lenses (Springer, 1992)

  24. [24]

    Chakraborty and S

    A. Chakraborty and S. Mukherjee, GLANCE – Grav- itational Lensing Authenticator using Non-modelled Cross-correlation Exploration of Gravitational Wave Sig- nals, Mon. Not. Roy. Astron. Soc. 532, 4842 (2024), arXiv:2403.03982 [gr-qc]

  25. [25]

    Chakraborty and S

    A. Chakraborty and S. Mukherjee, mu-GLANCE: A Novel Technique to Detect Chromatically and Achromat- ically Lensed Gravitational-wave Signals, Astrophys. J. 984, 107 (2025), arXiv:2410.06995 [gr-qc]

  26. [26]

    Chakraborty and S

    A. Chakraborty and S. Mukherjee, The first model- independent chromatic microlensing search: No evi- dence in the gravitational wave catalog of ligo-virgo-kagra (2025), arXiv:2503.16281 [gr-qc]

  27. [27]

    Martynov, E

    D. Martynov, E. Hall, B. Abbott, R. Abbott, T. Ab- bott, C. Adams, R. Adhikari, R. Anderson, S. Ander- son, K. Arai, M. Arain, S. Aston, L. Austin, S. Ballmer, M. Barbet, D. Barker, B. Barr, L. Barsotti, J. Bartlett, M. Barton, I. Bartos, J. Batch, A. Bell, I. Belopolski, J. Bergman, J. Betzwieser, G. Billingsley, J. Birch, S. Bis- cans, C. Biwer, E. Blac...

  28. [28]

    B. P. Abbott et al. (KAGRA, LIGO Scientific, Virgo, VIRGO), Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Ad- vanced Virgo and KAGRA, Living Rev. Rel.21, 3 (2018), arXiv:1304.0670 [gr-qc]

  29. [29]

    Buikema, C

    A. Buikema, C. Cahillane, G. Mansell, C. Blair, R. Abbott, C. Adams, R. Adhikari, A. Ananyeva, S. Appert, K. Arai, J. Areeda, Y. Asali, S. Aston, C. Austin, A. Baer, M. Ball, S. Ballmer, S. Banagiri, D. Barker, L. Barsotti, J. Bartlett, B. Berger, J. Bet- zwieser, D. Bhattacharjee, G. Billingsley, S. Biscans, R. Blair, N. Bode, P. Booker, R. Bork, A. Bram...

  30. [30]

    S.-S. Li, S. Mao, Y. Zhao, and Y. Lu, Gravitational lensing of gravitational waves: a statistical perspective, Monthly Notices of the Royal Astronomical Society 476, 2220 (2018)

  31. [31]

    Mukherjee, T

    S. Mukherjee, T. Broadhurst, J. M. Diego, J. Silk, and G. F. Smoot, Impact of astrophys- ical binary coalescence time-scales on the rate of lensed gravitational wave events, Monthly No- tices of the Royal Astronomical Society 506, 3751 (2021), https://academic.oup.com/mnras/article- pdf/506/3/3751/39456853/stab1980.pdf

  32. [32]

    Diego, T

    J. Diego, T. Broadhurst, and G. Smoot, Evidence for lensing of gravitational waves from ligo-virgo data, Phys- ical Review D 104, 10.1103/physrevd.104.103529 (2021)

  33. [33]

    Janquart, O

    J. Janquart, O. A. Hannuksela, K. Haris, and C. Van Den Broeck, A fast and precise methodology to search for and analyse strongly lensed gravitational-wave events, Monthly Notices of the Royal Astronomical Society 506, 5430 (2021), https://academic.oup.com/mnras/article- pdf/506/4/5430/39656596/stab1991.pdf

  34. [34]

    Wright and M

    M. Wright and M. Hendry, Gravelamps: Gravitational wave lensing mass profile model selection, The Astro- physical Journal 935, 68 (2022)

  35. [35]

    R. K. L. Lo and I. Maga˜ na Hernandez, Bayesian statistical framework for identifying strongly lensed gravitational-wave signals, Phys. Rev. D 107, 123015 (2023)

  36. [36]

    E. Seo, X. Shan, J. Janquart, O. A. Hannuksela, M. A. Hendry, and B. Hu, Residual Test to Search for Microlensing Signatures in Strongly Lensed Gravi- tational Wave Signals, Astrophys. J. 988, 159 (2025), arXiv:2503.02186 [gr-qc]

  37. [37]

    Abbottet al.(LIGO Scientific, KAGRA, VIRGO), Astrophys

    R. Abbott et al. (LIGO Scientific, KAGRA, VIRGO), Search for Gravitational-lensing Signatures in the Full Third Observing Run of the LIGO–Virgo Network, As- trophys. J. 970, 191 (2024), arXiv:2304.08393 [gr-qc]

  38. [38]

    Janquart, M

    J. Janquart et al. , Follow-up analyses to the O3 LIGO–Virgo–KAGRA lensing searches, Mon. Not. Roy. Astron. Soc. 526, 3832 (2023), arXiv:2306.03827 [gr-qc]

  39. [39]

    Karathanasis, B

    C. Karathanasis, B. Revenu, S. Mukherjee, and F. Stachurski, GWSim: Python package for creating mock GW samples for different astrophysical popula- tions and cosmological models of binary black holes, Astron. Astrophys. 677, A124 (2023), [Erratum: As- tron.Astrophys. 682, C1 (2024)], arXiv:2210.05724 [astro- ph.CO]

  40. [40]

    Bilby: A user-friendly Bayesian inference library for gravitational-wave astronomy

    G. Ashton et al. , BILBY: A user-friendly Bayesian infer- ence library for gravitational-wave astronomy, Astrophys. J. Suppl. 241, 27 (2019), arXiv:1811.02042 [astro-ph.IM]

  41. [41]

    D. W. Hogg, Distance measures in cosmology (2000), arXiv:astro-ph/9905116 [astro-ph]

  42. [42]

    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

  43. [43]

    Morr´ as, J

    G. Morr´ as, J. F. N. n. Siles, J. Garc´ ıa-Bellido, and E. R. Morales, False alarms induced by gaussian noise in gravitational wave detectors, Phys. Rev. D 107, 023027 (2023)

  44. [44]

    Keitel, False positives for gravitational lensing: the gravitational-wave perspective, Phil

    D. Keitel, False positives for gravitational lensing: the gravitational-wave perspective, Phil. Trans. Roy. Soc. Lond. A 383, 20240128 (2025), arXiv:2407.12974 [gr-qc]

  45. [45]

    Mukherjee, T

    S. Mukherjee, T. Broadhurst, J. M. Diego, J. Silk, and G. F. Smoot, Impact of astrophysical binary coales- cence time-scales on the rate of lensed gravitational wave events, Mon. Not. Roy. Astron. Soc. 506, 3751 (2021), arXiv:2106.00392 [gr-qc]

  46. [46]

    Tyson, Mapping dark matter with gravi- tational lenses, Physics Today 45, 24 (1992), https://pubs.aip.org/physicstoday/article- pdf/45/6/24/11176778/24 1 online.pdf

    A. Tyson, Mapping dark matter with gravi- tational lenses, Physics Today 45, 24 (1992), https://pubs.aip.org/physicstoday/article- pdf/45/6/24/11176778/24 1 online.pdf

  47. [47]

    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. D 93, 044006 (2016), arXiv:1508.07250 [gr-qc]

  48. [48]

    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. D 93, 044007 (2016), arXiv:1508.07253 [gr-qc]

  49. [49]

    collaboration, Noise curves used for simulations in the update of the observing scenarios paper (2022)

    L. collaboration, Noise curves used for simulations in the update of the observing scenarios paper (2022)

  50. [50]

    J. S. Speagle, dynesty: a dynamic nested sampling pack- age for estimating bayesian posteriors and evidences, Monthly Notices of the Royal Astronomical Society 493, 3132–3158 (2020)

  51. [51]

    S. R. Shukla, L. Pathak, and A. S. Sengupta, How i won- der where you are: pinpointing coalescing binary neutron star sources with the igwn, including ligo-aundha (2023), arXiv:2311.15695 [gr-qc]

  52. [52]

    E. D. Hall, Cosmic explorer: A next-generation ground-based gravitational-wave observatory, Galaxies 10, 10.3390/galaxies10040090 (2022)

  53. [53]

    Punturo et al

    M. Punturo et al. , The Einstein Telescope: A third- generation gravitational wave observatory, Class. Quant. Grav. 27, 194002 (2010)

  54. [54]

    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. D 107, 043029 (2023), arXiv:2206.02803 [astro-ph.CO]

  55. [55]

    C. R. Harris, K. J. Millman, S. J. van der Walt, R. Gom- mers, P. Virtanen, D. Cournapeau, E. Wieser, J. Tay- lor, S. Berg, N. J. Smith, R. Kern, M. Picus, S. Hoyer, M. H. van Kerkwijk, M. Brett, A. Haldane, J. F. del R´ ıo, M. Wiebe, P. Peterson, P. G´ erard-Marchant, K. Shep- pard, T. Reddy, W. Weckesser, H. Abbasi, C. Gohlke, and T. E. Oliphant, Array ...

  56. [56]

    Virtanen, R

    P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Pe- terson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙I. Po- 16 lat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. H...

  57. [57]

    J. D. Hunter, Matplotlib: A 2d graphics environment, Computing in Science & Engineering 9, 90 (2007)

  58. [58]

    A. Nitz, I. Harry, D. Brown, C. M. Biwer, J. Willis, T. D. Canton, C. Capano, T. Dent, L. Pekowsky, G. S. C. Davies, S. De, M. Cabero, S. Wu, A. R. Williamson, B. Machenschalk, D. Macleod, F. Pannarale, P. Kumar, S. Reyes, dfinstad, S. Kumar, M. Tapai, L. Singer, P. Ku- mar, veronica villa, maxtrevor, B. U. V. Gadre, S. Khan, S. Fairhurst, and A. Tolley, ...

  59. [59]

    D. M. Macleod, J. S. Areeda, S. B. Coughlin, T. J. Massinger, and A. L. Urban, GWpy: A Python pack- age for gravitational-wave astrophysics, SoftwareX 13, 100657 (2021)

  60. [60]

    LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration, LVK Algorithm Library - LAL- Suite, Free software (GPL) (2018)

  61. [61]

    Cosmic Star Formation History

    P. Madau and M. Dickinson, Cosmic Star Formation History, Ann. Rev. Astron. Astrophys. 52, 415 (2014), arXiv:1403.0007 [astro-ph.CO]. Appendix A: The rate of detectable lensed GW events: term-by term explanation In this appendix, we summarize the factors deciding the rate of detectable lensed events, presented in the equation 17:

  62. [62]

    Here ⃗θ consists all the source intrinsic properties that account for a GW

    The first term in the second integral, P (⃗θ) is the joint probability density of⃗θ. Here ⃗θ consists all the source intrinsic properties that account for a GW. Here we have considered the dependence only on the source masses and spins ⃗θ = (m1s, m2s, ⃗ χ1, ⃗ χ2). Note that, R P (⃗θ)d⃗θ = 1

  63. [63]

    The selection function has a maximum value of unity with a minimum of zero

    The second term in the second integral ,S(⃗θ, µ, z) is the selection function that decides whether an event is detectable or not given a lensing magnification of µ, source luminosity distance dl (which corresponds to a source redshift z) and the source masses and spins (m1s, m2s, ⃗ χ1, ⃗ χ2). The selection function has a maximum value of unity with a mini...

  64. [64]

    The first term in the first integral R(z) is the merger density rate (in units of number of events / unit time / unit volume) as a function of the red- shift. Here R0 is the merger density rate at z ≈ 0, zp is the peak redshift of the merger density andmα and mβ constants that governs the nature of the merger rate density at low redshift (upto z ≈ zp) and...

  65. [65]

    Here dV can be considered as the comoving volume of the spherical shell at redshiftz with thickness dz

    The second term of the first integral dV /dz is the differential comoving volume which tells how the volume of a sphere up to a redshift z varies with z. Here dV can be considered as the comoving volume of the spherical shell at redshiftz with thickness dz. Therefore by chain rule dV = dV dz dz produces the volume element to be integrated over in the firs...

  66. [66]

    The lensing optical depth de- scribes the probability of an GW at redshiftz being lensed with a magnification µ

    The third term in the first integral τ(≥ µ, z) tells about the lensing optical depth of the universe in variation with z. The lensing optical depth de- scribes the probability of an GW at redshiftz being lensed with a magnification µ. Therefore, performing the first integral after the second integral without τ(≥ µ, z) and µ = 1, we obtain the de- tectable...

  67. [67]

    For the BBH merger rate redshift evolution, we have used the merger-rate equation based on the empirical Madau-Dickinson star-formation rate [55] given by, R(z) = R0(1 + z)mα 1 + (1 +zp)−(mα+mβ) 1 + 1+z 1+zp (mα+mβ)

  68. [68]

    dVc dz is the differential comoving volume at a redshift z. 17

  69. [69]

    The distributions are given by, P (m1 | Mmin, Mmax, α) ∝ m−α 1 Mmin < m1 < Mmax 0 otherwise and, G (m1 | µ, σ) = 1 σ √ 2π exp −(m1 − µ)2 2σ2 !

    For the BBH mass-distribution model, we choose primary masses from a power-law plus Gaussian distribution given by 19, p (m1 | Mmin, Mmax, α, µ, σ, λ, δ) ∝ S (m1 | Mmin, δ) ×    (1 − λ) P (m1 | Mmax, Mmin, α) +λG (m1 | µ, σ) for Mmin < m1 < Mmax, 0 otherwise Here, P (m1 | Mmax, Mmin, α) is the truncated power law distribution and G (m1 | µ, σ) is a G...

  70. [70]

    Here z is the direction of the angu- lar momentum vector

    The magnitude of the BBH spins are chosen with the help of effective spin parameter χeff = m1χ1z+m1χ1z m1+m2 . Here z is the direction of the angu- lar momentum vector. Therefore, the effective spin χeff tells us the sum of the mass-weighted aligned spins of the binary black hole. The value of χeff are chosen from the uniform distribution given by, P (χef...

  71. [71]

    In the table B, we have mentioned our choice for differ- ent parameters for the generation of a realistic BBH pop- ulation using GWSIM

    The selection function S(⃗θ, µ, z) selects those events above magnified by a flux-magnification factor µ with a signal-to-noise ratio ( ρ) of 8 in at least two of the detectors. In the table B, we have mentioned our choice for differ- ent parameters for the generation of a realistic BBH pop- ulation using GWSIM. The figure 10, shows the distribu- tion for...