Codex @codex  0

The disk scale height measures the vertical spread of matter or pressure. For an isothermal, non-self-gravitating thin disk, $\rho (z)={\rho }_0\exp [-z^2/(2H^2)]$ and $H=c_s/{\Omega }_z$. Stronger vertical gravity confines the same gas more tightly; a higher sound speed makes it thicker. Here $H$ is a standard deviation, not a sharp surface; a polytropic surface uses a different convention.^[ref]

Assume a passive disk has local isothermal sound speed squared $c_s^2=aT$, disk scale height $H=c_s/{\Omega }_K$, and irradiated height $z_s=fH$ with constant $f$. Its irradiation grazing angle of a disk and Stefan–Boltzmann law give $q^{\prime }=C{q}^8/R^3$ for $q=z_s/R$, where $C=4\pi \sigma (GM/a{)}^4/(L_{\ast }{f}^8)$. Integrating gives $q^{-7}=R_0^2/R^2+B$, with $R_0^2=7C/2$. Where the first term dominates and the thin disk assumptions hold, $q\simeq (R/R_0{)}^{2/7}$, $H\propto R^{9/7}$ and $T\propto R^{-3/7}$. The law cannot be extrapolated into a geometrically thick region using the same grazing-angle approximation.

An accretion disk forms when infalling material has enough angular momentum to orbit before it can reach the central object. In a thin disc dominated by a central mass $M$, the orbital motion is approximately Keplerian rotation:Vertical hydrostatic equilibrium gives the disk scale height $H\sim c_s/{\Omega }_K$; a geometrically thin disc has $H/R\ll 1$. The surface density is $\Sigma =\int \rho dz$, and a steady inward mass flow rate is $\dot{M}=-2\pi R\Sigma v_R>0$.

Moving inward requires losing specific angular momentum. Internal stresses transport angular momentum outward; much of the mass can then accrete while an outer fraction spreads to larger radius. Molecular viscosity alone is generally inadequate for rapid astrophysical accretion. Turbulent and magnetic stresses, including those generated by the magnetorotational instability, provide effective transport. A common parametrization is $\nu =\alpha c_{s}H$ for an effective kinematic viscosity, with $\alpha$ a dimensionless stress parameter.

For a thin Keplerian accretion disk, the evolution equation isThis viscous evolution of an accretion disk occurs on the viscous timescale $t_{\nu }\sim R^2/\nu$, much longer than an orbital period in a thin disc with moderate effective stress.

The orbital specific energy is $-GM/(2R)$, so inward motion releases binding energy. In a steady, optically thick, radiatively efficient disc with zero torque at its inner edge, the emergent flux from each face isThe Stefan–Boltzmann law identifies the local effective temperature. Far from the inner boundary, $T_{\mathrm{eff}}\propto R^{-3/4}$, so the emitted spectrum combines hotter inner and cooler outer annuli. Integrating both faces gives $L_{\mathrm{disc}}=GM\dot{M}/(2R_{\mathrm{in}})$ under these assumptions. This is not a universal spectrum or efficiency: optical depth, inner torque, energy advection and outflows can change the radiation. The essential connection is outward angular-momentum transport, inward accretion, and release of binding energy.

A sufficiently self-gravitating disk can fragment directly, bypassing the slow assembly of a massive solid core. In the simplest thin gas-disk model the competition is between self-gravity, pressure and rotational shear. Axisymmetric perturbations havewhere $\kappa$ is the radial epicyclic frequency, not the opacity. In a nearly Keplerian disk, $\kappa \simeq \Omega$. This is Toomre's stability criterion; finite thickness and global/nonaxisymmetric modes modify the precise boundary. The fastest local wavelength is $\lambda =2c_s^2/(G\Sigma )$, of order the disk scale height times $2\pi$ when $Q_T\simeq 1$.

Instability alone does not ensure bound planets. Compression and shocks heat the gas, raising $Q_T$; a slowly cooling disk can maintain self-regulated spiral turbulence. Cooling-limited disk fragmentation requires energy loss on a sufficiently short dynamical timescale, conventionally expressed as $\beta =\Omega t_{\mathrm{cool}}$ below a model-dependent threshold. The often quoted $\beta \sim 3$ belongs to a particular local razor-thin calculation, not every real disk.^[ref] Cold massive outer disks are promising sites, but fragment contraction, tidal survival, further accretion and migration determine whether the result remains a planetary-mass body or grows into a brown dwarf.

