pith. machine review for the scientific record. sign in

arxiv: 2604.19073 · v1 · submitted 2026-04-21 · 🌌 astro-ph.HE · astro-ph.IM· gr-qc

Recognition: unknown

A comprehensive framework for phase-coherent mapping of the gravitational-wave sky with pulsar timing arrays

Authors on Pith no claims yet

Pith reviewed 2026-05-10 02:27 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.IMgr-qc
keywords pulsar timing arraysgravitational wavessky mappingphase-coherent analysisstochastic backgroundanisotropysource identificationpolarisation
0
0 comments X

The pith

Phase-coherent mapping preserves amplitude, phase and polarization of gravitational waves in every sky pixel from pulsar timing arrays.

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

The paper introduces a practical phase-coherent mapping technique for pulsar timing arrays that resolves the full complex polarisation state of the gravitational-wave sky as a function of direction and frequency. Unlike standard cross-correlation methods, this approach keeps the signal amplitude, phase and polarisation intact in each sky pixel. The resulting maps form a compact summary of the data from which characterisation of a stochastic background, searches for anisotropy and identification of individual sources can all be derived in one unified framework. The method remains compatible with established pulsar timing data analysis techniques and has been tested through simulations that vary array configurations, noise properties and signal types to show robust recovery of source amplitudes and locations.

Core claim

This framework resolves the full complex polarisation state of the gravitational-wave sky as a function of direction and frequency, producing maps that preserve amplitude, phase and polarisation in every sky pixel and serve as a minimally processed summary from which all subsequent analyses follow within a single framework.

What carries the argument

Phase-coherent mapping that resolves the full complex polarisation state of the gravitational-wave sky and preserves amplitude, phase and polarisation in every sky pixel.

If this is right

  • Characterisation of a stochastic background, searches for anisotropy and identification of individual sources can all be performed from the same set of maps.
  • The maps act as a compact, minimally processed summary so that later analyses do not require reprocessing the original pulsar timing data.
  • The technique remains fully compatible with established pulsar timing data analysis methods.
  • Recovery of source amplitudes and sky locations works across different array configurations, noise levels and signal types in simulations.

Where Pith is reading between the lines

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

  • The unified maps could allow direct combination of data from multiple pulsar timing arrays without separate reprocessing steps for each analysis type.
  • Direction-dependent sensitivity and polarisation leakage effects could be studied quantitatively by inspecting the maps pixel by pixel rather than through separate statistical tests.
  • Once produced, the maps might serve as input for machine-learning searches for rare or transient signals that current cross-correlation pipelines miss.

Load-bearing premise

The framework assumes that realistic simulations with varying array configurations, noise properties and signal types accurately represent real pulsar timing data and that robust recovery of amplitudes and locations remains possible despite polarization leakage and direction-dependent sensitivity.

What would settle it

A simulation in which the method fails to recover the known amplitude or sky location of an injected gravitational-wave signal in realistic pulsar timing data would show that the claimed robust recovery does not hold.

Figures

Figures reproduced from arXiv: 2604.19073 by Dawson S. Gaynor, Eric Thrane, Ma{\l}gorzata Cury{\l}o, Paul D. Lasky.

