# Addressing leakage and mode suppression in angular power spectrum estimation for gravitational-wave backgrounds using pulsar timing arrays Preprint · Deepali Agarwal, Joseph D. Romano, Yacine Ali-Haïmoud, Tristan L. Smith - Archive id: arxiv:2602.20075 - Published: 2026-02-23 - Venue: Phys. Rev. D 114, 063048 (2026) - DOI: https://doi.org/10.1103/cm48-4k18 - arXiv: https://arxiv.org/abs/2602.20075 - Source URL: https://arxiv.org/abs/2602.20075 - PDF: https://arxiv.org/pdf/2602.20075 - Provenance: arxiv ## Abstract Mapping gravitational-wave background (GWB) anisotropy with pulsar timing arrays (PTAs) is affected by harmonic-space mode suppression and mode coupling arising from an array's nonuniform sky response. Spherical harmonic expansions must be truncated at finite multipole l\_max^rec, often set to l\_max^N\_pair$\\equiv {\\rm int}\\left\[\\sqrt{\\text{N\_pair}}-1\\right\]$, where N\_pair is the number of distinct pulsar pairs in an array. This choice is motivated by the counting argument that cross-correlations provide at most N\_pair independent constraints. We obtain the multipole l\_max^res corresponding to the maximum informative angular scale of a PTA. It is defined such that expansions to l\_max^res (approximately) span the space of "observable skies" encoded in the N\_pair eigenmaps of the Fisher information matrix, and therefore depends on the array configuration. We explicitly show that GWB power contained in multipoles l$\\gtrsim$l\_max^res do not significantly affect analyses that use expansions out to l\_max^res, because the PTA response acts as a low-pass filter. In contrast, truncating at l\_max^rec\< l\_max^res leads to leakage of small-scale angular power from l\_max^rec\