Codex @codex  0

A hydrogen-burning shell converts hydrogen into helium in a layer surrounding a hydrogen-exhausted stellar core. It moves outward in enclosed mass as it supplies fresh helium to that core; its movement in radius need not have the same sign.

Let $m(r)$ be the background enclosed mass, so $g=Gm/r^2$ and $m^{\prime }=4\pi r^2\rho$. The background hydrostatic equilibrium is $p^{\prime }=-\rho g$. Use a Lagrangian fluid displacement $\xi ={\xi }_r$ of each spherical mass shell; $\Delta$ denotes its Lagrangian perturbation. Conservation of shell mass givesand the adiabatic pressure law gives the Lagrangian pressure perturbationThis retains perturbed self-gravity: a shell keeps its enclosed mass while its gravitational acceleration changes with its radius.

In mass coordinates, the exact radial shell equation is $\ddot{r}=-4\pi r^2\partial p/\partial m-Gm/r^2$. Its first-order perturbation isHydrostatic equilibrium gives ${\partial }_{m}p=-g/(4\pi r^2)$, and ${\partial }_m\Delta p=(\Delta p{)}^{\prime }/(4\pi r^2\rho )$. Thus the perturbation of the pressure-force geometry contributes another $2g\xi /r$, in addition to the gravity term:Put $\eta =\xi /r$, so ${\xi }^{\prime }+2\xi /r=r{\eta }^{\prime }+3\eta$. Expanding the derivative, using constant $\gamma$ and $p^{\prime }=-\rho g$, yields the radial stellar pulsation equation:The coefficient $4$ includes both geometric pressure and self-gravity effects; dropping the gravitational perturbation would not give this equation.

For a normal mode $\xi (r,t)=\xi (r)e^{-i\omega t}$, multiply by $\rho r^2{\xi }^{\ast }$ and integrate. Since ${\xi }^{\ast }/r={\eta }^{\ast }$, integration by parts givesRegularity at the center requires $\xi =O(r)$ with finite $\eta$, so the lower endpoint vanishes. For an isolated star take $p(R)=0$ and the free-surface condition $\Delta p=0$, with a regular finite-energy displacement; these eliminate the upper boundary term. HenceThese boundary conditions give the self-adjoint radial Sturm-Liouville problem behind the weighted stellar pulsation Rayleigh quotient.

Use the homologous trial fluid displacement ${\xi }_r=r$, or $\eta =1$. It is regular at the center, and $\Delta p=-3\gamma p$ vanishes at the free surface. The derivative term in the Rayleigh quotient is zero, leavingThe ratio is a positive weighted average of $g/r$. Since the central mass density is its maximum,For $\gamma >4/3$ multiplication preserves this inequality, proving the central-density upper bound for a radial stellar frequency:The full positive quadratic form also shows radial stability for constant $\gamma >4/3$ under these assumptions. A uniform-density equilibrium has constant $g/r$, making the homologous displacement an exact mode and saturating the bound.

If $\gamma <4/3$, the same trial value is strictly negative. Therefore ${\omega }_{\min }^2<0$: the star has an exponentially growing radial instability, with growth rate $\sqrt{-{\omega }_{\min }^2}$. One must use this negative trial quotient, rather than carrying the central-density inequality across the sign change unchanged. At the limiting value $\gamma =4/3$, the homologous displacement has zero restoring force and is a neutral radial mode in this constant-exponent model.

For a spherical stellar polytrope in hydrostatic equilibrium, combine $d{m}_r/dr=4\pi r^2\rho$ with $dP/dr=-G{m}_r\rho /r^2$ to eliminate the enclosed mass:For $n>0$, put $P=K{\rho }^{1+1/n}$, $\rho ={\rho }_c{\theta }^n$, $P=P_c{\theta }^{n+1}$ and $r=a\xi$. Since ${\rho }^{-1}dP/dr=(n+1)P_c{\theta }^{\prime }/(a{\rho }_c)$, chooseThe Lane-Emden equation is thenA regular center requires $\theta (0)=1$ and ${\theta }^{\prime }(0)=0$, with positive chosen central mass density and pressure. Locally $\theta =1-{\xi }^2/6+n{\xi }^4/120+\cdots$. The stellar surface is the first positive zero ${\xi }_1$, where the idealized external pressure is zero; retain the positive solution before it. Thus $R=a{\xi }_1$ and the Lane-Emden mass formula isThe surface condition selects where to stop a centrally regular solution, rather than replacing its central regularity conditions.

