Codex @codex  0

For a grey atmosphere with $T^4=\frac{3}{4}{T}_{\mathrm{eff}}^4(\tau +2/3)$ and $P^2=P_0^2\log (1+3\tau /2)$, logarithmic differentiation givesA monatomic perfect gas of uniform composition has adiabatic temperature gradient $2/5$. The Schwarzschild criterion therefore predicts the first marginally unstable layer at ${\tau }_c=\frac{2}{3}(e^{4/5}-1)$, and convection for larger optical depth. Beyond that onset, the assumed purely radiative profile is no longer a self-consistent transport solution.

For no incident radiation and outward intensity approximated as direction independent, the surface angular moments are $F=\pi I_0$ and $P_{\mathrm{rad}}(0)=2\pi I_0/(3c)$. Integrating the radiation momentum equation with the Eddington closure approximation yields the grey atmosphere law $T^4=3T_{\mathrm{eff}}^4(\tau +2/3)/4$. This is an approximate angular boundary condition, not isotropy across both hemispheres.

Assume an ideal-gas grey atmosphere with $P(0)=0$, constant surface gravity of a star $g$, $T^4=T_e^4(1+3\tau /2)/2$ and $\kappa ={\kappa }_0\rho T^{4\beta +1}$ for $\beta >1$. Using $\rho =\mu P/(\mathcal{R}T)$ in $dP/d\tau =g/\kappa$ and integrating givesDifferentiation yields $d\log T/d\log P=[(1+3\tau /2{)}^{\beta -1}-1]/[2(\beta -1)]$. Applying the Schwarzschild criterion with monatomic adiabatic gradient $2/5$ gives ${\tau }_b=\frac{2}{3}[((4\beta +1)/5{)}^{1/(\beta -1)}-1]$. Since ${\tau }_b$ depends only on $\beta$, matching to $P=K{T}^{5/2}$ fixes $K\propto g^{1/2}{T}_e^{-2\beta -5/2}$.

Take an ideal gas grey atmosphere with $T^4=T_e^4(3\tau +2)/4$, constant surface gravity of a star $g$, vanishing external pressure and $\kappa ={\kappa }_0{P}^{1/2}{T}^8$. Integrating hydrostatic equilibrium in optical depth gives $P^{3/2}=8g[1/2-(3\tau +2{)}^{-1}]/({\kappa }_0{T}_e^8)$. Its logarithmic temperature-pressure gradient is $9\tau /16$, so the Schwarzschild criterion for a monatomic gas is crossed at ${\tau }_c=32/45$. The radiative continuation beyond that depth is unstable and must be replaced by convective transport.

For a grey atmosphere with $T^4=T_e^4(2+3\tau )/4$, an ideal gas and $\kappa ={\kappa }_0\rho T^{17}$, hydrostatic equilibrium in optical depth integrates to $P^2=512\mathcal{R}g[1/8-(2+3\tau {)}^{-3}]/(9\mu {\kappa }_0{T}_e^{16})$, taking $P(0)=0$. The logarithmic radiative gradient is $[(2+3\tau {)}^3-8]/48$. Equating it to the monatomic adiabatic temperature gradient $2/5$ gives $(2+3{\tau }_b{)}^3=136/5$ and the displayed onset temperature. Beyond this point the radiative profile must be replaced by convective transport.

In a thin ideal-gas grey atmosphere with $T^4=T_e^4(1/2+3\tau /4)$, constant enclosed mass $M$ and luminosity $L$, and $\kappa ={\kappa }_0\rho T^{13}$, divide the stellar hydrostatic equation by radiative diffusion in a star. This gives $PdP/dT=16\pi acGM\mathcal{R}{T}^{-9}/(3{\kappa }_{0}L\mu )$. Integrating from $P=0$, $T^4=T_e^4/2$, yieldsLogarithmic derivative gives $\nabla =T^8/T_e^8-1/4$. For a monatomic gas gas, the Schwarzschild criterion reaches equality at $\nabla =2/5$, hence $T/T_e=(13/20{)}^{1/8}$. The purely radiative profile cannot be continued past this point as a convectively stable model.

