---
title: Magnetically Dominated, Thermally Stable Disk Models
url: https://www.emergentmind.com/topics/thermally-stable-magnetically-dominated-disk-models
type: topic
---

# Magnetically Dominated, Thermally Stable Disk Models

A thermally stable, magnetically dominated disk is an accretion flow in which magnetic pressure—typically from toroidal (azimuthal) fields—contributes a significant fraction (often the majority) of the total pressure in the disk, thereby modifying the classical relations between heating, cooling, and instability. These disks contrast with traditional gas- or radiation-pressure supported geometrically thin disks, in that their stability, structure, and variability properties are fundamentally altered by magnetic support and the associated dynamics.

## 1. Theoretical Framework and Governing Equations

Thermally stable magnetically dominated disks are governed by the full set of ideal magnetohydrodynamic (MHD) equations, sometimes including general relativistic and radiative transport effects, depending on context and domain (e.g., X-ray binaries, AGN, TDEs). The central physical variable is the midplane plasma beta,
\[
\beta \equiv \frac{P_{\rm gas}+P_{\rm rad}}{P_{\rm mag}},
\]
where $P_{\rm mag} = B^2/8\pi$ is the magnetic pressure, and $P_{\rm gas}, P_{\rm rad}$ are gas and radiation pressures, respectively. A key regime of interest is $\beta \lesssim 1$, and often $\beta \ll 1$ in “hyper-magnetized” disks [2310.04507, 2310.04506, 2407.00160, 2502.05268].

The equations of mass, momentum, induction, and energy (including radiative transfer) can be summarized, in conservative form, as:
- Continuity:
  \[
  \frac{\partial \rho}{\partial t} + \nabla\cdot(\rho \mathbf{v}) = 0
  \]
- Momentum (including magnetic stress):
  \[
  \frac{\partial (\rho\mathbf{v})}{\partial t} + \nabla\cdot\left[ \rho\mathbf{v}\mathbf{v} + P_{\rm tot}\mathbf{I} - \frac{\mathbf{B}\mathbf{B}}{4\pi} \right] = -\rho\nabla\Phi
  \]
- Induction:
  \[
  \frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B})
  \]
- Energy (if included):
  \[
  \frac{\partial E}{\partial t} + \nabla\cdot\left[ (E+P_{\rm tot})\mathbf{v} - \frac{(\mathbf{B}\cdot\mathbf{v})}{4\pi}\mathbf{B} + \mathbf{F}_{\rm rad} \right] = -\rho \mathbf{v}\cdot\nabla\Phi
  \]
Here, $E$ includes all relevant energy densities, $P_{\rm tot} = P_{\rm gas} + P_{\rm rad} + P_{\rm mag}$.

Angular momentum transport is set by the Maxwell and Reynolds stresses, not a prescribed $\alpha$-viscosity. In regions where magnetic pressure dominates, $P_{\rm tot} \approx P_{\rm mag}$, and vertical support becomes
\[
H \sim \frac{v_A}{\Omega}
\]
with $v_A = B/\sqrt{4\pi \rho}$ the Alfvén speed.

These equations can be closed with an appropriate energy equation for the disk midplane, or in some idealized models, a locally isothermal equation of state with fixed thermal scale height (used in convergence tests of MRI-sustained turbulence) [2505.12671].

## 2. Magnetically-Induced Thermal Stabilization: Physical Mechanisms

In classical radiation-pressure dominated Shakura–Sunyaev disks ($P_{\rm rad} \gg P_{\rm gas}$, $P_{\rm mag} \ll P_{\rm rad}$), the local viscous heating $Q^+ \propto P_{\rm rad}^2$ rises more steeply with temperature than the radiative cooling $Q^- \propto P_{\rm rad} \propto T^4$, leading to the well-known thermal instability [2209.03317, 1601.06785]. When magnetic pressure becomes non-negligible, this scaling flattens:
\[
\frac{\partial \ln Q^+}{\partial \ln P_{\rm rad}} = \frac{2 - \beta_r^{-1}}{1 - \beta_r^{-1}}
\]
with $\beta_r^{-1}=P_{\rm mag}/(P_{\rm rad}+P_{\rm mag})$ [2209.03317]. Therefore, for $P_{\rm mag} \gtrsim 0.5 P_{\rm tot}$ (i.e., $\beta_r^{-1} \gtrsim 0.5$), the heating response falls below the cooling slope, and the disk is stabilized against the classical thermal runaway.

Similarly, in optically thin, cooling-dominated two-temperature disks, linear analysis shows the critical magnetic pressure fraction is $\beta_m \equiv P_{\rm mag}/P \gtrsim 0.7$ for unconditional thermal stability [1501.01241]. The stabilizing influence is further enhanced by energy extraction associated with magnetically driven winds and X-ray coronae [2501.16087, 2309.06939, 1403.2305].

