Numerical simulations of winds driven by radiation force from the corona above a thin disk

Observations show that winds can be driven from the innermost region (inside a 50 Schwarschild radius) of a thin disk. It is interesting to study the winds launched from the innermost region. A hot corona above the black hole (BH) thin disk is irradiated by the disk. We perform two-dimensional hydrodynamical simulations to study the winds driven by radiation force from the corona in the innermost regions. The hard X-ray spectrum from active galactic nuclei (AGNs) suggests that the corona temperature is about $10^9$ K, so that we mainly analyze the properties of winds (or outflows) from the $10^9$ K corona. The disk luminosity plays an important role in driving the outflows. The more luminous the disk, the stronger the outflows. Mass outflow rate ($\dot{M}_{\rm out}$) at a 90 Schwarschild radius depends on disk luminosity, which can be described as $\dot{M}_{\rm out}\propto 10^{3.3 \Gamma}$ ($\Gamma$ is the ratio of the disk luminosity to the Eddington luminosity). In the case of high luminosity (e.g. $\Gamma=0.75$), the supersonic outflows with maximum speed $1.0 \times 10^4$ Km s$^{-1}$ are launched at $\sim17^{o}$ --$30^{o}$ and $\sim50^{o}$ --$80^{o}$ away from the pole axis. The Bernoulli parameter keeps increasing with the outward propagation of outflows. The radiation force keeps accelerating the outflows when outflows move outward. Therefore, we can expect the outflows to escape from the BH gravity and go to the galactic scale. The interaction between outflows and interstellar medium may be an important AGN feedback process.Comment: 9 pages, 12 figures, accepted for publication in Ap

Two dimensional numerical simulations of Supercritical Accretion Flows revisited

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{\theta\phi}$ component of the viscous stress and consider various values of the viscous parameter $\alpha$. We find that when $T_{\theta\phi}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $\alpha$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $\rho(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $\alpha$ ranging from $\alpha=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $\alpha$, the flow is marginally stable (when $\alpha$ is small) or unstable (when $\alpha$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.Comment: 16 pages, 13 figures, accepted for publication in Ap

Thermal wind from hot accretion flows at large radii

We study slowly rotating accretion flow at parsec and sub-parsec scale irradiated by a low luminosity active galactic nuclei. We take into account the Compton heating, photoionization heating by the central X-rays. The bremsstrahlung cooling, recombination and line cooling are also included. We find that due to the Compton heating, wind can be thermally driven. The power of wind is in the range $(10^{-6}-10^{-3}) L_{\rm Edd}$, with $L_{\rm Edd}$ being the Eddington luminosity. The mass flux of wind is in the range $(0.01-1) \dot M_{\rm Edd}$ ($\dot M_{\rm Edd}= L_{\rm Edd}/0.1c^2$ is the Eddington accretion rate, $c$ is speed of light). We define the wind generation efficiency as $\epsilon = P_W/\dot {M}_{\rm BH}c^2$, with $P_W$ being wind power, $\dot M_{\rm BH}$ being the mass accretion rate onto the black hole. $\epsilon$ lies in the rage $10^{-4}-1.18$. Wind production efficiency decreases with increasing mass accretion rate. The possible role of the thermally driven wind in the active galactic feedback is briefly discussed.Comment: 9 pages, 6 figures, accepted by MNRA

What is the real accretion rate onto a black hole for low angular momentum accretion?

Mass accretion rate is a key parameter in accretion disk theory. It determines black hole accretion mode. In large scale cosmological simulations studying galaxy formation and evolution, Bondi radius can at most be marginally resolved. In those simulations, Bondi accretion formula is always used to estimate black hole accretion rate. Bondi solution is too simple, which cannot represent the real accretion process. We perform simulations of hot accretion flow at parsec scale irradiated by a low luminosity active galactic nucleus (AGN). We perform 77 simulations with varying density and temperature at outer boundary (10 parsec). Our purpose is to find a formula to calculate real black hole accretion rate based on the gas density and temperature at parsec scale. We define Eddington accretion rate to be $\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ been Eddington luminosity and speed of light respectively. We set black hole mass to be $10^8M_{\odot}$, $M_{\odot}$ is solar mass. We find that black hole accretion rate can be expressed as $\frac{\dot M_{\rm BH} }{\dot M_{\rm Edd}}=10^{-3.11} (\rho_0/10^{-22}{\rm g cm}^{-3})^{1.36} (T_0/10^7 {\rm K})^{-1.9}$, with $\rho_0$ and $T_0$ being density and temperature at parsec scale, respectively. We find the formula can accurately predict the luminosity of observed low-luminosity AGNs (with black hole mass $\sim 10^8M_{\odot}$). This formula can be used in the sub-grid models in large scale cosmological simulations with a black hole mass of $\sim 10^8M_{\odot}$.Comment: 12 pages, 9 figures, accepted by MNRA