Mass Growth Histories Research Articles

The Mass Assembly History for Galaxies with MaNGA

How galaxies assemble masses through their own star formation or interaction with the external environment is still an important topic in the field of galaxy formation and evolution. We use Value Added Catalogs with galaxy features that are spatially and temporally resolved from Sloan Digital Sky Survey Data Release 17 to investigate the mass growth histories of early-type galaxies (ETGs) and late-type galaxies (LTGs). We find that the mass growth of ETGs is earlier than that of LTGs for massive galaxies (M * > 1010 M ⊙), while low-mass (M * ≤ 1010 M ⊙) ETGs have statistically similar mass assembly histories as low-mass LTGs. The stellar metallicity of all massive galaxies shows a negative gradient and basically does not change with time. However, in low-mass galaxies, the stellar metallicity gradient of elliptical galaxies is negative, and the stellar metallicity gradient of lenticular and spiral galaxies evolves from positive to negative. ETGs are not all in a high-density environment, but exhibit mass dependence. As the tidal strength increases, the star formation rate of low-mass ETGs rapidly decreases. These results support a picture where massive galaxies exhibit inside-out quenching mode, while low-mass galaxies show outside-in quenching mode. Environmental effects play an important role in regulating the mass assembly histories of low-mass ETGs.

From Halos to Galaxies. VII. The Connections between Stellar Mass Growth History, Quenching History, and Halo Assembly History for Central Galaxies

The assembly of galaxies over cosmic time is tightly connected to the assembly of their host dark matter halos. We investigate the stellar mass growth history and the chemical enrichment history of central galaxies in the Sloan Digital Sky Survey Mapping Nearby Galaxies at Apache Point Observatory. We find that the derived stellar metallicity of passive central galaxies is always higher than that of the star-forming ones. This stellar metallicity enhancement becomes progressively larger toward low-mass galaxies (at a given epoch) and earlier epochs (at a given stellar mass), which suggests strangulation as the primary mechanism for star formation quenching in central galaxies not only in the local Universe but also very likely at higher redshifts up to z ∼ 3. We show that at the same present-day stellar mass, passive central galaxies assembled half of their final stellar mass ∼2 Gyr earlier than star-forming central galaxies, which agrees well with the semi-analytic model. Exploring the semi-analytic model, we find that this is because passive central galaxies reside in, on average, more massive halos with a higher halo mass increase rate across cosmic time. As a consequence, passive central galaxies are assembled faster and also quenched earlier than their star-forming counterparts. While at the same present-day halo mass, different halo assembly history also produces a very different final stellar mass of the central galaxy within, and halos assembled earlier host more massive centrals with a higher quenched fraction, in particular around the “golden halo mass” at 1012 M ⊙. Our results call attention back to the dark matter halo as a key driver of galaxy evolution.

EDGE: The direct link between mass growth history and the extended stellar haloes of the faintest dwarf galaxies

ABSTRACT Ultra-faint dwarf galaxies (UFDs) are commonly found in close proximity to the Milky Way and other massive spiral galaxies. As such, their projected stellar ellipticity and extended light distributions are often thought to owe to tidal forces. In this paper, we study the projected stellar ellipticities and faint stellar outskirts of tidally isolated ultra-faints drawn from the ‘Engineering Dwarfs at Galaxy Formation’s Edge’ (EDGE) cosmological simulation suite. Despite their tidal isolation, our simulated dwarfs exhibit a wide range of projected ellipticities (0.03 < ε < 0.85), with many possessing anisotropic extended stellar haloes that mimic tidal tails, but owe instead to late-time accretion of lower mass companions. Furthermore, we find a strong causal relationship between ellipticity and formation time of a UFD, which is robust to a wide variation in the feedback model. We show that the distribution of projected ellipticities in our suite of simulated EDGE dwarfs matches well with a sample of 19 Local Group dwarf galaxies and a sample of 11 isolated dwarf galaxies. Given ellipticity in EDGE arises from an ex-situ accretion origin, the agreement in shape indicates the ellipticities of some observed dwarfs may also originate from a non-tidal scenario. The orbital parameters of these observed dwarfs further support that they are not currently tidally disrupting. If the baryonic content in these galaxies is still tidally intact, then the same may be true for their dark matter content, making these galaxies in our Local Group pristine laboratories for testing dark matter and galaxy formation models.

Galaxy Quenching with Mass Growth History of Galaxy Groups and Clusters: The Importance of Post-processing

