skip to main content


Title: Formation and evolution of globular clusters in cosmological simulations
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
Award ID(s):
1909063
NSF-PAR ID:
10166674
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Proceedings of the International Astronomical Union
Volume:
14
Issue:
S351
ISSN:
1743-9213
Page Range / eLocation ID:
34 to 39
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We present a suite of galaxy formation simulations that directly model star cluster formation and disruption. Starting from a model previously developed by our group, here we introduce several improvements to the prescriptions for cluster formation and feedback, then test these updates using a large suite of cosmological simulations of Milky Way mass galaxies. We perform a differential analysis with the goal of understanding how each of the updates affects star cluster populations. Two key parameters are the momentum boost of supernova feedback fboost and star formation efficiency per free-fall time ϵff. We find that fboost has a strong influence on the galactic star formation rate, with higher values leading to less star formation. The efficiency ϵff does not have a significant impact on the global star formation rate, but dramatically changes cluster properties, with increasing ϵff leading to a higher maximum cluster mass, shorter age spread of stars within clusters, and higher integrated star formation efficiencies. We also explore the redshift evolution of the observable cluster mass function, finding that most massive clusters have formed at high redshift z > 4. Extrapolation of cluster disruption to z = 0 produces good agreement with both the Galactic globular cluster mass function and age–metallicity relation. Our results emphasize the importance of using small-scale properties of galaxies to calibrate subgrid models of star cluster formation and feedback.

     
    more » « less
  2. Abstract The Fornax dwarf spheroidal galaxy has an anomalous number of globular clusters, five, for its stellar mass. There is a longstanding debate about a potential sixth globular cluster (Fornax 6) that has recently been “rediscovered” in DECam imaging. We present new Magellan/M2FS spectroscopy of the Fornax 6 cluster and Fornax dSph. Combined with literature data we identify ∼15–17 members of the Fornax 6 cluster, showing that this overdensity is indeed a star cluster and associated with the Fornax dSph. The cluster is significantly more metal-rich (mean metallicity of [ Fe / H ] ¯ = −0.71 ± 0.05) than the other five Fornax globular clusters (−2.5 < [Fe/H] < −1.4) and more metal-rich than the bulk of Fornax. We measure a velocity dispersion of 5.6 − 1.6 + 2.0 km s − 1 corresponding to an anomalously high mass-to-light of 15 < M / L < 258 at 90% confidence when calculated assuming equilibrium. Two stars inflate this dispersion and may be either Fornax field stars or as yet unresolved binary stars. Alternatively, the Fornax 6 cluster may be undergoing tidal disruption. Based on its metal-rich nature, the Fornax 6 cluster is likely younger than the other Fornax clusters, with an estimated age of ∼2 Gyr when compared to stellar isochrones. The chemodynamics and star formation history of Fornax shows imprints of major events such as infall into the Milky Way, multiple pericenter passages, star formation bursts, and/or potential mergers or interactions. Any of these events may have triggered the formation of the Fornax 6 cluster. 
    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

    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
  5. 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