skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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. 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
  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 The early growth of black holes (BHs) in high-redshift galaxies is likely feedback regulated. While radiative feedback has been extensively studied, the role of mechanical feedback has received less scrutiny to date. Here, we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on 100 M⊙ BH accretion in the dense, low-metallicity gas expected in early protogalaxies. As the jet propagates, it shocks the surrounding gas forming a jet cocoon. The cocoon consists of a rapidly cooling cold phase at the interface with the background gas and an overpressured subsonic phase of reverse shock-heated gas filling the interior. We vary the background gas density and temperature, BH feedback efficiency, and the jet model. We found that the width of the jet cocoon roughly follows a scaling derived by assuming momentum conservation in the jet-propagation direction and energy conservation in the lateral directions. Depending on the assumed gas and jet properties, the cocoon either stays elongated to large radii or isotropizes before reaching the Bondi radius, forming a nearly spherical bubble. Lower jet velocities and higher background gas densities result in self-regulation to higher momentum fluxes and elongated cocoons. In all cases, the outward cocoon momentum flux balances the inward inflowing gas momentum flux near the Bondi radius, which ultimately regulates BH accretion. The time-averaged accretion rate always remains below the Bondi rate, and exceeds the Eddington rate only if the ambient medium is dense and cold, and/or the jet is weak (low velocity and mass loading). 
    more » « less
  4. null (Ed.)
    We studied the accretion disc structure in the doubly imaged lensed quasar SDSS J1339+1310 using r -band light curves and UV-visible to near-IR spectra from the first 11 observational seasons after its discovery. The 2009−2019 light curves displayed pronounced microlensing variations on different timescales, and this microlensing signal permitted us to constrain the half-light radius of the 1930 Å continuum-emitting region. Assuming an accretion disc with an axis inclined at 60° to the line of sight, we obtained log( r 1/2 /cm) = 15.4 −0.4 +0.93 . We also estimated the central black hole mass from spectroscopic data. The width of the C  IV , Mg  II , and H β emission lines, and the continuum luminosity at 1350, 3000, and 5100 Å, led to log( M BH / M ⊙ ) = 8.6 ± 0.4. Thus, hot gas responsible for the 1930 Å continuum emission is likely orbiting a 4.0 × 10 8   M ⊙ black hole at an r 1/2 of only a few tens of Schwarzschild radii. 
    more » « less
  5. Abstract Multiwavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nucleus transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole (BH). However, the formation, structure, and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ∼35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semiopaque corona reach, respectively,Ti≳ 1010K andTe≳ 5 × 108K. The system produces radiation, powerful collimated jets, and broader winds at the total energy efficiency exceeding 90%, the highest ever energy extraction efficiency from a spinning BH by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state. 
    more » « less