We investigate the fraction of quenched satellite galaxies in host galaxy groups and clusters using TNG300 in the IllustrisTNG cosmological magnetohydrodynamical simulations. The simulations show that most satellites are quenched after they fall into their final hosts, and that post-processing is a more dominant mechanism of galaxy quenching than pre-processing. We find that the fraction of quenched satellites at z = 0 increases with host mass, which implies that more massive hosts have higher quenching efficiency because they have more massive groups infalling. Furthermore, we find that hosts that have many early-infall satellites show a higher fraction of quenched satellites at z = 0 than those that have many late-infall satellites, which results in a scatter of the quenched fraction of satellites in a given mass range of hosts at z = 0. Our results highlight the significance of the mass of hosts and the different infall times of satellites in understanding galaxy quenching.

A pair of early- and late-forming galaxy cluster samples: A novel way of studying halo assembly bias assisted by a constrained simulation

The halo assembly bias, a phenomenon referring to dependencies of the large-scale bias of a dark matter halo other than its mass, is a fundamental property of the standard cosmological model. First discovered in 2005 from the Millennium Run simulation, it has been proven very difficult to be detected observationally, with only a few convincing claims of detection so far. The main obstacle lies in finding an accurate proxy of the halo formation time. In this study, by utilizing a constrained simulation that can faithfully reproduce the observed structures larger than 2 Mpc in the local universe, for a sample of 634 massive clusters at z ≤ 0.12, we found their counterpart halos in the simulation and used the mass growth history of the matched halos to estimate the formation time of the observed clusters. This allowed us to construct a pair of early- and late-forming clusters, with a similar mass as measured via weak gravitational lensing, and large-scale biases differing at the ≈3σ level, suggestive of the signature of assembly bias, which is further corroborated by the properties of cluster galaxies, including the brightest cluster galaxy and the spatial distribution and number of member galaxies. Our study paves a way to further detect assembly bias based on cluster samples constructed purely on observed quantities.

A deep learning approach to halo merger tree construction

ABSTRACT A key ingredient for semi-analytic models of galaxy formation is the mass assembly history of haloes, encoded in a tree structure. The most commonly used method to construct halo merger histories is based on the outcomes of high-resolution, computationally intensive N-body simulations. We show that machine learning (ML) techniques, in particular Generative Adversarial Networks (GANs), are a promising new tool to tackle this problem with a modest computational cost and retaining the best features of merger trees from simulations. We train our GAN model with a limited sample of merger trees from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulation suite, constructed using two halo finders–tree builder algorithms: SUBFIND – D-TREES and ROCKSTAR – ConsistentTrees. Our GAN model successfully learns to generate well-constructed merger tree structures with high temporal resolution, and to reproduce the statistical features of the sample of merger trees used for training, when considering up to three variables in the training process. These inputs, whose representations are also learned by our GAN model, are mass of the halo progenitors and the final descendant, progenitor type (main halo or satellite), and distance of a progenitor to that in the main branch. The inclusion of the latter two inputs greatly improves the final learned representation of the halo mass growth history, especially for SUBFIND-like ML trees. When comparing equally sized samples of ML merger trees with those of the EAGLE simulation, we find better agreement for SUBFIND-like ML trees. Finally, our GAN-based framework can be utilized to construct merger histories of low- and intermediate-mass haloes, the most abundant in cosmological simulations.

EDGE: What shapes the relationship between H i and stellar observables in faint dwarf galaxies?

ABSTRACT We show how the interplay between feedback and mass-growth histories introduces scatter in the relationship between stellar and neutral gas properties of field faint dwarf galaxies ($M_{\star }\lessapprox 10^{6} \, \mbox{M}_\mathrm{\odot }$). Across a suite of cosmological, high-resolution zoomed simulations, we find that dwarf galaxies of stellar masses $10^5 \le M_{\star }\le 10^{6} \, \mbox{M}_\mathrm{\odot }$ are bimodal in their cold gas content, being either H i-rich or H i-deficient. This bimodality is generated through the coupling between (i) the modulation of H i contents by the background of ultraviolet radiation (UVB) at late times and (ii) the significant scatter in the stellar-to-halo mass relationship induced by reionization. Furthermore, our H i-rich dwarfs exhibit disturbed and time-variable neutral gas distributions primarily due to stellar feedback. Over the last four billion years, we observe order-of-magnitude changes around the median $M_{\mathrm{H\,\small {I} }}$, factor-of-a-few variations in H i spatial extents, and spatial offsets between H i and stellar components regularly exceeding the galaxies’ optical sizes. Time variability introduces further scatter in the $M_{\star }\!-\! M_{\mathrm{H\,\small {I} }}$ relation and affects a galaxy’s detectability in H i at any given time. These effects will need to be accounted for when interpreting observations of the population of faint, H i-bearing dwarfs by the combination of optical and radio wide, deep surveys.