For gravitational collapse of a settled particle layer, if “core collapse” instead denotes a solid-particle layer, its random-velocity dispersion replaces gas sound speed in an appropriate particle-disk stability condition. Settling and concentration can enable collapse into planetesimals; gas turbulence and particle stirring oppose it. This route makes solid seeds for later core accretion and is physically distinct from direct gas fragmentation. In both cases one must establish self-gravity against velocity support and tides, rather than infer collapse from high surface density alone.

A protoplanetary disk is a flattened, rotating reservoir of gas and solids around a young star. Its radial gravity is usually dominated by the host, making its rotation close to Keplerian rotation, ${\Omega }_K=(G{M}_{\star }/r^3{)}^{1/2}$. Gas pressure provides vertical support and a small radial correction to the orbital speed; it does not make the disk a nonrotating spherical atmosphere.

At height $|z|\ll r$, vertical hydrostatic equilibrium in the stellar gravitational field givesFor a locally vertically isothermal ideal gas, $P_g={\rho }_g{c}_s^2$ with $c_s^2=k_{B}T/(\mu m_H)$, where $\mu$ is mean molecular weight. Integration givesHere ${\Sigma }_g$ is the vertically integrated gas surface density of a disk, and $H$ is the isothermal disk scale height. The disk aspect ratio $H/r=c_s/v_K$ is small in a geometrically thin disk. For $T\propto r^{-q}$ it scales as $r^{(1-q)/2}$, so $q<1$ produces disk flaring. This illustrative local vertical model does not determine a unique radial surface-density law or temperature normalization from the information given.

The gas is predominantly molecular hydrogen and helium. A much smaller component consists of heavier elements, partitioned between vapor, solid refractory grains and condensed volatile ices. A gas-to-solid mass ratio of order $100$ is a common initial fiducial model, not a universal measurement applying at every disk radius or evolutionary stage. Solids can grow and settle while gas accretes or disperses, changing the local ratio substantially.^[ref]

Thermal structure organizes composition. In the hottest inner region, grains can sublimate and much of their material remains in the gas. Refractory silicate minerals and metallic grains survive at lower temperatures; farther out, water and then other sufficiently volatile species can condense when their local condensation conditions are met. The water snow line can raise the solid inventory available to planet formation. Its location depends on stellar irradiation, accretion heating, opacity and gas pressure, and evolves with time; it is not fixed merely by naming a star or adopting a solar-system distance.

Stellar radiation heats the exposed surface of a passive disk. Beneath it, shielding can leave a cooler, molecule-rich interior with a higher mass density. In the inner actively accreting disk, dissipation can instead make the optically thick midplane the hottest region. These are distinct heating regimes; a vertically isothermal profile is an analytic approximation rather than a claim that every real layer has the same temperature. Dust opacity controls how radiation is absorbed and re-emitted. Infrared and millimetre emission probe solid material, while spectral lines probe gas; neither observable alone supplies a model-free inventory of the whole disk.

Dust settling in an astrophysical disk makes larger, less tightly coupled particles more concentrated toward the midplane, while turbulence and mixing oppose that settling. For an aerodynamic stopping time $t_s$ short compared with ${\Omega }_K^{-1}$, terminal vertical balance gives $v_z\simeq -{\Omega }_K^2{t}_{s}z$. Small grains remain better coupled to the gas. Radial drift also changes the distribution of solids and their access to snow lines and embryos.

The disk evolves through accretion and angular momentum transport, potentially involving turbulent stresses, magnetorotational instability or magnetically driven winds where conditions permit. Protoplanetary disk dispersal eventually removes gas through accretion and mass loss, including photoevaporation. In standard evolutionary pictures this occurs on a scale of a few million years, with substantial system-to-system variation. That finite reservoir lifetime is central to whether an embryo acquires a massive gaseous envelope.^[ref]

The essential structure is a nearly Keplerian gas disk with vertical pressure support, radial and vertical thermal gradients, and an evolving solid population concentrated by settling, growth and drift. Its composition is set jointly by the original elemental inventory, local condensation conditions and transport.

Use the leading vertically averaged horizontal momentum equations of a thin disk. For an isothermal ideal gas, write $P=\int pdz=c_s^2\Sigma$, with constant isothermal sound speed $c_s$. The disk scale height is of order $H=c_s/{\Omega }_z$; this is exact for local isothermal vertical hydrostatic balance and remains the thickness scale for the estimate. Put $ϵ=H/r\ll 1$. The alpha disk prescription givesAt fixed radii of order $a$ inside the plunge and away from the marginal orbit, $|v_r|,v_{\varphi }$ and $r{\Omega }_z$ are fixed numerical multiples of $(GM/a{)}^{1/2}$. Hence the inflow time is dynamical, while the viscous timescale is longer by $1/(\alpha {ϵ}^2)$.

