pith. sign in

arxiv: 2510.05913 · v4 · submitted 2025-10-07 · 🌀 gr-qc · astro-ph.CO· astro-ph.IM

Pulsar timing array analysis in a Legendre polynomial basis

Bruce Allen , Arian L. von Blanckenburg , Ken D. Olum This is my paper

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

classification 🌀 gr-qc astro-ph.COastro-ph.IM
keywords pulsar timing arraysgravitational wave backgroundLegendre polynomialsHellings-Downs correlationpower-law spectraquadratic estimatortiming residuals
0
0 comments X

The pith

Switching to Legendre polynomials for pulsar timing data simplifies modeling of pulsar effects and delivers analytic expressions for the Hellings-Downs estimator when spectra are power laws.

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

The paper replaces the standard Fourier-mode basis with Legendre polynomials to analyze signals in pulsar timing arrays. This basis makes it straightforward to account for pulsar modeling, which removes constant, linear, and quadratic trends in the timing residuals and therefore sets the amplitudes of the first three Legendre polynomials to zero. The authors construct an optimal quadratic cross-correlation estimator for the Hellings and Downs correlation pattern and derive its variance following the Allen-Romano approach. When the gravitational-wave background and pulsar noise power spectra are power laws or sums of power laws, the Legendre basis supplies closed-form analytic results for these and related quantities. A reader cares because the closed forms can reduce computational cost and increase transparency when searching for a stochastic gravitational-wave background.

Core claim

By modeling pulsar timing residuals in a Legendre polynomial basis rather than trigonometric Fourier modes, the first three coefficients are automatically zeroed to incorporate pulsar modeling effects. This framework yields an optimal quadratic estimator mu-hat for the Hellings-Downs correlation together with an analytic expression for its variance sigma-squared when the gravitational-wave background and pulsar noise spectra are (sums of) power laws in frequency.

What carries the argument

The Legendre polynomial basis for timing residuals, which replaces Fourier modes and directly encodes the removal of low-order polynomial trends by zeroing the first three modes.

Load-bearing premise

The gravitational-wave background and pulsar noise power spectra must be sums of power laws in frequency.

What would settle it

Direct numerical evaluation of the estimator variance for a non-power-law spectrum should deviate from the closed-form expression derived in the Legendre basis.

Figures

Figures reproduced from arXiv: 2510.05913 by Arian L. von Blanckenburg, Bruce Allen, Ken D. Olum.

