Black hole mergers from dense star clusters with realistic binary populations

Abstract

We present a suite of 24 full-lifetime simulations of dense star clusters with the Cluster Monte Carlo (CMC) code, featuring updated input physics and a realistic distribution of initial binary systems. The latter encompasses a mass-dependent binary fraction, period distribution, and eccentricity distribution based on observations of well-studied stellar populations in the Solar neighborhood and nearby star-forming regions. We predict the cosmic rate, masses, and spins of binary black hole (BBH) mergers formed through dynamical assembly, primordial binary evolution, and hierarchical mergers within dense clusters. As with previous model grids with fewer binaries, dynamically assembled first-generation (1G) mergers dominate the rate of cluster-derived mergers, and the total merger rate is consistent with that inferred from LIGO-Virgo-KAGRA observations as of GWTC-5.0. Our models naturally reproduce key features of the inferred BBH population, including the broken-power-law behavior of the primary BH mass spectrum for $m_1 \gtrsim 20 M_\odot$, the shallower (steeper) slope of the secondary mass spectrum relative to the primary for $m_2 \lesssim 10 M_\odot$ ($m_2 \gtrsim 30 M_\odot$), and the shape of the mass-ratio distribution in the low- and high-mass domains. We predict broad distributions of the spin parameters $χ_{\mathrm{eff}}$ and $χ_{\mathrm{p}}$, consistent with previous studies of dynamical assembly in clusters. The merger rate from primordial binary systems within clusters is a small fraction of the total; however, their merger products are frequently involved in subsequent hierarchical mergers, with the result that the hierarchical merger rate evolves more steeply than the 1G dynamical merger rate with redshift.