Advection/Diffusion of Large-Scale B-Field in Accretion Disks

Activity of the nuclei of galaxies and stellar mass systems involving disk accretion to black holes is thought to be due to (1) a small-scale turbulent magnetic field in the disk (due to the magneto-rotational instability or MRI) which gives a large viscosity enhancing accretion, and (2) a large-scale magnetic field which gives rise to matter outflows and/or electromagnetic jets from the disk which also enhances accretion. An important problem with this picture is that the enhanced viscosity is accompanied by an enhanced magnetic diffusivity which acts to prevent the build up of a significant large-scale field. Recent work has pointed out that the disk’s surface layers are non-turbulent and thus highly conducting (or non-diffusive) because the MRI is suppressed high in the disk where the magnetic and radiation pressures are larger than the thermal pressure. Here, we calculate the vertical ($z$) profiles of the stationary accretion flows (with radial and azimuthal components), and the profiles of the large-scale, magnetic field taking into account the turbulent viscosity and diffusivity due to the MRI and the fact that the turbulence vanishes at the surface of the disk. The stationary solutions we find indicate that a weak ($\beta > 1$) large-scale field does not diffuse away as suggested by earlier work.

💡 Research Summary

The paper addresses a long‑standing paradox in accretion‑disk theory: while magnetorotational instability (MRI) generates small‑scale turbulence that enhances the effective viscosity and thus the accretion rate, the same turbulence also produces a turbulent magnetic diffusivity that should quickly erase any large‑scale magnetic field needed to launch winds or relativistic jets. Recent numerical and observational work has shown that the upper layers of a thin disk become magnetically and radiatively dominated, suppressing the MRI and creating a highly conducting, essentially non‑turbulent “surface layer.” The authors incorporate this layered structure into a semi‑analytic model that solves the steady‑state, axisymmetric MHD equations for a vertically stratified disk.

Key ingredients of the model are: (i) a turbulent viscosity ν(z) and magnetic diffusivity η(z) that follow the standard α‑prescription within the MRI‑active interior but decay smoothly to zero at the disk surface; (ii) a set of coupled ordinary differential equations derived from the continuity, momentum, and induction equations, retaining both radial (v_r) and azimuthal (v_φ) velocity components as well as the vertical (B_z) and toroidal (B_φ) magnetic‑field components; (iii) boundary conditions that fix B_z at the mid‑plane, enforce perfect conductivity (E_φ=0) at the surface, and require vanishing radial flow at the surface.

By solving these equations numerically, the authors obtain vertical profiles of v_r(z), v_φ(z), B_z(z), and B_φ(z). Inside the turbulent zone the radial inflow and azimuthal shear are strong, and B_φ is amplified by differential rotation. Near the surface, however, the vanishing ν and η cause the flow to stall and the toroidal field to decay sharply. Crucially, because η→0 in the surface layer, the large‑scale magnetic flux is effectively “locked” there and cannot diffuse outward. The calculations show that even a weak large‑scale field with plasma‑β>1 (i.e., magnetic pressure below the gas pressure) can survive for many accretion timescales, contrary to earlier analytic estimates that predicted rapid diffusion.

The study therefore revises the conventional view that MRI‑driven turbulence inevitably destroys large‑scale fields. Instead, the existence of a non‑turbulent, highly conducting surface sheath provides a natural barrier to magnetic diffusion, allowing the disk to retain the vertical flux needed for magnetocentrifugal winds or Blandford‑Znajek‑type jets. The authors conclude that any realistic model of jet‑producing accretion disks must include this layered conductivity structure. They also suggest that future work should explore how the thickness and physical properties of the non‑turbulent layer (e.g., radiation pressure, magnetic pressure, ionization state) affect field retention, and that fully three‑dimensional MHD simulations are required to validate the semi‑analytic predictions and to connect them with observable jet signatures.