More Like this

1. An elusive dark central mass in the globular cluster M4

   https://doi.org/10.1093/mnras/stad1068  

   Vitral, Eduardo ; Libralato, Mattia ; Kremer, Kyle ; Mamon, Gary A. ; Bellini, Andrea ; Bedin, Luigi R. ; Anderson, Jay ( May 2023 , Monthly Notices of the Royal Astronomical Society)

   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}$ ).more »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.

   « less
2. Great balls of FIRE II: The evolution and destruction of star clusters across cosmic time in a Milky Way-mass galaxy

   https://doi.org/10.1093/mnras/stad578  

   Rodriguez, Carl L. ; Hafen, Zachary ; Grudić, Michael Y. ; Lamberts, Astrid ; Sharma, Kuldeep ; Faucher-Giguère, Claude-André ; Wetzel, Andrew ( February 2023 , Monthly Notices of the Royal Astronomical Society)

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

   « less
3. Modeling Dense Star Clusters in the Milky Way and beyond with the Cluster Monte Carlo Code

   https://doi.org/10.3847/1538-4365/ac2edf  

   Rodriguez, Carl L. ; Weatherford, Newlin C. ; Coughlin, Scott C. ; Amaro-Seoane, Pau ; Breivik, Katelyn ; Chatterjee, Sourav ; Fragione, Giacomo ; Kıroğlu, Fulya ; Kremer, Kyle ; Rui, Nicholas Z. ; et al ( January 2022 , The Astrophysical Journal Supplement Series)

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

   « less
4. Compact Object Modeling in the Globular Cluster 47 Tucanae

   https://doi.org/10.3847/1538-4357/ac5b0b  

   Ye, Claire S. ; Kremer, Kyle ; Rodriguez, Carl L. ; Rui, Nicholas Z. ; Weatherford, Newlin C. ; Chatterjee, Sourav ; Fragione, Giacomo ; Rasio, Frederic A. ( May 2022 , The Astrophysical Journal)

   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
5. Stellar Escape from Globular Clusters. I. Escape Mechanisms and Properties at Ejection

   https://doi.org/10.3847/1538-4357/acbcc1  

   Weatherford, Newlin C. ; Kıroğlu, Fulya ; Fragione, Giacomo ; Chatterjee, Sourav ; Kremer, Kyle ; Rasio, Frederic A. ( April 2023 , The Astrophysical Journal)

   The theory of stellar escape from globular clusters (GCs) dates back nearly a century, especially the gradual evaporation of GCs via two-body relaxation coupled with external tides. More violent ejection can also occur via strong gravitational scattering, supernovae, gravitational wave-driven mergers, tidal disruption events, and physical collisions, but comprehensive study of the many escape mechanisms has been limited. Recent exquisite kinematic data from the Gaia space telescope has revealed numerous stellar streams in the Milky Way (MW) and traced the origin of many to specific MWGCs, highlighting the need for further examination of stellar escape from these clusters. In this study, the first of a series, we lay the groundwork for detailed follow-up comparisons between Cluster Monte Carlo GC models and the latest Gaia data on the outskirts of MWGCs, their tidal tails, and associated streams. We thoroughly review escape mechanisms from GCs and examine their relative contributions to the escape rate, ejection velocities, and escaper demographics. We show for the first time that three-body binary formation may dominate high-speed ejection from typical MWGCs, potentially explaining some of the hypervelocity stars in the MW. Due to their mass, black holes strongly catalyze this process, and their loss at themore »onset of observable core collapse, characterized by a steep central brightness profile, dramatically curtails three-body binary formation, despite the increased post-collapse density. We also demonstrate that even when born from a thermal eccentricity distribution, escaping binaries have significantly nonthermal eccentricities consistent with the roughly uniform distribution observed in the Galactic field.

   « less