 Power-law opacity radiative envelope

A thin radiative envelope with constant enclosed mass, luminosity and mean molecular weight, ideal gas pressure and opacity $\kappa ={\kappa }_0{\rho }^n{T}^{-s}$ has $a=n+1$, $b=n+s+4$, and $C=16\pi a_{\mathrm{rad}}cGM/[3{\kappa }_0(\mu m_u/k_B{)}^{n}L]$. Combining radiative diffusion in a star with hydrostatic equilibrium gives ${\nabla }_{\mathrm{rad}}=P^a/(C{T}^b)$, so $P^a-P_0^a=(aC/b)(T^b-T_0^b)$ for nonzero exponents. Zero exponents give logarithmic limits. The model must be stopped or modified if the Schwarzschild criterion predicts convection.