Scalable continuous gravitational wave detection in PTA data with non-parametric red noise suppression and optimal pulsar selection

TL;DR

This paper presents a fast, scalable frequentist method for detecting continuous gravitational waves in pulsar timing array data, reducing computational costs while maintaining accuracy comparable to Bayesian approaches.

Contribution

It introduces a non-parametric red noise suppression technique combined with optimal pulsar selection, enabling efficient analysis of large-scale PTAs.

Findings

Achieves similar detection accuracy to Bayesian methods with significantly less computation.

Completes analysis in less than 5 hours compared to 1-2 days for Bayesian methods.

Effective on simulated data based on NANOGrav 15-year dataset.

Abstract

Bayesian methods for the detection of continuous gravitational waves (CGWs) in Pulsar Timing Array (PTA) data incur substantial computational costs that grow rapidly due to the number of noise and signal parameters characterizing the fitted model being proportional to the size of the PTA. This computational burden limits the scalability of these methods for large-scale PTAs comprising hundreds of pulsars anticipated from next-generation radio astronomy facilities. In this work, we introduce a computationally efficient frequentist method designed to circumvent this challenge. This is achieved by combining an adaptive spline fitting algorithm that non-parametrically suppresses red noise, thereby eliminating the need for complex noise modeling inherent to Bayesian methods, with a novel scheme for optimizing the subsets of pulsars included in the search. We quantify the performance of our method on a simulated dataset based on the NANOGrav 15-year data release and find that it achieves a performance comparable to that of Bayesian analysis: for a CGW signal with a signal-to-noise ratio of $\approx 10$, our method yields a relative characteristic strain error of 1.0\% and a frequency error of 0.072\% from the injected values by using the optimal pulsar selections, while the same errors are 1.7\% and 0.16\%, respectively, for the standard Bayesian analysis. At the same time, our analysis completes in less than 5 hours, in contrast to the 1-2 days required by Bayesian methods. This allows us to perform a rigorous study of our method using multiple data realizations and signal parameters, establishing it as an efficient and scalable tool for CGW searches with large-scale PTAs.