For a polytrope of index zero, the density is constant and $({\xi }^2{\theta }^{\prime }{)}^{\prime }=-{\xi }^2$. Regularity givesThe pressure is $P_c{\theta }_0$ and $a^2=P_c/(4\pi G{\rho }_c^2)$, so $R^2=3P_c/(2\pi G{\rho }_c^2)$. Index zero is the structural incompressible limit: the expression $K{\rho }^{1+1/n}$ is not itself defined at $n=0$.

For a polytrope of index one, set $u=\xi \theta$. The equation becomes $u^{\prime \prime }+u=0$, while central regularity requires $u(0)=0$, $u^{\prime }(0)=1$. HenceThe mass is $4{\pi }^2{a}^3{\rho }_c$, so it can change with central mass density while the radius stays fixed.

The moment of inertia of a polytropic star about any axis through its center follows by integrating $r^2{\sin }^2\vartheta$ over spherical shells:For constant mass density this gives $I_0=2M{R}^2/5$. For index one, integration by parts gives ${\int }_0^{\pi }{\xi }^3\sin \xi d\xi =\pi ({\pi }^2-6)$, and thereforeThese are axial moments of inertia, not the scalar second mass moment $\int r^{2}dm$.

For a finite-radius centrally regular stellar polytrope with $0<n<5$ and $n\ne 1$, eliminate ${\rho }_c$ from $R=a{\xi }_1$ and $M=4\pi a^3{\rho }_c{\omega }_n$. The resulting polytropic mass-radius relation isHere ${\xi }_1,{\omega }_n$ are dimensionless functions of $n$. At $n=1$ no such single-valued mass-as-a-power-of-radius relation exists at fixed $K$: the radius is fixed instead. At $n=3$ the exponent is zero and $M=4\pi {\omega }_3(K/(\pi G){)}^{3/2}$ is independent of central mass density. For the incompressible case, $M=(4\pi {\rho }_c/3)R^3$ at fixed mass density. Regular $n\ge 5$ solutions have no finite zero-pressure surface, so the finite-radius formula does not apply to them.

In stellar homology, the dimensionless radial profiles are the same after scaling radius, enclosed mass, pressure, temperature and luminosity. At corresponding radii let $r^{\prime }=\xi r$, $m_r^{\prime }=m{m}_r$, $T^{\prime }=tT$, $P^{\prime }=pP$ and $L_r^{\prime }=l{L}_r$. The printed $R^{\prime }=\xi r$ is understood as this local radial scaling, with surface radius $R^{\prime }=\xi R$.

Let $q={\rho }^{\prime }/\rho$ and keep composition, opacity coefficient and nuclear coefficient fixed first. Mass conservation gives $q=m/{\xi }^3$, and hydrostatic equilibrium gives $p=mq/\xi =m^2/{\xi }^4$. The ideal gas law then gives $t=p/q=m/\xi$. Scaling the stellar energy-generation rate equation and the stellar radiative temperature gradient gives respectivelyThe second follows from $dT/dr\propto -\kappa \rho L_r/(r^2{T}^3)$, not from energy production. Equating the two luminosity scalings givesThis assumes the denominator is nonzero and a consistent homologous family exists. If the mean molecular weight differs by $u={\mu }^{\prime }/\mu$, then $t=um/\xi$ andRatios of the opacity and energy-generation coefficients multiply the right-hand side. Thus fixed composition is a real restriction, not an automatic property of every stellar sequence.

