An official website of the United States government
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
Stellar graveyards: clustering of compact objects in globular clusters NGC 3201 and NGC 6397
ABSTRACT We analyse Gaia EDR3 and re-calibrated HST proper motion data from the core-collapsed and non-core-collapsed globular clusters NGC 6397 and NGC 3201, respectively, with the Bayesian mass-orbit modelling code MAMPOSSt-PM. We use Bayesian evidence and realistic mock data sets constructed with Agama to select between different mass models. In both clusters, the velocities are consistent with isotropy within the extent of our data. We robustly detect a dark central mass (DCM) of roughly $$1000\, \rm M_\odot$$ in both clusters. Our MAMPOSSt-PM fits strongly prefer an extended DCM in NGC 6397, while only presenting a mild preference for it in NGC 3201, with respective sizes of a roughly one and a few per cent of the cluster effective radius. We explore the astrophysics behind our results with the CMC Monte Carlo N-body code, whose snapshots best matching the phase space observations lead to similar values for the mass and size of the DCM. The internal kinematics are thus consistent with a population of hundreds of massive white dwarfs in NGC 6397, and roughly 100 segregated stellar-mass black holes in NGC 3201, as previously found with CMC. Such analyses confirm the accuracy of both mass-orbit modelling and Monte Carlo N-body techniques, which together provide more robust predictions on the DCM of globular clusters (core-collapsed or not). This opens possibilities to understand a vast range of interesting astrophysical phenomena in clusters, such as fast radio bursts, compact object mergers, and gravitational waves.
more »
« less
- Award ID(s):
- 2001751
- PAR ID:
- 10421266
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 514
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 806-825
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract We estimate the absolute age of the globular cluster NGC 3201 using 10,000 sets of theoretical isochrones constructed through Monte Carlo simulation using the Dartmouth Stellar Evolution Program. These isochrones take into consideration the uncertainty introduced by the choice of stellar evolution parameters. We fit isochrones with three detached eclipsing binaries and obtained an age independent of distance. We also fit isochrones with differential reddening corrected Hubble Space Telescope photometry data utilizing two different Hess diagram-based fitting methods. Results from three different methods analyzing two different types of data agree to within 1σ, and we find the absolute age of NGC 3201 = 11.85 ± 0.74 Gyr. We also perform a variable importance analysis to study the uncertainty contribution from individual parameters, and we find the distance is the dominant source of uncertainty in photometry-based analysis, while total metallicity, helium abundance,α-element abundance, mixing length, and treatment of helium diffusion are an important source of uncertainties for all three methods.more » « less
-
ABSTRACT As an introduction of a kinematic survey of Magellanic Cloud (MC) star clusters, we report on the dynamical masses and mass-to-light ratios (M/L) of NGC 419 (Small Magellanic Cloud) and NGC 1846 (Large Magellanic Cloud). We have obtained more than one hundred high-resolution stellar spectra in and around each cluster using the multi-object spectrograph M2FS on the Magellan/Clay Telescope. Line-of-sight velocities and positions of the stars observed in each cluster were used as input to an expectation-maximization algorithm used to estimate cluster membership probabilities, resulting in samples of 46 and 52 likely members (PM ≥ 50 per cent) in NGC 419 and NGC 1846, respectively. This process employed single-mass King models constrained by the structural parameters of the clusters and provided self-consistent dynamical mass estimates for both clusters. Our best-fitting results show that NGC 419 has a projected central velocity dispersion of $$2.44^{+0.37}_{-0.21}$$ km s−1, corresponding to a total mass of $$7.6^{+2.5}_{-1.3}\times 10^4\ {\rm M}_{\odot }$$ and V-band M/L ratio of $$0.22^{+0.08}_{-0.05}$$ in solar units. For NGC 1846, the corresponding results are $$2.04^{+0.28}_{-0.24}$$ km s−1, $$5.4^{+1.5}_{-1.4}\times 10^4\ {\rm M}_{\odot }$$, and $$0.32^{+0.11}_{-0.11}$$. The mean metallicities of NGC 419 and NGC 1846 are found to be $$\rm [Fe/H]=-0.84\pm 0.19$$ and −0.70 ± 0.08, respectively, based on the spectra of likely cluster members. We find marginal statistical evidence of rotation in both clusters, though in neither cluster does rotation alter our mass estimates significantly. We critically compare our findings with those of previous kinematic studies of these two clusters in order to evaluate the consistency of our observational results and analytic tools.more » « less
-
Abstract We used high-resolution spectra acquired with the Magellan Telescope to measure radial and rotational velocities of approximately 200 stars in the Galactic globular cluster NGC 3201. The surveyed sample includes blue straggler stars (BSSs) and reference stars in different evolutionary stages (main-sequence turnoff, subgiant, red giant, and asymptotic giant branches). The average radial velocity value (〈Vr〉 = 494.5 ± 0.5 km s−1) confirms a large systemic velocity for this cluster and was used to distinguish 33 residual field interlopers. The final sample of member stars has 67 BSSs and 114 reference stars. Similarly to what is found in other clusters, the totality of the reference stars has negligible rotation (< 20 km s−1), while the BSS rotational velocity distribution shows a long tail extending up to ∼200 km s−1, with 19 BSSs (out of 67) spinning faster than 40 km s−1. This sets the percentage of fast-rotating BSSs to ∼28%. Such a percentage is roughly comparable to that measured in other loose systems (ωCentauri, M4, and M55) and significantly larger than that measured in high-density clusters (as 47 Tucanae, NGC 6397, NGC 6752, and M30). This evidence supports a scenario where recent BSS formation (mainly from the evolution of binary systems) is occurring in low-density environments. We also find that the BSS rotational velocity tends to decrease for decreasing luminosity and surface temperature, similarly to what is observed in main-sequence stars. Hence, further investigations are needed to understand the impact of BSS internal structure on the observed rotational velocities.more » « less
