No Arabic abstract
We set out to compare the age-metallicity relation (AMR) of massive clusters from Magellanic Cloud mass galaxies in the E-MOSAICS suite of numerical cosmological simulations with an amalgamation of observational data of massive clusters in the Large and Small Magellanic Clouds (LMC/SMC). We aim to test if: i) star cluster formation proceeds according to universal physical processes, suggestive of a common formation mechanism for young-massive clusters (YMCs), intermediate-age clusters (IACs), and ancient globular clusters (GCs); ii) massive clusters of all ages trace a continuous AMR; iii) the AMRs of smaller mass galaxies show a shallower relation when compared to more massive galaxies. Our results show that, within the uncertainties, the predicted AMRs of L/SMC-mass galaxies with similar star formation histories to the L/SMC follow the same relation as observations. We also find that the metallicity at which the AMR saturates increases with galaxy mass, which is also found for the field star AMRs. This suggests that relatively low-metallicity clusters can still form in dwarfs galaxies. Given our results, we suggest that ancient GCs share their formation mechanism with IACs and YMCs, in which GCs are the result of a universal process of star cluster formation during the early episodes of star formation in their host galaxies.
We present 25 cosmological zoom-in simulations of Milky Way-mass galaxies in the `MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) project. E-MOSAICS couples a detailed physical model for the formation, evolution, and disruption of star clusters to the EAGLE galaxy formation simulations. This enables following the co-formation and co-evolution of galaxies and their star cluster populations, thus realising the long-standing promise of using globular clusters (GCs) as tracers of galaxy formation and assembly. The simulations show that the age-metallicity distributions of GC populations exhibit strong galaxy-to-galaxy variations, resulting from differences in their evolutionary histories. We develop a formalism for systematically constraining the assembly histories of galaxies using GC age-metallicity distributions. These distributions are characterised through 13 metrics that we correlate with 30 quantities describing galaxy formation and assembly (e.g. halo properties, formation/assembly redshifts, stellar mass assembly time-scales, galaxy merger statistics), resulting in 20 statistically (highly) significant correlations. The GC age-metallicity distribution is a sensitive probe of the mass growth, metal enrichment, and minor merger history of the host galaxy. No such relation is found between GCs and major mergers, which play a sub-dominant role in GC formation for Milky Way-mass galaxies. Finally, we show how the GC age-metallicity distribution enables the reconstruction of the host galaxys merger tree, allowing us to identify all progenitors with masses $M_*gtrsim10^8$ M$_odot$ for redshifts $1leq zleq2.5$. These results demonstrate that cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies successfully unlock the potential of GCs as quantitative tracers of galaxy formation and assembly.
Linking globular clusters (GCs) to the assembly of their host galaxies is an overarching goal in GC studies. The inference of tight scaling relations between GC system properties and the mass of both the stellar and dark halo components of their host galaxies are indicative of an intimate physical connection, yet have also raised fundamental questions about how and when GCs form. Specifically, the inferred correlation between the mass of a GC system (Mgc) and the dark matter halo mass (Mhalo) of a galaxy has been posited as a consequence of a causal relation between the formation of dark matter mini-haloes and GC formation during the early epochs of galaxy assembly. We present the first results from a new simulation of a cosmological volume ($L=34.4$~cMpc on a side) from the E-MOSAICS suite, which includes treatments of the formation and evolution of GCs within the framework of a detailed galaxy formation model. The simulated Mgc-Mhalo relation is linear for halo masses $>5times10^{11}~Msun$, and is driven by the hierarchical assembly of galaxies. Below this halo mass, the simulated relation features a downturn, which we show is consistent with observations, and is driven by the underlying stellar mass-halo mass relation of galaxies. Our fiducial model reproduces the observed Mgc-Mstar relation across the full mass range, which we argue is more physically relevant than the Mgc-Mhalo relation. We also explore the physical processes driving the observed constant value of $Mgc / Mhalo sim 5times10^{-5}$ and find that it is the result of a combination of cluster formation physics and cluster disruption.
The formation histories of globular clusters (GCs) are a key diagnostic for understanding their relation to the evolution of the Universe through cosmic time. We use the suite of 25 cosmological zoom-in simulations of present-day Milky Way-mass galaxies from the E-MOSAICS project to study the formation histories of stars, clusters, and GCs, and how these are affected by the environmental dependence of the cluster formation physics. We find that the median lookback time of GC formation in these galaxies is ${sim}10.73~$Gyr ($z=2.1$), roughly $2.5~$Gyr earlier than that of the field stars (${sim}8.34~$Gyr or $z=1.1$). The epoch of peak GC formation is mainly determined by the time evolution of the maximum cluster mass, which depends on the galactic environment and largely increases with the gas pressure. Different metallicity subpopulations of stars, clusters and GCs present overlapping formation histories, implying that star and cluster formation represent continuous processes. The metal-poor GCs ($-2.5<[rm Fe/H]<-1.5$) of our galaxies are older than the metal-rich GC subpopulation ($-1.0<[rm Fe/H]<-0.5$), forming $12.13~$Gyr and $10.15~$Gyr ago ($z=3.7$ and $z=1.8$), respectively. The median ages of GCs are found to decrease gradually with increasing metallicity, which suggests different GC metallicity subpopulations do not form independently and their spatial and kinematic distributions are the result of their evolution in the context of hierarchical galaxy formation and evolution. We predict that proto-GC formation is most prevalent at $2lesssim z lesssim 3$, which could be tested with observations of lensed galaxies using JWST.
Stellar ages are a crucial component to studying the evolution of the Milky Way. Using Gaia DR2 distance estimates, it is now possible to estimate stellar ages for a larger volume of evolved stars through isochrone matching. This work presents [M/H]-age and [$alpha$/M]-age relations derived for different spatial locations in the Milky Way disc. These relations are derived by hierarchically modelling the star formation history of stars within a given chemical abundance bin. For the first time, we directly observe that significant variation is apparent in the [M/H]-age relation as a function of both Galactocentric radius and distance from the disc mid-plane. The [M/H]-age relations support claims that radial migration has a significant effect in the plane of the disc. Using the [M/H] bin with the youngest mean age at each radial zone in the plane of the disc, the present-day metallicity gradient is measured to be $-0.059 pm 0.010$ dex kpc$^{-1}$, in agreement with Cepheids and young field stars. We find a vertically flared distribution of young stars in the outer disc, confirming predictions of models and previous observations. The mean age of the [M/H]-[$alpha$/M] distribution of the solar neighborhood suggests that the high-[M/H] stars are not an evolutionary extension of the low-$alpha$ sequence. Our observational results are important constraints to Galactic simulations and models of chemical evolution.
In the local Universe, globular clusters (GCs) with metallicities $[{rm Fe}/{rm H}]<-2.5$ are extremely rare. In this Letter, the close connection between GC formation and galaxy evolution is used to show that this GC metallicity `floor results from the galaxy mass-metallicity relation of ultra low-luminosity galaxies (ULLGs) at high redshift, where the most metal-poor GCs must have formed. Galaxies with metallicities $[{rm Fe}/{rm H}]lesssim-2.5$ have too low masses to form GCs with initial masses $M_{rm i}gtrsim10^5~{rm M}_odot$, needed to survive for a Hubble time. This translates the galaxy mass-metallicity relation into a maximum initial cluster mass-metallicity relation for $[{rm Fe}/{rm H}]lesssim-1.8$, which naturally leads to the observed colour-magnitude relation of metal-poor GCs at $z=0$ (the `blue tilt). Its strength traces the slope of the gas phase mass-metallicity relation of ULLGs. Based on the observed blue tilt of GCs in the Virgo and Fornax Clusters, the galaxy mass-metallicity relation is predicted to have a slope of $alpha=0.4pm0.1$ for $10^5lesssim M_star/{rm M}_odotlesssim10^7$ at $zgtrsim2$. The GC metallicity floor implies a minimum host galaxy mass and a maximum redshift for GC formation. Any proto-GCs that may be detected at $z>9$ are most likely to end up in galaxies presently more massive than the Milky Way, whereas GCs in low-mass galaxies such as the Fornax dSph ($M_starapprox4times10^7~{rm M}_odot$) formed at $zlesssim3$.