Use the Stefan–Boltzmann law $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$ to locate the family on a Hertzsprung-Russell diagram. It implies $T_{\mathrm{eff}}\propto M^{(Y-2X)/4}$ andFor proton–proton chain burning take the usual local approximation $\eta =4$, with Kramers' opacity law $\lambda =1$, $\nu =-7/2$. ThenFor the CNO cycle with electron-scattering opacity, $\lambda =\nu =0$, soA conventional local choice $\eta =16$ gives $R\propto M^{15/19}$ and slope $76/9\simeq 8.44$. Choosing $\eta =17$ instead gives slope $60/7\simeq 8.57$. The source specifies no numerical nuclear exponents, and the effective exponent changes with temperature; the general expression is the unambiguous answer.

The Hertzsprung-Russell diagram places hotter stars to the left. Its branches rise toward higher luminosity and mass, and the CNO cycle/electron-scattering opacity branch has the steeper logarithmic slope. Their illustrative joining point and normalization are arbitrary because the proportional opacity and reaction laws do not specify absolute stellar scales. These are the fully radiative ideal-gas homology predictions, not exact observed main sequence relations: convection and increasing radiation support limit those assumptions in real stars.

Write $\mathcal{R}=k/(\mu H)$, where $H$ is the atomic mass unit used in the mean molecular weight convention, and let $\kappa$ be the Rosseland mean opacity per unit mass. For a static spherical star with negligible thermal storage, the stellar structure equations areHere $ϵ$ is the net stellar energy-generation rate per unit mass; $a$ is the radiation energy-density constant and $c$ the speed of light. The last differential equation follows from radiative diffusion, $F=-(c/3\kappa \rho )d(a{T}^4)/dr$, with $L=4\pi r^{2}F$. There is no gravothermal storage term in the specified stellar energy balance equation. Neglecting radiation pressure in force balance does not remove radiation as an energy carrier. Regular central conditions are $m(0)=L(0)=0$; take negligible external pressure, $P(R)=0$, for the analytic hydrostatic equilibrium model below.

Set $x=r/R$. Integrating the enclosed mass directly givesThusThe stellar hydrostatic equation then givesFactorizing makes the surface behaviour transparent:The ideal gas equation of state gives the temperature profileNo small-$x$ approximation was needed. In particular, close to the centre,These exact hydrostatic formulas reveal a limitation of the prescribed linear-density stellar model. Its mass density has a central cusp, and $T^{\prime }(0^{+})=T_c/R>0$. A positive opacity with positive outward luminosity would instead require $T^{\prime }<0$ in radiative diffusion. Moreover, regular positive nuclear heating gives $L=O(r^3)$ and therefore $T^{\prime }=O(r)$ at the centre. Thus the imposed profile is a hydrostatic toy model, not a complete positive-heating radiative-equilibrium solution of all the stellar structure equations. The requested algebraic profiles and nuclear-luminosity integral can still be computed consistently as properties of that toy model.

Write $\mathcal{R}=k/(\mu m_H)$ for the specific gas constant, with the printed $m_H$ denoting the atomic mass unit, and $a_{\mathrm{r}}=4\sigma /c$ for the radiation constant. In the low-mass envelope the enclosed mass is approximately $M$ and the outward luminosity is approximately $L$. The stellar hydrostatic equation, stellar radiative diffusion, and ideal gas equation giveThe Kramers' opacity law becomes $\kappa =({\kappa }_0/\mathcal{R})P{T}^{-9/2}$. Dividing the two differential equations eliminates radius and mass density:Integrating from a surface matching point $(P_s,T_s)$ givesDeep compared with the photosphere, the surface terms are negligible. The radiative-zero white-dwarf envelope therefore hasThis is the leading envelope approximation used in the later parts, rather than an exact atmospheric boundary condition at $T_{\mathrm{eff}}$.

