skip to main content

This content will become publicly available on December 1, 2024

Title: Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

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
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Date Published:
Journal Name:
The Astrophysical Journal Letters
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDEs). Motivated by recent hydrodynamical simulations, we assume that the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radiusRv∼ 1014cm (photosphere radius ∼1015cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contractionRvt−1, its temperature risingTefft1/2while its total luminosity remains roughly constant; the optical luminosity decays asνLνRv2Tefft3/2. Despite this similarity to the mass fallback rateṀfbt5/3, envelope heating from fallback accretion is subdominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days may account for the delayed onset of thermal X-rays, late-time radio flares, and high-energy neutrino generation, observed in some TDEs. We compare the model predictions to recent TDE light-curve correlation studies, finding both agreement and points of tension.

    more » « less
  2. 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
  3. 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)are9.35.4+4.6and4.22.0+1.9Mpc2(Kkms1)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(21)Σ0.5andαCO(10)Σ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
  4. 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
  5. Abstract

    We use ALMA observations of CO(2–1) in 13 massive (M*≳ 1011M) poststarburst galaxies atz∼ 0.6 to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The poststarburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (Data Release 14) based on their spectral shapes, as part of the Studying QUenching at Intermediate-z Galaxies: Gas, anguLarmomentum, and Evolution (SQuIGGLE) program. Early results showed that two poststarburst galaxies host large H2reservoirs despite their low inferred star formation rates (SFRs). Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, withMH2109M. Given their high stellar masses, this mass limit corresponds to an average gas fraction offH2MH2/M*7%or ∼14% using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with theDn4000 spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched ≲150 Myr prior to observation host detectable CO(2–1) emission, while older poststarburst galaxies are undetected. The large H2reservoirs and low SFRs in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100–200 Myr, theSQuIGGLEgalaxies require the additional and efficient heating or removal of cold gas to bring their low SFRs in line with standard H2scaling relations.

    more » « less