A neural-network surrogate for the log-likelihood gradients makes Hamiltonian Monte Carlo trajectories 30 times faster than relative-binning gradients and recovers LVK-consistent posteriors for two binary neutron star mergers.
A Hamiltonian Monte Carlo method for Bayesian Inference of Supermassive Black Hole Binaries
1 Pith paper cite this work. Polarity classification is still indexing.
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.
fields
gr-qc 1years
2025 1verdicts
CONDITIONAL 1representative citing papers
citing papers explorer
-
DeepHMC : a deep-neural-network acclerated Hamiltonian Monte Carlo algorithm for binary neutron star parameter estimation
A neural-network surrogate for the log-likelihood gradients makes Hamiltonian Monte Carlo trajectories 30 times faster than relative-binning gradients and recovers LVK-consistent posteriors for two binary neutron star mergers.