Fact-checked by Grok 2 weeks ago

Accretion disk

An accretion disk is a flattened, rotating structure composed of gas, dust, and plasma that forms around a massive central object, such as a black hole, neutron star, protostar, or white dwarf, when surrounding material with non-zero angular momentum spirals inward under gravity.^[1]^[2] Due to the conservation of angular momentum, infalling matter cannot plunge directly toward the center but instead orbits in a disk-like configuration, gradually losing energy through viscous processes that cause it to migrate inward while heating to extreme temperatures.^[3]^[1] These disks are ubiquitous in astrophysical environments and play a pivotal role in phenomena ranging from star formation to high-energy emissions observed in the universe. In protoplanetary disks around young stars, cooler material at outer radii can coalesce into planets, while inner regions supply gas to the star itself.^[1] Around compact objects like black holes in binary systems, accretion disks form when gas is stripped from a companion star, leading to rapid orbital decay and eventual infall.^[2] The frictional heating in these disks—driven by turbulence, magnetic fields, and differential rotation—raises temperatures to hundreds of thousands or even millions of Kelvin, particularly in the inner regions, resulting in intense thermal radiation across X-rays and other wavelengths.^[3]^[1] Accretion disks are highly efficient at converting gravitational potential energy into electromagnetic radiation, with efficiencies up to 30 times greater than nuclear fusion in stars, making them the dominant light sources for active galactic nuclei, quasars, and X-ray binaries.^[3] For supermassive black holes at galactic centers, disks can power jets of relativistic particles and outshine the entire host galaxy's stars combined.^[1] Observations reveal disk properties through spectral signatures, including Doppler shifts from rotation and gravitational redshift near the central object, while simulations incorporate magnetohydrodynamics to model viscosity and angular momentum transport.^[3]^[2]

Definition and Formation

Basic Principles

An accretion disk is a flattened, rotating structure composed of gas, dust, and plasma that orbits a central gravitating body, such as a black hole, neutron star, or young star, where material spirals inward by dissipating angular momentum through viscous or turbulent processes, thereby converting gravitational potential energy into thermal energy and radiation. This concept was first proposed in the context of powering quasars through accretion onto massive objects by Salpeter (1964) and extended to galactic nuclei as collapsed quasars involving supermassive black holes by Lynden-Bell (1969).^[4]^[5]^[6] The basic structure of an accretion disk features a radial extent from an inner edge—often truncated at the innermost stable circular orbit (ISCO) for relativistic systems like black holes—to an outer radius set by the material supply. Vertically, the disk is typically thin, with a scale height H much smaller than the radial distance r (H/r \ll 1) in standard models, and is primarily composed of ionized gas at temperatures ranging from thousands to millions of Kelvin. The orbital motion is predominantly Keplerian, with the angular velocity given by

\Omega = \sqrt{\frac{GM}{r^3}},

where G is the gravitational constant, M is the central mass, and r is the radial coordinate.^[7]^[4] Key parameters governing the disk's dynamics include the mass accretion rate \dot{M}, which quantifies the inward mass flow; the central mass M; the dimensionless viscosity parameter \alpha introduced by Shakura and Sunyaev (1973) to parameterize turbulent stresses ($0 < \alpha < 1); the surface mass density \Sigma; and the midplane temperature T(r), which decreases outward. In thin disk approximations, the radial drift velocity v_r is small and negative (inward), approximated as

v_r \approx -\frac{3}{2} \alpha \left( \frac{H}{r} \right)^2 v_\phi,

where v_\phi \approx \sqrt{GM/r} is the azimuthal velocity. These parameters determine the disk's luminosity and spectral properties, with \alpha typically inferred to be around 0.01–0.1 from observations and simulations.^[8]^[7]

Formation Processes

Accretion disks form through several primary mechanisms depending on the astrophysical context, such as the capture of interstellar medium by a moving compact object, mass transfer via Roche lobe overflow in binary systems, wind accretion described by the Bondi-Hoyle process, and the gravitational collapse of molecular clouds during star formation. In the case of interstellar medium capture, diffuse gas is drawn into the gravitational well of a compact object like a black hole or neutron star, where differential rotation imparts angular momentum leading to disk buildup. Roche lobe overflow occurs when a companion star in a binary system expands to fill its Roche lobe, transferring material equatorially toward the accretor, which then circularizes into a disk due to the conservation of angular momentum.^[8] The Bondi-Hoyle mechanism applies to accretion from stellar winds or ambient flows, where the relative motion between the accretor and gas creates a shocked accretion column that can form a disk if the impact parameter provides sufficient angular momentum. During star formation, rotating molecular cloud cores collapse under self-gravity, channeling infalling material into a centrifugally supported disk around the central protostar.^[9] The initial conditions for disk formation typically involve material in a roughly spherical or toroidal configuration approaching the central gravitating body, such as a protostar or compact object. As this material falls inward, angular momentum is conserved, causing it to flatten into a disk-like structure perpendicular to the rotation axis. The centrifugal barrier arises from this conserved angular momentum, halting radial infall and preventing direct accretion onto the central object once the material reaches radii where rotational support balances gravity. For disk formation to occur, the specific angular momentum l of the incoming material must exceed the Keplerian value at the innermost stable circular orbit (ISCO), ensuring the material cannot plunge directly into the central object. The specific angular momentum is given by

l = r^2 \Omega,

where r is the radial distance and \Omega is the angular velocity. In the case of a Schwarzschild black hole, the condition for disk formation is l > l_{\rm ISCO} \approx 2\sqrt{3} \, GM/c, where G is the gravitational constant, M is the black hole mass, and c is the speed of light. The formation timescale is on the order of the dynamical time t_{\rm dyn} = 1/\Omega, which corresponds to the orbital period at the initial radius of collapse or infall, typically ranging from hours to years depending on the system scale. Once formed, the disk grows through continuous infall of additional material, with the inner disk boundary set by the centrifugal barrier and angular momentum distribution. Viscosity plays a brief role post-formation in allowing gradual inward migration of material while transporting angular momentum outward. Specific examples illustrate these processes: protoplanetary disks emerge from the gravitational collapse of rotating molecular cloud cores, where conserved angular momentum halts infall at radii of ~100 AU around young stars of solar mass.^[9] In X-ray binaries, mass transfer from a Roche lobe-filling companion to a neutron star or black hole forms a disk with inner radius determined by the specific angular momentum of the transferred stream, often leading to luminosities exceeding 10^{37} erg/s.^[8]

