---
title: Viscous Alpha-Disk Prescription
url: https://www.emergentmind.com/topics/viscous-alpha-disk-prescription
type: topic
---

# Viscous Alpha-Disk Prescription

The viscous alpha-disk prescription is a phenomenological framework central to the modeling of angular momentum transport, mass accretion, and disk evolution in various astrophysical disks, including protoplanetary disks, Be star decretion disks, active galactic nuclei (AGN), X-ray binaries, and self-gravitating systems. Originating with the Shakura–Sunyaev ansatz, this prescription parametrizes the turbulent stress tensor via a dimensionless efficiency parameter, $\alpha$, encoding the effects of unresolved microphysics—typically magnetohydrodynamic (MHD) turbulence—allowing construction of tractable time-dependent or steady-state disk models. The alpha-disk formalism is implemented across a wide range of regimes, from classical, thin, Keplerian disks to fully self-gravitating and wind-influenced systems.

## 1. Foundations and Key Equations

The core assumption of the viscous alpha-disk framework is that the kinematic viscosity, $\nu(r,t)$, mediating radial angular momentum transport, can be written as
\[
\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)
\]
where $c_s(r, t)$ is the local isothermal or adiabatic sound speed and $H(r, t)$ is the local disk scale height. In typical thin-disk (non-self-gravitating) regimes, $H = c_s / \Omega_K$, with $\Omega_K$ the Keplerian angular velocity. The viscous stress tensor is parameterized as $T_{r\varphi} = -\alpha\,P$, directly relating the turbulent stress to local pressure [1203.6851].

The viscous evolution of the disk surface density, $\Sigma(r,t)$, follows the diffusion equation:
\[
\frac{\partial \Sigma}{\partial t} = \frac{3}{r} \frac{\partial}{\partial r} \left[ r^{1/2} \frac{\partial}{\partial r} \left(\nu(r)\, \Sigma(r,t)\, r^{1/2}\right)\right]
\]
This relation, first formalized for accretion disks by Lynden-Bell & Pringle, underpins time-dependent simulations and is widely adopted in 1D hydrodynamic evolution codes such as SINGLEBE [1702.06982].

## 2. Physical Assumptions and Parameter Regimes

The standard alpha-disk prescription is built upon the following typical assumptions [1203.6851, 1210.0470]:
- The disk is thin ($H/r\ll1$) and axisymmetric.
- Local vertical hydrostatic equilibrium governs the vertical structure.
- Turbulent stresses responsible for angular momentum transport are proportional to total (gas plus, if relevant, radiation) pressure.
- The dominant balance is between viscous heating and radiative cooling, with negligible radial advection for steady-state disks.
- The dimensionless $\alpha$ parameter is sub-unity to ensure turbulence remains subsonic and correlates with the disk's ability to transport angular momentum efficiently.

Empirical and simulation-guided values for $\alpha$ span a broad range: $\alpha\sim10^{-4}$–$10^{-1}$ in protoplanetary and Be star disks, and up to order unity during Be outbursts [1702.06982, 2009.03323]. In MHD simulations of relativistic disks, $\alpha$ typically increases toward the ISCO due to enhanced shear and mean-field effects [1211.0526].

## 3. Extensions: Spatial and Temporal Variability of Alpha

While the original prescription treated $\alpha$ as constant, recent GRMHD simulations and analytic work have demonstrated the need for radius-dependent models. In particular, the radial profile of $\alpha(r)$ can be decomposed as a sum of turbulent and mean-field contributions:
\[
\alpha(r) = \alpha_0\,\left[\frac{q(r)}{1.5}\right]^n - \alpha_1\,\frac{b_{\hat{r}} b_{\hat{\varphi}}}{\rho^\Gamma}
\]
Here, $q(r)$ is the dimensionless shear parameter, analytically specified in the Kerr metric, and $n\approx6$ encodes the sensitivity of MRI-driven turbulence to the local shear [1211.0526]. A widely used fit is
\[
\alpha(r) = 0.025\,[q(r)/1.5]^6
\]
which captures the observed rise of $\alpha$ in the inner relativistic disk compared to the Newtonian regime.

In time-dependent Be disk models, $\alpha$ is implemented as piecewise-constant: held fixed within each phase (e.g., outburst or quiescence) but allowed to jump between segments. Observational fits require larger $\alpha$ during disk buildup ($\alpha\sim0.8$–$1$) and lower values during dissipation ($\alpha\sim0.1$–$0.2$) [1702.06982].

## 4. Application to Special Disk Environments

### a) Self-Gravitating and Gravitoturbulent Disks

When disk self-gravity becomes important (e.g., outer protoplanetary disks, AGN feeding), vertical support and viscosity scaling are modified. The vertical scale height transitions from $H = c_s/\Omega$ to $H = c_s^2/(\pi G \Sigma)$, and the viscosity takes the form $\nu\propto \alpha\,c_s^3/(\pi G\Sigma)$ [1602.04069]. Gravitational instability (Toomre $Q\lesssim1$) leads to a gravitoturbulent regime, with effective $\alpha$ determined by equating viscous heating to radiative cooling under $Q=Q_0\sim1$:
\[
\alpha_{gt}^R(r, \Sigma) = \frac{8}{9}\frac{\sigma ( \pi G Q_0 )^6}{f(\tau_Q)}\left(\frac{\mu}{k_B}\right)^4 \frac{\Sigma^5}{\Omega^7}\left[ 1 - \left(\frac{T_{\rm irr}}{T_Q}\right)^4\right]
\]
This yields a unique, local viscous stress for disks maintaining marginal gravitational stability [1501.04980].

