Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 15th until 12:00 AM ET on Friday, October 16th due to maintenance. We apologize for the inconvenience.


Title: Great Balls of FIRE: III. Modelling black hole mergers from massive star clusters in simulations of galaxies
Context.Despite the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Aims.Here we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. Methods.To this end, we applied a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment, with a unique star formation history, that hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolutionCMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. Results.As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favouring the creation of black holes, low-mass galaxies contain few massive clusters and therefore make a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear, identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations.  more » « less
Award ID(s):
2310362
PAR ID:
10556484
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
EDP Sciences
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
686
ISSN:
0004-6361
Page Range / eLocation ID:
A106
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract GW231123, the most massive binary black hole (BBH) merger detected by LIGO–Virgo–KAGRA, highlights the need to understand the origins of massive, high-spin stellar black holes (BHs). Dense star clusters provide natural environments for forming such systems, beyond the limits of standard massive star evolution to core collapse. While repeated BBH mergers can grow BHs through dynamical interactions (the so-called “hierarchical merger” channel), most star clusters with masses ≲106M⊙have escape speeds too low to retain higher-generation BHs, limiting growth into or beyond the mass gap. In contrast, BH–star collisions with subsequent accretion of the collision debris can grow and retain BHs irrespective of the cluster escape speed. UsingN-body (Cluster Monte Carlo) simulations, we study BH growth and spin evolution through this process, and we find that accretion can drive BH masses up to at least ∼200M⊙, with spins set by the details of the growth history. BHs up to about 150M⊙can reach dimensionless spinsχ ≳ 0.7 via single coherent episodes, while more massive BHs form through multiple stochastic accretion events and eventually spin down toχ ≲ 0.4. These BHs later form binaries through dynamical encounters, producing BBH mergers that contribute up to ∼10% of all detectable events, comparable to predictions for the hierarchical channel. However, the two pathways predict distinct signatures: hierarchical mergers yield more unequal mass ratios, whereas accretion-grown BHs preferentially form near-equal-mass binaries. The accretion-driven channel allows dense clusters with low escape speeds, such as globular clusters, to produce highly spinning BBHs with both components in or above the mass gap, providing a natural formation pathway to GW231123-like systems. 
    more » « less
  2. Abstract Merging binary black holes (BBHs) formed dynamically in dense star clusters are expected to have uncorrelated spin–orbit orientations since they are assembled through many random interactions. However, measured effective spins in BBHs detected by LIGO/Virgo/KAGRA hint at additional physical processes that may introduce anisotropy. Here we address this question by exploring the impact of stellar collisions and accretion of collision debris on the spin–orbit alignment in merging BBHs formed in dense star clusters. Through hydrodynamic simulations, we study the regime where the disruption of a massive star by a BBH causes the stellar debris to form individual accretion disks bound to each black hole (BH). We show that these disks, which are randomly oriented relative to the binary orbital plane after the initial disruption of the star, can be reoriented by strong tidal torques in the binary near pericenter passages. Following accretion by the BHs on longer timescales, BBHs with small but preferentially positive effective spin parameters (χeff≲ 0.2) are formed. Our results indicate that BBH collisions in young massive star clusters could contribute to the observed trend toward small positiveχeff, and we suggest that the standard assumption often made that dynamically assembled BBHs should have isotropically distributed BH spins is not always justified. 
    more » « less
  3. Abstract With about one hundred mergers of binary black holes (BBHs) detected via gravitational waves by the LIGO-Virgo-KAGRA (LVK) Collaboration, our understanding of the darkest objects in the universe has taken unparalleled steps forward. While most of the events are expected to consist of black holes (BHs) directly formed from the collapse of massive stars, some may contain the remnants of previous BBH mergers. In the most massive globular clusters and in nuclear star clusters, successive mergers can produce second- (2G) or higher-generation BHs, and even form intermediate-mass BHs (IMBHs). Overall, we predict that up to ∼10%, ∼1%, or ∼0.1% of the BBH mergers have one component being a 2G, 3G, or 4G BH, respectively. Assuming that ∼500 BBH mergers will be detected in O4 by LVK, this means that ∼50, ∼5, or ∼0.5 events, respectively, will involve a 2G, 3G, or 4G BH, if most sources are produced dynamically in dense star clusters. With their distinctive signatures of higher masses and spins, such hierarchical mergers offer an unprecedented opportunity to learn about the BH populations in the densest stellar systems and to shed light on the elusive IMBHs that may form therein. 
    more » « less
  4. Abstract The active galactic nucleus (AGN) channel for the formation of binary black hole (BBH) mergers has been previously studied as a potential formation channel for the merging compact binaries observed by the LIGO–Virgo–KAGRA (LVK) scientific collaboration. The first two papers in this series explored theMcFACTScode for the evolution of black hole orbits in AGN accretion disks for individual galaxy models and described the characteristics of predicted BBH populations in realizations of those models (such as the correlation between mass ratio and aligned spin). In this work, we explore the impact of the properties of AGN host galaxies and assume an AGN lifetime and cosmological model for the density of AGN in a universe like our own. By sampling from an inferred population of AGN, we marginalize over galaxy mass to predict a population of BBH mergers observable by modern ground-based gravitational-wave observatories. We find that for reasonable assumptions, AGN disk environments may account for massive BBH mergers such as GW190521 and GW190929_012149. We find that the majority of observable BBH mergers from our simulation are expected to originate in galaxies with a supermassive black hole between 107M⊙and 109.4M⊙. We also find that if hierarchical mergers from AGN disks account for a substantial part of the LVK population, our current models require an AGN lifetime of 0.5–2.5 Myr. 
    more » « less
  5. The gravitational wave event GW231123, with component black hole masses lying within or above the pair-instability mass gap, poses a significant challenge to current stellar evolution models. In this work, we describe how we investigated its origin by coupling the galaxy formation modelGAMESHwith the cluster population synthesis codeRAPSTERand using two distinct binary population synthesis codes,SEVNandBSEEMP. This framework allowed us, for the first time, to reconstruct the life cycle of GW231123-like candidates within the same cosmological simulation, enabling a self-consistent comparison between different formation channels. Although both population synthesis codes can in principle produce black holes compatible with GW231123, we find that isolated binary evolution fails to reproduce the inferred merger redshift. InSEVN, massive black hole binaries form with semi-major axes > 103R⊙, preventing coalescences within a Hubble time. InBSEEMP, candidates arise only at extremely low metallicities (Z ≈ 10−10), which contribute negligibly to the star formation rate density in our overdense simulated volume. Our results therefore strongly support a dynamical hierarchical origin. The observed black hole masses are naturally reproduced through successive mergers in dense globular clusters. The high dimensionless spins reported by the LIGO-Virgo-KAGRA Collaboration are consistent with this hierarchical population. We found a local merger rate density of 0.78 Gpc−3yr−1, with a peak atz = 4 − 6, tracing the maximum formation rate of globular clusters in metal-poor environments (Z ≈ 0.006). Overall, GW231123 may represent a benchmark event for a robust population of hierarchical black holes formed in the early Universe. 
    more » « less