ABSTRACT
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} per cent 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 (a few 108 M⊙) and the GCs formerly associated with Canis Major (${\sim }10^9\, \mbox{M$_\odot $}$). The third satellite has no known associated relic and was likely accreted between z = 0.6 and 1.3. We name this enigmatic galaxy Kraken and propose that it is the most massive satellite ($M_*\sim 2\times 10^9\, \mbox{M$_\odot $}$) ever accreted by the MW. We predict that ∼40 per cent of the Galactic GCs formed ex situ (in galaxies with masses M* = 2 × 107–$2\times 10^9\, \mbox{M$_\odot $}$), with 6 ± 1 being former nuclear clusters.
Abstract
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 realizing 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 characterized 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_*\gtrsim 10^8\,{\rm M}_\odot$ for redshifts 1 ≤ z ≤ 2.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.
An EAGLE’s view of ex situ galaxy growth Davison, Thomas A; Norris, Mark A; Pfeffer, Joel L ...
Monthly Notices of the Royal Astronomical Society,
09/2020, Volume:
497, Issue:
1
Journal Article
Peer reviewed
Open access
ABSTRACT
Modern observational and analytical 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 field of view (FoV) Integral Field Unit (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 100 and 25 cMpc3 EAGLE (Evolution and Assembly of GaLaxies and their Environments) 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 × 107 to 1.9 × 1012 M⊙. For those galaxies of mass M* > 5 × 108 M⊙, we find that the total ex situ mass fraction is greater for more extended galaxies at fixed mass. When categorizing 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 counterintuitive result may be due to high passing velocities within large cluster haloes inhibiting efficient accretion on to individual galaxies.
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.
ABSTRACT
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 ($\mbox{$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 $\mbox{$r_{M_\star }$}=1$:$31^{+34}_{-16}$ and $\mbox{$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 organizing these accretion events into the most detailed reconstruction to date of the Milky Way’s merger tree.
ABSTRACT
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 MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE 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 ∼ 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.
The recent discovery of supermassive black holes (SMBHs) in high mass ultra-compact dwarf galaxies (UCDs) suggests that at least some UCDs are the nuclear star clusters of stripped galaxies. In this ...paper we present a new method to estimate how many UCDs host an SMBH 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 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 to predict which UCDs most likely host an SMBH. We 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 clusters alone there should be stripped nuclei with SMBHs. This analysis shows that stripped nuclei are almost as common in clusters as present-day galaxy nuclei. We estimate the SMBH number density caused by stripped nuclei to be (2-8) × 10−3 Mpc−3, which represents a significant fraction (8%-32%) of the SMBH density in the local universe. These SMBHs hidden in stripped nuclei 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.
ABSTRACT
The α-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 α/Fe–Fe/H distribution of field stars and GCs in 25 Milky Way–mass galaxies. The α/Fe–Fe/H distribution of GCs largely follows that of the field stars and can also therefore be used as tracers of the α/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 α/Fe at fixed Fe/H. Therefore, if a GC is observed to have low α/Fe for its Fe/H there is an increased possibility that this GC was accreted recently alongside a dwarf galaxy. There is a wide range of shapes for the field star α/Fe–Fe/H distribution, with a notable subset of galaxies exhibiting bimodal distributions, in which the high α/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 and 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.
ABSTRACT
The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the emp-Pathfindersuite of cosmological zoom-in ...Milky Way-mass simulations. These simulations contain a subgrid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside their host galaxies across cosmic time. As a key ingredient in these simulations, we include the physics of the multiphase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects star cluster formation and evolution. We consider two different star formation prescriptions: a constant star formation efficiency per free-fall time, as well as an environmentally dependent, turbulence-based prescription. We identify two key results drawn from these simulations. First, we find that the tidal shock-driven disruption caused by the graininess of the cold ISM produces old ($\tau \gt 10~\mbox{${\rm Gyr}$}$) stellar cluster populations with properties that are in excellent agreement with the observed populations in the Milky Way and M31. Importantly, the addition of the cold ISM addresses the areas of disagreement found in previous simulations that lacked the cold gas phase. Secondly, we find that the formation of stellar clusters is extremely sensitive to the baryonic physics that govern the properties of the cold, dense gas reservoir in the galaxy. This implies that the demographics of the stellar cluster population represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM.