Analytic series expansion of the overlap reduction function for gravitational wave search with pulsar timing arrays

In our previous paper A. Boîtier et al. [Phys. Rev. D 102, 064051 (2020), 10.1103/PhysRevD.102.064051], we derived a generic expression for the pulse redshift, the main observable for the Pulsar Timing Array (PTA) experiment for detection of gravitational waves for all possible polarizations induced by modifications of general relativity. In this work, we provide a generic expression of the overlap reduction function for PTA without using the short wavelength approximation for tensorial polarization. We find that when the overlap reduction function does not have the exponential terms, which is the case when using short wavelength approximation, it leads to discontinuities and poles. In this work, we provide a series expansion to calculate the integral exactly and investigate the behavior of the series for short wavelength values via numerical evaluation of the analytical series. We find a disagreement for the limit of colocated pulsars with the Hellings and Downs curve. Otherwise, our formalism agrees with the Hellings and Downs curve for a broad range of pulsars.