skip to main content

Title: The contribution of in situ and ex situ star formation in early-type galaxies: MaNGA versus IllustrisTNG

We compare stellar mass surface density, metallicity, age, and line-of-sight velocity dispersion profiles in massive ($M_*\ge 10^{10.5}\, \mathrm{M_\odot }$) present-day early-type galaxies (ETGs) from the MaNGA survey with simulated galaxies from the TNG100 simulation of the IllustrisTNG suite. We find an excellent agreement between the stellar mass surface density profiles of MaNGA and TNG100 ETGs, both in shape and normalization. Moreover, TNG100 reproduces the shapes of the profiles of stellar metallicity and age, as well as the normalization of velocity dispersion distributions of MaNGA ETGs. We generally also find good agreement when comparing the stellar profiles of central and satellite galaxies between MaNGA and TNG100. An exception is the velocity dispersion profiles of very massive ($M_*\gtrsim 10^{11.5}\, \mathrm{M_\odot }$) central galaxies, which, on average, are significantly higher in TNG100 than in MaNGA ($\approx 50\, \mathrm{km\, s^{-1}}$). We study the radial profiles of in situ and ex situ stars in TNG100 and discuss the extent to which each population contributes to the observed MaNGA profiles. Our analysis lends significant support to the idea that high-mass ($M_*\gtrsim 10^{11}\, \mathrm{M_\odot }$) ETGs in the present-day Universe are the result of a merger-driven evolution marked by major mergers that tend to homogenize the more » stellar populations of the progenitors in the merger remnant.

