Searching for distinctive signatures, which characterize different formation channels of binary black holes (BBHs), is a crucial step towards the interpretation of current and future gravitational wave detections. Here, we investigate the demography of merging BBHs in young star clusters (SCs), which are the nursery of massive stars. We performed 4 × 10 3 N-body simulations of SCs with metallicity Z = 0.002, initial binary fraction 0.4 and fractal initial conditions, to mimic the clumpiness of star forming regions. Our simulations include a novel population-synthesis approach based on the code MOBSE. We find that SC dynamics does not affect the merger rate significantly, but leaves a strong fingerprint on the properties of merging BBHs. More than 50 % of merging BBHs in young SCs form by dynamical exchanges in the first few Myr. Dynamically formed merging BBHs are significantly heavier than merging BBHs in isolated binaries: merging BBHs with total mass up to ∼ 120 M form in young SCs, while the maximum total mass of merging BBHs in isolated binaries with the same metallicity is only ∼ 70 M . Merging BBHs born via dynamical exchanges tend to have smaller mass ratios than BBHs in isolated binaries. Furthermore, SC dynamics speeds up the merger: the delay time between star formation and coalescence is significantly shorter in young SCs. In our simulations, massive systems such as GW170729 form only via dynamical exchanges. Finally ∼ 2 % of merging BBHs in young SCs have mass in the pair-instability mass gap (∼ 60 − 120 M ). This represents a unique fingerprint of merging BBHs in SCs.

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z = 0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius ρh ≥ 3.4 × 104 and ≥1.5 × 102 M⊙ pc−3 in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find ∼6, ∼2, and <1 % of BHs with mass mBH > 60 M⊙ at Z = 0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to ∼320 M⊙. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to ∼140 M⊙ (∼80 M⊙). The most massive BBH merger in our simulations has primary mass ∼88 M⊙, inside the pair-instability mass gap, and a mass ratio of ∼0.5. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density ∼110 and ∼55 Gpc−3 yr−1, if we assume that all stars form in loose and dense star clusters, respectively.