Efficient Bayesian inference and model selection for continuous gravitational waves in pulsar timing array data

Abstract

Finding and characterizing gravitational waves from individual supermassive black hole binaries is a central goal of pulsar timing array experiments, which will require analysis methods that can be efficient on our rapidly growing datasets. Here we present a novel approach built on three key elements: i) precalculating and interpolating expensive matrix operations; ii) semi-analytically marginalizing over the gravitational-wave phase at the pulsars; iii) numerically marginalizing over the pulsar distance uncertainties. With these improvements the recent NANOGrav 15yr dataset can be analyzed in minutes after an $\mathcal{O}(1\mathrm{hour})$ setup phase, instead of an analysis taking days-weeks with previous methods. The same setup can be used to efficiently analyze the dataset under any sinusoidal deterministic model. In particular, this will aid testing the binary hypothesis by allowing for efficient analysis of competing models (e.g. incoherent, monopolar, or dipolar sine wave model) and scrambled datasets for false alarm studies. The same setup can be updated in minutes for new realizations of the data, which enables large simulation studies.

arXiv comment19 pages, 3 figures, version matching published paper

References in corpus (51)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
• 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
• Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library
• The NANOGrav 15-year Data Set: Search for Signals from New Physics
• The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
• The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars
• The second data release from the European Pulsar Timing Array I. The dataset and timing analysis
• Gravitational waves from resolvable massive black hole binary systems and observations with Pulsar Timing Arrays
• 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
• The Parkes Pulsar Timing Array Third Data Release
• Gravitational wave astronomy of single sources with a pulsar timing array
• The Local Nanohertz Gravitational-Wave Landscape From Supermassive Black Hole Binaries
• The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• Nested sampling for physical scientists
• The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget
• An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set
• NANOGrav Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries in Circular Orbits
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• The Sensitivity of the Parkes Pulsar Timing Array to Individual Sources of Gravitational Waves
• 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
• Electromagnetic Counterparts to Massive Black Hole Mergers
• Optimal strategies for continuous gravitational wave detection in pulsar timing arrays
• All correlations must die: Assessing the significance of a stochastic gravitational-wave background in pulsar-timing arrays
• The second data release from the European Pulsar Timing Array V. Search for continuous gravitational wave signals
• Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array
• Towards Robust Gravitational Wave Detection with Pulsar Timing Arrays
• 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
• Eryn : A multi-purpose sampler for Bayesian inference
• Multi-messenger Time-domain Signatures Of Supermassive Black Hole Binaries
• Accelerated Bayesian model-selection and parameter-estimation in continuous gravitational-wave searches with pulsar-timing arrays
• Pulsar Timing Array Observations of Massive Black Hole Binaries
• Searching for gravitational waves from known pulsars using the $\mathcal{F}$ and $\mathcal{G}$ statistics
• The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
• Joint search for isolated sources and an unresolved confusion background in pulsar timing array data
• Pulsar timing array signals induced by black hole binaries in relativistic eccentric orbits
• Multi-Messenger Astrophysics with Pulsar Timing Arrays
• Constraining alternative polarization states of gravitational waves from individual black hole binaries using pulsar timing arrays
• The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B
• Exploring the Capabilities of Gibbs Sampling in Pulsar Timing Arrays
• Post-Newtonian-accurate pulsar timing array signals induced by inspiralling eccentric binaries: accuracy, computational cost, and single-pulsar search
• Implementation of an efficient Bayesian search for gravitational wave bursts with memory in pulsar timing array data
• Iterative time-domain method for resolving multiple gravitational wave sources in Pulsar Timing Array data
• Rapid parameter estimation for pulsar-timing-array datasets with variational inference and normalizing flows

Cited by in corpus (2)

• Separating deterministic and stochastic gravitational wave signals in realistic pulsar timing array datasets
• Realistic assessment of a single gravitational wave source localization taking into account precise pulsar distances with pulsar timing arrays