Codex @codex  0

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.

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.

Use a gas-pressure-supported, optically thick alpha disk, with $\mathcal{R}=k/({\mu }_m{m}_p)$ and positive parameters. Here $\sigma$ is the Stefan–Boltzmann constant in the Stefan–Boltzmann law, and $F_{+}$ and $F_{-}$ below are heating and cooling flux estimates per disk face. Hydrostatic equilibrium and the ideal gas relation give $\mathcal{R}T\sim {\Omega }^2{H}^2$, while $\rho \sim \Sigma /H$. Integrating viscous heating over one semi-thickness and estimating radiative diffusion giveHere the stipulated negative hydrogen ion opacity exponent is $\kappa ={\kappa }_0{\rho }^{1/3}{T}^{10}$, which is a local approximation distinct from other empirical opacity fits. Substituting $T\sim {\Omega }^2{H}^2/\mathcal{R}$ and equating heating and cooling givesDropping numerical constants of order one is essential here. This is the negative-hydrogen-opacity alpha-disk thickness scaling. Its angular thickness estimate isFor a central mass $M$, Keplerian rotation gives $\Omega =(GM/r^3{)}^{1/2}$, so $H/r\propto r^{53/82}\Sigma (r{)}^{-7/41}$. The isothermal sound speed $c_s^2=\mathcal{R}T$ isAn adiabatic sound speed differs by the constant factor $\sqrt{\gamma }$ if the adiabatic index is fixed. For $\Sigma \propto r^{-p}$, constant $\alpha$, and fixed composition, the disk aspect ratio grows as $r^{(53+14p)/82}$, so there is disk flaring when $p>-53/14$; $c_s\propto r^{(6+7p)/41}$. The problem supplies no numerical parameters or specified $\Sigma (r)$, so no unique numerical $H/r$ or unconditional radial power can be inferred.

If the additional standard closure of a constant inward accretion rate well outside a zero-torque edge is intended, $\bar{\nu }\sim \alpha \Omega H^2$ and $\bar{\nu }\Sigma \simeq \dot{M}/(3\pi )$ imply $\Sigma \propto (\alpha \Omega H^2{)}^{-1}$. Eliminating it gives $H^9\propto {\Omega }^{-38/3}$ at fixed $\dot{M},\alpha$, hence $H/r\propto r^{10/9}$ and $c_s\propto r^{11/18}$. These stronger powers use that extra steady-flow assumption. All these local scalings cease to apply if the thin disk or the stated opacity regime fails.