The essential stabilizing features are:
- Maxwell stress replaces $\alpha$-prescriptions, with $Q^+ \propto P_{\rm mag}$.
- Magnetic pressure is less sensitive to temperature than radiation pressure.
- Additional angular momentum and energy extraction via winds/coronae increases the net cooling slope [2501.16087, 2309.06939].

## 3. Numerical and Analytic Models: Global Simulations and Similarity Solutions

Thermally stable, magnetically dominated disks have been realized in both semi-analytic steady-state solutions and global GR(M)HD simulations [2310.04506, 2310.04507, 2209.03317, 1601.06785, 2502.05268, 2505.12671]. Their main features are summarized in the table:

| Model Type           | Pressure Dominance     | Key Diagnostic $\beta$           | Typical $H/R$      |
|----------------------|-----------------------|-----------------------------------|--------------------|
| Shakura–Sunyaev      | $P_{\rm rad}$, $P_{\rm gas}$ | $\gg 1$                        | $0.01-0.03$        |
| Magnetically dominated (“Flux-frozen”) | $P_{\rm mag}$            | $10^{-4} - 1$                  | $0.03-0.3$         |
| Puffy disk           | $P_{\rm mag} \sim P_{\rm rad}$ | $\sim 1$                    | $0.1-1$            |

In the “hyper-magnetized” or “flux-frozen” regime, radial accretion brings ISM magnetic flux into the disk, amplified to $B_\phi$ fields sufficient for $P_{\rm mag} \gg P_{\rm gas}, P_{\rm rad}$ across most radii [2310.04506, 2407.00160]. Turbulence is trans-Alfvénic ($v_{\rm turb} \sim v_A$), and the turbulent Mach number $\mathcal{M}_s \gg 1$ [2502.05268].

Key simulation results include:
- Magnetically dominated disks remain thin ($H/R \sim 0.1$) but are “puffed-up” compared to pure $\alpha$-disks due to magnetic support.
- Maxwell stresses drive accretion at $\alpha \sim 0.1$.
- Surface density, temperature, $B$-field, and other quantities follow analytically predicted scaling laws (e.g., $B_\phi^2/8\pi \propto R^{-2/3}$) [2310.04507, 2407.00160].
- Thermal timescale $t_{\rm th} \ll t_{\rm visc}$; $Q^+ \approx Q^-$ in global energy balance.
- Magnetic field topology is crucial: net vertical or quadrupolar fields are required for robust, sustained magnetic pressure; pure dipole or small-scale-loops are less effective [2209.03317].

The global stability properties are robust over $0.01\lesssim\dot M/\dot M_{\rm Edd}\lesssim1000$ and black hole masses $10 - 10^{10}M_\odot$ [2310.04506, 2407.00160, 2512.09017].

## 4. Effects of Magnetized Winds and Coronae

Magnetically driven winds and coronae play a dual role: they extract angular momentum and mass, and carry off a fraction of the dissipated energy, further steepening the effective cooling curve in the $Q^+$-$Q^-$ plane.

- The mass-loss rate and the wind angular momentum parameter (mass-loading $\mu$) are critical in one-dimensional simulations; stability is achieved when the wind cooling term $Q_{\rm w}^-$ grows rapidly enough with $T$ [2309.06939, 1403.2305].
- The coronal energy loss fraction $f_{\rm cor}$ acts analogously: as $f_{\rm cor}$ increases, the stabilizing effect is enhanced [2501.16087].
- Stability boundaries in the $f_{\rm cor}$–$s$ (wind index) plane are quantified: typically, $f_{\rm cor}\gtrsim0.2$, $s\gtrsim0.2$, and $\beta_{\rm mag}\gtrsim0.3$ suffice to stabilize the hot, magnetically dominated branch [2501.16087].

This mechanism provides a natural explanation for the absence of radiation-pressure-driven variability (“limit cycles”) in luminous X-ray binaries and AGN [1403.2305, 2309.06939].

## 5. Radial and Vertical Structure

On the magnetically dominated branch:
- Vertical support: $H/R \sim v_A/v_K \sim 0.03$–$0.3$, regulated by $P_{\rm mag}$ [2310.04507, 2407.00160].
- Surface density: $\Sigma \propto R^{-5/6}$ [2310.04507].
- Temperature: midplane $T$ grows inward but is regulated by rapid cooling; profiles remain shallow compared to radiative-pressure-dominated disks.
- Plasma $\beta$ never rises above unity except at the innermost radii or in the midplane after flux escapes [2505.12671].
- Turbulence is trans- to super-Alfvénic, with high Reynolds and Maxwell stresses ($\alpha \sim 0.1$).
- At high $\dot{m}$, the disk remains magnetically supported and thermally stable; radiation pressure saturates at $< P_{\rm mag}$ due to efficient vertical mixing [2502.05268].

