GW190521 formation via three-body encounters in young massive star clusters

GW190521 is the most massive binary black hole (BBH) merger observed to date, and its primary component lies in the pair-instability (PI) mass gap. Here, we investigate the formation of GW190521-like systems via three-body encounters in young massive star clusters. We performed 2× 105 simulations of binary-single interactions between a BBH and a massive ≥60 M⊙ black hole (BH), including post-Newtonian terms up to the 2.5 order and a prescription for relativistic kicks. In our initial conditions, we take into account the possibility of forming BHs in the PI mass gap via stellar collisions. If we assume that first-generation BHs have low spins, $\sim {0.17}{{\ \rm per\ cent}}$ of all the simulated BBH mergers have component masses, effective and precessing spin, and remnant mass and spin inside the $90{{\ \rm per\ cent}}$ credible intervals of GW190521. Seven of these systems are first-generation exchanged binaries, while five are second-generation BBHs. We estimate a merger rate density $\mathcal {R}_{\rm GW190521}\sim {0.03}\,$Gpc−3 yr−1 for GW190521-like binaries formed via binary-single interactions in young star clusters. This rate is extremely sensitive to the spin distribution of first-generation BBHs. Stellar collisions, second-generation mergers and dynamical exchanges are the key ingredients to produce GW190521-like systems in young star clusters.

Mass and Rate of Hierarchical Black Hole Mergers in Young, Globular and Nuclear Star Clusters

Hierarchical mergers are one of the distinctive signatures of binary black hole (BBH) formation through dynamical evolution. Here, we present a fast semi-analytic approach to simulate hierarchical mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs). Hierarchical mergers are more common in NSCs than they are in both GCs and YSCs because of the different escape velocity. The mass distribution of hierarchical BBHs strongly depends on the properties of first-generation BBHs, such as their progenitor’s metallicity. In our fiducial model, we form black holes (BHs) with masses up to ∼103 M⊙ in NSCs and up to ∼102 M⊙ in both GCs and YSCs. When escape velocities in excess of 100 km s−1 are considered, BHs with mass >103 M⊙ are allowed to form in NSCs. Hierarchical mergers lead to the formation of BHs in the pair instability mass gap and intermediate-mass BHs, but only in metal-poor environments. The local BBH merger rate in our models ranges from ∼10 to ∼60 Gpc−3 yr−1; hierarchical BBHs in NSCs account for ∼10−2–0.2 Gpc−3 yr−1, with a strong upper limit of ∼10 Gpc−3 yr−1. When comparing our models with the second gravitational-wave transient catalog, we find that multiple formation channels are favored to reproduce the observed BBH population.

New insights on binary black hole formation channels after GWTC-2: young star clusters versus isolated binaries

With the recent release of the second gravitational-wave transient catalogue (GWTC-2), which introduced dozens of new detections, we are at a turning point of gravitational wave astronomy, as we are now able to directly infer constraints on the astrophysical population of compact objects. Here, we tackle the burning issue of understanding the origin of binary black hole (BBH) mergers. To this effect, we make use of state-of-the-art population synthesis and N-body simulations, to represent two distinct formation channels: BBHs formed in the field (isolated channel) and in young star clusters (dynamical channel). We then use a Bayesian hierarchical approach to infer the distribution of the mixing fraction f, with f = 0 (f = 1) in the pure dynamical (isolated) channel. We explore the effects of additional hyper-parameters of the model, such as the spread in metallicity σZ and the parameter σsp, describing the distribution of spin magnitudes. We find that the dynamical model is slightly favoured with a median value of f = 0.26, when σsp = 0.1 and σZ = 0.4. Models with higher spin magnitudes tend to strongly favour dynamically formed BBHs (f ≤ 0.1 if σsp = 0.3). Furthermore, we show that hyper-parameters controlling the rates of the model, such as σZ, have a large impact on the inference of the mixing fraction, which rises from 0.18 to 0.43 when we increase σZ from 0.2 to 0.6, for a fixed value of σsp = 0.1. Finally, our current set of observations is better described by a combination of both formation channels, as a pure dynamical scenario is excluded at the $99{{\ \rm per\ cent}}$ credible interval, except when the spin magnitude is high.

Intermediate-mass black holes from stellar mergers in young star clusters

ABSTRACT Intermediate-mass black holes (IMBHs) in the mass range $10^2\!-\!10^5\, \mathrm{M_{\odot }}$ bridge the gap between stellar black holes (BHs) and supermassive BHs. Here, we investigate the possibility that IMBHs form in young star clusters via runaway collisions and BH mergers. We analyse 104 simulations of dense young star clusters, featuring up-to-date stellar wind models and prescriptions for core collapse and (pulsational) pair instability. In our simulations, only nine IMBHs out of 218 form via binary BH mergers, with a mass ∼100–140 M⊙. This channel is strongly suppressed by the low escape velocity of our star clusters. In contrast, IMBHs with masses up to ∼438 M⊙ efficiently form via runaway stellar collisions, especially at low metallicity. Up to ∼0.2 per cent of all the simulated BHs are IMBHs, depending on progenitor’s metallicity. The runaway formation channel is strongly suppressed in metal-rich (Z = 0.02) star clusters, because of stellar winds. IMBHs are extremely efficient in pairing with other BHs: ∼70 per cent of them are members of a binary BH at the end of the simulations. However, we do not find any IMBH–BH merger. More massive star clusters are more efficient in forming IMBHs: ∼8 per cent (∼1 per cent) of the simulated clusters with initial mass 104–3 × 104 M⊙ (103–5 × 103 M⊙) host at least one IMBH.