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.

19 pages, 3 figures, version matching published paper

References in corpus (14)

• Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library
• The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
• 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
• The Local Nanohertz Gravitational-Wave Landscape From Supermassive Black Hole Binaries
• Nested sampling for physical scientists
• 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
• Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array
• 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
• Multi-Messenger Astrophysics with Pulsar Timing Arrays
• Constraining alternative polarization states of gravitational waves from individual black hole binaries using pulsar timing arrays
• 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