Assembling GW231123 in star clusters through the combination of stellar binary evolution and hierarchical mergers

Abstract: GW231123 is the most massive binary black hole (BBH) merger detected to date by the LIGO-Virgo-KAGRA collaboration. With at least one black hole (BH) in the upper-mass gap and both BHs exhibiting high spins ($\chi_{1,2} \gtrsim 0.8$), this event challenges standard isolated binary evolution models. A compelling alternative is a dynamical origin in star clusters, where stellar binaries and hierarchical mergers may both contribute to the formation of similar BBHs. In this work, we investigate the formation of GW231123-like events in different cluster environments using the B-POP semi-analytic population synthesis code. We find that low-metallicity environments ($Z \lesssim 0.002$) are ideal for producing BBH mergers similar to GW231123. In young and globular clusters, these BBHs have components formed in stellar binaries, whilst in nuclear clusters there is also a significant contribution from BHs built-up via hierarchical mergers. Natal spins of BHs formed in stellar binaries are crucial to find GW231123 analogs. In particular, our models suggest that BHs from stellar binaries are likely characterized by high-spins. Simulated GW231123-like systems exhibit short delay times, $t_\mathrm{del} \sim 0.1 - 1$ Gyr, which suggests their progenitors formed close to the inferred merger redshift ($z = 0.39^{+0.27}_{-0.24}$). We argue that star clusters in metal-poor dwarf galaxies or Milky Way-like galaxies are ideal nurseries, inferring an upper limit to the local merger rate of $\mathcal{R} \sim 1.6\times10^{-3} - 0.16$ yr$^{-1}$ Gpc$^{-3}$ for nuclear clusters, $\sim 0.036 - 0.72$ yr$^{-1}$ Gpc$^{-3}$ for globular clusters, and $4\times10^{-4}-0.041$ yr$^{-1}$ Gpc$^{-3}$ for young clusters.