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 more »

Authors:
; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10361831
Journal Name:
The Astrophysical Journal Supplement Series
Volume:
258
Issue:
2
Page Range or eLocation-ID:
Article No. 22
ISSN:
0067-0049
Publisher:
DOI PREFIX: 10.3847
National Science Foundation
##### More Like this
1. 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–120M). 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 500Min 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 thatallof 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 themore »

2. 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 ≳ 10 6 ) and high density. Here we present detailed and self-consistent 47 Tuc models computed with the Cluster 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 in CMC 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 binariesmore »
3. (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 populationmore »
4. Abstract

We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly modeled sample to date. Using theProspectorstellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of$z=0.64−0.32+0.83$(68% confidence), including nine photometric redshifts atz≳ 1. We further find a median mass-weighted age oftm=$0.8−0.53+2.71$Gyr, stellar mass of log(M*/M) =$9.69−0.65+0.75$, star formation rate of SFR =$1.44−1.35+9.37$Myr−1, stellar metallicity of log(Z*/Z) =$−0.38−0.42+0.44$, and dust attenuation of$AV=0.43−0.36+0.85$mag (68% confidence). Overall, the majority of short GRB hosts are star-forming (≈84%), with small fractions that are either transitioning (≈6%) or quiescent (≈10%); however, we observe a much larger fraction (≈40%) of quiescent and transitioning hosts atz≲ 0.25, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies as the general galaxy population, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution,more »

5. Abstract Close encounters between neutron stars and main-sequence stars occur in globular clusters and may lead to various outcomes. Here we study encounters resulting in the tidal disruption of the star. Using N -body models, we predict the typical stellar masses in these disruptions and the dependence of the event rate on the host cluster properties. We find that tidal disruption events occur most frequently in core-collapsed globular clusters and that roughly 25% of the disrupted stars are merger products (i.e., blue straggler stars). Using hydrodynamic simulations, we model the tidal disruptions themselves (over timescales of days) to determine the mass bound to the neutron star and the properties of the accretion disks formed. In general, we find roughly 80%–90% of the initial stellar mass becomes bound to the neutron star following disruption. Additionally, we find that neutron stars receive impulsive kicks of up to about 20 km s −1 as a result of the asymmetry of unbound ejecta; these kicks place these neutron stars on elongated orbits within their host cluster, with apocenter distances well outside the cluster core. Finally, we model the evolution of the (hypercritical) accretion disks on longer timescales (days to years after disruption) to estimatemore »