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.

arXiv comment16 pages, 11 figures, matches published version

References in corpus (40)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background
• Data analysis of gravitational-wave signals from spinning neutron stars. I. The signal and its detection
• The NANOGrav 11-year Data Set: High-precision timing of 45 Millisecond Pulsars
• A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals
• 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
• European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries
• Gravitational wave astronomy of single sources with a pulsar timing array
• The generalized F-statistic: multiple detectors and multiple GW pulsars
• The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• The NANOGrav 12.5 yr Data Set: Observations and Narrowband Timing of 47 Millisecond Pulsars
• An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset
• Astrophysics Milestones For Pulsar Timing Array Gravitational Wave Detection
• Resolving multiple supermassive black hole binaries with pulsar timing arrays
• Single Sources in the Low-Frequency Gravitational Wave Sky: properties and time to detection by pulsar timing arrays
• Optimal strategies for continuous gravitational wave detection in pulsar timing arrays
• A Review of Multiple Try MCMC algorithms for Signal Processing
• All correlations must die: Assessing the significance of a stochastic gravitational-wave background in pulsar-timing arrays
• Towards Robust Gravitational Wave Detection with Pulsar Timing Arrays
• The Search for Massive Black Hole Binaries with LISA
• Detecting eccentric supermassive black hole binaries with pulsar timing arrays: Resolvable source strategies
• A Bayesian analysis pipeline for continuous GW sources in the PTA band
• 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
• Interacting Multiple Try Algorithms with Different Proposal Distributions
• Joint search for isolated sources and an unresolved confusion background in pulsar timing array data
• Umbrella sampling: a powerful method to sample tails of distributions
• Pulsar timing array signals induced by black hole binaries in relativistic eccentric orbits
• Multi-messenger Approaches to Supermassive Black Hole Binary Detection and Parameter Estimation: Implications for Nanohertz Gravitational Wave Searches with Pulsar Timing Arrays
• Associating Host Galaxy Candidates to Massive Black Hole Binaries resolved by Pulsar Timing Arrays
• Constraining alternative polarization states of gravitational waves from individual black hole binaries using pulsar timing arrays
• Towards a unified treatment of gravitational-wave data analysis
• Bayesian search for gravitational wave bursts in pulsar timing array data
• A Parallelized Bayesian Approach To Accelerated Gravitational-Wave Background Characterization