Physical Mechanisms

Angular Momentum and Viscosity

In accretion disks, material in near-Keplerian rotation around a central gravitating body possesses substantial specific angular momentum, on the order of \sqrt{GM r} per unit mass, where M is the central mass and r the radial distance. Without a mechanism to redistribute this angular momentum outward, the material cannot spiral inward to accrete, as conservation of angular momentum would otherwise maintain or increase the orbital radius. This "angular momentum problem" necessitates efficient transport processes to enable sustained accretion rates observed in astrophysical systems.^[10] Viscosity plays a central role in resolving this issue by generating torques from the shear arising between differentially rotating annuli in the disk. Molecular viscosity alone is insufficient due to its extremely low value in astrophysical plasmas, so an effective kinematic viscosity \nu must arise from turbulent motions. The Shakura-Sunyaev \alpha-disk prescription parameterizes this as \nu = \alpha c_s H, where c_s is the isothermal sound speed, H is the vertical scale height of the disk, and \alpha (typically $0.01 \lesssim \alpha \lesssim 0.1) quantifies the strength of the turbulence relative to orderly rotation. The turbulence is commonly sourced by the magnetorotational instability in magnetized disks. This formulation allows angular momentum to be transferred outward while permitting inward radial drift of the accreting material.^[8] The key viscous stress is the azimuthal-radial component of the stress tensor, t_{r\phi} = -\alpha P, where P is the total (gas plus radiation) pressure; this stress exerts torques that drive the transport. In steady-state conditions, angular momentum conservation integrates to yield the mass accretion rate \dot{M} = 3\pi \nu \Sigma, where \Sigma is the surface mass density, assuming \nu is independent of radius and the inner boundary torque is negligible. A brief derivation begins with the azimuthal angular momentum equation from the Navier-Stokes equations in cylindrical coordinates, balancing viscous torques with the advective flux of angular momentum. For a thin disk, this leads to the radial velocity

v_r = -\frac{3}{2} \frac{\nu}{r} \frac{d \ln (\nu \Sigma r^{1/2})}{d \ln r},

which, when combined with mass conservation \dot{M} = -2\pi r \Sigma v_r, recovers the accretion rate expression for constant \nu \Sigma. This inward velocity is small compared to the Keplerian speed, ensuring the disk remains geometrically thin.^[10] Although viscosity dominates in standard models, non-viscous mechanisms can contribute to angular momentum transport in certain regimes. Gravitational torques, for instance, arise in warped disks where Lense-Thirring precession from a rotating central black hole induces differential tilting, leading to the Bardeen-Petterson effect that aligns the inner disk while torquing the warp to redistribute momentum. Global modes, such as density waves excited by self-gravity or external perturbations, can also drive torques without relying on local viscosity. These alternatives are particularly relevant in dense, self-gravitating disks or those with significant misalignment.^[11]

Energy Generation and Radiative Processes

In accretion disks, the primary energy source is the release of gravitational potential energy as matter spirals inward toward the central object. For a thin disk around a black hole, the total gravitational energy liberated per unit mass accreted is approximately E_{\rm grav} \approx \frac{GM \dot{M}}{2 R_{\rm in}}, where M is the central mass, \dot{M} is the accretion rate, and R_{\rm in} is the inner radius of the disk. This corresponds to an accretion efficiency \eta \approx 1/12 (or about 8%) in the Newtonian approximation for thin disks, assuming the inner edge is at roughly three times the Schwarzschild radius. This gravitational energy is converted into thermal energy through viscous dissipation within the differentially rotating disk. The viscous heating rate per unit surface area is given by Q_+ = \frac{9}{4} \nu \Sigma \Omega^2, where \nu is the kinematic viscosity, \Sigma is the surface density, and \Omega is the Keplerian angular velocity. In steady-state conditions, this heating is balanced by radiative cooling, leading to a local thermal equilibrium. Cooling primarily occurs through radiation from the disk surfaces. In optically thick regions, the cooling rate per unit area is Q_- = 2 \sigma T_{\rm eff}^4, where \sigma is the Stefan-Boltzmann constant and T_{\rm eff} is the effective temperature; the factor of 2 accounts for emission from both sides of the disk. The spectrum in these regions approximates a blackbody, modified by the radial temperature gradient. For optically thinner cases or diffusion approximations, the cooling can be expressed as Q_- \approx 2 \sigma T_{\rm eff}^4 / \tau, with \tau the vertical optical depth.^[12] The effective temperature profile in a steady-state thin disk follows T(r) \propto r^{-3/4}, reflecting the radial dependence of the dissipation rate. More precisely, T_{\rm eff}^4 = \frac{3 G M \dot{M}}{8 \pi \sigma r^3} \left(1 - \sqrt{\frac{r_{\rm in}}{r}}\right), which accounts for the boundary condition of no torque at the inner edge. This results in a multicolour blackbody spectrum, with outer regions emitting in the infrared and inner regions dominated by ultraviolet and X-ray emission due to higher temperatures. Key opacity sources influencing the disk's radiative transfer include electron scattering, which dominates in ionized regions, and free-free absorption, important in cooler, denser parts of the disk. These processes determine the optical depth \tau and thus the transition between optically thick and thin regimes across different radial zones.

Theoretical Models

Thin Disk Models