At the base, temperature, mass density and pressure match the isothermal degenerate core. Let $K_{\mathrm{NR}}=K_1{\mu }_e^{-5/3}$, where ${\mu }_e$ is the mean molecular weight per electron. Equating the gas and electron degeneracy pressure givesHence $P_b={\mu }_e^{5/2}{\mathcal{R}}^{5/2}{K}_1^{-3/2}{T}_c^{5/2}$. Substituting this into the envelope relation at $T=T_c$ yieldsThe matching of a white-dwarf envelope to a degenerate core thus gives $L/M\propto T_c^{7/2}$, the envelope relation underlying Mestel's cooling law. It depends on pressure as well as mass density continuity at the interface; otherwise the supplied core equation would not determine the matching mass density.

For a slowly evolving spherical star, use a Lagrangian enclosed mass coordinate $m$. Let $u$ be specific thermal internal energy and $D/Dt$ follow a mass element. The stellar energy balance equation isThe outward luminosity is energy transported through the mass shell, ${ϵ}_{\mathrm{nuc}}$ is nuclear rest-mass energy released per unit mass and time, and ${ϵ}_{\nu }$ is escaping neutrino power per unit mass. Define the gravothermal stellar energy generation by ${ϵ}_{\mathrm{grav}}=-Du/Dt-PD(1/\rho )/Dt$. At fixed composition it is $-TDS/Dt$ by the first law of thermodynamics. During composition changes the full internal-energy derivative must use the appropriate equation of state, with nuclear rest-mass changes counted only once. If a tabulated nuclear rate already subtracts reaction-neutrino energy, that loss must not be subtracted a second time in ${ϵ}_{\nu }$.

In a nearly steady main sequence star, stellar nuclear fusion supplies most of the luminosity. Before stable hydrogen ignition, Kelvin-Helmholtz contraction releases gravitational energy. In a gas-supported hydrostatic star the stellar virial theorem gives $U\simeq -\Omega /2$, so about half the change in gravitational binding energy heats the star and half is available for radiation. The associated timescale is $G{M}^2/(RL)$. After nuclear burning ends, a white dwarf can shine by losing stored thermal energy; crystallization releases latent heat, and composition separation can add gravitational energy. These are heat and gravitational reservoirs rather than sustained hydrogen fusion.

Stellar nuclear fusion is possible because lighter nuclei can combine into products with larger nuclear binding energy. The energy per reaction is $Q=(m_{\mathrm{initial}}-m_{\mathrm{final}})c^2$, including the relevant particles consistently. Although thermal energies are below the Coulomb barrier, quantum tunnelling permits fusion. The thermonuclear reaction rate averages a cross-section over the thermal relative-speed distribution:For nonresonant charged-particle reactions the energy weighting contains $\exp [-E/(kT)-\sqrt{E_G/E}]$; the compromise between the thermal tail and penetration produces the Gamow peak. A resonance can greatly enhance the rate. Density, abundances and temperature therefore all affect the stellar energy-generation rate, and power laws such as $T^4$ or $T^{16}$ are local approximations.

For the proton–proton chain, first produce deuterium and helium-3:The initial reaction is slow because it converts a proton to a neutron through the weak interaction. It is the bottleneck that allows long hydrogen-burning lifetimes. The alternative pep reaction $p+e^{-}+p\to {}^2\mathrm{H}+{\nu }_e$ also feeds the same chain. There are three principal proton-proton chain branches. In pp I, two helium-3 nuclei terminate the chain:In pp II, helium-3 first captures an existing helium-4 nucleus, then electron capture and proton capture finish the branch:In pp III, the beryllium-7 instead captures a proton:Relative branch weights change with temperature and composition. All branches convert four protons into a net helium-4 nucleus, with two weak conversions and two neutrinos, allowing for the electron consumed in pp II. After including positron annihilation, the atomic-mass energy budget is approximately $26.7\mathrm{MeV}$ per net helium-4 nucleus, but the deposited heat is smaller by the branch-dependent escaping-neutrino energy. Near ordinary low-mass main-sequence conditions, ${ϵ}_{pp}\propto \rho X_H^2{T}^4$ is a useful local approximation.

