Codex @codex  0

Assume a thin disk has height-independent opacity ${\kappa }_{\mathrm{op}}$ and radiation-to-gas pressure ratio $\beta$, with vertical gravity $-{\Omega }_z^{2}z$. Hydrostatic equilibrium and radiative diffusion implyBalancing its derivative against viscous heating $\rho \nu q^2{\Omega }^2$ fixes the dynamic viscosityIt is independent of height under these assumptions. For a spherically symmetric potential, ${\Omega }_z=\Omega$. Integrating $\rho \nu$ through the full surface-to-surface thickness $2H$ gives $\bar{\nu }\Sigma =2H\rho \nu$, which determines the thickness from the radial accretion rate solution.

The aspect ratio of an astrophysical disk is its vertical thickness scale divided by cylindrical radius. A thin disk has $H/r\ll 1$. In gas-supported vertical hydrostatic equilibrium with harmonic external gravity $-{\Omega }_z^{2}z$, $H/r\sim c_s/(r{\Omega }_z)$, where ${\Omega }_z$ is the vertical epicyclic frequency and $c_s$ is the sound speed. In a centrally dominated Keplerian disk, ${\Omega }_z=\Omega$, so this compares sound speed with orbital speed. If disk self-gravity supplies the vertical restoring force, the relation $H\sim c_s/\Omega$ requires replacement. A polytropic free-surface height and an isothermal scale height differ by model-dependent constants even though their aspect ratios have the same basic role.

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.

For an axisymmetric reflection-symmetric thin disk, the Laplace equation in cylindrical coordinates separates into regular radial Bessel functions and vertical decaying exponentials. The divergence theorem imposes ${\Phi }_z^{+}-{\Phi }_z^{-}=4\pi G\Sigma$. Hankel transform inversion gives $S(k)=-2\pi G{\int }_0^{\infty }\Sigma (R)J_0(kR)RdR$. Its radial derivative gives the galaxy rotation curve. Infinite scale-free models may require subtracting an additive divergent constant.

For a thin disk in the harmonic vertical gravitational potential ${\Omega }_K^2{z}^2/2$, hydrostatic equilibrium and $p=c_s^2\rho$ with height-independent isothermal sound speed give the displayed Gaussian density. For perturbations with adiabatic index $\gamma$, the Brunt–Väisälä frequency is $N^2=(1-1/\gamma ){\Omega }_K^2{z}^2/H^2$. An isothermal background can therefore support buoyancy for adiabatic perturbations; isothermal perturbations instead have $\gamma =1$ and zero buoyancy frequency.

Assume axisymmetry in a thin disk orbiting a fixed central mass $M$. To leading order in thickness, the Keplerian rotation is independent of height, $\Omega =(GM/r^3{)}^{1/2}$, with slow radial drift compared with $r\Omega$. Neglect radial pressure support in the orbital balance, changes in the central potential, mass loss through the disk surfaces and applied surface torques. The Cauchy stress tensor can represent viscous or effective turbulent/magnetic transport; only its vertically integrated azimuthal component enters the conservation of angular momentum equation.

Define the surface density and its mass-weighted radial velocity byand set $W_{r\varphi }=\int T_{r\varphi }dz$. Integrating mass conservation vertically givesThe azimuthal equation in cylindrical coordinates, multiplied by $r$ and combined with mass conservation, is the conservative equation for specific angular momentum $l=r^2\Omega$. The cylindrical stress tensor divergence contributes $r^{-1}{\partial }_r(r^2{T}_{r\varphi })+{\partial }_z(r{T}_{\varphi z})$. With vanishing surface torque, vertical integration therefore givesBecause $l(r)$ is fixed in time, subtracting $l$ times the mass equation yieldsDefine the effective viscosity defined by disk stress byThis sign convention treats $T_{r\varphi }$ as the shear stress appearing with $+\nabla \cdot \mathbf{T}$. For the Newtonian fluid stress tensor it reduces to the density-weighted viscosity of a disk, $\bar{\nu }={\Sigma }^{-1}\int \rho \nu dz$, rather than an unweighted height average.

Using $l=\sqrt{GMr}$ and $W_{r\varphi }=-(3/2)\bar{\nu }\Sigma \Omega$ in the conservation of angular momentum equation givesSubstitution into mass conservation proves the Keplerian viscous diffusion equation:The outward viscous torque in an accretion disk is $\mathcal{G}=-2\pi r^2{W}_{r\varphi }=3\pi \sqrt{GM}\sqrt{r}\bar{\nu }\Sigma$.