Thin disk models provide analytic descriptions of geometrically thin and optically thick accretion disks, where viscous heating is balanced locally by radiative cooling, applicable to steady-state accretion at sub-Eddington rates onto compact objects such as black holes or neutron stars. These models assume Keplerian orbital motion in the disk, with the disk's scale height H satisfying H/r \ll 1, ensuring a slim geometry that allows efficient photon diffusion outward. Viscosity is parameterized via the \alpha-prescription, where the kinematic viscosity \nu = \alpha c_s H and $0 < \alpha < 1, representing turbulent stresses proportional to the total pressure; energy transport is dominated by radiation without significant radial advection.^[13]^[14] The Shakura-Sunyaev framework divides the disk into three radial regions based on the dominant pressure and opacity sources, reflecting variations in temperature and density with radius. In the inner region, radiation pressure dominates with electron-scattering opacity (\kappa = \kappa_{es} \approx 0.4 cm²/g). The middle region features gas pressure support with free-free opacity (\kappa_{ff}). The outer region is characterized by gas pressure and Kramers opacity (\kappa \propto \rho T^{-7/2}), where bound-free and free-free processes prevail at lower temperatures. These divisions arise from solving the vertical hydrostatic equilibrium and opacity expressions, with transitions occurring at specific radii depending on the accretion rate \dot{M} and central mass M.^[13]^[14] Key scaling relations emerge from balancing angular momentum transport and energy dissipation. In the middle region, the surface density scales as \Sigma \propto \dot{M}^{3/5} r^{3/5} \alpha^{-4/5}, reflecting the interplay of viscous diffusion and radiative transfer. The total disk luminosity follows L \propto \dot{M}, with an efficiency \eta \approx 0.06 for non-spinning black holes, determined by the binding energy released near the innermost stable circular orbit (ISCO). These scalings are derived under the assumptions of steady accretion and local thermodynamic equilibrium.^[13]^[15] The \alpha-disk solution involves integrating the full vertical structure to obtain density \rho(z) and temperature T(z) profiles, assuming hydrostatic equilibrium \frac{dP}{dz} = -\rho \Omega^2 z and radiative diffusion for energy transport. Boundary conditions include zero stress at the disk midplane and matching to the no-torque condition at the ISCO radius (r = 6GM/c^2 for Schwarzschild black holes), ensuring angular momentum conservation. The central temperature T_c and \Sigma profiles are determined by equating the vertically integrated viscous heating rate per unit area Q_+ to the radiative cooling rate Q_-:

Q_+ = Q_-

Here, Q_+ \approx \frac{9}{8} \nu \Sigma \Omega_K^2 for Keplerian angular velocity \Omega_K = \sqrt{GM/r^3}, while Q_- = \frac{8 \sigma T_c^4}{3 \tau} with optical depth \tau = \kappa \Sigma / 2; solving this yields T_c \propto (\dot{M} \alpha^{-1} r^{-3})^{1/4} in simplified forms, enabling predictions of spectral energy distributions.^[13]^[14] These models break down at low accretion rates (\dot{M} \lesssim 0.01 \dot{M}_{Edd}), where radiative cooling becomes inefficient and advection dominates energy transport, leading to hotter, geometrically thicker flows. At high rates (\dot{M} \gtrsim 0.4 \dot{M}_{Edd}), radiation pressure triggers thermal instabilities, and photon trapping reduces radiative efficiency, invalidating the thin-disk approximation.^[14]

Thick and Advection-Dominated Models

Thick accretion disks, where the vertical scale height is comparable to the radial extent, arise in regimes where radiative cooling is suppressed, leading to hot, puffed-up structures supported by gas or radiation pressure. In these models, energy generated by viscous dissipation is primarily transported inward via advection rather than radiated locally, resulting in low radiative efficiencies. This contrasts with thinner disks by emphasizing entropy advection as the dominant cooling mechanism, particularly in low-density environments around compact objects. A prominent example is the advection-dominated accretion flow (ADAF) model, developed by Narayan and Yi, which describes optically thin, hot flows reaching virial temperatures (~10^9-10^{12} K) with sub-Keplerian rotation velocities. In ADAFs, the plasma is collisionless and turbulent, driven by magnetic fields or other instabilities, and the advection parameter f_{\rm adv} \approx 1 - \epsilon_{\rm rad} quantifies the fraction of viscous heating advected versus radiated, where \epsilon_{\rm rad} is typically much less than 1 due to inefficient bremsstrahlung and synchrotron cooling.^[16] The "Polish doughnut" model, introduced by Abramowicz et al., represents an early analytic framework for such thick disks, featuring toroidal, pressure-supported configurations with constant specific angular momentum distribution. These structures maintain hydrostatic and vertical equilibrium without self-gravity, remaining stable even under super-Eddington accretion rates where thin-disk assumptions break down.^[17] For moderately high accretion rates, slim disk models by Abramowicz et al. extend these ideas to mildly thick, optically thick flows where advection competes with radiation pressure. Here, photon trapping becomes significant, reducing the effective optical depth and allowing super-Eddington luminosities without excessive inflation, as trapped photons are advected inward alongside the flow.^[18] In all these models, the radial advection of entropy balances the energy budget, expressed as the advection rate

q_{\rm adv} = f \Sigma T v_r \frac{ds}{dr},

where \Sigma is the surface density, T the temperature, v_r the radial velocity, s the specific entropy, and f a geometric factor (~1-2); this term equals the difference between viscous heating Q_+ and radiative cooling Q_-. Model applicability depends on accretion rate: ADAFs prevail at low \dot{M} < 0.01 \dot{M}_{\rm Edd}, yielding dim, hard spectra, while slim disks emerge at higher \dot{M} > 0.1 \dot{M}_{\rm Edd}, bridging to brighter, softer emissions. ADAFs, in particular, explain the underluminous spectra of low-luminosity active galactic nuclei and the Galactic center black hole Sgr A*.

Dynamical Instabilities and Magnetism

Magnetorotational Instability

The magnetorotational instability (MRI) arises in differentially rotating, weakly magnetized plasmas, such as those in accretion disks, where a weak magnetic field couples with the shear to drive axisymmetric perturbations that grow exponentially. This mechanism destabilizes the flow by magnetic tension acting like a spring, pulling radially displaced fluid elements in the azimuthal direction while the differential rotation winds up field lines, leading to instability on the dynamical timescale of approximately the orbital frequency \Omega^{-1}.^[19]^[20] The instability criterion, derived from ideal magnetohydrodynamics (MHD), requires that the angular velocity decreases outward (d\Omega / d \ln r < 0) in the presence of a weak magnetic field; this condition is satisfied in nearly all astrophysical accretion disks due to their approximate Keplerian rotation profiles where \Omega \propto r^{-3/2}.^[19] In linear theory, the dispersion relation for perturbations with wavenumber k is given by

\omega^4 - \omega^2 \left[ \kappa^2 + 2 (k \cdot v_A)^2 \right] + (k \cdot v_A)^2 \left[ (k \cdot v_A)^2 + \frac{d \Omega^2}{d \ln r} \right] = 0,

