2026ApJ...998..138M

The LIGO/Virgo/Kagra (LVK) Collaboration has detected numerous binary black hole (BBH) mergers with properties that challenge standard binary evolution scenarios, such as component masses above the pair-instability gap and high spin magnitudes. Dense stellar environments such as globular clusters provide a natural channel for producing such systems through hierarchical mergers, where black hole (BH) remnants formed in earlier mergers are retained in the cluster and undergo successive mergers. However, gravitational-wave recoil kicks often eject merger remnants from typical globular clusters, which limits hierarchical growth. Massive clusters with deeper potential wells, such as those found in giant elliptical galaxies like M87, may overcome this barrier, but direct simulations of such massive globular clusters remain computationally challenging. In this study, we present a 10-million-body cluster simulation performed with the Cluster Monte Carlo (CMC) code, referred to as colossus, which serves as a proxy for the most massive low-metallicity globular clusters observed in the local Universe. This simulation demonstrates that extended chains of hierarchical mergers can occur in massive globular clusters, producing BHs up to the fifth generation with masses approaching 250 M_⊙, comparable to the most massive LVK events observed to date (e.g., GW231123). Combining the colossus simulation with the previous CMC Cluster Catalog, we develop a framework to extrapolate BBH merger predictions for the thousands of globular clusters seen in the Virgo Supercluster.