Pulsar timing array analysis in a Legendre polynomial basis

We use Legendre polynomials, previously employed in this context by Lee et al., van Haasteren and Levin, and Pitrou and Cusin, to model signals in pulsar timing arrays. These replace the (Fourier mode) basis of trigonometric functions normally used for data analysis. The Legendre basis makes it simpler to incorporate pulsar modeling effects, which remove constant-, linear-, and quadratic-in-time terms from pulsar timing residuals. In the Legendre basis, this zeroes the amplitudes of the the first three Legendre polynomials. We use this basis to construct an optimal quadratic cross-correlation estimator $\widehat{\mu}$ of the Hellings and Downs (HD) correlation and compute its variance $\sigma^2_{\widehat{\mu}}$ in the way described by Allen and Romano. Remarkably, if the gravitational-wave background (GWB) and pulsar noise power spectra are (sums of) power laws in frequency, then in this basis one obtains analytic closed forms for many quantities of interest.