Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

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↗
DOE OSTI2025

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↗
DOE OSTI2016

The final fates of accreting supermassive stars

The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have a less concentrated structure than a fully convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit ≳10{sup 5} M {sub ⊙} derived for the fully convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with M-dot ≲0.1 M{sub ⊙} yr{sup −1}. With M-dot ≃0.3--1 M{sub ⊙} yr{sup −1}, the star becomes GR unstable during the helium-burning stage at M ≃ 2–3.5 × 10{sup 5} M {sub ⊙}. In an extreme case with 10 M{sub ⊙} yr{sup −1}, the star does not collapse until the mass reaches ≃8.0 × 10{sup 5} M {sub ⊙}, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Global hydromagnetic simulations of a planet embedded in a dead zone: Gap opening, gas accretion, and formation of a protoplanetary jet

We present global hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (1) the accretion flow in the Hill sphere is intrinsically three-dimensional for HD and MHD models. Net inflow toward the planet is dominated by high-latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (2) The opening of a gap magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere, and a CPD whose structure displays high levels of variability. (3) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch, a highly collimated, one-sided jet is observed. (4) The CPD's surface density is ∼30 g cm{sup −2}, small enough for significant ionization and turbulence to develop. (5) The accretion rate onto the planet in the MHD simulation reaches a steady value 8 × 10{sup –3} M {sub ⊕} yr{sup –1} and is similar in the viscous HD runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with a dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Variable accretion processes in the young binary-star system UY Aur

We present new K-band spectroscopy of the UY Aur binary star system. Our data are the first to show H{sub 2} emission in the spectrum of UY Aur A and the first to spectrally resolve the Brγ line in the spectrum of UY Aur B. We see an increase in the strength of the Brγ line in UY Aur A and a decrease in Brγ and H{sub 2} line luminosity for UY Aur B compared to previous studies. Converting Brγ line luminosity to accretion rate, we infer that the accretion rate onto UY Aur A has increased by 2 × 10{sup –9} M {sub ☉} yr{sup –1} per year since a rate of zero was observed in 1994. The Brγ line strength for UY Aur B has decreased by a factor of 0.54 since 1994, but the K-band flux has increased by 0.9 mag since 1998. The veiling of UY Aur B has also increased significantly. These data evince a much more luminous disk around UY Aur B. If the lower Brγ luminosity observed in the spectrum of UY Aur B indicates an intrinsically smaller accretion rate onto the star, then UY Aur A now accretes at a higher rate than UY Aur B. However, extinction at small radii or mass pile-up in the circumstellar disk could explain decreased Brγ emission around UY Aur B even when the disk luminosity implies an increased accretion rate. In addition to our scientific results for the UY Aur system, we discuss a dedicated pipeline we have developed for the reduction of echelle-mode data from the ARIES spectrograph.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Anomalous accretion activity and the spotted nature of the DQ Tau binary system

We report the detection of an anomalous accretion flare in the tight eccentric pre-main-sequence binary system DQ Tau. In a multi-epoch survey consisting of randomly acquired low- to moderate-resolution near-infrared spectra obtained over a period of almost 10 yr, we detect a significant and simultaneous brightening of four standard accretion indicators (Ca II infrared triplet, the Paschen and Brackett series H I lines, and He I 1.083 μm), on back-to-back nights (φ = 0.372 and 0.433) with the flare increasing in strength as the system approached apastron (φ = 0.5). The mass accretion rate measured for the anomalous flare is nearly an order of magnitude stronger than the average quiescent rate. While previous observations established that frequent, periodic accretion flares phased with periastron passages occur in this system, these data provide evidence that orbitally modulated accretion flares occur near apastron, when the stars make their closest approach to the circumbinary disk. The timing of the flare suggests that this outburst is due to interactions of the stellar cores (or the highly truncated circumstellar disks) with material in non-axisymmetric structures located at the inner edge of the circumbinary disk. We also explore the optical/infrared spectral type mismatch previously observed for T Tauri stars (TTSs) and successfully model the shape of the spectra from 0.8 to 1.0 μm and the strengths of the TiO and FeH bands as manifestations of large cool spots on the surfaces of the stellar companions in DQ Tau. These findings illustrate that a complete model of near-infrared spectra of many TTSs must include parameters for spot filling factors and temperatures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Prestellar core formation, evolution, and accretion from gravitational fragmentation in turbulent converging flows

