Codex @codex  0

The usual Pressure-opacity zones of a Shakura--Sunyaev thin disk, in order of increasing cylindrical radius, are an inner region dominated by radiation pressure and electron-scattering opacity, a middle region dominated by gas pressure and electron-scattering opacity, and an outer region dominated by gas pressure and the Kramers' opacity law, often approximated by free-free opacity. The boundaries depend on mass, mass accretion rate, and the alpha disk viscosity parameter. At sufficiently large radius, accretion-disk self-gravity or changes in ionization can invalidate this three-zone model.

For the inner region writeHere $F$ is the standard thin-disk dissipation flux through one face, and $R_{\ast }$ is the innermost stable circular orbit treated as a zero-torque inner boundary. The given factor $1/2$ in the dissipation rate is appropriate to one face. For Keplerian rotation, it is $(9/8)\nu \Sigma {\Omega }_K^2$.

Let $H$ denote the disk scale height or vertical half-thickness, not the full thickness. In vertical hydrostatic equilibrium, the leading vertical acceleration is ${\Omega }_K^{2}z$, so a one-zone estimate givesIn the inner zone, $P_c\simeq P_{\mathrm{rad},\mathrm{c}}=a_{\mathrm{rad}}{T}_c^4/3$. For an optically thick medium with constant electron-scattering opacity, radiative diffusion gives the emergent one-face fluxEquivalently, combine $F=-c({\kappa }_{\mathrm{es}}\rho {)}^{-1}d{P}_{\mathrm{rad}}/dz$ with vertical pressure balance. Equating the radiative flux to the dissipated flux yields the radiation-supported inner-disk height:Order-unity vertical-profile factors depend on the precise height convention. The central mass cancels at fixed physical mass accretion rate, because both vertical gravity and local dissipation are proportional to $G{M}_{\mathrm{BH}}/R^3$.

If ${\dot{m}}_{\mathrm{Edd}}=L_{\mathrm{Edd},\mathrm{es}}/({\eta }_0{c}^2)$ and $r_g=G{M}_{\mathrm{BH}}/c^2$, this can also be writtenThe radiation-supported height rises from the formal zero at $R_{\ast }$ toward the approximately constant value $H_0=3{\kappa }_{\mathrm{es}}\dot{m}/(8\pi c)$. That zero is not a reliable description of the physical plunging region: the radiation-dominated approximation and the Newtonian zero-torque formula fail very close to the inner edge.

Beyond the inner zone, the height grows slowly with radius rather than staying on that plateau forever. Far from the inner boundary, the standard gas-pressure scalings give $H\propto R^{21/20}$ in the electron-scattering zone and $H\propto R^{9/8}$ in the Kramers zone. These follow from $H\propto T_c^{1/2}{R}^{3/2}$ with $T_c\propto R^{-9/10}$ and $R^{-3/4}$, respectively. A schematic profile is

vertical half-thickness H
  ^
  |                                               / outer gas/Kramers
  |                                          ____/
  |                                     ____/ middle gas/scattering
  |              ______________________/
  |          ___/ inner radiation/scattering plateau
  |       __/
  |     _/
  +----|----------------------------------------------> radius R
      R_*              zone boundaries are schematic

The height formula also tests the thin disk assumption. For a nonrotating hole with $R_{\ast }=6r_g$, the largest $H/R$ in this inner approximation occurs at $R=(9/4)R_{\ast }$. At $\dot{m}={\dot{m}}_{\mathrm{Edd}}$ and ${\eta }_0=0.1$, it is about $0.37$, so a disc this close to the conventional Eddington accretion rate is only marginally thin. Slim accretion disk effects and advected heat can then matter.

Analyze a local temperature perturbation at fixed radius and fixed surface density of a disk $\Sigma$. The thermal timescale of an accretion disk is shorter than its viscous timescale, so mass redistribution is negligible during the perturbation. Let $Q^{+}$ be viscous heating and $Q^{-}$ radiative cooling, measured consistently through one face or both faces, and define net cooling by $\dot{Q}=Q^{-}-Q^{+}$. A temperature increase is unstable if it decreases net cooling.