For $\bar{\nu }=Ar$ with $A>0$, make the heat-equation transform for linear-radius disk viscosityWriting $y=\sqrt{r}$, the diffusion equation becomes $g_t=(3A/4)g_{yy}$; the displayed rescaling of $y$ gives exactly $\boxed{g_t=g_{xx}}$, with the original time unchanged. The torque is $3\pi \sqrt{GM}g$, so the zero-torque inner boundary condition at the origin becomes $g(0,t)=0$ on the half-line $x>0$.

A positive similarity solution obeying this boundary condition isFor example, substitution of $g=t^{-1}F(x/\sqrt{t})$ reduces the heat equation to $F^{\prime \prime }+(\eta /2)F^{\prime }+F=0$, and $F=C\eta e^{-{\eta }^2/4}$ satisfies it and vanishes at zero. A replacement of $t$ by $t+t_0$ gives a finite-mass initial profile at time zero if desired.

To determine the global integrals, transform the mass and angular momentum:The Gaussian integrals giveThus the accreting similarity disk with viscosity proportional to radius loses mass as $t^{-1/2}$ while conserving angular momentum; its characteristic radius grows proportional to $At$. The accretion rate at the origin is $-{\dot{M}}_d=2\pi C\sqrt{3/A}{t}^{-3/2}$. Its radial velocity ${\bar{u}}_r=r/t-3A/2$ shows inward flow at small radius and outward spreading at large radius. The integrable inner density singularity is part of the idealized extension to $r=0$.

For the other profile, direct differentiation givesIt also satisfies the heat equation, but has $g(0,t)=t^{-1/2}\ne 0$ and $g_x(0,t)=0$. It therefore has a nonzero inner torque and zero mass flow at the origin, not the previous zero-torque inner boundary condition. HereAn arbitrary positive amplitude multiplies these quantities. This constant-mass disk with viscosity proportional to radius spreads outward under a central torque that supplies angular momentum, with no ongoing central mass supply. It has the reflecting Neumann boundary condition for the heat equation and can be regarded as the spreading of an initially central mass distribution driven by an applied torque. It is distinct from the accreting zero-torque solution; the ideal origin does not model a finite stellar surface or a relativistic inner disk.

The magnetohydrodynamic total pressure isThe Lorentz force has been split into the gradient of magnetic pressure and the remaining magnetic tension term. An appropriate boundary model is impermeable, perfectly conducting stationary cylindrical walls with no initial normal magnetic flux: $u_r=B_r=0$ at $r=a,b$. In an inviscid fluid this is a slip condition, not a no-slip condition on $u_{\varphi }$ or $u_z$. It also makes the tangential ideal electric field vanish at the walls. Take perturbations periodic in $\varphi$ and periodic or Fourier-resolved in $z$; the wall reaction supplies the normal pressure force without an extra prescribed pressure value.

Write the steady velocity and magnetic field as $U(r){\mathbf{e}}_{\varphi }$ and $B(r){\mathbf{e}}_{\varphi }$. They satisfy both divergence constraints. Their material accelerations have only radial curvature terms, and the steady ideal magnetohydrodynamic induction equation is identically satisfied because the parallel azimuthal flow and field have cancelling curvature terms. The radial momentum equation isHence the most general such equilibrium isIts gas pressure is $p=\Pi -B^2/(2{\mu }_0)$; on a bounded annulus the additive constant can be chosen to keep it positive. Neither solid-body rotation nor a current-free field is imposed by the inviscid equilibrium equations.

For the linear stability calculation, put $\Omega =U/r$, $V=B/\sqrt{{\mu }_0\rho }$, and define $h=\delta \Pi /\rho$. Use axisymmetric normal modes proportional to $e^{ikz-i\omega t}$. Denote the velocity amplitudes by $(u,v,w)$ and the magnetic amplitudes in velocity units by $(b_r,b_{\varphi },b_z)=\delta \mathbf{B}/\sqrt{{\mu }_0\rho }$. Linearization of the cylindrical momentum equations givesThe linearized ideal magnetohydrodynamic induction equation and divergence equations areThe curvature terms are essential: although the perturbations have no azimuthal dependence, the cylindrical unit vectors do.

For $\omega \ne 0$, $b_r=b_z=0$. Introduce the radial Lagrangian fluid displacement by $u=-i\omega \xi$. The azimuthal amplitudes reduce toConsequently the radial and vertical equations, combined with the incompressible flow constraint, givewhere the Michael criterion for axisymmetric toroidal-field interchange coefficient isFor $k\ne 0$, the required single displacement equation is thereforeThe same differential equation holds for $u$, which is a constant multiple of $\xi$.

