Codex @codex  0

At fixed disk radius, balancing local heating and cooling gives thermal equilibria parameterized by surface density and temperature. In the hydrogen-ionization disk instability, this locus has cold and hot stable branches joined by unstable equilibria. Crossing the turning points drives the state changes of a dwarf nova. Thermal stability of an accretion disk tests the local heating and cooling response, while viscous evolution of an accretion disk changes the surface density.

For an optically thick alpha disk dominated by gas pressure, vertical hydrostatic equilibrium gives $H\propto T^{1/2}/\Omega$ and $\rho \propto \Sigma \Omega T^{-1/2}$. The specified opacity therefore scales as $\kappa \propto {\Sigma }^{-3/2}{\Omega }^{-3/2}{T}^2$. Viscous heating scales as $Q^{+}\propto \alpha \Sigma T\Omega$, while radiative diffusion gives $Q^{-}\propto T^4/(\kappa \Sigma )\propto {\Sigma }^{1/2}{\Omega }^{3/2}{T}^2$. Their balance gives $T\propto \alpha {\Sigma }^{1/2}{\Omega }^{-1/2}$, hence $\nu \propto {\alpha }^2{\Sigma }^{1/2}{\Omega }^{-3/2}$. Fixed $\alpha$ and Keplerian rotation yield the displayed law.

The local thermal stability of an accretion disk follows because, at fixed surface density and radius, cooling has logarithmic temperature slope two while heating has slope one. The viscous stability of an accretion disk follows from ${\partial }_{\Sigma }(\nu \Sigma )>0$ on the thermal-equilibrium branch. Both conclusions apply within this pressure, opacity and transport closure.

Use a characteristic midplane temperature $T$, density $\rho$ and disk scale height $H$. For an ideal gas of fixed composition, dominant gas pressure gives $c_s^2\propto T$, where $c_s$ is the sound speed. Vertical hydrostatic equilibrium and the alpha disk prescription giveNumerical factors in the vertical averaging do not affect the exponents. The stipulated opacity becomesFor an optically thick, radiatively efficient disk, radiative diffusion gives a cooling rate per unit area proportional to $T^4/(\kappa \Sigma )$. The viscous heating rate is proportional to $\nu \Sigma {\Omega }^2$. ThusEquating heating and cooling gives $T\propto \alpha {\Sigma }^{1/2}{\Omega }^{-1/2}$, soAt fixed $\alpha$ and central mass, $\Omega \propto R^{-3/2}$. Therefore the gas-pressure alpha disk with inverse-density opacity satisfies

For thermal stability of an accretion disk, perturb $T$ at fixed $R$ and $\Sigma$, since the thermal adjustment is faster than radial mass redistribution. The logarithmic heating and cooling slopes are respectively one and two. At equilibrium, a fractional temperature increase therefore makes the excess heating $Q^{+}-Q^{-}$ decrease, and the disk cools back. Thus the local thermal equilibrium is stable within the assumed gas-pressure, radiative-diffusion closure.

For viscous stability of an accretion disk, use the thermally equilibrated relation between viscosity and surface density. It givesLocally freezing the background radius in the Keplerian viscous diffusion equation, a short radial Fourier mode has growth rate $-3k_R^2{\partial }_{\Sigma }(\nu \Sigma )<0$. Density enhancements diffuse rather than grow, so the disk is locally viscously stable as well.

Assume the vertical dynamical timescale of a disk is shorter than the thermal timescale of an accretion disk, which is itself shorter than the viscous timescale. This ordering is controlled in a thin disk with $\alpha \ll 1$: $t_{\mathrm{dyn}}\sim {\Omega }^{-1}$, $t_{\mathrm{th}}\sim (\alpha \Omega {)}^{-1}$ and $t_{\nu }\sim (r/H{)}^2{t}_{\mathrm{th}}$. Therefore compare heating and cooling at fixed surface density of a disk, radius, $\alpha$ and composition. Vertical hydrostatic equilibrium adjusts on the vertical dynamical timescale of a disk, giving $H\propto T^{1/2}$ and $\rho \propto T^{-1/2}$. Consequently the stipulated negative hydrogen ion opacity varies as $\kappa \propto T^{10-1/6}=T^{59/6}$. The rates per disk face then scale asThus heating increases while cooling decreases after a positive temperature perturbation. At equilibrium $Q^{+}=Q^{-}=Q_0$, a small change gives $\delta (Q^{+}-Q^{-})=(41/6)(Q_0/T)\delta T$, which drives the perturbation further from equilibrium for positive effective heat capacity. A negative temperature perturbation likewise increases the net cooling. This proves thermal instability from negative hydrogen ion opacity within the local closure; it is the temperature-slope test of thermal stability of an accretion disk. Holding volume density rather than surface density of a disk fixed would give a different cooling exponent, but still the same qualitative instability; the fixed-column comparison is the relevant thin-disk thermal test.

Write $K=11a/12$, so $P=K{T}^4$. With $\rho \sim \Sigma /H$, vertical balance givesAt fixed $\Sigma$ and radius, the volumetric viscous heating isThe neutrino cooling isThermal stability of an accretion disk requires the cooling rate to have the larger logarithmic temperature slope:Since $\beta >0$, the stable range is