pith. machine review for the scientific record.
sign in

arxiv: 2510.04537 · v2 · submitted 2025-10-06 · 🌌 astro-ph.HE · gr-qc

Finding Supermassive Black Hole Binary Mergers in Pulsar Timing Array Data

Pith reviewed 2026-05-18 09:52 UTC · model grok-4.3

classification 🌌 astro-ph.HE gr-qc
keywords supermassive black hole binariespulsar timing arraysgravitational wave memorywaveform modelingparameter estimationmerger signalscontinuous gravitational waves
0
0 comments X

The pith

A complete waveform model covering the full merger process lets pulsar timing arrays recover parameters of supermassive black hole binary mergers from simulated data.

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

The paper develops a search method for rare supermassive black hole binary mergers in pulsar timing array observations by using one waveform that follows the binary from its slow inspiral all the way through merger, ringdown, and the permanent gravitational-wave memory shift. This single model treats both the oscillating continuous signal and the non-oscillatory memory effect together instead of handling them separately. Tests on simulated pulsar timing datasets show that strong enough signals can be recovered with log Bayes factors above 10, yielding constraints on chirp mass and luminosity distance while also producing sky localizations of a few degrees. The same tests reveal that the simpler memory-burst approximations still used in the field return biased amplitude and parameter values even when tuned to match the full signal.

Core claim

A physically complete waveform model including inspiral, merger, ringdown, and gravitational-wave memory enables a unified treatment of continuous emission and memory signals in PTA data, allowing recovery of SMBHB parameters with log Bayes factors greater than 10 on simulated datasets and exposing biases in common memory-burst approximations.

What carries the argument

The physically complete SMBHB waveform model that includes inspiral, merger, ringdown, and gravitational-wave memory, used to perform unified parameter estimation on PTA timing residuals.

If this is right

  • Strong signals yield constraints on chirp mass and luminosity distance despite their degeneracy.
  • Sky localization to a few degrees opens the possibility of electromagnetic follow-up observations.
  • Common memory-burst approximations produce biased strain amplitudes and source parameters relative to the full model.
  • A pathway exists for systematic searches of SMBHB mergers in PTA data that uses complete waveform models rather than piecemeal approximations.

Where Pith is reading between the lines

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

  • Real PTA searches using this model could be combined with electromagnetic surveys to confirm mergers through multi-messenger detections.
  • The bias findings suggest that existing upper limits on memory signals may need revision once full waveforms are applied to archival data.
  • Extending the approach to eccentric or precessing binaries would test how robust the recovery remains for more general systems.

Load-bearing premise

The simulated PTA datasets accurately capture the noise properties, timing residuals, and array configuration of real observations.

What would settle it

Applying both the full waveform model and a memory-burst approximation to the same real PTA dataset that contains a candidate merger signal and finding that the two methods return statistically inconsistent values for chirp mass or strain amplitude.

Figures

Figures reproduced from arXiv: 2510.04537 by Boris Goncharov, Rutger van Haasteren, Sharon Mary Tomson.

