2023MNRAS.520.5259A

The current interpretation of LIGO-Virgo-KAGRA data suggests that the primary mass function of merging binary black holes (BBHs) at redshift z ≲ 1 contains multiple structures, while spins are relatively low. Theoretical models of BBH formation in different environments can provide a key to interpreting the population of observed mergers, but they require the simultaneous treatment of stellar evolution and dynamics, galaxy evolution, and general relativity. We present B-POP, a population synthesis tool to model BBH mergers formed in the field or via dynamical interactions in young, globular, and nuclear clusters. Using B-POP, we explore how black hole (BH) formation channels, star cluster evolution, hierarchical mergers, and natal BH properties affect the population of BBH mergers. We find that the primary mass distribution of BBH mergers extends beyond $M_1 \simeq 200\, {}$ M_⊙, and the effective spin parameter distribution hints at different natal spins for single and binary BHs. Observed BBHs can be interpreted as members of a mixed population comprised of $\sim 34{{\ \rm per\ cent}} \, {}(66{{\ \rm per\ cent}})$ isolated (dynamical) BBHs, with the latter likely dominating at redshift z > 1. Hierarchical mergers constitute the $4.6\!-\!7.9{{\ \rm per\ cent}}$ of all mergers in the reference model, dominating the primary mass distribution beyond $M_1 \gt 65\, {}$ M_⊙. The inclusion of cluster mass-loss and expansion causes an abrupt decrease in the probability for mergers beyond the third generation to occur. Considering observational biases, we find that $2.7\!-\!7.5{{\ \rm per\ cent}}$ of mock mergers involve intermediate-mass black hole (IMBH) seeds formed via stellar collisions. Comparing this percentage to observed values will possibly help us to constrain IMBH formation mechanisms.