A Trifecta of Modelling Tools: Navigating COMPAS, Shark and Bayesian Inference for Binary Black Hole Model Selection
Authors:
Liana Rauf,
Cullan Howlett,
Simon Stevenson,
Jeff Riley
Abstract:
Gravitational waves (GWs) have revealed surprising properties of binary black hole (BBH) populations, but there is still mystery surrounding how these compact objects evolve. We apply Bayesian inference and an efficient method to calculate the BBH merger rates in the Shark host galaxies, to determine the combination of COMPAS parameters that outputs a population most like the GW sources from the L…
▽ More
Gravitational waves (GWs) have revealed surprising properties of binary black hole (BBH) populations, but there is still mystery surrounding how these compact objects evolve. We apply Bayesian inference and an efficient method to calculate the BBH merger rates in the Shark host galaxies, to determine the combination of COMPAS parameters that outputs a population most like the GW sources from the LVK transient catalogue. For our COMPAS models, we calculate the likelihood with and without the dependence on the predicted number of BBH merger events. We find strong correlations between hyper-parameters governing the specific angular momentum (AM) of mass lost during mass transfer, the mass-loss rates of Wolf--Rayet stars via winds and the chemically homogeneous evolution (CHE) formation channel. We conclude that analysing the marginalised and unmarginalised likelihood is a good indicator of whether the population parameters distribution and number of observed events reflect the LVK data. In doing so, we see that the majority of the models preferred in terms of the population-level parameters of the BBHs greatly overpredict the number of events we should have observed to date. Looking at the smaller number of models which perform well with both likelihoods, we find that those with no CHE, AM loss occurring closer to the donor during the first mass-transfer event, and/or higher rates of mass-loss from Wolf--Rayet winds are generally preferred by current data. We find these conclusions to be robust to our choice of selection criteria.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
Exploring binary black hole mergers and host galaxies with Shark and COMPAS
Authors:
Liana Rauf,
Cullan Howlett,
Tamara M. Davis,
Claudia D. P. Lagos
Abstract:
We explore the connection between the gravitational wave (GW) merger rates of stellar-mass binary black holes (BBH) and galaxy properties. We do this by generating populations of stars using the binary population synthesis code COMPAS and evolving them in galaxies from the semi-analytic galaxy formation model Shark, to determine the number of mergers occurring in each simulation time-step. We find…
▽ More
We explore the connection between the gravitational wave (GW) merger rates of stellar-mass binary black holes (BBH) and galaxy properties. We do this by generating populations of stars using the binary population synthesis code COMPAS and evolving them in galaxies from the semi-analytic galaxy formation model Shark, to determine the number of mergers occurring in each simulation time-step. We find that metal-rich and massive galaxies with star formation rate (SFR) greater than $1M_{\odot}/ \rm yr$ are 10 times more likely to have GW events compared to younger, less massive and metal poor galaxies. Our simulation with the default input parameters predicts a higher local merger rate density compared to the third GW transient catalogue (GWTC-3) prediction from LIGO, VIRGO and KAGRA, due to short coalescence times, low metallicities and a high SFR at low redshift in the simulation, which produces more BBHs that merge within the age of the Universe compared to observations. We identify alternate remnant mass models that more accurately reproduce the volumetric rate and provide updated fits to the merger rate as a function of redshift. We then investigate the relative fraction of GW events in our simulation that are in observable host galaxies from upcoming galaxy surveys, determining which of those are ideal for tracing host galaxies with high merger rates. The implications of this work can be utilised for constraining stellar evolution models, better informing follow-up programs, and placing informative priors on host galaxies when measuring cosmological parameters such as the Hubble constant.
△ Less
Submitted 12 June, 2023; v1 submitted 16 February, 2023;
originally announced February 2023.