where \omega is the frequency, \kappa is the epicyclic frequency, and v_A is the Alfvén speed; instability occurs for wavenumbers satisfying k v_A / \Omega < \sqrt{ - d \ln \Omega / d \ln r }, with the maximum growth rate approximately (3/4) \Omega for Keplerian shear.^[19]^[20] In the nonlinear regime, the MRI saturates into a turbulent state characterized by correlated radial-azimuthal fluctuations that transport angular momentum outward, quantified by the stress \langle t_{r\phi} \rangle \approx 0.01 \rho v_A^2. This turbulence enables an effective viscosity parameter \alpha \sim 0.01 - 0.1 in the Shakura-Sunyaev prescription, consistent with observed accretion rates in protoplanetary, stellar, and active galactic nucleus disks.^[21]^[20] Numerical simulations confirm that the MRI is essential for sustaining accretion by generating the required turbulent transport, as laminar disks would accrete far too slowly to match observations.^[21]

Magnetic Fields and Outflows

Large-scale magnetic fields play a crucial role in the dynamics of accretion disks, particularly in enabling the efficient extraction and ejection of energy and angular momentum through outflows. These fields are typically weak initially but can be amplified through dynamo processes driven by the magnetorotational instability (MRI), which acts as a seed mechanism for turbulence that sustains field growth. Amplification occurs via dynamo action, where turbulent motions and differential rotation in the disk wind up and strengthen the fields until they reach equipartition, characterized by a plasma β ≈ 1, where β = P_gas / P_mag balances gas and magnetic pressures.^[22] This process allows magnetic fields to exert significant influence on disk evolution, transporting angular momentum outward and facilitating the launching of powerful outflows. Two primary mechanisms explain how these amplified fields drive relativistic jets from accretion disks around black holes. The Blandford-Payne mechanism involves magneto-centrifugal winds launched directly from the disk surface, where rotating magnetic field lines anchored to the disk accelerate plasma along open field lines inclined at more than 30° to the disk plane, extracting rotational energy from the Keplerian motion. In contrast, the Blandford-Znajek process extracts energy from the black hole's spin by threading the ergosphere with twisted magnetic fields, converting rotational energy into electromagnetic energy via frame-dragging effects. Both mechanisms rely on strong poloidal fields that become toroidal due to rotation, enabling efficient energy transfer. Collimation of these jets into narrow, relativistic beams is primarily achieved through the hoop stress provided by the toroidal magnetic field component, which acts like a magnetic tension confining the outflow against expansion, similar to the tension in a wire loop.^[23] This self-confinement allows jets to propagate over vast distances while maintaining their structure. Observationally, such jets are evident in active galactic nuclei (AGN) and Galactic microquasars, where they reach speeds up to ~0.9c and carry powers on the order of ~0.1 times the accretion luminosity L_acc, powering phenomena like radio lobes and synchrotron emission. The power in these magnetically dominated jets is often carried predominantly by the Poynting flux, with the jet power approximated as

P_\mathrm{jet} \approx \frac{B_\phi^2}{8\pi} \, 4\pi r^2 v,

where B_φ is the toroidal magnetic field strength, r is the cylindrical radius, and v is the outflow speed; this can achieve efficiencies up to 100% of the accreted energy in ideal cases.^[24] Recent advances in mean-field dynamo models have improved predictions of field strength evolution, incorporating nonlinear saturation effects and shear-driven α-Ω dynamos to better describe how large-scale fields persist and interact with disk turbulence over long timescales.^[25]

Observational Evidence

Manifestations in Stellar Systems

Accretion disks around compact objects and young stars in stellar systems produce observable signatures across X-ray, ultraviolet, optical, and infrared wavelengths, revealing their structure, dynamics, and evolution. These manifestations include thermal emission, variability patterns, and molecular line profiles that trace disk properties such as temperature, density, and rotation. In systems with stellar-mass black holes or neutron stars, such as X-ray binaries, the disks exhibit state-dependent spectra and rapid fluctuations linked to accretion instabilities. Cataclysmic variables, involving white dwarf accretors, show outburst cycles driven by thermal-viscous instabilities in the disk. Protostellar disks surrounding T Tauri stars display infrared excesses and kinematic tracers indicative of ongoing mass inflow and orbital motion. In X-ray binaries, accretion disks transition between high/soft and low/hard spectral states, reflecting changes in accretion rate and disk geometry. During the high/soft state, the spectrum is dominated by thermal blackbody emission from the inner disk, where the disk extends close to the innermost stable circular orbit, as observed in the black hole binary Cygnus X-1.^[26] In contrast, the hard state features a truncated disk with non-thermal Comptonized emission prevailing, reducing the thermal component's contribution.^[27] Quasi-periodic oscillations (QPOs) in these systems, often detected in the soft state of sources like Cygnus X-1, arise from instabilities in the inner accretion flow, such as relativistic precession or magnetorotational turbulence, with frequencies ranging from millihertz to kilohertz.^[28] Variability timescales as short as milliseconds probe the inner disk regions, reflecting dynamical processes near the compact object.^[29] Spectral analysis commonly employs the diskbb model in XSPEC to fit the multitemperature blackbody emission, providing estimates of disk temperature and inner radius that align with thin disk theory predictions.^[30] Cataclysmic variables, binary systems with white dwarf accretors, exhibit accretion disks that undergo dramatic brightness outbursts, particularly in dwarf novae subtypes. These outbursts, recurring on timescales of weeks to months, are explained by the disk instability model (DIM), where thermal-viscous instabilities cause rapid heating and enhanced accretion when the disk surface density exceeds a critical value.^[31] During quiescence, the disk cools and accumulates material until instability triggers a sudden rise in viscosity, leading to a factor of 100 or more increase in luminosity.^[32] The emission peaks in the ultraviolet and optical bands, arising from the heated disk atmosphere, with UV delays relative to optical rises supporting the inward propagation of heating fronts in the DIM framework.^[32] Protostellar accretion disks around T Tauri stars, low-mass pre-main-sequence objects, are identified through infrared excesses beyond the stellar photosphere, resulting from dust grains heated by stellar radiation and viscous dissipation.^[33] These excesses, prominent from near- to mid-infrared wavelengths, indicate flared disk geometries that intercept and reprocess stellar light. Molecular line observations, such as CO rotational transitions, trace the disk's rotation via Doppler-shifted profiles, revealing Keplerian velocity fields and inclinations out to radii of hundreds of AU.^[34] Typical disk lifetimes span 1–10 million years, during which mass accretion rates onto the star average around \dot{M} \sim 10^{-8} \, M_\odot \, \mathrm{yr}^{-1}, sustaining the observed excesses and line emission.^[35]^[36] Recent James Webb Space Telescope (JWST) observations as of 2025 have provided high-resolution mid-infrared spectra and images, revealing detailed disk chemistry, nested wind structures, and extended geometries that enhance models of accretion dynamics and planet formation.^[37]^[38]

