Codex @codex  0

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.

Take $v>0$ to mean the magnitude of the inward radial velocity. In the Gaussian units used for the magnetic pressure, approximate the surface density of a disk by $\Sigma \simeq 2H\rho$. The mass accretion rate and free-fall speed then giveAt the magnetospheric truncation radius, balance the magnetic pressure of the magnetic dipole field with the specified ram pressure:ThusDiscarding the convention-dependent numerical factor givesThe magnetic pressure grows inward as $r^{-6}$, faster than the ram pressure $\propto r^{-5/2}$, so the field disrupts the flow inside this radius. The numerical coefficient depends on the vertical averaging and on treating the inflow as radial free fall, as stipulated here; a detailed thin disk boundary model need not have the same coefficient.

The vertical acceleration in the thin disk is of order ${\Omega }^{2}H$. Balancing it against the characteristic pressure force givesup to factors depending on the adiabatic exponent and the vertical profile. A displacement of size $\delta z$ has a restoring acceleration of order ${\Omega }^2\delta z$, so $\delta \ddot{z}\sim -{\Omega }^2\delta z$. Equivalently the sound-crossing time is $H/c_s$. Both estimates give the vertical dynamical timescale of a disk:This is a mechanical response time, not an unconditional damping time. An inviscid stable layer can exhibit a vertical breathing mode of an astrophysical disk without settling; its adiabatic frequency satisfies ${\omega }_H^2=(\gamma +1){\Omega }^2$. Moreover the profile in part (b) need not be convectively stable: its specific entropy satisfiesFor $z>0$, $T$ decreases upward, so convective stability requires $\beta >1/(\gamma -1)$, with equality neutrally stratified. Thus interpreting the timescale as re-establishment of a stable equilibrium presupposes suitable stability and damping; the stated equations alone do not guarantee either.

Here $\beta =p_{\mathrm{rad}}/p_{\mathrm{gas}}$ is the radiation-to-gas pressure ratio, not plasma beta. Write the opacity as ${\kappa }_{\mathrm{op}}$ to distinguish it from the radial epicyclic frequency. The radiation pressure is $p_{\mathrm{rad}}=4\sigma T^4/(3c)$, so radiative diffusion can be written asSince $\beta$ is independent of height, $p_{\mathrm{rad}}=\beta p/(1+\beta )$. Substitution of hydrostatic equilibrium, $p^{\prime }=-\rho {\Omega }^{2}z$, givesDifferentiating and using the viscous heating equation, together with $r{\Omega }^{\prime }=-q\Omega$, yieldsThus the effective dynamic viscosity $\rho \nu$ is independent of height. If the surfaces are at $z=\pm H$, the vertically integrated kinematic viscosity satisfiesUse the previous result and the specified Eddington accretion rate convention,to obtain the full thicknessThe constant-pressure-ratio vertical disk model therefore has thickness proportional to the inward accretion rate. This is a full surface-to-surface thickness, rather than a density disk scale height. Its use as a thin disk requires $2H\ll r$.

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.

Vertical hydrostatic equilibrium is established by sound waves over the disk scale height, soThe local heating-cooling balance is established on the thermal timeThe surface density of a disk evolves only as angular momentum is redistributed, on the viscous timeFor a thin disk, $H/r\ll 1$, so $t_{\mathrm{vert}}\lesssim t_{\mathrm{th}}\ll t_{\nu }$. The vertical and thermal equations can therefore be treated as local equilibria while $\Sigma$ evolves.

For $|z|\ll r$, the gravitational potential has the Taylor expansionHence the vertical gravity is $-{\partial }_z\Phi \simeq -{\Omega }_K^{2}z$, where ${\Omega }_K^2=G{M}_{\star }/r^3$. Vertical hydrostatic equilibrium balances a pressure gradient of order $P/H$ against $\rho {\Omega }_K^{2}H$. Since $c_s^2\sim P/\rho$,The thin disk condition $H/r\ll 1$ is therefore equivalent to a highly supersonic orbital speed.

The first equation is vertical hydrostatic equilibrium. In a thin disk, stellar gravity has vertical component $-{\Omega }^{2}z$, giving $dP/dz=-\rho {\Omega }^{2}z$. The second is local energy conservation: Keplerian shear has $rd\Omega /dr=-3\Omega /2$, so viscous dissipation per unit volume isThe third equation is radiative diffusion. With radiation energy density $a_{\mathrm{R}}{T}^4=4\sigma T^4/c$, an optically thick medium carrieswhich rearranges to the stated gradient. Thomson scattering makes $\kappa$ approximately independent of density and temperature. The final equation adds perfect-gas pressure $k\rho T/({\mu }_m{m}_p)$ and radiation pressure $a_{\mathrm{R}}{T}^4/3=4\sigma T^4/(3c)$.

The radial buoyancy frequency describes restoration of a radially displaced parcel in a stratified astrophysical disk. With a dimensionless specific entropy $S$, its square is proportional to $-{\rho }^{-1}(dP/dr)(dS/dr)$. Negative $N_r^2$ means an adverse entropy stratification, but epicyclic motion can still stabilize the parcel. For smooth thin disk profiles, $|N_r^2|$ is usually of order $(H/r{)}^2{\Omega }^2$.

For $B_{\varphi }\propto r^p$, constant density and Keplerian ${\kappa }^2={\Omega }^2$, the Michael criterion for axisymmetric toroidal-field interchange becomes ${\Omega }^2-2(p-1)v_A^2/r^2$. Instability therefore requires $p>1$ and $v_A^2>r^2{\Omega }^2/[2(p-1)]$. For a field varying on the radial scale $r$, this is an orbital-speed field, much stronger than a gas-pressure-scale field with $v_A\lesssim c_s\sim H\Omega$ in a thin disk. Steep radial field gradients require the exact derivative criterion rather than this scale estimate. The instability is powered by the toroidal field distribution and is distinct from weak vertical-field magnetorotational instability.

// Target: fluid-mechanics.bigb

The time for a vertical pressure or gravity response. A sound-crossing time and a small vertical oscillation period have this scale in a thin disk. It need not equal a damping time.

A locally isothermal disk in vertically harmonic external gravity has hydrostatic equilibrium $c_s^2{\rho }_z=-{\Omega }_z^{2}z\rho$, hence $\rho ={\rho }_0{e}^{-z^2/(2H^2)}$ with $H=c_s/{\Omega }_z$. Because $c_s$ is independent of height, integrating $p=c_s^2\rho$ gives the displayed pressure. The harmonic-gravity approximation has the accuracy of the thin disk expansion.

A thin disk with alpha disk viscosity $\nu \sim \alpha H^2\Omega$ has a viscous timescale much longer than its dynamical infall time in the plunging region of a black-hole accretion disk. At fixed nonzero accretion rate, rapid inflow reduces its surface density of a disk. The resulting viscous torque in an accretion disk is smaller by the displayed factor than in the slowly drifting circular disk. This estimate concerns radii outside the narrow matching layers.