Past exam of the mathematics course of the University of Cambridge / 2018 / iii / Paper 317 / 2

In the thin radiative envelope, take enclosed mass $m(r)\simeq M$, luminosity $L_r\simeq L$, and constant mean molecular weight $\mu$. Combining radiative diffusion in a star with the hydrostatic pressure support equation givesThe ideal gas relation is $\rho =AP/T$, where $A=\mu m_u/k_B$, so the opacity law impliesHere $a_{\mathrm{rad}}$ is the radiation constant, distinct from the exponent $a$. Since $\nabla =(P/T)dT/dP$, separation gives $P^{a-1}dP=C{T}^{b-1}dT$. With photospheric $T=T_{\mathrm{eff}}$, the power-law opacity radiative envelope therefore hasThis expression assumes $a,b\ne 0$ and a positive bracket on the physical branch. If $b=0$, replace $(T^b-T_{\mathrm{eff}}^b)/b$ by $\log (T/T_{\mathrm{eff}})$; if $a=0$, the integrated relation is $\log (P/P_{\mathrm{phot}})=C(T^b-T_{\mathrm{eff}}^b)/b$, with the same logarithmic limit when $b=0$.

At the photosphere, use the stated boundary pressure and the effective temperature relation $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$. The radiation constant obeys $a_{\mathrm{rad}}c=4\sigma$ by the Stefan–Boltzmann law. Substitution into the radiative temperature gradient givesThis uses the approximate stellar surface boundary condition at optical depth $2/3$; it is not a universal photospheric gradient for arbitrary atmospheres.

For negative hydrogen ion opacity, $n=1/2$ and $s=-9$, so the exponents of the power-law opacity radiative envelope are $a=3/2$ and $b=-9/2$. Dividing the integrated pressure relation by $C{T}^b$ givesSince $a/b=-1/3$ and ${\nabla }_{\mathrm{phot}}=1/8$,It starts below the monatomic adiabatic temperature gradient and rises as the temperature increases inward.

For a monatomic ideal gas of fixed composition, the adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=2/5$. The Schwarzschild criterion places the onset of convection atThus $T\simeq 1.11T_{\mathrm{eff}}$. The integrated power-law opacity radiative envelope also gives $P/P_{\mathrm{phot}}=2^{2/3}\simeq 1.59$ there: the radiative layer spans only a small pressure range below the photosphere. Its geometrical depth is of order $H_P\log (1.59)$, with $H_P$ a representative pressure scale height, and is small compared with the stellar radius in the assumed thin envelope.

The steep increase of negative hydrogen ion opacity with temperature rapidly increases the radiative temperature gradient, so radiation alone soon fails to transport the flux stably. Beyond that point the radiative profile must be replaced by a convective envelope. Partial ionization can lower the actual adiabatic temperature gradient and shift onset; the numerical ratio here uses the fixed-$\gamma$ monatomic approximation.