We use the age-metallicity distribution of 96 Galactic globular clusters (GCs) to infer the formation and assembly history of the Milky Way (MW), culminating in the reconstruction of its merger tree. Based on a quantitative comparison of the Galactic GC population to the 25 cosmological zoom-in simulations of MW-mass galaxies in the E-MOSAICS project, which self-consistently model the formation and evolution of GC populations in a cosmological context, we find that the MW assembled quickly for its mass, reaching {25, 50}% of its present-day halo mass already at z = {3, 1.5} and half of its present-day stellar mass at z = 1.2. We reconstruct the MW's merger tree from its GC age-metallicity distribution, inferring the number of mergers as a function of mass ratio and redshift. These statistics place the MW's assembly rate among the 72th-94th percentile of the E-MOSAICS galaxies, whereas its integrated properties (e.g. number of mergers, halo concentration) match the median of the simulations. We conclude that the MW has experienced no major mergers (mass ratios >1:4) since z ∼ 4, sharpening previous limits of z ∼ 2. We identify three massive satellite progenitors and constrain their mass growth and enrichment histories. Two are proposed to correspond to Sagittarius (few 10 8 M ) and the GCs formerly associated with Canis Major (∼ 10 9 M ). The third satellite has no known associated relic and was likely accreted between z = 0.6-1.3. We name this enigmatic galaxy Kraken and propose that it is the most massive satellite (M * ∼ 2 × 10 9 M ) ever accreted by the MW. We predict that ∼ 40% of the Galactic GCs formed ex-situ (in galaxies with masses M * = 2 × 10 7 -2 × 10 9 M ), with 6 ± 1 being former nuclear clusters.

Globular clusters (GCs) formed when the Milky Way experienced a phase of rapid assembly. We use the wealth of information contained in the Galactic GC population to quantify the properties of the satellite galaxies from which the Milky Way assembled. To achieve this, we train an artificial neural network on the E-MOSAICS cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies. The network uses the ages, metallicities, and orbital properties of GCs that formed in the same progenitor galaxies to predict the stellar masses and accretion redshifts of these progenitors. We apply the network to Galactic GCs associated with five progenitors: Gaia-Enceladus, the Helmi streams, Sequoia, Sagittarius, and the recently discovered, ‘low-energy’ GCs, which provide an excellent match to the predicted properties of the enigmatic galaxy ‘Kraken’. The five galaxies cover a narrow stellar mass range [M⋆ = (0.6 − 4.6) × 108 M⊙], but have widely different accretion redshifts (${z_{\rm acc}}=0.57{-}2.65$). All accretion events represent minor mergers, but Kraken likely represents the most major merger ever experienced by the Milky Way, with stellar and virial mass ratios of ${r_{M_\star }}=1$:$31^{+34}_{-16}$ and ${r_{M_{\rm h}}}=1$:$7^{+4}_{-2}$, respectively. The progenitors match the z = 0 relation between GC number and halo virial mass, but have elevated specific frequencies, suggesting an evolution with redshift. Even though these progenitors likely were the Milky Way’s most massive accretion events, they contributed a total mass of only log (M⋆, tot/M⊙) = 9.0 ± 0.1, similar to the stellar halo. This implies that the Milky Way grew its stellar mass mostly by in-situ star formation. We conclude by organising these accretion events into the most detailed reconstruction to date of the Milky Way’s merger tree.

It has recently been established that the properties of young star clusters (YSCs) can vary as a function of the galactic environment in which they are found. We use the cluster catalogue produced by the Legacy Extragalactic UV Survey (LEGUS) collaboration to investigate cluster properties in the spiral galaxy M51. We analyse the cluster population as a function of galactocentric distance and in arm and interarm regions. The cluster mass function exhibits a similar shape at all radial bins, described by a power law with a slope close to −2 and an exponential truncation around 10 5 M . While the mass functions of the YSCs in the spiral arm and interarm regions have similar truncation masses, the inter-arm region mass function has a significantly steeper slope than the one in the arm region; a trend that is also observed in the giant molecular cloud mass function and predicted by simulations. The age distribution of clusters is dependent on the region considered, and is consistent with rapid disruption only in dense regions, while little disruption is observed at large galactocentric distances and in the inter-arm region. The fraction of stars forming in clusters does not show radial variations, despite the drop in the H 2 surface density measured as function of galactocentric distance. We suggest that the higher disruption rate observed in the inner part of the galaxy is likely at the origin of the observed flat cluster formation efficiency radial profile.

No abstract

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 ∼10.73 Gyr (z = 2.1), roughly 2.5 Gyr earlier than that of the field stars (∼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 < [Fe/H] < −1.5) of our galaxies are older than the metal-rich GC subpopulation (−1.0 < [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 2 z 3, which could be tested with observations of lensed galaxies using JWST.