Normalizing Flows for Hierarchical Bayesian Analysis: A Gravitational Wave Population Study

Abstract

We propose parameterizing the population distribution of the gravitational wave population modeling framework (Hierarchical Bayesian Analysis) with a normalizing flow. We first demonstrate the merit of this method on illustrative experiments and then analyze four parameters of the latest LIGO/Virgo data release: primary mass, secondary mass, redshift, and effective spin. Our results show that despite the small and notoriously noisy dataset, the posterior predictive distributions (assuming a prior over the parameters of the flow) of the observed gravitational wave population recover structure that agrees with robust previous phenomenological modeling results while being less susceptible to biases introduced by less flexible models. Therefore, the method forms a promising flexible, reliable replacement for population inference distributions, even when data is highly noisy.

Cited by in corpus (4)

• Gravitational wave populations and cosmology with neural posterior estimation
• Is there an excess of black holes around $20M_{\odot }$? Optimising the complexity of population models with the use of reversible jump MCMC
• Applications of machine learning in gravitational wave research with current interferometric detectors
• Towards a robust and reliable deep learning approach for detection of compact binary mergers in gravitational wave data