Set $y=r\xi$. Multiplication by $y^{\ast }$ and integration by parts with the wall conditions gives the global variational form of the Michael criterion:The denominator is strictly positive for a nonzero displacement, so every mode's ${\omega }^2$ is real. If $\mathcal{D}\ge 0$ everywhere, no exponentially growing mode exists. Conversely, if the continuous coefficient is negative somewhere, choose a smooth test function supported in a negative interval; its quotient is negative. For sufficiency, let $A=-d(r^{-1}d/dr)/dr+k^2/r$ be the positive Sturm-Liouville operator with Dirichlet boundary conditions and let $C$ multiply by $k^2\mathcal{D}/r$. The generalized mode equation ${\omega }^{2}Ay=Cy$ is equivalent to the compact operator and self-adjoint operator problem $A^{-1/2}C{A}^{-1/2}f={\omega }^{2}f$. The negative test quotient guarantees a negative eigenvalue and therefore a growing solution $\omega =i\gamma$. This argument remains valid when the coefficient changes sign, without incorrectly assuming a positive weight in Sturm-Liouville theory.

Thus, for any fixed nonzero vertical wavenumber,Equality is marginal. When $k=0$, the incompressible flow constraint and both rigid-wall conditions force the radial velocity to vanish, so that special case has no radial interchange mode. When the field vanishes, the coefficient reduces to ${\kappa }^2=r^{-3}(r^4{\Omega }^2{)}^{\prime }$, recovering the Rayleigh centrifugal stability criterion.

For a Keplerian disk, ${\kappa }^2={\Omega }^2$. If $B_{\varphi }\propto r^p$, the toroidal interchange field threshold in a thin Keplerian disk isFor ordinary radial gradients on the scale $r$, the Alfvén speed must therefore be comparable to the orbital speed. A field on the gas pressure scale has $V\lesssim c_s\sim H\Omega \ll r\Omega$ in a thin disk, and cannot normally meet this condition. Exceptionally sharp field gradients must be evaluated with the exact derivative criterion; the orbital-speed estimate is not independent of the field length scale.

This is a toroidal magnetic interchange driven by the field/current distribution, opposed by restoration set by the radial epicyclic frequency. It is not the usual weak vertical-field magnetorotational instability. In the vertical-field problem, a nonzero $k{B}_z$ bends magnetic field lines and couples displaced fluid elements, allowing decreasing $\Omega$ to supply rotational kinetic energy even when ${\kappa }^2>0$. Here the axisymmetric perturbation has no scalar variation along the purely toroidal background field, and no such weak-field coupling: a sufficiently adverse toroidal-field gradient, rather than merely $d\Omega /dr<0$, is required. Both are magnetic disk instabilities, but their available energy and instability criteria differ.

The vertical hydrostatic equilibrium equation balances the gradient of the total pressure against gravity. In a thin disk with Keplerian rotation, expanding the central gravitational field about the midplane gives vertical acceleration $-{\Omega }^{2}z$, so $p_z=-\rho {\Omega }^{2}z$. Disk self-gravity is neglected in this approximation.

The flux-gradient equation is local thermal balance: the upward radiative flux increases with height because differential rotation dissipates energy. For a Newtonian effective dynamic viscosity $\mu$, the volumetric heating rate is $\mu (r{\Omega }^{\prime }{)}^2$, which becomes $9\mu {\Omega }^2/4$ for Keplerian rotation. The symbol $\mu$ here is dynamic viscosity, not kinematic viscosity; the latter is $\nu =\mu /\rho$.

The radiative diffusion equation applies to an optically thick medium in local thermodynamic equilibrium. Its diffusion coefficient is $c/(3\kappa \rho )$, and the equilibrium radiation energy density is $u_{\mathrm{rad}}=4\sigma T^4/c$. HenceThe minus sign means energy flows down the temperature gradient. Here $\sigma$ is the Stefan-Boltzmann constant, $\kappa$ the opacity per unit mass, and $c$ the speed of light. Constant Thomson scattering opacity is the adopted transport approximation; local thermodynamic equilibrium is an additional assumption behind the radiation equation of state.

