skip to main content


Title: Tidal disruption of star clusters in galaxy formation simulations
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
Award ID(s):
1909063
NSF-PAR ID:
10369110
Author(s) / Creator(s):
;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
515
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
p. 1065-1077
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We study environmental quenching using the spatial distribution of current star formation and stellar population ages with the full SAMI Galaxy Survey. By using a star formation concentration index [C-index, defined as log10(r50, H α/r50, cont)], we separate our sample into regular galaxies (C-index ≥−0.2) and galaxies with centrally concentrated star formation (SF-concentrated; C-index <−0.2). Concentrated star formation is a potential indicator of galaxies currently undergoing ‘outside-in’ quenching. Our environments cover ungrouped galaxies, low-mass groups (M200 ≤ 1012.5M⊙), high-mass groups (M200 in the range 1012.5–14 M⊙) and clusters (M200 > 1014M⊙). We find the fraction of SF-concentrated galaxies increases as halo mass increases by 9 ± 2 per cent, 8 ± 3 per cent, 19 ± 4 per cent, and 29 ± 4 per cent for ungrouped galaxies, low-mass groups, high-mass groups, and clusters, respectively. We interpret these results as evidence for ‘outside-in’ quenching in groups and clusters. To investigate the quenching time-scale in SF-concentrated galaxies, we calculate light-weighted age (AgeL) and mass-weighted age (AgeM) using full spectral fitting, as well as the Dn4000 and HδA indices. We assume that the average galaxy age radial profile before entering a group or cluster is similar to ungrouped regular galaxies. At large radius (1–2 Re), SF-concentrated galaxies in high-mass groups have older ages than ungrouped regular galaxies with an age difference of 1.83 ± 0.38 Gyr for AgeL and 1.34 ± 0.56 Gyr for AgeM. This suggests that while ‘outside-in’ quenching can be effective in groups, the process will not quickly quench the entire galaxy. In contrast, the ages at 1–2 Re of cluster SF-concentrated galaxies and ungrouped regular galaxies are consistent (difference of 0.19 ± 0.21 Gyr for AgeL, 0.40 ± 0.61 Gyr for AgeM), suggesting the quenching process must be rapid.

     
    more » « less
  2. 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
  3. ABSTRACT

    The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modelling cluster formation in galaxy simulations: mapping giant molecular clouds (GMCs) formed self-consistently in a FIRE-2 magnetohydrodynamic galaxy simulation on to a cluster population according to a GMC-scale cluster formation model calibrated to higher resolution simulations, obtaining detailed properties of the galaxy’s star clusters in mass, metallicity, space, and time. We find $\sim 10{{\ \rm per\ cent}}$ of all stars formed in the galaxy originate in gravitationally bound clusters overall, and this fraction increases in regions with elevated Σgas and ΣSFR, because such regions host denser GMCs with higher star formation efficiency. These quantities vary systematically over the history of the galaxy, driving variations in cluster formation. The mass function of bound clusters varies – no single Schechter-like or power-law distribution applies at all times. In the most extreme episodes, clusters as massive as 7 × 106 M⊙ form in massive, dense clouds with high star formation efficiency. The initial mass–radius relation of young star clusters is consistent with an environmentally dependent 3D density that increases with Σgas and ΣSFR. The model does not reproduce the age and metallicity statistics of old ($\gt 11\rm Gyr$) globular clusters found in the Milky Way, possibly because it forms stars more slowly at z > 3.

     
    more » « less
  4. ABSTRACT

    Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the starforge simulation suite – a set of radiation-magnetohydrodynamic simulations of star-forming turbulent clouds that include all key stellar feedback processes inherent to star formation. We use nbody6++gpu to follow the evolution of the young stellar systems after gas removal. We use HDBSCAN to define stellar groups and analyse the stellar kinematics to identify the true bound star clusters. The conditions modeled by the simulations, i.e. global cloud surface densities below 0.15 g cm−2, star formation efficiencies below 15 per cent, and gas expulsion time-scales shorter than a free fall time, primarily produce expanding stellar associations and small clusters. The largest star clusters, which have ∼1000 bound members, form in the densest and lowest velocity dispersion clouds, representing ∼32 and 39 per cent of the stars in the simulations, respectively. The cloud’s early dynamical state plays a significant role in setting the classical star formation efficiency versus bound fraction relation. All stellar groups follow a narrow mass-velocity dispersion power-law relation at 10 Myr with a power-law index of 0.21. This correlation result in a distinct mass–size relationship for bound clusters. We also provide valuable constraints on the gas dispersal time-scale during the star formation process and analyse the implications for the formation of bound systems.

     
    more » « less
  5. ABSTRACT

    At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star-forming main sequence, ΔSFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios μ ≥ 0.25), minor (0.1 ≤ μ < 0.25), and mini mergers (0.01 ≤ μ < 0.1) in driving this connection amongst star-forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated IllustrisTNG (TNG) galaxies over 0.1 ≤ z ≤ 0.7 with log (M⋆/M⊙) ≥ 9 (∼750 k images). Using their instantaneous star formation rates (SFRs), known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between ΔSFMS and asymmetry and (2) a strikingly similar trend emerges between ΔSFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long time-scales (at least ∼3 Gyr after coalescence, on average) – in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only ∼1 Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by μ ≥ 0.01 mergers since z = 0.7 in TNG50 SFGs, we show that mini mergers are responsible for (i) 55 per cent of all merger-driven star formation and (ii) 70 per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost time-scales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs.

     
    more » « less