Accurate pulsar timing array residual variances and correlation of the stochastic gravitational wave background

Abstract: Pulsar timing arrays have reported a compelling evidence of a nanohertz stochastic gravitational wave background. However, the origin of the signal remains undetermined, largely because its spectrum is bluer for an astrophysical source and can be explained by cosmological models. In this letter, we derive theoretically accurate expressions for the Fourier bin variances and correlation of pulsar timing residuals, and demonstrate their outstanding agreement with astrophysical simulations. In contrast, we show that a common power law traditionally used to interpret a stochastic gravitational wave background spectrum is generally faced with systematic errors, one of which is the illusion of a bluer signal. This hints at a conservative solution, supportive of an astrophysical source, to the observed correlated common spectrum process in pulsar timing arrays.