# Cosmic Variance in Anisotropy Searches at Pulsar Timing Arrays

## TL;DR

This paper investigates how cosmic variance affects anisotropy measurements in pulsar timing array data, demonstrating that statistical inference remains robust and can distinguish between astrophysical and cosmological gravitational wave backgrounds.

## Contribution

It provides a detailed analysis of cosmic variance impacts on PTA anisotropy searches and confirms the robustness of statistical inference in this context.

## Key findings

- Anisotropy measurements are robust despite cosmic variance.
- Maximum multipole scales linearly with the number of pulsars.
- PTA data can effectively discriminate between astrophysical and cosmological GWBs.

## Abstract

Recent pulsar timing array (PTA) analyses show evidence for a gravitational wave background (GWB) with angular correlations consistent with the Hellings-Downs curve. Anisotropies are a key discriminator of the origin of this GWB, as they are expected to be at 1--20\% for astrophysical sources, but suppressed for cosmological GWBs. However, contrary to gravitational wave detectors at higher frequencies, PTAs only take a few independent measurements of a GWB and consequently are highly sensitive to cosmic variance, which induces apparent anisotropies in individual realizations of an isotropic GWB. We demonstrate explicitly that statistical inference nevertheless remains robust, i.e., measurements are consistent with the underlying assumption of isotropy. This confirms that searches for anisotropies will be able to robustly discriminate astrophysical from cosmological GWBs. En route, we demonstrate that the maximum multipole constrained by a PTA dataset scales linearly with the number of pulsars $\ell_{\rm max} \sim N_p$.