In a grey atmosphere with $T^4=T_0^4(1+3\tau /2)$ and $\kappa ={\kappa }_0{P}^{\alpha -1}{T}^{4-4\beta }$, hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. With $\alpha ,\beta >0$ and $P(0)=0$, integration gives the displayed pressure-temperature relation. At the matching surface $T=T_e$, the Eddington surface boundary condition gives $T_e^4=2T_0^4$, hence $P\kappa /g=4\alpha (1-2^{-\beta })/(3\beta )$.

The Hopf function encodes the departure of a plane-parallel grey atmosphere from a purely linear source profile. With physical bolometric radiative flux $F$ and standard mean intensity $J$, it is defined by $J=3F[\tau +q(\tau )]/(4\pi )$. In thermal radiative equilibrium, $T^4=3T_{\mathrm{eff}}^4[\tau +q(\tau )]/4$. The Eddington closure approximation with its approximate surface condition gives constant $q=2/3$; exact angular transfer generally gives a depth-dependent function.

Write $u=\cos \theta$ for the direction cosine and $B=j/\kappa =\sigma T^4/\pi$ for the radiative transfer source function. Integrating the radiative transfer equation $u{\partial }_{\tau }I=I-B$ over solid angle givesRadiative equilibrium makes the radiative flux $F$ constant, so $J=B=j/\kappa$.

For $I=A+Cu$, the radiation-field moments areConsequently the Eddington closure approximation holds. Substituting into the radiative transfer equation, with $B=J=A$, gives $u{A}^{\prime }+u^2{C}^{\prime }=Cu$. Equality for every $u$ requires $A^{\prime }=C$ and $C^{\prime }=0$, with $C=3F/(4\pi )$. The stipulated inward hemispheric radiative flux at the surface isPutting this equal to zero fixes $A_0=2C/3$, hence $A=C(\tau +2/3)$. Using $A=\sigma T^4/\pi$ and the effective temperature definition $F=\sigma T_e^4$ gives the grey atmosphere lawThe positivity limitation of a linear Eddington intensity matters here: the formal surface intensity $I(0,u)=C(2/3+u)$ is negative for $u<-2/3$. Thus this angular ansatz gives the requested Eddington surface boundary condition for approximate moments; vanishing integrated inward flux is not an exact, nonnegative, pointwise no-incoming radiation field.

Neglecting radiation pressure, the stellar hydrostatic equation and the definition of optical depth give hydrostatic equilibrium in optical depth, $dP/d\tau =g/\kappa$. Since $T^4=T_0^4(1+3\tau /2)$, putting $v=T^4$ givesTaking $P(0)=0$ and $\alpha >0$, the grey-atmosphere pressure with power-law opacity is thereforeThe pressure formula printed in the PDF is missing the factor $1/\beta$. Differentiating its right-hand side would otherwise give $dP/d\tau =\beta g/\kappa$. For $\beta =0$ the correct limit is $P^{\alpha }=2\alpha g\log (T^4/T_0^4)/(3{\kappa }_0{T}_0^4)$.

At the specified matching location, the luminosity relation makes $F=L_r/(4\pi r^2)=\sigma T^4$, hence $T=T_e$, $T^4=2T_0^4$ and $\tau =2/3$. Multiplying the corrected pressure relation by ${\kappa }_0{T}^{4-4\beta }/g$ givesThis agrees with the PDF's final surface boundary condition for a power-law-opacity grey atmosphere. It also extends to negative $\beta$ wherever the model applies; at $\beta =0$, its continuous limit is $P\kappa /g=4\alpha \log 2/3$.

