Viscous Alpha-Disk Prescription

Updated 25 January 2026

Viscous alpha-disk prescription is a framework that parameterizes turbulent stresses in astrophysical disks using a dimensionless alpha parameter.
It is applied across various systems, including protoplanetary, AGN, and Be star disks, to model angular momentum transport and mass accretion.
Modern models incorporate spatial/temporal variations, disk winds, and gravitational instabilities to address limitations of the original, ad hoc formulation.

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 (Montesinos, 2012).

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 (Ghoreyshi et al., 2017).

2. Physical Assumptions and Parameter Regimes

The standard alpha-disk prescription is built upon the following typical assumptions (Montesinos, 2012, Abbassi et al., 2012):

The disk is thin ( $\nu(r,t)$ 0) 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 $\nu(r,t)$ 1 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 $\nu(r,t)$ 2 span a broad range: $\nu(r,t)$ 3– $\nu(r,t)$ 4 in protoplanetary and Be star disks, and up to order unity during Be outbursts (Ghoreyshi et al., 2017, Ribas et al., 2020). In MHD simulations of relativistic disks, $\nu(r,t)$ 5 typically increases toward the ISCO due to enhanced shear and mean-field effects (Penna et al., 2012).

3. Extensions: Spatial and Temporal Variability of Alpha

While the original prescription treated $\nu(r,t)$ 6 as constant, recent GRMHD simulations and analytic work have demonstrated the need for radius-dependent models. In particular, the radial profile of $\nu(r,t)$ 7 can be decomposed as a sum of turbulent and mean-field contributions: $\nu(r,t)$ 8 Here, $\nu(r,t)$ 9 is the dimensionless shear parameter, analytically specified in the Kerr metric, and $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 0 encodes the sensitivity of MRI-driven turbulence to the local shear (Penna et al., 2012). A widely used fit is

$\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 1

which captures the observed rise of $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 2 in the inner relativistic disk compared to the Newtonian regime.

In time-dependent Be disk models, $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 3 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 $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 4 during disk buildup ( $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 5– $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 6) and lower values during dissipation ( $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 7– $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 8) (Ghoreyshi et al., 2017).

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 $\nu(r, t) = \alpha\, c_s(r, t)\, H(r, t)$ 9 to $c_s(r, t)$ 0, and the viscosity takes the form $c_s(r, t)$ 1 (Kubsch et al., 2016). Gravitational instability (Toomre $c_s(r, t)$ 2) leads to a gravitoturbulent regime, with effective $c_s(r, t)$ 3 determined by equating viscous heating to radiative cooling under $c_s(r, t)$ 4: $c_s(r, t)$ 5 This yields a unique, local viscous stress for disks maintaining marginal gravitational stability (Rafikov, 2015).

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: $c_s(r, t)$ 6 with the effective, vertically averaged $c_s(r, t)$ 7 set by the fractional mass of MRI-active gas ( $c_s(r, t)$ 8 common in T Tauri systems) (Landry et al., 2013).

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 $c_s(r, t)$ 9, in tension with direct turbulent constraints from ALMA line measurements ( $H(r, t)$ 0). This suggests that disk winds, not turbulence, may dominate angular momentum extraction in many systems (Ribas et al., 2020).

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 $H(r, t)$ 1 and $H(r, t)$ 2 via vertical structure, radiative equilibrium, and energy balance equations.
Using $H(r, t)$ 3 along with the continuity and angular momentum equations (or diffusion equation) to evolve $H(r, t)$ 4.
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 (Ghoreyshi et al., 2017, Abbassi et al., 2012).
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 $H(r, t)$ 5, $H(r, t)$ 6, dust properties, and disk geometry (Ribas et al., 2020).

In addition, variations such as the delayed-heating prescription introduce a timescale $H(r, t)$ 7 between the pressure and stress response. This can stabilize radiation-pressure-dominated regions otherwise thermally unstable under the instantaneous alpha law (Ciesielski et al., 2011).

6. Observational and Theoretical Constraints

Empirical constraints on $H(r, t)$ 8 arise from diverse sources:

Light curve modeling in Be stars requires $H(r, t)$ 9 variations of an order of magnitude between build-up and dissipation phases ( $H = c_s / \Omega_K$ 0– $H = c_s / \Omega_K$ 1 for $H = c_s / \Omega_K$ 2 CMa) (Ghoreyshi et al., 2017).
Classical accretion rates and SED modeling in T Tauri disks yield $H = c_s / \Omega_K$ 3, but direct turbulence measurements from line broadening suggest much lower values ( $H = c_s / \Omega_K$ 4) (Ribas et al., 2020).
Disk outburst models (e.g., dwarf novae, FU Ori) require switching between low and high $H = c_s / \Omega_K$ 5 values to match observed viscous timescales (Montesinos, 2012).
In MHD simulations, $H = c_s / \Omega_K$ 6 rises in relativistic (ISCO/proximal) regions, with fits such as $H = c_s / \Omega_K$ 7 reproducing simulation results to within $H = c_s / \Omega_K$ 8 (Penna et al., 2012).

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 $H = c_s / \Omega_K$ 9 in the original framework. Major limitations include (Montesinos, 2012):

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 $\Omega_K$ 0-disk (Reynolds-number regulated, $\Omega_K$ 1) or direct MHD closures using Maxwell and Reynolds stresses, offer physically distinct approaches but at the expense of analytic simplicity. Hybrid models combine viscous ( $\Omega_K$ 2) and wind-driven ( $\Omega_K$ 3 or explicit magnetic field) effects to more completely describe angular momentum extraction in astrophysical disks (Rafikov, 2015, Abbassi et al., 2012).