Quality over Quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals

Abstract

The search for gravitational waves using Pulsar Timing Arrays (PTAs) is a computationally expensive complex analysis that involves source-specific noise studies. As more pulsars are added to the arrays, this stage of PTA analysis will become increasingly challenging. Therefore, optimizing the number of included pulsars is crucial to reduce the computational burden of data analysis. Here, we present a suite of methods to rank pulsars for use within the scope of PTA analysis. First, we use the maximization of the signal-to-noise ratio as a proxy to select pulsars. With this method, we target the detection of stochastic and continuous gravitational wave signals. Next, we present a ranking that minimizes the coupling between spatial correlation signatures, namely monopolar, dipolar, and Hellings & Downs correlations. Finally, we also explore how to combine these two methods. We test these approaches against mock data using frequentist and Bayesian hypothesis testing. For equal-noise pulsars, we find that an optimal selection leads to an increase in the log-Bayes factor two times steeper than a random selection for the hypothesis test of a gravitational wave background versus a common uncorrelated red noise process. For the same test but for a realistic EPTA dataset, a subset of 25 pulsars selected out of 40 can provide a log-likelihood ratio that is $89\%$ of the total, implying that an optimally selected subset of pulsars can yield results comparable to those obtained from the whole array. We expect these selection methods to play a crucial role in future PTA data combinations.

References in corpus (15)

• Array Programming with NumPy
• TEMPO2, a new pulsar timing package. I: Overview
• Tempo2, a new pulsar timing package. II: The timing model and precision estimates
• The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
• On the evidence for a common-spectrum process in the search for the nanohertz gravitational-wave background with the Parkes Pulsar Timing Array
• Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search
• Assessing the Role of Spin Noise in the Precision Timing of Millisecond Pulsars
• Expected properties of the first gravitational wave signal detected with pulsar timing arrays
• On measuring the gravitational-wave background using Pulsar Timing Arrays
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data
• A pulsar-based timescale from the International Pulsar Timing Array
• Massive black hole binary systems and the NANOGrav 12.5 year results
• Common-spectrum process versus cross-correlation for gravitational-wave searches using pulsar timing arrays
• Ephemeris Errors and the Gravitational Wave Signal: Harmonic Mode Coupling in Pulsar Timing Array Searches