No Arabic abstract
We present three-dimensional calculations of spherically symmetric Bondi accretion onto a stationary supermassive black hole (SMBH) of mass $10^{8}$ $M_{odot}$ within a radial range of $0.02-10$ pc, using a modified version of the smoothed particle hydrodynamics (SPH) gad sp code, which ensures approximate first-order consistency (i.e., second-order accuracy) for the particle approximation. First-order consistency is restored by allowing the number of neighbours, $n_{rm neigh}$, and the smoothing length, $h$, to vary with the total number of particles, $N$, such that the asymptotic limits $n_{rm neigh}toinfty$ and $hto 0$ hold as $Ntoinfty$. The ability of the method to reproduce the isothermal ($gamma =1$) and adiabatic ($gamma =5/3$) Bondi accretion is investigated with increased spatial resolution. In particular, for the isothermal models the numerical radial profiles closely match the Bondi solution, except near the accretor, where the density and radial velocity are slightly underestimated. However, as $n_{rm neigh}$ is increased and $h$ is decreased, the calculations approach first-order consistency and the deviations from the Bondi solution decrease. The density and radial velocity profiles for the adiabatic models are qualitatively similar to those for the isothermal Bondi accretion. Steady-state Bondi accretion is reproduced by the highly resolved consistent models with a percent relative error of $lesssim 1$% for $gamma =1$ and $sim 9$% for $gamma =5/3$, with the adiabatic accretion taking longer than the isothermal case to reach steady flow. The performance of the method is assessed by comparing the results with those obtained using the standard Gadget and the Gizmo codes.
In this paper, we present the classical Bondi accretion theory for the case of non-isothermal accretion processes onto a supermassive black hole (SMBH), including the effects of X-ray heating and the radiation force due to electron scattering and spectral lines. The radiation field is calculated by considering an optically thick, geometrically thin, standard accretion disk as the emitter of UV photons and a spherical central object as a source of X-ray emission. In the present analysis, the UV emission from the accretion disk is assumed to have an angular dependence, while the X-ray/central object radiation is assumed to be isotropic. This allows us to build streamlines in any angular direction we need to. The influence of both types of radiation is evaluated for different flux fractions of the X-ray and UV emissions with and without the effects of spectral line driving. We find that the radiation emitted near the SMBH interacts with the infalling matter and modifies the accretion dynamics. In the presence of line driving, a transition resembles from pure type 1 & 2 to type 5 solutions (see Fig2.1 of Frank etal. 2002), which takes place regardless of whether or not the UV emission dominates over the X-ray emission. We compute the radiative factors at which this transition occurs, and discard type 5 solution from all our models. Estimated values of the accretion radius and accretion rate in terms of the classical Bondi values are also given. The results are useful for the construction of proper initial conditions for time-dependent hydrodynamical simulations of accretion flows onto SMBH at the centre of galaxies.
Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead predict that accretion can be dominated by a stochastic, clumpy distribution of very cold molecular clouds - a departure from the hot mode accretion model - although unambiguous observational support for this prediction remains elusive. Here we report observations that reveal a cold, clumpy accretion flow towards a supermassive black hole fuel reservoir in the nucleus of the Abell 2597 Brightest Cluster Galaxy (BCG), a nearby (redshift z=0.0821) giant elliptical galaxy surrounded by a dense halo of hot plasma. Under the right conditions, thermal instabilities can precipitate from this hot gas, producing a rain of cold clouds that fall toward the galaxys centre, sustaining star formation amid a kiloparsec-scale molecular nebula that inhabits its core. The observations show that these cold clouds also fuel black hole accretion, revealing shadows cast by the molecular clouds as they move inward at about 300 kilometres per second towards the active supermassive black hole in the galaxy centre, which serves as a bright backlight. Corroborating evidence from prior observations of warmer atomic gas at extremely high spatial resolution, along with simple arguments based on geometry and probability, indicate that these clouds are within the innermost hundred parsecs of the black hole, and falling closer towards it.
We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $rsim R_{rm B}$, where the density follows $rho propto (1+R_{rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $rho propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($dot{M}propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $dot{M}_{rm B}$ The net accretion rate depends on the viscous strength, following $dot{M}/dot{M}_{rm B}propto (alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $dot{M}_{rm B}/dot{M}_{rm Edd}< 10^{-3}$ having minimal radiation cooling, where $dot{M}_{rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $dot{M}/dot{M}_{rm B} sim 10^{-3}-10^{-2}$. This rate is as low as $dot{M}/dot{M}_{rm Edd} sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.
Gas undergoing Bondi accretion onto a supermassive black hole (SMBH) becomes hotter toward smaller radii. We searched for this signature with a Chandra observation of the hot gas in NGC 3115, which optical observations show has a very massive SMBH. Our analysis suggests that we are resolving, for the first time, the accretion flow within the Bondi radius of an SMBH. We show that the temperature is rising toward the galaxy center as expected in all accretion models in which the black hole is gravitationally capturing the ambient gas. There is no hard central point source that could cause such an apparent rise in temperature. The data support that the Bondi radius is at about 4 arcsec-5 arcsec (188-235 pc), suggesting an SMBH of 2 x 10^9 M_sun that is consistent with the upper end of the optical results. The density profile within the Bondi radius has a power-law index of 1.03^{+0.23}_{-0.21} which is consistent with gas in transition from the ambient medium and the accretion flow. The accretion rate at the Bondi radius is determined to be {dot M}_B = 2.2 x 10^{-2} M_sun yr^{-1}. Thus, the accretion luminosity with 10% radiative efficiency at the Bondi radius (10^{44} erg s^{-1}) is about six orders of magnitude higher than the upper limit of the X-ray luminosity of the nucleus.
The fully analytical solution for isothermal Bondi accretion on a black hole (MBH) at the center of two-component Jaffe (1983) galaxy models is presented. In a previous work we provided the analytical expressions for the critical accretion parameter and the radial profile of the Mach number in the case of accretion on a MBH at the center of a spherically symmetric one-component Jaffe galaxy model. Here we apply this solution to galaxy models where both the stellar and total mass density distributions are described by the Jaffe profile, with different scale-lengths and masses, and to which a central MBH is added. For such galaxy models all the relevant stellar dynamical properties can also be derived analytically (Ciotti & Ziaee Lorzad 2018). In these new models the hydrodynamical and stellar dynamical properties are linked by imposing that the gas temperature is proportional to the virial temperature of the galaxy stellar component. The formulae that are provided allow to evaluate all flow properties, and are then useful for estimates of the scale-radius and the mass flow rate when modeling accretion on massive black holes at the center of galaxies. As an application, we quantify the departure from the true mass accretion rate of estimates obtained using the gas properties at various distances from the MBH, under the hypothesis of classical Bondi accretion.