The CNO cycle provides an alternative catalyzed hydrogen burning route in hotter cores. The CNO-I cycle, also called the CN cycle, consists ofThe carbon seed is regenerated: the net reaction again consumes four protons and produces one helium-4 nucleus, two positrons and two neutrinos. The slow ${}^{14}\mathrm{N}(p,\gamma ){}^{15}\mathrm{O}$ reaction governs the ordinary cycle and causes nitrogen-14 to accumulate. The CNO-II cycle branches through ${}^{15}\mathrm{N}(p,\gamma ){}^{16}\mathrm{O}$, ${}^{16}\mathrm{O}(p,\gamma ){}^{17}\mathrm{F}$, ${}^{17}\mathrm{F}\to {}^{17}\mathrm{O}+e^{+}+{\nu }_e$, and ${}^{17}\mathrm{O}(p,\alpha ){}^{14}\mathrm{N}$, after which the main cycle continues. Thus CNO nuclei are catalytic in the closed cycles, although their relative abundances change. A local rate is ${ϵ}_{\mathrm{CNO}}\propto \rho X_H{X}_{\mathrm{CNO}}{T}^{\eta }$ with a much steeper exponent than the PP chain; the preceding homology problem specifies $\eta =16$. The concentration of this heating near the centre favors a convective core in hotter, more massive hydrogen-burning stars.

Once hydrogen is exhausted in a core, contraction can raise its temperature enough for core helium burning. The Triple-alpha process overcomes the absence of stable mass-5 and mass-8 nuclei by maintaining a tiny transient beryllium-8 population:The Hoyle state, an excited carbon-12 resonance near the three-alpha threshold, makes the second capture efficient enough despite beryllium-8's very short lifetime. The net conversion is $3{}^4\mathrm{He}\to {}^{12}\mathrm{C}$, releasing about $7.27\mathrm{MeV}$. Because three alpha particles are required, the specific rate scales as ${\rho }^2{Y}^3$; a schematic resonant rate isIts local temperature exponent is $-3+44/T_8$, of order forty near $T_8=1$. A degenerate core cannot expand promptly to regulate this increase, giving a helium flash; in a nondegenerate core the stellar thermostat permits stable helium burning. The competing ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$ reaction determines much of the eventual carbon-oxygen mixture.

If the core becomes hot enough, carbon burning follows, for example through ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{20}\mathrm{Ne}+\alpha$ and ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{23}\mathrm{Na}+p$. Neon burning begins with photodisintegration, ${}^{20}\mathrm{Ne}+\gamma \to {}^{16}\mathrm{O}+\alpha$, followed by alpha capture such as ${}^{20}\mathrm{Ne}+\alpha \to {}^{24}\mathrm{Mg}+\gamma$. The first step consumes heat, but the combined rearrangement can release net energy. Oxygen burning includes ${}^{16}\mathrm{O}+{}^{16}\mathrm{O}\to {}^{28}\mathrm{Si}+\alpha$ and channels producing phosphorus and other nearby nuclei. At still higher temperatures, silicon burning proceeds through photodisintegration and particle captures in a reaction network approaching quasi-statistical equilibrium, producing iron-group nuclei. It is not simply direct silicon-plus-silicon fusion; the products depend on the electron fraction and weak-interaction timescale.

The increasing nuclear binding energy per nucleon supplies progressively less energy as products approach the iron group. Further fusion past that region cannot provide sustained net heat to support a core. Massive stars can therefore end with an unstable iron-group core; lower-mass stars do not attain all these burning stages and instead leave remnants such as carbon-oxygen white dwarfs. Heavy-element neutron captures can synthesize nuclei beyond the iron group without being the principal hydrostatic power source.

