Efficient Gravitational Wave Searches with Pulsar Timing Arrays using Hamiltonian Monte Carlo

Abstract

Pulsar timing arrays (PTAs) detect low-frequency gravitational waves (GWs) by looking for correlated deviations in pulse arrival times. Current Bayesian searches use Markov Chain Monte Carlo (MCMC) methods, which struggle to sample the large number of parameters needed to model the PTA and GW signals. As the data span and number of pulsars increase, this problem will only worsen. An alternative Monte Carlo sampling method, Hamiltonian Monte Carlo (HMC), utilizes Hamiltonian dynamics to produce sample proposals informed by first-order gradients of the model likelihood. This in turn allows it to converge faster to high dimensional distributions. We implement HMC as an alternative sampling method in our search for an isotropic stochastic GW background, and show that this method produces equivalent statistical results to similar analyses run with standard MCMC techniques, while requiring 100-200 times fewer samples. We show that the speed of HMC sample generation scales as $\mathcal{O}(N_{\mathrm{psr}}^{5/4})$ where $N_{\mathrm{psr}}$ is the number of pulsars, compared to $\mathcal{O}(N_{\mathrm{psr}}^2)$ for MCMC methods. These factors offset the increased time required to generate a sample using HMC, demonstrating the value of adopting HMC techniques for PTAs.

arXiv comment9 pages, 5 figures. Matches version accepted in Physical Review D

References in corpus (32)

• MCMC using Hamiltonian dynamics
• The NANOGrav 11-year Data Set: High-precision timing of 45 Millisecond Pulsars
• The dynamical evolution of massive black hole binaries - I. Hardening in a fixed stellar background
• High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array
• The NANOGrav 11-year Data Set: Pulsar-timing Constraints On The Stochastic Gravitational-wave Background
• The number of cosmic string loops
• The Population of Viscosity- and Gravitational Wave-Driven Supermassive Black Hole Binaries Among Luminous AGN
• The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background
• Gravitational wave stochastic background from cosmic (super)strings
• Low-frequency gravitational radiation from coalescing massive black hole binaries in hierarchical cosmologies
• A General Framework for the Parametrization of Hierarchical Models
• The International Pulsar Timing Array: Second data release
• The Astrophysics of Nanohertz Gravitational Waves
• Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data
• Expected properties of the first gravitational wave signal detected with pulsar timing arrays
• Systematic investigation of the expected gravitational wave signal from supermassive black hole binaries in the pulsar timing band
• Interaction of massive black hole binaries with their stellar environment: I. Ejection of hypervelocity stars
• Gravitational-wave cosmology across 29 decades in frequency
• On measuring the gravitational-wave background using Pulsar Timing Arrays
• The gravitational wave signal from massive black hole binaries and its contribution to the LISA data stream
• Gas driven massive black hole binaries: signatures in the nHz gravitational wave background
• Understanding and analysing time-correlated stochastic signals in pulsar timing
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• The Parkes Pulsar Timing Array Project: Second data release
• New advances in the Gaussian-process approach to pulsar-timing data analysis
• Gravitational-wave memory and pulsar timing arrays
• Prospects for gravitational-wave detection and supermassive black hole astrophysics with pulsar timing arrays
• The Indian Pulsar Timing Array: First data release
• Constraining the Solution to the Last Parsec Problem with Pulsar Timing
• A Bayesian analysis pipeline for continuous GW sources in the PTA band
• An Efficient Approximation to the Likelihood for Gravitational Wave Stochastic Background Detection Using Pulsar Timing Data
• Taming outliers in pulsar-timing datasets with hierarchical likelihoods and Hamiltonian sampling

Cited by in corpus (9)

• The NANOGrav 15-year Gravitational-Wave Background Methods
• Accurate characterization of the stochastic gravitational-wave background with pulsar timing arrays by likelihood reweighting
• 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
• Spatial and Spectral Characterization of the Gravitational-wave Background with the PTA Optimal Statistic
• Exploring the Capabilities of Gibbs Sampling in Pulsar Timing Arrays
• Solving the PTA Data Analysis Problem with a Global Gibbs Scheme
• Bayesian pulsar timing and noise analysis with Vela.jl: an overview
• Realistic assessment of a single gravitational wave source localization taking into account precise pulsar distances with pulsar timing arrays
• Rapid inference for individual binaries and a stochastic background with pulsar timing array data