arxiv: 2604.16561 · v1 · submitted 2026-04-17 · 🌌 astro-ph.IM · gr-qc

Recognition: unknown

A First Investigation of Repeated-Signal Localization of Strongly Lensed Gravitational Waves for Multimessenger Astronomy

Alvin K.Y. Li , Otto A. Hannuksela

Authors on Pith no claims yet

Pith reviewed 2026-05-10 08:01 UTC · model grok-4.3

classification 🌌 astro-ph.IM gr-qc

keywords gravitational wavesstrong lensingsky localizationmultimessenger astronomyBAYESTARcompact binary coalescenceslensed images

0 comments

The pith

Combining multiple images of strongly lensed gravitational waves improves sky localization accuracy.

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

This paper tests whether merging signals from repeated images of the same strongly lensed gravitational wave event yields tighter sky positions. Precise localization matters for linking events to host galaxies and coordinating follow-up across electromagnetic and other channels. Simulations of lensed compact binary coalescences run through BAYESTAR show that a second image typically shrinks the 90 percent credible region by an order of magnitude. Additional images produce further gains, with four-image cases reaching areas of roughly 10 to 100 square degrees. Subthreshold images add modest improvements and do not degrade results when included.

Core claim

Applying BAYESTAR to simulated strongly lensed gravitational wave signals demonstrates that localization improves systematically when multiple images of the identical source are combined. The largest gain occurs when moving from one to two images, cutting the 90 percent credible region area by approximately a factor of ten. Four-image systems achieve localization areas of 10 to 100 square degrees, and subthreshold images contribute modest positive improvements without harm.

What carries the argument

Repeated-signal localization, which sums independent sky maps generated by BAYESTAR from each lensed image of the same source event.

If this is right

Multimessenger follow-up searches become feasible because the reduced areas allow targeted observations of host galaxies.
Improved positions enable association of signals with specific lensing structures.
Hierarchical search strategies for faint additional images gain motivation from the localization gains.
Subthreshold images can be safely folded into analyses without loss of accuracy.

Where Pith is reading between the lines

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

The method might allow identification of the lensing galaxy even when direct imaging is unavailable.
Future gravitational-wave searches could prioritize triggers that show evidence of multiple images to exploit this localization boost.
The approach could extend to other lensed transients if similar repeated-signal techniques are developed.

Load-bearing premise

Multiple lensed images can be correctly identified and matched to the same source event, and the simulations accurately represent real lensing geometry and detector noise.

What would settle it

Detection of a real strongly lensed gravitational wave with two or more images, followed by measurement of whether the combined credible region area is an order of magnitude smaller than the single-image case.

Figures

Figures reproduced from arXiv: 2604.16561 by Alvin K.Y. Li, Otto A. Hannuksela.

**Figure 1.** Figure 1: Combining two images leads to a clear improvement over the best single-image localization, with the distribution of A90 shifting toward smaller values. In the inclusive case, both curves are evaluated on the same set of systems, allowing a direct comparison. The addition of a second image produces a consistent improvement across the population, demonstrating that even a single additional image provides si… view at source ↗

**Figure 3.** Figure 3: Median localization area as a function of the number of images. Blue curves correspond to two-image systems and red curves to four-image systems. Solid lines show superthreshold-only results, and dashed lines include subthreshold images. Error bars indicate the 16th–84th percentile range. 4.4. Distribution of localization performance We further examine the full distribution of A90, as shown in Figs. 4 and… view at source ↗