Manifestations in Supermassive Black Holes

In active galactic nuclei (AGN) and quasars, accretion disks around supermassive black holes manifest through prominent spectral features driven by the disk's thermal and reprocessed emission. The intense ultraviolet (UV) radiation from the hot inner disk photoionizes surrounding gas in the broad-line region, producing broad emission lines in optical and UV spectra with widths of thousands of km/s, characteristic of high-velocity motions near the black hole. These lines, such as Hβ and Lyα, serve as key diagnostics for the disk's ionizing flux and geometry. A hallmark of quasar spectra is the "big blue bump," a broad excess in the UV continuum attributed to thermal blackbody emission from the disk's surface, peaking around 1000–2000 Å and contributing significantly to the bolometric luminosity. In X-rays, the disk reflects coronal emission, producing a Compton-reflected spectrum with a prominent iron Kα fluorescence line at 6.4 keV, whose relativistic broadening and redshift reveal the disk's inner extent down to a few gravitational radii from the black hole. Bright quasars typically operate near the Eddington accretion rate, with luminosities L / L_Edd ≈ 1 and total output reaching ~10^{46} erg/s, powering their immense brightness across cosmic distances. In low-luminosity systems like the Galactic center black hole Sgr A*, the accretion disk is dim and sub-Eddington (L / L_Edd ≪ 1), modeled as an advection-dominated accretion flow (ADAF) that is hot, geometrically thick, optically thin, and inefficient at radiating its gravitational energy. Observations of the Fermi bubbles—giant gamma-ray structures extending from the Galactic center—provide evidence for past energetic outflows possibly driven by disk winds or jets from enhanced accretion episodes around Sgr A*. The Event Horizon Telescope (EHT) has directly imaged these disks: the 2019 M87* observation revealed a bright ring of emission encircling the black hole shadow, interpreted as synchrotron radiation from a hot plasma in the inner disk, with the jet base closely connected to disk magnetic fields; subsequent 2024-2025 analyses confirmed a persistent shadow and provided insights into jet-disk coupling.^[39] Similarly, 2022 EHT images of Sgr A* show a comparable ring-like structure from orbiting disk material, and 2024 polarized observations revealed strong, spiraling magnetic fields at the edge of the accretion region, confirming the presence of a compact accretion flow with organized magnetism despite its low luminosity.^[40] Temporal variability in AGN light curves, occurring on timescales of hours to days, constrains the emitting region's size to the inner disk, roughly scaling with the black hole mass and accretion rate. Reverberation mapping, which measures time delays between UV continuum variations and broad-line responses, estimates disk sizes as R ≈ (G M / c²) (Ṁ / Ṁ_Edd)^{1/2}, linking the ionized broad-line region to the disk's illumination. Additionally, linear polarization in the UV and optical spectra of AGN arises from electron scattering in the disk atmosphere, providing insights into the disk's inclination and structure.

Advanced Developments

Numerical Simulations

Numerical simulations play a crucial role in elucidating the dynamics of accretion disks, capturing nonlinear phenomena such as turbulence and outflows that analytic models cannot fully resolve. These computations typically employ magnetohydrodynamic (MHD) frameworks to incorporate magnetic fields, which are essential for angular momentum transport. Local simulations using the shearing box approximation examine small-scale physics within a co-rotating frame of the disk, particularly the magnetorotational instability (MRI), while global simulations address the full radial and vertical structure, including relativistic effects near black holes. A key method for studying the MRI involves local shearing box simulations, which isolate a rectangular patch of the disk to focus on differential rotation and magnetic instabilities without curvature effects. These simulations reveal how weak seed magnetic fields amplify into turbulent states that enable accretion. For global relativistic disks around black holes, general relativistic MHD (GRMHD) codes such as HARM and Athena++ are widely used; HARM solves the equations in a 3+1 formalism for stationary spacetimes, while Athena++ employs advanced Riemann solvers and constrained transport to maintain numerical stability. Both codes evolve the MHD equations in conservative form:

\frac{\partial U}{\partial t} + \nabla \cdot \mathbf{F} = \mathbf{S}

where U represents the vector of conserved variables (including density, momenta, energy, and magnetic field components), \mathbf{F} the fluxes, and \mathbf{S} the source terms accounting for gravitational forces and relativistic effects.^[41]^[42] Simulations have yielded several pivotal insights. The MRI saturates at an effective Shakura-Sunyaev viscosity parameter \alpha \approx 0.01, quantifying the turbulent stress that drives accretion at rates consistent with observations of low-luminosity systems. Disk winds emerge ubiquitously from MRI turbulence, expelling material along magnetic field lines and contributing substantially to mass and angular momentum loss. In three-dimensional global MHD simulations, powerful jets are launched from the disk corona via magneto-centrifugal acceleration, demonstrating the viability of the Blandford-Payne mechanism in realistic geometries.^[43]^[44]^[45] Despite these advances, challenges remain in accurately resolving the multi-scale nature of disk turbulence. Insufficient grid resolution can suppress small-scale MRI modes, leading to underestimated transport rates, while numerical viscosity from discretization schemes introduces artificial diffusion that mimics physical effects. High-resolution runs, often requiring adaptive mesh refinement, are essential to mitigate these artifacts.^[46] Post-2018 developments have leveraged GPU acceleration to overcome computational limits, with codes like cuHARM enabling simulations of warped or long-term disk evolution at unprecedented scales. Additionally, integrating radiation transport into GRMHD frameworks has facilitated realistic modeling of slim disks, where photon trapping and advection dominate, revealing how radiative cooling influences disk thickness and outflow properties.^[47] Notable applications include GRMHD model comparisons with Event Horizon Telescope (EHT) images of M87* in the 2020s, where magnetically arrested disk configurations align with observed polarization patterns near the event horizon, constraining black hole spin and plasma properties. In hot accretion flows, simulations indicate outflow mass-loading rates \dot{M}_{\rm out} / \dot{M}_{\rm acc} \approx 1, implying that outflows remove comparable mass to what accretes, explaining the radial decline in inflow rates.^[48]^[49]