We investigate prestellar core formation and accretion based on three-dimensional hydrodynamic simulations. Our simulations represent local ∼1 pc regions within giant molecular clouds where a supersonic turbulent flow converges, triggering star formation in the post-shock layer. We include turbulence and self-gravity, applying sink particle techniques, and explore a range of inflow Mach numbers M=2−16. Two sets of cores are identified and compared: t{sub 1} cores are identified from a time snapshot in each simulation and represent dense structures in a single cloud map; t{sub coll} cores are identified at their individual time of collapse and represent the initial mass reservoir for accretion. We find that cores and filaments form and evolve at the same time. At the stage of core collapse, there is a well-defined, converged characteristic mass for isothermal fragmentation that is comparable to the critical Bonnor–Ebert mass at the post-shock pressure. The core mass functions (CMFs) of t{sub coll} cores show a deficit of high-mass cores (≳7 M{sub ⊙}) compared to the observed stellar initial mass function (IMF). However, the CMFs of t{sub 1} cores are similar to the observed CMFs and include many low-mass cores that are gravitationally stable. The difference between t{sub 1} cores and t{sub coll} cores suggests that the full sample from observed CMFs may not evolve into protostars. Individual sink particles accrete at a roughly constant rate throughout the simulations, gaining one core mass per freefall time even after the initial mass reservoir is accreted. High-mass sinks gain proportionally more mass at later times than low-mass sinks. There are outbursts in accretion rates, resulting from clumpy density structures falling into the sinks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The effect of core-collapse supernova accretion phase turbulence on neutrino flavor evolution

During the accretion phase of a core-collapse supernovae, large amplitude turbulence is generated by the combination of the standing accretion shock instability and convection driven by neutrino heating. The turbulence directly affects the dynamics of the explosion, but there is also the possibility of an additional, indirect, feedback mechanism due to the effect turbulence can have upon neutrino flavor evolution and thus the neutrino heating. Here in this paper we consider the effect of turbulence during the accretion phase upon neutrino evolution, both numerically and analytically. Adopting representative supernova profiles taken from the accretion phase of a supernova simulation, we find the numerical calculations exhibit no effect from turbulence. We explain this absence using two analytic descriptions: the stimulated transition model and the distorted phase effect model. In the stimulated transition model turbulence effects depend upon six different lengthscales, and three criteria must be satisfied between them if one is to observe a change in the flavor evolution due to stimulated transition. We further demonstrate that the distorted phase effect depends upon the presence of multiple semi-adiabatic MSW resonances or discontinuities that also can be expressed as a relationship between three of the same lengthscales. When we examine the supernova profiles used in the numerical calculations we find the three stimulated transition criteria cannot be satisfied, independent of the form of the turbulence power spectrum, and that the same supernova profiles lack the multiple semi-adiabatic MSW resonances or discontinuities necessary to produce a distorted phase effect. Thus, we conclude that even though large amplitude turbulence is present in supernova during the accretion phase, it has no effect upon neutrino flavor evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

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↗
DOE OSTI2021

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↗
DOE OSTI2021

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↗
DOE OSTI2023

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↗
DOE OSTI2024

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↗
DOE OSTI2019

Nucleosynthetic yields from neutron stars accreting in binary common envelopes

Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion on to the neutron star. Accretion on to these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron-rich isotopes of low Z material at the most extreme conditions, and proton-rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modelling of the accretion phase. Finally, we show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

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↗
DOE OSTI2022

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↗
DOE OSTI2025

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 (↗
DOE OSTI2025

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↗