Accretion
The Shakura-Sunyaev Alpha-Disk: How Turbulent Viscosity Sets the Accretion Rate
In 1973, two Soviet physicists compressed the impossibly messy problem of a swirling disk of gas falling onto a black hole into a single Greek letter. That letter, α (alpha), hides all the unknown turbulence in a number that most astronomers set to about 0.01–0.1 — and with it you can predict the temperature, brightness, and lifetime of an accretion disk from a white dwarf binary to a billion-solar-mass quasar.
The Shakura-Sunyaev alpha-disk is the standard analytic model for a geometrically thin, optically thick accretion disk. Its central insight is a viscosity prescription — ν = α c_s H — that parametrizes the unknown turbulent stress transporting angular momentum outward, allowing gas to spiral inward and release gravitational energy as radiation. It remains the workhorse framework of accretion astrophysics more than fifty years later.
- TypeGeometrically thin, optically thick accretion disk model
- IntroducedShakura & Sunyaev 1973 (A&A, vol. 24, p. 337)
- Key equationν = α c_s H (turbulent viscosity prescription)
- Typical alphaα ≈ 0.01–0.1 (MRI simulations give ~0.01; dwarf novae need ~0.1–0.3)
- RegimeSub-Eddington, radiatively efficient (η ≈ 0.06–0.42)
- Observed inCataclysmic variables, X-ray binaries, AGN/quasars, protoplanetary disks
Interactive visualization
Press play, or step through manually. The visualization is yours to drive — try it before reading on.
Watch the 60-second explainer
A condensed visual walkthrough — narrated, captioned, under a minute.
What the Alpha-Disk Is: The Angular-Momentum Problem
Gas in orbit around a compact object cannot simply fall in — it carries angular momentum. To accrete, a gas element must shed that angular momentum to material farther out. In a Keplerian disk (where orbital velocity ∝ R^(−1/2)) inner rings rotate faster than outer rings, so there is a shear. Friction across that shear transports angular momentum outward and lets mass drift inward, converting orbital energy into heat and light.
The problem Nikolai Shakura and Rashid Sunyaev faced in 1973 is that ordinary molecular viscosity in ionized disk gas is laughably weak — it would take longer than the age of the universe to accrete a disk. The real transport must come from turbulence, whose microphysics was unknown.
- Their solution was pragmatic: don't derive the turbulence, parametrize it.
- They assumed the disk is geometrically thin (vertical scale height H ≪ radius R) and optically thick, so it radiates locally as a blackbody.
- All the microphysical ignorance is bundled into one dimensionless number, α.
The Mechanism: The ν = α c_s H Prescription
The turbulent viscous stress cannot exceed the disk's pressure, and turbulent eddies cannot move faster than the sound speed (supersonic motion would shock and dissipate) nor be larger than the disk thickness. Shakura and Sunyaev captured both bounds by writing the kinematic viscosity as:
ν = α c_s H
where c_s is the local sound speed, H is the vertical scale height, and 0 < α ≲ 1. Equivalently, the vertically-integrated shear stress is τ = α P — the stress is a fixed fraction α of the pressure.
- Because H ≈ c_s / Ω (hydrostatic balance) and the viscous timescale scales as R²/ν, larger α means faster inflow.
- The steady-state accretion rate Ṁ is set by matching viscous heating to radiative cooling; α controls how quickly angular momentum diffuses.
- A landmark result is that the emitted flux is independent of the viscosity's details: dissipation per unit area D(R) ∝ Ṁ Ω², so the disk temperature profile T(R) ∝ R^(−3/4) does not depend on α.
Alpha sets the rate and structure (density, thickness), but the broadband spectral shape is remarkably robust.
Key Quantities and a Worked Example
The disk temperature follows the multicolor blackbody law T(R) = [3 G M Ṁ / (8π σ R³)]^(1/4) · f (with f accounting for the inner boundary), giving T ∝ R^(−3/4). Consider a stellar-mass black hole:
- M = 10 M_sun, accreting near Ṁ ≈ 10^(−8) M_sun/yr (~10% Eddington).
- The inner disk (a few Schwarzschild radii, ~90 km) reaches T ≈ 10^7 K (~1 keV) — hence these systems shine in soft X-rays.
- Total luminosity L = η Ṁ c², with efficiency η ≈ 0.06 (Schwarzschild) to 0.42 (maximal Kerr).
For a supermassive black hole (M = 10^8 M_sun), the same relation gives peak inner temperatures of only ~10^5 K, so quasar disks radiate in the ultraviolet/optical — the famous "big blue bump." With α = 0.1 the viscous inflow time across a CV disk is days to weeks; with α = 0.01 an AGN disk viscous time reaches thousands of years, which is why quasar disks appear steady while dwarf-nova disks visibly cycle.
How It's Observed and Where It Appears
The alpha-disk is not directly imaged (except, arguably, in the Event Horizon Telescope images of M87* and Sgr A*), but its signatures are ubiquitous:
- Cataclysmic variables (white-dwarf accretors): the model reproduces their thermal-viscous limit cycles. Dwarf-nova outbursts require the disk to switch between low-α (~0.01–0.03) quiescence and high-α (~0.1–0.3) outburst — a direct empirical constraint on α.
- X-ray binaries: the soft/thermal state shows a multicolor-disk-blackbody spectrum matching T ∝ R^(−3/4); disk-fitting yields black-hole spin via the inner disk radius.
- Active galactic nuclei / quasars: the UV "big blue bump" is the integrated alpha-disk spectrum; microlensing of lensed quasars measures disk sizes consistent (though somewhat larger) than predicted.
- Protoplanetary disks: α ~ 10^(−4)–10^(−2) inferred from ALMA turbulence limits and accretion rates onto young stars.
Observationally inferred α values cluster around 0.1–0.4 in ionized disks — notably higher than the ~0.01 that local MRI simulations produce, a persistent tension.
Related Regimes: Slim Disks, ADAFs, and the MRI
The alpha-disk is one member of a family, valid only when the disk is thin and radiatively efficient:
- Slim disks: near or above the Eddington limit, radiation pressure puffs the disk (H/R ~ 1) and photons are advected inward before escaping ("photon trapping"), lowering efficiency. Relevant to ultraluminous X-ray sources and tidal disruption events.
- ADAF/RIAF (advection-dominated flows): at very low Ṁ (≪ 0.01 Ṁ_Edd) the gas can't radiate efficiently and carries its heat into the black hole; the flow becomes hot, quasi-spherical, and dim. This describes Sgr A*.
- The MRI: the deepest development came in 1991, when Steven Balbus and John Hawley showed the magnetorotational instability — a weak magnetic field plus outward-decreasing angular velocity — spontaneously generates the turbulence α was standing in for. MHD simulations recover an effective α ~ 0.01–0.1, giving the phenomenological parameter a physical foundation.
Alpha is thus best seen as an emergent transport coefficient, not a fundamental constant.
Significance, Famous Cases, and Open Questions
Few equations in astrophysics have been as productive as ν = α c_s H. Shakura & Sunyaev (1973) is one of the most-cited papers in the field, and its framework underpins interpretations of essentially every accreting system, including:
- Cygnus X-1, the first confirmed black-hole candidate, whose soft-state disk spectrum is a textbook alpha-disk fit.
- Quasar continuum modeling and black-hole mass/accretion-rate estimates across cosmic time.
Yet major questions remain unresolved:
- The α discrepancy: dwarf novae demand α ~ 0.1–0.4, while zero-net-flux MRI boxes give ~0.01. Net magnetic flux, convection, and disk winds may bridge the gap.
- Radiation-pressure instability: in luminous inner disks, dominance by radiation pressure makes the classic α-disk thermally and viscously unstable — yet observed disks look stable, a long-standing puzzle.
- Is α even constant? Simulations suggest it varies with radius, magnetization, and thermodynamic state, motivating variable-α(R) prescriptions.
The alpha-disk endures not because it is exact, but because it is the right level of abstraction — a single knob that captures the essential physics while gracefully containing our ignorance of turbulence.
| Model / regime | Geometry (H/R) | Accretion rate (ṁ = Ṁ/Ṁ_Edd) | Radiative efficiency | Where it applies |
|---|---|---|---|---|
| Shakura-Sunyaev thin disk | H/R ≪ 1 (~0.001–0.1) | ~0.01 to ~1 | High (~5–40%) | CVs, soft-state X-ray binaries, AGN |
| Slim disk | H/R ~ 0.1–1 | ~1 to few (super-Eddington) | Reduced (photon trapping) | ULXs, tidal disruption events |
| ADAF / RIAF | H/R ~ 1 (quasi-spherical) | ≪ 0.01 | Very low (advected into hole) | Sgr A*, low-luminosity AGN, quiescent XRBs |
| Radiation-dominated inner disk | H/R ~ 0.1 (thermally unstable) | ≳ 0.1 | High but variable | Inner regions of luminous XRBs/AGN |
Frequently asked questions
What does the alpha parameter physically represent?
Alpha is the ratio of the turbulent shear stress to the local pressure in the disk (τ = α P). It bundles all the unknown microphysics of turbulence into one dimensionless number between 0 and about 1. A larger α means more efficient angular-momentum transport and therefore faster accretion; the viscosity is written ν = α c_s H, where c_s is the sound speed and H the disk scale height.
Who were Shakura and Sunyaev?
Nikolai Shakura and Rashid Sunyaev are Soviet/Russian astrophysicists who introduced the model in a 1973 Astronomy & Astrophysics paper (volume 24, page 337). Sunyaev is also co-discoverer of the Sunyaev-Zeldovich effect. Their alpha-disk paper is among the most-cited works in all of astrophysics and defined the standard accretion-disk framework.
Why doesn't the disk temperature depend on alpha?
The local energy dissipation rate in a steady disk depends only on the accretion rate Ṁ and orbital frequency Ω, not on how the viscosity is parametrized — this is because in steady state the same Ṁ must flow through every radius regardless of α. So the temperature profile T(R) ∝ R^(−3/4) and the emitted spectrum are set by M, Ṁ, and R. Alpha instead controls the disk's density, thickness, and the timescale on which it evolves.
What is the connection between the alpha-disk and the MRI?
When Shakura and Sunyaev wrote their model in 1973, the source of the required turbulence was unknown. In 1991 Steven Balbus and John Hawley showed that the magnetorotational instability (MRI) — a weak magnetic field in a differentially rotating disk — naturally drives the turbulence. MHD simulations of the MRI produce an effective α of roughly 0.01–0.1, giving the phenomenological alpha parameter a first-principles physical origin.
What typical value does alpha take?
Numerical MRI simulations with zero net magnetic flux typically give α ≈ 0.01. However, observations of dwarf-nova outbursts require α ≈ 0.1–0.3 in the hot ionized state, and protoplanetary disks suggest α ~ 10^(−4)–10^(−2). This roughly order-of-magnitude discrepancy between simulations and observations is an active and unresolved research topic.
When does the standard alpha-disk model break down?
The thin alpha-disk assumes the disk is geometrically thin and radiatively efficient. It breaks down near or above the Eddington limit, where radiation pressure inflates the disk into a 'slim disk' with photon trapping, and at very low accretion rates, where the flow becomes a hot, quasi-spherical, advection-dominated flow (ADAF/RIAF), as seen in Sgr A*. Radiation-pressure-dominated inner regions are also formally thermally and viscously unstable in the classic model.