Masers and Broad-Line Mapping Favor Magnetically-Dominated AGN Accretion Disks

TL;DR

This work confronts the long-standing assumption that AGN outer accretion disks are thermally or radiation-pressure dominated by using direct kinematic probes from maser and BLR observations. By linking the observed near-Keplerian rotation and constrained disk masses to pressure support via $V_c^2/K$, the authors show that thermal, radiation, cosmic-ray, and simple turbulence-dominated disks predict disk masses and rotation curves that conflict with data, and would require unphysical temperatures or luminosities. In contrast, magnetic-pressure dominated (flux-frozen) disks naturally satisfy the dynamical constraints, align with maser Zeeman-field limits, and reproduce BLR/maser gas properties, implying outer disks are in a hyper-magnetized state with $P_{ m mag} \,\gg\, P_{ m therm}$ but $M_{\rm disk}(<R) \ll M_{\rm BH}$. The findings significantly shift the favored picture of SMBH fueling and disk structure, with magnetic stresses driving accretion in the outer regions and reducing the need for large, self-gravitating disks. Together, these results refine accretion disk theory and offer concrete observational discriminants for future high-resolution kinematic studies.

Abstract

We present a novel and powerful constraint on the physics of supermassive black hole (BH) accretion disks. We show that in the outer disk (radii $R \gtrsim 0.01\,$pc or $\gtrsim 1000\,R_{G}$), models supported by thermal or radiation pressure predict disk masses which are much larger than the BH mass and increase with radius - i.e. rapidly-rising, extremely non-Keplerian rotation curves. More generally, we show that any observational upper limit to the deviation from Keplerian potentials at these radii directly constrains the physical form of the pressure in disks. We then show that existing maser and broad line region (BLR) kinematic observations immediately rule out the classic thermal-pressure-dominated Shakura Sunyaev-like $α$-disk model, and indeed rule out any thermal or radiation (or cosmic-ray) pressure-dominated disk, as the required temperatures and luminosities of the gas at large radii would exceed those observed by orders of magnitude. We show that models where the pressure comes entirely from turbulence (without thermal, radiation, or magnetic sources) could in principle be viable but would require turbulent Toomre $Q \gtrsim 100$, far larger than predicted by self gravitating/gravito-turbulent models. However, recently proposed models of magnetic pressure-dominated disks agree with all of the observational constraints. These magnetically-dominated models also appear to agree better with constraints on maser magnetic fields, compared to the other possibilities. Observations appear to strongly favor the hypothesis that the outer regions of BH accretion disks are in the 'hyper-magnetized' state.

Figures (5)