Recent Observations and Open Questions

Recent observations from the James Webb Space Telescope (JWST), launched in 2021 and operational since 2022, have provided unprecedented mid-infrared imaging of active galactic nuclei (AGN) tori and associated accretion disks, revealing detailed structures in obscured quasars at high redshifts. For instance, JWST's Mid-Infrared Instrument (MIRI) has detected hot dusty tori surrounding supermassive black holes in quasars at z > 7.^[50] In 2025, JWST MIRI-MRS spectroscopy of the highest-redshift luminous quasars has measured rest-frame optical/IR properties, advancing models of early universe black hole growth.^[51] Additionally, JWST observations have imposed limits on slim disk emission in high-redshift objects, supporting super-Eddington accretion scenarios where geometrically thick, radiation-pressure-dominated flows explain the lack of strong mid-IR excesses in luminous quasars.^[52] JWST's Near-Infrared Spectrograph (NIRSpec) has further illuminated disk dynamics through high-resolution spectra of quasars. These 2023–2024 spectra, with resolutions up to R ≈ 2700, trace metal enrichment in intervening absorbers along quasar sightlines, linking feedback processes to galaxy evolution at cosmic dawn.^[53] The Event Horizon Telescope (EHT) has advanced imaging of the Sagittarius A* (Sgr A*) accretion disk with polarized light observations published in 2024, revealing strong, ordered magnetic fields spiraling around the event horizon. The 2022 polarized images show a prominent spiral electric vector polarization angle pattern in the emission ring, consistent with twisted magnetic fields threading the disk and supporting general relativistic magnetohydrodynamic models.^[54] Follow-up EHT data analyzed through 2024 have detected variability patterns consistent with the magnetorotational instability (MRI), where turbulent magnetic fields amplify disk accretion and produce observed flux variations on hourly timescales.^[55] Despite these advances, several open questions persist in accretion disk physics. The origin of quasi-periodic oscillations (QPOs) remains debated, with models invoking relativistic precession of orbital hotspots in the inner disk competing against parametric resonance mechanisms where epicyclic frequencies couple to drive instabilities, as neither fully explains the observed frequency ratios in black hole systems. The efficiency of disk-jet coupling is another unresolved issue, with estimates suggesting only 5–20% of the disk's bolometric luminosity converts to jet power, but the exact role of black hole spin and magnetic field geometry in this process lacks a unified framework across spectral states.^[56] In ultraluminous X-ray sources (ULXs), the limits of super-Eddington accretion are unclear, as these systems exceed the Eddington luminosity by factors of 10–100 without evident pulsations in many cases, raising questions about photon trapping, wind geometries, and whether intermediate-mass black holes or beamed emission resolve the discrepancy. Finally, the near-universal value of the Shakura-Sunyaev viscosity parameter α ≈ 0.01 across diverse accreting systems—from cataclysmic variables to AGN—lacks a fundamental theoretical explanation, as it empirically fits observations but varies slightly without a clear physical origin tied to turbulence or magnetic processes.^[43]

Excretion Disks

Excretion disks, also referred to in some contexts as decretion disks or viscous spreading disks, describe astrophysical structures where angular momentum is transported radially outward, leading to disk expansion, in contrast to the primary inward mass flow in standard accretion disks.^[57] This outward angular momentum transport can occur through mechanisms such as stellar torque or injection of material with excess specific angular momentum, causing the inner material to lose angular momentum and spiral inward while the disk overall expands. For specific cases of net outward mass transport from the central body, see the decretion disks subsection below. Key mechanisms driving such disk evolution include viscous spreading, where internal friction diffuses angular momentum outward, resulting in disk growth over time. In protoplanetary disks, this viscous evolution leads to an increase in disk size as material migrates to larger radii, while net mass accretes onto the central star.^[58] The structure of these disks is similar to accretion disks in terms of rotational support and thin geometry, but features regions with positive radial velocity v_r > 0 in the outer parts, directing some mass flow away from the center. For self-similar solutions in viscous evolution scenarios, the surface density \Sigma decreases with time as \Sigma \propto t^{-3/2} (for power-law index \gamma = 1), reflecting the disk's expansion while conserving total mass.^[58]^[59] The dynamics are governed by a modified form of the viscous mass transport equation, where the mass flux \dot{M}(r) for outward flow is obtained by flipping the sign of the standard accretion expression:

\dot{M}(r) = 3\pi \nu \Sigma \frac{d \ln (\nu \Sigma r^{1/2}) }{d \ln r}

Here, \nu is the kinematic viscosity, \Sigma is the surface density, and r is the radial distance; positive \dot{M} denotes outward transport, analogous to the inward accretion case but with reversed angular momentum gradient. Representative examples include the viscous evolution of debris disks around main-sequence stars, such as \beta Pictoris, where low-density gas diffuses outward under viscous forces, shaping the disk's extended structure. In planet formation contexts, rings within protoplanetary disks experience viscous spreading driven by dust-induced instabilities, facilitating outward material transport and influencing ring widths.^[60]^[61]

Decretion Disks