Figure 1
Figure 1. Figure 1: FIG. 1. The Legendre polynomials [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The transmission function [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Comparison of sinc [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Starting from timing residuals [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
read the original abstract

We use Legendre polynomials, previously employed in this context by Lee et al., van Haasteren and Levin, and Pitrou and Cusin, to model signals in pulsar timing arrays. These replace the (Fourier mode) basis of trigonometric functions normally used for data analysis. The Legendre basis makes it simpler to incorporate pulsar modeling effects, which remove constant-, linear-, and quadratic-in-time terms from pulsar timing residuals. In the Legendre basis, this zeroes the amplitudes of the the first three Legendre polynomials. We use this basis to construct an optimal quadratic cross-correlation estimator $\widehat{\mu}$ of the Hellings and Downs (HD) correlation and compute its variance $\sigma^2_{\widehat{\mu}}$ in the way described by Allen and Romano. Remarkably, if the gravitational-wave background (GWB) and pulsar noise power spectra are (sums of) power laws in frequency, then in this basis one obtains analytic closed forms for many quantities of interest.

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 proposes replacing the conventional Fourier basis with Legendre polynomials for modeling signals in pulsar timing arrays. This choice naturally incorporates the effects of pulsar timing models by zeroing the amplitudes of the first three Legendre polynomials, which absorb constant, linear, and quadratic terms in the timing residuals. The authors construct an optimal quadratic cross-correlation estimator for the Hellings-Downs correlation and derive its variance following the Allen-Romano formalism, obtaining analytic closed-form expressions for several quantities when the gravitational-wave background and pulsar noise power spectra are (sums of) power laws in frequency.

Significance. If the closed-form results hold under the stated power-law assumptions, the Legendre basis offers a computationally efficient alternative for PTA analyses that avoids numerical integration over frequency for common spectral shapes. This builds directly on prior Legendre-basis work cited in the abstract and could simplify variance calculations for the HD correlation estimator without introducing additional free parameters.

major comments (2)
  1. [§4] §4 (variance derivation): The claim of analytic closed forms for σ²_μ̂ relies on the power-law assumption for the spectra, but the explicit evaluation of the frequency integrals over the Legendre basis functions (which should reduce to Gamma functions or similar) is not shown; without this step the central simplification cannot be verified independently.
  2. [§5] §5 (validation): No comparison of the closed-form variance against numerical integration or simulated PTA data is presented, even for a simple single power-law case; this leaves the practical accuracy of the analytic expressions untested within the manuscript.
minor comments (2)
  1. [Abstract] Abstract: the phrase 'many quantities of interest' is vague; explicitly listing the quantities (e.g., the estimator variance, the HD overlap reduction function projection) would improve clarity.
  2. [§3] Notation: the definition of the quadratic estimator μ̂ should be written with an explicit sum over pulsar pairs and Legendre coefficients to make the zeroing of the first three modes immediately visible.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable feedback on our manuscript. We address each of the major comments below.

read point-by-point responses

  1. Referee: [§4] §4 (variance derivation): The claim of analytic closed forms for σ²_μ̂ relies on the power-law assumption for the spectra, but the explicit evaluation of the frequency integrals over the Legendre basis functions (which should reduce to Gamma functions or similar) is not shown; without this step the central simplification cannot be verified independently.

    Authors: We appreciate this observation. The manuscript asserts that closed-form expressions are obtained for power-law spectra, but we concur that the explicit computation of the integrals over the Legendre basis functions was not presented in detail. In the revised manuscript, we will add an appendix detailing the evaluation of these frequency integrals, showing their reduction to expressions involving Gamma functions for the assumed power-law forms of the spectra. This will allow independent verification of the analytic results. revision: yes

  2. Referee: [§5] §5 (validation): No comparison of the closed-form variance against numerical integration or simulated PTA data is presented, even for a simple single power-law case; this leaves the practical accuracy of the analytic expressions untested within the manuscript.

    Authors: We acknowledge the lack of explicit validation in the current version. To strengthen the manuscript, we will include in the revised §5 a comparison between the analytic closed-form variance and numerical integration results for a representative single power-law spectrum. If feasible, we will also present results from simulated PTA data to further test the expressions. This addition will demonstrate the practical accuracy of the derived formulas. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's central results are analytic closed-form expressions for the quadratic estimator and its variance, obtained by substituting the Legendre basis into the standard cross-correlation formalism and integrating against power-law spectra. This reduction is explicitly conditional on the external power-law assumption for GWB and noise spectra (stated in the abstract), which is not derived or fitted within the paper. The Legendre basis itself is adopted from prior literature (Lee et al., van Haasteren & Levin, Pitrou & Cusin) and correctly nulls the first three modes to absorb the quadratic timing model; the citation to Allen & Romano merely describes the general quadratic-estimator construction rather than supplying a load-bearing uniqueness theorem or ansatz. No step equates a prediction to its own fitted input or renames a known result as a derivation. The derivation chain is therefore self-contained once the stated power-law assumption is granted.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the orthogonality and completeness properties of Legendre polynomials together with the domain assumption that timing models remove the lowest three modes and the power-law spectral assumption that enables closed forms.

axioms (2)
  • standard math Legendre polynomials form an orthogonal basis on the relevant interval and can represent the same time-series information as Fourier modes
    Invoked to replace the trigonometric basis while preserving completeness.
  • domain assumption Pulsar timing models subtract constant, linear, and quadratic terms, which corresponds exactly to zeroing the first three Legendre polynomial amplitudes
    Central simplification stated in the abstract.

pith-pipeline@v0.9.0 · 5703 in / 1432 out tokens · 40964 ms · 2026-05-18T09:26:14.403206+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. Constraints on the Primordial Black Hole Abundance using Pulsar Parameter Drifts

    astro-ph.CO 2026-04 unverdicted novelty 8.0

    The first search for scalar-induced gravitational waves via pulsar parameter drifts yields f_PBH < 10^{-10} (95% CL) for PBH masses 0.3 to 4e4 solar masses, strongly disfavoring a primordial black hole origin for LVK ...

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

    astro-ph.CO 2025-12 conditional novelty 6.0

    Combined five-PTA dataset yields posterior on SGWB power-law amplitude and index consistent with nonzero signal but below 5-sigma significance, with reconstructed angular correlations matching the Hellings-Downs prediction.

  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

21 extracted references · 21 canonical work pages · cited by 3 Pith papers

  1. [1]

    K. J. Lee, C. G. Bassa, G. H. Janssen, R. Karuppusamy, M. Kramer, R. Smits, and B. W. Stappers, The optimal schedule for pulsar timing array observations, Mon. Not. Roy. Astron. Soc.423, 2642 (2012), see Appendix A

    work page 2012

  2. [2]

    van Haasteren and Y

    R. van Haasteren and Y. Levin, Understanding and analysing time-correlated stochastic signals in pulsar tim- ing, Mon. Not. Roy. Astron. Soc.428, 1147 (2013)

    work page 2013

  3. [3]

    Pitrou and G

    C. Pitrou and G. Cusin, Mitigating cosmic variance in the Hellings-Downs curve: A cosmic microwave background analogy, Phys. Rev. D111, 083544 (2025)

    work page 2025

  4. [4]

    Allen and J

    B. Allen and J. D. Romano, Optimal reconstruction of the Hellings and Downs correlation, Phys. Rev. Lett. 134, 031401 (2025)

    work page 2025

  5. [5]

    Antoniadis, P

    J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A.-S. B. Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau,et al.(EPTA and InPTA Collaborations), The second data release from the European Pulsar Timing Array: III. Search for gravitational wave signals, Astronomy & ...

    work page 2023

  6. [6]

    Agazie, A

    G. Agazie, A. Anumarlapudi, A. M. Archibald, Z. Ar- zoumanian, P. T. Baker, B. Becsy, L. Blecha, A. Brazier, P. R. Brook, S. Burke-Spolaor, R. Burnette, R. Case, M. Charisi, S. Chatterjee, K. Chatziioannou,et al. (NANOGrav Collaboration), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background, The Astrophysical Journal Letters951, L8 (2023)

    work page 2023

  7. [7]

    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,et al.(CPTA Collab- oration), Searching for the Nano-Hertz Stochastic Gravi- tational Wave Background with the Chinese Pulsar Tim- ing Array Data Release I, Research in Astronomy and 17 Astrophysics23, 075024 (2023)

    work page 2023

  8. [8]

    A. Zic, D. J. Reardon, A. Kapur, G. Hobbs, R. Mandow, M. Cury lo, R. M. Shannon, J. Askew, M. Bailes, N. D. R. Bhat, A. Cameron, Z.-C. Chen, S. Dai, V. D. Marco, Y. Feng,et al., The Parkes Pulsar Timing Array third data release, Publ. Astron. Soc. Austral.40, e049 (2023)

    work page 2023

  9. [9]

    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,et al., The MeerKAT Pulsar Timing Array: the first search for grav- itational waves with the MeerKAT radio telescope, Mon. Not. Roy. Astron. Soc.536, 1489 (2025)

    work page 2025

  10. [10]

    R. W. Hellings and G. S. Downs, Upper limits on the isotropic gravitational radiation background from pulsar timing analysis, Astrophys. J.265, L39 (1983)

    work page 1983

  11. [11]

    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

  12. [12]

    K. A. Gersbach, S. R. Taylor, P. M. Meyers, and J. D. Romano, Spatial and spectral characterization of the gravitational-wave background with the PTA optimal statistic, Phys. Rev. D111, 023027 (2025)

    work page 2025

  13. [13]

    J. S. Hazboun, J. D. Romano, and T. L. Smith, Realistic sensitivity curves for pulsar timing arrays, Phys. Rev. D 100, 104028 (2019)

    work page 2019

  14. [14]

    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, record ascl:1509.002 (2015)

    work page 2015

  15. [15]

    Hobbs, R

    G. Hobbs, R. Edwards, and R. Manchester, TEMPO2, a new pulsar timing package. 1. Overview, Mon. Not. Roy. Astron. Soc.369, 655 (2006)

    work page 2006

  16. [16]

    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)

    work page 2021

  17. [17]

    Gradshteyn and I

    I. Gradshteyn and I. Ryzhik,Table of Integrals, Series and Products, 7th Edition(Academic Press, 2007)

    work page 2007

  18. [18]

    F. W. Olver, A. B. O. Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain, NIST Digital Library of Mathematical Functions,https:// dlmf.nist.gov/10.60.E12(2025), equation (10.60.E12)

    work page 2025

  19. [19]

    J. Nay, K. K. Boddy, T. L. Smith, and C. M. F. Min- garelli, Harmonic analysis for pulsar timing arrays, Phys. Rev. D110, 044062 (2024)

    work page 2024

  20. [20]

    Allen and J

    B. Allen and J. D. Romano, Harmonic spectrum of pulsar timing array angular correlations (2024), arXiv:2412.14852 [astro-ph.IM]

    work page arXiv 2024

  21. [21]

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

    work page 2024