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 zeroangularmomentum 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,
Collapsing stars constitute the main black hole (BH) formation channel, and are occasionally associated with the launch of relativistic jets that power
 NSFPAR ID:
 10524410
 Publisher / Repository:
 IOP Publishing
 Date Published:
 Journal Name:
 The Astrophysical Journal Letters
 Volume:
 952
 Issue:
 2
 ISSN:
 20418205
 Page Range / eLocation ID:
 L32
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

Abstract R _{B}= 10^{3}R _{g}. The BH develops dynamically important largescale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideR _{B}and suppress BH accretion to 1.5% of the Bondi rate, . Thus, lowangularmomentum accretion in the MAD state can form largescale 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 ${\stackrel{\u0307}{M}}_{\mathrm{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% ${\stackrel{\u0307}{M}}_{\mathrm{B}}$ . Because the BAD and RAD states tangle up the jets and destroy them well inside ${\stackrel{\u0307}{M}}_{\mathrm{B}}$R _{B}, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies. 
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 multizone 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 nonspinning SMBH from an external medium with a Bondi radius of
R _{B}≈ 2 × 10^{5}GM _{•}/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 ^{−1}rather thanr ^{−3/2}and the accretion rate is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\stackrel{\u0307}{M}$ . We find continuous energy feedback from the accretion flow to the external medium at a level of ${\stackrel{\u0307}{M}}_{\mathrm{B}}$ . 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. $\sim {10}^{2}\stackrel{\u0307}{M}{c}^{2}\sim 5\phantom{\rule{0.15em}{0ex}}\times \phantom{\rule{0.15em}{0ex}}{10}^{5}{\stackrel{\u0307}{M}}_{\mathrm{B}}{c}^{2}$ 
Abstract We investigate the stellar mass–black hole mass (
) relation with type 1 active galactic nuclei (AGNs) down to ${\mathit{\ue239}}_{*}\u2013{\mathit{\ue239}}_{\mathrm{BH}}$ , corresponding to a ≃ −21 absolute magnitude in restframe ultraviolet, at ${\mathit{\ue239}}_{\mathrm{BH}}={10}^{7\phantom{\rule{0.11em}{0ex}}}{M}_{\odot}$z = 2–2.5. Exploiting the deep and largearea spectroscopic survey of the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX), we identify 66 type 1 AGNs with ranging from 10^{7}–10^{10} ${\mathit{\ue239}}_{\mathrm{BH}}$M _{⊙}that are measured with singleepoch virial method using Civ emission lines detected in the HETDEX spectra. of the host galaxies are estimated from optical to nearinfrared photometric data taken with Spitzer, the Widefield Infrared Survey Explorer, and groundbased 4–8 m class telescopes by ${\mathit{\ue239}}_{*}$CIGALE spectral energy distribution (SED) fitting. We further assess the validity of SED fitting in two cases by hostnuclear decomposition performed through surface brightness profile fitting on spatially resolved host galaxies with the James Webb Space Telescope/NIRCam CEERS data. We obtain the relation covering the unexplored lowmass ranges of ${\mathit{\ue239}}_{*}\u2013{\mathit{\ue239}}_{\mathrm{BH}}$ , and conduct forward modeling to fully account for the selection biases and observational uncertainties. The intrinsic ${\mathit{\ue239}}_{\mathrm{BH}}\phantom{\rule{0.25em}{0ex}}\sim \phantom{\rule{0.25em}{0ex}}{10}^{7}\u2013{10}^{8}\phantom{\rule{0.25em}{0ex}}{M}_{\odot}$ relation at ${\mathit{\ue239}}_{*}\u2013{\mathit{\ue239}}_{\mathrm{BH}}$z ∼ 2 has a moderate positive offset of 0.52 ± 0.14 dex from the local relation, suggestive of more efficient black hole growth at higher redshift even in the lowmass regime of . Our ${\mathit{\ue239}}_{\mathrm{BH}}\phantom{\rule{0.25em}{0ex}}\sim \phantom{\rule{0.25em}{0ex}}{10}^{7}\u2013{10}^{8}\phantom{\rule{0.25em}{0ex}}{M}_{\odot}$ relation is inconsistent with the ${\mathit{\ue239}}_{*}\u2013{\mathit{\ue239}}_{\mathrm{BH}}$ suppression at the low ${\mathit{\ue239}}_{\mathrm{BH}}$ regime predicted by recent hydrodynamic simulations at a 98% confidence level, suggesting that feedback in the lowmass systems may be weaker than those produced in hydrodynamic simulations. ${\mathit{\ue239}}_{*}$ 
Abstract The polluted white dwarf (WD) system SDSS J122859.93+104032.9 (SDSS J1228) shows variable emission features interpreted as originating from a solid core fragment held together against tidal forces by its own internal strength, orbiting within its surrounding debris disk. Estimating the size of this orbiting solid body requires modeling the accretion rate of the polluting material that is observed mixing into the WD surface. That material is supplied via sublimation from the surface of the orbiting solid body. The sublimation rate can be estimated as a simple function of the surface area of the solid body and the incident flux from the nearby hot WD. On the other hand, estimating the accretion rate requires detailed modeling of the surface structure and mixing in the accreting WD. In this work, we present MESA WD models for SDSS J1228 that account for the thermohaline instability and mixing in addition to heavy element sedimentation to constrain accurately the sublimation and accretion rate necessary to supply the observed pollution. We derive a total accretion rate of
, several orders of magnitude higher than the ${\stackrel{\u0307}{M}}_{\mathrm{acc}}=1.8\times {10}^{11}\phantom{\rule{0.25em}{0ex}}\mathrm{g}\phantom{\rule{0.25em}{0ex}}{\mathrm{s}}^{1}$ estimate obtained in earlier efforts. The larger mass accretion rate implies that the minimum estimated radius of the orbiting solid body is ${\stackrel{\u0307}{M}}_{\mathrm{acc}}=5.6\times {10}^{8}\phantom{\rule{0.25em}{0ex}}\mathrm{g}\phantom{\rule{0.25em}{0ex}}{\mathrm{s}}^{1}$ = 72 km, which, although significantly larger than prior estimates, still lies within the upper bounds (a few hundred kilometers) for which the internal strength could no longer withstand the tidal forces from the gravity of the WD. ${r}_{\mathrm{min}}$ 
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 at
z > 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 volumeaveraged neutral fraction of the IGM is either relatively low ( ) or close to unity ( ${\overline{x}}_{\mathrm{H}\phantom{\rule{0.25em}{0ex}}\mathrm{I}}\lesssim {10}^{3}$ ). In particular, the neutral fraction evolution of the IGM at the critical redshift range of ${\overline{x}}_{\mathrm{H}\phantom{\rule{0.25em}{0ex}}\mathrm{I}}\sim 1$z = 6–7 is poorly constrained. We present new constraints on at ${\overline{x}}_{\mathrm{H}\phantom{\rule{0.25em}{0ex}}\mathrm{I}}$z ∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z < 7.09. We derive modelindependent 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 modelindependent constraints on the IGM neutral hydrogen fraction atz ∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of (1 ${\overline{x}}_{\mathrm{H}\phantom{\rule{0.25em}{0ex}}\mathrm{I}}(z=6.3)<0.79\pm 0.04$σ ), (1 ${\overline{x}}_{\mathrm{H}\phantom{\rule{0.25em}{0ex}}\mathrm{I}}(z=6.5)<0.87\pm 0.03$σ ), and (1 ${\overline{x}}_{\mathrm{H}\phantom{\rule{0.25em}{0ex}}\mathrm{I}}(z=6.7)<{0.94}_{0.09}^{+0.06}$σ ). The dark pixel fractions atz > 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.