Radiative AGN feedback on a moving mesh: the impact of the galactic disc and dust physics on outflow properties
ABSTRACT Feedback from accreting supermassive black holes (BHs), active galactic nuclei (AGNs), is now a cornerstone of galaxy formation models. In this work, we present radiation-hydrodynamic simulations of radiative AGN feedback using the novel arepo-rt code. A central BH emits radiation at a constant luminosity and drives an outflow via radiation pressure on dust grains. Utilizing an isolated Navarro–Frenk–White (NFW) halo we validate our set-up in the single- and multiscattering regimes, with the simulated shock front propagation in excellent agreement with the expected analytic result. For a spherically symmetric NFW halo, an examination of the simulated outflow properties with radiation collimation demonstrates a decreasing mass outflow rate and momentum flux, but increasing kinetic power and outflow velocity with decreasing opening angle. We then explore the impact of a central disc galaxy and the assumed dust model on the outflow properties. The contraction of the halo during the galaxy’s formation and modelling the production of dust grains result in a factor 100 increase in the halo’s optical depth. Radiation then couples momentum more efficiently to the gas, driving a stronger shock and producing a mass-loaded $\sim \!10^{3}\, \mathrm{M}_{\odot }\, \mathrm{yr}^{-1}$ outflow with a velocity of $\sim \!2000\, \mathrm{km}\, \mathrm{s}^{-1}$. However, the more »
Authors:
; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10104821
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
494
Issue:
1
Page Range or eLocation-ID:
1143 to 1164
ISSN:
0035-8711
1. ABSTRACT The galaxy size–stellar mass and central surface density–stellar mass relationships are fundamental observational constraints on galaxy formation models. However, inferring the physical size of a galaxy from observed stellar emission is non-trivial due to various observational effects, such as the mass-to-light ratio variations that can be caused by non-uniform stellar ages, metallicities, and dust attenuation. Consequently, forward-modelling light-based sizes from simulations is desirable. In this work, we use the skirt  dust radiative transfer code to generate synthetic observations of massive galaxies ($M_{*}\sim 10^{11}\, \rm {M_{\odot }}$ at z = 2, hosted by haloes of mass $M_{\rm {halo}}\sim 10^{12.5}\, \rm {M_{\odot }}$) from high-resolution cosmological zoom-in simulations that form part of the Feedback In Realistic Environments project. The simulations used in this paper include explicit stellar feedback but no active galactic nucleus (AGN) feedback. From each mock observation, we infer the effective radius (Re), as well as the stellar mass surface density within this radius and within $1\, \rm {kpc}$ (Σe and Σ1, respectively). We first investigate how well the intrinsic half-mass radius and stellar mass surface density can be inferred from observables. The majority of predicted sizes and surface densities are within a factor of 2 of the intrinsic values.more »
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 »
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 »
5. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factorsmore »