Thermal and dynamical properties of gas accreting onto a supermassive black hole in an active galactic nucleus

Abstract

We study stability of gas accretion in active galactic nuclei (AGNs). Our grid-based simulations cover a radial range from 0.1 to 200 pc, which may enable linking the galactic/cosmological simulations with small-scale black hole (BH) accretion models within a few hundreds of Schwarzschild radii. Here, as in previous studies by our group, we include gas radiative cooling as well as heating by a sub-Eddington X-ray source near the central supermassive BH of 108 M⊙. Our theoretical estimates and simulations show that for the X-ray luminosity, LX ∼ 0.008 L Edd, the gas is thermally and convectively unstable within the computational domain. In the simulations, we observe that very tiny fluctuations in an initially smooth, spherically symmetric, accretion flow, grow first linearly and then nonlinearly. Consequently, an initially one-phase flow relatively quickly transitions into a two-phase/cold-hot accretion flow. For LX = 0.015 L Edd or higher, the cold clouds continue to accrete but in some regions of the hot phase, the gas starts to move outward. For LX < 0.015 L Edd, the cold phase contribution to the total mass accretion rate only moderately dominates over the hot phase contribution. This result might have some consequences for cosmological simulations of the so-called AGN feedback problem. Our simulations confirm the previous results of Barai et al. who used smoothed particle hydrodynamic (SPH) simulations to tackle the same problem. Here, however, because we use a grid-based code to solve equations in one dimension and two dimensions, we are able to follow the gas dynamics at much higher spacial resolution and for longer time compared with the three-dimensional SPH simulations. One of the new features revealed by our simulations is that the cold condensations in the accretion flow initially form long filaments, but at the later times, those filaments may break into smaller clouds advected outward within the hot outflow. Therefore, these simulations may serve as an attractive model for the so-called narrow-line region in AGNs. © 2013. The American Astronomical Society. All rights reserved.