A Hamiltonian Monte Carlo method for Bayesian Inference of Supermassive Black Hole Binaries

Abstract

We investigate the use of a Hamiltonian Monte Carlo to map out the posterior density function for supermassive black hole binaries. While previous Markov Chain Monte Carlo (MCMC) methods, such as Metropolis-Hastings MCMC, have been successfully employed for a number of different gravitational wave sources, these methods are essentially random walk algorithms. The Hamiltonian Monte Carlo treats the inverse likelihood surface as a "gravitational potential" and by introducing canonical positions and momenta, dynamically evolves the Markov chain by solving Hamilton's equations of motion. We present an implementation of the Hamiltonian Markov Chain that is faster, and more efficient by a factor of approximately the dimension of the parameter space, than the standard MCMC.

16 pages, 8 figures

References in corpus (15)

• Bayesian model selection for testing the no-hair theorem with black hole ringdowns
• Tests of Bayesian Model Selection Techniques for Gravitational Wave Astronomy
• The Gravitational Universe
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• A Solution to the Galactic Foreground Problem for LISA
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• An algorithm for detection of extreme mass ratio inspirals in LISA data
• MCMC Exploration of Supermassive Black Hole Binary Inspirals
• Catching Super Massive Black Hole Binaries Without a Net
• A Three-Stage Search for Supermassive Black Hole Binaries in LISA Data
• The search for spinning black hole binaries in mock LISA data using a genetic algorithm
• Bayesian inference on EMRI signals using low frequency approximations
• Inference on white dwarf binary systems using the first round Mock LISA Data Challenges data sets
• Eccentric Massive Black Hole Binaries in LISA I : The Detection Capabilities of Circular Templates
• A Bayesian approach to the study of white dwarf binaries in LISA data: The application of a reversible jump Markov chain Monte Carlo method

Cited by in corpus (11)

• Exploring the Bayesian parameter estimation of binary black holes with LISA
• Gravitational-wave surrogate models powered by artificial neural networks: The ANN-Sur for waveform generation
• GPU-accelerated massive black hole binary parameter estimation with LISA
• Machine Learning Gravitational Waves from Binary Black Hole Mergers
• Fisher vs. Bayes : A comparison of parameter estimation techniques for massive black hole binaries to high redshifts with eLISA
• Bayesian inference for binary neutron star inspirals using a Hamiltonian Monte Carlo Algorithm
• Accelerated parameter estimation of supermassive black hole binaries in LISA using a meshfree approximation
• User-specified random sampling of quantum channels and its applications
• A fully-automated end-to-end pipeline for massive black hole binary signal extraction from LISA data
• Source Confusion of Massive Black Hole Binaries for the Taiji Mission
• The Sequential Monte Carlo goes NUTS: Boosting Gravitational-Wave Inference