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
Galaxy cold gas contents in modern cosmological hydrodynamic simulations
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
- NSF-PAR ID:
- 10390149
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 497
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 146 to 166
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
ABSTRACT In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a comparative analysis of the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $z\approx 0$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($M_{\rm 200c}\lesssim 10^{11.5}\, \mathrm{M}_{\odot }$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $2\!-\!3\times R_{\rm 200c}$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the interstellar medium, recycle within the circumgalactic medium (within $R_{\rm 200c}$). We find that active galactic nucleus (AGN)-driven outflows in SIMBA are notably potent, reaching several times $R_{\rm 200c}$ even at halo masses up to $M_{\rm 200c}\approx 10^{13.5}\, \mathrm{M}_{\odot }$. In both TNG and EAGLE, AGN feedback can eject gas beyond $R_{\rm 200c}$ at this mass scale, but seldom beyond $2\!-\!3\times R_{\rm 200c}$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $R_{\rm 200c}$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform subgrid feedback models in the next generation of simulations.
-
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. -
ABSTRACT We present a study of the stellar host galaxy, CO (1–0) molecular gas distribution and AGN emission on 50–500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300–500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5–20 kpc in size. The CO (1–0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1–0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation. We find CO–H2 conversion factors of around αCO = 1.5 ± 0.5 (K km s−1 pc2)−1, molecular gas masses of >3 × 1010 M⊙, dynamical masses of ∼1011 M⊙, and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa1894
