We present a new method of quantifying a galaxy's accretion history from its integrated spectrum alone. Using full spectral fitting and calibrated regularization techniques we show how we can accurately derive a galaxy's mass distribution in age-metallicity space and further separate this into stellar populations from different chemical enrichment histories. By exploiting the fact that accreted lower mass galaxies will exhibit an offset to lower metallicities at fixed age compared to the in-situ stellar population, we quantify the fraction of light that comes from past merger events, that are long since mixed in phase-space and otherwise indistinguishable. Empirical age-metallicity relations (AMRs) parameterized for different galaxy masses are used to identify the accreted stellar populations and link them back to the progenitor galaxy's stellar mass. This allows us to not only measure the host galaxy's total ex-situ mass fraction (f acc ), but also quantify the relative amount of accreted material deposited by satellite galaxies of different masses, i.e. the accreted satellite mass function in analogy to the subhalo mass function. Using mock spectra of simulated, present-day galaxies from the EAGLE suite we demonstrate that our method can recover the total accreted fraction to within ≈ 12 %, the stellar mass of the most massive accreted subhalo to within ≈ 26 % and the slope of the accreted satellite mass function to within ≈ 16 % of the true values from the EAGLE merger trees. Future application of this method to observations could potentially provide us accretion histories of hundreds of individual galaxies, for which deep integrated light spectroscopy is available.

Recent photometric observations revealed a massive, extended (MGC ≳ 105 M⊙; Rh ∼ 14 pc) globular cluster (GC) in the central region (D3D ≲ 100 pc) of the low-mass (M* ∼ 5 × 106 M⊙) dwarf irregular galaxy Pegasus. This massive GC offers a unique opportunity to study star cluster inspiral as a mechanism for building up nuclear star clusters, and the dark matter (DM) density profile of the host galaxy. Here, we present spectroscopic observations indicating that the GC has a systemic velocity of ΔV = 3 ± 8 km s−1 relative to the host galaxy, and an old, metal-poor stellar population. We run a suite of orbital evolution models for a variety of host potentials (cored to cusped) and find that the GC’s observed tidal radius (which is ∼3 times larger than the local Jacobi radius), relaxation time, and relative velocity are consistent with it surviving inspiral from a distance of Dgal ≳ 700 pc (up to the maximum tested value of Dgal = 2000 pc). In successful trials, the GC arrives to the galaxy centre only within the last ∼1.4 ± 1 Gyr. Orbits that arrive in the centre and survive are possible in DM haloes of nearly all shapes, however to satisfy the GC’s structural constraints a galaxy DM halo with mass MDM ≃ 6 ± 2 × 109 M⊙, concentration c ≃ 13.7 ± 0.6, and an inner slope to the DM density profile of −0.9 ≤ γ ≤ −0.5 is preferred. The gas densities necessary for its creation and survival suggest the GC could have formed initially near the dwarf’s centre, but then was quickly relocated to the outskirts where the weaker tidal field permitted an increased size and relaxation time – with the latter preserving the former during subsequent orbital decay.

Current instruments and spectral analysis programs are now able to decompose the integrated spectrum of a stellar system into distributions of ages and metallicities. The reliability of these methods has rarely been tested on nearby systems with resolved stellar ages and metallicities. Here we derive the age–metallicity distribution of M54, the nucleus of the Sagittarius dwarf spheroidal galaxy, from its integrated Multi-Unit Spectroscopic Explorer (MUSE) spectrum. We find a dominant old (8–14 Gyr), metal-poor (−1.5 dex) component and a young (1 Gyr), metal-rich (+0.25 dex) component—consistent with the complex stellar populations measured from individual stars in the same MUSE data set. There is excellent agreement between the (mass-weighted) average age and metallicity of the resolved and integrated analyses. Differences are only 3% in age and 0.2 dex metallicity. By co-adding individual stars to create M54's integrated spectrum, we show that the recovered age–metallicity distribution is insensitive to the magnitude limit of the stars or the contribution of blue horizontal branch stars—even when including additional blue wavelength coverage from the WiFeS Atlas of Galactic Globular cluster Spectra survey. However, we find that the brightest stars can induce the spurious recovery of an old (>8 Gyr), metal-rich (+0.25 dex) stellar population, which is otherwise not expected from our understanding of chemical enrichment in M54. The overall derived stellar mass-to-light ratio of M54 is M/L V = 1.46 with a scatter of 0.22 across the field of view, which we attribute to the stochastic contribution of a young, metal-rich component. These findings provide strong evidence that complex stellar population distributions can be reliably recovered from integrated spectra of extragalactic systems.

Galaxy mergers are instrumental in dictating the final mass, structure, stellar populations, and kinematics of galaxies. Cosmological galaxy simulations indicate that the most massive galaxies at z = 0 are dominated by high fractions of ‘ex-situ’ stars, which formed first in distinct independent galaxies, and then subsequently merged into the host galaxy. Using spatially resolved MUSE spectroscopy we quantify and map the ex-situ stars in thirteen massive early-type galaxies. We use full spectral fitting together with semi-analytic galaxy evolution models to isolate the signatures in the galaxies’ light which are indicative of ex-situ populations. Using the large MUSE field of view we find that all galaxies display an increase in ex-situ fraction with radius, with massive and more extended galaxies showing a more rapid increase in radial ex-situ fraction (reaching values between ∼30 per cent and 100 per cent at 2 effective radii) compared to less massive and more compact sources (reaching between ∼5 per cent and 40 per cent ex-situ fraction within the same radius). These results are in line with predictions from theory and simulations which suggest ex-situ fractions should increase significantly with radius at fixed mass for the most massive galaxies.

Context. Modern astronomical observations give unprecedented access to the physical properties of nearby galaxies, including spatially resolved stellar populations. However, observations can only give a present-day view of the Universe, whereas cosmological simulations give access to the past record of the processes that galaxies have experienced in their evolution. To connect the events that happened in the past with galactic properties as seen today, simulations must be taken to a common ground before being compared to observations. Therefore, a dedicated effort is needed to forward-model simulations into the observational plane. Aims. We emulate data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which is the largest integral field spectroscopic galaxy survey to date with its 10 000 nearby galaxies of all types. For this, we use the latest hydro-cosmological simulations IllustrisTNG to generate MaNGIA (Mapping Nearby Galaxies with IllustrisTNG Astrophysics), a mock MaNGA sample of similar size that emulates observations of galaxies for stellar population analysis. Methods. We chose TNG galaxies to match the MaNGA sample selection in terms of mass, size, and redshift in order to limit the impact of selection effects. We produced MaNGA-like datacubes from all simulated galaxies, and processed them with the stellar population analysis code pyPipe3D. This allowed us to extract spatially resolved maps of star formation history, age, metallicity, mass, and kinematics, following the same procedures used as part of the official MaNGA data release. Results. This first paper presents the approach used to generate the mock sample and provides an initial exploration of its properties. We show that the stellar populations and kinematics of the simulated MaNGIA galaxies are overall in good agreement with observations. Specific discrepancies, especially in the age and metallicity gradients in low- to intermediate-mass regimes and in the kinematics of massive galaxies, require further investigation, but are likely to uncover new physical understanding. We compare our results to other attempts to mock similar observations, all of smaller datasets. Conclusions. Our final dataset is released with this publication, consisting of ≳10 000 post-processed datacubes analysed with pyPipe3D, along with the codes developed to create it. Future work will employ modern machine learning and other analysis techniques to connect observations of nearby galaxies to their cosmological evolutionary past.