Title: Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD
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 M ̇ is consequently suppressed by over 2 orders of magnitude below the Bondi rate M ̇ B . We find continuous energy feedback from the accretion flow to the external medium at a level of 10 2 M ̇ c 2 5 × 10 5 M ̇ B c 2 . 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
Award ID(s):
1743747
PAR ID:
10486083
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
ADS
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
959
Issue:
2
ISSN:
2041-8205
Page Range / eLocation ID:
L22
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We measure the CO-to-H2conversion factor (αCO) in 37 galaxies at 2 kpc resolution, using the dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metal ratio. In total, we have ∼790 and ∼610 independent measurements ofαCOfor CO (2–1) and (1–0), respectively. The mean values forαCO (2–1)andαCO (1–0)are 9.3 5.4 + 4.6 and 4.2 2.0 + 1.9 M pc 2 ( K km s 1 ) 1 , respectively. The CO-intensity-weighted mean is 5.69 forαCO (2–1)and 3.33 forαCO (1–0). We examine howαCOscales with several physical quantities, e.g., the star formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ( U ¯ ). Among them, U ¯ , ΣSFR, and the integrated CO intensity (WCO) have the strongest anticorrelation with spatially resolvedαCO. We provide linear regression results toαCOfor all quantities tested. At galaxy-integrated scales, we observe significant correlations betweenαCOandWCO, metallicity, U ¯ , and ΣSFR. We also find thatαCOin each galaxy decreases with the stellar mass surface density (Σ) in high-surface-density regions (Σ≥ 100Mpc−2), following the power-law relations α CO ( 2 1 ) Σ 0.5 and α CO ( 1 0 ) Σ 0.2 . The power-law index is insensitive to the assumed dust-to-metal ratio. We interpret the decrease inαCOwith increasing Σas a result of higher velocity dispersion compared to isolated, self-gravitating clouds due to the additional gravitational force from stellar sources, which leads to the reduction inαCO. The decrease inαCOat high Σis important for accurately assessing molecular gas content and star formation efficiency in the centers of galaxies, which bridge “Milky Way–like” to “starburst-like” conversion factors. 
    more » « less
  2. Abstract The growth of supermassive black holes is strongly linked to their galaxies. It has been shown that the population mean black hole accretion rate ( BHAR ¯ ) primarily correlates with the galaxy stellar mass (M) and redshift for the general galaxy population. This work aims to provide the best measurements of BHAR ¯ as a function ofMand redshift over ranges of 109.5<M< 1012Mandz< 4. We compile an unprecedentedly large sample with 8000 active galactic nuclei (AGNs) and 1.3 million normal galaxies from nine high-quality survey fields following a wedding cake design. We further develop a semiparametric Bayesian method that can reasonably estimate BHAR ¯ and the corresponding uncertainties, even for sparsely populated regions in the parameter space. BHAR ¯ is constrained by X-ray surveys sampling the AGN accretion power and UV-to-infrared multiwavelength surveys sampling the galaxy population. Our results can independently predict the X-ray luminosity function (XLF) from the galaxy stellar mass function (SMF), and the prediction is consistent with the observed XLF. We also try adding external constraints from the observed SMF and XLF. We further measure BHAR ¯ for star-forming and quiescent galaxies and show that star-forming BHAR ¯ is generally larger than or at least comparable to the quiescent BHAR ¯
    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, M ̇ 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 M ̇ B reach 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% M ̇ 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 Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ( x ¯ H I 10 3 ) or close to unity ( x ¯ H I 1 ). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints on x ¯ H I atz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of x ¯ H I ( z = 6.3 ) < 0.79 ± 0.04 (1σ), x ¯ H I ( z = 6.5 ) < 0.87 ± 0.03 (1σ), and x ¯ H I ( z = 6.7 ) < 0.94 0.09 + 0.06 (1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018. 
    more » « less
  5. Abstract We have observed electron impact fluorescence from CO2to excite the Cameron bands (CBs), CO (a3Π →X1Σ+; 180–280 nm), the first-negative group (1NG) bands, CO+(B2Σ+→X2Σ+; 180–320 nm), the fourth-positive group (4PG) bands, CO (A1Π →X1Σ+; 111–280 nm), and the UV doublet, CO2+( B ˜ 2 Σ u + X ˜ 2 Π g ; 288.3 and 289.6 nm) in the ultraviolet (UV). This wavelength range matches the spectral region of past and present spacecraft equipped to observe UV dayglow and aurora emissions from the thermospheres (100–300 km) of Mars and Venus. Our large vacuum system apparatus is able to measure the emission cross sections of the strongest optically forbidden UV transitions found in planetary spectra. Based on our cross-sectional measurements, previous CB emission cross-sectional errors exceed a factor of 3. The UV doublet lifetime is perturbed through B ˜ 2 Σ u + A ˜ 2 Π u spin–orbit coupling. Forward modeling codes of the Mars dayglow have not been accurate in the mid-UV due to systematic errors in these two emission cross sections. We furnish absolute emission cross sections for several band systems over electron energies 20–100 eV for CO2. We present a CB lifetime, which together with emission cross sections, furnish a set of fundamental physical constants for electron transport codes such as AURIC (Atmospheric Ultraviolet Radiance Integrated Code). AURIC and Trans-Mars are used in the analysis of UV spectra from the Martian dayglow and aurora. 
    more » « less