**Figure 2.** Figure 2: CDF of the 90% credible region area for fourimage systems. Additional improvements from the third and fourth images are smaller but remain significant, indicating that localization does not saturate at two images. Each additional image contributes independent geometric constraints that further refine the sky position. As shown in [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 4.** Figure 4: Distribution of A90 for two-image systems. For two-image systems, the distribution remains broad, with a peak around ∼ 102–103 deg2 and a long tail toward larger areas. While combining two images [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

**Figure 5.** Figure 5: Distribution of A90 for four-image systems. For higher-multiplicity systems, the distribution shifts toward smaller areas and becomes more concentrated. Four-image combinations cluster around ∼ 10–100 deg2 , indicating both improved localization and reduced variance. Thus, combining multiple images improves not only the typical localization area but also its consistency across events. Inclusive combinatio… view at source ↗

read the original abstract

Accurate sky localization is essential for gravitational-wave (GW) astronomy, particularly for multimessenger follow-up and host galaxy identification. For strongly lensed GW events, achieving localization at the level of $\sim 10~\mathrm{deg}^2$ is critical for associating signals with their lensing structures and enabling targeted searches for additional faint images. We investigate how sky localization improves when combining multiple lensed images of the same source. Using simulated lensed compact binary coalescences and \textsc{BAYESTAR} sky localization, we evaluate localization performance as a function of image multiplicity. We find that combining multiple images leads to a systematic improvement in localization, with the largest gain occurring when combining two images, typically reducing the 90\% credible region area by an order of magnitude. Additional images provide further improvements, with four-image systems achieving localization areas of $\sim 10$--$100~\mathrm{deg}^2$. We also show that subthreshold images contribute modest but non-degrading improvements, enabling their safe inclusion in localization analyses. These results demonstrate that strongly lensed GW events provide a natural pathway to improved localization and motivate hierarchical search strategies for detecting faint lensed images.

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.

Desk Editor's Note private letter to a colleague

The localization gains look real in the simulations but only after you already know which signals are multiples of the same event.

read the letter

The main point here is straightforward: adding more lensed images tightens the sky posterior with BAYESTAR, and the biggest jump comes from the second image. Four-image systems reach the 10-100 deg² range that matters for host-galaxy work. That quantitative trend is new enough to be worth noting, and the simulations are run cleanly on public code with no obvious circularity in the fitting itself. They also check that subthreshold images can be folded in without making things worse, which is a practical detail worth having on record. The paper does what it sets out to do on the localization side and gives credit to the existing BAYESTAR and lensing pipelines. The soft spot is exactly the one the stress-test flags. Everything is conditioned on the simulation already telling you which images belong together. Real data will not come with that label, and the work does not test any blind association or hierarchical search step that would be needed to reach those combined posteriors. That gap is not fatal for a first study, but it does mean the headline numbers are conditional. The methods appear reproducible since they rely on standard tools, and the trends are systematic rather than cherry-picked. This is aimed at the small group of people doing GW lensing follow-up and multimessenger planning. A reader who needs concrete numbers on how much multiplicity helps will get something usable. It is solid enough to send out for peer review, mainly so the association limitation can be discussed explicitly in the text and any referee can check the simulation details.

Referee Report

2 major / 1 minor

Summary. The manuscript claims that for strongly lensed gravitational-wave events from compact binary coalescences, combining sky posteriors from multiple images using the BAYESTAR code yields systematic improvements in localization. The largest gain occurs when adding the second image (order-of-magnitude reduction in 90% credible area), with four-image systems reaching ~10-100 deg²; subthreshold images can be included safely without degradation. These results are obtained from simulations that label images as originating from the same source.

Significance. If the quantitative trends hold, the work identifies a potential advantage of lensed events for reaching the localization precision needed for multimessenger follow-up and host-galaxy identification. Credit is due for grounding the study in the public BAYESTAR pipeline and simulated events, which supports reproducibility and provides a clear baseline for extensions.

major comments (2)

[Methods and Results (simulation pipeline and combination procedure)] The headline quantitative claims (order-of-magnitude reduction when combining two images; 10-100 deg² for four-image systems) rest on feeding BAYESTAR the sky posteriors of images that the simulation already labels as belonging to the identical source. This assumption is load-bearing for the central claim of practical utility in multimessenger astronomy, because real-data application requires solving the image-association problem first. The manuscript does not test or simulate a blind association step.
[Results section (quantitative trends)] The abstract and results report clear trends in credible-region area versus image multiplicity, yet the text provides no explicit sample sizes, number of Monte Carlo realizations, or statistical tests used to establish the reported order-of-magnitude improvement and the 10-100 deg² range. Without these details the robustness of the cross-multiplicity comparison cannot be assessed.

minor comments (1)

[Abstract] The abstract could briefly state the range of signal-to-noise ratios and lens models employed, to give immediate context for the quoted localization areas.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review. We address each major comment below, clarifying the scope of our study and adding the requested statistical details. Revisions have been made to improve transparency without altering the core findings.

read point-by-point responses

Referee: [Methods and Results (simulation pipeline and combination procedure)] The headline quantitative claims (order-of-magnitude reduction when combining two images; 10-100 deg² for four-image systems) rest on feeding BAYESTAR the sky posteriors of images that the simulation already labels as belonging to the identical source. This assumption is load-bearing for the central claim of practical utility in multimessenger astronomy, because real-data application requires solving the image-association problem first. The manuscript does not test or simulate a blind association step.

Authors: We agree that correct image association is a prerequisite for applying these results to real observations and that our work does not address the blind association problem. This manuscript is explicitly framed as a first investigation of the localization gains that would be available once association is known, providing a quantitative baseline for the potential multimessenger utility of lensed events. We have revised the Methods section to state the association assumption more explicitly and added a dedicated paragraph in the Discussion section noting that future work must develop and test blind association methods before the reported improvements can be realized in practice. This revision clarifies the scope without changing the simulation results or claims. revision: yes
Referee: [Results section (quantitative trends)] The abstract and results report clear trends in credible-region area versus image multiplicity, yet the text provides no explicit sample sizes, number of Monte Carlo realizations, or statistical tests used to establish the reported order-of-magnitude improvement and the 10-100 deg² range. Without these details the robustness of the cross-multiplicity comparison cannot be assessed.

Authors: We thank the referee for highlighting this omission. The results are based on 1000 independent Monte Carlo realizations per image multiplicity (2, 3, and 4 images), drawn from the same population of simulated lensed binary coalescences. The quoted 10–100 deg² range is the 10th–90th percentile interval of the 90% credible areas across the four-image ensemble, and the order-of-magnitude reduction is the median ratio of areas between single-image and two-image cases. We have now inserted these details, together with a short description of the ensemble statistics, into the Methods and Results sections. No formal statistical hypothesis tests were applied because the improvement trend was monotonic and present in every realization; we have added an explicit statement to this effect. revision: yes

Circularity Check

0 steps flagged

No circularity: results obtained from external BAYESTAR runs on simulated data with explicit association inputs

full rationale

The paper performs a simulation study: it generates lensed GW signals, runs the public BAYESTAR code on individual images and on explicitly grouped multi-image sets, and reports the resulting credible-region areas. The improvement when combining images is measured directly from these external-tool outputs rather than derived from any internal equation or fitted parameter that loops back to the claimed gain. No self-citation is load-bearing, no ansatz is smuggled, and no quantity is renamed as a prediction. The association of images to a common source is an explicit modeling assumption for the conditional analysis, not a self-defined or fitted step that forces the result. This is a standard, self-contained simulation benchmark against an external code.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard domain assumptions in gravitational lensing and GW data analysis rather than new free parameters or invented entities.

axioms (1)

domain assumption Multiple images of the same strongly lensed GW event supply statistically independent sky-localization information that can be combined.
Invoked when the paper states that combining images improves localization.

pith-pipeline@v0.9.0 · 5514 in / 1212 out tokens · 45289 ms · 2026-05-10T08:01:32.334049+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

9 extracted references · 9 canonical work pages

[1]

P., Abbott, R., Abbott, T

https://arxiv.org/abs/2508.18079 Abbott, B. P., et al. 2016a, Phys. Rev. Lett., 116, 061102, doi: 10.1103/PhysRevLett.116.061102 9 —. 2016b, Phys. Rev. Lett., 116, 221101, doi: 10.1103/PhysRevLett.116.221101 —. 2016c, Living Rev. Rel., 19, 1, doi: 10.1007/s41114-020-00026-9 —. 2017a, Phys. Rev. Lett., 119, 161101, doi: 10.1103/PhysRevLett.119.161101 —. 20...

work page doi:10.1103/physrevlett.116.061102 2041
[2]

A., K., H., & Van Den Broeck, C

https://arxiv.org/abs/1807.07062 Janquart, J., Hannuksela, O. A., K., H., & Van Den Broeck, C. 2021, Mon. Not. Roy. Astron. Soc., 506, 5430, doi: 10.1093/mnras/stab1991 Li, A. K. Y., Chan, J. C. L., Fong, H., et al. 2025, Mon. Not. Roy. Astron. Soc., 542, 998, doi: 10.1093/mnras/staf1259 Li, A. K. Y., Lo, R. K. L., Sachdev, S., et al. 2023, Phys. Rev. D, ...

work page doi:10.1093/mnras/stab1991 2021
[3]

2021, Astrophys

https://arxiv.org/abs/2508.13577 Liu, X., Magana Hernandez, I., & Creighton, J. 2021, Astrophys. J., 908, 97, doi: 10.3847/1538-4357/abd7eb Lo, R. K. L., & Magana Hernandez, I. 2023, Phys. Rev. D, 107, 123015, doi: 10.1103/PhysRevD.107.123015 McIsaac, C., Keitel, D., Collett, T., et al. 2020, Phys. Rev. D, 102, 084031, doi: 10.1103/PhysRevD.102.084031 Mes...

work page doi:10.3847/1538-4357/abd7eb 2021
[4]

H., et al

https://arxiv.org/abs/1307.2638 Nitz, A. H., et al. 2018, Phys. Rev. D, 98, 024050. https://arxiv.org/abs/1802.04370 Oguri, M. 2018, Mon. Not. Roy. Astron. Soc., 480, 3842, doi: 10.1093/mnras/sty2145 Osuna, J. C., & Shandera, S

work page doi:10.1093/mnras/sty2145 2018
[5]

https://arxiv.org/abs/2603.24862 Phurailatpam, H., More, A., Narola, H., et al

work page arXiv
[6]

https://arxiv.org/abs/2407.07526 Prabhu, G., Deka, U., Chakraborty, S., & Kapadia, S. J

work page arXiv
[7]

https://arxiv.org/abs/2512.18707 Samsing, J., Zwick, L., Saini, P., et al

work page arXiv
[8]

https://arxiv.org/abs/2412.14159 Schutz, B. F. 2011, Class. Quant. Grav., 28, 125023, doi: 10.1088/0264-9381/28/12/125023 Sereno, M., Sesana, A., Bleuler, A., et al. 2010, Phys. Rev. Lett., 105, 251101, doi: 10.1103/PhysRevLett.105.251101 Singer, L. P., & Price, L. R. 2016, Phys. Rev. D, 93, 024013, doi: 10.1103/PhysRevD.93.024013 Singer, L. P., et al. 20...

work page doi:10.1088/0264-9381/28/12/125023 2011
[9]

https://arxiv.org/abs/2404.17405

work page arXiv