bajes: Bayesian inference of multimessenger astrophysical data, methods and application to gravitational-waves
Abstract
We present bajes, a parallel and lightweight framework for Bayesian inference of multimessenger transients. bajes is a Python modular package with minimal dependencies on external libraries adaptable to the majority of the Bayesian models and to various sampling methods. We describe the general workflow and the parameter estimation pipeline for compact-binary-coalescence gravitational-wave transients. The latter is validated against injections of binary black hole and binary neutron star waveforms, including confidence interval tests that demonstrates the inference is well-calibrated. Binary neutron star postmerger injections are also studied using a network of five detectors made of LIGO, Virgo, KAGRA and Einstein Telescope. Postmerger signals will be detectable for sources at 80\,Mpc, with Einstein Telescope contributing over 90\% of the total signal-to-noise ratio. As a full scale application, we re-analyze the GWTC-1 black hole transients using the effective-one-body TEOBResumS approximant, and reproduce selected results with other approximants. bajes inferences are consistent with previous results; the direct comparison of bajes and bilby analyses of GW150914 shows a maximum Jensen-Shannon divergence of 5.2×10-4. GW170817 is re-analyzed using TaylorF2 with 5.5PN point-mass and 7.5PN tides, TEOBResumSPA, and IMRPhenomPv2\NRTidal with different cutoff-frequencies of 1024\,Hz and 2048\,Hz. We find that the former choice minimizes systematics on the reduced tidal parameter, while a larger amount of tidal information is gained with the latter choice. bajes can perform these analyses in about 1~day using 128 CPUs.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.