Modern observational and analytic techniques now enable the direct measurement of star formation histories and the inference of galaxy assembly histories. However, current theoretical predictions of assembly are not ideally suited for direct comparison with such observational data. We therefore extend the work of prior examinations of the contribution of ex-situ stars to the stellar mass budget of simulated galaxies. Our predictions are specifically tailored for direct testing with a new generation of observational techniques by calculating ex-situ fractions as functions of galaxy mass and morphological type, for a range of surface brightnesses. These enable comparison with results from large FoV IFU spectrographs, and increasingly accurate spectral fitting, providing a look-up method for the estimated accreted fraction. We furthermore provide predictions of ex-situ mass fractions as functions of galaxy mass, galactocentric radius and environment. Using $z=0$ snapshots from the 100cMpc$^3$ and 25cMpc$^3$ EAGLE simulations we corroborate the findings of prior studies, finding that ex-situ fraction increases with stellar mass for central and satellite galaxies in a stellar mass range of 2$\times$10$^{7}$ - 1.9$\times$10$^{12}$ M$_{\odot}$. For those galaxies of mass M$_*$>5$\times$10$^{8}$M$_{\odot}$, we find that the total ex-situ mass fraction is greater for more extended galaxies at fixed mass. When categorising satellite galaxies by their parent group/cluster halo mass we find that the ex-situ fraction decreases with increasing parent halo mass at fixed galaxy mass. This apparently counter-intuitive result may be the result of high passing velocities within large cluster halos inhibiting efficient accretion onto individual galaxies.

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 (M_gc) and the dark matter halo mass (M_halo) 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.4cMpc 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 M_gc-M_halo relation is linear for halo masses $>5\times10^{11}$ Msun, and is driven by the hierarchical assembly of galaxies, in agreement with previous studies. 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 (M_star)-halo mass relation of galaxies. Our fiducial model reproduces the observed M_gc-M_star relation across the full mass range, which we argue is more physically relevant than the M_gc-M_halo relation. We also explore the physical processes driving the observed constant value of M_gc / M_halo $\sim 5\times10^{-5}$ and find that it is the result of a combination of cluster formation physics and cluster disruption.

Globular clusters (GCs) are typically old, with most having formed at z >~ 2. This makes understanding their birth environments difficult, as they are typically too distant to observe with sufficient angular resolution to resolve GC birth sites. Using 25 cosmological zoom-in simulations of Milky Way-like galaxies from the E-MOSAICS project, with physically-motivated models for star formation, feedback, and the formation, evolution, and disruption of GCs, we identify the birth environments of present-day GCs. We find roughly half of GCs in these galaxies formed in-situ (52.0 +/- 1.0 per cent) between z ~ 2 - 4, in turbulent, high-pressure discs fed by gas that was accreted without ever being strongly heated through a virial shock or feedback. A minority of GCs form during mergers (12.6 +/- 0.6 per cent in major mergers, and 7.2 +/- 0.5 per cent in minor mergers), but we find that mergers are important for preserving the GCs seen today by ejecting them from their natal, high density interstellar medium (ISM), where proto-GCs are rapidly destroyed due to tidal shocks from ISM substructure. This chaotic history of hierarchical galaxy assembly acts to mix the spatial and kinematic distribution of GCs formed through different channels, making it difficult to use observable GC properties to distinguish GCs formed in mergers from ones formed by smooth accretion, and similarly GCs formed in-situ from those formed ex-situ. These results suggest a simple picture of GC formation, in which GCs are a natural outcome of normal star formation in the typical, gas-rich galaxies that are the progenitors of present-day galaxies.

