We study the dynamics of stellar-mass black holes (BH) in star clusters with particular attention to the formation of BH–BH binaries, which are inter

The study of stellar-remnant black holes (BH) in dense stellar clusters is now in the spotlight, especially due to their intrinsic ability to form binary black holes (BBH) through dynamical encounters, which potentially coalesce via gravitational-wave (GW) radiation. In this work, which is a continuation from a recent study (Paper I), additional models of compact stellar clusters with initial masses ≲ 10⁵ M_⊙ and also those with small fractions of primordial binaries (≲ 10 per cent) are evolved for long term, applying the direct N-body approach, assuming state-of-the-art stellar-wind and remnant-formation prescriptions. That way, a substantially broader range of computed models than that in Paper I is achieved. As in Paper I, the general-relativistic BBH mergers continue to be mostly mediated by triples that are bound to the clusters rather than happen among the ejected BBHs. In fact, the number of such in situ BBH mergers, per cluster, tends to increase significantly with the introduction of a small population of primordial binaries. Despite the presence of massive primordial binaries, the merging BBHs, especially the in situ ones, are found to be exclusively dynamically assembled and hence would be spin–orbit misaligned. The BBHs typically traverse through both the LISA's and the LIGO's detection bands, being audible to both instruments. The ‘dynamical heating’ of the BHs keeps the electron-capture-supernova (ECS) neutron stars (NS) from effectively mass segregating and participating in exchange interactions; the dynamically active BHs would also exchange into any NS binary within ≲1 Gyr. Such young massive and open clusters have the potential to contribute to the dynamical BBH merger detection rate to a similar extent as their more massive globular-cluster counterparts.

Context. As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. Aims. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N-body evolution program NBODY7. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. Methods. All the above newly implemented schemes are also shared with the standalone versions of SSE and BSE. All these demonstrations are performed with both the updated standalone BSE and the updated NBODY7/BSE. Results. When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing. Conclusions. Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.

Stellar-remnant black holes (BH) in dense stellar clusters have always drawn attention due to their potential in a number of phenomena, especially the dynamical formation of binary black holes (BBH), which potentially coalesce via gravitational-wave radiation. This study presents a preliminary set of evolutionary models of compact stellar clusters with initial masses ranging over 1.0 × 10⁴–5.0 × 10⁴ M_⊙, and half-mass radius of 2 or 1 pc, which is typical for young massive and starburst clusters. They have metallicities between 0.05 Z_⊙ and Z_⊙. Including contemporary schemes for stellar wind and remnant formation, such model clusters are evolved, for the first time, using the state-of-the-art direct N-body evolution program nbody7, until their dissolution or at least for 10 Gyr. That way, a self-regulatory behaviour in the effects of dynamical interactions among the BHs is demonstrated. In contrast to earlier studies, the BBH coalescences obtained in these models show a prominence in triple-mediated coalescences while being bound to the clusters, compared to those occurring among the BBHs that are dynamically ejected from the clusters. A broader mass spectrum of the BHs and lower escape velocities of the clusters explored here might cause this difference, which is yet to be fully understood. Among the BBH coalescences obtained here, there are ones that resemble the detected GW151226, LVT151012 and GW150914 events and also ones that are even more massive. A preliminary estimate suggests few 10–100 s of BBH coalescences per year, originating due to dynamics in stellar clusters that can be detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) at its design sensitivity.

I present a set of long-term, direct, relativistic many-body computations of model dense stellar clusters with up-to-date stellar-evolutionary, supernova (SN), and remnant natal-kick models, including pair instability and pulsation pair instability supernova (PSN and PPSN), using an updated version of ${\rm{\small NBODY7}}$ N-body simulation program. The N-body model also includes stellar evolution-based natal spins of black holes (BHs) and treatments of binary black hole (BBH) mergers based on numerical relativity. These, for the first time in a direct N-body simulation, allow for second-generation BBH mergers. The set of 65 evolutionary models have initial masses $10^4{\!-\!}10^5\, \mathrm{M}_{\odot }$, sizes 1–3 pc, metallicity 0.0001–0.02, with the massive stars in primordial binaries and they represent young massive clusters (YMC) and moderately massive open clusters (OC). Such models produce dynamically paired BBH mergers that agree well with the observed masses, mass ratios, effective spin parameters, and final spins of the LVC O1/O2 merger events, provided BHs are born with low or no spin but spin-up after undergoing a BBH merger or matter accretion on to it. In particular, the distinctly higher mass, effective spin parameter, and final spin of GW170729 merger event is naturally reproduced, as also the mass asymmetry of the O3 event GW190412. The computed models produce intermediate-mass, $\sim 100\, \mathrm{M}_{\odot }$ BBH mergers with primary mass within the ‘PSN gap’ and also yield mergers involving remnants in the ‘mass gap’. They also suggest that YMCs and OCs produce persistent, Local-Universe GW sources detectable by LISA. Such clusters are also capable of producing eccentric LIGO-Virgo mergers.