# Pulsar timing array analysis in a Legendre polynomial basis

Preprint · Bruce Allen, Arian L. von Blanckenburg, Ken D. Olum

- Archive id: arxiv:2510.05913
- Published: 2025-10-07
- Venue: PRD 113, 102001 (2026)
- DOI: https://doi.org/10.1103/2dft-4rjj
- arXiv: https://arxiv.org/abs/2510.05913
- Source URL: https://arxiv.org/abs/2510.05913
- PDF: https://arxiv.org/pdf/2510.05913
- Provenance: arxiv

## 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μ$ of the Hellings and Downs (HD) correlation and compute its variance $σ^2\_{\\widehatμ}$ 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.

_Abstract reproduced under CC0-1.0 (arXiv descriptive metadata)._

---
Space Kitty Archive · https://space-kitty.com · API: https://space-kitty.com/api/v1 · MCP: https://space-kitty.com/mcp · Premium: $0.20/request via x402 (Base Sepolia (testnet)). Metadata: arXiv (CC0), OpenAlex (CC0), Europe PMC (CC BY/CC0 works), ClinicalTrials.gov.