A Parallelized Bayesian Approach To Accelerated Gravitational-Wave Background Characterization

Abstract

The characterization of nanohertz-frequency gravitational waves (GWs) with pulsar-timing arrays requires a continual expansion of datasets and monitored pulsars. Whereas detection of the stochastic GW background is predicated on measuring a distinctive pattern of inter-pulsar correlations, characterizing the background's spectrum is driven by information encoded in the power spectra of the individual pulsars' time series. We propose a new technique for rapid Bayesian characterization of the stochastic GW background that is fully parallelized over pulsar datasets. This Factorized Likelihood (FL) technique empowers a modular approach to parameter estimation of the GW background, multi-stage model selection of a spectrally-common stochastic process and quadrupolar inter-pulsar correlations, and statistical cross-validation of measured signals between independent pulsar sub-arrays. We demonstrate the equivalence of this technique's efficacy with the full pulsar-timing array likelihood, yet at a fraction of the required time. Our technique is fast, easily implemented, and trivially allows for new data and pulsars to be combined with legacy datasets without re-analysis of the latter.

14 pages, 6 figures. Matches version accepted by PRD

References in corpus (11)

• The stochastic gravitational-wave background from massive black hole binary systems: implications for observations with Pulsar Timing Arrays
• 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
• Assessing the Role of Spin Noise in the Precision Timing of Millisecond Pulsars
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data
• Massive black hole binary systems and the NANOGrav 12.5 year results
• Low-rank approximations for large stationary covariance matrices, as used in the Bayesian and generalized-least-squares analysis of pulsar-timing data
• Common-spectrum process versus cross-correlation for gravitational-wave searches using pulsar timing arrays
• Bayesian cross validation for gravitational-wave searches in pulsar-timing array data