## 6. Observational Implications and Applications

Thermally stable, magnetically dominated disks predict and are consistent with a variety of observed phenomena:
- Absence of thermal/viscous limit cycles in X-ray binary high/soft states even at $L\sim0.1-0.5\,L_{\rm Edd}$ [2209.03317, 1403.2305, 1601.06785, 2512.09017].
- Predicts broad UV/optical plateaus in TDEs and LFBOTs, with luminosities and durations explained by stable magnetically dominated disk spreading, not by classical gas- or radiation-pressure disks [2503.03811, 2512.09017].
- Provides a physical mechanism for the broad-line region, warm corona, and flared torus in AGN, all as different radial (and vertical) phases in a unified “multi-phase” magnetically dominated accretion structure [2407.00160, 2310.04506, 2502.05268].
- Jet powers via Blandford-Znajek-like mechanisms are set by the magnetic flux regulated in these disks, producing both radio-quiet and jetted TDEs within the observed range [2211.00704].
- Magnetically dominated disks remain optically thick and can reproduce the spectral hardening seen in AGN and TDEs, with characteristic color correction factors $f_{\rm col}\sim 1.4-2$ [2211.00704].

## 7. Methodological and Model-Building Considerations

Construction and diagnosis of thermally stable, magnetically dominated disks require:
- Seeding the disk with sufficient net vertical (or large-scale quadrupolar) magnetic flux to enable MRI amplification to $P_{\rm mag}/P_{\rm tot} \gtrsim 0.5$ [2209.03317, 1601.06785].
- Sufficient vertical resolution in simulations ($N_H \gtrsim 4$ cells per scale height) to capture MRI, Parker & buoyancy instabilities, and vertical escape [2505.12671].
- Accounting for flux escape, local dynamo, and boundary-driven advection for field replenishment—radial advection alone is insufficient inside the circularization radius, requiring vertical flux for sustained magnetization [2505.12671].
- Including energy loss via magnetically powered winds and/or coronae, with appropriate parameterizations for $f_{\rm wind}, f_{\rm cor}$, and wind/orbital parameters [2501.16087, 2309.06939, 1403.2305].

Stability is robust over a wide parameter space: $P_{\rm mag}\gtrsim 0.3\,P_{\rm tot}$, $f_{\rm wind}+f_{\rm cor}\gtrsim 0.1$, and moderate to strong net flux topologies guarantee stability across physically realistic conditions for AGN, XRBs, TDEs, and beyond [2501.16087, 2211.00704, 2310.04506].

---

**References**:  
- "The Role of Strong Magnetic Fields in Stabilizing Highly Luminous, Thin Disks" [2209.03317]  
- "Thin accretion disks are stabilized by a strong magnetic field" [1601.06785]  
- "Idealized Global Models of Accretion Disks with Strong Toroidal Magnetic Fields" [2505.12671]  
- "Thermal instability of thin disk in the presence of wind and corona" [2501.16087]  
- "Zooming In On The Multi-Phase Structure of Magnetically-Dominated Quasar Disks" [2502.05268]  
- "Multi-Phase Thermal Structure & The Origin of the Broad-Line Region, Torus, and Corona in Magnetically-Dominated Accretion Disks" [2407.00160]  
- "FORGE'd in FIRE II: The Formation of Magnetically-Dominated Quasar Accretion Disks from Cosmological Initial Conditions" [2310.04506]  
- "An Analytic Model For Magnetically-Dominated Accretion Disks" [2310.04507]  
- "Magnetically Dominated Disks in Tidal Disruption Events and Quasi-Periodic Eruptions" [2211.00704]  
- "Time-dependent global simulations of a thin accretion disc: the effects of magnetically-driven winds on thermal instability" [2309.06939]  
- "Thermal stability of thin disk with magnetically driven winds" [1403.2305]  
- "Thermal Stability of Magnetized, Optically Thin, Radiative Cooling-dominated Accretion Disks" [1501.01241]  
- "Puffy accretion disks: sub-Eddington, optically thick, and stable" [1908.08396]  
- "Late-Time Evolution of Magnetized Disks in Tidal Disruption Events" [2503.03811]  
- "Viscously Spreading Accretion Disks around Black Holes: Implications for TDEs, LFBOTs and other Transients" [2512.09017]

Source: https://www.emergentmind.com/topics/thermally-stable-magnetically-dominated-disk-models