We present a detailed comparison of the Milky Way (MW) globular cluster (GC) kinematics with the 25 Milky Way-mass cosmological simulations from the E-MOSAICS project. While the MW falls within the kinematic distribution of GCs spanned by the simulations, the relative kinematics of its metal-rich ($[\rm{Fe/H}]>-1.2$) versus metal-poor ($[\rm{Fe/H}]<-1.2$), and inner ($r<8\rm{kpc}$) versus outer ($r>8\rm{kpc}$) populations are atypical for its mass. To understand the origins of these features, we perform a comprehensive statistical analysis of the simulations, and find 18 correlations describing the assembly of $L^*$ galaxies and their dark matter haloes based on their GC population kinematics. The correlations arise because the orbital distributions of accreted and in-situ GCs depend on the masses and accretion redshifts of accreted satellites, driven by the combined effects of dynamical fraction, tidal stripping, and dynamical heating. Because the kinematics of in-situ/accreted GCs are broadly traced by the metal-rich/metal-poor and inner/outer populations, the observed GC kinematics are a sensitive probe of galaxy assembly. We predict that relative to the population of $L^*$ galaxies, the MW assembled its dark matter and stellar mass rapidly through a combination of in-situ star formation, more than a dozen low-mass mergers, and $1.4\pm1.2$ early ($z=3.1\pm1.3$) major merger. The rapid assembly period ended early, limiting the fraction of accreted stars. We conclude by providing detailed quantitative predictions for the assembly history of the MW.

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: {\it 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_\star=(0.6{-}4.6)\times10^8~{\rm M}_\odot$], 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_{\rm \star,tot}/{\rm M}_\odot)}=9.0\pm0.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.