Decretion disks are circumstellar structures formed by the ejection of material from a central star, resulting in the outward radial transport of mass and angular momentum, in direct contrast to the inward flow characteristic of accretion disks.^[62] These disks typically develop around rapidly rotating stars where equatorial material gains sufficient angular momentum to orbit in a quasi-Keplerian configuration, driven by mechanisms such as stellar rotation or torque.^[63] In classical Be stars—rapidly rotating B-type main-sequence stars with equatorial emission lines—decretion disks form a defining feature, often extending to tens of stellar radii and exhibiting gaseous, Keplerian rotation profiles.^[63] The material originates from the stellar photosphere, particularly the equator, where rotation velocities approach 70–80% of the Keplerian break-up speed, leading to mass injection rates that sustain disk buildup.^[64] Recent models highlight boundary layer effects at the star-disk interface, where a geometrically thick inner disk (aspect ratio H/R \approx 0.1) facilitates steady-state solutions even below critical rotation, enabling torque that slows the star while supporting observed decretion rates of approximately $10^{-9} to $10^{-10} M_\odot yr^{-1}. As of 2025, high-resolution simulations have further explored decretion disk formation in binary systems, confirming the role of boundary layers in smooth transitions from stellar rotation to disk flow.^[64]^[65]^[66] The viscous decretion disk (VDD) model provides a foundational framework for understanding their evolution, incorporating \alpha-viscosity prescriptions to describe angular momentum transport via turbulent diffusion.^[63] In this model, disks are geometrically thin and optically thick within about 20 stellar radii, transitioning to transonic outflows at sonic points often beyond 100 stellar radii, influenced by radiative line forces and viscosity parameters (\alpha \approx 0.1–0.6).^[63] During phases of disk growth, higher \alpha values promote rapid expansion, while dissipation phases feature slower viscous spreading and eventual truncation by stellar winds or binary interactions.^[63] Beyond Be stars, decretion-accretion diffusional disks appear in broader astrophysical contexts, such as evolving binary systems or young stellar objects, where nonstationary diffusion governs the interplay between inward accretion and outward decretion, parameterized by turbulent velocities and mean free paths.^[67] These hybrid structures highlight the diffusive nature of disk evolution, with analytical solutions for mass and angular momentum distributions underscoring their role in stellar mass loss and angular momentum regulation across diverse environments.^[67]

References