Buoyancy and the convective criterion. Consider a small fluid parcel displaced outward by $\delta r$. Assume that sound waves keep its pressure equal to the ambient pressure, while heat exchange is slow enough that its displacement is an adiabatic process. Composition is uniform, so the parcel and its surroundings have the same mean molecular weight. For an ideal gas, $P\propto \rho T$, and an adiabatic parcel obeys $P\propto {\rho }^{\gamma }$. If $\delta \ln P<0$ is the pressure change on moving outward, thenThe density excess of the parcel at its new location is thereforeIf $\nabla >(\gamma -1)/\gamma$, this excess is negative: the outward-displaced parcel is lighter and accelerates farther outward. A downward displacement similarly continues downward. Conversely, a smaller gradient gives a restoring buoyancy force. The Schwarzschild criterion is thus stability for $\nabla <{\nabla }_{\mathrm{ad}}$ and instability for $\nabla >{\nabla }_{\mathrm{ad}}$, where the adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=(\gamma -1)/\gamma$. For a monatomic gas,A composition gradient would require the additional density effect in the Ledoux criterion; it is absent in this calculation.

Opacity and pressure. The supplied opacity is characteristic of a cool hydrogen-rich stellar atmosphere dominated by negative hydrogen ion opacity. Mostly neutral hydrogen can attach electrons supplied by partially ionized metals; bound-free and free-free absorption involving the resulting $H^{-}$ ions produces a steep positive temperature dependence. The weak density dependence and strong temperature sensitivity can be represented by local power-law fits. At fixed composition the given pressure form is equivalent to $\kappa \propto {\rho }^{1/2}{T}^{17/2}$, because $P\propto \rho T$. This is an approximate fit over the relevant cool-atmosphere regime, not an exact opacity law at all temperatures; molecular opacity or hydrogen ionization eventually changes the regime.

For a geometrically thin atmosphere take the surface gravity of a star to be constant, $g_s=GM/R^2$. Combining $dP/dr=-\rho g_s$ with $d\tau /dr=-\kappa \rho$ gives hydrostatic equilibrium in optical depth,The grey atmosphere relation gives $T^8=T_e^8(3\tau +2{)}^2/16$. With negligible external pressure, $P(0)=0$, integration yieldsAn imposed nonzero external pressure would add $P(0{)}^{3/2}$; the boundary assumption is what fixes the integration constant here.

The factor in parentheses is $3\tau /[2(3\tau +2)]$. HenceThe grey-atmosphere convection onset with eighth-power temperature opacity occurs when this reaches $2/5$, soThe radiative continuation below this boundary would be unstable. Efficient convection instead keeps the deeper star approximately adiabatic.

Matching to the convective interior. Write $P=K{\rho }^{5/3}$. Eliminating density using the ideal gas law givesat fixed composition. Since ${\tau }_c$ is a fixed number, the matching temperature $T_b$ is proportional to $T_e$, while the matching pressure satisfies $P_b^{3/2}\propto g_s{T}_e^{-8}$. ConsequentlyThe polytropic mass-radius relation for index $3/2$ gives $K\propto R{M}^{1/3}$. Raising the matching relation to the ninth power and using $g_s\propto M/R^2$ therefore givesFinally the Stefan–Boltzmann law $L=4\pi R^2{\sigma }_{\mathrm{SB}}{T}_e^4$ yields the Hayashi relation with eighth-power temperature opacity,Composition and the opacity normalization are fixed in this comparison. The entropy coefficient $K$ is spatially constant in each convective model, but need not be the same for stars with different masses or radii.

For the Schwarzschild criterion, displace a small gas parcel upwards. Assume its composition is unchanged, it remains in pressure balance with its surroundings, and it moves quickly enough that heat exchange is negligible. The parcel obeys the adiabatic equation of state $P{\rho }^{-\gamma }=\mathrm{constant}$. Combining this with the ideal gas equation gives $T\propto P^{(\gamma -1)/\gamma }$, so its adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=(\gamma -1)/\gamma =2/5$.

Let the surrounding temperature gradient be $\nabla =d\log T/d\log P$. For an upward displacement $\delta \log P<0$, the parcel's temperature relative to its new surroundings isIf $\nabla <{\nabla }_{\mathrm{ad}}$, this is negative: the parcel is cooler and, at equal pressure and composition, denser than its surroundings. It therefore sinks back; a downward displacement likewise gives a restoring buoyancy force. If the inequality is reversed, the displaced parcel continues to rise. The stability condition is thusEquality is the neutral boundary. There is no composition-gradient term because the ideal-gas parcel and the background have the same mean molecular weight.

