Analytical Estimates of Gravitational Wave Background Anisotropies from Shot Noise and Large-Scale Structure in Pulsar Timing Arrays

An important next step for pulsar timing arrays (PTAs) is to measure anisotropies in the gravitational wave background (GWB) at $\sim$ nano-Hz frequencies. We calculate the expected GWB anisotropies using empirically calibrated models for the merger rates of supermassive black hole binaries (SMBHBs). The anisotropies reflect both shot-noise in the discrete SMBHB populations while also tracing, in part, the large-scale structure (LSS) of the universe. The shot-noise term is sensitive to the high-mass end of the merging SMBH mass function, depends somewhat on the low-redshift tail of the merger distribution, and is a strong function of observing frequency. The precise frequency dependence provides a test of SMBHB residence times. In our models, the mean shot-noise anisotropy typically lies close to or above the broad frequency-band NANOGrav upper limits. Consequently, near-future PTA data, and potentially re-analyses of existing measurements using frequency-dependent shot-noise anisotropy templates, should be capable of detecting this signal or placing meaningful constraints on SMBHB merger models. A full interpretation, however, will require modeling the probability distribution of shot-noise amplitudes rather than relying solely on ensemble-averaged predictions. The LSS-induced anisotropies are at least two to three orders of magnitude smaller. Although the LSS contribution contains valuable information regarding the redshift distribution and clustering bias of the merging SMBHBs, detecting this component will be challenging.

💡 Research Summary

The paper presents a comprehensive analytical framework for predicting anisotropies in the nanohertz stochastic gravitational‑wave background (GWB) as measured by pulsar timing arrays (PTAs). The authors identify two distinct sources of anisotropy: (i) shot‑noise arising from the discrete, Poissonian nature of the supermassive black‑hole binary (SMBHB) population, and (ii) large‑scale‑structure (LSS) fluctuations because SMBHB mergers trace the underlying matter density field.

Starting from the definition of the strain‑squared field h²(f, Ω̂), they construct the fractional fluctuation δh²(Ω̂) and its two‑point correlation function ωₕ²(θ). By expanding δh² in spherical harmonics they obtain the angular power spectrum Cℓ,ₕ², which naturally separates into a shot‑noise term C_SN,ℓ and an LSS term C_LSS,ℓ. The shot‑noise amplitude is expressed through an effective source density N_eff = ⟨h⁴⟩/(4π⟨h²⟩²). Both ⟨h²⟩ and ⟨h⁴⟩ are derived from integrals over the SMBHB chirp‑mass, mass‑ratio, and redshift distributions, using the standard GW energy‑spectrum for circular, GW‑driven binaries. Because ⟨h⁴⟩ is weighted by M⁵, the anisotropy is highly sensitive to the high‑mass tail of the SMBH mass function (M ≳ 10⁹ M⊙).

A key result is the steep frequency scaling of the shot‑noise power: C_SN,ℓ ∝ f⁸⁄³ for pure GW‑driven inspirals. This scaling follows from the residence‑time relation dt_r/dln f ∝ f^(-β), with β = 8/3 in the GW‑only case. If additional mechanisms (gas drag, stellar scattering) shorten the residence time at low frequencies (β < 8/3), the frequency dependence flattens, providing a direct observational probe of SMBHB dynamics.

For the LSS contribution, the merger rate density is modulated by a linear bias factor b_BH(M,q,z) and the linear matter overdensity δ_lin(Ω̂,z). After integrating over mass and mass‑ratio, the authors define a redshift‑dependent effective bias ⟨b_BH(z)⟩ and a weighting function P_{h²}(z) that captures how much each redshift slice contributes to the total strain. The resulting angular power spectrum C_LSS,ℓ is proportional to the matter power spectrum evaluated at k = ℓ/χ(z) and is suppressed by roughly two to three orders of magnitude relative to the shot‑noise term.

Concrete numerical examples employ empirically calibrated SMBHB merger‑rate models from recent literature (e.g., Refs.