skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Early-type galaxy density profiles from IllustrisTNG – II. Evolutionary trend of the total density profile
ABSTRACT We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from z = 0 to 4 and measure the power-law slope γ′ of the total density profile for their main progenitors. We find that their slopes γ′ steepen on average during z ∼ 4–2, then becoming shallower until z = 1, after which they remain almost constant, aside from a residual trend of becoming shallower towards z = 0. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since z = 3, and after z ≈ 1, they remain nearly invariant with a mild increase towards z = 0. These evolutionary trends are mass dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analysing the impact of mergers and AGN feedback on the progenitors’ evolution, we conjecture a multiphase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen γ′ from z = 4 to 2, whereas subsequent AGN feedback (especially in the kinetic mode) makes γ′ shallower again from z = 2 to 1. Afterwards, passive evolution from z = 1 to 0, mainly through gas-poor mergers, mildly decreases γ′ and maintains the overall mass distribution close to isothermal.  more » « less
Award ID(s):
1814053 1814259
PAR ID:
10104818
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
490
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5722 to 5738
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 counterparts 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. 
    more » « less
  2. ABSTRACT Rotation curves of galaxies probe their total mass distributions, including dark matter. Dwarf galaxies are excellent systems to investigate the dark matter density distribution, as they tend to have larger fractions of dark matter compared to higher mass systems. The core-cusp problem describes the discrepancy found in the slope of the dark matter density profile in the centres of galaxies (β*) between observations of dwarf galaxies (shallower cores) and dark matter-only simulations (steeper cusps). We investigate β* in six nearby spiral dwarf galaxies for which high-resolution CO J = 1–0 data were obtained with ALMA (Atacama Large Millimeter/submillimeter Array). We derive rotation curves and decompose the mass profile of the dark matter using our CO rotation curves as a tracer of the total potential and 4.5 $$\mu$$m photometry to define the stellar mass distribution. We find 〈β*〉 = 0.6 with a standard deviation of ±0.1 among the galaxies in this sample, in agreement with previous measurements in this mass range. The galaxies studied are on the high stellar mass end of dwarf galaxies and have cuspier profiles than lower mass dwarfs, in agreement with other observations. When the same definition of the slope is used, we observe steeper slopes than predicted by the FIRE and NIHAO simulations. This may signal that these relatively massive dwarfs underwent stronger gas inflows towards their centres than predicted by these simulations, that these simulations overpredict the frequency of accretion or feedback events, or that a combination of these or other effects are at work. 
    more » « less
  3. null (Ed.)
    ABSTRACT The James Webb Space Telescope (JWST) is expected to observe galaxies at z > 10 that are presently inaccessible. Here, we use a self-consistent empirical model, the universemachine, to generate mock galaxy catalogues and light-cones over the redshift range z = 0−15. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy–halo relationships, and galaxy–galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at z < 10. We predict that Cycle 1 JWST surveys will very likely detect galaxies with M* > 107 M⊙ and/or M1500 < −17 out to at least z ∼ 13.5. Number density uncertainties at z > 12 expand dramatically, so efforts to detect z > 12 galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits (M* < 107 M⊙ or M1500 > −17). For reionization models, extrapolating luminosity functions with a constant faint-end slope from M1500 = −17 down to M1500 = −12 gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to z ∼ 12. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and light-cones for common fields are available online. 
    more » « less
  4. ABSTRACT 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 in 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. 
    more » « less
  5. ABSTRACT We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, namely SIMBA, EAGLE, and IllustrisTNG, versus observations from z ∼ 0 to 2. These simulations all rely on similar subresolution prescriptions to model cold interstellar gas that they cannot represent directly, and qualitatively reproduce the observed z ≈ 0 H i and H2 mass functions (HIMFs and H2MFs, respectively), CO(1–0) luminosity functions (COLFs), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density μ*, with some quantitative differences. To compare to the COLF, we apply an H2-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in αCO and significant differences between models. Using this, predicted z = 0 COLFs agree better with data than predicted H2MFs. Out to z ∼ 2, EAGLE’s and SIMBA’s HIMFs and COLFs strongly increase, while IllustrisTNG’s HIMF declines and COLF evolves slowly. EAGLE and simba reproduce high-LCO(1–0) galaxies at z ∼ 1–2 as observed, owing partly to a median αCO(z = 2) ∼ 1 versus αCO(z = 0) ∼ 3. Examining H i, H2, and CO scaling relations, their trends with M* are broadly reproduced in all models, but EAGLE yields too little H i in green valley galaxies, IllustrisTNG and SIMBA overproduce cold gas in massive galaxies, and SIMBA overproduces molecular gas in small systems. Using SIMBA variants that exclude individual active galactic nucleus (AGN) feedback modules, we find that SIMBA’s AGN jet feedback is primarily responsible by lowering cold gas contents from z ∼ 1 → 0 by suppressing cold gas in $$M_*\gtrsim 10^{10}{\rm \,M}_\odot$$ galaxies, while X-ray feedback suppresses the formation of high-μ* systems. 
    more » « less