skip to main content


Title: Toward Horizon-scale Accretion onto Supermassive Black Holes in Elliptical Galaxies
Abstract We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new GPU-accelerated version of the Athena++ AMR code, and they span more than six orders of magnitude in radius, reaching scales similar to that of the black hole horizon. The key physical ingredients are radiative cooling and a phenomenological heating model. We find that the accretion flow takes the form of multiphase gas at radii less than about a kpc. The cold gas accretion includes two dynamically distinct stages: the typical disk stage in which the cold gas resides in a rotationally supported disk, and relatively rare chaotic stages (≲10% of the time) in which the cold gas inflows via chaotic streams. Though cold gas accretion dominates the time-averaged accretion rate at intermediate radii, accretion at the smallest radii is dominated by hot virialized gas at most times. The accretion rate scales with radius as M ̇ ∝ r 1 / 2 when hot gas dominates, and we obtain M ̇ ≃ 10 − 4 – 10 − 3 M ⊙ yr − 1 near the event horizon, similar to what is inferred from EHT observations. The orientation of the cold gas disk can differ significantly on different spatial scales. We propose a subgrid model for accretion in lower-resolution simulations in which the hot gas accretion rate is suppressed relative to the Bondi rate by ∼ ( r g / r Bondi ) 1 / 2 . Our results can also provide more realistic initial conditions for simulations of black hole accretion at the event horizon scale.  more » « less
Award ID(s):
2107872
NSF-PAR ID:
10459080
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
946
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
26
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process, which is suitable for analytical and numerical studies. With radiusrscaled to Schwarzschild units and coronal mass accretion rateṁcto Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and an active galactic nucleus. The corona solution consists of two power-law segments separated at a break radiusrb∼ 103(α/0.3)−2, whereαis the viscosity parameter. Gas evaporates from the disk to the corona forr>rb, and condenses back forr<rb. Atrb,ṁcreaches its maximum,ṁc,max0.02(α/0.3)3. If atrrbthe thin disk accretes withṁd<ṁc,max, then the disk evaporates fully before reachingrb, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model that includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, andrbremains unchanged. This model predicts strong coronal winds forr>rb, and aT∼ 5 × 108K Compton-cooled corona forr<rb. Two-temperature effects are ignored, but may be important at small radii.

     
    more » « less
  2. Abstract X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up ≲10% of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole (SMBH), SMBH merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the SMBH sphere of influence of radius R B to the SMBH of gravitational radius R g at the largest scale separation, R B / R g = 10 3 , to date. Initially, our axisymmetric system of constant-density hot gas contains a weak vertical magnetic field and rotates in the equatorial plane of a rapidly spinning SMBH. We seed the gas with small-scale 2% level pressure perturbations. Infalling gas forms an accretion disk, and the SMBH launches relativistically magnetized collimated jets reaching well outside R B . Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the SMBH, becomes dynamically important, and leads to a magnetically arrested disk state. The SMBH accretes at 2% of the Bondi rate ( M ̇ ≃ 2.4 × 10 − 3 M ⊙ yr − 1 for M87*) and launches twin jets at η = 150% efficiency. These jets are powerful enough ( P jets ≃ 2 × 10 44 erg s −1 ) to escape along the SMBH spin axis and end the short-lived intermittent jet state, whose transient nature can account for the rarity of XRGs. 
    more » « less
  3. Abstract

    Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius,RB= 103Rg. The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideRBand suppress BH accretion to 1.5% of the Bondi rate,ṀB. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% ofṀBreach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2%ṀB. Because the BAD and RAD states tangle up the jets and destroy them well insideRB, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.

     
    more » « less
  4. Abstract

    Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius ofRB≈ 2 × 105GM/c2, whereMis the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼r−1rather thanr−3/2and the accretion rateṀis consequently suppressed by over 2 orders of magnitude below the Bondi rateṀB. We find continuous energy feedback from the accretion flow to the external medium at a level of102Ṁc25×105ṀBc2. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

     
    more » « less
  5. null (Ed.)
    ABSTRACT We examine the 1.4 GHz radio luminosities of galaxies arising from star formation and active galactic nuclei (AGNs) within the state-of-the-art cosmological hydrodynamic simulation Simba. Simba grows black holes via gravitational torque limited accretion from cold gas and Bondi accretion from hot gas, and employs AGN feedback including jets at low Eddington ratios. We define a population of radio loud AGNs (RLAGNs) based on the presence of ongoing jet feedback. Within RLAGN, we define high and low excitation radio galaxies (HERGs and LERGs) based on their dominant mode of black hole accretion: torque limited accretion representing feeding from a cold disc, or Bondi representing advection-dominated accretion from a hot medium. Simba predicts good agreement with the observed radio luminosity function (RLF) and its evolution, overall as well as separately for HERGs and LERGs. Quiescent galaxies with AGN-dominated radio flux dominate the RLF at $\gtrsim 10^{22-23}$ W Hz−1, while star formation dominates at lower radio powers. Overall, RLAGNs have higher black hole accretion rates and lower star formation rates than non-RLAGN at a given stellar mass or velocity dispersion, but have similar black hole masses. Simba predicts an LERG number density of 8.53 Mpc−3, ∼10× higher than for HERGs, broadly as observed. While LERGs dominate among most massive galaxies with the largest black holes and HERGs dominate at high specific star formation rates, they otherwise largely populate similar-sized dark matter haloes and have similar host galaxy properties. Simba thus predicts that deeper radio surveys will reveal an increasing overlap between the host galaxy demographics of HERGs and LERGs. 
    more » « less