Rapid inference for individual binaries and a stochastic background with pulsar timing array data

Abstract

The analysis of pulsar timing array data has provided evidence for a gravitational wave background in the nanohertz band. This raises the question of what is the source of the signal, is it astrophysical or cosmological in origin? If the signal originates from a population of supermassive black hole binaries, as is generally assumed, we can expect to see evidence for both anisotropy and to be able to resolve signals from individual binaries as more data are collected. The anisotropy and resolvable systems are caused by a small number of loud signals that stand out from the crowd. Here we focus on the joint detection of individual signals and a stochastic background. While methods have previously been developed to perform such an analysis, they are currently held back by the cost of computing the joint likelihood function. Each individual source is described by $N=8+2N_p$ parameters, where $N_p$ are the number of pulsars in the array. With the latest combined datasets having over one hundred pulsars, the parameter space is very large, and consequently, it takes a large number of likelihood evaluations to explore these models. Here we present a new approach that extends the Fourier basis method, previously introduced to accelerate analyses for stochastic signals, to also include deterministic signals. Key elements of the method are that the likelihood evaluations are per-pulsar, avoiding expensive operations on large matrices, and the templates for individual binaries can be computed analytically or using fast Fourier methods on a sparsely sampled grid of time samples. This analysis method scales better than quadratically with the size of the dataset, while the approach currently being used in most analyses scales quartically or worse with the number of data points. As datasets grow with more observations, this analysis will be orders of magnitude faster than previous approaches.

19 pages, 11 figures

References in corpus (29)

• The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background
• Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array
• Searching for the nano-Hertz stochastic gravitational wave background with the Chinese Pulsar Timing Array Data Release I
• The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals
• The NANOGrav 15-year Data Set: Search for Signals from New Physics
• Detection methods for stochastic gravitational-wave backgrounds: a unified treatment
• The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars
• A General Framework for the Parametrization of Hierarchical Models
• Assessing the Role of Spin Noise in the Precision Timing of Millisecond Pulsars
• On measuring the gravitational-wave background using Pulsar Timing Arrays
• The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• Understanding and analysing time-correlated stochastic signals in pulsar timing
• The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Resolving multiple supermassive black hole binaries with pulsar timing arrays
• A Bayesian analysis pipeline for continuous GW sources in the PTA band
• Low-rank approximations for large stationary covariance matrices, as used in the Bayesian and generalized-least-squares analysis of pulsar-timing data
• The NANOGrav 15-year Gravitational-Wave Background Methods
• Joint search for isolated sources and an unresolved confusion background in pulsar timing array data
• Targeting supermassive black hole binaries and gravitational wave sources for the pulsar timing array
• Towards a unified treatment of gravitational-wave data analysis
• Accelerating pulsar timing data analysis
• Bayesian search for gravitational wave bursts in pulsar timing array data
• Optimal reconstruction of the Hellings and Downs correlation
• Exploring the Capabilities of Gibbs Sampling in Pulsar Timing Arrays
• Efficient Gravitational Wave Searches with Pulsar Timing Arrays using Hamiltonian Monte Carlo
• Solving the PTA Data Analysis Problem with a Global Gibbs Scheme
• Regularizing the Pulsar Timing Array likelihood: A path towards Fourier Space