In the thin upper stellar atmosphere, take a plane-parallel approximation with constant $g$, optical depth increasing inward, and zero overlying pressure $P(0)=0$. Hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. The ideal gas equation converts the stipulated opacity intoDefine $u=1+3\tau /2$, so the given grey atmosphere relation is $T^4=T_e^{4}u/2$. Integrating from the surface givesThis proves the required pressure law. In particular $P^2\propto g{T}_e^{-4\beta }$ times a function only of $u$ and $\beta$.

For the grey-atmosphere convection onset with density-linear opacity, differentiate the two profiles to obtain the radiative temperature gradient:It increases monotonically from zero, since $\beta >1$, and first reaches the adiabatic temperature gradient at $u_b^{\beta -1}=(4\beta +1)/5$. ThereforeBelow this point the radiative stratification is unstable and convection must transport some of the outward energy flux.

At the outer convection boundary, $u_b$ depends only on $\beta$, so $T_b=T_e(u_b/2{)}^{1/4}\propto T_e$ and $P_b\propto g^{1/2}{T}_e^{-2\beta }$ for fixed $\mu ,{\kappa }_0,\beta$. The matching constant in $P=K{T}^{5/2}$ is $K=P_b/T_b^{5/2}$. ConsequentlyThis temperature form of the polytropic relation follows from an ideal monatomic adiabat; its $K$ is not the density-polytrope coefficient in Question 3.

For the deep radiative boundary beneath a power-law-opacity convective envelope, negligible envelope mass and absent local energy generation give $m\simeq M$ and $L_r\simeq L$. Divide the radiative temperature equation by the pressure equation to get the stellar radiative temperature gradientThis is the gradient that would be needed if radiation alone carried $L$, and it is the quantity to compare with $2/5$ even while the actual region is convective. The deeper Kramers' opacity law gives $\kappa ={\kappa }_1\mu P{T}^{-9/2}/\mathcal{R}$. Substituting $P=K{T}^{5/2}$ yieldsAs $T$ rises inward this decreases, so the transition back to a stable radiative interior is atThus $T_{\mathrm{in}}^{7/2}\propto K^{2}L/M$. Using the outer matching and the Stefan–Boltzmann law,The mass and radius cancel. Hence, at fixed opacity coefficients and composition,This uses the same adiabat through the convection zone and distinguishes its inner transition temperature from its outer onset temperature.

Let $\Omega =-{\int }_0^{M}Gmdm/r$ be the gravitational energy. Multiply the stellar hydrostatic equation by $4\pi r^3$ and integrate through the star. Integration by parts and zero surface pressure givewhile its hydrostatic right-hand side is $-{\int }_0^{R}4\pi Gmr\rho dr=\Omega$. Hence the stellar virial theorem is

There is a genuine reference-sheet error to resolve. A monatomic ideal gas with adiabatic exponent $5/3$ has specific internal energy of an ideal gas $u=P/[\rho (5/3-1)]=3P/(2\rho )$, not the printed $3P/\rho$ on PDF page 2. The latter would give total internal energy $-\Omega$ and zero total energy, inconsistent with radiating contraction. The following derivation uses the consistent monatomic value.

A fully convective star is approximately isentropic and has an adiabatic stellar polytrope with index $3/2$. The given polytropic gravitational-energy formula yields $\Omega =-6G{M}^2/(7R)$. The virial relation gives total internal energy $U=-\Omega /2$, so total energy is $U+\Omega =-3G{M}^2/(7R)$. At constant mass with no nuclear heating, Kelvin-Helmholtz contraction obeys $d(U+\Omega )/ds=-L$, whenceHere $s$ denotes physical time.

In the plane-parallel atmosphere, the inward optical depth satisfies $d\tau =-\kappa \rho dr$, so hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. Neglecting radiation pressure, $\rho =\mu P/(\mathcal{R}T)$ and the specified opacity becomes $\kappa =({\kappa }_0\mu /\mathcal{R})P{T}^{16}$. ThereforeThe grey atmosphere relation and the boundary $P(0)=0$ integrate toTo apply the Schwarzschild criterion, put $w=2+3\tau$. Logarithmic differentiation givesThis increases inward, starting from zero. Marginal convection has ${\nabla }_{\mathrm{rad}}=2/5$, so $w_b^3=136/5$. Since $T_b^{12}/T_e^{12}=w_b^3/64$, the grey-atmosphere convection onset with seventeenth-power opacity isThe purely radiative atmosphere is used only up to this boundary.

