Shakura–Sunyaev

Shakura–Sunyaev
Name	Nikolai Shakura and Rashid Sunyaev
Field	Astrophysics, Theoretical Physics
Known for	Alpha disk model, accretion theory, X-ray binaries
Notable works	1973 alpha viscosity paper
Awards	Lenin Prize, Crafoord Prize

Shakura–Sunyaev.

Shakura–Sunyaev is the eponymous description of a seminal accretion disk model introduced by Nikolai Shakura and Rashid Sunyaev, which provided a practical prescription for angular momentum transport in thin disks around compact objects such as black holes, neutron stars, and white dwarfs. The model established a phenomenological viscosity parameter that linked turbulent stress to local pressure, influencing theoretical work across studies of Soviet astrophysics, Roger Penrose, Subrahmanyan Chandrasekhar, John Wheeler, Kip Thorne, Martin Rees, and later numerical projects at institutions like the NASA Jet Propulsion Laboratory, European Space Agency, and observatories including Chandra X-ray Observatory and XMM-Newton. Its influence extends to interpretation of spectra, timing, and stability in systems studied by collaborations at Harvard-Smithsonian Center for Astrophysics, Max Planck Institute for Astrophysics, and Cambridge University.

Introduction

The Shakura–Sunyaev construct emerged in the early 1970s to resolve discrepancies in modeling luminous disks in systems such as Cygnus X-1, SS 433, and active nuclei like 3C 273 and Messier 87. Shakura and Sunyaev introduced a dimensionless parameter, alpha (α), to parametrize turbulent stresses without specifying microscopic mechanisms, enabling connection to observations from missions like Uhuru, Einstein Observatory, EXOSAT, and later ROSAT. The framework influenced theoretical programs at groups led by figures such as Evgeny Lifshitz, Lev Landau, Andrei Sakharov, and modern numerical efforts at Princeton University, Caltech, and MIT.

Shakura–Sunyaev α-disk Model

The α-disk model assumes a geometrically thin, optically thick disk in Keplerian rotation, balancing viscous heating, radiative cooling, and vertical hydrostatic equilibrium. In practice the model expresses the viscous shear stress as τ_{rφ} = α P, linking local pressure P to stress through α, a scaling inspired by turbulent analogies used by researchers associated with Ludwig Prandtl in fluid dynamics and later adapted by magnetohydrodynamic studies informed by work at Los Alamos National Laboratory and Lawrence Livermore National Laboratory. The model yields radial profiles for surface density, temperature, and emitted flux used to interpret continuum spectra from sources studied with Hubble Space Telescope, Very Large Telescope, Keck Observatory, and radio arrays like Very Large Array.

Mathematical Formulation and Assumptions

The formalism uses conservation equations for mass, angular momentum, and energy in a thin-disk approximation around a central mass M, often applied to environments exemplified by Sagittarius A*, stellar remnants cataloged by the International Astronomical Union, and quasars from surveys by Sloan Digital Sky Survey. Key assumptions include axisymmetry, local viscous torque representation via α, vertical hydrostatic balance, and radiative diffusion as the dominant cooling mechanism. Shakura and Sunyaev combined the Navier–Stokes approach with Keplerian angular velocity Ω_K from classical potentials used in works by Isaac Newton and refined in relativistic settings by Albert Einstein and Roy Kerr, later extended using general relativistic corrections developed in studies at University of Chicago and KIPAC. The α prescription reduces complex magnetohydrodynamic turbulence—later linked to the magnetorotational instability discovered by Steven Balbus and John Hawley—to a single parameter that can be constrained from variability and spectral fitting used in analyses by teams at Columbia University and University of California, Berkeley.

Applications to Accretion Systems

Practitioners apply the model to interpret emission from X-ray binaries such as GX 339-4, cataclysmic variables like SS Cygni, and active galactic nuclei exemplified by NGC 1068 and Seyfert galaxies. It underpins theoretical descriptions of spectral states, outburst cycles modeled in frameworks developed by groups at Max Planck Institute for Astronomy and INAF, and informs estimates of mass accretion rates for objects in catalogs compiled by the European Southern Observatory and National Radio Astronomy Observatory. The α parameter has been empirically inferred across systems observed by RXTE, NuSTAR, and multiwavelength campaigns organized by collaborations including ESO, NOAO, and Keck Observatory.

Observational Tests and Limitations

Comparisons with timing and spectral data from missions such as RXTE, XMM-Newton, and Chandra X-ray Observatory reveal successes and tensions: the model explains thermal continuum emission in high-accretion-rate states of sources like LMC X-3 but struggles with low-luminosity, radiatively inefficient flows studied in Sgr A* and low-hard states of black hole binaries like V404 Cygni. Limitations include the ad hoc nature of α, uncertainties in vertical structure influenced by magnetic fields explored at Princeton Plasma Physics Laboratory, and departures from local thermal equilibrium considered in work at Instituto de Astrofísica de Canarias and National Astronomical Observatory of Japan.

Extensions and Alternative Models

Extensions incorporate magnetohydrodynamic turbulence via the magnetorotational instability, leading to global simulations by groups at University of Illinois Urbana-Champaign, University of Oxford, and Yale University. Alternative models include advection-dominated accretion flows developed by Ramesh Narayan and Insituto de Astrofísica de Canarias collaborators, slim disks formulated in studies by Wuhan University and P. J. Armitage's groups, and radiatively inefficient accretion flows used to describe systems analyzed by Harvard University and Stanford University. Modern efforts combine general relativistic magnetohydrodynamics from teams at Max Planck Institute for Gravitational Physics and observational constraints from the Event Horizon Telescope collaboration to refine or replace the original prescription.

Category:Astrophysics