The Number Densities and Stellar Populations of Massive Galaxies at 3 < z < 6: A Diverse, Rapidly Forming Population in the Early Universe

Abstract We present the census of massive (log(M */M ⊙) > 11) galaxies at 3 < z < 6 identified over the COSMOS/UltraVISTA Ultra-Deep field stripes: consisting of ≈100 and ≈20 high-confidence candidates at 3 < z < 4 and 4 < z < 6, respectively. The 3 < z < 4 population is comprised of post-starburst, UV-star-forming, and dusty star-forming galaxies in roughly equal fractions, while UV-star-forming galaxies dominate at 4 < z < 6 . We account for various sources of biases in the spectral energy distribution (SED) modeling, finding that the treatment of emission line contamination is essential for understanding the number densities and mass growth histories of massive galaxies at z > 3. The significant increase in observed number densities at z ∼ 4 (> × 5 in ≲600 Myr) implies that this is the epoch at which log(M */M ⊙) > 11 galaxies emerge in significant numbers, with stellar ages (≈500–900 Myr) indicating rapid formation epochs as early as z ∼ 7. Leveraging ancillary multiwavelength data sets, we perform panchromatic SED modeling to constrain the total star formation activity of the sample. The star formation activity of the sample is generally consistent with being on the star formation main sequence at the considered redshifts, with ≈15%–25% of the population showing evidence of suppressed star formation rates, indicating that quenching mechanisms are already at play by z ∼ 4. We stack the available Hubble Space Telescope imaging, confirming their compact nature (r e ≲ 2.2 kpc), consistent with expected sizes of high-z star-forming galaxies. Finally, we discuss how our results are in-line with the early formation epochs and short formation timescales inferred from the fossil records of the most massive galaxies in the universe.

Turbulence and its connection to episodic accretion in binary YSOs

ABSTRACT We report signatures of episodic accretion in young stellar objects (YSOs) that emerge in protobinary configurations in a gravoturbulent gas collapse. We find in most of these protobinary systems strong accretion bursts between the two companions with a recurrence time-scale of about 1 kyr. The accretion rate on to the secondary star typically exceeds that on to the primary with a peak value of 2 × 10−2 M⊙ yr−1 for the former and 6 × 10−3 M⊙ yr−1 for the latter. We propose that the secondary companion, which remains more active in its episodes of accretion bursts, especially for the gas cores with subsonic velocity dispersion, may provide observational opportunities to find traces of episodic accretion in the surrounding gas of the embedded YSOs that are in a binary configuration. Also, protostars evolving as single objects in the same environment show fewer accretion bursts and all together a more steady mass growth history. The prestellar cores with subsonic velocity dispersion exhibit an order of magnitude more intense accretion bursts than in the case of cores with supersonic velocity dispersions. The latter shows the formation of some of the protobinaries in which the primary acts as a more actively accreting companion. This can support these binaries to become systems of extreme mass ratio. Moreover, the YSOs in binary configurations with small semimajor axis a ≈ 50 au and high mass ratio q > 0.7 support phases of intense episodic accretion. The eccentricity, however, seems to play no significant role in the occurrence of accretion bursts.

First Structure Formation under the Influence of Gas–Dark Matter Streaming Velocity and Density: Impact of the “Baryons Trace Dark Matter” Approximation

Abstract The impact of streaming between baryons and dark matter on the first structures has been actively explored by recent studies. We investigate how the key results are affected by two popular approximations. One is to implement the streaming by accounting for only the relative motion while assuming “baryons trace dark matter” spatially at the initialization of simulation. This neglects the smoothing on the gas density taking place before the initialization. In our simulation initialized at z i = 200, it overestimates the gas density power spectrum by up to 40% at k ≈ 102 h Mpc−1 at z = 20. Halo mass (M h ) and baryonic fraction in halos ( ) are also overestimated, but the relation between the two remains unchanged. The other approximation tested is to artificially amplify the density/velocity fluctuations in the cosmic mean density to simulate the first minihalos that form in overdense regions. This gives a head start to the halo growth while the subsequent growth is similar to that in the mean density. The growth in a true overdense region, on the other hand, is accelerated gradually in time. For example, raising σ 8 by 50% effectively transforms in the halo mass growth history while, at 2σ overdensity, the growth is accelerated by a constant in redshift: . As a result, halos have grown more massive in the former than in the latter before z ≈ 27 and vice versa after. The –M h relation is unchanged in those cases as well, suggesting that the Population III formation rate for a given M h is insensitive to the tested approximations.