Finally, stellar neutrino energy loss competes with every heat source. Reaction neutrinos accompany the proton–proton and CNO chains. Hot or dense material also emits neutrino pairs through $e^{-}+e^{+}\to \nu +\bar{\nu }$, plasmon decay, $e^{-}+\gamma \to e^{-}+\nu +\bar{\nu }$, and electron-ion bremsstrahlung. These neutrinos usually escape, unlike the photons whose transport is diffusive; neutrino trapping requires much more extreme collapse conditions. Late burning has smaller fuel-energy reservoirs and strong neutrino cooling, so its duration is much shorter than core hydrogen burning. The surface luminosity is the escaping photon power, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$, not the sum of photons and unobserved neutrino power. Radiative diffusion in a star, convection and electron conduction redistribute energy, while the energy-balance equation distinguishes true sources, losses and storage. Accretion or mass loss adds boundary energy and mechanical work when present, rather than changing the nuclear reaction budget.

Assume a Newtonian, spherically symmetric, nonrotating body in hydrostatic equilibrium, supported by a polytropic equation of state $P=K{\rho }^{1+1/n}$ with constant $K>0$ and $n>0$. The enclosed mass and hydrostatic pressure support equation areEliminate $m$ by first writing $r^2{\rho }^{-1}dP/dr=-Gm$ and then differentiating:Let ${\rho }_c$ be central mass density, and introduce the Lane-Emden variables for a stellar polytropewhereThen $(1/\rho )dP/dr=(n+1)K{\rho }_c^{1/n}d\theta /dr$. Substitution cancels the dimensional factors and yields the Lane-Emden equationThe central conditions enforce the chosen central mass density and regular spherical symmetry. The local regular expansion is $\theta =1-{\xi }^2/6+n{\xi }^4/120+\cdots$. If a finite first zero ${\xi }_1$ exists and surface pressure is negligible, the physical radius is $R=a{\xi }_1$, and the Lane-Emden mass formula is $M=4\pi a^3{\rho }_c[-{\xi }_1^2{\theta }^{\prime }({\xi }_1)]$. For $0<n<5$ the standard isolated solution has a finite surface; $n=5$ has infinite extent, so a finite surface must not be assumed for all indices. Irradiation, composition stratification and non-polytropic equations of state require more general structure equations.