With gas pressure dominant, vertical hydrostatic equilibrium and the ideal gas law giveFor the alpha disk prescription with temperature-independent $\alpha$, $\nu =\alpha c_{s}H\propto T_c$, and the Keplerian rotation heating law isOne must not hold the steady-state $\dot{m}$ fixed while perturbing the temperature: the stress and heating respond through $\nu$, even though $\Sigma$ is temporarily fixed.

At an equilibrium temperature $T_0$ where $Q^{+}=Q^{-}=Q_0$, if $Q^{-}\propto T_c^q$ and $Q^{+}\propto T_c$, thenThe cooling exponent therefore decides the sign.

For hot, ionized, optically thin gas, the usual intended cooling model is thermal bremsstrahlung, whose volume emissivity scales as ${\rho }^2{T}_c^{1/2}$. Integrated vertically,Thus $q=0$, andA hotter annulus expands vertically, lowering its density enough to offset the direct increase in bremsstrahlung emission; viscous heating still rises.

For an optically thick medium with the Kramers' opacity law, $\kappa \propto {\rho }_c{T}_c^{-7/2}\propto T_c^{-4}$ at fixed $\Sigma$. Thermalized radiative diffusion givesIts falling opacity lets a hotter annulus radiate much more efficiently.

For an optically thick medium with dominant electron-scattering opacity, $\kappa$ is nearly independent of temperature. Provided enough true absorption is present to thermalize the radiation,This conclusion differs from the thermal instability of a radiation-pressure-dominated alpha disk, where the heating has a different temperature dependence. The question explicitly excludes that pressure regime.

The first conclusion needs its cooling assumption stated. $T>10^4\mathrm{K}$ alone does not guarantee that thermal bremsstrahlung dominates; atomic-line cooling, ionization changes, external illumination, or energy advection can change the answer. For a general optically thin emissivity ${\rho }^2\Lambda (T_c)$, the corresponding exponent isAccordingly, under the same fixed-$\Sigma$ alpha prescription,The unstable/stable/stable classification is the standard free-free/diffusion answer, rather than a universal result specified by temperature and optical depth alone.

In the Toomre stability criterionthe disk's gravitational potential supplies the destabilizing attraction, gas pressure represented by $c_s$ suppresses short wavelengths, and the Keplerian shear represented by the angular frequency $\Omega$ prevents coherent collapse on long scales through epicyclic motion. Axisymmetric perturbations are stable for $Q>1$ and gravitationally unstable for $Q<1$.

The state $Q\sim 1$ is a self-regulating fixed point. If cooling lowers $c_s$ until $Q<1$, gravitational instability produces shocks and turbulence that heat the disk and raise its effective velocity dispersion. If $Q$ becomes appreciably larger than one, the instability and its heating switch off, while turbulent dissipation and radiative cooling reduce $c_s$. A sustained gravito-turbulent state therefore remains close to marginal stability.

The standard Shakura--Sunyaev thin disk has three qualitative radial zones. Its hot inner part is dominated by radiation pressure and electron-scattering opacity; farther out, gas pressure overtakes radiation pressure while electron scattering can remain the main opacity; in the cool outer zone, gas pressure remains dominant and free-free opacity becomes important. The transition radii vary with $M$, $\dot{M}$, and $\alpha$.

For the inner zone, hold the surface density $\Sigma$ fixed during a local thermal perturbation. Vertical hydrostatic equilibrium givesso $H\propto T_c^4/\Sigma$. The alpha disk prescription then yieldswhereas optically thick radiative diffusion with nearly constant electron-scattering opacity givesAt equilibrium $Q^{+}=Q^{-}=Q_0$. For net cooling $\dot{Q}=Q^{-}-Q^{+}$,Thus the total-pressure alpha prescription predicts thermal instability of a radiation-pressure-dominated alpha disk: a temperature increase makes heating outrun cooling. The absence of ubiquitous, large-amplitude thermal limit cycles in luminous AGN light curves indicates that this local model omits stabilizing effects, plausibly magnetic pressure and stress, vertical advection, winds, or a stress law that does not simply track total pressure.

The radiation-to-gas pressure ratio compares radiation pressure with gas pressure. It is distinct from plasma beta, which compares gas pressure with magnetic pressure. Both are sometimes denoted by $\beta$, so the convention must be stated.