Multi-Phase Thermal Structure & The Origin of the Broad-Line Region, Torus, and Corona in Magnetically-Dominated Accretion Disks

TL;DR

This paper develops an analytic model for hyper-magnetized, flux-frozen accretion disks around supermassive black holes, extending from near-horizon to galactic scales. By enforcing magnetic pressure dominance (β ≪ 1) and trans-Alfvénic turbulence, the authors derive a self-similar disk structure in Keplerian potentials and compute the thermo-chemical states across multiple zones. Remarkably, the model naturally yields the BLR, dusty torus, warm Comptonizing skin, and hard X-ray corona as integral disk components without introducing new physics, instead as emergent properties of the disk’s geometry and opacities. The work highlights strong, observable predictions for covering factors, densities, ionization parameters, and reprocessing fractions, offering a cohesive framework that links accretion-disk physics to the AGN ecosystem and guiding future radiative-transfer studies and observations.

Abstract

Recent simulations have demonstrated the formation of 'flux-frozen' and hyper-magnetized disks, qualitatively distinct from both classical $α$ disks and magnetically-arrested disks, as a natural consequence of fueling gas to supermassive black holes in galactic nuclei. We previously showed that the dynamical structure of said disks can be approximated by simple analytic similarity models. Here we study the thermal properties of these models over a wide range of physical scales and accretion rates (from highly sub-critical to super-critical). We show there are several characteristic zones: a dusty torus-like region, a multi-phase neutral and then multi-phase ionized, broad line-emitting region interior to the sublimation radius, before finally a transition to a thermal accretion disk with a warm Comptonizing layer. The disks are strongly-flared with large scale heights, and reprocess and/or scatter an order-one fraction of the central disk emission. As a result, this simple accretion disk model predicts phenomena including the existence of a dusty torus and its covering factor, geometry, clumpiness, and dust temperatures; a broad-line-region (BLR) with its characteristic sizes and luminosities and ionization properties; extended scattering/reprocessing surfaces producing cooler disk continuum and apparently large observed disk sizes; and existence of warm Comptonizing layers and hard coronal gas. Remarkably, these properties emerge without our having to introduce new components or parameters: they are all part of the accretion flow if the disks are in the hyper-magnetized limit.

Figures (14)