« less
; ; ; ; ; ; ; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 5651-5670
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this

    Early-type galaxies (ETGs) possess total density profiles close to isothermal, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at 1.5Reff (effective radius) that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on these kinematic features through mock ETGs in IllustrisTNG. We create mock Integral-Field-Unit observations to extract projected stellar kinematic features for 207 ETGs with stellar mass $M_{\ast }\geqslant 10^{11} \, \mathrm{M_{\odot}}$ in TNG100-1. The mock observations reproduce the key outer (1.5Reff) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope γouter and the kurtosis h4’s gradient. We find that h4 is uncorrelated with stellar orbital anisotropy beyond Reff; instead, we find that the variations in γouter and outer h4 (a good proxy for h4 gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the origin of non-Gaussian outer kinematic structure inmore »ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates the two quantities.

    « less
  2. ABSTRACT We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at z = 0 with stellar masses $10^{10.7}\, \mathrm{M}_{\odot } \leqslant M_{\ast } \leqslant 10^{11.9}\, \mathrm{M}_{\odot }$, the total power-law slope has a mean of 〈γ′〉 = 2.011 ± 0.007 and a scatter of $\sigma _{\gamma ^{\prime }} = 0.171$ over the radial range 0.4–4 times the stellar half-mass radius. Several correlations between γ′ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction, and in situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, γ′ is almost constant with redshift below z = 2. Through analysing IllustrisTNG model variations, we show that black hole kinetic winds are crucial to lowering γ′ and matching observed galaxy correlations. The effects of stellar winds on γ′ are subdominant compared to active galactic nucleus (AGN) feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter haloes are well described by steeper than NFW profiles, and they are steeper in the full physics (FP) run than their counterpartsmore »in the dark matter-only (DMO) run. Their inner density slopes anticorrelate (remain constant) with the halo mass in the FP (DMO) run, and anticorrelate with the halo concentration parameter c200 in both the types of runs. The dark matter haloes of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.« less
  3. ABSTRACT We perform a consistent comparison of the mass and mass profiles of massive (M⋆ > 1011.4 M⊙) central galaxies at z ∼ 0.4 from deep Hyper Suprime-Cam (HSC) observations and from the Illustris, TNG100, and Ponos simulations. Weak lensing measurements from HSC enable measurements at fixed halo mass and provide constraints on the strength and impact of feedback at different halo mass scales. We compare the stellar mass function (SMF) and the Stellar-to-Halo Mass Relation (SHMR) at various radii and show that the radius at which the comparison is performed is important. In general, Illustris and TNG100 display steeper values of α where $M_{\star } \propto M_{\rm vir}^{\alpha }$. These differences are more pronounced for Illustris than for TNG100 and in the inner rather than outer regions of galaxies. Differences in the inner regions may suggest that TNG100 is too efficient at quenching in situ star formation at Mvir ≃ 1013 M⊙ but not efficient enough at Mvir ≃ 1014 M⊙. The outer stellar masses are in excellent agreement with our observations at Mvir ≃ 1013 M⊙, but both Illustris and TNG100 display excess outer mass as Mvir ≃ 1014 M⊙ (by ∼0.25 and ∼0.12 dex, respectively). We argue that reducing stellar growth at earlymore »times in $M_\star \sim 10^{9-10} \, \mathrm{M}_{\odot }$ galaxies would help to prevent excess ex-situ growth at this mass scale. The Ponos simulations do not implement AGN feedback and display an excess mass of ∼0.5 dex at r < 30 kpc compared to HSC which is indicative of overcooling and excess star formation in the central regions. The comparison of the inner profiles of Ponos and HSC suggests that the physical scale over which the central AGN limits star formation is r ≲ 20 kpc. Joint comparisons between weak lensing and galaxy stellar profiles are a direct test of whether simulations build and deposit galaxy mass in the correct dark matter haloes and thereby provide powerful constraints on the physics of feedback and galaxy growth. Our galaxy and weak lensing profiles are publicly available to facilitate comparisons with other simulations.« less
  4. Abstract This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, α-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at z ≤ 0.08 from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to stellar population synthesis models. This requires spectra with high signal-to-noise which we achieve by stacking in bins of luminosity (Lr) and central velocity dispersion (σ0). We find that: 1) There is a bimodality in the S0 population: S0s more massive than 3 × 1010M⊙ show stronger velocity dispersion and age gradients (age and σr decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above 2 × 1011M⊙ the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. 2) S0s have steeper velocity dispersion profiles than fast rotating elliptical galaxies (E-FRs) of the same luminositymore »and velocity dispersion. The kinematic profiles and stellar population gradients of E-FRs are both more similar to those of slow rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. 3) At fixed σ0, more luminous S0s and E-FRs are younger, more metal rich and less α-enhanced. Evidently for these galaxies, the usual statement that ‘massive galaxies are older’ is not true if σ0 is held fixed.« less
  5. ABSTRACT We present estimates of stellar population (SP) gradients from stacked spectra of slow rotator (SR) and fast rotator (SR) elliptical galaxies from the MaNGA-DR15 survey. We find that (1) FRs are ∼5 Gyr younger, more metal rich, less α-enhanced and smaller than SRs of the same luminosity Lr and central velocity dispersion σ0. This explains why when one combines SRs and FRs, objects which are small for their Lr and σ0 tend to be younger. Their SP gradients are also different. (2) Ignoring the FR/SR dichotomy leads one to conclude that compact galaxies are older than their larger counterparts of the same mass, even though almost the opposite is true for FRs and SRs individually. (3) SRs with σ0 ≤ 250 km s−1 are remarkably homogeneous within ∼Re: they are old, α-enhanced, and only slightly supersolar in metallicity. These SRs show no gradients in age and M*/Lr, negative gradients in metallicity, and slightly positive gradients in [α/Fe] (the latter are model dependent). SRs with σ0 ≥ 250 km s−1 are slightly younger and more metal rich, contradicting previous work suggesting that age increases with σ0. They also show larger M*/Lr gradients. (4) Self-consistently accounting for M*/L gradients yields Mdyn ≈ M* because gradients reducemore »Mdyn by ∼0.2 dex while only slightly increasing the M* inferred using a Kroupa (not Salpeter) initial mass function. (5) The SR population starts to dominate the counts above $M_*\ge 3\times 10^{11}\, \mathrm{M}_\odot$; this is the same scale at which the size–mass correlation and other scaling relations change. Our results support the finding that this is an important mass scale that correlates with the environment and above which mergers matter.« less