Bayesian inference for binary neutron star inspirals using a Hamiltonian Monte Carlo Algorithm

Abstract

The coalescence of binary neutron stars are one of the main sources of gravitational waves for ground-based gravitational wave detectors. As Bayesian inference for binary neutron stars is computationally expensive, more efficient and faster converging algorithms are always needed. In this work, we conduct a feasibility study using a Hamiltonian Monte Carlo algorithm (HMC). The HMC is a sampling algorithm that takes advantage of gradient information from the geometry of the parameter space to efficiently sample from the posterior distribution, allowing the algorithm to avoid the random-walk behaviour commonly associated with stochastic samplers. As well as tuning the algorithm's free parameters specifically for gravitational wave astronomy, we introduce a method for approximating the gradients of the log-likelihood that reduces the runtime for a $10^6$ trajectory run from ten weeks, using numerical derivatives along the Hamiltonian trajectories, to one day, in the case of non-spinning neutron stars. Testing our algorithm against a set of neutron star binaries using a detector network composed of Advanced LIGO and Advanced Virgo at optimal design, we demonstrate that not only is our algorithm more efficient than a standard sampler, but a $10^6$ trajectory HMC produces an effective sample size on the order of $10^4-10^5$ statistically independent samples.

16 pages, 8 figures. Submitted to PRD

References in corpus (15)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Advanced Virgo: a 2nd generation interferometric gravitational wave detector
• Multi-messenger Observations of a Binary Neutron Star Merger
• GW170817: Measurements of Neutron Star Radii and Equation of State
• GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence
• Properties of the binary neutron star merger GW170817
• GW170608: Observation of a 19-solar-mass Binary Black Hole Coalescence
• Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library
• Comparison of post-Newtonian templates for compact binary inspiral signals in gravitational-wave detectors
• The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo
• Parameter estimation for binary neutron-star coalescences with realistic noise during the Advanced LIGO era
• Parameter estimation of spinning binary inspirals using Markov-chain Monte Carlo
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• Catching Super Massive Black Hole Binaries Without a Net
• Fisher vs. Bayes : A comparison of parameter estimation techniques for massive black hole binaries to high redshifts with eLISA

Cited by in corpus (4)

• The Science of the Einstein Telescope
• Gravitational-wave surrogate models powered by artificial neural networks: The ANN-Sur for waveform generation
• Binary Neutron Stars Gravitational Wave Detection Based on Wavelet Packet Analysis And Convolutional Neural Networks
• Bayesian inference in single-line spectroscopic binaries with a visual orbit