[1]
Glossary term: Accretion Disk
Accretion disks are found in a variety of astrophysical situations such as around supermassive black holes, stellar remnants, gamma ray bursts, or protostars.
[2]
Black Hole Anatomy - NASA Science
The main light source from a black hole is a structure called an accretion disk. Black holes grow by consuming matter, a process scientists call accretion.Anatomy of a Black Hole · Event Horizon · Accretion Disk · Doppler Beaming
[3]
Luminous disks: How black holes light up their surroundings
But once a particle of infalling matter plunges into an accretion disk – and possibly earlier – the particle's motion is disturbed. Due to frequent collisions ...
[4]
Accretion discs - Scholarpedia
Aug 12, 2014 · Accretion discs are flattened astronomical objects made of rapidly rotating gas which slowly spirals onto a central gravitating body.
[5]
https://ui.adsabs.harvard.edu/abs/1964ApJ...140..796S/abstract
[6]
Galactic Nuclei as Collapsed Old Quasars - Nature
LYNDEN-BELL, D. Galactic Nuclei as Collapsed Old Quasars. Nature 223, 690–694 (1969). https://doi.org/10.1038/223690a0. Download citation. Received: 08 July ...
[7]
Foundations of Black Hole Accretion Disk Theory | Living Reviews in ...
Jan 14, 2013 · This leads us to the four primary accretion disk models that we review: Polish doughnuts (thick disks), Shakura-Sunyaev (thin) disks, slim disks ...
[8]
https://ui.adsabs.harvard.edu/abs/1973A&A....24..337S/abstract
[9]
From Molecular Cloud Cores to Protostars and Protoplanetary Disks
Abstract. This review describes realistic evolution of magnetic field and rotation of the protostars, dynamics of outflows and jets, and the formation and ...
[10]
https://ui.adsabs.harvard.edu/abs/1981ARA&A..19..137P/abstract
[11]
https://ui.adsabs.harvard.edu/abs/1975ApJ...195L..65B/abstract
[12]
[PDF] accretion disks - MPA Garching
Abstract. In this lecture the basic theory of accretion disks is reviewed, with emphasis on aspects relevant for X-ray binaries and Cataclysmic Vari-.
[13]
https://doi.org/10.1007/BF01228565
[14]
https://doi.org/10.12942/lrr-2013-1
[15]
https://doi.org/10.1086/152966
[16]
https://ui.adsabs.harvard.edu/abs/1994ApJ...428L..13N/abstract
[17]
https://ui.adsabs.harvard.edu/abs/1978A&A....63..221A/abstract
[18]
https://ui.adsabs.harvard.edu/abs/1988ApJ...332..646A/abstract
[19]
https://ui.adsabs.harvard.edu/abs/1991ApJ...376..214B/abstract
[20]
https://ui.adsabs.harvard.edu/abs/1998RvMP...70....1B/abstract
[21]
https://ui.adsabs.harvard.edu/abs/1995ApJ...440..742H/abstract
[22]
Magnetic field amplification via protostellar disc dynamos
It has been shown that the MRI can amplify and sustain magnetic fields in accretion discs, despite dissipative effects. ... Differential disc rotation causes this ...
[23]
Jet Collimation by Small-Scale Magnetic Fields - IOPscience
The magnetic "hoop stress" is balanced either by the gas pressure of the jet or by centrifugal force if the jet is spinning. Since the length scale of the ...
[24]
Poynting Jets from Accretion Disks - IOPscience
We find solutions with magnetically collimated Poynting jets in which there is a continuous outflow of energy, angular momentum, and toroidal magnetic flux ...Missing: formula | Show results with:formula
[25]
Advances in mean-field dynamo theory and applications to ...
Aug 7, 2018 · Ensemble averaging obeys the Reynolds rules and could describe non-axisymmetric mean fields, but the practical meaning of such averaging is ...
[26]
Modelling hard and soft states of Cygnus X-1 with propagating mass ...
Aug 19, 2017 · We present a timing analysis of three Rossi X-ray Timing Explorer observations of the black hole binary Cygnus X-1 with the propagating mass ...Missing: disks | Show results with:disks
[27]
The relation between radio and X-ray luminosity of black hole binaries
If the accretion disk reaches inward towards the last stable orbit, the spectrum is dominated by a strong thermal component, the state is known as the high/soft ...
[28]
Detection of a Low-frequency Quasi-periodic Oscillation in the Soft ...
We compare QPOs at similar frequencies that have been previously detected in other persistent high-mass X-ray binaries in the soft state; we speculate that such ...
[29]
Accretion disc variability in the hard state of black hole X-ray binaries
We use this technique to show that the disc blackbody component varies on all time-scales, but varies more, relative to the power law, on longer time-scales.
[30]
disk: accretion disk, black body - HEASARC
The spectrum from an accretion disk, where the opacities are dominated by free-free absorption, i.e., the so-called blackbody disk model. Not correct for a disk ...Missing: fitting variability timescales
[31]
[PDF] A review of the disc instability model for dwarf novae, soft X-ray ...
Oct 4, 2019 · Dwarf novae (DNe) are cataclysmic variables (CV) that undergo regular outbursts (see Warner, 2003, for an ency- clopaedic review), in contrast ...
[32]
Testing the disk instability model of cataclysmic variables
The disk instability model (DIM) attributes the outbursts of dwarf novae ... optical magnitudes do for CVs, which in outburst emit primarily UV radiation.
[33]
Optical and Near-infrared Excesses are Correlated in T Tauri Stars
Apr 4, 2022 · Accretion produces a UV and optical excess, while re-radiated emission at the inner edge of the dust component of the disk produces a near- ...
[34]
Mapping Protoplanetary Disk Vertical Structure with CO ...
May 8, 2023 · ... CO rotational lines are tracing the cold gas in the outer disk. Thus, as the line emission is originating from similar heights, we expect ...
[35]
LIFETIMES AND ACCRETION RATES OF PROTOPLANETARY DISKS
Mar 16, 2016 · lifetime of protoplanetary disks is in the range of 1–10 Myr ... systems have either no disk or disks with very short lifetimes. (less than 5. 10 ...
[36]
The Mass Accretion Rates of Intermediate-Mass T Tauri Stars
We find that the average mass accretion rate for IMTTSs is ∼3 × 10-8 M⊙ yr-1, a factor of ∼5 higher than that for their low-mass counterparts. The new ...
[37]
Not found - ADS
**Summary:**
[38]
An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-mesh Constrained Transport
- **Confirmation**: Yes, Athena++ is a GRMHD (General Relativistic Magnetohydrodynamics) code used for black hole accretion disks.
[39]
The viscosity parameter α and the properties of accretion disc ...
We analyse a set of observations of dwarf-nova and AM CVn star outbursts and from the measured decay times determine the hot-disc viscosity parameter α h.Missing: structure | Show results with:structure
[40]
Shearing box simulations of accretion disk winds
The turbulence driven by MRI leads to internal stresses that redistribute angular momentum as would a high viscosity, thus driving accretion. MRI in disks have ...
[41]
Three-dimensional MHD Simulations of Jets from Accretion Disks
Aug 18, 2005 · We report the results of 3-dimensional magnetohydrodynamic (MHD) simulations of a jet formation by the interaction between an accretion disk and ...Missing: launching | Show results with:launching
[42]
Measuring the numerical viscosity in simulations of protoplanetary ...
This information can be used to make informed decisions on how to measure and minimize numerical diffusion in hydrodynamics simulations of accretion disks, and ...
[43]
Radiation Transport Two-temperature GRMHD Simulations of ...
Feb 22, 2023 · We present the first radiative two-temperature general relativistic magnetohydrodynamic (GRMHD) simulation of a strongly tilted (65°) accretion disk around an ...
[44]
On the Comparison of AGN with GRMHD Simulations: II. M87
### Summary of Key Points on EHT Comparisons with GRMHD Models for M87*
[45]
[1401.0586] Hot Accretion Flows Around Black Holes - arXiv
Jan 3, 2014 · Hot accretion flows, the topic of this review, are virially hot and optically thin. They occur at lower mass accretion rates, and are described ...
[46]
(PDF) Super-Eddington Accretion in Quasars - ResearchGate
Feb 10, 2025 · PDF | This review provides an observational perspective on the fundamental properties of super-Eddington accretion onto supermassive black ...
[47]
Metal enrichment and evolution in four z > 6.5 quasar sightlines ...
In this paper, we present high signal-to-noise ratio R ≃ 2700 JWST/NIRSpec spectra of four well-studied quasars at 6.5 < z < 7.5. We focus our analysis ...Missing: disk | Show results with:disk
[48]
First Sagittarius A* Event Horizon Telescope Results. VII ...
The linear polarization images demonstrate that the emission ring is highly polarized, exhibiting a prominent spiral electric vector polarization angle pattern ...
[49]
Astronomers unveil strong magnetic fields spiraling at the edge of ...
Mar 27, 2024 · The EHT has conducted several observations since 2017 and is scheduled to observe Sgr A* again in April 2024. Each year, the images improve as ...
[50]
The contribution of spin to jet-disk coupling in black holes
We estimate that 5–20% of the bolometric disk luminosity is converted to jet power. Most of the scatter to the optical-radio correlation is due to environment; ...
[51]
https://arxiv.org/abs/2511.02902
[52]
Observed sizes of planet-forming disks trace viscous spreading
In the viscous disk theory, conservation of angular momentum ensures that some parts of the disk move outward. As a result, disk sizes should grow with time.
[53]
Revisiting viscous transonic decretion disks of Be stars
We re-studied the viscous transonic decretion disk model of these stars. This model is driven by a radiative force due to an ensemble of optically thin lines ...
[54]
Self-similar spreading of a viscous, star-forming galactic disc. - ADS
The self-similar solution matches the galactic gas distribution indeed well. To represent the stellar profile would require a radially more extended gas disc ...Missing: Σ t^
[55]
self-consistent model for the evolution of the gas produced in the ...
To understand the distribution of gas in β Pictoris, we develop a numerical model of viscous diffusion of gas in the low-density regime expected in debris discs ...
[56]
Dust-driven viscous ring-instability in protoplanetary disks
From viscous disk theory it is known that a disk tends to increase its surface density in regions where the viscosity (i.e. turbulence) is low. Local maxima in ...
[57]
Astrophysical discs - DAMTP
Rapidly rotating early-type main-sequence stars, known as Be stars, are surrounded by an equatorial decretion disc that is expelled through the action of torque ...
[58]
Formation of Be star decretion discs through boundary layer effects
Mar 3, 2025 · We find one-dimensional steady state decretion disc solutions that smoothly transition from a rapidly rotating star that is in hydrostatic balance.
[59]
Models of astrophysical decretion-accretion diffusional disks
The formation of gaseous diffusional accretion-decretion disks is an important stage in the evolution of numerous astronomical objects. Matter is accreted.