skip to main content

Explore Research Products in the NSF-PAR

Title: Dust dynamics in AGN winds: a new mechanism for multiwavelength AGN variability
ABSTRACT

Partial dust obscuration in active galactic nuclei (AGNs) has been proposed as a potential explanation for some cases of AGN variability. The dust–gas mixture present in AGN tori is accelerated by radiation pressure, leading to the launching of an AGN wind. Dust under these conditions has been shown to be unstable to a generic class of fast-growing resonant drag instabilities (RDIs). In this work, we present the first numerical simulations of radiation-driven outflows that explicitly include dust dynamics in conditions resembling AGN winds. We investigate the implications of RDIs on the torus morphology, AGN variability, and the ability of radiation to effectively launch a wind. We find that the RDIs rapidly develop, reaching saturation at times much shorter than the global time-scales of the outflows, resulting in the formation of filamentary structure on box-size scales with strong dust clumping and super-Alfvénic velocity dispersions. The instabilities lead to fluctuations in dust opacity and gas column density of 10–20  per cent when integrated along mock observed lines of sight to the quasar accretion disc. These fluctuations occur over year to decade time-scales and exhibit a red-noise power spectrum commonly observed for AGNs. Additionally, we find that the radiation effectively couples with the dust–gas mixture, launching highly supersonic winds that entrain 70–90  per cent of the gas, with a factor of ≲3 photon momentum loss relative to the predicted multiple-scattering momentum loading rate. Therefore, our findings suggest that RDIs play an important role in driving the clumpy nature of AGN tori and generating AGN variability consistent with observations.

  more » « less
Award ID(s):
2009234
NSF-PAR ID:
10450094
Author(s) / Creator(s):
;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
525
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2668-2689
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT Observations of emission lines in active galactic nuclei (AGNs) often find fast (∼1000 km s−1) outflows extending to kiloparsec scales, seen in ionized, neutral atomic and molecular gas. In this work we present radiative transfer calculations of emission lines in hydrodynamic simulations of AGN outflows driven by a hot wind bubble, including non-equilibrium chemistry, to explore how these lines trace the physical properties of the multiphase outflow. We find that the hot bubble compresses the line-emitting gas, resulting in higher pressures than in the ambient interstellar medium or that would be produced by the AGN radiation pressure. This implies that observed emission line ratios such as [O iv]$_{25 \, \rm {\mu m}}$ / [Ne ii]$_{12 \, \rm {\mu m}}$, [Ne v]$_{14 \, \rm {\mu m}}$ / [Ne ii]$_{12 \, \rm {\mu m}}$, and [N iii]$_{57 \, \rm {\mu m}}$ / [N ii]$_{122 \, \rm {\mu m}}$ constrain the presence of the bubble and hence the outflow driving mechanism. However, the line-emitting gas is under-pressurized compared to the hot bubble itself, and much of the line emission arises from gas that is out of pressure, thermal and/or chemical equilibrium. Our results thus suggest that assuming equilibrium conditions, as commonly done in AGN line emission models, is not justified if a hot wind bubble is present. We also find that ≳50 per cent of the mass outflow rate, momentum flux, and kinetic energy flux of the outflow are traced by lines such as [N ii]$_{122 \, \rm {\mu m}}$ and [Ne iii]$_{15 \, \rm {\mu m}}$ (produced in the 10$^{4} \, \rm {K}$ phase) and [C ii]$_{158 \, \rm {\mu m}}$ (produced in the transition from 10$^{4} \, \rm {K}$ to 100 K). 
    more » « less
  2. ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    The ultraviolet (UV) bright accretion disc in active galactic nuclei (AGNs) should give rise to line driving, producing a powerful wind that may play an important role in AGN feedback as well as in producing structures like the broad-line region. However, coupled radiation-hydrodynamic codes are complex and expensive, so we calculate the winds instead using a non-hydrodynamical approach (the qwind framework). The original qwind model assumed the initial conditions in the wind, and had only simple radiation transport. Here, we present an improved version that derives the wind initial conditions and has significantly improved ray tracing to calculate the wind absorption self-consistently, given the extended nature of the UV emission. We also correct the radiation flux for relativistic effects and assess the impact of this on the wind velocity. These changes mean the model is more physical, so its predictions are more robust. We find that, even when accounting for relativistic effects, winds can regularly achieve velocities ≃(0.1−0.5)c, and carry mass-loss rates that can be up to 80 per cent of the accreted mass for black hole masses of 107−9 M⊙, and mass accretion rates of 50 per cent of the Eddington rate. Overall, the ratio of kinetic power carried by the wind to bolometric luminosity increases with mass accretion rate at a given black hole mass, unlike the constant fraction generally assumed in current cosmological simulations that include AGN feedback. The updated code, qwind3, is publicly available in GitHub.1

     
    more » « less
  3. ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    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
  4. ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    E+A galaxies are believed to be a short phase connecting major merger ultraluminous infrared galaxies (ULIRGs) with red and dead elliptical galaxies. Their optical spectrum suggests a massive starburst that was quenched abruptly, and their bulge-dominated morphologies with tidal tails suggest that they are merger remnants. Active galactic nucleus (AGN)-driven winds are believed to be one of the processes responsible for the sudden quenching of star formation and for the expulsion and/or destruction of the remaining molecular gas. Little is known about AGN-driven winds in this short-lived phase. In this paper, we present the first and unique sample of post-starburst galaxy candidates with AGNs that show indications of ionized outflows in their optical emission lines. Using Infrared Astronomical Satellite–far infrared (IRAS–FIR) observations, we study the star formation in these systems and find that many systems selected to have post-starburst signatures in their optical spectrum are in fact obscured starbursts. Using SDSS spectroscopy, we study the stationary and outflowing ionized gas. We also detect neutral gas outflows in 40 per cent of the sources with mass outflow rates 10–100 times more massive than in the ionized phase. The mean mass outflow rate and kinetic power of the ionized outflows in our sample ($\dot{M}\sim 1\, \mathrm{M_{\odot }\, yr^{-1}}$, $\dot{E}\sim 10^{41}\, \mathrm{erg\, s}^{-1}$) are larger than those derived for active galaxies of similar AGN luminosity and stellar mass. For the neutral outflow ($\dot{M}\sim 10\, \mathrm{M_{\odot }\, yr^{-1}}$, $\dot{E}\sim 10^{42}\, \mathrm{erg\, s}^{-1}$), their mean is smaller than that observed in (U)LIRGs with and without AGN.

     
    more » « less
  5. ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    In dusty cool-star outflow or ejection events around asymptotic giant branch (AGB) or R Coronae Borealis or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. It has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulations of the RDIs operating on a spectrum of dust grain sizes subject to radiative acceleration (allowing for different grain optical properties), coupled to the gas with a realistic drag law, including or excluding the effects of magnetic fields and charged grains, and calculate for the first time how the RDIs could contribute to observed variability. We show that the RDIs naturally produce significant variations (spatially and temporally) ($\sim 10\!-\!20{{\ \rm per\ cent}}$ 1 σ-level) in the extinction, corresponding to $\sim 0.1\!-\!1\,$mag level in the stellar types above, on time-scales of order months to a year. The fluctuations are surprisingly robust to the assumed size of the source as they are dominated by large-scale modes, which also means their spatial structure could be resolved in some nearby systems. We also quantify how this produces variations in the line-of-sight grain size-distribution. All of these variations are similar to those observed, suggesting that the RDIs may play a key role driving observed spatial and temporal variability in dust extinction within dusty outflow/ejection events around cool stars. We further propose that the measured variations in grain sizes could directly be used to identify the presence of the RDIs in close by systems with observations.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email