Recent studies of nearby globular clusters have discovered excess dark mass in their cores, apparently in an extended distribution, and simulations indicate that this mass is composed mostly of white dwarfs (respectively stellar-mass black holes) in clusters that are core collapsed (respectively with a flatter core). We perform mass-anisotropy modelling of the closest globular cluster, M4, with intermediate slope for the inner stellar density. We use proper motion data from Gaia Early Data Release 3 (EDR3) and from observations by the Hubble Space Telescope. We extract the mass profile employing Bayesian Jeans modelling, and check our fits with realistic mock data. Our analyses return isotropic motions in the cluster core and tangential motions (β ≈ −0.4 ± 0.1) in the outskirts. We also robustly measure a dark central mass of roughly $800\pm 300 \, \rm M_\odot$ , but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), and a dark population of stellar remnants of extent ${\approx} 0.016\, {\rm pc} \simeq 3300\, {\rm au}$ . However, when removing a high-velocity star from the cluster centre, the same mass excess is found, but more extended (${\sim} 0.034\, {\rm pc} \approx 7000\, {\rm au}$ ). We use Monte Carlo N-body models of M4 to interpret the second outcome, and find that our excess mass is not sufficiently extended to be confidently associated with a dark population of remnants. Finally, we discuss the feasibility of these two scenarios (i.e. IMBH versus remnants), and propose new observations that could help to better grasp the complex dynamics in M4’s core.
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
-
more » « less
ABSTRACT -
more » « less
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.CMC treats 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 ofCMC is pinned directly to theCOSMIC population 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 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 compact stellar remnants in the mass range 2–5
M ⊙has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5M ⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte CarloN -body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc−3yr−1. -
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.