The convective relation $P=K{T}^{5/2}$ is spatially uniform in its entropy coefficient $K$ at each instant; $K$ can change as the star radiates and contracts. At the convection boundary $w_b$ is a fixed number, so $P_b\propto g^{1/2}{T}_e^{-8}$ at fixed composition and opacity normalization. ThusOn the other hand, the central scalings for the $n=3/2$ stellar polytrope give $P_c\propto G{M}^2/R^4$ and $T_c=\mu P_c/(\mathcal{R}{\rho }_c)\propto \mu GM/(\mathcal{R}R)$. Hence, at fixed composition,Equating the interior and atmosphere expressions yields $T_e^{21/2}\propto M{R}^{1/2}$. Applying the Stefan–Boltzmann law then proves the Hayashi contraction with seventeenth-power opacity relation

At fixed mass write $L=A{M}^{8/21}{R}^{\alpha }$, with constant $A$ and $\alpha =46/21$. The contraction equation is $\dot{R}=-C{R}^{\alpha +2}$, where $C=7A{M}^{8/21}/(3G{M}^2)$. Direct integration givesMultiplying through by $R^{\alpha +1}$ and substituting the current luminosity givesFor the stated very large initial radius, or after that initial-radius term becomes negligible,This is the luminosity-age relation within the contraction model, rather than an assertion that a physical star begins at literally infinite radius.

Write $u=\cos \theta$, $S=j/\kappa$ and use axial symmetry. The radiative flux is $F=2\pi {\int }_{-1}^{1}uIdu$. Integrating the radiative transfer equation over directions givesThere are no energy sources in the atmosphere, so radiative equilibrium makes $F$ constant. Consequently

For the linear angular approximation $I=A+Cu$, the radiation-field moments areThus it satisfies the Eddington closure approximation, andSince $F$ is constant, $C^{\prime }=0$. Substitution in the radiative transfer equation, using $S=A$, gives $u{A}^{\prime }=Cu$, hence $\boxed{A^{\prime }=C}$.

In this moment approximation, the inward hemispheric flux at the surface isSetting this boundary moment to zero gives $A_0=2C/3$. ThereforeCombining it with the thermal source function and $F=\sigma T_e^4$ gives the grey atmosphere temperature law:

The surface condition is an approximate angular-moment condition. If the affine intensity were interpreted as an exact nonnegative boundary radiation field, no incident radiation would require $I(0,u)=0$ for every $u<0$; a linear polynomial vanishing on that interval would have $A_0=C=0$ and no outward flux either. The derived formal profile indeed has $I(0,-1)=-F/(4\pi )<0$. This positivity limitation of a linear Eddington intensity is why one uses the Eddington surface boundary condition as a moment closure, not the affine ansatz as an exact emergent angular intensity. An outward direction-independent intensity with no inward intensity gives the same boundary moments.

For the gas-pressure-dominated atmosphere, combine $dP/dz=-\rho g$ with $d\tau /dz=-\kappa \rho$. Thus the hydrostatic equilibrium in optical depth is $dP/d\tau =g/\kappa$. Let $w=T^4$. The grey atmosphere law implies $dw/d\tau =3T_0^4/2$. With the specified power-law opacity,At the outer vacuum limit take $P=0$ and $T=T_0$. Since $\alpha ,\beta >0$, integration givesThis is the grey-atmosphere pressure with power-law opacity relation.

At the matching radius the local radiative flux is $L_r/(4\pi r^2)=\sigma T^4$. In the thin plane-parallel atmosphere approximation this equals $F=\sigma T_e^4$, so $T=T_e$ and $T^4=2T_0^4$. Multiplying the pressure relation by ${\kappa }_0{T}^{4-4\beta }/g$ yields the required surface pressure-opacity relationHere $P\kappa$ is a product, not a pressure subscript.