Figure 1
Figure 1. Figure 1: FIG. 1: Gravitational wave strain and timing residuals from a merger of a non-spinning supermassive black hole [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: The figure illustrates the posterior distributions from our simulation recovery studies with 25 pulsars. Simu [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Sky maps of the posterior distributions of simulated supermassive black hole binary mergers, shown in [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Timing residuals from a pulsar located at (ra,dec) = (258 [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Lower limits on luminosity distance as a function [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Posterior distributions of SMBHB merger with null memory parameters from simulation studies on 25 pulsars. [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Posterior distribution of the SMBHB merger from simulation studies on 25 pulsars. The dataset used is [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Posterior distributions of SMBHB merger with null memory parameters from simulation studies on 25 pulsars. [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Posterior distribution of the SMBHB merger from simulation studies on 25 pulsars. Parameters of the binary [PITH_FULL_IMAGE:figures/full_fig_p019_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Posterior distributions of SMBHB merger with null memory parameters from simulation studies on 25 [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Posterior distribution of the SMBHB merger from simulation studies on 25 pulsars. Parameters of the binary [PITH_FULL_IMAGE:figures/full_fig_p021_12.png] view at source ↗
read the original abstract

Galaxy observations suggest that mergers of supermassive black hole binaries (SMBHBs) are rare events, with rates of order one per decade across the observable Universe. We present a framework to search for merging SMBHBs in pulsar timing array (PTA) data using a physically complete waveform model including inspiral, merger, ringdown, and gravitational-wave memory. This enables a unified treatment of continuous emission and the non-oscillatory memory signal. Using simulated PTA datasets, we demonstrate parameter estimation for representative systems with chirp masses of $10^8$ and $10^{10}~M_\odot$ at distances of $3$ Mpc to $100$ Mpc respectively. For sufficiently strong signals, we recover binaries with log Bayes factors >10 and constrain chirp mass and luminosity distance, subject to their characteristic degeneracy. Sky localization uncertainties of a few degrees could potentially enable electromagnetic follow-up and multi-messenger observations of SMBHB mergers. We further demonstrate that commonly used memory burst approximations lead to biased strain amplitudes and inferred source parameters when compared to the full SMBHB waveform, even when optimally tuned. These results establish a pathway for searching for SMBHB mergers with PTAs using complete waveform models.

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

1 major / 2 minor

Summary. The manuscript presents a framework for searching for supermassive black hole binary (SMBHB) mergers in pulsar timing array (PTA) data using a physically complete waveform model that includes inspiral, merger, ringdown, and gravitational-wave memory phases. This enables unified treatment of continuous emission and memory signals. On simulated PTA datasets for representative systems (chirp masses 10^8 and 10^10 M_⊙ at distances 3–100 Mpc), the authors demonstrate Bayesian parameter recovery with log Bayes factors >10 for strong signals, constraints on chirp mass and luminosity distance (subject to degeneracy), sky localization to a few degrees, and biases in common memory-burst approximations even when optimally tuned. The work positions this as a pathway for future real-data searches and multi-messenger follow-up.

Significance. If the simulation results hold under realistic conditions, the unified waveform approach would be a valuable advance for PTA analyses of rare SMBHB mergers, improving upon separate treatments of continuous waves and memory bursts. The explicit bias comparison to approximations is a concrete strength, highlighting the need for complete models. The simulation-based demonstration of parameter recovery and Bayes-factor thresholds provides a reproducible starting point, though the absence of real-data application or full error-budget analysis limits immediate observational impact.

major comments (1)
  1. The central results on log Bayes factors >10 and bias in memory approximations rest on the fidelity of the simulated PTA noise model (red noise, DM variations, pulsar terms, array configuration). The manuscript does not provide a quantitative comparison of these simulated residuals to actual PTA datasets (e.g., NANOGrav 15-year or EPTA DR2), which is load-bearing for translating the recovery performance and degeneracy structure to real observations.
minor comments (2)
  1. The abstract states 'for sufficiently strong signals' without a specific SNR or strain threshold; this should be quantified in the results section with reference to the simulated injections.
  2. Clarify in the methods whether the waveform model is implemented in the time domain or frequency domain for PTA timing residuals, and specify the sampling of the memory step function.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive and detailed review of our manuscript. We address the major comment below and have made revisions to strengthen the connection between our simulations and real PTA observations.

read point-by-point responses
  1. Referee: The central results on log Bayes factors >10 and bias in memory approximations rest on the fidelity of the simulated PTA noise model (red noise, DM variations, pulsar terms, array configuration). The manuscript does not provide a quantitative comparison of these simulated residuals to actual PTA datasets (e.g., NANOGrav 15-year or EPTA DR2), which is load-bearing for translating the recovery performance and degeneracy structure to real observations.

    Authors: We thank the referee for this important observation. Our simulated datasets employ red-noise spectra, DM variations, and pulsar-term contributions drawn from standard models in the PTA literature that are calibrated to reproduce the statistical properties of current arrays. We agree that an explicit quantitative comparison would improve the manuscript's utility for interpreting real-data searches. In the revised version we will add a dedicated subsection (and accompanying figure) in the Methods that directly compares the power spectral densities, residual RMS values, and autocorrelation properties of our simulated noise realizations against published noise parameters from the NANOGrav 15-year and EPTA DR2 datasets. This addition will clarify the fidelity of the simulations while preserving the paper's focus on the complete waveform model. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results from forward simulation and Bayesian recovery on independent mock data.

full rationale

The paper's central claims rest on generating simulated PTA timing residuals for representative SMBHB systems and then performing Bayesian parameter estimation and model comparison using a complete waveform model. These steps are forward simulations followed by recovery; they do not involve fitting parameters to a subset of data and then relabeling the fit as a prediction, nor do they reduce any derived quantity to a self-definition or self-citation chain. The comparison to memory-burst approximations is an internal consistency check on the same simulated datasets and does not create a tautology. No load-bearing uniqueness theorems, ansatzes smuggled via prior self-citations, or renamings of known results appear in the provided abstract or described methodology. The derivation chain is therefore self-contained against external benchmarks (simulated data with stated noise properties).

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The method relies on standard general-relativity waveform models and PTA noise assumptions that are not derived in the paper; no new free parameters are introduced beyond the usual source parameters (chirp mass, distance, sky location) that are fitted to the simulated data.

free parameters (1)
  • chirp mass and luminosity distance
    Fitted parameters recovered from the simulated signals; their degeneracy is noted but not removed by additional priors in the abstract.
axioms (2)
  • domain assumption Standard general-relativity inspiral-merger-ringdown-memory waveform is an accurate description of SMBHB signals in the PTA band.
    Invoked when constructing the complete waveform model used for injection and recovery.
  • domain assumption PTA noise is stationary and can be modeled sufficiently well in simulations to test detection and parameter estimation.
    Required for the simulated datasets to be representative.

pith-pipeline@v0.9.0 · 5755 in / 1488 out tokens · 28327 ms · 2026-05-18T09:52:03.933427+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

76 extracted references · 76 canonical work pages · 33 internal anchors

  1. [1]

    It also cor- responds to the frequency of gravitational wave signals

    Pulsar noise It is convenient to express pulsar noise in terms of the power spectral density of residuals,P(f), wherefis the Fourier frequency of the timing residuals. It also cor- responds to the frequency of gravitational wave signals. White noise (WN) corresponds to the constant power spectral density across frequency, and it is assumed to be known in ...

  2. [2]

    Namely, the stochastic gravitational wave background (GWB)

    Gravitational Wave Background In some simulations, we test the robustness of our parameter estimation for SMBHB mergers to the pres- ence of an additional gravitational wave signal. Namely, the stochastic gravitational wave background (GWB). It is modeled as a Gaussian process with the same power spectrum in all pulsars and additional inter-pulsar corre- ...

  3. [3]

    We construct a multivari- ate Gaussian likelihood

    Likelihood In this work, we perform an analysis of simulated PTA data using Bayesian inference. We construct a multivari- ate Gaussian likelihood. The likelihood of the vector of timing residualsδtis a multivariate Gaussian function [50, 58, 59] L(δt|θ) = 1p det(2πC) exp h − 1 2(δt−µ) TC −1(δt−µ) i , (2.13) whereµrepresents the time series model for an SM...

  4. [4]

    Simulationwn+rn+mem: white noise, red noise whose hyperparameters (ARN, γRN) are drawn from the uniform distributions [−15,−14] and [1,4], re- spectively, , and SMBHB merger

  5. [6]

    This differs from the smoother, continuous buildup of memory in circular binaries

    presents a derivation of the null memory from bina- ries on elliptical orbits and shows that the memory flux becomes sharply peaked near each pericenter passage as eccentricity increases. This differs from the smoother, continuous buildup of memory in circular binaries. While the total memory evolves toward the circular case as the orbit circularizes, sys...

  6. [7]

    R. S. Foster and D. C. Backer, Astrophys. J.361, 300 (1990)

  7. [8]

    M. A. McLaughlin, Classical and Quantum Gravity30, 224008 (2013), arXiv:1310.0758 [astro-ph.IM]

  8. [9]

    R. N. Manchester, G. Hobbs, M. Bailes, W. A. Coles, W. van Straten, M. J. Keith, R. M. Shannon, N. D. R. Bhat, A. Brown, S. G. Burke-Spolaor, D. J. Cham- pion, A. Chaudhary, R. T. Edwards, G. Hampson, A. W. Hotan, A. Jameson, F. A. Jenet, M. J. Kesteven, J. Khoo, J. Kocz, K. Maciesiak, S. Oslowski, V. Ravi, J. R. Reynolds, J. M. Sarkissian, J. P. W. Verbi...

  9. [10]

    High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array

    G. Desvignes, R. N. Caballero, L. Lentati, J. P. W. Ver- biest, D. J. Champion, B. W. Stappers, G. H. Janssen, P. Lazarus, S. Os lowski, S. Babak, C. G. Bassa, P. Brem, M. Burgay, I. Cognard, J. R. Gair, E. Graikou, L. Guille- mot, J. W. T. Hessels, A. Jessner, C. Jordan, R. Karup- pusamy, M. Kramer, A. Lassus, K. Lazaridis, K. J. Lee, K. Liu, A. G. Lyne,...

  10. [11]

    Tarafdar, K

    P. Tarafdar, K. Nobleson, P. Rana, J. Singha, M. A. Kr- ishnakumar, B. C. Joshi, A. K. Paladi, N. Kolhe, N. D. Batra, N. Agarwal, A. Bathula, S. Dandapat, S. Desai, L. Dey, S. Hisano, P. Ingale, R. Kato, D. Kharbanda, T. Kikunaga, P. Marmat, B. A. Pandian, T. Prabu, A. Srivastava, M. Surnis, S. C. Susarla, A. Susobhanan, K. Takahashi, P. Arumugam, M. Bagc...

  11. [12]

    K. J. Lee, inFrontiers in Radio Astronomy and F AST Early Sciences Symposium 2015, Astronomical Society of the Pacific Conference Series, Vol. 502, edited by L. Qain and D. Li (2016) p. 19

  12. [13]

    M. T. Miles, R. M. Shannon, M. Bailes, D. J. Rear- don, M. J. Keith, A. D. Cameron, A. Parthasarathy, M. Shamohammadi, R. Spiewak, W. van Straten, S. Buchner, F. Camilo, M. Geyer, A. Karastergiou, M. Kramer, M. Serylak, G. Theureau, and V. Venkatra- man Krishnan, Monthly Notices of the Royal Astronom- ical Society519, 3976 (2023), arXiv:2212.04648 [astro-...

  13. [14]

    J. P. W. Verbiest, L. Lentati, G. Hobbs, R. van Haasteren, P. B. Demorest, G. H. Janssen, J. B. Wang, G. Desvignes, R. N. Caballero, M. J. Keith, D. J. Cham- pion, Z. Arzoumanian, S. Babak, C. G. Bassa, N. D. R. Bhat, A. Brazier, P. Brem, M. Burgay, S. Burke-Spolaor, S. J. Chamberlin, S. Chatterjee, B. Christy, I. Cognard, J. M. Cordes, S. Dai, T. Dolch, ...

  14. [15]

    Arzoumanian, P

    Z. Arzoumanian, P. T. Baker, H. Blumer, B. B´ ecsy, A. Brazier, P. R. Brook, S. Burke-Spolaor, S. Chatterjee, S. Chen, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromartie, M. E. Decesar, P. B. Demorest, T. Dolch, J. A. Ellis, E. C. Ferrara, W. t. Fiore, and Nanograv Collaboration, Astrophysical Journal Letters905, L34 (2020), arXiv:2009.04496 [astro-ph.HE]

  15. [16]

    Goncharov, R

    B. Goncharov, R. M. Shannon, D. J. Reardon, G. Hobbs, A. Zic, M. Bailes, M. Cury lo, S. Dai, M. Kerr, M. E. Lower, R. N. Manchester, R. Mandow, H. Middleton, M. T. Miles, A. Parthasarathy, E. Thrane, N. Thyagara- jan, X. Xue, X. J. Zhu, A. D. Cameron, Y. Feng, R. Luo, C. J. Russell, J. Sarkissian, R. Spiewak, S. Wang, J. B. Wang, L. Zhang, and S. Zhang, A...

  16. [17]

    Goncharov, E

    B. Goncharov, E. Thrane, R. M. Shannon, J. Harms, N. D. R. Bhat, G. Hobbs, M. Kerr, R. N. Manch- ester, D. J. Reardon, C. J. Russell, X.-J. Zhu, and A. Zic, Astrophysical Journal Letters932, L22 (2022), arXiv:2206.03766 [gr-qc]

  17. [18]

    S. Chen, R. N. Caballero, Y. J. Guo, A. Chalumeau, K. Liu, G. Shaifullah, K. J. Lee, S. Babak, G. Desvignes, A. Parthasarathy, H. Hu, E. van der Wateren, J. Anto- niadis, A. S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Burgay, D. J. Champion, I. Cognard, M. Falxa, R. D. Ferdman, P. C. C. Freire, J. R. Gair, E. Graikou, L. Guillemot, J. Jang, G. H. Janss...

  18. [20]

    The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

    G. Agazie, A. Anumarlapudi, A. M. Archibald, Z. Ar- zoumanian, P. T. Baker, B. B´ ecsy, L. Blecha, A. Brazier, P. R. Brook, S. Burke-Spolaor, R. Burnette, R. Case, M. Charisi, S. Chatterjee, K. Chatziioannou, B. D. Cheeseboro, S. Chen, T. Cohen, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromartie, K. Crowter, C. J. Cutler, M. E. Decesar, D. Degan, P...

  19. [21]

    The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals

    EPTA Collaboration, InPTA Collaboration, J. An- toniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Cham- pion, S. Chanlaridis, S. Chen, I. Cognard, S. Danda- pat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-B...

  20. [22]

    D. J. Reardon, A. Zic, R. M. Shannon, G. B. Hobbs, M. Bailes, V. Di Marco, A. Kapur, A. F. Rogers, E. Thrane, J. Askew, N. D. R. Bhat, A. Cameron, M. Cury lo, W. A. Coles, S. Dai, B. Goncharov, M. Kerr, A. Kulkarni, Y. Levin, M. E. Lower, R. N. Manchester, R. Mandow, M. T. Miles, R. S. Nathan, S. Os lowski, C. J. Russell, R. Spiewak, S. Zhang, and X.-J. Z...

  21. [23]

    M. T. Miles, R. M. Shannon, D. J. Reardon, M. Bailes, D. J. Champion, M. Geyer, P. Gitika, K. Grunthal, M. J. Keith, M. Kramer, A. D. Kulkarni, R. S. Nathan, A. Parthasarathy, J. Singha, G. Theureau, E. Thrane, F. Abbate, S. Buchner, A. D. Cameron, F. Camilo, B. E. Moreschi, G. Shaifullah, M. Shamohammadi, A. Possenti, and V. V. Krishnan, Monthly Notices ...

  22. [24]

    H. Xu, S. Chen, Y. Guo, J. Jiang, B. Wang, J. Xu, Z. Xue, R. N. Caballero, J. Yuan, Y. Xu, J. Wang, L. Hao, J. Luo, K. Lee, J. Han, P. Jiang, Z. Shen, M. Wang, N. Wang, R. Xu, X. Wu, R. Manchester, L. Qian, X. Guan, M. Huang, C. Sun, and Y. Zhu, Re- search in Astronomy and Astrophysics23, 075024 (2023), arXiv:2306.16216 [astro-ph.HE]

  23. [25]

    The NANOGrav 15-year Data Set: Search for Signals from New Physics

    A. Afzal, G. Agazie, A. Anumarlapudi, A. M. Archibald, Z. Arzoumanian, P. T. Baker, B. B´ ecsy, J. J. Blanco- Pillado, L. Blecha, K. K. Boddy, A. Brazier, P. R. Brook, S. Burke-Spolaor, R. Burnette, R. Case, M. Charisi, S. Chatterjee, K. Chatziioannou, B. D. Cheeseboro, S. Chen, T. Cohen, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromartie, K. Crowt...

  24. [26]

    An- toniadis, P

    EPTA Collaboration, InPTA Collaboration, J. An- toniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Cham- pion, S. Chanlaridis, S. Chen, I. Cognard, S. Danda- pat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-B...

  25. [27]

    Ellis, F

    J. Ellis, F. Jenet, and X. Siemens, inAmerican Astro- nomical Society Meeting Abstracts #219, American As- 14 tronomical Society Meeting Abstracts, Vol. 219 (2012) p. 146.20

  26. [28]

    Gravitational wave bursts from cosmic strings

    T. Damour and A. Vilenkin, Phys. Rev. Lett.85, 3761 (2000), arXiv:gr-qc/0004075 [gr-qc]

  27. [29]

    Y. Xia, J. Wang, S. Kuroyanagi, W. Yan, Y. Wen, A. Ka- pur, J. Zou, Y. Feng, V. Di Marco, S. Mishra, C. J. Russell, S. Wang, D. Zhao, and X. Zhu, Universe11, 81 (2025), arXiv:2502.21069 [gr-qc]

  28. [30]

    Turner, Astrophys

    M. Turner, Astrophys. J.216, 610 (1977)

  29. [31]

    The gravitational-wave memory effect

    M. Favata, Classical and Quantum Gravity27, 084036 (2010), arXiv:1003.3486 [gr-qc]

  30. [32]

    Gravitational-wave memory and pulsar timing arrays

    R. van Haasteren and Y. Levin, Monthly Notices of the Royal Astronomical Society401, 2372 (2010), arXiv:0909.0954 [astro-ph.IM]

  31. [33]

    Dandapat, A

    S. Dandapat, A. Susobhanan, L. Dey, A. Gopakumar, P. T. Baker, and P. Jetzer, Phys. Rev. D109, 103018 (2024), arXiv:2402.03472 [astro-ph.HE]

  32. [34]

    The NANOGrav 15 yr Data Set: Targeted Searches for Supermassive Black Hole Binaries

    N. Agarwal, G. Agazie, A. Anumarlapudi, A. M. Archibald, Z. Arzoumanian, J. G. Baier, P. T. Baker, B. Becsy, L. Blecha, A. Brazier, P. R. Brook, S. Burke- Spolaor, R. Burnette, R. Case, J. A. Casey-Clyde, Y.-T. Chang, M. Charisi, S. Chatterjee, T. Cohen, P. Coppi, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromar- tie, K. Crowter, M. E. DeCesar, P. B...

  33. [35]

    J. A. Taylor, R. Burnette, B. B´ ecsy, and N. J. Cornish, Phys. Rev. D111, 022006 (2025), arXiv:2408.07864 [gr- qc]

  34. [36]

    Dandapat, M

    S. Dandapat, M. Ebersold, A. Susobhanan, P. Rana, A. Gopakumar, S. Tiwari, M. Haney, H. M. Lee, and N. Kolhe, Phys. Rev. D108, 024013 (2023), arXiv:2305.19318 [gr-qc]

  35. [37]

    Epstein, Astrophys

    R. Epstein, Astrophys. J.223, 1037 (1978)

  36. [38]

    D. R. Madison, D. F. Chernoff, and J. M. Cordes, Phys. Rev. D96, 123016 (2017), arXiv:1710.04974 [astro- ph.GA]

  37. [39]

    N. Sago, K. Ioka, T. Nakamura, and R. Yamazaki, Phys. Rev. D70, 104012 (2004), arXiv:gr-qc/0405067 [gr-qc]

  38. [40]

    Y. B. Zel’dovich and A. G. Polnarev, Soviet Astronomy 18, 17 (1974)

  39. [41]

    V. B. Braginsky and K. S. Thorne, Nature (London)327, 123 (1987)

  40. [42]

    Christodoulou, Phys

    D. Christodoulou, Phys. Rev. Lett.67, 1486 (1991)

  41. [43]

    Nonlinear gravitational-wave memory from binary black hole mergers

    M. Favata, Astrophysical Journal Letters696, L159 (2009), arXiv:0902.3660 [astro-ph.SR]

  42. [44]

    Search for Memory and Inspiral Gravitational Waves from Super-Massive Binary Black Holes with Pulsar Timing Arrays

    N. Seto, Monthly Notices of the Royal Astronomical So- ciety400, L38 (2009), arXiv:0909.1379 [astro-ph.CO]

  43. [45]

    M. S. Pshirkov, D. Baskaran, and K. A. Postnov, Monthly Notices of the Royal Astronomical Society402, 417 (2010), arXiv:0909.0742 [astro-ph.CO]

  44. [46]

    J. B. Wang, G. Hobbs, W. Coles, R. M. Shannon, X. J. Zhu, D. R. Madison, M. Kerr, V. Ravi, M. J. Keith, R. N. Manchester, Y. Levin, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, S. Dai, S. Os lowski, W. van Straten, L. Toomey, N. Wang, and L. Wen, Monthly Notices of the Royal Astronomical Society446, 1657 (2015), arXiv:1410.3323 [astro-ph.GA]

  45. [47]

    Aggarwal, Z

    K. Aggarwal, Z. Arzoumanian, P. T. Baker, A. Bra- zier, P. R. Brook, S. Burke-Spolaor, S. Chatterjee, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromartie, K. Crowter, M. DeCesar, P. B. Demorest, T. Dolch, J. A. Ellis, R. D. Ferdman, E. C. Ferrara, E. Fonseca, N. Garver-Daniels, P. Gentile, and D. t. Good, The As- trophysical Journal889, 38 (2020)

  46. [48]

    Agazie, Z

    G. Agazie, Z. Arzoumanian, P. T. Baker, B. B´ ecsy, L. Blecha, H. Blumer, A. Brazier, P. R. Brook, S. Burke- Spolaor, R. Burnette, R. Case, J. A. Casey-Clyde, M. Charisi, S. Chatterjee, T. Cohen, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cromartie, M. E. Decesar, D. Degan, P. B. Demorest, T. Dolch, B. Drachler, J. A. Ellis, R. D. Ferdman, E. C. Ferr...

  47. [49]

    Agazie, A

    G. Agazie, A. Anumarlapudi, A. M. Archibald, Z. Arzou- manian, J. G. Baier, P. T. Baker, B. B´ ecsy, L. Blecha, A. Brazier, P. R. Brook, S. Burke-Spolaor, R. Burnette, J. A. Casey-Clyde, M. Charisi, S. Chatterjee, T. Cohen, J. M. Cordes, N. J. Cornish, F. Crawford, H. T. Cro- martie, K. Crowter, M. E. DeCesar, P. B. Demorest, H. Deng, L. Dey, T. Dolch, E....

  48. [50]

    J. Yoo, K. Mitman, V. Varma, M. Boyle, S. E. Field, N. Deppe, F. H´ ebert, L. E. Kidder, J. Moxon, H. P. Pfeiffer, M. A. Scheel, L. C. Stein, S. A. Teukolsky, W. Throwe, and N. L. Vu, Phys. Rev. D108, 064027 (2023), arXiv:2306.03148 [gr-qc]

  49. [51]

    Cauchy-characteristic extraction in numerical relativity

    N. Bishop, R. Gomez, L. Lehner, and J. Winicour, arXiv e-prints , gr-qc/9705033 (1997), arXiv:gr-qc/9705033 [gr- qc]

  50. [52]

    Moxon, M

    J. Moxon, M. A. Scheel, S. A. Teukolsky, N. Deppe, N. Vu, F. H´ ebert, L. E. Kidder, and W. Throwe, Phys. Rev. D107, 064013 (2023)

  51. [53]

    Mitman, J

    K. Mitman, J. Moxon, M. A. Scheel, S. A. Teukolsky, M. Boyle, N. Deppe, L. E. Kidder, and W. Throwe, Phys. Rev. D102, 104007 (2020), arXiv:2007.11562 [gr-qc]

  52. [54]

    J. W. McKee, G. H. Janssen, B. W. Stappers, A. G. Lyne, R. N. Caballero, L. Lentati, G. Desvignes, A. Jessner, C. A. Jordan, R. Karuppusamy, M. Kramer, I. Cognard, D. J. Champion, E. Graikou, P. Lazarus, S. Os lowski, 15 D. Perrodin, G. Shaifullah, C. Tiburzi, and J. P. W. Ver- biest, Monthly Notices of the Royal Astronomical Society 461, 2809 (2016), arX...

  53. [55]

    M. T. Miles, R. M. Shannon, D. J. Reardon, M. Bailes, D. J. Champion, M. Geyer, P. Gitika, K. Grunthal, M. J. Keith, M. Kramer, A. D. Kulkarni, R. S. Nathan, A. Parthasarathy, N. K. Porayko, J. Singha, G. Theureau, F. Abbate, S. Buchner, A. D. Cameron, F. Camilo, B. E. Moreschi, G. Shaifullah, M. Shamo- hammadi, and V. V. Krishnan, Monthly Notices of the ...

  54. [56]

    On measuring the gravitational-wave background using Pulsar Timing Arrays

    R. van Haasteren, Y. Levin, P. McDonald, and T. Lu, Monthly Notices of the Royal Astronomical Society395, 1005 (2009), arXiv:0809.0791 [astro-ph]

  55. [57]

    Pulsar timing analysis in the presence of correlated noise

    W. Coles, G. Hobbs, D. J. Champion, R. N. Manchester, and J. P. W. Verbiest, Monthly Notices of the Royal Astronomical Society418, 561 (2011), arXiv:1107.5366 [astro-ph.IM]

  56. [58]

    R. M. Shannon and J. M. Cordes, Astrophys. J.725, 1607 (2010), arXiv:1010.4794 [astro-ph.SR]

  57. [59]

    Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data

    L. Lentati, P. Alexander, M. P. Hobson, S. Taylor, J. Gair, S. T. Balan, and R. van Haasteren, Phys. Rev. D87, 104021 (2013), arXiv:1210.3578 [astro-ph.IM]

  58. [60]

    Low-rank approximations for large stationary covariance matrices, as used in the Bayesian and generalized-least-squares analysis of pulsar-timing data

    R. van Haasteren and M. Vallisneri, Monthly Notices of the Royal Astronomical Society446, 1170 (2015), arXiv:1407.6710 [astro-ph.IM]

  59. [61]

    Allen and J

    B. Allen and J. D. Romano, Phys. Rev. Lett.134, 031401 (2025), arXiv:2407.10968 [gr-qc]

  60. [62]

    Crisostomi, R

    M. Crisostomi, R. van Haasteren, P. M. Meyers, and M. Vallisneri, arXiv e-prints , arXiv:2506.13866 (2025), arXiv:2506.13866 [astro-ph.IM]

  61. [63]

    R. W. Hellings and G. S. Downs, Astrophysical Journal Letters265, L39 (1983)

  62. [64]

    The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

    Z. Arzoumanian, A. Brazier, S. Burke-Spolaor, S. J. Chamberlin, S. Chatterjee, B. Christy, J. M. Cordes, N. J. Cornish, K. Crowter, P. B. Demorest, X. Deng, T. Dolch, J. A. Ellis, R. D. Ferdman, E. Fonseca, N. Garver-Daniels, M. E. Gonzalez, F. Jenet, G. Jones, M. L. Jones, V. M. Kaspi, M. Koop, M. T. Lam, T. J. W. Lazio, L. Levin, A. N. Lommen, D. R. Lor...

  63. [65]

    The International Pulsar Timing Array second data release: Search for an isotropic gravitational wave background,

    J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A. S. Bak Nielsen, P. T. Baker, C. G. Bassa, B. B´ ecsy, A. Berthereau, M. Bonetti, A. Brazier, P. R. Brook, M. Burgay, S. Burke-Spolaor, R. N. Caballero, J. A. Casey-Clyde, and C.et al., Monthly Notices of the Royal Astronomical Society510, 4873 (2022), arXiv:2201.03980 [astro-ph.HE]

  64. [66]

    Resolving multiple supermassive black hole binaries with pulsar timing arrays II: genetic algorithm implementation

    A. Petiteau, S. Babak, A. Sesana, and M. de Ara´ ujo, Phys. Rev. D87, 064036 (2013), arXiv:1210.2396 [astro- ph.CO]

  65. [67]

    Ellis and R

    J. Ellis and R. van Haasteren, PTMCMCSampler: Paral- lel tempering MCMC sampler package written in Python, Astrophysics Source Code Library, record ascl:1912.017 (2019)

  66. [68]

    The stochastic gravitational-wave background from massive black hole binary systems: implications for observations with Pulsar Timing Arrays

    A. Sesana, A. Vecchio, and C. N. Colacino, Monthly No- tices of the Royal Astronomical Society390, 192 (2008), arXiv:0804.4476 [astro-ph]

  67. [69]

    Gravitational memory in binary black hole mergers

    D. Pollney and C. Reisswig, Astrophysical Journal Let- ters732, L13 (2011), arXiv:1004.4209 [gr-qc]

  68. [70]

    The gravitational-wave memory from eccentric binaries

    M. Favata, Phys. Rev. D84, 124013 (2011), arXiv:1108.3121 [gr-qc]

  69. [71]
  70. [72]

    E. M. Guti´ errez, L. Combi, S. C. Noble, M. Campan- elli, J. H. Krolik, F. L´ opez Armengol, and F. Garc´ ıa, Astrophys. J.928, 137 (2022), arXiv:2112.09773 [astro- ph.HE]

  71. [73]

    Bogdanovi´ c, M

    T. Bogdanovi´ c, M. C. Miller, and L. Blecha, Living Re- views in Relativity25, 3 (2022), arXiv:2109.03262 [astro- ph.HE]

  72. [74]

    Cardinal Tremblay, B

    J. Cardinal Tremblay, B. Goncharov, R. van Haasteren, N. D. R. Bhat, Z.-C. Chen, V. Di Marco, S. Iguchi, A. Ka- pur, W. Ling, R. Mandow, S. Mishra, D. J. Reardon, R. M. Shannon, H. Sudou, J. Wang, S.-Y. Zhao, X.-J. Zhu, and A. Zic, arXiv e-prints , arXiv:2508.20007 (2025), arXiv:2508.20007 [astro-ph.HE]

  73. [75]

    B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Aber- nathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai, M. C. Araya, C. C. Arceneaux, J. S. Areeda, N....

  74. [76]

    Field, V

    S. Field, V. Varma, J. Blackman, B. Gadre, C. Gal- ley, T. Islam, K. Mitman, M. P¨ urrer, A. Ravichan- dran, M. Scheel, L. Stein, and J. Yoo, The Journal of Open Source Software10, 7073 (2025), arXiv:2504.08839 [astro-ph.IM]

  75. [77]

    J. A. Ellis, M. Vallisneri, S. R. Taylor, and P. T. Baker, ENTERPRISE: Enhanced Numerical Toolbox Enabling a Robust PulsaR Inference SuitE (2020). 16 Appendix A: Appendix

  76. [78]

    These results provide a complete view of the parameter estimation performance of our SMBHB merger model across different simulated scenarios

    Posterior distributions The appendix presents the full posterior distributions for all simulation runs discussed in the main text, including the recovered sky locations, masses, distances, and other model parameters. These results provide a complete view of the parameter estimation performance of our SMBHB merger model across different simulated scenarios...