pith. sign in

arxiv: 2512.08666 · v3 · submitted 2025-12-09 · 🌌 astro-ph.CO · gr-qc

Stochastic gravitational-wave background search using data from five pulsar timing arrays

Wang-Wei Yu , Bruce Allen This is my paper

Pith reviewed 2026-05-17 00:02 UTC · model grok-4.3

classification 🌌 astro-ph.CO gr-qc
keywords stochastic gravitational wave backgroundpulsar timing arraysHellings-Downs correlationnanohertz gravitational wavesBayesian posterior analysisdirect combination methodmulti-PTA data merging
0
0 comments X

The pith

Combined data from five pulsar timing arrays shows a stochastic gravitational-wave background signal consistent with predictions but below the 5-sigma detection threshold.

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

The paper combines public pulse time-of-arrival data from five pulsar timing arrays into a dataset of 121 pulsars to search for a stationary, isotropic, unpolarized stochastic gravitational-wave background at nanohertz frequencies. It employs a direct combination method to merge data for pulsars observed by multiple arrays, avoiding reconciliation of different timing models. The analysis yields posterior distributions for the amplitude and spectral index of a power-law SGWB energy-density spectrum that are consistent with a nonzero signal. However, the statistical significance, assessed through Bayesian odds ratios and noise-marginalized false-alarm probabilities, remains below the conventional 5 sigma threshold for detection. The reconstructed angular correlation of timing residuals between pulsar pairs matches the Hellings-Downs prediction expected for an isotropic SGWB.

Core claim

Using a combined 121-pulsar dataset from five PTAs that is approximately four times larger than any individual array, the analysis finds posterior distributions for the SGWB power-law amplitude A_gw and exponent gamma_gw that are consistent with a nonzero amplitude, with the inter-pulsar timing-residual correlations matching the Hellings-Downs curve, though the detection significance falls short of 5 sigma.

What carries the argument

The direct combination method for merging astrophysical timing models and data from multiple PTAs for shared pulsars, together with the power-law model for the SGWB energy-density spectrum.

If this is right

  • The larger combined dataset increases the sensitivity of the SGWB search compared to individual PTA analyses.
  • The results support the presence of a stochastic background whose angular correlation follows the Hellings-Downs pattern.
  • Bayesian odds ratios and p-values for detection statistics indicate the signal has not yet reached conventional discovery thresholds.
  • Future improvements in PTA data volume and precision could push the significance higher if the signal is real.

Where Pith is reading between the lines

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

  • Extending the direct combination approach to future PTAs with overlapping observations could accelerate detection timelines.
  • The consistency with Hellings-Downs strengthens the case for an astrophysical origin from supermassive black hole binaries.
  • Similar combination techniques might be applied to other multi-experiment gravitational wave searches to boost sensitivity.
  • Additional years of timing data could provide the observations needed to cross the 5 sigma threshold.

Load-bearing premise

The direct combination method for merging data from pulsars observed by multiple PTAs correctly preserves the astrophysical timing residuals and noise properties without introducing unaccounted biases or correlations.

What would settle it

A significant deviation of the reconstructed inter-pulsar timing-residual correlation function from the Hellings-Downs prediction, or a posterior distribution for A_gw that becomes consistent with zero upon inclusion of more data or refined noise modeling.

Figures

Figures reproduced from arXiv: 2512.08666 by Bruce Allen, Wang-Wei Yu.