### b) MRI-Driven and Layered Accretion

In protostellar disks with strong non-ideal MHD effects, MRI activity is restricted to surface layers, leading to a highly stratified alpha profile:
\[
\alpha(R,z) = 
\begin{cases}
1/(2\beta(R,z)), & \Lambda\geq1,\ \beta>\beta_{\rm min}(Am)\\
\alpha_{\min}\ll1, & \text{otherwise}
\end{cases}
\]
with the effective, vertically averaged $\alpha$ set by the fractional mass of MRI-active gas ($\Sigma_{\rm active}\sim 10\,{\rm g\,cm}^{-2}$ common in T Tauri systems) [1305.0770].

### c) Inclusion of Disk Winds

Strong magnetocentrifugal winds or magnetothermal outflows are not captured by a pure alpha prescription. Observational SED modeling finds that fitting far-IR fluxes often requires $\alpha\sim10^{-2}$, in tension with direct turbulent constraints from ALMA line measurements ($\alpha\lesssim10^{-3}$). This suggests that disk winds, not turbulence, may dominate angular momentum extraction in many systems [2009.03323].

## 5. Implementation in Disk Evolution Models

The alpha-disk prescription is broadly implemented in both steady-state and time-dependent models. A typical workflow includes:
- Solving for self-consistent $T(r)$ and $H(r)$ via vertical structure, radiative equilibrium, and energy balance equations.
- Using $\nu(r)=\alpha\,c_s\,H$ along with the continuity and angular momentum equations (or diffusion equation) to evolve $\Sigma(r, t)$.
- For global simulations, accommodating boundary conditions: specified mass-flux at the inner rim, zero-torque or outflow at outer boundaries, or mass-loss due to winds [1702.06982, 1210.0470].
- For fits to observational data (e.g., SEDs), embedding forward models within Bayesian inference frameworks, often with acceleration by artificial neural networks, to capture degeneracies between $\alpha$, $\dot{M}$, dust properties, and disk geometry [2009.03323].

In addition, variations such as the delayed-heating prescription introduce a timescale $\tau$ between the pressure and stress response. This can stabilize radiation-pressure-dominated regions otherwise thermally unstable under the instantaneous alpha law [1106.2335].

## 6. Observational and Theoretical Constraints

Empirical constraints on $\alpha$ arise from diverse sources:
- Light curve modeling in Be stars requires $\alpha$ variations of an order of magnitude between build-up and dissipation phases ($\alpha=0.1$–$1.0$ for $\omega$ CMa) [1702.06982].
- Classical accretion rates and SED modeling in T Tauri disks yield $\alpha\sim10^{-2}$, but direct turbulence measurements from line broadening suggest much lower values ($\alpha\lesssim10^{-3}$) [2009.03323].
- Disk outburst models (e.g., dwarf novae, FU Ori) require switching between low and high $\alpha$ values to match observed viscous timescales [1203.6851].
- In MHD simulations, $\alpha$ rises in relativistic (ISCO/proximal) regions, with fits such as $\alpha(r) = 0.025[q(r)/1.5]^6$ reproducing simulation results to within $\sim20\%$ [1211.0526].

These findings support both the flexibility and the limitations of the alpha prescription, motivating continuous development of more physically-grounded viscosity closures and hybrid models that can accommodate observed phenomena beyond pure turbulence, such as disk winds and gravito- or magneto-turbulent transport.

## 7. Limitations and Alternative Prescriptions

Despite its ubiquity, the alpha-disk model is recognized as fundamentally ad hoc, with no microscopic derivation of $\alpha$ in the original framework. Major limitations include [1203.6851]:
- Inability to self-consistently model transition to or the properties of dead zones, laminar flows, and regions of steep radial pressure gradients.
- Predicted thermal and viscous instabilities in radiation-pressure-dominated disks, mitigated only by additional physical ingredients (e.g., delayed heating).
- Lack of incorporation of MHD wind torques and non-local transport, both of which are increasingly favored for explaining observed disk behavior.
- In self-gravitating outer disks, the local alpha law may overestimate the transport rate compared to global torques.

Alternative viscosity prescriptions, such as the $\beta$-disk (Reynolds-number regulated, $\nu = \beta\,r^2\,\Omega$) or direct MHD closures using Maxwell and Reynolds stresses, offer physically distinct approaches but at the expense of analytic simplicity. Hybrid models combine viscous ($\alpha$) and wind-driven ($\beta$ or explicit magnetic field) effects to more completely describe angular momentum extraction in astrophysical disks [1501.04980, 1210.0470].

---

**References:**  
[1702.06982], [2009.03323], [1211.0526], [1305.0770], [1203.6851], [1501.04980], [1106.2335], [1210.0470], [1602.04069]

Source: https://www.emergentmind.com/topics/viscous-alpha-disk-prescription