Figure 1
Figure 1. Figure 1: Radiometer (left) and clean (right) maps of the real components of two polarisation states h+ (top) and h× (bottom). The injected continuous GW source is defined to emit entirely in Re(h+) and any strong signal in other channels constitutes leakage (bottom-left). The black hexagon indicates the injected source position, whereas the dark circle shows the recovered high S/N area. trix, which accounts for the… view at source ↗
Figure 3
Figure 3. Figure 3: Total power radiometer and clean maps for the MPTA simulation, demonstrating the need for regularisation when localising the source. The dark circle in both maps shows the highest S/N region (here it is again a single pixel) with true source position denoted with a black hexagon. In the case of a highly non-isotropic PTA array such as the MPTA, the position of the continuous GW source is critical, which we… view at source ↗
Figure 2
Figure 2. Figure 2: Total power maps for the IPTA simulation with a continuous GW signal. The dark circle in each plot cor￾responds to highest S/N region (in this case single pixel) identified in the clean total power map. The black hexagon indicated the true position of the injected GW source. (a) Radiometer S/N. (b) Clean S/N [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: Clean maps of the real component of h+ polarisation state for the MPTA simulation. Each map shows the same continuous GW source injected in different sky locations, demonstrating differences in directional sensitivity of the array. sky localization of continuous GW sources by effectively turning each Earth-pulsar baseline into an independent directional measurement (see e.g., Tsai et al. 2025). In our pipe… view at source ↗
Figure 5
Figure 5. Figure 5: Clean S/N maps of the real h+ polarisation for three types of signals: a) noise only, b) GW background, and c) two identical continuous GW sources. spanning multiple frequency bins. We will present de￾tailed studies of spectral leakage and a continuous GW parameter estimation in coherently combined maps in an upcoming follow up paper. 5.3. Connection to the isotropic analysis In the standard isotropic anal… view at source ↗
Figure 6
Figure 6. Figure 6: Clean S/N maps for an isotropic GW background at the first Fourier-frequency bin for the IPTA simulation. REFERENCES Abadie, J., Abbott, B. P., Abbott, R., et al. 2011, PhRvL, 107, 271102, doi: 10.1103/PhysRevLett.107.271102 Abbott, B., Abbott, R., Adhikari, R., et al. 2008, PhRvD, 77, 022001, doi: 10.1103/PhysRevD.77.022001 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016a, PhRvL, 116, 061102, doi: 1… view at source ↗
Figure 7
Figure 7. Figure 7: Impact of the unmodelled pulsar term on the clean map S/N for two sky locations in the MPTA-like array. Top: best-covered location (bottom-right), where the coherent earth term dominates and the source is correctly recovered. Bottom: poorly-covered location (upper-left), where the incoherent pulsar term power is comparable to the Earth-term signal, degrading both localisation and amplitude recovery [PITH_… view at source ↗
read the original abstract

We present a practical implementation of a phase-coherent mapping technique for pulsar timing arrays that resolves the full complex polarisation state of the gravitational-wave sky as a function of direction and frequency. Unlike standard cross-correlation methods, this approach preserves the amplitude, phase, and polarisation of the signal in every sky pixel. The resulting maps constitute a compact, minimally processed summary of the data from which all subsequent analyses -- characterisation of a stochastic background, searches for anisotropy, and identification of individual sources -- can be derived within a single unified framework. Our implementation is fully compatible with established pulsar timing data analysis methods. We validate the framework through a series of realistic simulations with varying array configurations, noise properties, and signal types. We demonstrate robust recovery of source amplitudes and sky locations across different scenarios, and discuss the impact of polarisation leakage, noise, and direction-dependent array sensitivity on the recovery of astrophysical signals.

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 presents a practical implementation of a phase-coherent mapping technique for pulsar timing arrays that resolves the full complex polarisation state of the gravitational-wave sky as a function of direction and frequency. Unlike standard cross-correlation methods, the approach constructs sky maps by solving for complex amplitudes in each pixel from timing residuals, preserving amplitude, phase, and polarisation information. These maps are claimed to form a compact, minimally processed summary of the data from which stochastic background characterisation, anisotropy searches, and individual source identification can all be derived within a single unified framework. The implementation is compatible with established PTA analysis methods and is validated through realistic simulations with varying array configurations, noise properties, and signal types, demonstrating robust recovery of source amplitudes and sky locations while discussing polarisation leakage, noise, and direction-dependent sensitivity.

Significance. If the invertibility of the mapping holds under realistic conditions, this framework would represent a meaningful advance for PTA gravitational-wave astronomy by supplying a single, information-rich data product that supports multiple downstream science cases without repeated reprocessing of the raw timing residuals. It could improve consistency and efficiency in analyses of upcoming large PTA datasets, particularly for detecting anisotropic signals or resolving individual sources, provided the retained phase and polarisation information survives the finite-array degeneracies.

major comments (2)
  1. [Abstract and §5] Abstract and §5 (validation): the central claim that the maps constitute a 'compact, minimally processed summary' from which all subsequent analyses can be derived rests on the assumption that the mapping operator is invertible without significant information loss. With a finite number of pulsars the response matrix is necessarily rank-deficient; the simulations demonstrate recovery inside the chosen ensemble but supply no analytic bound on the condition number of the operator or on residual mixing between pixels or Stokes parameters due to polarisation leakage and direction-dependent sensitivity.
  2. [§3] §3 (mapping construction): the procedure of solving for complex amplitudes per pixel does not appear to include an explicit treatment of the null-space projection or regularisation; without this, it is unclear whether phase and polarisation information that is lost in the initial mapping can be recovered in later analyses such as individual-source searches.
minor comments (2)
  1. The abstract would be strengthened by including at least one quantitative metric (e.g., fractional amplitude recovery error or sky-location uncertainty) from the simulations rather than qualitative statements of 'robust recovery'.
  2. Notation for the complex amplitudes, polarisation basis, and frequency binning should be defined explicitly at first use to aid readers unfamiliar with the PTA mapping literature.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on the manuscript. We address each major comment below, providing clarifications and indicating the revisions we will make.

read point-by-point responses
  1. Referee: [Abstract and §5] Abstract and §5 (validation): the central claim that the maps constitute a 'compact, minimally processed summary' from which all subsequent analyses can be derived rests on the assumption that the mapping operator is invertible without significant information loss. With a finite number of pulsars the response matrix is necessarily rank-deficient; the simulations demonstrate recovery inside the chosen ensemble but supply no analytic bound on the condition number of the operator or on residual mixing between pixels or Stokes parameters due to polarisation leakage and direction-dependent sensitivity.

    Authors: We agree that the absence of an analytic bound on the condition number leaves the invertibility claim somewhat qualitative. Deriving a general closed-form bound is difficult because the response matrix depends on the specific pulsar positions, frequencies, and noise realizations. In the revised manuscript we will add a dedicated paragraph in §5 that reports the numerically computed condition numbers for all simulated array configurations, together with quantitative measures of pixel-to-pixel and Stokes-parameter mixing recovered from the simulations. These numerical diagnostics will be used to qualify the statement that the maps form a compact summary suitable for downstream analyses, while explicitly noting the information loss inherent to finite arrays. revision: yes

  2. Referee: [§3] §3 (mapping construction): the procedure of solving for complex amplitudes per pixel does not appear to include an explicit treatment of the null-space projection or regularisation; without this, it is unclear whether phase and polarisation information that is lost in the initial mapping can be recovered in later analyses such as individual-source searches.

    Authors: The referee is correct that the current description in §3 does not explicitly discuss null-space handling. The linear solve is performed via a truncated singular-value decomposition in which singular values below a chosen threshold (set by the noise level) are discarded; this implicitly projects out the null space. We will revise §3 to state this procedure clearly, specify the threshold criterion, and add a short discussion of how the retained modes preserve the recoverable phase and polarisation information. We will also include a brief demonstration, using the existing simulation suite, that individual-source searches performed on the resulting maps recover the injected signals with fidelity comparable to direct timing-residual analyses. revision: yes

Circularity Check

0 steps flagged

No circularity: mapping framework is a direct implementation validated by simulation

full rationale

The paper introduces a phase-coherent mapping technique that solves for complex amplitudes per sky pixel from timing residuals. This construction is presented as a new implementation compatible with standard PTA methods, with validation performed through independent realistic simulations covering array configurations, noise, and signal types. No load-bearing step reduces by definition to a fitted parameter, self-citation chain, or ansatz imported from the authors' prior work. The claim that the maps form a compact summary for downstream analyses follows directly from the pixel-wise amplitude recovery and is not shown to be equivalent to the input data by construction. The provided abstract and context contain no equations or citations that exhibit self-definitional or fitted-input circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The framework rests on standard domain assumptions in pulsar timing analysis and simulation-based validation; no explicit free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Pulsar timing residuals contain detectable gravitational-wave signals that can be modeled with known noise properties
    Invoked implicitly when stating compatibility with established methods and validation via simulations of noise and signals.

pith-pipeline@v0.9.0 · 5476 in / 1333 out tokens · 49726 ms · 2026-05-10T02:27:24.944862+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

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

  1. [1]

    P., Abbott, R., et al

    Abadie, J., Abbott, B. P., Abbott, R., et al. 2011, PhRvL, 107, 271102, doi: 10.1103/PhysRevLett.107.271102

  2. [2]

    2008, PhRvD, 77, 022001, doi: 10.1103/PhysRevD.77.022001

    Abbott, B., Abbott, R., Adhikari, R., et al. 2008, PhRvD, 77, 022001, doi: 10.1103/PhysRevD.77.022001

  3. [3]

    P., Abbott, R., Abbott, T

    Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016a, PhRvL, 116, 061102, doi: 10.1103/PhysRevLett.116.061102

  4. [4]

    P., Abbott, R., Abbott, T

    Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016b, PhRvL, 116, 241102, doi: 10.1103/PhysRevLett.116.241102

  5. [5]

    M., et al

    Agazie, G., Anumarlapudi, A., Archibald, A. M., et al. 2023a, ApJL, 951, L8, doi: 10.3847/2041-8213/acdac6

  6. [6]

    2023, Astrophys

    Agazie, G., Anumarlapudi, A., Archibald, A. M., et al. 2023b, ApJL, 956, L3, doi: 10.3847/2041-8213/acf4fd Ali-Ha¨ ımoud, Y., Smith, T. L., & Mingarelli, C. M. F. 2021, PhRvD, 103, 042009, doi: 10.1103/PhysRevD.103.042009

  7. [7]

    Allen, B., & Ottewill, A. C. 1997, PhRvD, 56, 545, doi: 10.1103/PhysRevD.56.545

  8. [8]

    Ashton, G., Thrane, E., & Smith, R. J. E. 2019, PhRvD, 100, 123018, doi: 10.1103/PhysRevD.100.123018

  9. [9]

    Ballmer, S. W. 2006, Classical and Quantum Gravity, 23, S179, doi: 10.1088/0264-9381/23/8/S23

  10. [10]

    2021, MNRAS, 507, 5451, doi: 10.1093/mnras/stab2479

    Banagiri, S., Criswell, A., Kuan, T., et al. 2021, MNRAS, 507, 5451, doi: 10.1093/mnras/stab2479

  11. [11]

    R., Charisi, M., et al

    Burke-Spolaor, S., Taylor, S. R., Charisi, M., et al. 2019, A&A Rv, 27, 5, doi: 10.1007/s00159-019-0115-7

  12. [12]

    Caprini, C., & Figueroa, D. G. 2018, Classical and Quantum Gravity, 35, 163001, doi: 10.1088/1361-6382/aac608

  13. [13]

    2026, PhRvD, 113, 043042, doi: 10.1103/czxp-zrd6

    Chen, Y., Zhao, S.-Y., Peng, Z.-Z., et al. 2026, PhRvD, 113, 043042, doi: 10.1103/czxp-zrd6

  14. [14]

    J., & van Haasteren, R

    Cornish, N. J., & van Haasteren, R. 2014, arXiv e-prints, arXiv:1406.4511, doi: 10.48550/arXiv.1406.4511

  15. [15]

    Vallisneri, M. 2025, arXiv e-prints, arXiv:2506.13866, doi: 10.48550/arXiv.2506.13866 18 Figure 7.Impact of the unmodelled pulsar term on the clean map S/N for two sky locations in the MPTA-like array. Top: best-covered location (bottom-right), where the coherent earth term dominates and the source is correctly recovered. Bottom: poorly-covered location (...

  16. [16]

    Monthly Notices of the Royal Astronomical Society , author =

    Cusin, G., Pitrou, C., Pijnenburg, M., & Sesana, A. 2025, MNRAS, 541, 2884, doi: 10.1093/mnras/staf1074 Di Marco, V., Zic, A., Shannon, R. M., Thrane, E., &

  17. [17]

    Kulkarni, A. D. 2025, ApJ, 990, 85, doi: 10.3847/1538-4357/adf108

  18. [18]

    A., Vallisneri, M., Taylor, S

    Ellis, J. A., Vallisneri, M., Taylor, S. R., & Baker, P. T. 2020, ENTERPRISE: Enhanced Numerical Toolbox Enabling a Robust PulsaR Inference SuitE,, Zenodo doi: 10.5281/zenodo.4059815 EPTA Collaboration, InPTA Collaboration, Antoniadis, J., et al. 2023, A&A, 678, A50, doi: 10.1051/0004-6361/202346844

  19. [19]

    2025, A&A, 703, A86, doi: 10.1051/0004-6361/202556035

    Fastidio, F., Bortolas, E., Gualandris, A., et al. 2025, A&A, 703, A86, doi: 10.1051/0004-6361/202556035

  20. [20]

    D., Taylor, S., & Mingarelli, C

    Gair, J., Romano, J. D., Taylor, S., & Mingarelli, C. M. F. 2014, PhRvD, 90, 082001, doi: 10.1103/PhysRevD.90.082001

  21. [21]

    A., Taylor, S

    Gersbach, K. A., Taylor, S. R., B´ ecsy, B., et al. 2025, arXiv e-prints, arXiv:2509.07090, doi: 10.48550/arXiv.2509.07090 G´ orski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759, doi: 10.1086/427976

  22. [22]

    J., Thrane, E., et al

    Grunthal, K., Champion, D. J., Thrane, E., et al. 2026a, arXiv e-prints, arXiv:2601.13957. https://arxiv.org/abs/2601.13957

  23. [23]

    J., & Kramer, M

    Grunthal, K., Porayko, N., Champion, D. J., & Kramer, M. 2026b, A&A, 706, A299, doi: 10.1051/0004-6361/202555394

  24. [24]

    S., Thrane, E., et al

    Grunthal, K., Nathan, R. S., Thrane, E., et al. 2025, MNRAS, 536, 1501, doi: 10.1093/mnras/stae2573

  25. [25]

    S., Simon, J., Siemens, X., & Romano, J

    Hazboun, J. S., Simon, J., Siemens, X., & Romano, J. D. 2020, ApJL, 905, L6, doi: 10.3847/2041-8213/abca92

  26. [26]

    2022, MNRAS, 509, 3488, doi:10.1093/mnras/stab3239 28 LILA Black Hole Binary Detection Landscape&Science CasesNguyenet al.2026

    Izquierdo-Villalba, D., Sesana, A., Bonoli, S., & Colpi, M. 2022, MNRAS, 509, 3488, doi: 10.1093/mnras/stab3239

  27. [27]

    2026, PhRvD, 113, 022001, doi: 10.1103/4fyl-fzt8

    Kato, R., & Takahashi, K. 2026, PhRvD, 113, 022001, doi: 10.1103/4fyl-fzt8

  28. [28]

    Z., Blecha, L., & Hernquist, L

    Kelley, L. Z., Blecha, L., & Hernquist, L. 2017, MNRAS, 464, 3131, doi: 10.1093/mnras/stw2452

  29. [29]

    2025, arXiv e-prints, arXiv:2509.07074, doi: 10.48550/arXiv.2509.07074

    Konstandin, T., Lemke, A.-M., Mitridate, A., & Perboni, E. 2025, arXiv e-prints, arXiv:2509.07074, doi: 10.48550/arXiv.2509.07074

  30. [30]

    T., Shannon, R

    Miles, M. T., Shannon, R. M., Reardon, D. J., et al. 2025a, MNRAS, 536, 1489, doi: 10.1093/mnras/stae2571

  31. [31]

    T., Shannon, R

    Miles, M. T., Shannon, R. M., Reardon, D. J., et al. 2025b, MNRAS, 536, 1467, doi: 10.1093/mnras/stae2572

  32. [32]

    Mingarelli, C. M. F., Sidery, T., Mandel, I., & Vecchio, A. 2013, PhRvD, 88, 062005, doi: 10.1103/PhysRevD.88.062005

  33. [33]

    Mingarelli, C. M. F., Casey-Clyde, J. A., Chang, Y. T., et al. 2026, arXiv e-prints, arXiv:2603.13643, doi: 10.48550/arXiv.2603.13643

  34. [34]

    2008, PhRvD, 77, 042002, doi: 10.1103/PhysRevD.77.042002

    Mitra, S., Dhurandhar, S., Souradeep, T., et al. 2008, PhRvD, 77, 042002, doi: 10.1103/PhysRevD.77.042002

  35. [35]

    2020, Phys

    Payne, E., Banagiri, S., Lasky, P. D., & Thrane, E. 2020, PhRvD, 102, 102004, doi: 10.1103/PhysRevD.102.102004

  36. [36]

    R., & Romano, J

    Pol, N., Taylor, S. R., & Romano, J. D. 2022, ApJ, 940, 173, doi: 10.3847/1538-4357/ac9836

  37. [37]

    Ravi, V., Wyithe, J. S. B., Hobbs, G., et al. 2012, ApJ, 761, 84, doi: 10.1088/0004-637X/761/2/84

  38. [38]

    J., Zic, A., Shannon, R

    Reardon, D. J., Zic, A., Shannon, R. M., et al. 2023, ApJL, 951, L6, doi: 10.3847/2041-8213/acdd02

  39. [39]

    Semenzato, F., Bellomo, N., Raccanelli, A., & Mingarelli, C. M. F. 2025, arXiv e-prints, arXiv:2510.24857, doi: 10.48550/arXiv.2510.24857

  40. [40]

    2010, PhRvD, 81, 104008, doi: 10.1103/PhysRevD.81.104008

    Sesana, A., & Vecchio, A. 2010, PhRvD, 81, 104008, doi: 10.1103/PhysRevD.81.104008

  41. [41]

    2021, Physical Review Research, 3, 043049, doi: 10.1103/PhysRevResearch.3.043049

    Talbot, C., Thrane, E., Biscoveanu, S., & Smith, R. 2021, Physical Review Research, 3, 043049, doi: 10.1103/PhysRevResearch.3.043049

  42. [42]

    Taylor, S. R. 2021, arXiv e-prints, arXiv:2105.13270, doi: 10.48550/arXiv.2105.13270

  43. [43]

    R., & Gair, J

    Taylor, S. R., & Gair, J. R. 2013, PhRvD, 88, 084001, doi: 10.1103/PhysRevD.88.084001

  44. [44]

    R., van Haasteren, R., & Sesana, A

    Taylor, S. R., van Haasteren, R., & Sesana, A. 2020, PhRvD, 102, 084039, doi: 10.1103/PhysRevD.102.084039

  45. [45]

    R., Mingarelli, C

    Taylor, S. R., Mingarelli, C. M. F., Gair, J. R., et al. 2015, PhRvL, 115, 041101, doi: 10.1103/PhysRevLett.115.041101

  46. [46]

    D., et al

    Thrane, E., Ballmer, S., Romano, J. D., et al. 2009, PhRvD, 80, 122002, doi: 10.1103/PhysRevD.80.122002

  47. [47]

    2016, MNRAS, 455, 4339, doi: 10.1093/mnras/stv2143

    Tiburzi, C., Hobbs, G., Kerr, M., et al. 2016, MNRAS, 455, 4339, doi: 10.1093/mnras/stv2143

  48. [48]

    Reaching diffraction-limited localization with coherent PTAs

    Tsai, A. C., Jow, D. L., & Pen, U.-L. 2025, arXiv e-prints, arXiv:2512.10795, doi: 10.48550/arXiv.2512.10795

  49. [49]

    2020, libstempo: Python wrapper for Tempo2,, Astrophysics Source Code Library, record ascl:2002.017 http://ascl.net/2002.017 van Haasteren, R., Allen, B., & Romano, J

    Vallisneri, M. 2020, libstempo: Python wrapper for Tempo2,, Astrophysics Source Code Library, record ascl:2002.017 http://ascl.net/2002.017 van Haasteren, R., Allen, B., & Romano, J. D. 2025, arXiv e-prints, arXiv:2509.06489, doi: 10.48550/arXiv.2509.06489

  50. [50]

    S., et al

    Wright, D., Wayt, K., Hazboun, J. S., et al. 2025, arXiv e-prints, arXiv:2512.14807, doi: 10.48550/arXiv.2512.14807

  51. [51]

    2023, Research in Astronomy and Astrophysics, 23, 075024, doi: 10.1088/1674-4527/acdfa5

    Xu, H., Chen, S., Guo, Y., et al. 2023, Research in Astronomy and Astrophysics, 23, 075024, doi: 10.1088/1674-4527/acdfa5

  52. [52]

    M., et al

    Zic, A., Hobbs, G., Shannon, R. M., et al. 2022, MNRAS, 516, 410, doi: 10.1093/mnras/stac2100