skip to main content

Title: Modeling Dense Star Clusters in the Milky Way and beyond with the Cluster Monte Carlo Code

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 » recently reported merger with component masses in the pulsational pair-instability mass gap).

« less
; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
2009916 2108624 2001751
Publication Date:
Journal Name:
The Astrophysical Journal Supplement Series
Page Range or eLocation-ID:
Article No. 22
DOI PREFIX: 10.3847
Sponsoring Org:
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 »peak IMBH merger rate to be2Gpc3yr1at redshiftz≈ 2.

    « less
  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 »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.« less
  3. 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 populationmore »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.« less
  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 ofz=0.640.32+0.83(68% confidence), including nine photometric redshifts atz≳ 1. We further find a median mass-weighted age oftm=0.80.53+2.71Gyr, stellar mass of log(M*/M) =9.690.65+0.75, star formation rate of SFR =1.441.35+9.37Myr−1, stellar metallicity of log(Z*/Z) =0.380.42+0.44, and dust attenuation ofAV=0.430.36+0.85mag (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 »with a fast-merging channel atz> 1 and a decreased neutron star binary formation efficiency from high to low redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the Broadband Repository for Investigating Gamma-ray burst Host Traits website.

    « less
  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 »the accretion rate onto the neutron stars and accompanying spin-up. As long as ≳1% of the bound mass accretes onto the neutron star, millisecond spin periods can be attained. We argue the growing numbers of isolated millisecond pulsars observed in globular clusters may have formed, at least in part, through this mechanism. In the case of significant mass growth, some of these neutron stars may collapse to form low-mass (≲3 M ⊙ ) black holes.« less