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 a quasistatic spherical star with fixed total mass, the stellar structure equations in radius areHere $s$ is the specific entropy of a moving mass shell, and the entropy term accounts for release or storage of thermal and gravitational energy. In thermal equilibrium it vanishes. An equation of state, opacity, energy-generation law and composition evolution close the equations. The absence of mass loss fixes the outer mass coordinate; it does not require the star to be in thermal equilibrium.

The radiative temperature gradient needed to carry the full luminosity isHere the actual temperature gradient is $\nabla =d\log T/d\log P$. For uniform composition the Schwarzschild criterion says that a radiative layer is stable if ${\nabla }_{\mathrm{rad}}<{\nabla }_{\mathrm{ad}}$, with marginal stability at equality, where the adiabatic temperature gradient is $(\partial \log T/\partial \log P{)}_s$. Use $\nabla ={\nabla }_{\mathrm{rad}}$ there. If ${\nabla }_{\mathrm{rad}}>{\nabla }_{\mathrm{ad}}$, convection carries some of the flux. Efficient convection gives $\nabla \simeq {\nabla }_{\mathrm{ad}}$; inefficient surface convection needs a transport prescription and can be superadiabatic. For composition gradients, use the Ledoux criterion, with threshold ${\nabla }_{\mathrm{ad}}+\varphi {\nabla }_{{\mu }_{\mathrm{mol}}}/{\delta }_{\mathrm{th}}$ instead. Here ${\nabla }_{{\mu }_{\mathrm{mol}}}=d\log {\mu }_{\mathrm{mol}}/d\log P$, $\varphi =(\partial \log \rho /\partial \log {\mu }_{\mathrm{mol}}{)}_{P,T}$ and ${\delta }_{\mathrm{th}}=-(\partial \log \rho /\partial \log T{)}_{P,{\mu }_{\mathrm{mol}}}$, using the mean molecular weight ${\mu }_{\mathrm{mol}}$.

Above the thin burning shell, take constant luminosity $L$, constant mean molecular weight ${\mu }_{\mathrm{mol}}$, ideal-gas pressure $P=\mathcal{R}T\rho$ with $\mathcal{R}=k_B/({\mu }_{\mathrm{mol}}{m}_u)$, and the Kramers' opacity law $\kappa ={\kappa }_0\rho T^{-7/2}$. Then the temperature equation isTo implement the printed assumption that $P,\rho$ and $M_r$ are all power laws while retaining mass conservation, write $M_r\propto r^m$, $\rho \propto r^b$, $P\propto r^p$ and $T\propto r^d$. The mass equation gives $m=b+3$; the ideal gas law gives $d=p-b$; hydrostatic equilibrium gives $d=m-1$. The diffusion equation gives $d-1=2b-2-13d/2$. Solving this linear system gives the self-gravitating Kramers power-law envelope:The coefficients can also be matched consistently. If $r_c$ is the core boundary and $M_c$ its enclosed mass,Radiative transfer fixes the constant luminosity throughwhose right-hand side is independent of radius for these exponents. Thus this is a solution of all four envelope equations, not just a dimensional estimate.

Matching the boundary temperature to the isothermal core and extrapolating the idealized power law to the specified photosphere givesFor reference $M_R/M_c=(R/r_c{)}^{1/11}=2500^{1/10}\simeq 2.19$, so the envelope mass is not negligible in this particular solution. The given $M_c\simeq M_{\odot }$ fixes the core-radius normalization $r_c=(11/42)G{M}_c/(\mathcal{R}10^7\mathrm{K})$; a numerical absolute radius needs the unspecified molecular weight. The exponent ratio also gives $\nabla =d/p=5/21<2/5$, so the fully ionized monatomic version of this idealization is convectively stable.

A different common approximation neglects the envelope's self-gravity and sets $M_r\simeq M_c$ in the force equation. It cannot obey the exact mass equation with a nonzero density and exactly constant $M_r$. Under that additional approximation, the same calculation instead givesThis is the constant-core-mass power-law opacity radiative envelope limit, not the full power-law solution with changing enclosed mass. Stating which approximation is used resolves the otherwise different numerical radius ratios.

A compact nearly isothermal core, a thin luminosity-producing shell and a much larger cool envelope describe a shell-burning red giant. The full power-law profile is a deliberately simple model. The Kramers' opacity law and an ionized ideal gas are not reliable at a $4000\mathrm{K}$ photosphere; partial ionization, other opacity sources and a convective envelope usually modify real cool giants. The quoted radius ratios are extrapolations within the stipulated model.

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.