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 of
We present an exact
- Award ID(s):
- 2110425
- NSF-PAR ID:
- 10472067
- Publisher / Repository:
- IOP Science - Institute of Physics
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 945
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 76
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract R B≈ 2 × 105GM •/c 2, whereM •is 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 . We find continuous energy feedback from the accretion flow to the external medium at a level of . 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. -
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,
R B= 103R g. The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideR Band suppress BH accretion to 1.5% of the Bondi rate, . 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 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% . Because the BAD and RAD states tangle up the jets and destroy them well insideR B, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies. -
Abstract We present theoretical expectations for infall toward supercluster-scale cosmological filaments, motivated by the Arecibo Pisces–Perseus Supercluster Survey (APPSS) to map the velocity field around the Pisces–Perseus Supercluster (PPS) filament. We use a minimum spanning tree applied to dark matter halos the size of galaxy clusters to identify 236 large filaments within the Millennium simulation. Stacking the filaments along their principal axes, we determine a well-defined, sharp-peaked velocity profile function that can be expressed in terms of the maximum infall rate
V maxand the distanceρ maxbetween the location of maximum infall and the principal axis of the filament. This simple, two-parameter functional form is surprisingly universal across a wide range of linear mass densities.V maxis positively correlated with the halo mass per length along the filament, andρ maxis negatively correlated with the degree to which the halos are concentrated along the principal axis. We also assess an alternative, single-parameter method usingV 25, the infall rate at a distance of 25 Mpc from the axis of the filament. Filaments similar to the PPS have 116 km s−1, Mpc, andV 25= 329 ± 68 km s−1. We create mock observations to model uncertainties associated with viewing angle, lack of three-dimensional velocity information, limited sample size, and distance uncertainties. Our results suggest that it would be especially useful to measure infall for a larger sample of filaments to test our predictions for the shape of the infall profile and the relationships among infall rates and filament properties. -
Abstract We confirm the planetary nature of TOI-5344 b as a transiting giant exoplanet around an M0-dwarf star. TOI-5344 b was discovered with the Transiting Exoplanet Survey Satellite photometry and confirmed with ground-based photometry (the Red Buttes Observatory 0.6 m telescope), radial velocity (the Habitable-zone Planet Finder), and speckle imaging (the NN-Explore Exoplanet Stellar Speckle Imager). TOI-5344 b is a Saturn-like giant planet (
ρ = 0.80 g cm−3) with a planetary radius of 9.7 ± 0.5R ⊕(0.87 ± 0.04R Jup) and a planetary mass of (0.42 ). It has an orbital period of days and an orbital eccentricity of . We measure a high metallicity for TOI-5344 of [Fe/H] = 0.48 ± 0.12, where the high metallicity is consistent with expectations from formation through core accretion. We compare the metallicity of the M-dwarf hosts of giant exoplanets to that of M-dwarf hosts of nongiants (≲8R ⊕). While the two populations appear to show different metallicity distributions, quantitative tests are prohibited by various sample caveats. -
Abstract We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star,
r p, is well within the tidal radius of the black hole,r t, i.e., whenβ ≡r t/r p≫ 1. We find that shocks form forβ ≳ 3, but they are weak (with Mach numbers ∼1) for allβ , and that they reach the center of the star prior to the time of maximum adiabatic compression forβ ≳ 10. The maximum density and temperature reached during the TDE follow much shallower relations withβ than the previously predicted and scalings. Belowβ ≃ 10, this shallower dependence occurs because the pressure gradient is dynamically significant before the pressure is comparable to the ram pressure of the free-falling gas, while aboveβ ≃ 10, we find that shocks prematurely halt the compression and yield the scalings and . We find excellent agreement between our results and high-resolution simulations. Our results demonstrate that, in the Newtonian limit, the compression experienced by the star is completely independent of the mass of the black hole. We discuss our results in the context of existing (affine) models, polytropic versus non-polytropic stars, and general relativistic effects, which become important when the pericenter of the star nears the direct capture radius, atβ ∼ 12.5 (2.7) for a solar-like star disrupted by a 106M ⊙(107M ⊙) supermassive black hole.