skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Great balls of FIRE II: The evolution and destruction of star clusters across cosmic time in a Milky Way-mass galaxy
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
Award ID(s):
2108230 1715216 1652522 2045928 2310362
PAR ID:
10428539
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
521
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
124 to 147
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The properties of globular clusters (GCs) contain valuable information of their host galaxies and dark-matter halos. In the remarkable example of ultra-diffuse galaxy, NGC5846-UDG1, the GC population exhibits strong radial mass segregation, indicative of dynamical-friction-driven orbital decay, which opens the possibility of using imaging data alone to constrain the dark-matter content of the galaxy. To explore this possibility, we develop a semianalytical model of GC evolution, which starts from the initial mass, structural, and spatial distributions of the GC progenitors, and follows the effects of dynamical friction, tidal evolution, and two-body relaxation. Using Markov Chain Monte Carlo, we forward-model the GCs in a UDG1-like potential to match the observed GC statistics, and to constrain the profile of the host halo and the origin of the GCs. We find that, with the assumptions of zero mass segregation when the star clusters were born, UDG1 is relatively dark-matter-poor compared to what is expected from stellar-to-halo–mass relations, and its halo concentration is lower than the cosmological average, irrespective of having a cuspy or a cored profile. Its GC population has an initial spatial distribution more extended than the smooth stellar distribution. We discuss the results in the context of scaling laws of galaxy–halo connections, and warn against naively using the GC-abundance–halo–mass relation to infer the halo mass of ultra-diffuse galaxies. Our model is generally applicable to GC-rich dwarf galaxies, and is publicly available. 
    more » « less
  2. ABSTRACT The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest haloes as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyse the catalogues of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass haloes with M200 ≲ 1011.2 M⊙. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive haloes is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intracluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than 10 GC tracers are particularly prone to severe uncertainties. Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive haloes suggested by GC kinematics around some UDGs. 
    more » « less
  3. Abstract In a series of three papers, we introduced a novel cluster formation model that describes the formation, growth, and disruption of star clusters in high-resolution cosmological simulations. We tested this model on a Milky Way-sized galaxy and found that various properties of young massive clusters, such as the mass function and formation efficiency, are consistent with observations in the local universe. Interestingly, most massive clusters – globular cluster candidates – are preferentially formed during major merger events. We follow the dynamical evolution of clusters in the galactic tidal field. Due to tidal disruption, the cluster mass function evolves from initial power law to a peaked shape. The surviving clusters at z = 0 show a broad range of metallicity [Fe/H] from -3 to -0.5. A robust prediction of the model is the age–metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3 Gyr. 
    more » « less
  4. Abstract Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we createN-body models that reproduce the ratio of single-to-binary pulsars in Milky Way–like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play dominant roles: (1) tidal disruption of main-sequence stars by neutron stars and (2) gravitational collapse of heavy white dwarf binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC. 
    more » « less
  5. Some ultra diffuse galaxies (UDGs) reveal many more globular clusters (GCs) than classical dwarf galaxies of the same stellar mass. These UDGs, with a mass in their GC system (M_GC) approaching 10 per cent of their host galaxy stellar mass (M_*), are also inferred to have high halo mass to stellar mass ratios (M_halo/M_*). They have been dubbed Failed Galaxies. It is unknown what role high GC formation efficiencies and/or low destruction rates play in determining the high M_GC/M_* ratios of some UDGs. Here we present a simple model, which is informed by recent JWST observations of lensed galaxies and by a simulation in the literature of GC mass loss and tidal disruption in dwarf galaxies. With this simple model, we aim to constrain the effects of GC efficiency/destruction on the observed GC richness of UDGs and their variation with the integrated stellar populations of UDGs. We assume no ongoing star formation (i.e. quenching at early times) and that the disrupted GCs contribute their stars to those of the host galaxy. We find that UDGs, with high M_GC/M_* ratios today, are most likely the result of very high GC formation efficiencies combined with modest rates of GC destruction. The current data loosely follow the model that ranges from the mean stellar population of classical dwarfs to that of metal-poor GCs as M_GC/M_* increases. As more data becomes available for UDGs, our simple model can be refined and tested further. 
    more » « less