The MeerKAT Pulsar Timing Array: Maps of the gravitational-wave sky with the 4.5 year data release

Abstract: In an accompanying publication, the MeerKAT Pulsar Timing Array (MPTA) collaboration reports tentative evidence for the presence of a stochastic gravitational-wave background, following observations of similar signals from the European and Indian Pulsar Timing Arrays, NANOGrav, the Parkes Pulsar Timing Array and the Chinese Pulsar Timing Array. If such a gravitational-wave background signal originates from a population of inspiraling supermassive black-hole binaries, the signal may be anisotropically distributed on the sky. In this Letter we evaluate the anisotropy of the MPTA signal using a spherical harmonic decomposition. We discuss complications arising from the covariance between pulsar pairs and regularisation of the Fisher matrix. Applying our method to the 4.5 yr dataset, we obtain two forms of sky maps for the three most sensitive MPTA frequency bins between 7 -21 nHz. Our "clean maps'' estimate the distribution of gravitational-wave strain power with minimal assumptions. Our radiometer maps answer the question: is there a statistically significant point source? We find a noteworthy hotspot in the 7 nHz clean map with a $p$-factor of $p=0.015$ (not including trial factors). Future observations are required to determine if this hotspot is of astrophysical origin.

• The paper presents an extensive analysis of the gravitational-wave background using timing data from 83 millisecond pulsars observed over 4.5 years.
• It employs spherical harmonic decomposition and overlap reduction functions to map anisotropy across multiple nanohertz frequency bins.
• The detection of a potential hotspot in the 7 nanohertz map hints at anisotropic distributions, emphasizing the need for further observational studies.

Anisotropy in the Gravitational-Wave Background: Insights from the MeerKAT Pulsar Timing Array

The paper "The MeerKAT Pulsar Timing Array: Maps of the gravitational-wave sky with the 4.5 year data release" presents a detailed analysis of the gravitational-wave background using data obtained from the MeerKAT Pulsar Timing Array (MPTA). This research is an extension of the efforts of the International Pulsar Timing Array (IPTA) and associated entities, aimed at understanding the characteristics and potential anisotropy of gravitational waves (GWs) emanating from cosmological and astrophysical sources, particularly supermassive black-hole binaries.

The authors utilize timing data from 83 millisecond pulsars observed by the MeerKAT radio telescope over four and a half years. The MPTA is instrumental in probing nanohertz frequency ranges, which are expected to be populated by GWs from inspiraling supermassive black-hole binaries. Identifying GW-induced anisotropy could provide critical insights into the distribution and nature of these binaries, as well as other GW sources such as cosmic strings or primordial processes in the early universe.

Methodology

The investigation employs a spherical harmonic decomposition of the GW signal to evaluate anisotropy. This is operationalized through the computation of overlap reduction functions for pulsar pairs, which are sensitive to angular patterns in the GW background. The methodology is informed by established frameworks employed in both the LIGO community and PTA analyses, albeit with adjustments to account for PTAs' unique sensitivity.

The response of each pulsar to GWs is modeled, and covariance matrices are constructed to understand the correlations between pulsars in response to the GW background. The core analysis involves deriving clean maps and radiometer maps using statistical approaches optimized for both broad-scale anisotropy and point source detection, respectively. These maps are regularized and analyzed to identify significant deviations from isotropy, which would suggest anisotropic distribution of GW sources across the sky.

Results

The analysis results in sky maps of GW power for three frequency bins: 7, 14, and 21 nanohertz. The study finds a potential hotspot in the 7 nanohertz map with a low p-value of 0.015, hinting at possible anisotropy or a point source. However, the authors maintain a cautious interpretation of these results due to the statistical insignificance and suggest that future observations may clarify the nature of this anomaly.

Implications and Future Work

This research contributes to the growing body of work aimed at mapping the GW background and understanding its sources. The ability to detect anisotropy in the GW background could significantly advance our understanding of supermassive black hole mergers and other processes that might contribute to the GW background. The detection of hotspots could eventually tie these GW signals to specific astrophysical events or structures in the universe.

The authors acknowledge the need for further analysis, particularly concerning cosmic variance and the self-consistent treatment of GW noise in the analysis framework. The ongoing enhancement of PTAs, both in terms of the number of pulsars observed and the precision of the measurements, will be crucial for improving the sensitivity and resolution of these sky maps. As the field progresses, integrating data from multiple PTAs and refining analytical methods will be essential for deciphering the complexities of the gravitational-wave sky.