The Paczyński-Wiita ballistic plunge obeys $v_r{v}_r^{\prime }-v_{\varphi }^2/r=-{\Phi }^{\prime }$ and $v_r{j}_0^{\prime }=0$ exactly. For smoothly varying disk quantities on scale $r$, the pressure correction to radial acceleration has magnitude $c_s^2/r$, and the viscous correction is of order $\nu v/r^2$. Relative to the orbital acceleration scale $r{\Omega }_z^2$, these are respectivelyFinite-height corrections to the mid-plane gravitational force are also $O({ϵ}^2)$ by reflection symmetry. Thus the leading horizontal flow is ${\mathbf{u}}_h=\mathbf{v}+O({ϵ}^{2}v)$, with specific angular momentum differing from $j_0$ by a fraction $O(\alpha {ϵ}^2)$ over a dynamical transit, for $\alpha \lesssim 1$. Vertical structure and the smaller vertical velocity must be supplied by the thin-disk expansion; they are not specified by the planar ballistic expression.

This estimate is not uniform at the marginal orbit. In the ballistic solution, steady mass conservation givesConsequently $d\ln \Sigma /dr=[2(r-a){]}^{-1}+3[2(3a-r){]}^{-1}$. Writing $\delta =3a-r\ll a$, the pressure acceleration scales as $c_s^2/\delta$, while the net ballistic radial force scales as $(GM/a^2)(\delta /a{)}^2$. Their ratio is $O[{ϵ}^2(a/\delta {)}^3]$, evaluated with the characteristic thickness near $3a$. A transonic matching layer of width $\delta /a=O({ϵ}^{2/3})$ is therefore needed. The approximation concerns the bulk plunging flow outside that layer, not a uniform relative approximation arbitrarily close to $3a$.

For the torque comparison, the steady angular-momentum flux is $\mathcal{G}-\dot{M}j=\mathrm{constant}$, where $\mathcal{G}=-2\pi \nu \Sigma r^3{\Omega }^{\prime }$ is the viscous torque in an accretion disk. Inside the plunge, $\Omega =j_0/r^2$, soOutside, the approximately circular disk has $\mathcal{G}=\dot{M}[j(r)-j_0]+O(\alpha {ϵ}^2\dot{M}{j}_0)$, since the inner plunge torque is small. At $4a$,a fixed number independent of disk thickness. At $2a$, radial infall is dynamical and the inner torque is $O(\alpha {ϵ}^2\dot{M}{j}_0)$. Therefore the weak viscous torque in a thin-disk plunging region estimate givesThe numerical prefactor need not be small; the suppression is an asymptotic statement in the thin-disk parameter. Rapid infall and its low surface density explain why a nearly zero inner torque is a useful leading boundary condition for the outer accretion disk.

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.

A passive disk is heated predominantly by intercepted stellar radiation, rather than by its own viscous dissipation. A rough luminosity criterion, ignoring geometric factors, isThe intrinsic accretion disk luminosity is of order $GM\dot{M}/(2R_{\ast })$, with additional accretion energy potentially released at the stellar surface. The local test is more precise: absorbed irradiation must dominate the viscous flux at the radius being considered. For an ordinary steady Keplerian accretion disk, sufficiently far outside its inner edge,where $\alpha$ is the irradiation grazing angle of a disk. Reflection, very small incidence angles and shadowing can invalidate a simple total-luminosity comparison. The following calculation assumes an illuminated, absorbing surface.

At a point $(R,z_s)$ on the upper surface, the outward unit normal in the meridional plane is proportional to $(-z_s^{\prime },1)$, and a light ray from the point source is directed along $(R,z_s)$. Its flux through the upper face is the incident inverse-square flux times minus their scalar product:For a thin disk with small surface slope this becomesThe tangent slope and ray slope differ by $\alpha \simeq z_s^{\prime }-z_s/R$. Positive incidence requires disk flaring, $d(z_s/R)/dR>0$, and an unobstructed ray from the star. A negative value of the formula would mean that direct illumination of this face is absent, not that it receives a negative heating flux. The hierarchy $z_s\gg R_{\ast }$ makes the finite angular size of the star negligible here.

