Fast Bayesian analysis of individual binaries in pulsar timing array data

Abstract

Searching for gravitational waves in pulsar timing array data is computationally intensive. The data is unevenly sampled, and the noise is heteroscedastic, necessitating the use of a time-domain likelihood function with attendant expensive matrix operations. The computational cost is exacerbated when searching for individual supermassive black hole binaries, which have a large parameter space due to the additional pulsar distance, phase offset and noise model parameters needed for each pulsar. We introduce a new formulation of the likelihood function which can be used to make the Bayesian analysis significantly faster. We divide the parameters into projection and shape parameters. We then accelerate the exploration of the projection parameters by more than four orders of magnitude by precomputing the expensive inner products for each set of shape parameters. The projection parameters include nuisance parameters such as the gravitational wave phase offset at each pulsar. In the new scheme, these troublesome nuisance parameters are efficiently marginalized over using multiple-try Markov chain Monte Carlo sampling as part of a Metropolis-within-Gibbs scheme. The acceleration provided by our method will become increasingly important as pulsar timing datasets rapidly grow. Our method also makes sophisticated analyses more tractable, such as searches for multiple binaries, or binaries with non-negligible eccentricities.

16 pages, 11 figures, matches published version

References in corpus (15)

• The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background
• 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
• Expected properties of the first gravitational wave signal detected with pulsar timing arrays
• An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set
• Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset
• Astrophysics Milestones For Pulsar Timing Array Gravitational Wave Detection
• Multi-messenger Time-domain Signatures Of Supermassive Black Hole Binaries
• Common-spectrum process versus cross-correlation for gravitational-wave searches using pulsar timing arrays
• Accelerated Bayesian model-selection and parameter-estimation in continuous gravitational-wave searches with pulsar-timing arrays
• Multi-messenger Approaches to Supermassive Black Hole Binary Detection and Parameter Estimation: Implications for Nanohertz Gravitational Wave Searches with Pulsar Timing Arrays
• Constraining alternative polarization states of gravitational waves from individual black hole binaries using pulsar timing arrays
• Bayesian search for gravitational wave bursts in pulsar timing array data
• A Parallelized Bayesian Approach To Accelerated Gravitational-Wave Background Characterization