At the stellar surface boundary condition, match the interior to a stellar atmosphere, rather than setting the temperature and pressure to zero at an arbitrary radius. At a photosphere of radius $R$, the enclosed mass is $M$, the luminosity is $L$, and the outward radiative flux defines the effective temperature through $F=L/(4\pi R^2)={\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4$. With no external illumination the incoming specific intensity satisfies $I_{\nu }(0,\mu )=0$ for $\mu <0$. A deeper boundary is matched to the nearly isotropic radiative diffusion field. A grey atmosphere commonly matches near $\tau =2/3$; under the Eddington surface boundary condition this gives $T=T_{\mathrm{eff}}$. For roughly constant gravity and opacity, the total pressure rises by $g\tau /\kappa$ relative to its outer boundary value; if radiation supplies appreciable support, the gas-pressure gradient uses the effective gravity $g-\kappa F/c$ instead.

Take $z$ increasing outward and $\mu$ the outward direction cosine. If $j_{\nu }$ is emission per unit path length and solid angle, the radiative transfer equation isDefine inward-increasing optical depth by $d{\tau }_{\nu }=-\rho {\kappa }_{\nu }dz$ and the source function by $S_{\nu }=j_{\nu }/(\rho {\kappa }_{\nu })$. Then $\mu {\partial }_{{\tau }_{\nu }}{I}_{\nu }=I_{\nu }-S_{\nu }$. This states the emission-coefficient convention; if $j_{\nu }$ is defined per unit mass instead, the emission term is $\rho j_{\nu }$.

For frequency-integrated intensity, introduce the radiation-field momentsThe first angular approximation writes $I\simeq J+3\mu H$. Its odd part does not contribute to $K$, soThis is the Eddington closure approximation, not an exact description of the escaping angular distribution. An approximate surface condition takes the outgoing hemisphere to have direction-independent intensity $I_0$ and the incoming hemisphere to be dark. Its moments give $F=\pi I_0$ and $J(0)=I_0/2=F/(2\pi )$. Integrating the moment equations with the closure then gives $J=3F(\tau +2/3)/(4\pi )$; in thermal equilibrium, $J=a_{\mathrm{rad}}c{T}^4/(4\pi )$ gives $T^4=(3/4)T_{\mathrm{eff}}^4(\tau +2/3)$. This derives the usual approximate photospheric matching condition rather than imposing a zero surface temperature. For conservative coherent isotropic scattering, $S=J$. Integrating the transfer equation gives $dH/d\tau =0$ and $dK/d\tau =H$, the moment conditions for radiative equilibrium.

The exponential ansatz in the PDF has a genuine consistency problem. Compute its moments before making any approximation. For a nonsingular real intensity on all $-1\le \mu \le 1$, take $|ϵ|<1$ and defineThe continuous zero limits are $A_0=1$, $A_1=0$, $A_2=1/3$. Direct integration yieldsFor a decaying nontrivial correction, $\delta >0$ and $\gamma \ne 0$. Exact constant flux then requires $A_1=0$. But $\mathrm{atanh}ϵ/ϵ>1$ for every nonzero real $|ϵ|<1$, so $A_1=0$ forces $ϵ=0$. Thus exact radiative equilibrium cannot imply a small nonzero $ϵ$ for this ansatz.

Substitution into the full transfer equation makes the obstruction sharper. The exponential coefficient must satisfy, for every $\mu$,Comparing the constant and linear coefficients gives $A_0=1$ and $\delta =A_0ϵ=ϵ$. Hence regular real coefficients require $\delta =ϵ=0$, eliminating any decaying correction. This is the obstruction to a single smooth exponential mode in conservative grey transfer. There is no exact nontrivial real solution of the stipulated form with a small nonzero decay parameter. Choosing, for example, $ϵ=\delta =0.1$ and $\gamma \ne 0$ supplies an explicit counterexample to the claimed exact flux constancy, because $A_1\ne 0$.

Several requested conclusions remain useful as formal approximations, provided they are identified as such. If $\delta >0$, the large-depth exponential vanishes, giving $F\to 4\pi \alpha /3$, hence $\alpha =3F/(4\pi )$ for the asymptotic physical flux. It also gives $c{P}_{\mathrm{rad}}\to (4\pi /3)J$. Since $A_0=1+{ϵ}^2/3+O({ϵ}^4)$, neglecting second-order angular errors yields the formal relation $\delta \simeq ϵ$. However $A_1=-ϵ/3+O({ϵ}^3)$, so the finite-depth flux still varies at first order. The formal relation is not an exact transfer or constant-flux solution.

There is also a flux-normalization mismatch in the PDF's last displayed relation. With the same physical $F$ used above, the Hopf function of a grey atmosphere is defined byThe PDF's coefficient $3F/4$ would require redefining its $F$ as physical flux divided by $\pi$; it cannot simultaneously use the earlier $\alpha =3F/(4\pi )$ for the same physical flux. With the consistent physical-flux normalization, formal moment matching givesIn the intended small-$ϵ$ approximation this isThe approximate outgoing surface ansatz then gives the limb darkening ratioFor finite $ϵ$, $\gamma =-0.14/A_0$ gives a parameter-dependent ratio, not a unique exact value. The two Hopf endpoint values alone do not determine the exact emergent intensity: the formal solution of the radiative transfer equation uses the entire source profile, $I(0,\mu )={\int }_0^{\infty }J(t)e^{-t/\mu }dt/\mu$ for $\mu >0$. Nor does the simple linear-plus-exponential ansatz satisfy the exact no-incoming-radiation condition for every negative $\mu$. Thus $2.75$ is a clearly labelled approximate disk-centre-to-limb result, not an exact consequence of inconsistent premises.

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.