Codex @codex  0

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.

The required pressure integral uses ${\int }_0^1(1-x^2{)}^{4}dx=128/315$. With $\Sigma =(32/35){\rho }_{0}H$ and $p_0/{\rho }_0={\Omega }^2{H}^2/8$,The vertically averaged alpha closure for a mixed-pressure layer therefore givesThe height-integrated dynamic viscosity is $2H\mu =8cH\lambda /[9\kappa (1+\lambda )]$. Equate it to $\bar{\nu }\Sigma$ to obtainThese are integrated equalities, not an imposed pointwise relation $\mu (z)\propto p(z)$.

Now keep $\alpha$, the opacity and molecular constants independent of radius, and use Keplerian rotation $\Omega \propto r^{-3/2}$. In the gas-dominated limit, $\lambda \ll 1$ but positive, the last equality gives $\lambda \propto H\Omega \Sigma$. The column-density relation gives $\lambda \Sigma \propto H^7{\Omega }^6$. Eliminating $\lambda$ yields $H^6\propto {\Sigma }^2{\Omega }^{-5}$, whenceThis is the gas-pressure branch of a mixed-pressure alpha disk. In the radiation-dominated limit, $\lambda /(1+\lambda )\to 1$, so $H\propto (\Omega \Sigma {)}^{-1}$ andFor viscous stability of an accretion disk, linearize the Keplerian viscous diffusion equation at fixed radius: a local mass density perturbation has diffusion coefficient $3{\partial }_{\Sigma }(\bar{\nu }\Sigma )$. The gas branch has $\bar{\nu }\Sigma \propto {\Sigma }^{5/3}$ and a positive derivative, so it smooths perturbations. The radiation branch has $\bar{\nu }\Sigma \propto {\Sigma }^{-1}$ and a negative derivative, giving radiation-pressure viscous instability. The gas-pressure branch is viscously stable; the radiation-pressure branch is viscously unstable in this closure. This concerns radial mass-transport stability on wavelengths where the vertically averaged thin-disc model applies.