skip to main content

Title: 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 inclusion of dust destruction mechanisms, like thermal sputtering, leads to the rapid destruction of dust grains within the outflow, reducing its properties below the initial NFW halo. We conclude that radiative AGN feedback can drive outflows, but a thorough numerical and physical treatment is required to assess its true impact.  more » « less
Award ID(s):
1814053 1814259
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
1143 to 1164
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    Radio-loud active galactic nuclei (RLAGNs) are a unique AGN population and were thought to be preferentially associated with supermassive black holes (SMBHs) at low accretion rates. They could impact the host galaxy evolution by expelling cold gas through the jet-mode feedback. In this work, we studied CO(6−5) line emission and continuum emission in a high-redshift radio galaxy, MRC 0152−209, at z = 1.92 using ALMA (Atacama Large Millimeter/submillimeter Array) up to a 0.024″ resolution (corresponding to ∼200 pc at z = 1.92). This system is a starburst major merger comprising two galaxies: the north-west (NW) galaxy hosting the RLAGN with jet kinetic power Ljet ≳ 2 × 1046  erg s−1 and the other galaxy to the south-east (SE). Based on the spectral energy distribution fitting for the entire system (NW+SE galaxies), we find an AGN bolometric luminosity LAGN, bol ∼ 3 × 1046  erg s−1 with a lower limit of ∼0.9 × 1046  erg s−1 for the RLAGN. We estimate the black hole mass through MBH–M⋆ scaling relations and find an Eddington ratio of λEdd ∼ 0.07–4 conservatively by adopting the lower limit of LAGN, bol and considering the dispersion of the scaling relation. These results suggest that the RLAGN is radiatively efficient and the powerful jets could be launched from a super-Eddington accretion disc. ALMA Cycle 6 observations further reveal a massive (${M}_\mathrm{H_2}=(1.1-2.3)\times 10^9\ \rm M_\odot$), compact (∼500 pc), and monopolar molecular outflow perpendicular to the jet axis. The corresponding mass outflow rate ($1200^{+300}_{-300}-2600^{+600}_{-600}\ \mathrm{M_\odot }\ \rm yr^{-1}$) is comparable with the star formation rate of at least $\sim 2100\ \mathrm{M_\odot }\ \rm yr^{-1}$. Depending on the outflowing molecular gas mass, the outflow kinetic power/LAGN, bol ratio of ∼0.008–0.02, and momentum boost factor of ∼3–24 agree with a radiative-mode AGN feedback scenario. On the other hand, the jets can also drive the molecular outflow within its lifetime of ∼2 × 105 yr without additional energy supply from AGN radiation. The jet-mode feedback is then capable of removing all cold gas from the host galaxy through the long-term, episodic launching of jets. Our study reveals a unique object where starburst activity, powerful jets, and rapid BH growth co-exist, which may represent a fundamental stage of AGN-host galaxy co-evolution.

    more » « less

    Radiation pressure-driven outflows from luminous accreting supermassive black holes are an important part of active galactic nucleus (AGN) feedback. The effective Eddington limit, based on absorption of radiation by dust, not electron scattering, is revealed in the plane of AGN absorption column density NH as a function of Eddington fraction λEdd = Lbol/LEdd, where a lack of objects is seen in the region where the effective limit is exceeded. Here, we conduct radiation simulation using the cloudy code to deduce the radiative force applied on to dusty gas at the nucleus and compare to the gravitational force to reveal the outflow region and its boundary with long-lived absorption clouds. We also investigate how the outflow condition is affected by various AGN and dust properties and distribution. As expected, the dust abundance has the largest effect on the NH–λEdd diagram since the higher the abundance, the more effective the radiative feedback, while the impact of the inner radius of the dusty gas shell, the shell width, and the AGN spectral shape are relatively negligible. The presence of other central masses, such as a nuclear star cluster, can also make the feedback less effective. The AGN spectral energy distribution depends on the mass of the black hole and its spin. Though the effects of the AGN spectral energy distribution on the diagram are relatively small, the fraction of ionizing ultraviolet photons from the blackbody accretion disc is affected more by black hole mass than spin, and can influence the efficiency of radiation pressure.

    more » « less
  3. null (Ed.)
    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. We then compare our predictions to the observed size–mass relationship and the Σ1−M⋆ and Σe−M⋆ relationships. At z ≳ 2, the simulated massive galaxies are in general agreement with observational scaling relations. At z ≲ 2, they evolve to become too compact but still star forming, in the stellar mass and redshift regime where many of them should be quenched. Our results suggest that some additional source of feedback, such as AGN-driven outflows, is necessary in order to decrease the central densities of the simulated massive galaxies to bring them into agreement with observations at z ≲ 2. 
    more » « less
  4. 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

    Without active galactic nucleus (AGN) feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z ≲ 2. In this paper, we perform high-resolution re-simulations of a massive ($M_{\star }\sim 10^{11}\, \rm {{\rm M}_{\odot }}$) galaxy at z ∼ 2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multiwavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ∼1 to $\sim 0.1\, \rm {kpc}$ in $\lesssim 40\, \rm {Myr}$ when AGN winds are not included, but increases to $\sim 2\, \rm {kpc}$ when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over $35\, \rm {Myr}$, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations.

    more » « less