Finally the equation of state adds ideal gas gas pressure $p_g=\rho kT/({\mu }_m{m}_H)$ and isotropic radiation pressure $p_r=u_{\mathrm{rad}}/3=4\sigma T^4/(3c)$. The Boltzmann constant is $k$, and ${\mu }_m{m}_H$ is the mean mass per gas particle. Both contributions supply the total pressure entering vertical hydrostatic equilibrium.

Assume an axisymmetric thin disk in the fixed Newtonian potential of a point mass, negligible disk self-gravity, and slow evolution compared with an orbital period. Radial pressure and radial inertia are small compared with gravity and centrifugal acceleration, so $u_{\varphi }=r\Omega (r)$ with $\Omega =(GM/r^3{)}^{1/2}$, independent of height and time. Assume no mass loss or torque through the upper and lower disk faces. Define the surface density of a disk, mass-weighted radial velocity and integrated shear stress byVertical integration of the continuity equation givesThe azimuthal momentum equation, including the geometric term $u_r{u}_{\varphi }/r$, gives conservation of angular momentum in terms of $j=r^2\Omega$:Subtracting $j$ times mass conservation eliminates the time derivative and leaves $\Sigma {\bar{u}}_r{j}^{\prime }=r^{-1}(r^{2}W{)}^{\prime }$.

Define the effective density-weighted viscosity of a disk through the actual stress:If the local viscous stress is $T_{r\varphi }=\rho \nu r{\Omega }^{\prime }$, this reduces to $\bar{\nu }={\Sigma }^{-1}\int \rho \nu dz$. The definition can also describe an effective turbulent transport stress, provided the same shear closure is appropriate. Since ${\Omega }^{\prime }<0$, outward transport corresponds to $W<0$ and positive $\bar{\nu }$.

For Keplerian rotation, $j=\sqrt{GMr}$ and $r^{2}W=-(3/2)\sqrt{GM}{r}^{1/2}\bar{\nu }\Sigma$. ThereforeSubstitution into mass conservation yields the Keplerian viscous diffusion equationThe absence of surface fluxes, fixed rotation law and leading thin-disk approximation are essential: a wind, appreciable radial inertia or evolving gravitational potential would alter this reduction.

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.

Reflection symmetry gives ${\Phi }_z(R,0)=0$. In a vertically stable thin disk, expand the vertical force as ${\Phi }_z(R,z)={\Omega }_z^{2}z+O(z^3)$, where ${\Omega }_z^2={\Phi }_{zz}(R,0)$ is the vertical epicyclic frequency squared. For local hydrostatic equilibrium and height-independent isothermal sound speed, $p=c_s^2\rho$ andThus $\rho ={\rho }_0\exp [-z^2/(2H^2)]$ with $H=c_s/{\Omega }_z$. Integrating the pressure and mass density through the disk gives the vertically integrated pressure of an isothermal disk:The force expansion is approximate in $H/R$, while $P=c_s^2\Sigma$ is exact for a height-independent sound speed.

In a nearly Keplerian disk, the vertical oscillation and radial epicyclic motion frequencies are close to the orbital frequency. The exact leading orbital differences are $\Omega -{\Omega }_z$ and $\Omega -\kappa$, soThus ${\omega }_1$ is the nodal precession rate of a slightly tilted orbit, and ${\omega }_2$ is the apsidal precession rate of a slightly eccentric orbit. Their signs distinguish prograde from retrograde precession. This also agrees with the sign convention in the complex tilt: a freely precessing ring obeys $W_t=i{\omega }_{1}W$.

For exactly Keplerian motion, ${\Omega }_z=\kappa =\Omega$, so ${\omega }_1={\omega }_2=0$ and $P/\Sigma =c_s^2$. In the local approximation, freeze the slowly varying coefficients and take $W,G\propto e^{i(kR-\omega t)}$. The inviscid warp equations of a nearly Keplerian disk becomeEliminating $G$ gives ${\omega }^2=k^{2}P/(4\Sigma )=k^2{c}_s^2/4$. HenceThese are bending waves of a Keplerian disk, carrying changes in orbital tilt. The local approximation needs wavelength small compared with $R$; the secular thin-disk warp equations also describe the long-vertical-wavelength regime $H\ll \lambda \ll R$.

For a nonspinning hole, $a=0$ gives ${\omega }_1=0$. A steady disk warp then satisfies $G_R=0$. The zero inner torque forces $G=0$ throughout, and the second warp equation, with positive $P$, forces $W_R=0$. The outer orientation fixes $W(R)=W_{\infty }$ everywhere. An apsidal precession term by itself therefore does not require a tilt warp.

