An official website of the United States government
Abstract Active galactic nuclei (AGNs) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) andSimbacosmological hydrodynamic simulations to show that the low-redshift Lyαforest can provide constraints on the impact of AGN feedback. We show that TNG overpredicts the number density of absorbers at column densitiesNHI< 1014cm−2compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducialSimbamodel, which includes AGN jet feedback, is the preferred fit to the observed CDD of thez= 0.1 Lyαforest across 5 orders of magnitude in column density. We show that theSimbaresults with jets produce a quantitatively better fit to the observational data than theSimbaresults without jets, even when the ultraviolet background is left as a free parameter. AGN jets inSimbaare high speed, collimated, weakly interacting with the interstellar medium (via brief hydrodynamic decoupling), and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the intergalactic medium when compared to TNG’s kinetic feedback mode, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low-redshift Lyαforest provides plausible evidence for long-range AGN jet feedback.
more »
« less
An Exploration of AGN and Stellar Feedback Effects in the Intergalactic Medium via the Low-redshift Lyα Forest
Abstract We explore the role of galactic feedback on the low-redshift Lyα(Lyα) forest (z≲ 2) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba subgrid models. We find that both AGN and stellar feedback in Simba play a role in setting the Lyαforest column density distribution function (CDD) and the Doppler width (b-value) distribution. The Simba AGN jet feedback mode is able to efficiently transport energy out to the diffuse IGM, causing changes in the shape and normalization of the CDD and a broadening of theb-value distribution. We find that stellar feedback plays a prominent role in regulating supermassive black hole growth and feedback, highlighting the importance of constraining stellar and AGN feedback simultaneously. In IllustrisTNG, the AGN feedback variations explored in CAMELS do not affect the Lyαforest, but varying the stellar feedback model does produce subtle changes. Our results imply that the low-zLyαforest can be sensitive to changes in the ultraviolet background, stellar and black hole feedback, and that AGN jet feedback in particular can have a strong effect on the thermal state of the IGM.
more »
« less
- PAR ID:
- 10473123
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astronomical Journal
- Volume:
- 166
- Issue:
- 6
- ISSN:
- 0004-6256
- Format(s):
- Medium: X Size: Article No. 228
- Size(s):
- Article No. 228
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Most diffuse baryons, including the circumgalactic medium (CGM) surrounding galaxies and the intergalactic medium (IGM) in the cosmic web, remain unmeasured and unconstrained. Fast radio bursts (FRBs) offer an unparalleled method to measure the electron dispersion measures (DMs) of ionized baryons. Their distribution can resolve the missing baryon problem and constrain the history of feedback theorized to impart significant energy to the CGM and IGM. We analyze the Cosmology and Astrophysics with Machine Learning Simulations using three suites, IllustrisTNG, SIMBA, and Astrid, each varying six parameters (two cosmological and four astrophysical feedback), for a total of 183 distinct simulation models. We find significantly different predictions between the fiducial models of the suites owing to their different implementations of feedback. SIMBA exhibits the strongest feedback, leading to the smoothest distribution of baryons and reducing the sight-line-to-sight-line variance in DMs betweenz= 0 and 1. Astrid has the weakest feedback and the largest variance. We calculate FRB CGM measurements as a function of galaxy impact parameter, with SIMBA showing the weakest DMs due to aggressive active galactic nucleus (AGN) feedback and Astrid the strongest. Within each suite, the largest differences are due to varying AGN feedback. IllustrisTNG shows the most sensitivity to supernova feedback, but this is due to the change in the AGN feedback strengths, demonstrating that black holes, not stars, are most capable of redistributing baryons in the IGM and CGM. We compare our statistics directly to recent observations, paving the way for the use of FRBs to constrain the physics of galaxy formation and evolution.more » « less
-
Abstract The absorption by neutral hydrogen in the intergalactic medium (IGM) produces the Ly α forest in the spectra of quasars. The Ly α forest absorbers have a broad distribution of neutral hydrogen column density N H I and Doppler b parameter. The narrowest Ly α absorption lines (of lowest b ) with neutral hydrogen column density above ∼10 13 cm −2 are dominated by thermal broadening, which can be used to constrain the thermal state of the IGM. Here we constrain the temperature-density relation T = T 0 ( ρ / ρ ¯ ) γ − 1 of the IGM at 1.6 < z < 3.6 by using N H I and b parameters measured from 24 high-resolution and high-signal-to-noise quasar spectra and by employing an analytic model to model the N H I -dependent low- b cutoff in the b distribution. In each N H I bin, the b cutoff is estimated using two methods, one non-parametric method from computing the cumulative b distribution and a parametric method from fitting the full b distribution. We find that the IGM temperature T 0 at the mean gas density ρ ¯ shows a peak of ∼1.5 × 10 4 K at z ∼ 2.7–2.9. At redshift higher than this, the index γ approximately remains constant, and it starts to increase toward lower redshifts. The evolution in both parameters is in good agreement with constraints from completely different approaches, which signals that He ii reionization completes around z ∼ 3.more » « less
-
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
-
ABSTRACT We quantify the cosmological spread of baryons relative to their initial neighbouring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighbouring distribution owing to chaotic gravitational dynamics on spatial scales comparable to their host dark matter halo. In contrast, gas in hydrodynamic simulations spreads much further from the initial neighbouring dark matter owing to feedback from supernovae (SNe) and active galactic nuclei (AGN). We show that large-scale baryon spread is very sensitive to model implementation details, with the fiducial simba model spreading ∼40 per cent of baryons >1 Mpc away compared to ∼10 per cent for the IllustrisTNG and astrid models. Increasing the efficiency of AGN-driven outflows greatly increases baryon spread while increasing the strength of SNe-driven winds can decrease spreading due to non-linear coupling of stellar and AGN feedback. We compare total matter power spectra between hydrodynamic and paired N-body simulations and demonstrate that the baryonic spread metric broadly captures the global impact of feedback on matter clustering over variations of cosmological and astrophysical parameters, initial conditions, and (to a lesser extent) galaxy formation models. Using symbolic regression, we find a function that reproduces the suppression of power by feedback as a function of wave number (k) and baryonic spread up to $$k \sim 10\, h$$ Mpc−1 in SIMBA while highlighting the challenge of developing models robust to variations in galaxy formation physics implementation.more » « less
