An official website of the United States government
Abstract The globular cluster 47 Tucanae (47 Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars (N≳ 106) and high density. Here we present detailed and self-consistent 47 Tuc models computed with theCluster Monte Carlocode (CMC). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47 Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions (as in, e.g., Kroupa 2001), and an initial Elson profile (Elson et al. 1987) that is overfilling the cluster’s tidal radius. We include new prescriptions inCMCfor the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47 Tuc; our best-fit model contains ∼50 millisecond pulsars, 70% of which are formed through giant collisions and tidal captures. Our models also suggest that 47 Tuc presently contains up to ∼200 stellar-mass black holes, ∼5 binary black holes, ∼15 low-mass X-ray binaries, and ∼300 cataclysmic variables.
more »
« less
Modeling Dense Star Clusters in the Milky Way and beyond with the Cluster Monte Carlo Code
Abstract We describe the public release of the Cluster Monte Carlo (CMC) code, a parallel, star-by-starN-body code for modeling dense star clusters.CMCtreats collisional stellar dynamics using Hénon’s method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release ofCMCis pinned directly to theCOSMICpopulation synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest (N= 108) star-by-starN-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap).
more »
« less
- PAR ID:
- 10361831
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Supplement Series
- Volume:
- 258
- Issue:
- 2
- ISSN:
- 0067-0049
- Format(s):
- Medium: X Size: Article No. 22
- Size(s):
- Article No. 22
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
ABSTRACT The current generation of galaxy simulations can resolve individual giant molecular clouds, the progenitors of dense star clusters. But the evolutionary fate of these young massive clusters, and whether they can become the old globular clusters (GCs) observed in many galaxies, is determined by a complex interplay of internal dynamical processes and external galactic effects. We present the first star-by-star N-body models of massive (N ∼ 105–107) star clusters formed in a FIRE-2 MHD simulation of a Milky Way-mass galaxy, with the relevant initial conditions and tidal forces extracted from the cosmological simulation. We select 895 (∼30 per cent) of the YMCs with >6 × 104 M⊙ from Grudić et al. 2022 and integrate them to z = 0 using the cluster Monte Carlo code, CMC. This procedure predicts a MW-like system with 148 GCs, predominantly formed during the early, bursty mode of star formation. Our GCs are younger, less massive, and more core-collapsed than clusters in the Milky Way or M31. This results from the assembly history and age-metallicity relationship of the host galaxy: Younger clusters are preferentially born in stronger tidal fields and initially retain fewer stellar-mass black holes, causing them to lose mass faster and reach core collapse sooner than older GCs. Our results suggest that the masses and core/half-light radii of GCs are shaped not only by internal dynamical processes, but also by the specific evolutionary history of their host galaxies. These results emphasize that N-body studies with realistic stellar physics are crucial to understanding the evolution and present-day properties of GC systems.more » « less
-
Abstract We study the evolution of populations of binary stars within massive cluster-forming regions. We simulate the formation of young massive star clusters within giant molecular clouds with masses ranging from 2 × 104to 3.2 × 105M⊙. We use Torch, which couples stellar dynamics, magnetohydrodynamics, star and binary formation, stellar evolution, and stellar feedback through the Amuseframework. We find that the binary fraction decreases during cluster formation at all molecular cloud masses. The binaries’ orbital properties also change, with stronger and quicker changes in denser, more massive clouds. Most of the changes we see can be attributed to the disruption of binaries wider than 100 au, although the close binary fraction also decreases in the densest cluster-forming region. The binary fraction for O stars remains above 90%, but exchanges and dynamical hardening are ubiquitous, indicating that O stars undergo frequent few-body interactions early during the cluster formation process. Changes to the populations of binaries are a by-product of hierarchical cluster assembly: most changes to the binary population take place when the star formation rate is high, and there are frequent mergers between subclusters in the cluster-forming region. A universal primordial binary distribution based on observed inner companions in the Galactic field is consistent with the binary populations of young clusters with resolved stellar populations, and the scatter between clusters of similar masses could be explained by differences in their formation history.more » « less