Let $a=\mathcal{R}/\mu$, so the ideal gas equation of state gives the isothermal sound speed $c_s^2=a{T}_c$. In the stellar vertical gravitational potential, hydrostatic equilibrium has vertical acceleration $-{\Omega }_K^{2}z$, where ${\Omega }_K^2=GM/R^3$. Integration of $c_s^2{\partial }_z\rho =-\rho {\Omega }_K^{2}z$ gives a Gaussian atmosphere with disk scale heightConsequently $q=z_s/R=fH/R$ satisfies $T_c=GM{q}^2/(a{f}^{2}R)$. Substitute this into the stipulated Stefan–Boltzmann law $F_d=\sigma T_c^4$ and the irradiation flux:Here $f$ specifies surface height in units of disk scale height and is unrelated to the wind-energy fraction in the preceding solution. This separable differential equation integrates exactly toIf the thin-disc solution is formally matched to $q(R_M)=1$, then $B=1-R_0^2/R_M^2$, givingIt follows that the flaring passive disk asymptotic isThis applies in the thin intermediate region $R\ll \min (R_0,R_M)$, in addition to the point-source and stellar-gravity assumptions. The corrections are of relative order $R^2/R_0^2+R^2/R_M^2$. In particular, the printed hierarchy $R\ll R_M$ must still be accompanied by the initial thin-disc condition $q\ll 1$: if $R_M\gg R_0$, radii larger than $R_0$ cannot satisfy that approximation. The matching at $q=1$ is an extrapolation used to set an integration constant; the derivation is not a valid treatment of a thick surface at $R_M$ itself.

Identify the stellar surface temperature with its effective radiating temperature, so that luminosity of a spherical blackbody gives $L_{\ast }=4\pi R_{\ast }^2\sigma T_{\ast }^4$. The given stellar surface scale height has $H_{\ast }^2=a{T}_{\ast }{R}_{\ast }^3/(GM)$, hence $GM/(a{T}_{\ast }{R}_{\ast })=(R_{\ast }/H_{\ast }{)}^2$. Substituting into $R_0$ yieldsThus a small stellar surface scale height leaves a wide interval above $R_{\ast }$ and below $R_0$ in which the thin disk flaring law can apply. The point-source approximation imposes the stronger lower cutoff obtained from $z_s\gg R_{\ast }$:This cutoff is still much smaller than $R_0$ when $R_0\gg R_{\ast }$. The full interval is therefore $R_{\ast }^{7/9}{R}_0^{2/9}\ll R\ll \min (R_0,R_M)$, whenever these bounds overlap. An assumed thickening radius must lie above the lower cutoff for all of the original approximations to hold simultaneously. The stellar luminosity-temperature relation is essential to the final comparison; independent arbitrary values of $L_{\ast }$ and $T_{\ast }$ would not imply it.

For a razor-thin disk, integrate the delta function in height to getThe corresponding horizontal force kernel is proportional to $(\mathbf{R}-{\mathbf{R}}^{\prime })/(|\mathbf{R}-{\mathbf{R}}^{\prime }{|}^2+{ϵ}^2{)}^{3/2}$. It removes the short-distance point-force singularity and smooths forces on separations of order $ϵ$. This is height-evaluation gravitational softening.

To represent finite vertical thickness, choose $ϵ$ of order the disc disk scale height $H$. There is no universal exact numerical choice: a true vertical mass density profile produces the Fourier reduction factor ${\Sigma }^{-1}\int \rho (z)e^{-k|z|}dz$, which is not generally $e^{-kϵ}$. Matching the long-wave term gives $ϵ=⟨|z|{⟩}_{\rho }$, the vertical-profile softening match. For an exponential vertical profile this mean is $H$, while its full reduction factor is $(1+kH{)}^{-1}$.

For a slowly rotating central star, the accretion-disk boundary layer radiates a luminosity $L_{\mathrm{BL}}\sim G{M}_{\ast }\dot{M}/(2r_{\ast })$. A layer of thickness comparable to the disk scale height has emitting area $A_{\mathrm{BL}}\sim r_{\ast }H$ times a geometric constant. Applying the Stefan–Boltzmann law to this blackbody area givesComparing with $\sigma T_{\mathrm{in}}^4=3G{M}_{\ast }\dot{M}/(8\pi r_{\ast }^3)$ gives the boundary-layer temperature scalingOnly the scaling is fixed: for a two-faced annulus of area $4\pi r_{\ast }H$, for example, $T_{\mathrm{BL}}^4/T_{\mathrm{in}}^4=r_{\ast }/(3H)$. A surface belt gives a different order-one coefficient.