EDGE: from quiescent to gas-rich to star-forming low-mass dwarf galaxies

ABSTRACT We study how star formation is regulated in low-mass field dwarf galaxies ($10^5 \le M_{\star } \le 10^6 \, \mbox{M}_\mathrm{\odot }$), using cosmological high-resolution ($3 \, \mathrm{pc}$) hydrodynamical simulations. Cosmic reionization quenches star formation in all our simulated dwarfs, but three galaxies with final dynamical masses of $3 \times 10^{9} \, \mbox{M}_\mathrm{\odot }$ are subsequently able to replenish their interstellar medium by slowly accreting gas. Two of these galaxies reignite and sustain star formation until the present day at an average rate of $10^{-5} \, \mbox{M}_\mathrm{\odot } \, \text{yr}^{-1}$, highly reminiscent of observed low-mass star-forming dwarf irregulars such as Leo T. The resumption of star formation is delayed by several billion years due to residual feedback from stellar winds and Type Ia supernovae; even at z = 0, the third galaxy remains in a temporary equilibrium with a large gas content but without any ongoing star formation. Using the ‘genetic modification’ approach, we create an alternative mass growth history for this gas-rich quiescent dwarf and show how a small $(0.2\, \mathrm{dex})$ increase in dynamical mass can overcome residual stellar feedback, reigniting star formation. The interaction between feedback and mass build-up produces a diversity in the stellar ages and gas content of low-mass dwarfs, which will be probed by combining next-generation H i and imaging surveys.

Corona-heated Accretion-disk Reprocessing: A Physical Model to Decipher the Melody of AGN UV/Optical Twinkling