• Figure 1: Heuristic illustration of the accretion disk properties for a magnetically-dominated, flux-frozen disk, with distinct "zones" with different thermo-chemical properties. Black line shows the disk scale-height $z=H$ versus cylindrical radius $R$. We label critical radii including: the gravitational radius/horizon $R_{g}$, ISCO $R_{\rm ISCO}$, radius interior to which the midplane is thermalized and self-ionized $R_{\rm therm,\,ion}$, radius at which the atomic disk is thermalized $R_{\rm therm,\,atomic}$, radius of dust sublimation in the shielded midplane $R_{{\rm sub},\,T_{\rm eff}}$, radius interior to which there is a Thompson-thick electron-scattering intercepting surface illumination $R_{\rm es,\,ion}$, radius interior to which the surface/illuminated layers have dust sublimated $R_{\rm sub,\,ext}$, and BH radius of influence $R_{\rm BHROI}$. We divide the system into distinct "zones" labeled (§ \ref{['sec:zone1.ism']}-\ref{['sec:zone8.warm.scattering']}), described in Fig. \ref{['fig:cartoon.quadrants']}.
• Figure 2: Illustration as Fig. \ref{['fig:cartoon.definitions']} of a flux-frozen disk, with descriptions of each zone. Top right: Zone label (§ \ref{['sec:zone1.ism']}-\ref{['sec:zone8.warm.scattering']}; as Fig. \ref{['fig:cartoon.definitions']}). Black line shows disk scale-height $z=H$, versus cylindrical radius $R$. Bottom right: Phase of the gas in each zone (e.g. "dusty" outside of the sublimation radii, or atomic, or ionized). Bottom left: Whether each phase is effectively optically-thin or thick to absorption, and multiphase or thermalized. Top left: Phenomenological structures corresponding to each zone. Zones include the galactic ISM (1; § \ref{['sec:zone1.ism']}), dust torus (2; § \ref{['sec:zone2.torus']}), broad-line region (3; § \ref{['sec:zone3.blr']}), neutral (4; § \ref{['sec:zone4.neutral']}) and multi-phase (5; § \ref{['sec:zone5.multiphase.disk']}) optically-thin disks, central thermalized and blackbody-emitting disk (6; § \ref{['sec:zone6.disk']}), corona (7; § \ref{['sec:zone7.corona']}), and extended scattering layers and surfaces (8; § \ref{['sec:zone8.warm.scattering']}).
• Figure 3: Predicted disk surface densities $\Sigma_{\rm gas}$ versus cylindrical radius $R$ (in units of gravitational radius $R_{g}$ on the bottom axis, or light-days on the top axis). Each panel shows a single disk example (value of BH mass $M_{\rm BH}$ and accretion rate $\dot{m}$), plotted from the ISCO ($R_{\rm ISCO}$) to the BHROI ($R_{\rm BHROI}$R), for the flux-frozen magnetized ($\beta \ll 1$) models here ( black solid line). We contrast the prediction for the standard SS73 thermal-pressure-dominated $\alpha$-disk model ( pink dotted line). We compare different parameters: ${\rm M}_{\rm BH}=10^{7}\,{\rm M}_{\odot}$, $\dot{m}=1$ ( top left); ${\rm M}_{\rm BH}=10^{9}\,{\rm M}_{\odot}$, $\dot{m}=1$ ( top right); ${\rm M}_{\rm BH}=10^{7}\,{\rm M}_{\odot}$, $\dot{m}=0.01$ ( bottom left); ${\rm M}_{\rm BH}=10^{7}\,{\rm M}_{\odot}$, $\dot{m}=100$ ( bottom right). Shaded bars along the horizontal axis represent the locations of different zones (colored as Figs. \ref{['fig:cartoon.definitions']}-\ref{['fig:cartoon.quadrants']}) for each case, with their name ( bottom right) and dividing radii ( bottom left) labeled. Their $y$-axis values are not meaningful (they only highlight the range of radii for each zone). Bottom row shows the midplane ($|z|\ll H$) zones: thermalized disk (6a), ionized (5) and neutral (4) optically-thin disks, and obscuring torus (2; with 2b the thermalized-dust subregion). The row above corresponds to the disk illuminated surface ($|z|\sim H$) zones: warm comptonizing skin (6b), the scattering/reprocessing (8) and optically-thin ionized illuminated disk (3) BLR-like region; and illuminated warm dust-reprocessing torus (2c). Above this we show the range of radii of the coronal gas ($|z| \gtrsim H$), zone (7). Galactic ISM (Zone (1)) resides to the right of the plot ($R > R_{\rm BHROI}$) and is not modeled here. At most radii, masses, and accretion rates, the proportionally much stronger Maxwell stresses in the flux-frozen disks translate to lower $\Sigma_{\rm gas}$ compared to a thermal-pressure-dominated disk. Note the declining central densities in SS73 at high $\dot{m}$ are a consequence of the inner disks becoming radiation-pressure dominated. This only influences the models here weakly, as there is a small change in slope in the $\dot{m}=100$ flux-frozen case at $x_{g} \equiv R/R_{g} \sim 1500$ interior to which the disk becomes radiatively inefficient (§ \ref{['sec:super.eddington']}).
• Figure 4: Predicted midplane 3D gas densities $n\equiv \rho/m_{p}$, as Fig. \ref{['fig:surface.densities']}. We compare the same BH mass and accretion rates (panels), range of radii (ISCO to BHROI), zone locations (shaded horizontal bars), and flux-frozen disk ( thick black) versus thermal-pressure-dominated (SS73; thin pink) prediction. SS73-like models assume disks are weakly turbulent/laminar, so there is only one midplane density to plot. Flux-frozen disks are supersonically turbulent and multi-phase, so we plot our estimate of the volume-weighted mean midplane density $\langle \rho \rangle_{\rm vol} \sim \Sigma_{\rm gas}/(2\,H)$,; mass-weighted mean midplane density $\langle \rho \rangle_{\rm mass} \equiv M^{-1} \int \rho dm \approx \langle \rho \rangle_{\rm vol} C^{1/2}$ where $C$ is the "clumping factor"; and mass-weighted $\pm 2\sigma$ range of densities ( grey shaded). We assume a lognormal density distribution with the standard variance-Mach number relation $S \approx \ln[1 + (\mathcal{M}_{s}/3)^{2}]$ for supersonic turbulence to compute these konstantin:mach.compressive.relation. The discontinuities in $\langle \rho \rangle_{\rm mass}$ appear at midplane zone transitions because we use the analytic approximations from § \ref{['sec:zone1.ism']}-\ref{['sec:zone8.warm.scattering']} for the temperature/phase structure and dominant opacities, hence sound speed and $\mathcal{M}_{s}$. The mean profile crudely follows $\rho \propto R^{-2}$, with a broad range of densities (i.e. "clumpy"/inhomogeneous structure) at all radii. In general the densities are orders-of-magnitude lower than in an SS73 disk (owing both to lower $\Sigma_{\rm gas}$, Fig. \ref{['fig:surface.densities']}, and thicker disks with $H/R \sim 0.1-1$), except for interior regions at large $\dot{m}$ where radiation pressure modifies the magnetized solutions much more weakly. Note the similarity of the densities for flux-frozen disks to observational inferences for the BLR (§ \ref{['sec:zone3.blr']}) and dusty torus (§ \ref{['sec:zone2.torus']}) at corresponding radii (where the SS73 model is many orders-of-magnitude more dense).
• Figure 5: Total enclosed disk gas mass $M_{\rm gas,\,disk}(<R) \equiv \int_{0}^{R} 2\pi\Sigma_{\rm gas}\,R\,dR$ versus radii $R$ for the flux-frozen models ( thick black) or thermal-pressure-dominated models ( thin pink) with different BH masses and accretion rates (linestyles labeled). Per Fig. \ref{['fig:surface.densities']}, the flux-frozen disks predict much lower-disk masses. The gas masses required for thermal pressure-dominated disks become enormous at large radii, which in turn raises many theoretical challenges goodman:qso.disk.selfgrav, and even comparable to the BH mass (so Keplerian approximations break down, the disk cannot stabilize, and orbits become strongly perturbed) at $\gtrsim 1000\,R_{g}$. In contrast the flux-frozen disks require far less disk mass for the same BH mass and accretion rate/luminosity.