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
The impact of quenching on galaxy profiles in the simba simulation
ABSTRACT We study specific star formation rate (sSFR) and gas profiles of star-forming (SF) and green valley (GV) galaxies in the simba cosmological hydrodynamic simulation. SF galaxy half-light radii (Rhalf) at z = 0 and their evolution (∝(1 + z)−0.78) agree with observations. Passive galaxy Rhalf agree with observations at high redshift, but by z = 0 are too large, owing to numerical heating. We compare simbaz = 0 sSFR radial profiles for SF and GV galaxies to observations. simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in GV galaxies and massive SF systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, simba yields GV galaxies with strongly dropping sSFR profiles beyond ≳Rhalf, regardless of active galactic nucleus feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. simba’s satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At z = 2, simba does not show central depressions in massive SF galaxies, suggesting simba’s X-ray feedback should be more active at high-z. High-resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in z = 0 GV galaxies represents a unique success of simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms.
more »
« less
- NSF-PAR ID:
- 10185176
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 494
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 6053 to 6071
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
ABSTRACT We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\lesssim$4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $\gt t_{\rm H}^{-1}$ to $\lt 0.2t_{\rm H}^{-1}$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\lesssim 1{{\ \rm per\ cent}}$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the H2 content, but for $M_*\gtrsim 10^{10.5} \,\mathrm{ M}_{\odot }$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift, ∝(1 + z)3.5, but the quenching rate evolves much more slowly, ∝(1 + z)0.9; there are insufficient mergers to explain the quenching rate at $z\lesssim 1.5$. simba first quenches galaxies at $z\gtrsim 3$, with a number density in good agreement with observations. The quenching time-scales τq are strongly bimodal, with ‘slow’ quenchings (τq ∼ 0.1tH) dominating overall, but ‘fast’ quenchings (τq ∼ 0.01tH) dominating in M* ∼ 1010–1010.5 M$\odot$ galaxies, likely induced by simba’s jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.
-
more » « less
Abstract Current models of galaxy formation require strong feedback from active galactic nuclei (AGN) to explain the observed lack of star formation in massive galaxies since
z ≈ 2, but direct evidence of this energy input is limited. We use the SIMBA cosmological galaxy formation simulations to assess the ability of thermal Sunyaev–Zel’dovich (tSZ) measurements to provide such evidence, by mapping the pressure structure of the circumgalactic medium around massivez ≈ 0.2–1.5 galaxies. We undertake a stacking approach to calculate the total tSZ signal and its radial profile in simulations with varying assumptions of AGN feedback, and we assess its observability with current and future telescopes. By convolving our predictions with the 2.′1 beam of the Atacama Cosmology Telescope, we show that current observations atz ≈ 1 are consistent with SIMBA’s fiducial treatment of AGN feedback and inconsistent with SIMBA models without feedback. Atz ≈ 0.5, observational signals lie between SIMBA run with and without AGN feedback, suggesting AGN in SIMBA may inject too much energy at late times. By convolving our data with a 9.″5 beam corresponding to the TolTEC camera on the Large Millimeter Telescope Alfonso Serrano, we predict a unique profile for AGN feedback that can be distinguished with future higher-resolution measurements. Finally, we explore a novel approach to quantify the nonspherically symmetric features surrounding our galaxies by plotting radial profiles representing the component of the stack with m-fold symmetry. -
more » « less
ABSTRACT We investigate the origin of rare star formation in an otherwise red-and-dead population of S0 galaxies, using spatially resolved spectroscopy. Our sample consists of 120 low redshift (z < 0.1) star-forming S0 (SF-S0) galaxies from the SDSS-IV MaNGA DR15. We have selected this sample after a visual inspection of deep images from the DESI Legacy Imaging Surveys DR9 and the Subaru/HSC-SSP survey PDR3 to remove contamination from spiral galaxies. We also construct two control samples of star-forming spirals (SF-Sps) and quenched S0s (Q-S0s) to explore their evolutionary link with the star-forming S0s. To study star formation at resolved scales, we use dust-corrected H α luminosity and stellar density (Σ⋆) maps to construct radial profiles of star formation rate (SFR) surface density (ΣSFR) and specific SFR (sSFR). Examining these radial profiles, we find that star formation in SF-S0s is centrally dominated as opposed to disc-dominated star formation in spirals. We also compared various global (size–mass relation, bulge-to-total luminosity ratio) and local (central stellar velocity dispersion) properties of SF-S0s to those of the control sample galaxies. We find that SF-S0s are structurally similar to the quenched S0s and are different from star-forming spirals. We infer that SF-S0s are unlikely to be fading spirals. Inspecting stellar and gas velocity maps, we find that more than $50{{\ \rm per\ cent}}$ of the SF-S0 sample shows signs of recent galaxy interactions such as kinematic misalignment, counter-rotation, and unsettled kinematics. Based on these results, we conclude that in our sample of SF-S0s, star formation has been rejuvenated, with minor mergers likely to be a major driver.
-
ABSTRACT We compare the star-forming main sequence (SFMS) of galaxies – both integrated and resolved on 1 kpc scales – between the high-resolution TNG50 simulation of IllustrisTNG and observations from the 3D-HST slitless spectroscopic survey at z ∼ 1. Contrasting integrated star formation rates (SFRs), we find that the slope and normalization of the star-forming main sequence in TNG50 are quantitatively consistent with values derived by fitting observations from 3D-HST with the Prospector Bayesian inference framework. The previous offsets of 0.2–1 dex between observed and simulated main-sequence normalizations are resolved when using the updated masses and SFRs from Prospector. The scatter is generically smaller in TNG50 than in 3D-HST for more massive galaxies with M*> 1010 M⊙, by ∼10–40 per cent, after accounting for observational uncertainties. When comparing resolved star formation, we also find good agreement between TNG50 and 3D-HST: average specific star formation rate (sSFR) radial profiles of galaxies at all masses and radii below, on, and above the SFMS are similar in both normalization and shape. Most noteworthy, massive galaxies with M*> 1010.5 M⊙, which have fallen below the SFMS due to ongoing quenching, exhibit a clear central SFR suppression, in both TNG50 and 3D-HST. In contrast, the original Illustris simulation and a variant TNG run without black hole kinetic wind feedback, do not reproduce the central SFR profile suppression seen in data. In TNG, inside-out quenching is due to the supermassive black hole (SMBH) feedback model operating at low accretion rates.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa1169