Abstract Active galactic nuclei (AGNs) have long been observed to “twinkle” (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here, we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter α, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about 0.01–0.2. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black-hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy.

Predicting fully self-consistent satellite richness, galaxy growth and starformation rates from the STastical sEmi-Empirical modeL steel.

Abstract Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, steel, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST.

Not So Heavy Metals: Black Hole Feedback Enriches the Circumgalactic Medium

Abstract We examine the effects of supermassive black hole (SMBH) feedback on the circumgalactic medium (CGM) using a cosmological hydrodynamic simulation (Romulus25) and a set of four zoom-in “genetically modified” Milky-Way–mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the CGM, we show that O vi is a sensitive indicator of SMBH feedback. First, we calculate the column densities of O vi in simulated Milky-Way–mass galaxies and compare them with observations from the COS-Halos Survey. Our simulations show column densities of O vi in the CGM consistent with those of COS-Halos star-forming and quenched galaxies. These results contrast with those from previous simulation studies which typically underproduce CGM column densities of O vi. We determine that a galaxy’s star formation history and assembly record have little effect on the amount of O vi in its CGM. Instead, column densities of O vi are closely tied to galaxy halo mass and BH growth history. The set of zoom-in, genetically modified Milky-Way–mass galaxies indicates that the SMBH drives highly metal-enriched material out into its host galaxy’s halo, which in turn elevates the column densities of O vi in the CGM.

Investigating the Stellar Mass Growth Histories of Satellite Galaxies as a Function of Infall Time Using Phase-space

Abstract We compile a large sample of nearby galaxies that are satellites of hosts using a well-known Sloan Digital Sky Survey (SDSS) group catalog. From this sample, we create an “ancient infallers” and “recent infallers” subsample, based on the mean infall time predicted from cosmological simulations for galaxies with their location in phase-space. We compare the stellar mass growth histories of the galaxies in these two subsamples, as determined from multiwavelength spectral energy distribution fitting that uses a comprehensive library of star formation history shapes derived from cosmological simulations. By simultaneously controlling for satellite stellar mass and host halo mass, we can clearly see the impact of time spent in their hosts. As we might predict, the ancient infaller population shows clear signs of earlier quenching, especially for lower mass satellites in more massive hosts. More importantly, we find that the effects are not limited to massive hosts. We find that hosts that might be considered low-mass groups (halo masses ∼1013 M ⊙) can significantly alter their satellites, even for massive satellites (stellar masses ∼1010 M ⊙). Intriguingly, we see changes in the mass growth history of the satellites of clusters as early as 8 or 9 Gyr ago, when they had not yet entered the virial radius of their current host. We propose that this could be the result of galaxies being preprocessed in low-mass substructures in the protocluster outskirts, prior to infall.

The Inflow and Outflow Rate Evolution of Local Milky Way–mass Star-forming Galaxies since z = 1.3

Abstract We study the gas inflow rate (ζ inflow) and outflow rate (ζ outflow) evolution of local Milky Way–mass star-forming galaxies (SFGs) since z = 1.3. The stellar mass growth history of Milky Way–mass progenitor SFGs is inferred from the evolution of the star formation rate (SFR)−stellar mass (M *) relation, and the gas mass (M gas) is derived using the recently established gas-scaling relations. With the growth curve, the net inflow rate κ is quantified at each cosmic epoch. At z ∼ 1.3, κ is comparable with the SFR, whereas it rapidly decreases to ∼0.15 × SFR at z = 0. We then constrain the average outflow rate ζ outflow of progenitor galaxies by modeling the evolution of their gas-phase metallicity. The best-fit ζ outflow is found to be (0.5–0.8) × SFR. Combining κ and ζ outflow, we finally investigate the evolution of ζ inflow since z = 1.3. We find that ζ inflow rapidly decreases by ∼80% from z = 1.3 to z = 0.5. At z < 0.5, ζ inflow continuously decreases but with a much lower decreasing rate. Implications of these findings on galaxy evolution are discussed.

The E-MOSAICS project: tracing galaxy formation and assembly with the age–metallicity distribution of globular clusters

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.

Optical integral field spectroscopy observations applied to simulated galaxies: testing the fossil record method

By means of post-processed cosmological Hydrodynamics simulations we explore the ability of the fossil record method to recover the stellar mass, age gradients, and global/radial star formation and mass growth histories of galaxies observed with an optical Integral Field Unity. We use two simulations of representative Milky Way-sized galaxies, one disk dominated and another bulge-dominated. We generate data cubes emulating MaNGA/SDSS-IV instrumental and observational conditions, and analyze them using the Pipe3D pipeline. For optimal MaNGA instrumental/observational setups, the global masses, ages, and dust extinctions of both galaxies are well recovered within the uncertainties, while for the radial distributions, the trend is to get slightly flatter age gradients and weaker inside-out growth modes than the true ones. The main bias is in recovering younger stellar populations in the inner, older regions, and slightly older stellar populations in the outer, younger regions, as well as lower global stellar masses. These trends are enhanced as lower is the spatial sampling/resolution (less number of fibers), more inclined and dust-attenuated is the observed galaxy, lower is the signal-to-noise ratio, and worse is the seeing. The recovered properties and histories are more affected for the disk-dominated galaxy with a prominent inside-out growth mode than for the bulge-dominated galaxy assembled earlier and more coherently. We make publicly available the code and the mock IFU datacubes used in this paper.

The formation and assembly history of the Milky Way revealed by its globular cluster population

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.

Field spheroid-dominated galaxies in a Λ-CDM Universe

Context. Understanding the formation and evolution of early-type, spheroid-dominated galaxies is an open question within the context of the hierarchical clustering scenario, particularly in low-density environments. Aims. Our goal is to study the main structural, dynamical, and stellar population properties and assembly histories of field spheroid-dominated galaxies formed in a Λ-cold dark matter (Λ-CDM) scenario to assess to what extent they are consistent with observations. Methods. We selected spheroid-dominated systems from a Λ-CDM simulation that includes star formation (SF), chemical evolution, and supernova feedback. The sample is made up of 18 field systems with MStar ≲ 6 × 1010M⊙ that are dominated by the spheroid component. For this sample we estimated the fundamental relations of ellipticals and compared them with current observations. Results. The simulated spheroid galaxies have sizes that are in good agreement with observations. The bulges follow a Sersic law with Sersic indexes that correlate with the bulge-to-total mass ratios. The structural-dynamical properties of the simulated galaxies are consistent with observed Faber–Jackson, fundamental plane, and Tully–Fisher relations. However, the simulated galaxies are bluer and with higher star formation rates (SFRs) than the observed isolated early-type galaxies. The archaeological mass growth histories show a slightly delayed formation and more prominent inside-out growth mode than observational inferences based on the fossil record method. Conclusions. The main structural and dynamical properties of the simulated spheroid-dominated galaxies are consistent with observations. This is remarkable since our simulation has not been calibrated to match them. However, the simulated galaxies are blue and star-forming, and with later stellar mass growth histories compared to observational inferences. This is mainly due to the persistence of extended discs in the simulations. The need for more efficient quenching mechanisms able to avoid further disc growth and SF is required in order to reproduce current observational trends.