Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accretion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Rapid accretion state transitions following the tidal disruption event AT2018fyk

Following a tidal disruption event (TDE), the accretion rate can evolve from quiescent to near- Eddington levels and back over months{years timescales. This provides a unique opportunity to study the formation and evolution of the accretion flow around supermassive black holes (SMBHs). We present two years of multi-wavelength monitoring observations of the TDE AT2018fyk at X-ray, UV, optical and radio wavelengths. We identify three distinct accretion states and two state transitions between them. These appear remarkably similar to the behaviour of stellar-mass black holes in outburst. The X-ray spectral properties show a transition from a soft (thermal-dominated) to a hard (power-law dominated) spectral state around Lbol ~few x 10^(-2) LEdd, and the strengthening of the corona over time ~100-200 days after the UV/optical peak. Contemporaneously, the spectral energy distribution (in particular, the UV-to-X-ray spectral slope alphaox ) shows a pronounced softening as the outburst progresses. The X-ray timing properties also show a marked change, initially dominated by variability at long (>day) timescales while a high frequency (~10-3 Hz) component emerges after the transition into the hard state. At late times (~500 days after peak), a second accretion state transition occurs, from the hard into the quiescent state, as identi fied by the sudden collapse of the bolometric (X-ray+UV) emission to levels below 10-3.4 LEdd. Our fi ndings illustrate that TDEs can be used to study the scale (in)variance of accretion processes in individual SMBHs. Consequently, they provide a new avenue to study accretion states over seven orders of magnitude in black hole mass, removing limitations inherent to commonly used ensemble studies.

T. Wevers↗

Electromagnetic Time and Wavelength Domain Predictions for Accreting Binary Black Holes

Accreting massive binary black holes in gaseous environments are keymulti-messenger sources for the LISA mission. We will present some ofthe latest predictions of the observational signatures of accreting binary black holes across the electromagnetic spectrum from the fieldand provide some of our own. In particular, we will provide a briefsummary of the progress made in our group to understand these systemstheoretically using high-performance simulations of generalrelativistic magnetohydrodynamics (MHD) and dynamic spacetimes.Simulating these systems realistically is challenging asradiation-coupled MHD must be considered over large dynamic ranges inspace and time. Constraining our view to the circumbinary disk regionfirst, we will show how the evolution duration, accretion disk size,mass ratio, and thermodynamics model affect the structure andvariability of the accretion flow. We particularly emphasize how theseparameters influence the over-density feature, or "lump", which orbitsthe binary near the edge of the cavity, since it is responsible formost of the electromagnetic emission's variability---a key signatureof a system being a binary. Extending to smaller length scales, wewill report on simulations following accretion all the down to theevent horizons so that we may begin to investigate how black hole spinaffects mini-disk dynamics, accretion rate, and jet power. We willalso report on our radiative transfer predictions for how black holespin alters the spectral energy distribution of these systems.Wrapping up, we will itemize the remaining key theoreticaltechnologies and directions needed to be accomplished to betterunderstand these sources, and assess strategies for searching for themin the coming decade.

Scott C Noble↗

Bremsstrahlung in GRMHD Models of Accreting Black Holes

Here, the role of bremsstrahlung in the emission from hot accretion flows around slowly accreting supermassive black holes is not thoroughly understood. In order to appraise the importance of bremsstrahlung relative to other radiative processes, we compute spectral energy distributions (SEDs) of accretion disks around slowly accreting supermassive black holes including synchrotron radiation, inverse Compton scattering, and bremsstrahlung. We compute SEDs for (i) four axisymmetric radiative general relativistic magnetohydrodynamics (RadGRMHD) simulations of 10 8 M ⊙ black holes with accretion rates between 10 -8 $\dot{M}$ Edd and 10 -5 $\dot{M}$ Edd , (ii) four axisymmetric RadGRMHD simulations of M87* with varying dimensionless spin a$_*$ and black hole mass, and (iii) a 3D GRMHD simulation scaled for Sgr A*. At 10 -8 $\dot{M}$ Edd , most of the luminosity is synchrotron radiation, while at 10 -5 $\dot{M}$ Edd the three radiative processes have similar luminosities. In most models, bremsstrahlung dominates the SED near 512 keV. In the M87* models, bremsstrahlung dominates this part of the SED if a$_*$ = 0.5, but inverse Compton scattering dominates if a$_*$ = 0.9375. Since scattering is more variable than bremsstrahlung, this result suggests that 512 keV variability could be a diagnostic of black hole spin. In the appendix, we compare some bremsstrahlung formulæ found in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a "multizone" technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH ($a\ast$ = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r –1 inside the Bondi radius R B , which is located at R B = 2 × 10 5 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ~ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is ≈1% of the analytical Bondi accretion rate ${\dot{M}}_{{\rm{B}}};$ and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate $\approx 0.02\dot{M}{c}^{2}$. Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ~ 10 2 –10 7 r g , and find that $\dot{M}/{\dot{M}}_{{B}}\approx {({R}_{{B}}/6\,{r}_{g})}^{-0.5}$.

79 ASTRONOMY AND ASTROPHYSICS↗

River bank erosion and accretion rates, and planform metrics of the East River, downstream of Gothic Colorado over intervals between 1955 and 2015

This dataset provides the tabular summary of analysis of an alluvial floodplain reach of the East River, downstream of Gothic, CO near Crested Butte. The measurements include bank erosion and accretion rates, channel width, bank curvature, and the aspect/orientation of the river banks. The measurements were derived from binary masks of the location of the river channel from aerial and satellite imagery collected on seven dates over a sixty-year timespan: 1955, 1973, 1983, 1990, 2001, 2011, and 2015. The masks were analyzed using the Spatially Continuous Riverbank Erosion and Accretion Measurements (SCREAM) software detailed in Rowland et al. 2016. The masks used in this analysis can be found an accompanying dataset (DOI: ). Rates of change along the East River were measured over a total of 8-time intervals: 1955-1973; 1955-2015; 1973-1983; 1983-1990; 1990-2001; 2001-2011; 2001-2015; and 2011-2015.In files with “summary” in the name, the data is provided at a pixel level, where each mapped bank pixel has an associated erosion or accretion value, a channel width, a curvature value, and an aspect each river and time period will have an individual file. Files with “Segments” in name provide data that is averaged along segments of the rivers. These segments are approximately 10 channel widths in length. In addition to erosion and accretion rates, the segment-based results include area measurements of erosion and accretion, islands, and channels. The number of islands is also included.

54 ENVIRONMENTAL SCIENCES↗

Accretion Properties of PDS 70b with MUSE

We report a new evaluation of the accretion properties of PDS 70b obtained with the Very Large Telescope/Multi Unit Spectroscopic Explorer. The main difference from the previous studies of Haffert et al. and Aoyama & Ikoma is in the mass accretion rate. Simultaneous multiple line observations, such as Hα and Hβ, can better constrain the physical properties of an accreting planet. While we clearly detected Hα emissions from PDS 70b, no Hβ emissions were detected. We estimate the line flux of Hβ with a 3σ upper limit to be 2.3 × 10{sup −16} erg s{sup −1} cm{sup −2}. The flux ratio F {sub Hβ}/F {sub Hα} for PDS 70b is <0.28. Numerical investigations by Aoyama et al. suggest that F {sub Hβ}/F {sub Hα} should be close to unity if the extinction is negligible. We attribute the reduction of the flux ratio to the extinction, and estimate the extinction of Hα (A {sub Hα}) for PDS 70b to be >2.0 mag using the interstellar extinction value. By combining with the Hα linewidth and the dereddening line luminosity of Hα, we derive the PDS 70b mass accretion rate to be ≳5 × 10{sup −7} M {sub Jup} yr{sup −1}. The PDS 70b mass accretion rate is an order of magnitude larger than that of PDS 70. We found that the filling factor f {sub f} (the fractional area of the planetary surface emitting Hα) is ≳0.01, which is similar to the typical stellar value. The small value of f {sub f} indicates that the Hα emitting areas are localized at the surface of PDS 70b.

79 ASTRONOMY AND ASTROPHYSICS↗

Keck/NIRC2 L’-band Imaging of Jovian-mass Accreting Protoplanets around PDS 70

We present L’-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wave front sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semimajor axis of 20{sub −4}{sup +3} au and PDS 70 c to have a semimajor axis of 34{sub −6}{sup +12} au (95% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2–3 R {sub Jup}. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 M {sub Jup} and a mean mass accretion rate between 3 × 10{sup −7} and 8 × 10{sup −7} M {sub Jup}/yr. For PDS 70 c, we computed a mass between 1 and 3 M {sub Jup} and mean mass accretion rate between 1 × 10{sup −7} and 5 × 10{sup −7} M {sub Jup}/yr. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets’ SEDs.

79 ASTRONOMY AND ASTROPHYSICS↗

Spin and Accretion Rate Dependence of Black Hole X-Ray Spectra

We present a survey of how the spectral features of black hole X-ray binary systems depend on spin, accretion rate, viewing angle, and Fe abundance when predicted on the basis of first-principles physical calculations. The power-law component hardens with increasing spin. The thermal component strengthens with increasing accretion rate. The Compton bump is enhanced by higher accretion rate and lower spin. The Fe Kα equivalent width grows sublinearly with Fe abundance. Strikingly, the Kα profile is more sensitive to accretion rate than to spin because its radial surface brightness profile is relatively flat, and higher accretion rate extends the production region to smaller radii. The overall radiative efficiency is at least 30%–100% greater than as predicted by the Novikov–Thorne model.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiation GRMHD Models of Accretion onto Stellar-mass Black Holes. I. Survey of Eddington Ratios

We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

The chemical compositions of accreted and in situ galactic globular clusters according to SDSS/APOGEE

ABSTRACT Studies of the kinematics and chemical compositions of Galactic globular clusters (GCs) enable the reconstruction of the history of star formation, chemical evolution, and mass assembly of the Galaxy. Using the latest data release (DR16) of the SDSS/APOGEE survey, we identify 3090 stars associated with 46 GCs. Using a previously defined kinematic association, we break the sample down into eight separate groups and examine how the kinematics-based classification maps into chemical composition space, considering only α (mostly Si and Mg) elements and Fe. Our results show that (i) the loci of both in situ and accreted subgroups in chemical space match those of their field counterparts; (ii) GCs from different individual accreted subgroups occupy the same locus in chemical space. This could either mean that they share a similar origin or that they are associated with distinct satellites which underwent similar chemical enrichment histories; (iii) the chemical compositions of the GCs associated with the low orbital energy subgroup defined by Massari and collaborators is broadly consistent with an in situ origin. However, at the low-metallicity end, the distinction between accreted and in situ populations is blurred; (iv) regarding the status of GCs whose origin is ambiguous, we conclude the following: the position in Si–Fe plane suggests an in situ origin for Liller 1 and a likely accreted origin for NGC 5904 and NGC 6388. The case of NGC 288 is unclear, as its orbital properties suggest an accretion origin, its chemical composition suggests it may have formed in situ.

Horta, Danny↗

Observing correlations between dark matter accretion and galaxy growth: II. testing the impact of galaxy mass, star formation indicator, and neighbour colours

A crucial question in galaxy formation is what role new accretion has in star formation. Theoretical models have predicted a wide range of correlation strengths between halo accretion and galaxy star formation. Previously, we presented a technique to observationally constrain this correlation strength for isolated Milky Way mass galaxies at z ~ 0.12, based on the correlation between halo accretion and the density profile of neighbouring galaxies. By applying this technique to both observational data from the Sloan Digital Sky Survey and simulation data from the UNIVERSEMACHINE , where we can test different correlation strengths, we ruled out positive correlations between dark matter accretion and recent star formation activity. In this work, we expand our analysis by (1) applying our technique separately to red and blue neighbouring galaxies, which trace different infall populations, (2) correlating dark matter accretion rates with D n 4000 measurements as a longer-term quiescence indicator than instantaneous star-formation rates, and (3) analysing higher-mass isolated central galaxies with 10 11.0 < M * /M ⊙ < 10 11.5 out to z ~ 0.18. Finally, in all cases, our results are consistent with non-positive correlation strengths with ≳ 85 percent confidence, which is most consistent with models where processes such as gas recycling dominate star formation in massive z = 0 galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Using action space clustering to constrain the recent accretion history of Milky Way-like galaxies

ABSTRACT In the currently favoured cosmological paradigm galaxies form hierarchically through the accretion of satellites. Since a satellite is less massive than the host, its stars occupy a smaller volume in action space. Actions are conserved when the potential of the host halo changes adiabatically, so stars from an accreted satellite would remain clustered in action space as the host evolves. In this paper, we identify recently disrupted accreted satellites in three Milky Way-like disc galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellites through simulation snapshots. We try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the z = 0 snapshot. Even with completely error-free mock data we find that only 35 per cent (14/39) satellites are well recovered while the rest (25/39) are poorly recovered (i.e. either contaminated or split up). Most (10/14 ∼70 per cent) of the well-recovered satellites have infall times <7.1 Gyr ago and total mass >4 × 108M⊙ (stellar mass more than 1.2 × 106 M⊙, although our upper mass limit is likely to be resolution dependent). Since cosmological simulations predict that stellar haloes include a population of in situ stars, we test our ability to recover satellites when the data include 10–50 per cent in situ contamination. We find that most previously well-recovered satellites stay well recovered even with 50 per cent contamination. With the wealth of 6D phase space data becoming available we expect that cluster analysis in action space will be useful in identifying the majority of recently accreted and moderately massive satellites in the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

Testing AGN outflow and accretion models with C IV and He II emission line demographics in z ≈ 2 quasars

Using ≈190 000 spectra from the 17th data release of the Sloan Digital Sky Survey (SDSS), we investigate the ultraviolet emission line properties in z ≈ 2 quasars. Specifically, we quantify how the shape of C IV λ1549 and the equivalent width (EW) of He II λ1640 depend on the black hole mass and Eddington ratio inferred from Mg II λ2800. Above L/L Edd ≳ 0.2, there is a strong mass dependence in both C IV blueshift and He II EW. Large C IV blueshifts are observed only in regions with both high mass and high accretion rate. Including X-ray measurements for a subsample of 5000 objects, we interpret our observations in the context of AGN accretion and outflow mechanisms. The observed trends in He II and 2 keV strength are broadly consistent with theoretical qsosed models of AGN spectral energy distributions (SEDs) for low spin black holes, where the ionizing SED depends on the accretion disc temperature and the strength of the soft excess. High spin models are not consistent with observations, suggesting SDSS quasars at z ≈ 2 may in general have low spins. We find a dramatic switch in behaviour at L/L Edd ≲ 0.1: the ultraviolet emission properties show much weaker trends, and no longer agree with qsosed predictions, hinting at changes in the structure of the broad line region. Overall, the observed emission line trends are generally consistent with predictions for radiation line driving where quasar outflows are governed by the SED, which itself results from the accretion flow and hence depends on both the SMBH mass and accretion rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Establishing accretion flares from supermassive black holes as a source of high-energy neutrinos

ABSTRACT The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

Astronomy & Astrophysics↗

WASP-107b’s Density Is Even Lower: A Case Study for the Physics of Planetary Gas Envelope Accretion and Orbital Migration

With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet’s low surface gravity and the star’s brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4 yr Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and provide a detailed study of the physics governing the accretion of the gas envelope of WASP-107b. We reveal that WASP-107b’s mass is only 1.8 Neptune masses (M {sub b} = 30.5 ± 1.7 M {sub ⊕}). The resulting extraordinarily low density suggests that WASP-107b has a H/He envelope mass fraction of >85% unless it is substantially inflated. The corresponding core mass of <4.6 M {sub ⊕} at 3σ is significantly lower than what is traditionally assumed to be necessary to trigger massive gas envelope accretion. We demonstrate that this large gas-to-core mass ratio most plausibly results from the onset of accretion at ≳1 au onto a low-opacity, dust-free atmosphere and subsequent migration to the present-day a {sub b} = 0.0566 ± 0.0017 au. Beyond WASP-107b, we also detect a second, more massive planet (M{sub c}sin i=0.36±0.04M{sub J}) on a wide eccentric orbit (e {sub c} = 0.28 ± 0.07) that may have influenced the orbital migration and spin–orbit misalignment of WASP-107b. Overall, our new RV observations and envelope accretion modeling provide crucial insights into the intriguing nature of WASP-107b and the system’s formation history. Looking ahead, WASP-107b will be a keystone planet to understand the physics of gas envelope accretion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MaNGA 8313-1901: Gas Accretion Observed in a Blue Compact Dwarf Galaxy?

Gas accretion is an important process in the evolution of galaxies, but it has limited direct observational evidences. In this paper, we report the detection of a possible ongoing gas accretion event in a blue compact dwarf (BCD) galaxy, MaNGA 8313-1901, observed by the Mapping Nearby Galaxies and Apache Point Observatory (MaNGA) program. This galaxy has a distinct off-centered blue clump to the northeast (the NE clump) that shows low metallicity and enhanced star formation. The kinematics of the gas in the NE clump also seems to be detached from the host BCD galaxy. Together with the metallicity drop of the NE clump, it suggests that the NE clump likely has an external origin, such as gas accretion or galaxy interaction, rather than an internal origin, such as an H II complex in the disk. After removing the underlying host component, we find that the spectrum of the "pure" clump can match very well with a modeled spectrum containing a stellar population of the young stars (≤7 Myr) only. This may imply that the galaxy is experiencing an accretion of cold gas, instead of a merger event involving galaxies with significant preexisting old stars. We also find signs of another clump (the SW clump) at the southwest corner of the host galaxy, and the two clumps may share the same origin of gas accretion.

79 ASTRONOMY AND ASTROPHYSICS↗

Joint Modeling of Quasar Variability and Accretion Disk Reprocessing Using Latent Stochastic Differential Equations

Quasars are bright active galactic nuclei powered by the accretion of matter around supermassive black holes at the center of galaxies. Their stochastic brightness variability depends on the physical properties of the accretion disk and black hole. The upcoming Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to observe tens of millions of quasars, so there is a need for efficient techniques like machine learning that can handle the large volume of data. Quasar variability is believed to be driven by an X-ray corona, which is reprocessed by the accretion disk and emitted as UV/optical variability. We are the first to introduce an auto-differentiable simulation of the accretion disk and reprocessing. We use the simulation as a direct component of our neural network to jointly model the driving variability and reprocessing, trained with supervised learning on simulated LSST-like 10 yr quasar light curves. We encode the light curves using a transformer encoder, and the driving variability is reconstructed using latent stochastic differential equations, a physically motivated generative deep learning method that can model continuous-time stochastic dynamics. By embedding the physical processes of the driving signal and reprocessing into our network, we achieve a model that is more robust and interpretable. We demonstrate that our model outperforms a Gaussian process regression baseline and can infer accretion disk parameters and time delays between wave bands, even for out-of-distribution driving signals. Our approach provides a powerful framework that can be adapted to solve other inverse problems in multivariate time series.

Fagin, Joshua [City Univ. of New York (CUNY), NY (↗

Bridging Scales in Black Hole Accretion and Feedback: Relativistic Jet Linking the Horizon to the Host Galaxy

Simulating black hole (BH) accretion and feedback from the BH horizon to galactic scales is extremely challenging, as it involves a vast range of scales. Recently, our multizone method has successfully achieved global dynamical steady states of hot accretion flows in 3D general relativistic magnetohydrodynamic simulations by tracking the bidirectional interaction between a nonspinning BH and its host galaxy. In this paper, we present technical improvements to the method and apply it to spin a * = 0.9 BHs, which power relativistic jets. We first test the new multizone setup with a smaller Bondi radius, R B ≈ 400 r g , where r g is the gravitational radius. The strongly magnetized accretion launches a relativistic jet with an intermediate feedback efficiency η ∼ 30%, in between that of a prograde (η ∼ 100%) and retrograde (η ∼ 10%) torus. Interestingly, both prograde and retrograde simulations also eventually converge to the same intermediate efficiency when evolved long enough, as accumulated magnetic fields remove gas rotation. We then extend strongly magnetized simulations to larger Bondi radii, R B ≈ 2 × 10 3 , 2 × 10 4 , 2 × 10 5 r g . We find that the BH accretion rate $\dot{M}$ is suppressed with respect to the Bondi rate as $\dot{M}_{\textrm{B}}$ as $\dot{M}/\dot{M}_{\textrm{B}} ∝ R_{\textrm{B}}^{-1/2}$. However, despite some variability, the time-averaged feedback efficiency remains at η ∼ 30%, independent of R B . This suggests that BH feedback efficiency in hot accretion flows is mainly governed by the BH spin (a * ) rather than by the galactic properties (R B ). From these first-principles simulations, we provide a feedback subgrid prescription for cosmological simulations: $\dot{E}_{\textrm{fb}} = 2$ x $10^{-3}[R_{\textrm{B}}/(2$ x $10^5 r_g)]^{-1/2}$ $\dot{M}_{\textrm{B}}c^2$ for BH spin a * = 0.9.

79 ASTRONOMY AND ASTROPHYSICS↗