skip to main content


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
NSF-PAR ID:
10400092
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
521
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
p. 124-147
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. ABSTRACT

    Galaxy mergers are known to host abundant young massive cluster (YMC) populations, whose formation mechanism is still not well-understood. Here, we present a high-resolution galaxy merger simulation with explicit star formation and stellar feedback prescriptions to investigate how mergers affect the properties of the interstellar medium and YMCs. Compared with a controlled simulation of an isolated galaxy, the mass fraction of dense and high-pressure gas is much higher in mergers. Consequently, the mass function of both molecular clouds and YMCs becomes shallower and extends to higher masses. Moreover, cluster formation efficiency is significantly enhanced and correlates positively with the star formation rate surface density and gas pressure. We track the orbits of YMCs and investigate the time evolution of tidal fields during the course of the merger. At an early stage of the merger, the tidal field strength correlates positively with YMC mass, λtid ∝ M0.71, which systematically affects the shape of the mass function and age distribution of the YMCs. At later times, most YMCs closely follow the orbits of their host galaxies, gradually sinking into the centre of the merger remnant due to dynamical friction, and are quickly dissolved via efficient tidal disruption. Interestingly, YMCs formed during the first passage, mostly in tidal tails and bridges, are distributed over a wide range of galactocentric radii, greatly increasing their survivability because of the much weaker tidal field in the outskirts of the merger system. These YMCs are promising candidates for globular clusters that survive to the present day.

     
    more » « less
  5. ABSTRACT

    We investigate the evolution of the tidal field experienced by massive star clusters using cosmological simulations of Milky Way-sized galaxies. Clusters in our simulations experience the strongest tidal force in the first few hundred Myr after formation, when the maximum eigenvalue of the tidal tensor reaches several times 104 Gyr−2. After about 1 Gyr the tidal field plateaus at a lower value, with the median λm ∼ 3 × 103 Gyr−2. The fraction of time clusters spend in high tidal strength (λm > 3 × 104 Gyr−2) regions also decreases with their age from ∼20 per cent immediately after formation to less than 1 per cent after 1 Gyr. At early ages both the in situ and ex situ clusters experience similar tidal fields, while at older ages the in situ clusters in general experience stronger tidal field due to their lower orbits in host galaxy. This difference is reflected in the survival of clusters: we looked into cluster disruption calculated in simulation runtime and found that ex situ star clusters of the same initial mass typically end up with higher bound fraction at the last available simulation snapshot than the in situ ones.

     
    more » « less