Figure 1
Figure 1. Figure 1: FIG. 1. A Mollweide projection showing sky locations of [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The posteriors for the amplitude and exponent of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Histograms of NPMV detection statistic [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Reconstruction of the mean correlation between [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

Using public pulse time-of-arrival data from five pulsar timing arrays (PTAs), we search for a stationary, isotropic, and unpolarized nHz stochastic gravitational-wave background (SGWB). This analysis is more sensitive than previous individual PTA searches because the combined 121-pulsar dataset is about four times larger than any single PTA's. For pulsars observed by multiple PTAs, we employ a new "direct combination" method to merge their astrophysical models and data. This avoids the challenge of reconciling different PTA timing models to obtain a single "best" model. A central result of our analysis is the posterior distribution of the amplitude $A_{gw}$ and exponent $\gamma_{gw}$ of the putative SGWB energy-density spectrum, modeled as a power law in frequency. While these results are consistent with a nonzero SGWB amplitude $A_{gw}$, the statistical significance-assessed via a Bayesian odds ratio and noise-marginalized false-alarm probabilities ($p$-values) for three detection statistics-remains below the conventional $5\sigma$ threshold for a confident detection. The inter-pulsar timing-residual correlation, reconstructed as a function of angle $\theta$ between the pulsar lines of sight, matches the Hellings and Downs (HD) prediction.

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 paper reports a search for a stationary isotropic unpolarized nHz stochastic gravitational-wave background using public pulse time-of-arrival data from five pulsar timing arrays. The combined dataset contains 121 pulsars; for objects observed by more than one PTA the authors introduce a direct combination procedure to merge the timing residuals and noise models without constructing a single reconciled timing solution. Posterior distributions are obtained for the power-law parameters A_gw and γ_gw of the SGWB energy-density spectrum; the results are stated to be consistent with a nonzero amplitude but remain below the conventional 5σ threshold according to Bayesian odds ratios and noise-marginalized p-values. The reconstructed angular correlation of timing residuals is reported to match the Hellings-Downs prediction.

Significance. A validated direct-combination analysis of the public multi-PTA data would constitute a meaningful increase in effective array size and sensitivity relative to any single PTA. The explicit reconstruction of the angular correlation function and the use of standard Bayesian methods on publicly available residuals are strengths that support reproducibility. The work therefore has the potential to tighten constraints on the nHz SGWB once the methodological validation is completed.

major comments (1)
  1. [Section describing the direct combination method] The direct combination method for pulsars observed by multiple PTAs is load-bearing for the central claim of a 121-pulsar analysis. The manuscript presents the procedure but does not report recovery tests on simulated data sets that inject a known Hellings-Downs signal together with realistic PTA-specific noise; consequently the quantitative effect of any residual mismatch on the reported Bayesian odds ratio and noise-marginalized p-values remains unquantified (see the section describing the direct combination method and the results in §5).
minor comments (2)
  1. [Abstract] The abstract refers to “three detection statistics” without naming them; a brief parenthetical list would improve clarity for readers who do not reach the methods section.
  2. [Throughout the text and equations] Notation for the spectral index is given as both γ_gw and γ in different equations; a single consistent symbol should be used throughout.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive review and for recognizing the potential of the combined multi-PTA analysis. We address the single major comment below in a point-by-point manner.

read point-by-point responses

  1. Referee: [Section describing the direct combination method] The direct combination method for pulsars observed by multiple PTAs is load-bearing for the central claim of a 121-pulsar analysis. The manuscript presents the procedure but does not report recovery tests on simulated data sets that inject a known Hellings-Downs signal together with realistic PTA-specific noise; consequently the quantitative effect of any residual mismatch on the reported Bayesian odds ratio and noise-marginalized p-values remains unquantified (see the section describing the direct combination method and the results in §5).

    Authors: We agree that the direct combination procedure is central to the 121-pulsar claim and that its validation via signal recovery tests is necessary to quantify any residual effects on the reported statistics. The current manuscript describes the method and applies it to the public data but does not include the requested injection-and-recovery tests on simulated datasets that incorporate PTA-specific noise. In the revised manuscript we will add these tests, injecting a Hellings-Downs correlated signal into realistic simulated timing residuals for each PTA and recovering the power-law parameters and detection statistics with the direct-combination pipeline. The results will be presented in a new subsection, with explicit quantification of any bias or degradation in the Bayesian odds ratio and noise-marginalized p-values. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is data-driven and self-contained.

full rationale

The paper's chain proceeds from public PTA timing residuals through a described direct-combination merging procedure for overlapping pulsars, followed by standard Bayesian posterior sampling on the power-law SGWB parameters A_gw and gamma_gw and reconstruction of the angular correlation function. These quantities are compared to the independently known Hellings-Downs curve. No equation reduces by construction to a fitted input, no parameter is renamed as a prediction, and no load-bearing step rests on a self-citation or author-specific uniqueness theorem. The analysis therefore remains externally falsifiable against the data and the pre-existing theoretical correlation pattern.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the power-law parametrization of the SGWB spectrum and the assumption that the direct combination method introduces no systematic bias in the recovered correlations.

free parameters (2)
  • A_gw
    Amplitude of the power-law SGWB energy-density spectrum, determined from the posterior.
  • gamma_gw
    Spectral exponent of the SGWB, determined from the posterior.
axioms (2)
  • domain assumption The stochastic gravitational-wave background is stationary, isotropic, and unpolarized.
    Explicitly stated as the model used for the search.
  • standard math The Hellings-Downs angular correlation pattern is the correct signature for an isotropic unpolarized SGWB.
    Standard general-relativity prediction invoked for comparison.

pith-pipeline@v0.9.0 · 5522 in / 1458 out tokens · 42419 ms · 2026-05-17T00:02:22.032590+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.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Prospects for multi-messenger discovery of the gravitational-wave background anisotropies via cross-correlation with galaxies

    astro-ph.CO 2026-05 unverdicted novelty 6.0

    New simulations show that cross-correlating gravitational wave background anisotropies with galaxy distributions can enable discovery at angular scales of 4-6 degrees with next-generation observatories.

  2. BB plot: A Tool for Accurate Model Selection Using Bayes factors

    gr-qc 2026-05 unverdicted novelty 6.0

    The BB plot is a new diagnostic that links Bayes factors to their expected distributions under competing hypotheses, enabling validation of calculations and low-cost background estimation for model selection in GW studies.

  3. Stochastic problems in pulsar timing

    astro-ph.HE 2026-04 unverdicted novelty 5.0

    Analytical solutions to Langevin equations for red noise and GWB in pulsars show that an Ornstein-Uhlenbeck spin frequency model is inconsistent with stationary signals, while an overdamped oscillator model and a two-...

Reference graph

Works this paper leans on

67 extracted references · 67 canonical work pages · cited by 3 Pith papers · 16 internal anchors

  1. [1]

    ref- erence PTA

    For pulsars observed by more than one PTA, a “ref- erence PTA” is selected via the arbitrary order- ing NG>EPTA>PPTA>MPTA. (InPTA is not listed because its pulsars are observed by other PTAs.) AnyotherPTAs that have data for this pulsar are called “target PTAs”

    work page

  2. [2]

    Since some pulsars have both J- and B-style names, the name is copied from the reference.parfile to the target.parfiles

    work page

  3. [3]

    Any.parfile that uses the TCB time standard is converted to TDB usingTempo2[40] tools

    work page

  4. [4]

    The parameter DM0 is initialized to the mean DMX value, and DM1 and DM2 are initialized to zero

    If the reference.parfile has a DMX model for DM fluctuations (NG and InPTA only), then it is re- placed by a Gaussian process (DMGP) model. The parameter DM0 is initialized to the mean DMX value, and DM1 and DM2 are initialized to zero

    work page

  5. [5]

    If the reference.parfile has a DMGP model, then the DM0 value is retained, and DM1 and DM2 are initialized to zero

    work page

  6. [6]

    JUMP” for time offsets and “FDx

    The astrophysical modelβ[1, Eq. 7.2], including sky position, proper motion, spin frequency, spin- down parameters, binary orbital parameters, and DM values, is copied from the reference.parfile to the target.parfiles. Detector-specific parameters in the.parfiles, for exam- ple, “JUMP” for time offsets and “FDx” for narrow-band frequency offsets, are unch...

    work page

  7. [7]

    S. R. Taylor,Nanohertz Gravitational Wave Astronomy, 1st ed. (CRC Press, Taylor & Francis Group, Boca Raton, 2021) see also arXiv:2105.13270 [astro-ph.HE], https://arxiv.org/abs/2105.13270

    work page arXiv 2021

  8. [8]

    J. D. Romano and B. Allen,Answers to frequently asked questions about the pulsar timing array Hellings and Downs curve, Classical and Quantum Gravity41, 175008 (2024), arXiv:2308.05847

    work page arXiv 2024

  9. [9]

    Agazie et al

    G. Agazie, A. Anumarlapudi, A. M. Archibald, P. T. Baker, B. B´ ecsy, L. Blecha, A. Bonilla, A. Brazier, P. R. Brook, S. Burke-Spolaor,et al.,The NANOGrav 15 yr Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational-wave Background, The Astrophysical Journal Letters952, L37 (2023), arXiv:2306.16220

    work page arXiv 2023

  10. [10]

    The Astrophysics of Nanohertz Gravitational Waves

    S. Burke-Spolaor, S. R. Taylor, M. Charisi, T. Dolch, J. S. Hazboun, A. M. Holgado, L. Z. Kelley, T. J. W. Lazio, D. R. Madison, N. McMann,et al.,The astro- physics of nanohertz gravitational waves, The Astronomy and Astrophysics Review27, 5 (2019), arXiv:1811.08826 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  11. [11]

    A Practical Theorem on Gravitational Wave Backgrounds

    E. Phinney,A Practical Theorem on Gravitational Wave Backgrounds, arXiv preprint astro-ph/0108028 (2001)

    work page internal anchor Pith review Pith/arXiv arXiv 2001

  12. [12]

    D. J. Reardonet al.(PPTA Collaboration),Search for an Isotropic Gravitational-wave Background with the Parkes Pulsar Timing Array, The Astrophysical Journal Letters 951, L6 (2023), arXiv:2306.16215

    work page internal anchor Pith review Pith/arXiv arXiv 2023

  13. [13]

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

    J. Antoniadiset al.(EPTA and InPTA Collaborations), The second data release from the European Pulsar Timing Array: III. Search for gravitational wave signals, Astron- omy & Astrophysics678, A50 (2023), arXiv:2306.16214

    work page internal anchor Pith review Pith/arXiv arXiv 2023

  14. [14]

    H. Xuet al.(CPTA Collaboration),Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I, Research in Astronomy and Astrophysics23, 075024 (2023), arXiv:2306.16216

    work page internal anchor Pith review Pith/arXiv arXiv 2023

  15. [15]

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

    G. Agazieet al.(NANOGrav Collaboration), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background, The Astrophysical Journal Letters951, L8 (2023), arXiv:2306.16213

    work page internal anchor Pith review Pith/arXiv arXiv 2023

  16. [16]

    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,et al.,The MeerKAT Pul- sar Timing Array: The first search for gravitational waves with the MeerKAT radio telescope, Monthly No- tices of the Royal Astronomical Society536, 1489 (2025), arXiv:2412.01153

    work page arXiv 2025

  17. [17]

    van Haasteren,Use model averaging instead of model selection in pulsar timing, Monthly Notices of the Royal Astronomical Society: Letters537, L1 (2024), arXiv:2409.06050

    R. van Haasteren,Use model averaging instead of model selection in pulsar timing, Monthly Notices of the Royal Astronomical Society: Letters537, L1 (2024), arXiv:2409.06050

    work page arXiv 2024

  18. [18]

    Allen, S

    B. Allen, S. Dhurandhar, Y. Gupta, M. McLaughlin, P. Natarajan, R. M. Shannon, E. Thrane, and A. Vec- chio, The International Pulsar Timing Array checklist for the detection of nanohertz gravitational waves (2023), arXiv:2304.04767 [astro-ph.IM]

    work page arXiv 2023

  19. [19]

    J. M. Goldstein, J. Veitch, A. Sesana, and A. Vecchio, Null-stream analysis of pulsar timing array data: local- ization of resolvable gravitational wave sources, Monthly Notices of the Royal Astronomical Society477, 5447 (2018), arXiv:1712.03975

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  20. [20]

    Perera and others (IPTA Consortium),The International Pulsar Timing Array: Second data release,Monthly Notices of the Royal Astronomical Society490(2019) 4666 [1909.04534]

    B. Perera, M. DeCesar, P. Demorest, M. Kerr, L. Lentati, D. Nice, S. Os lowski, S. Ransom, M. Keith, Z. Arzou- manian,et al.,The International Pulsar Timing Array: second data release, Monthly Notices of the Royal Astro- nomical Society490, 4666 (2019), arXiv:1909.04534

    work page arXiv 2019

  21. [21]

    J. G. Baier, J. S. Hazboun, and J. D. Romano,A sensi- tivity curve approach to tuning a pulsar timing array in the detection era, Class. Quant. Grav.42, 075008 (2025), arXiv:2409.00336 [astro-ph.HE]

    work page arXiv 2025

  22. [22]

    EPTA Collaboration, EPTA data set DR2new, 27.4.2023, commit be91c6be8623df1ab8b66d5e36281005c6ca6851https: //gitlab.in2p3.fr/epta/epta-dr2/-/tree/master/ EPTA-DR2/DR2new. As this draft was being completed, we found that there is a later commit dated 23.8.2025 with hash 1506123cae38bd2fb88df53d28f87cdb87c2f389, but have not yet investigated the differences

    work page 2023

  23. [23]

    InPTA Collaboration, InPTA data release 1, commit 2c400d51428abd59d6cf842cd8fe7c840e819d0d,https:// github.com/inpta/InPTA.DR1

    work page

  24. [24]

    MPTA Collaboration, MPTA 4.5-yr data set,https://docs.datacentral.org.au/ meerkat-pulsar-timing-array/45-year/ accessing-the-data/

    work page

  25. [25]

    NANOGrav Collaboration, NANOGrav data set 15yr, v1.0.0,https://zenodo.org/records/7967585

    work page arXiv

  26. [26]

    PPTA Collaboration, PPTA data release 3, https://data.csiro.au/collection/csiro%3A59381v1 (2023), cSIRO Data Collection

    work page 2023

  27. [27]

    van Haasteren and W.-W

    R. van Haasteren and W.-W. Yu, The direct com- bination method for PTA data analysis (2025), work in progress,https://github.com/vhaasteren/ metapulsar/blob/main/docs/METHOD_DESCRIPTION.md

    work page 2025

  28. [28]

    van Haasteren and W.-W

    R. van Haasteren and W.-W. Yu, MetaPulsar (2025), zenodo,https://doi.org/10.5281/zenodo.17626664

    work page doi:10.5281/zenodo.17626664 2025

  29. [29]

    J. A. Ellis, M. Vallisneri, S. R. Taylor, and P. T. Baker, Enterprise: Enhanced Numerical Toolbox Enabling a Ro- bust PulsaR Inference SuitE, Astrophysics Source Code Library , ascl (2019)

    work page 2019

  30. [30]

    NANOGrav Discovery Pipeline,https://github.com/ nanograv/discovery(2025), accessed: 2025-11-10

    work page 2025

  31. [31]

    R. W. Hellings, Pulsar timing and gravitational waves, inRelativistic Gravitational Experiments in Space: pro- ceedings of NASA Workshop, Annapolis, Maryland, June 28-30, 1988, NASA Conference Publication, Vol. 3046 (1989) pp. 93–97

    work page 1988

  32. [32]

    Allen,Variance of the Hellings-Downs corre- lation, Physical Review D107, 043018 (2023), arXiv:2205.05637

    B. Allen,Variance of the Hellings-Downs corre- lation, Physical Review D107, 043018 (2023), arXiv:2205.05637

    work page arXiv 2023

  33. [33]

    van Haasteren, B

    R. van Haasteren, B. Allen, and J. D. Romano, Opti- mal robust detection statistics for pulsar timing arrays (2025), arXiv:2509.06489 [astro-ph.IM]

    work page arXiv 2025

  34. [34]

    S. J. Vigeland, K. Islo, S. R. Taylor, and J. A. Ellis,Noise-marginalized optimal statistic: A robust hybrid frequentist-Bayesian statistic for the stochastic gravitational-wave background in pulsar timing arrays, Physical Review D98, 044003 (2018), arXiv:1805.12188

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  35. [35]

    Vallisneri, P

    M. Vallisneri, P. M. Meyers, K. Chatziioannou, and A. J. Chua,Posterior predictive checking for gravitational- wave detection with pulsar timing arrays. I. the opti- mal statistic, Physical Review D108, 123007 (2023), arXiv:2306.05558. 8

    work page arXiv 2023

  36. [36]

    N. J. Cornish and L. Sampson,Towards robust gravita- tional wave detection with pulsar timing arrays, Physical Review D93, 104047 (2016), arXiv:1512.06829

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  37. [37]

    All correlations must die: Assessing the significance of a stochastic gravitational-wave background in pulsar-timing arrays

    S. Taylor, L. Lentati, S. Babak, P. Brem, J. Gair, A. Sesana, and A. Vecchio,All correlations must die: As- sessing the significance of a stochastic gravitational-wave background in pulsar timing arrays, Physical Review D 95, 042002 (2017), arXiv:1606.09180

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  38. [38]

    van Haasteren,Calculation of p-values for quadratic statistics in pulsar timing arrays, Physical Review D112, 103009 (2025), arXiv:2506.10811 [astro-ph.IM]

    R. van Haasteren,Calculation of p-values for quadratic statistics in pulsar timing arrays, Physical Review D112, 103009 (2025), arXiv:2506.10811 [astro-ph.IM]

    work page arXiv 2025

  39. [39]

    J. S. Hazboun, P. M. Meyers, J. D. Romano, X. Siemens, and A. M. Archibald,Analytic distribu- tion of the optimal cross-correlation statistic for stochas- tic gravitational-wave-background searches using pulsar timing arrays, Physical Review D108, 104050 (2023), arXiv:2305.01116. [34]The nanograv 15 yr dataset: Posterior predictive checks for gravitationa...

    work page arXiv 2023

  40. [40]

    The dataset and timing analysis,Astron

    J. Antoniadis, S. Babak, A.-S. B. Nielsen, C. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. Brook, M. Burgay, R. Caballero,et al.,The second data release from the European Pulsar Timing Array-I. The dataset and timing analysis, Astronomy & Astrophysics678, A48 (2023), arXiv:2306.16224

    work page arXiv 2023

  41. [41]

    Tarafdar, K

    P. Tarafdar, K. Nobleson, P. Rana, J. Singha, M. Kr- ishnakumar, B. C. Joshi, A. K. Paladi, N. Kolhe, N. D. Batra, N. Agarwal,et al.,The Indian Pulsar Timing Ar- ray: First data release, Publications of the Astronomical Society of Australia39, e053 (2022), arXiv:2206.09289

    work page arXiv 2022

  42. [42]

    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,et al.,The MeerKAT Pulsar Timing Array: the 4.5-yr data release and the noise and stochas- tic signals of the millisecond pulsar population, Monthly Notices of the Royal Astronomical Society536, 1467 (2025), arXiv:2412.01148

    work page arXiv 2025

  43. [43]

    Agazie et al

    G. Agazie, M. F. Alam, A. Anumarlapudi, A. M. Archibald, Z. Arzoumanian, P. T. Baker, L. Blecha, V. Bonidie, A. Brazier, P. R. Brook,et al.,The NANOGrav 15 yr data set: observations and timing of 68 millisecond pulsars, The Astrophysical Journal Let- ters951, L9 (2023), arXiv:2306.16217

    work page arXiv 2023

  44. [44]

    A. Zic, D. J. Reardon, A. Kapur, G. Hobbs, R. Mandow, M. Cury lo, R. M. Shannon, J. Askew, M. Bailes, N. R. Bhat,et al.,The Parkes Pulsar Timing Array third data release, Publications of the Astronomical Society of Aus- tralia40, e049 (2023), arXiv:2306.16230

    work page arXiv 2023

  45. [45]

    TEMPO2, a new pulsar timing package. I: Overview

    G. Hobbs, R. Edwards, and R. Manchester,TEMPO2, a new pulsar timing package. 1. Overview, Mon. Not. R. Astron. Soc.369, 655 (2006), arXiv:astro-ph/0603381

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  46. [46]

    Goncharov, S

    B. Goncharov, S. Sardana, A. Sesana, S. M. Tomson, J. Antoniadis, A. Chalumeau, D. J. Champion, S. Chen, E. F. Keane, K. Liu,et al.,Reading signatures of su- permassive binary black holes in pulsar timing array observations, Nature Communications16, 9692 (2025), arXiv:2409.03627

    work page arXiv 2025

  47. [47]

    Efficient and Modular Im- plicit Differentiation

    M. Blondel, Q. Berthet, M. Cuturi, R. Frostig, S. Hoyer, F. Llinares-L´ opez, F. Pedregosa, and J.-P. Vert,Effi- cient and modular implicit differentiation, Advances in neural information processing systems35, 5230 (2022), arXiv:2105.15183

    work page arXiv 2022

  48. [48]

    Crisostomi, R

    M. Crisostomi, R. van Haasteren, P. M. Meyers, and M. Vallisneri, Beyond diagonal approximations: im- proved covariance modeling for pulsar timing array data analysis (2025), arXiv:2506.13866 [astro-ph.IM]

    work page arXiv 2025

  49. [49]

    Pulsar timing array analysis in a Legendre polynomial basis

    B. Allen, A. L. von Blanckenburg, and K. D. Olum, Pul- sar timing array analysis in a Legendre polynomial basis (2025), arXiv:2510.05913 [gr-qc]

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  50. [50]

    D. Nice, P. Demorest, I. Stairs, R. Manchester, J. Tay- lor, W. Peters, J. Weisberg, A. Irwin, N. Wex, and Y. Huang, TEMPO: Pulsar timing data analysis, Astro- physics Source Code Library, ascl:1509.002 (2015)

    work page 2015

  51. [51]

    J. Luo, S. Ransom, P. Demorest, P. S. Ray, A. M. Archibald, M. Kerr, R. J. Jennings, M. Bachetti, R. van Haasteren, C. A. Champagne, J. Colen, C. Phillips, J. Zimmerman, K. Stovall, M. T. Lam,et al.,PINT: A modern software package for pulsar timing, Astrophys. J. 911, 45 (2021), arXiv:2012.00074 [astro-ph.IM]

    work page arXiv 2021

  52. [52]

    Bradbury, R

    J. Bradbury, R. Frostig, P. Hawkins, M. J. Johnson, C. Leary, D. Maclaurin, G. Necula, A. Paszke, J. Van- derPlas, S. Wanderman-Milne, and Q. Zhang, JAX: com- posable transformations of Python+NumPy programs (2018)

    work page 2018

  53. [53]

    R. M. Nealet al.,MCMC using Hamiltonian dynamics, Handbook of Markov Chain Monte Carlo2, 2 (2011), arXiv:1206.1901

    work page arXiv 2011

  54. [54]

    M. D. Hoffman, A. Gelman,et al.,The No-U-Turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo., J. Mach. Learn. Res.15, 1593 (2014), arXiv:1111.4246

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  55. [55]

    D. Phan, N. Pradhan, and M. Jankowiak,Compos- able effects for flexible and accelerated probabilistic pro- gramming in numpyro, arXiv preprint arXiv:1912.11554 (2019)

    work page internal anchor Pith review Pith/arXiv arXiv 1912

  56. [56]

    Pyro: Deep Universal Probabilistic Programming

    E. Bingham, J. P. Chen, M. Jankowiak, F. Obermeyer, N. Pradhan, T. Karaletsos, R. Singh, P. A. Szerlip, P. Horsfall, and N. D. Goodman,Pyro: Deep universal probabilistic programming, J. Mach. Learn. Res.20, 28:1 (2019), arXiv:1810.09538

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  57. [57]

    A Conceptual Introduction to Hamiltonian Monte Carlo

    M. Betancourt,A conceptual introduction to Hamiltonian Monte Carlo, arXiv preprint arXiv:1701.02434 (2017)

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  58. [58]

    Ellis and R

    J. Ellis and R. van Haasteren, jellis18/ptmcmcsampler: Official release (2017)

    work page 2017

  59. [59]

    Allen and J

    B. Allen and J. D. Romano,Hellings and Downs corre- lation of an arbitrary set of pulsars, Physical Review D 108, 043026 (2023), arXiv:2208.07230

    work page arXiv 2023

  60. [60]

    Agazieet al.[International Pulsar Timing Array],Comparing recent PTA results on the nanohertz stochastic gravitational wave background,Astrophys

    G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzouma- nian, J. Askew, S. Babak, M. Bagchi, and et al. (Inter- national Pulsar Timing Array),Comparing Recent Pul- sar Timing Array Results on the Nanohertz Stochastic Gravitational-wave Background, Astrophys. J.966, 105 (2024), arXiv:2309.00693 [astro-ph.HE]

    work page arXiv 2024

  61. [61]

    Pernet,Null hypothesis significance testing: a short tutorial, F1000Research4, 621 (2016)

    C. Pernet,Null hypothesis significance testing: a short tutorial, F1000Research4, 621 (2016)

    work page 2016

  62. [62]

    Jeffreys,Theory of Probability, 3rd ed., International Series of Monographs on Physics (Clarendon Press, Ox- ford, 1961)

    H. Jeffreys,Theory of Probability, 3rd ed., International Series of Monographs on Physics (Clarendon Press, Ox- ford, 1961)

    work page 1961

  63. [63]

    E. M. Zahraoui, P. Maturana-Russel, W. van Straten, R. Meyer, and S. Gulyaev,Generalized stepping-stone sampling: efficient marginal likelihood estimation in gravitational wave analysis of pulsar timing array data, Monthly Notices of the Royal Astronomical Society540, 3818 (2025), arXiv:2411.14736 [astro-ph.IM]

    work page arXiv 2025

  64. [64]

    E. M. Zahraoui, P. Maturana-Russel, W. van Straten, R. Meyer, and S. Gulyaev, 9 https://github.com/apokryphav1/gss- estimator/tree/main (2024)

    work page 2024

  65. [65]

    The code we used in the paper will become a pull request inDiscovery

    work page

  66. [66]

    Allen and J

    B. Allen and J. D. Romano,Optimal reconstruction of the Hellings and Downs correlation, Physical Review Letters 134, 031401 (2025), arXiv:2407.10968

    work page arXiv 2025

  67. [67]

    NANOGrav Collaboration, NANOGrav 15-year stochas- tic analysis repository,https://github.com/nanograv/ 15yr_stochastic_analysis/blob/main/tutorials/ optimal_statistic_covariances.py(2023), file: tutorials/optimal statistic covariances.py, com- mit: 14ed645b7d26be788340036ded29b3c66286e29c

    work page 2023