• Figure 1: Expected scaling of gas circular velocity with radius and its deviation from Keplerian motion for different accretion disk models. Lines correspond to standard accretion disk models assuming a $M_{\rm BH}=10^{7}\,{\rm M_{\odot}}$ accreting near-Eddington ($m_{7}=\dot{m}=1$) with an effective accretion-stress-to-total-pressure ratio $\alpha=0.1$. We assume the disk total pressure is primarily: (1) thermal (a standard SS73-like or $\alpha$-disk; § \ref{['sec:thermal']}); (2) magnetic (a flux-frozen or hyper-magnetized disk; § \ref{['sec:magnetic']}); (3) radiation (a slim or radiation-pressure-supported disk; § \ref{['sec:radiation']}); (4) turbulent (constant-$Q_{0}$ or gravito-turbulent, with $Q_{0}=0.04$ as required to fit some NGC 1068 observations with this model; § \ref{['sec:turb']}). For each, we plot the predicted $V_{c}^{2}(R) \equiv G M_{\rm enc}(<R)/R = G\,[M_{\rm BH} + M_{\rm disk}(<R)]/R$. We compare to Keplerian (in absolute units at left or relative right). Per § \ref{['sec:theory']}, different assumptions about what dominates the pressure in the outer disk, for a given BH mass and luminosity, give very different predictions for the disk mass, with larger disk masses producing large deviations from Keplerian circular velocity curves (even rapidly-rising curves). Attempting to "re-tune" the thermal or radiation models by arbitrarily changing the predicted temperatures or opacities to suppress the deviations from Keplerian lead to other, immediately ruled-out predictions (like $\sim 10$ orders of magnitude larger luminosities).
• Figure 2: Observations (color; § \ref{['sec:obs']}), compared to the models (grayscale) per Fig. \ref{['fig:vc']}. We compile constraints from kinematics/rotation-curve-fitting from compilations of masers (shaded); GRAVITY optical interferometry of the BLR (hatched); microlensing response in broad-line wings (circles); plus individual resolved neutral-gas imaging and gas mass measurements in Circinus (star) and coordinated reverberation mapping of different BLR lines (squares). We compare model predictions assuming different midplane pressure sources as Fig. \ref{['fig:vc']}, for $\alpha=0.1$, $m_{7}=1$ and $\dot{m}=1$ (thick) or $\dot{m}=0.1$ (thin). Left: Deviation from Keplerian circular velocities. Right: Logarithmic slope of the circular velocity curve, fitted to each over the range of measured velocities. The best masers are within $\sim 1\%$ of Keplerian, while the best BLR measurements constrain the profile to within $\sim 5-10\%$ of Keplerian, sufficient to rule out radiation or thermal-pressure dominated models, or turbulence-dominated models with $Q \lesssim 1$ at $r \gtrsim 0.003\,$pc, but consistent with magnetic models at these radii. There is a hint of steeper-than-Keplerian curves in some BLR and maser systems, which only occurs in the magnetic models.
• Figure 3: Corresponding predictions (as Fig. \ref{['fig:vc']}) and lower/upper limits (Eqs. \ref{['eqn:HRmin']}-\ref{['eqn:rhomax']}) purely from the kinematic constraints in Fig. \ref{['fig:vc.obs']} on the central accretion disk scale-height $H/R$ (top) and volume-averaged midplane density $n_{\rm mid}$ (bottom). These constraints do not depend on the tracers being "part of" the disk, only that they feel the gravitational potential of the disk (and therefore constrain its mass). For the magnetically-dominated disk, since the disks are predicted to be highly multi-phase, we show the expected range of densities (shaded). We also show (ellipses) the approximate typical height and radius range of the BLR and "torus" regions, from typical covering factors and resolved imaging, as well as the typical densities of BLR clouds and of masing gas (from spectral modeling). The kinematic constraints clearly favor the magnetized disk model at $\gtrsim 0.003\,$pc. The limits on $H/R$ and density imply that the BLR, masers, and torus do not "sit above" a much denser, thinner disk containing most of the mass, but rather are fundamentally part of the same structure.
• Figure 4: Corresponding predictions as (Fig. \ref{['fig:vc']}) for the midplane in-plane magnetic field strengths $\langle | B_{\phi}|\rangle$ (assuming a Maxwell stress comparable to Reynolds/accretion stress for each model). We overplot upper limits from maser Zeeman constraints (Table \ref{['tbl:Bfields']}). Magnetically-dominated disks actually predict the weakest absolute $|{\bf B}|$ of any of the models, owing to their much lower overall densities (Fig. \ref{['fig:HR.rho']}). This appears to be clearly favored by the maser constraints at $\gtrsim 0.1\,$pc.
• Figure 5: Observational limits and predictions (as Fig. \ref{['fig:vc.obs']}) from different models for the minimum Toomre $Q(r)$ at a given radius $r$ in the disk (§ \ref{['sec:turb:grav']}), from kinematics alone (Eq. \ref{['eqn:Qmin']}). We plot the "turbulent" line at $Q=Q_{0}\sim1$, the usual prediction for gravito-turbulent or self-regulating models. Below $\sim$ pc scales, maser data and some BLR constraints require a minimum $Q \gtrsim 20-1000$, much larger than gravitoturbulent predictions, and many orders of magnitude larger than predicted for thermal or radiation-pressure supported disks, but consistent with magnetically-supported disk models.