The MeerKAT Pulsar Timing Array: The first search for gravitational waves with the MeerKAT radio telescope

TL;DR

This paper reports the first gravitational wave search using the MeerKAT Pulsar Timing Array, analyzing 4.5 years of data from 83 pulsars, and finds potential correlations consistent with a stochastic gravitational wave background.

Contribution

It presents the initial results of gravitational wave background searches with MeerKAT, demonstrating sensitivity and the importance of noise modeling in detecting spatial correlations.

Findings

Potential Hellings-Downs correlations at 3-3.4σ significance

Measured characteristic strain amplitude around 7.5×10^{-15} at spectral index -0.26

Results depend on noise assumptions and pulsar proximity

Abstract

Pulsar Timing Arrays search for nanohertz-frequency gravitational waves by regularly observing ensembles of millisecond pulsars over many years to look for correlated timing residuals. Recently the first evidence for a stochastic gravitational wave background has been presented by the major Arrays, with varying levels of significance ($\sim$2-4$\sigma$). In this paper we present the results of background searches with the MeerKAT Pulsar Timing Array. Although of limited duration (4.5 yr), the $\sim$ 250,000 arrival times with a median error of just $3 \mu$s on 83 pulsars make it very sensitive to spatial correlations. Detection of a gravitational wave background requires careful modelling of noise processes to ensure that any correlations represent a fit to the underlying background and not other misspecified processes. Under different assumptions about noise processes we can produce either what appear to be compelling Hellings-Downs correlations of high significance (3-3.4$\sigma$) with a spectrum close to that which is predicted, or surprisingly, under slightly different assumptions, ones that are insignificant. This appears to be related to the fact that many of the highest precision MeerKAT Pulsar Timing Array pulsars are in close proximity and dominate the detection statistics. The sky-averaged characteristic strain amplitude of the correlated signal in our most significant model is $h_{c, {\rm yr}} = 7.5^{+0.8}_{-0.9} \times 10^{-15}$ measured at a spectral index of $\alpha=-0.26$, decreasing to $h_{c, {\rm yr}} = 4.8^{+0.8}_{-0.9} \times 10^{-15}$ when assessed at the predicted $\alpha=-2/3$. These data will be valuable as the International Pulsar Timing Array project explores the significance of gravitational wave detections and their dependence on the assumed noise models.