Full analytic expressions of overlap reduction functions for anisotropies of the stochastic gravitational-wave background with pulsar timing arrays
Phys. Rev. D 114, 044003 – Published 3 August, 2026
DOI: https://doi.org/10.1103/34y2-pn7d
Abstract
Recent findings collectively suggest that pulsar timing arrays (PTAs) have detected the stochastic gravitational-wave background in the nanohertz frequency band, opening new avenues for testing theories of gravity and cosmology. PTAs serve as a unique tool for probing the nanohertz stochastic gravitational-wave background, enabling studies of supermassive black hole evolution and early-Universe physics. PTA data analysis fundamentally relies on cross-correlating timing residuals from pulsar pairs, for which the overlap reduction functions (ORFs) are critical in determining the sensitivity of such analyses. Conventional ORF calculations based on the “short-wavelength approximation” break down when handling the scalar longitudinal polarization mode and fail to accurately describe frequency-dependent effects and anisotropies. This work presents a significant advance by rigorously deriving the full response functions within an analytical cross-correlation framework. We systematically reveal, for the first time, the intrinsic symmetries and mathematical relationships among the anisotropic ORF integrals for different polarization modes. Through variable substitutions and coordinate rotations, we transform the complex integrals into analytically tractable forms, resolving the divergence difficulties encountered in the vector and scalar longitudinal modes. Building on these theoretical advances, we establish a universal framework capable of deriving fully analytical expressions for anisotropic ORFs to arbitrary order for all six polarization modes. As a direct application of this framework, we provide complete analytical expressions for anisotropic ORFs up to the order of . This universal framework is free from approximations; its (0,0) component naturally corresponds to the isotropic case and recovers the well-known Hellings-Downs curve. Compared to numerical integration, our results offer broader applicability, significantly faster computation, and higher precision. They provide a valuable theoretical and computational foundation for conducting high-precision anisotropic sky mapping, polarization-mode separation, and searches for new physics using PTA data.