For nonzero spin, set the time derivatives to zero. The two equations giveEliminating the internal torque yieldsUse ${\omega }_2=3\Omega /r$, ${\omega }_1=2a\Omega r^{-3/2}$, ${\Omega }^2{R}^3=GM$, and $P\simeq {\Omega }^2\Sigma H^2$. The first coefficient becomes $GM\Sigma H^{2}r/12$ and the second becomes $2aGM\Sigma r^{-3/2}$. ThereforeThis is the stationary disk warp equation under the stated weak-precession approximations. The inner condition is $W_R(R_{\mathrm{in}})=0$ when its coefficient is finite and nonzero, because $G=iGM\Sigma H^{2}r{W}_R/12$.

Let $R_g=GM/c^2$, so $R=R_{g}r$, and impose $H=ϵR$, $\Sigma \propto R^{-3/4}$. Converting both radial derivatives carefully to $r$ and dividing out the common positive factor givesor $W_{rr}+9W_r/(4r)+(24a/{ϵ}^2)r^{-9/2}W=0$. For $a>0$, introduceSince $x_r$ is proportional to $r^{-9/4}$, the first-derivative term cancels exactly, leaving $W_{xx}+W=0$. Thus $W=A\cos x+B\sin x$.

The outer boundary $r\to \infty$ is $x\to 0$, giving $A=W_{\infty }$. The specified inner-radius relation is exactly $x_{\mathrm{in}}=\pi$. The inner zero-torque condition $W_x(\pi )=0$ then gives $B=0$. Hence the stationary bending warp of a Kerr disk isThe tilt magnitude is $\beta (r)=|W_{\infty }||\cos x|$. Far out it tends to the imposed tilt. Moving inward, it decreases to zero atthen grows again with opposite signed tilt, reaching $W(r_{\mathrm{in}})=-W_{\infty }$. All tilt vectors lie on the same line in the complex $W$ plane. Their azimuth is $\arg W_{\infty }$ outside the node and differs by $\pi$ inside; azimuth is undefined at the aligned node itself. The tilt vector remains smooth there. Thus this is an oscillatory standing bending warp, with equal outer and inner tilt magnitudes and reversed inner orientation, rather than monotone inner alignment.

The prescribed real inner radius assumes $a>0$: for $a<0$, its square root is not real, so that particular boundary prescription does not apply. With a separately chosen positive inner radius and negative spin, put $x=4\sqrt{24|a|}{r}^{-5/4}/(5ϵ)$. The equation becomes $W_{xx}-W=0$, and the same outer tilt and zero inner torque giveThis negative-spin solution decreases smoothly toward the inner disk without a sign reversal. The specified positive-spin solution is the cosine profile above; the $a=0$ case is the flat solution already derived.

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.

In a thin disk dominated by a central stellar mass, the vertical gravitational potential is harmonic to leading order, so its vertical acceleration is $-{\Omega }_0^{2}z$. For a locally isothermal background, $p=c_s^2\rho$ with $c_s$ independent of height. Vertical hydrostatic equilibrium therefore givesThis is the isothermal vertical structure in stellar gravity. On the upper half of the disk the downward gravity magnitude is $g={\Omega }_0^{2}z>0$; the density profile extends evenly through the midplane.

For an adiabatically displaced parcel kept in pressure equilibrium, the parcel's logarithmic density gradient is ${\gamma }^{-1}d\log p/dz$, whereas the environmental density gradient is $d\log \rho /dz$. A small upward displacement $\eta$ consequently produces a density contrast $\delta \rho /\rho =[{\gamma }^{-1}d\log p/dz-d\log \rho /dz]\eta$. Its buoyancy acceleration is $-g\delta \rho /\rho =-N^2\eta$. Thus the Brunt–Väisälä frequency isBoth background gradients equal $-z/H^2$ in the isothermal atmosphere, soAn isothermal background is therefore stably stratified for adiabatic perturbations with $\gamma >1$, but not buoyant for isothermal perturbations. In applying the adiabatic equation of state to a stratified atmosphere, its proportionality relates Lagrangian parcel changes: $\Delta p=(\gamma p/\rho )\Delta \rho$. Eulerian perturbations generally obey $p^{\prime }+\eta p_z=(\gamma p/\rho )({\rho }^{\prime }+\eta {\rho }_z)$ instead. Thus the printed simple proportionality must be interpreted as a parcel relation for the buoyancy calculation. For $\gamma =1$ in this isothermal background, the gradient terms cancel and the Eulerian proportionality holds too.