The ages and metallicities of globular clusters (GCs) are known to be powerful tracers of the properties of their progenitor galaxies, enabling their use in determining the merger histories of galaxies. However, while useful in separating GCs into individual accretion events, the orbits of GC groups themselves have received less attention as probes of their progenitor galaxy properties. In this work, we use simulations of galaxies and their GC systems from the E-MOSAICS project to explore how the present-day orbital properties of GCs are related to the properties of their progenitor galaxies. We find that the orbits of GCs deposited by accretion events are sensitive to the mass and merger redshift of the satellite galaxy. Earlier mergers and larger galaxy masses deposit GCs at smaller median apocentres and lower total orbital energy. The orbital properties of accreted groups of GCs can therefore be used to infer the properties of their progenitor galaxy, though there exists a degeneracy between galaxy mass and accretion time. Combining GC orbits with other tracers (GC ages, metallicities) will help to break the galaxy mass/accretion time degeneracy, enabling stronger constraints on the properties of their progenitor galaxy. In situ GCs generally orbit at lower energies (small apocentres) than accreted GCs, however they exhibit a large tail to high energies and even retrograde orbits (relative to the present-day disc), showing significant overlap with accreted GCs. Applying the results to Milky Way GCs groups suggests a merger redshift $z \sim 1.5$ for the Gaia Sausage/Enceladus and $z>2$ for the `low-energy'/Kraken group, adding further evidence that the Milky Way had two significant mergers in its past.

We summarise recent results from analysis of APOGEE/Gaia data for stellar populations in the Galactic halo, disk, and bulge, leading to constraints on the contribution of dwarf galaxies and globular clusters to the stellar content of the Milky Way halo. Interpretation of the extant data in light of cosmological numerical simulations suggests that the Milky Way has been subject to an unusually intense accretion history at z >~ 1.5.

The {\alpha}-element abundances of the globular cluster (GC) and field star populations of galaxies encode information about the formation of each of these components. We use the E-MOSAICS cosmological simulations of ~L* galaxies and their GCs to investigate the [{\alpha}/Fe]-[Fe/H] distribution of field stars and GCs in 25 Milky Way-mass galaxies. The [{\alpha}/Fe]-[Fe/H] distribution go GCs largely follows that of the field stars and can also therefore be used as tracers of the [{\alpha}/Fe]-[Fe/H] evolution of the galaxy. Due to the difference in their star formation histories, GCs associated with stellar streams (i.e. which have recently been accreted) have systematically lower [{\alpha}/Fe] at fixed [Fe/H]. Therefore, if a GC is observed to have low [{\alpha}/Fe] for its [Fe/H] there is an increased probability that this GC was accreted recently alongside a dwarf galaxy. There is a wide range of shapes for the field star [{\alpha}/Fe]-[Fe/H] distribution, with a notable subset of galaxies exhibiting bimodal distributions, in which the high [{\alpha}/Fe] sequence is mostly comprised of stars in the bulge, a high fraction of which are from disrupted GCs. We calculate the contribution of disrupted GCs to the bulge component of the 25 simulated galaxies and find values between 0.3-14 per cent, where this fraction correlates with the galaxy's formation time. The upper range of these fractions is compatible with observationally-inferred measurements for the Milky Way, suggesting that in this respect the Milky Way is not typical of L* galaxies, having experienced a phase of unusually rapid growth at early times.

Globular clusters (GCs) have been posited, alongside dwarf galaxies, as significant contributors to the field stellar population of the Galactic halo. In order to quantify their contribution, we examine the fraction of halo stars formed in stellar clusters in the suite of 25 present-day Milky Way-mass cosmological zoom simulations from the E-MOSAICS project. We find that a median of $2.3$ and $0.3$ per cent of the mass in halo field stars formed in clusters and GCs, defined as clusters more massive than $5\times 10^3$ and $10^5~M_{\odot}$, respectively, with the $25$-$75$th percentiles spanning $1.9$-$3.0$ and $0.2$-$0.5$ per cent being caused by differences in the assembly histories of the host galaxies. Under the extreme assumption that no stellar cluster survives to the present day, the mass fractions increase to a median of $5.9$ and $1.8$ per cent. These small fractions indicate that the disruption of GCs plays a sub-dominant role in the build-up of the stellar halo. We also determine the contributed halo mass fraction that would present signatures of light-element abundance variations considered to be unique to GCs, and find that clusters and GCs would contribute a median of $1.1$ and $0.2$ per cent, respectively. We estimate the contributed fraction of GC stars to the Milky Way halo, based on recent surveys, and find upper limits of $2$-$5$ per cent (significantly lower than previous estimates), suggesting that models other than those invoking strong mass loss are required to describe the formation of chemically enriched stellar populations in GCs.

We present an analysis of young star clusters (YSCs) that form in the E-MOSAICS cosmological, hydrodynamical simulations of galaxies and their star cluster populations. Through comparisons with observed YSC populations, this work aims to test models for YSC formation and obtain an insight into the formation processes at work in part of the local galaxy population. We find that the models used in E-MOSAICS for the cluster formation efficiency and high-mass truncation of the initial cluster mass function ($M_\mathrm{c,\ast}$) both quantitatively reproduce the observed values of cluster populations in nearby galaxies. At higher redshifts ($z \geq 2$, near the peak of globular cluster formation) we find that, at a constant star formation rate (SFR) surface density, $M_\mathrm{c,\ast}$ is larger than at $z=0$ by a factor of four due to the higher gas fractions in the simulated high-redshift galaxies. Similar processes should be at work in local galaxies, offering a new way to test the models. We find that cluster age distributions may be sensitive to variations in the cluster formation rate (but not SFR) with time, which may significantly affect their use in tests of cluster mass loss. By comparing simulations with different implementations of cluster formation physics, we find that (even partially) environmentally-independent cluster formation is inconsistent with the brightest cluster-SFR and specific luminosity-$\Sigma_\mathrm{SFR}$ relations, whereas these observables are reproduced by the fiducial, environmentally-varying model. This shows that models in which a constant fraction of stars form in clusters are inconsistent with observations.

Globular cluster ages provide both an important test of models of globular cluster formation and a powerful method to constrain the assembly history of galaxies. Unfortunately, measuring the ages of unresolved old stellar populations has proven challenging. Here, we present a novel technique that combines optical photometry with metallicity constraints from near-infrared spectroscopy in order to measure ages. After testing the method on globular clusters in the Milky Way and its satellite galaxies, we apply our technique to three massive early-type galaxies using data from the SLUGGS Survey. The three SLUGGS galaxies and the Milky Way show dramatically different globular cluster age and metallicity distributions, with NGC 1407 and the Milky Way showing mostly old globular clusters while NGC 3115 and NGC 3377 show a range of globular ages. This diversity implies different galaxy formation histories and that the globular cluster optical colour-metallicity relation is not universal as is commonly assumed in globular cluster studies. We find a correlation between the median age of the metal rich globular cluster populations and the age of the field star populations, in line with models where globular cluster formation is a natural outcome of high intensity star formation.

We present the evolution of the rest-frame ultraviolet (UV) properties of the globular cluster (GC) populations and their host galaxies formed in the E-MOSAICS suite of cosmological hydrodynamical simulations. We compute the luminosities of all clusters associated with 25 simulated Milky Way-mass galaxies, discussed in previous works, in the rest-frame UV and optical bands by combining instantaneous cluster properties (age, mass, metallicity) with simple stellar population models, from redshifts $z=0$ to 10. Due to the rapid fading of young stellar populations in the UV, most of the simulated galaxies do not host GCs bright enough to be individually identified in deep Hubble Space Telescope (HST) observations, even in highly magnified systems. The median age of the most UV-luminous GCs is $<10$ Myr (assuming no extinction), increasing to $\gtrsim 100$ Myr for red optical filters. We estimate that these GCs typically only contribute a few per cent of the total UV luminosity of their host galaxies at any epoch. We predict that the number density of UV-bright proto-GCs (or cluster clumps) will peak between redshifts $z=1-3$. In the main progenitors of Milky Way-mass galaxies, $10$-$20$ per cent of the galaxies at redshifts $1 \lesssim z \lesssim 3$ have clusters brighter than $M_{\rm UV} < -15$, and less than $10$ per cent at other epochs. The brightest cluster in the galaxy sample at $z>2$ is typically $M_{\rm UV} \sim -16$, consistent with the luminosities of observed compact, high-redshift sources.

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 $2\lesssim z \lesssim 3$, which could be tested with observations of lensed galaxies using JWST.

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 galaxy's merger tree, allowing us to identify all progenitors with masses $M_*\gtrsim10^8$ M$_\odot$ for redshifts $1\leq z\leq2.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.

The recent discovery of super-massive black holes (SMBHs) in the centers of high-mass ultra compact dwarf galaxies (UCDs) suggests that at least some UCDs are the stripped nuclear star clusters of lower mass galaxies. Tracing these former nuclei of stripped galaxies provides a unique way to track the assembly history of a galaxy or galaxy cluster. In this paper we present a new method to estimate how many UCDs host an SMBH in their center and thus are stripped galaxy nuclei. We revisit the dynamical mass measurements that suggest many UCDs have more mass than expected from stellar population estimates, which recent observations have shown is due to the presence of an SMBH. We revise the stellar population mass estimates using a new empirical relation between the mass-to-light ratio (M/L) and metallicity, and use this to predict which UCDs are most likely to host an SMBH. This enables us to calculate the fraction of UCDs that host SMBHs across their entire luminosity range for the first time. We then apply the SMBH occupation fraction to the observed luminosity function of UCDs and estimate that in the Fornax and Virgo cluster alone there should be $69^{+32}_{-25}$ stripped nuclei with SMBHs. This analysis shows that stripped nuclei with SMBHs are almost as common in clusters as present-day galaxy nuclei. We estimate the local SMBH number density in stripped nuclei to $3-8\times10^{-3}Mpc^{-3}$, which represents a significant fraction (10-40\%) of the SMBH density in the local Universe. These SMBHs hidden in stripped nuclei will increase expected event rates for tidal disruption events and SMBH-SMBH and SMBH-BH mergers. The existence of numerous stripped nuclei with SMBHs are a direct consequence of hierarchical galaxy formation, but until now their impact on the SMBH density had not been quantified.

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 $\textit{Gaia}$-DR2 kinematics of halo stars. We show that $\sim$2/3 of nearby halo stars have high orbital eccentricities ($e \gtrsim 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]$\sim-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}\lesssim M_*\lesssim 10^{9}\mathrm{M_\odot}$. We show that the median orbital eccentricities of debris are largely unchanged since merger time, implying that this accretion event likely happened at $z\lesssim1.5$. The exact nature of the low $e$ population is unclear, but we hypothesise that it is a combination of $\textit{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.

Stellar haloes encode a fossil record of a galaxy's accretion history, generally in the form of structures of low surface brightness, such as stellar streams. While their low surface brightness makes it challenging to determine their age, metallicity, kinematics and spatial structure, the infalling galaxies also deposit globular clusters (GCs) in the halo, which are bright and therefore easier to observe and characterise. To understand how GCs associated with stellar streams can be used to estimate the stellar mass and the infall time of their parent galaxy, we examine a subset of 15 simulations of galaxies and their star clusters from the E-MOSAICS project. E-MOSAICS is a suite of hydrodynamical simulations incorporating a sub-grid model for GC formation and evolution. We find that more massive accreted galaxies typically contribute younger and more metal rich GCs. This lower age results from a more extended cluster formation history in more massive galaxies. In addition, at fixed stellar mass, galaxies that are accreted later host younger clusters, because they can continue to form GCs without being subjected to environmental influences for longer. This explains the large range of ages observed for clusters associated with the Sagittarius dwarf galaxy in the halo of the Milky Way compared to clusters which are thought to have formed in satellites accreted early in the Milky Way's formation history. Using the ages of the GCs associated with the Sagittarius dwarf, we estimate a virial radius crossing lookback time (infall time) of $9.3 \pm 1.8 Gyr$.

Several models have been advanced to explain the multiple stellar populations observed in globular clusters (GCs). Most models necessitate a large initial population of unenriched stars that provide the pollution for an enriched population, and which are subsequently lost from the cluster. This scenario generally requires clusters to lose $>90$ per cent of their birth mass. We use a suite of 25 cosmological zoom-in simulations of present-day Milky Way-mass galaxies from the \emosaics project to study whether dynamical disruption by evaporation and tidal shocking provides the necessary mass loss. We find that GCs with present-day masses $M>10^5~M_{\odot}$ were only $2$-$4$ times more massive at birth, in conflict with the requirements of the proposed models. This factor correlates weakly with metallicity, gas pressure at birth, or galactocentric radius, but increases towards lower GC masses. To reconcile our results with observational data, either an unphysically steep cluster mass-size relation must be assumed, or the initial enriched fractions must be similar to their present values. We provide the required relation between the initial enriched fraction and cluster mass. Dynamical cluster mass loss cannot reproduce the high observed enriched fractions nor their trend with cluster mass.

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 $\textit{rate}$ among the 72th-94th percentile of the E-MOSAICS galaxies, whereas its $\textit{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\sim4$, sharpening previous limits of $z\sim2$. We identify three massive satellite progenitors and constrain their mass growth and enrichment histories. Two are proposed to correspond to Sagittarius (few $10^8~{\rm M}_\odot$) and the GCs formerly associated with Canis Major ($\sim10^9~{\rm M}_\odot$). The third satellite has no known associated relic and was likely accreted between $z=0.6$-$1.3$. We name this enigmatic galaxy $\textit{Kraken}$ and propose that it is the most massive satellite ($M_*\sim2\times10^9~{\rm M}_\odot$) ever accreted by the MW. We predict that $\sim40\%$ of the Galactic GCs formed ex-situ (in galaxies with masses $M_*=2\times10^7$-$2\times10^9~{\rm M}_\odot$), with $6\pm1$ being former nuclear clusters.

The metal-poor sub-population of globular cluster (GC) systems exhibits a correlation between the GC average colour and luminosity, especially in those systems associated with massive elliptical galaxies. More luminous (more massive) GCs are typically redder and hence more metal-rich. This 'blue tilt' is often interpreted as a mass-metallicity relation stemming from GC self-enrichment, whereby more massive GCs retain a greater fraction of the enriched gas ejected by their evolving stars, fostering the formation of more metal-rich secondary generations. We examine the E-MOSAICS simulations of the formation and evolution of galaxies and their GC populations, and find that their GCs exhibit a colour-luminosity relation similar to that observed in local galaxies, without the need to invoke mass-dependent self-enrichment. We find that the blue tilt is most appropriately interpreted as a dearth of massive, metal-poor GCs: the formation of massive GCs requires high interstellar gas surface densities, conditions that are most commonly fostered by the most massive, and hence most metal rich, galaxies, at the peak epoch of GC formation. The blue tilt is therefore a consequence of the intimate coupling between the small-scale physics of GC formation and the evolving properties of interstellar gas hosted by hierarchically-assembling galaxies.

We discuss some of the key open questions regarding the formation and evolution of globular clusters (GCs) during galaxy formation and assembly within a cosmological framework. The current state-of-the-art for both observations and simulations is described, and we briefly mention directions for future research. The oldest GCs have ages $\ge$ 12.5 Gyr and formed around the time of reionisation. Resolved colour-magnitude diagrams of Milky Way GCs and direct imaging of lensed proto-GCs at z $\sim$ 6 with JWST promise further insight. Globular clusters are known to host multiple populations of stars with variations in their chemical abundances. Recently, such multiple populations have been detected in $\sim$2 Gyr old compact, massive star clusters. This suggests a common, single pathway for the formation of GCs at high and low redshift. The shape of the initial mass function for GCs remains unknown, however for massive galaxies a power-law mass function is favoured. Significant progress has been made recently modelling GC formation in the context of galaxy formation, with success in reproducing many of the observed GC-galaxy scaling relations.

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 $\sim L^\star$ 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\gtrsim1-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. (Abridged)

We compare the predictions of a semi-analytic model for ultra-compact dwarf galaxy (UCD) formation by tidal stripping to the observed properties of globular clusters (GCs) and UCDs in the Fornax and Virgo clusters. For Fornax we find the predicted number of stripped nuclei agrees very well with the excess number of GCs$+$UCDs above the GC luminosity function. GCs$+$UCDs with masses $>10^{7.3}$ M$_\odot$ are consistent with being entirely formed by tidal stripping. Stripped nuclei can also account for Virgo UCDs with masses $>10^{7.3}$ M$_\odot$ where numbers are complete by mass. For both Fornax and Virgo, the predicted velocity dispersions and radial distributions of stripped nuclei are consistent with that of UCDs within $\sim$50-100 kpc but disagree at larger distances where dispersions are too high and radial distributions too extended. Stripped nuclei are predicted to have radially biased anisotropies at all radii, agreeing with Virgo UCDs at clustercentric distances larger than 50 kpc. However, ongoing disruption is not included in our model which would cause orbits to become tangentially biased at small radii. We find the predicted metallicities and central black hole masses of stripped nuclei agree well with the metallicities and implied black hole masses of UCDs for masses $>10^{6.5}$ M$_\odot$. The predicted black hole masses also agree well with that of M60-UCD1, the first UCD with a confirmed central black hole. These results suggest that observed GC$+$UCD populations are a combination of genuine GCs and stripped nuclei, with the contribution of stripped nuclei increasing toward the high-mass end.

We use the Millennium II cosmological simulation combined with the semi-analytic galaxy formation model of Guo et al. (2011) to predict the contribution of galactic nuclei formed by the tidal stripping of nucleated dwarf galaxies to globular cluster (GC) and ultra-compact dwarf galaxy (UCD) populations of galaxies. We follow the merger trees of galaxies in clusters back in time and determine the absolute number and stellar masses of disrupted galaxies. We assume that at all times nuclei have a distribution in nucleus-to-galaxy mass and nucleation fraction of galaxies similar to that observed in the present day universe. Our results show stripped nuclei follow a mass function $N(M) \sim M^{-1.5}$ in the mass range $10^6 < M/M_\odot < 10^8$, significantly flatter than found for globular clusters. The contribution of stripped nuclei will therefore be most important among high-mass GCs and UCDs. For the Milky Way we predict between 1 and 3 star clusters more massive than $10^5 M_\odot$ come from tidally disrupted dwarf galaxies, with the most massive cluster formed having a typical mass of a few times $10^6 M_\odot$, like omega Centauri. For a galaxy cluster with a mass $7 \times 10^{13} M_\odot$, similar to Fornax, we predict $\sim$19 UCDs more massive than $2\times10^6 M_\odot$ and $\sim$9 UCDs more massive than $10^7 M_\odot$ within a projected distance of 300 kpc come from tidally stripped dwarf galaxies. The observed number of UCDs are $\sim$200 and 23, respectively. We conclude that most UCDs in galaxy clusters are probably simply the high mass end of the GC mass function.

Ultra Compact Dwarf Galaxies (UCDs) and dwarf galaxy nuclei have many common properties, such as internal velocity dispersions and colour-magnitude trends, suggesting tidally stripped dwarf galaxies as a possible UCD origin. However, UCDs typically have sizes more than twice as large as nuclei at the same luminosity. We use a GPU-enabled version of the particle-mesh code \textsc{superbox} to study the possibility of turning nucleated dwarf galaxies into UCDs by tidally stripping them in a Virgo-like galaxy cluster. We find that motion in spherical potentials, where close passages happen many times, leads to the formation of compact ($r_h \lesssim 20$ pc) star clusters/UCDs. In contrast, orbital motion where close passages happen only once or twice leads to the formation of extended objects which are large enough to account for the full range of observed UCD sizes. For such motion, we find that dwarf galaxies need close pericentre passages with distances less than 10 kpc to undergo strong enough stripping so that UCD formation is possible. As tidal stripping produces objects with similar properties to UCDs, and our estimates suggest dwarf galaxies have been destroyed in sufficient numbers to explain the observed number of UCDs in M87, we consider tidal stripping to be a likely origin of UCDs. However, comparison with cosmological simulations is needed to determine if the number and spatial distribution of UCDs formed by tidal stripping matches the observations of UCDs in galaxy clusters.