For a thin disk, $H/r_{\ast }\ll 1$, so a comparable luminosity emerges from a smaller emitting area at a higher effective temperature. The boundary-layer emission is harder than the disk emission, with its Planck law peak shifted to higher frequency. It is also closer to a single-temperature component than the broad multitemperature blackbody disk spectrum, under the adopted uniform-temperature approximation. Rapid stellar rotation weakens the heating and this conclusion's temperature contrast.

Take $\gamma >1$ for the ideal gas. At fixed cylindrical radius, put $\mathcal{R}=k/(\mu m_p)$, where $k$ is the Boltzmann constant, $\mu$ the mean molecular weight and $m_p$ the proton mass. Also write $a=9\alpha \Omega /4$. Vertical hydrostatic equilibrium, thermal equilibrium and the ideal gas law requireConsequently, with $K=A/a$,Integrating from the midplane yieldsThe polytropic vertical structure in stellar gravity therefore has, for $|z|\le H$,whereThis is a polytropic equation of state $P\propto {\rho }^{1+1/\beta }$ with polytropic index $\beta$. For $|z|>H$, set $\rho =P=0$; temperature has no physical gas value in that vacuum. Here $H$ is the finite semi-thickness, rather than an exponential disk scale height.

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.

Vertical hydrostatic equilibrium givesDefine vertically integrated pressure and the density-weighted disk scale height byMultiply the hydrostatic equation by $z$ and integrate. Since $zp\to 0$ at both boundaries, integration by parts givesTherefore

The Shakura--Sunyaev thin disk model replaces poorly resolved turbulent angular-momentum transport by a stress proportional to pressure,Equivalently, an eddy viscosity has $\nu \sim v_{\mathrm{eddy}}{\ell }_{\mathrm{eddy}}$. Turbulent motions much faster than the sound speed $c_s$ would shock, while eddies much larger than the disk scale height $H$ would not fit within the disk. Writing $v_{\mathrm{eddy}}{\ell }_{\mathrm{eddy}}=\alpha c_{s}H$ therefore gives the alpha disk prescriptionwhere vertical hydrostatic balance gives $H\simeq c_s/{\Omega }_K$. The parameter $\alpha$ summarizes the correlation and efficiency of turbulent or magnetic stresses; it is not a molecular viscosity.

The scale $M_{\mathrm{th}}=h^3{M}_{\star }$ makes the Bondi radius comparable with the disk scale height. It organizes when planet-induced disturbances become nonlinear; it is not identical to an exact universal gap-opening mass.

For a fixed-mass planet in a circular orbit of radius $a$, let $q=M_p/M_{\star }$ be the planet-to-star mass ratio, $\Sigma$ the local surface density, and $h=H/a$ the disk aspect ratio defined by disk scale height $H$. The local angular velocity is $\Omega =\sqrt{G{M}_{\star }/a^3}$. Suppose the torque on the planet iswhere the dimensionless coefficient $C$ can have either sign and depends on the disk model.

The planet's orbital angular momentum is $L=M_p{a}^2\Omega =M_p\sqrt{G{M}_{\star }a}$. Its derivative gives $\dot{L}=(M_{p}a\Omega /2)\dot{a}=\Gamma$, hencePositive $\Gamma$ in this convention produces outward planetary migration; negative $\Gamma$ produces inward planetary migration. The torque exerted on the disk has the opposite sign. The radial evolution timescale is half $L/|\Gamma |$, because $L\propto \sqrt{a}$. This expression assumes the low-mass Type I migration regime and the stated torque law, rather than a universal direction or rate.

The vertical structure of an astrophysical disk describes how mass density, pressure and temperature vary away from its midplane. For a thin, centrally dominated gas layer, hydrostatic equilibrium gives $dp/dz=-\rho {\Omega }_z^{2}z$. An isothermal ideal gas has $p=c_s^2\rho$, producing a Gaussian density profile with disk scale height $H=c_s/{\Omega }_z$. Disk self-gravity or magnetic support changes that balance.^[ref]

The negative-diffusion criterion for viscous disk instability is ${\partial }_{\Sigma }(\bar{\nu }\Sigma )<0$. A surface density enhancement then transports angular momentum less effectively and grows. Freezing the Keplerian viscous diffusion equation locally gives growth proportional to the square of radial wavenumber, within the thin-disc transport range; wavelengths must still be sufficiently long compared with the disk scale height for that closure. The formal short-wave limit is a backward heat equation.