# Beyond diagonal approximations: improved covariance modeling for pulsar timing array data analysis

Preprint · Marco Crisostomi, Rutger van Haasteren, Patrick M. Meyers, Michele Vallisneri

- Archive id: arxiv:2506.13866
- Published: 2025-06-16
- arXiv: https://arxiv.org/abs/2506.13866
- Source URL: https://arxiv.org/abs/2506.13866
- PDF: https://arxiv.org/pdf/2506.13866
- Provenance: arxiv

## Abstract

Pulsar Timing Array (PTA) searches for nHz gravitational-wave backgrounds (GWBs) typically model time-correlated noise by assuming a diagonal covariance in Fourier space, neglecting inter-frequency correlations introduced by the finite observation window. We show that this diagonal approximation can produce method-dependent shifts in estimates of spectral parameters, especially for pulsars with shallow intrinsic red noise spectra. To address these limitations, we present a method that (i) computes the time-domain autocorrelation on a coarse grid using a fast Fourier transform (FFT), (ii) interpolates it accurately to the unevenly sampled observation times, and (iii) incorporates it into a low-rank likelihood via the Sherman--Morrison--Woodbury identity. Using both analytic covariance comparisons and end-to-end simulations inspired by the NANOGrav 15-year dataset, we demonstrate that our method captures frequency correlations faithfully, avoids Gibbs ringing, and yields estimates closer to the injected values in the realization studied with modest computational cost. As PTA datasets increase in sensitivity and complexity, our approach offers a practical and scalable path to fully accurate covariance modeling for current and future analyses.

_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.