Outflows from inflows: the nature of Bondi-like accretion

The classic Bondi solution remains a common starting point both for studying black hole growth across cosmic time in cosmological simulations and for smaller scale simulations of active galactic nuclei (AGN) feedback. In nature, however, there will be inhomogeneous distributions of rotational velocity and density along the outer radius (R-o) marking the sphere of influence of a black hole. While there have been many studies of how the Bondi solution changes with a prescribed angular momentum boundary condition, they have all assumed a constant density at R-o. In this Letter, we show that a non-uniform density at R-o causes a meridional flow and due to conservation of angular momentum, the Bondi solution qualitatively changes into an inflow-outflow solution. Using physical arguments, we analytically identify the critical logarithmic density gradient vertical bar partial derivative ln rho/partial derivative theta vertical bar above which this change of the solution occurs. For realistic R-o, this critical gradient is less than 0.01 and tends to 0 as R-o -> infinity. We show using numerical simulations that, unlike for solutions with an imposed rotational velocity, the accretion rate for solutions under an inhomogeneous density boundary condition remains constant at nearly the Bondi rate (M) over dot(B), while the outflow rate can greatly exceed (M) over dot(B).