# Accurate pulsar timing array likelihoods: better low-rank covariance for gravitational wave background

Preprint · Heling Deng, Ken D. Olum, Rutger van Haasteren

- Archive id: arxiv:2610.33444
- Published: 2026-10-24
- arXiv: https://arxiv.org/abs/2610.33444
- Citations: 0
- Source URL: https://arxiv.org/abs/2610.33444
- PDF: https://arxiv.org/pdf/2610.33444
- Provenance: arxiv

## Abstract

In pulsar timing array analyses, the stochastic process describing the gravitational wave background is typically modeled using a Fourier basis at discrete frequencies with an interval $\\sim 0.1~\\mathrm{yr}^{-1}$, with a diagonal Fourier-domain covariance determined by the power spectral density at the corresponding frequencies. Although this rank-reduced modeling significantly accelerates likelihood evaluation and provides a good approximation in many scenarios, it can lead to biased parameter inference because it is insensitive to spectral features on scales smaller than $\\sim 0.1~\\mathrm{yr}^{-1}$. We introduce a better construction of the low-rank covariance for low-frequency processes by maximizing the log-likelihood averaged over an ensemble of realizations drawn from the true model likelihood. In this approach, the basis-domain covariance matrix is sensitive to the entire model spectrum and is generally non-diagonal. Using two mock datasets, we demonstrate improved likelihood accuracy and more faithful recovery of spectral features compared to the conventional Fourier-diagonal model. Benchmark tests show that this new model requires a computational cost comparable to that of the conventional approach.

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

Topics: [Gravitational-Wave Astronomy](https://space-kitty.com/topics/gravitational-waves)

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