We introduce the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) project. E-MOSAICS incorporates models describing the formation, evolution and disruption of star clusters into the EAGLE galaxy formation simulations, enabling the examination of the co-evolution of star clusters and their host galaxies in a fully cosmological context. A fraction of the star formation rate of dense gas is assumed to yield a cluster population; this fraction, and the population's initial properties, are governed by the physical properties of the natal gas. The subsequent evolution and disruption of the entire cluster population is followed accounting for two-body relaxation, stellar evolution, and gravitational shocks induced by the local tidal field. This introductory paper presents a detailed description of the model and initial results from a suite of 10 simulations of ∼ L galaxies with disc-like morphologies at z = 0. The simulations broadly reproduce key observed characteristics of young star clusters and globular clusters (GCs), without invoking separate formation mechanisms for each population. The simulated GCs are the surviving population of massive clusters formed at early epochs (z 1 − 2), when the characteristic pressures and surface densities of star-forming gas were significantly higher than observed in local galaxies. We examine the influence of the star formation and assembly histories of galaxies on their cluster populations, finding that (at similar present-day mass) earlier-forming galaxies foster a more massive and disruption-resilient cluster population, while galaxies with late mergers are capable of forming massive clusters even at late cosmic epochs. We find that the phenomenological treatment of interstellar gas in EAGLE precludes the accurate modelling of cluster disruption in low-density environments, but infer that simulations incorporating an explicitly-modelled cold interstellar gas phase will overcome this shortcoming.

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.

Recent work indicates that the nearby Galactic halo is dominated by the debris from a major accretion event. We confirm that result from an analysis of APOGEE-DR14 element abundances and Gaia-DR2 kinematics of halo stars. We show that ∼ 2/3 of nearby halo stars have high orbital eccentricities (e 0.8), and abundance patterns typical of massive Milky Way dwarf galaxy satellites today, characterised by relatively low [Fe/H], [Mg/Fe], [Al/Fe], and [Ni/Fe]. The trend followed by high e stars in the [Mg/Fe]-[Fe/H] plane shows a change of slope at [Fe/H]∼ −1.3, which is also typical of stellar populations from relatively massive dwarf galaxies. Low e stars exhibit no such change of slope within the observed [Fe/H] range and show slightly higher abundances of Mg, Al and Ni. Unlike their low e counterparts, high e stars show slightly retrograde motion, make higher vertical excursions and reach larger apocentre radii. By comparing the position in [Mg/Fe]-[Fe/H] space of high e stars with those of accreted galaxies from the EAGLE suite of cosmological simulations we constrain the mass of the accreted satellite to be in the range 10 8.5 M * 10 9 M . We show that the median orbital eccentricities of debris are largely unchanged since merger time, implying that this accretion event likely happened at z 1.5. The exact nature of the low e population is unclear, but we hypothesise that it is a combination of in situ star formation, high |z| disc stars, lower mass accretion events, and contamination by the low e tail of the high e population. Finally, our results imply that the accretion history of the Milky Way was quite unusual.

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 timescales, 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 galaxy's merger tree, allowing us to identify all progenitors with masses M * 10 8 M for redshifts 1 z 2.5. These results demonstrate that cosmological simulations of the coformation and co-evolution of GCs and their host galaxies successfully unlock the potential of GCs as quantitative tracers of galaxy formation and assembly.

We report evidence from APOGEE for the presence of a new metal-poor stellar structure located within ∼4 kpc of the Galactic Centre. Characterized by a chemical composition resembling those of low-mass satellites of the Milky Way, this new inner Galaxy structure (IGS) seems to be chemically and dynamically detached from more metal-rich populations in the inner Galaxy. We conjecture that this structure is associated with an accretion event that likely occurred in the early life of the Milky Way. Comparing the mean elemental abundances of this structure with predictions from cosmological numerical simulations, we estimate that the progenitor system had a stellar mass of ∼5 × 108 M⊙, or approximately twice the mass of the recently discovered Gaia-Enceladus/Sausage system. We find that the accreted:in situ ratio within our metal-poor ([Fe/H] < –0.8) bulge sample is somewhere between 1:3 and 1:2, confirming predictions of cosmological numerical simulations by various groups.