Rotating Accretion Flows: From Infinity to the Black Hole

Accretion onto a supermassive black hole of a rotating inflow is a particularly difficult problem to study because of the wide range of length scales involved. There have been broadly utilized analytic and numerical treatments of the global properties of accretion flows, but detailed numerical simulations are required to address certain critical aspects. We use the ZEUS code to run hydrodynamical simulations of rotating, axisymmetric accretion flows with Bremsstrahlung cooling, considering solutions for which the centrifugal balance radius significantly exceeds the Schwarzschild radius, with and without viscous angular momentum transport. Infalling gas is followed from well beyond the Bondi radius down to the vicinity of the black hole. We produce a continuum of solutions with respect to the single parameter Mdot_Bondi/Mdot_Edd, and there is a sharp transition between two general classes of solutions at an Eddington ratio of Mdot_Bondi/Mdot_Edd ~ few x 10^(-2). Our high inflow solutions are very similar to the standard Shakura & Sunyaev (1973) results. But our low inflow results are to zeroth order the stationary Papaloizou and Pringle (1984) solution, which has no accretion. To next order in the small, assumed viscosity they show circulation, with disk and conical wind outflows almost balancing inflow. These solutions are characterized by hot, vertically extended disks, and net accretion proceeds at an extremely low rate, only of order alpha times the inflow rate. Our simulations have converged with respect to spatial resolution and temporal duration, and they do not depend strongly on our choice of boundary conditions.

💡 Research Summary

The paper tackles the long‑standing problem of how rotating gas accretes onto a super‑massive black hole (SMBH) when the centrifugal balance radius is much larger than the Schwarzschild radius. Using the ZEUS hydrodynamics code, the authors perform axisymmetric two‑dimensional simulations that follow gas from well beyond the Bondi radius down to the immediate vicinity of the event horizon. The gas is allowed to cool via bremsstrahlung emission, and angular momentum transport is modeled with an α‑prescription for viscosity; runs are performed both with and without this viscous term.

A single dimensionless control parameter, the ratio of the Bondi accretion rate to the Eddington accretion rate (ṁ_Bondi/ṁ_Edd), is varied to generate a continuous family of solutions. The key result is a sharp transition at an Eddington ratio of order a few × 10⁻². Below this threshold the flow resembles the stationary Papaloizou‑Pringle (1984) torus: the centrifugal force nearly balances gravity, producing a hot, vertically thick configuration that, in the zeroth‑order limit, does not accrete at all. When a small viscosity (α ≪ 1) is introduced, the torus develops a weak circulation pattern and a conical wind that almost cancels the inflow. The net mass accretion onto the black hole is then suppressed to roughly α × ṁ_inflow, i.e., an extremely low fraction of the material supplied at large radii.

Above the critical ṁ_Bondi/ṁ_Edd, the flow switches to a thin, radiatively efficient disk that is essentially identical to the classic Shakura‑Sunyaev (1973) solution. Viscous angular‑momentum transport efficiently drives matter inward, and the accretion rate onto the black hole approaches the Bondi supply rate. The disk is cool, geometrically thin, and the associated wind is weak compared with the inflow.

The authors verify numerical convergence by increasing spatial resolution and extending the integration time; key diagnostics (accretion rate, disk scale height, wind velocity) change by less than a few percent. They also test the sensitivity to boundary conditions by varying the outer inflow properties and the inner absorbing radius, finding that the qualitative behavior of both low‑ and high‑Ṁ regimes is robust.

From a physical perspective, the study provides a clean explanation for the dichotomy observed in active galactic nuclei. Low‑luminosity systems (e.g., LINERs) can be interpreted as operating in the low‑Ṁ regime where a hot, thick torus plus weak viscous circulation limits the net inflow to a tiny fraction of the available gas. High‑luminosity quasars and Seyferts correspond to the high‑Ṁ regime where a thin, efficient disk transports most of the supplied mass onto the SMBH. Importantly, the work shows that even in the absence of magnetohydrodynamic effects, a modest viscous stress is sufficient to generate outflows and circulation that dramatically alter the accretion efficiency.

The paper concludes that the transition between these two modes is primarily governed by the dimensionless accretion rate, and that the results are insensitive to the details of the numerical setup. Future extensions should incorporate three‑dimensional MHD, radiation transfer, and direct comparison with multi‑wavelength observations to test the applicability of the two‑mode picture to real AGN populations.