Let $z$ increase outward and define the inward optical depth by $d\tau /dz=-\bar{\kappa }\rho$. With constant outward internal radiative flux $F=\sigma T_{\mathrm{int}}^4$, radiative diffusion givesThe constant $q$ is a boundary condition; diffusion alone does not determine it. For an unirradiated grey atmosphere with the Eddington closure approximation and Eddington surface boundary condition, $q=2/3$, soThe internal effective temperature of a planet is defined by its intrinsic cooling flux, not by its incident stellar heating. Since $dT/d\tau >0$ and $d\tau /dz<0$, temperature decreases outward. Towards the thin upper layers, $T\to 2^{-1/4}{T}_{\mathrm{int}}\simeq 0.84T_{\mathrm{int}}$; if the density and optical-depth gradient vanish there, $dT/dz\to 0$. This is the upper nearly isothermal atmosphere. The diffusion approximation itself fails at small optical depth: the grey boundary closure supplies the approximate continuation. Small irradiation perturbs this intrinsic-flux solution.

Young, self-luminous gas giants at wide orbital separations can have intrinsic cooling dominate their photospheric budget. They are favorable for exoplanet direct imaging, especially in the infrared, where their own thermal radiation is easier to separate from the host light. A strongly irradiated hot Jupiter instead has a large stable outer radiative region set mainly by stellar heating; it can be nearly isothermal over a broad pressure range or develop an atmospheric thermal inversion if stellar light is absorbed sufficiently high. Inversion is not inevitable for every strongly irradiated planet. Deep convection begins below its radiative-convective boundary.

The intrinsic curve follows the grey formula. The irradiated curves illustrate possible shapes only; they are not solutions for a specified opacity model. Smaller optical depth corresponds to greater altitude.

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.

Taking the logarithmic derivative of the grey atmosphere temperature law, $4\log T=\log (\tau +2/3)+\text{constant}$, gives

The outward radiative flux is $F=L/(4\pi r^2)$. A photon of energy $E$ carries momentum $E/c$, so the momentum transferred per unit volume and time to matter is $\kappa \rho F/c$, when $\kappa$ is the appropriate transport opacity. Momentum conservation makes the radiation stress decrease outward by this amount. In the thin atmosphere, where the radiation stress is treated by its scalar radiation pressure,This is the first moment of the radiative transfer equation. The grey atmosphere approximation uses a frequency-independent transport opacity; in an optically thick diffusion regime it becomes the Rosseland mean opacity. Near the surface the plane-parallel atmosphere approximation removes geometric stress-divergence corrections.

The stellar adiabatic exponent ${\Gamma }_2$ is defined at fixed specific entropy and composition byCombining parts (iii) and (iv), the radiative temperature gradient of the constant-opacity grey atmosphere isFor $1<{\Gamma }_2<4/3$, the adiabatic temperature gradient is strictly below $1/4$, and the Schwarzschild criterion is eventually violated. If ${\Gamma }_2$ is constant, solving for the crossing givesAt ${\Gamma }_2=4/3$ the radiative gradient approaches marginality from below, so there is no finite crossing in this model. If ${\Gamma }_2$ varies, the conclusion requires its subcritical value to persist in the deep layers. Uniform composition is needed for this use of the Schwarzschild criterion; a composition gradient requires the Ledoux criterion.

The Eddington closure approximation with the standard surface boundary condition gives the grey atmosphere profileTogether with the specified pressure law, logarithmic differentiation givesConsequently the radiative temperature gradient isFor a monatomic perfect gas, $\gamma =5/3$ and the adiabatic temperature gradient is $(\gamma -1)/\gamma =2/5$. Uniform composition makes the Schwarzschild criterion applicable. The first marginal layer therefore satisfies $\log (1+3{\tau }_c/2)=4/5$, giving the convection onset in a logarithmic-pressure grey atmosphereThe assumed profile becomes convectively unstable for $\tau >{\tau }_c$; further inward, a self-consistent model must account for convective transport.