skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on January 1, 2026

Title: Zooming In On The Multi-Phase Structure of Magnetically-Dominated Quasar Disks: Radiation From Torus to ISCO Across Accretion Rates
Recent radiation-thermochemical-magnetohydrodynamic simulations resolved formation of quasar accretion disks from cosmological scales down to ~300 gravitational radii R g , arguing they were ‘hyper-magnetized’ (plasma β 1 supported by toroidal magnetic fields) and distinct from traditional α -disks. We extend these, refining to 3 R g around a BH with multi-channel radiation and thermochemistry, and exploring a factor of 1000 range of accretion rates ( m ̇ 0.01 20 ). At smaller scales, we see the disks maintain steady accretion, thermalize and self-ionize, and radiation pressure grows in importance, but large deviations from local thermodynamic equilibrium and single-phase equations of state are always present. Trans-Alfvenic and highly-supersonic turbulence persists in all cases, and leads to efficient vertical mixing, so radiation pressure saturates at levels comparable to fluctuating magnetic and turbulent pressures even for m ̇ 1 . The disks also become radiatively inefficient in the inner regions at high m ̇ . The midplane magnetic field remains primarily toroidal at large radii, but at super-Eddington m ̇ we see occasional transitions to a poloidal-field dominated state associated with outflows and flares. Large-scale magnetocentrifugal and continuum radiation-pressure-driven outflows are weak at m ̇ < 1 , but can be strong at m ̇ 1 . In all cases there is a scattering photosphere above the disk extending to 1000 R g at large m ̇ , and the disk is thick and flared owing to magnetic support (with H / R nearly independent of m ̇ ), so the outer disk is strongly illuminated by the inner disk and most of the inner disk continuum scatters or is reprocessed at larger scales, giving apparent emission region sizes as large as .  more » « less
Award ID(s):
2307327 2108230
PAR ID:
10625000
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
OJAp
Date Published:
Journal Name:
The Open Journal of Astrophysics
Volume:
8
ISSN:
2565-6120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 m ̇ c to 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, m ̇ c reaches its maximum, m ̇ c , max 0.02 ( α / 0.3 ) 3 . If atr≫rbthe thin disk accretes with m ̇ d < m ̇ 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
  2. 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, M ̇ 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 M ̇ 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% M ̇ 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
  3. Abstract The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds,that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration (1.2 × 105rg/c) simulation of a rapidly rotating (a= 0.9) black hole feeding from a thick (H/r∼ 0.3), adiabatic, magnetically arrested disk (MAD), whose dynamically important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (i.e., winds and jets, respectively). By carefully disentangling the various angular momentum transport processes within the system, we demonstrate the novel result that disk winds and disk turbulence both extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates of L ̇ r 0.9 and M ̇ r 0.4 , respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma-ray bursts. 
    more » « less
  4. Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process (r-process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We presentνh-amr, a novel neutrino-transport general relativistic magnetohydrodynamic (νGRMHD) code, which we use to perform the first 3D globalνGRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ 0.1 1 M s - 1 , do not produce significant neutron-rich (Ye < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ few M s - 1 (e.g., for more centrally concentrated progenitors), can unbindMej ≲ Mof neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflowYe; however, the finalr-process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines becomer-process factories. 
    more » « less
  5. 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 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 M ̇ is consequently suppressed by over 2 orders of magnitude below the Bondi rate M ̇ B . We find continuous energy feedback from the accretion flow to the external medium at a level of 10 2 M ̇ c 2 5 × 10 5 M ̇ B c 2 . 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