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 (
We describe the public release of the Cluster Monte Carlo (
- NSF-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
-
more » « less
Abstract N ≳ 106) and high density. Here we present detailed and self-consistent 47 Tuc models computed with theCluster Monte Carlo code (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 inCMC for 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
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 The existence of black holes (BHs) with masses in the range between stellar remnants and supermassive BHs has only recently become unambiguously established. GW190521, a gravitational wave signal detected by the LIGO/Virgo Collaboration, provides the first direct evidence for the existence of such intermediate-mass BHs (IMBHs). This event sparked and continues to fuel discussion on the possible formation channels for such massive BHs. As the detection revealed, IMBHs can form via binary mergers of BHs in the “upper mass gap” (≈40–120
M ⊙). Alternatively, IMBHs may form via the collapse of a very massive star formed through stellar collisions and mergers in dense star clusters. In this study, we explore the formation of IMBHs with masses between 120 and 500M ⊙in young, massive star clusters using state-of-the-art Cluster Monte Carlo models. We examine the evolution of IMBHs throughout their dynamical lifetimes, ending with their ejection from the parent cluster due to gravitational radiation recoil from BH mergers, or dynamical recoil kicks from few-body scattering encounters. We find thatall of the IMBHs in our models are ejected from the host cluster within the first ∼500 Myr, indicating a low retention probability of IMBHs in this mass range for globular clusters today. We estimate the peak IMBH merger rate to bez ≈ 2. -
more » « less
Abstract Close encounters between stellar-mass black holes (BHs) and stars occur frequently in dense star clusters and in the disks of active galactic nuclei. Recent studies have shown that in highly eccentric close encounters, the star can be tidally disrupted by the BH in a microtidal disruption event (microTDE), resulting in rapid mass accretion and possibly bright electromagnetic signatures. Here we consider a scenario in which the star might approach the stellar-mass BH in a gradual, nearly circular inspiral, under the influence of dynamical friction in a circum-binary gas disk or three-body interactions in a star cluster. We perform hydrodynamics simulations of this scenario using the smoothed particle hydrodynamics code
PHANTOM . We find that under certain circumstances (for initial eccentricitye 0≳ 0.4 and penetration factorβ = 1, ore 0< 0.4 andβ ≲ 0.67), the mass of the star is slowly stripped away by the BH. We call this gradual tidal disruption a “tidal-peeling event.” Additionally, we discover that some low-eccentricity microTDEs (e 0< 0.4 andβ = 1) are a new form of fast luminous transients similar to parabolic microTDEs. Depending on the initial distance and eccentricity of the encounter, these low-eccentricity microTDEs might exhibit significant accretion rates and orbital evolution distinct from those of a typical (eccentric) microTDE. -
null (Ed.)ABSTRACT The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled magnetohydrodynamics and direct N-body star cluster formation code torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the q ≥ 0.1 detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semimajor axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.more » « less
