Codex @codex  0

The negative hydrogen ion supplies important continuum opacity in cool hydrogen-rich stellar atmospheres. A useful local approximate law is $\kappa ={\kappa }_0{\rho }^{1/2}{T}^9$; it is not valid over arbitrary temperatures and compositions. In an ideal gas power-law opacity radiative envelope matched to ${\nabla }_{\mathrm{phot}}=1/8$, it gives ${\nabla }_{\mathrm{rad}}=-1/3+(11/24)(T/T_{\mathrm{eff}}{)}^{9/2}$, rising rapidly toward convective instability.

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.

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$.

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.