$\texttt{BilbyFlow}$: user-friendly neural posterior estimation for gravitational-wave astronomy

Abstract

Bayesian inference plays a central role in the new field of gravitational-wave astronomy. However, traditional Bayesian inference with stochastic samplers is computationally expensive, taking hours to days per event. Transformative changes are therefore required to enable the science of next-generation observatories whose event rates and signal-to-noise ratios will increase significantly over the current generation. Recent work has shown that neural posterior estimation (NPE) is a promising path forward. A neural net is trained to approximate the posterior distribution of gravitational-wave parameters, allowing generation of posterior samples in a fraction of the time required by stochastic samplers. In this work, we introduce $\texttt{BilbyFlow}$, which harnesses the power of NPE in the popular $\texttt{Bilby}$ code suite. We use $\texttt{BilbyFlow}$ to analyze a subset of 38 high-mass events from the third LIGO-Virgo-KAGRA Gravitational-Wave Transient Catalog (GWTC-3). For 29 events (76\%), we obtained an importance-sampling efficiency $>$1%, allowing us to produce reliable posterior distributions within 3 min - 1.5 hours. For the other events, with importance-sampling efficiency $\ll$1%, the run time can be as long as 35 hours. We achieve a median importance-sampling efficiency of 7%, which is roughly comparable to the $\texttt{DINGO}$ package. We aim to significantly improve this efficiency with further development to make the runtime more reliably $O(\text{min})$. $\texttt{BilbyFlow}$ is open source and $\texttt{pip}$-installable.