Hierarchical binary black hole mergers in globular clusters: Mass function and evolution with redshift

Hierarchical black hole (BH) mergers are one of the most straightforward mechanisms producing BHs inside and above the pair-instability mass gap. We investigated the impact of globular cluster (GC) evolution on hierarchical mergers, accounting for the uncertainties related to BH mass pairing functions on the predicted primary BH mass, mass ratio, and spin distribution. We find that the evolution of the host GC quenches the hierarchical BH assembly at the third generation, mainly due to cluster expansion powered by a central BH subsystem. Hierarchical mergers match the primary BH mass distribution from GW events for m1 > 50 M⊙ regardless of the assumed BH pairing function. At lower masses, however, different pairing functions lead to dramatically different predictions on the primary BH mass merger-rate density. We find that the primary BH mass distribution evolves with redshift, with a larger contribution from mergers with m1 ≥30 M⊙ for z ≥2. Finally, we calculate the mixing fraction of binary black holes (BBHs) from GCs and isolated binary systems. Our predictions are very sensitive to the spins, which favor a large fraction (> 0.6) of BBHs born in GCs in order to reproduce misaligned spin observations.