Predictions of a simple parametric model of hierarchical black hole mergers

Hierarchical mergers of black holes are proposed as a mechanism to explain the observations of binary black holes with component masses between $\sim 50M_{\odot}\hbox{--}130M_{\odot}$ by LIGO/Virgo, often referred to as "upper mass gap". We study the efficiency with which hierarchical mergers can produce higher and higher masses using a simple model of the forward evolution of binary black hole populations in gravitationally bound systems like stellar clusters. The model relies on pairing probability and initial mass functions for the black hole population, along with numerical relativity fitting formulas for the mass, spin, and kick speed of the merger remnant. We carry out an extensive comparison of the predictions of our model with clusterBHBdynamics (cBHBd) model, a fast method for the evolution of star clusters and black holes therein. For this comparison, we consider three different pairing functions of black holes and consider simulations from high- and low-metallicity cluster environments from cBHBd. We find good agreements between our model and the cBHBd results when the pairing probability of binaries depends on both total mass and mass ratio. We also assess the efficiency of hierarchical mergers as a function of merger generation and derive the mass distribution of black holes using our model. We find that the multi-modal features in the observed binary black hole mass spectrum -- revealed by the non-parametric population models -- can be interpreted by invoking the hierarchical merger scenario in dense, metal-rich, stellar environments. Further, the two subdominant peaks in the GWTC-3 component mass spectrum are consistent with second and third-generation mergers in metal-rich, dense environments. With more binary black hole detections, our model could be used to infer the black hole initial mass function and pairing probability exponent.