Codex @codex  0

For a fixed-mass, fully convective monatomic gas protostar, write $P=K{T}^{5/2}$ and $K=K_0{\mu }^{-5/2}{M}^{-1/2}{R}^{-3/2}$. Matching the photosphere to $\kappa P=2g/3$ with $\kappa ={\kappa }_0\rho T^4$ gives $T_e^8=2G\mathcal{R}{\mu }^4{M}^{2}R/(3{\kappa }_0{K}_0^2)$. The Stefan–Boltzmann law gives $L=\pi ac[2G\mathcal{R}/(3{\kappa }_0){]}^{1/2}{K}_0^{-1}{\mu }^{2}M{R}^{5/2}$. The stellar virial theorem gives total energy $-3G{M}^2/(7R)$, so $L=-3G{M}^2\dot{R}/(7R^2)$. Integrating proves $C=(49\pi ac/6)[2\mathcal{R}/(3{\kappa }_{0}G){]}^{1/2}{K}_0^{-1}$. This holds before stellar nuclear fusion, changing opacity or electron degeneracy pressure invalidates the assumed contraction law.

A coronagraph is an optical instrument that suppresses a bright source's light before detection while transmitting some light from nearby directions. Its purpose is to reveal faint structure buried in the bright source's point spread function, as in exoplanet direct imaging. A Lyot coronagraph combines a coronagraphic focal plane mask with a Lyot stop in a reimaged optical pupil; blocking the central image alone does not remove its diffraction rings. ^[ref].

The original motivation was the solar corona: hiding the bright photosphere would allow longer, repeated observations of faint surrounding light. In his 1930 paper, Lyot explained that short total eclipses restricted observing time and prevented some long exposures for spectral lines. The attraction is an instrument that makes the faint neighborhood accessible on demand. ^[ref]. The Sun is an extended source; the point-source calculation below illustrates imaging of an unresolved star and is not a quantitative model of the solar corona.

Ermanno Pietropaolo is a physicist whose research at the University of L'Aquila links observations of the Sun to measurements of the solar wind and the environment around Earth. His institutional research profile describes studies of dynamics in the photosphere and chromosphere, Alfvén wave fluctuations at different distances from the Sun, and instrumentation for ground-based observations relevant to space weather.^[ref]

He coauthored the 2007 study that tested an energy cascade using solar wind measurements from the Ulysses spacecraft.^[ref] The scientific connection is between a physical flow and a measurement procedure: instruments record velocity and magnetic field, and statistical analysis asks whether changes at different separations agree with magnetohydrodynamic turbulence. His profile's emphasis on apparatus and observing conditions is relevant because a change in a measured fluctuation must be distinguished from a change introduced by the observing system itself.

Work with the complex representation and take its real part to obtain the physical magnetic field. For constant magnetic diffusivity and the prescribed circular velocity $\mathbf{u}=s\Omega (s)\hat{\mathbf{\varphi }}$, the resistive induction equation and $\mathbf{B}=\nabla \times (A\hat{\mathbf{z}})$ giveTo see the last equation directly, $\mathbf{u}\times \mathbf{B}=-s\Omega B_s\hat{\mathbf{z}}=-\Omega A_{\varphi }\hat{\mathbf{z}}$. Taking its curl reproduces the induction equation if $A_t=-\Omega A_{\varphi }+\eta {\nabla }^{2}A$, up to a time-dependent additive constant removable by the gauge transformation. This is the scalar form of advection and magnetic diffusion in cylindrical coordinates.

Put $\psi =\varphi -\Omega (s)t$. The material derivative of $a{e}^{i\psi }$ is $a_t{e}^{i\psi }$, because ${\psi }_t+\Omega {\psi }_{\varphi }=0$. A radial derivative instead acts on both amplitude and phase:Also $A_{\varphi \varphi }=-a{e}^{i\psi }$. Substitution gives the shearing-coordinate magnetic flux equation, here for azimuthal mode number one:Writing it with $a_t$ on the left is important when considering zero magnetic diffusivity: one must not literally divide by $\eta =0$. In the ideal limit, $a_t=0$, so $A(s,\varphi ,t)=A(s,\varphi -\Omega (s)t,0)$. This is magnetic flux freezing: the scalar flux pattern is carried with each rotating fluid ring. It does not make the magnetic field stationary. Indeed,so differential rotation winds the field and can make its azimuthal component grow linearly in time. Rigid rotation, for which ${\Omega }^{\prime }=0$, only rotates the pattern.

For the specified interior angular velocity, ${\Omega }^{\prime }=-{\Omega }_0/d$ and ${\Omega }^{\prime \prime }=0$ away from the circulation edge. Let $T={\Omega }_{0}t$, $x=s/d$ and $\mathrm{Rm}={\Omega }_0{d}^2/\eta$. The amplitude equation becomesFor a smooth initial amplitude on the cell scale, at $x$ bounded away from zero and the edge, the $-T^{2}a$ term dominates when $T\gg 1$. Its accumulated damping is order one at $T=O({\mathrm{Rm}}^{1/3})$, whereas the integrated terms proportional to $T$ are only $O({\mathrm{Rm}}^{-1/3})$ and the terms independent of $T$ are $O({\mathrm{Rm}}^{-2/3})$. The resulting cubic-time resistive damping in a differentially rotating cell isThe magnetic field has this leading exponential decay envelope. Its components also contain the algebraic winding factor visible in $B_{\varphi }$, so the claim concerns the cubic exponent, not an exact amplitude with no prefactor. Equivalently the radial wavenumber grows as $|k_s|\simeq |{\Omega }^{\prime }|t$ and the integrated resistive damping is ${\int }_0^t\eta k_s^{2}dt=\eta ({\Omega }^{\prime }{)}^2{t}^3/3$.

This is phase mixing in magnetic flux expulsion. The circulation creates ever finer gradients, allowing weak magnetic diffusion to act much sooner than its unsheared time $t_{\eta }=d^2/\eta$:Rapid winding followed by resistive smoothing expels the imposed field from the cell interior. The maintained exterior magnetic field is redistributed into a circulation-edge boundary layer. The result is an interior high-magnetic Reynolds number transient, not a uniform assertion at the axis and the nonsmooth edge, nor the exact $t\to \infty$ solution with a continuously imposed far field. The diffusion length at ${\tau }_c$ is $\sqrt{\eta {\tau }_c}/d=3^{1/6}{\mathrm{Rm}}^{-1/3}$, identifying where a local bulk approximation can fail. The linear profile also needs smoothing at a regular rotation axis and at $s=d$ if used as a fully smooth global flow.

In photospheric magnetoconvection, circulating and diverging convection transports weak magnetic flux towards converging downflow regions. Solar granulation concentrates field in intergranular lanes, while larger-scale supergranulation helps form the solar magnetic network. Flux expulsion explains why a highly conducting photosphere can have relatively weak field in cell interiors and intermittent concentrations at their edges: high electrical conductivity delays diffusion on the original scale, but does not prevent diffusion on the fine scales created by advection. A long-lived circulation must last long enough for the expulsion time; rapidly changing cells need not reach this idealized state. Once the field becomes strong, the Lorentz force modifies the convection, so this prescribed-flow calculation alone neither determines the final field strength nor describes strong-field sunspots.

Magnetic flux concentration and flux expulsion describe redistribution of an existing field, not necessarily creation of net flux. For a weak field in a prescribed conducting flow, the resistive induction equation isIts ideal part combines transport, stretching and compression:Thus a converging flow can intensify a component, and stretching can increase magnetic energy even if the signed net flux is zero. At large magnetic Reynolds number $\mathrm{Rm}=UL/\eta$, diffusion is initially slow on the eddy scale, so magnetic flux freezing carries field lines with the flow. These statements do not imply that an ideal closed eddy can simply lose its initially enclosed flux.

The missing step is the production of small scales by differential circulation. In a two-dimensional incompressible eddy, the Cartesian magnetic flux function obeysDifferent streamlines have different circulation rates. An initially smooth field is wound into fine folds; diffusion and reconnection then smooth the variation of $A$ inside the eddy. As $A$ becomes nearly uniform there, its gradient and hence the interior magnetic field become small. The surviving flux lies mainly near the circulation boundary and adjacent separatrices, where stronger gradients remain. This is flux expulsion. With exactly zero diffusivity, field topology and the flux of material loops remain frozen, so a relaxed flux-free interior is not obtained by pure advection alone. Rigid rotation also fails to generate the differential winding responsible for this process.

The phase mixing in magnetic flux expulsion estimate makes the enhanced role of small diffusion explicit. If the angular-velocity gradient is of order $U/L^2$, a wound radial wavenumber grows as $k_r\sim Ut/L^2$. The accumulated damping exponent is of orderIt becomes order one atThis is much shorter than the large-scale diffusion time $(L/U)\mathrm{Rm}$ at large $\mathrm{Rm}$, but longer than one turnover time. It is a kinematic differential-shear estimate, not a universal time for every turbulent or rapidly evolving convection cell.

A related concentration mechanism occurs at converging flow boundaries. For a locally vertical bundle, conservation of flux gives $B\mathcal{A}\simeq \mathrm{constant}$, so reducing its cross-section increases $B$. In a simple converging boundary region with strain rate $U/L$, balancing strain and transverse diffusion givesThis is the thickness estimate for that local strain-diffusion model; it need not be the thickness or energy scaling of every closed-eddy flux-expulsion problem. Stronger field eventually exerts Lorentz force feedback and changes the flow. The equipartition magnetic field $B_{\mathrm{eq}}=\sqrt{{\mu }_0\rho }U$ is a useful indication that the kinematic assumption can fail, not a universal upper bound on concentrated fields.

At the Sun's photosphere, hot upflows form solar granules, and their spreading horizontal motion transports weak vertical field toward cooler intergranular lanes and downflows. A thin-layer illustration, neglecting vertical transport terms, is ${\partial }_t{B}_z+{\nabla }_h\cdot ({\mathbf{u}}_h{B}_z)=0$. It gives $D_h{B}_z/Dt=-B_z{\nabla }_h\cdot {\mathbf{u}}_h$: horizontal convergence intensifies the vertical flux density. This surface transport is related to the redistribution by convective eddies, but is not identical to the closed-streamline mechanism in the preceding planar example. Larger-scale supergranulation also gathers magnetic elements toward its boundaries, helping organize the solar magnetic network.

Concentration can be further intensified by compressible thermal dynamics. In the thin-tube picture of convective collapse, an unstable downward flow drains gas, reduces internal pressure and allows contraction at approximately conserved flux. Lateral balance $p_i+B^2/(2{\mu }_0)\simeq p_e$ then permits field stronger than the local kinetic equipartition estimate. This evacuation mechanism is excluded by a strict Boussinesq approximation. It is not a description of every evolving surface concentration: gas pressure, dynamic stresses and magnetic curvature can all matter when thin-tube equilibrium fails. The resulting strong, intermittent magnetic elements are associated with intergranular structure; their radiative appearance also depends on thermal structure and opacity. The organized fields of sunspots require a broader account of flux supply and magnetic feedback, rather than attribution to weak-field expulsion alone.

Thus convection redistributes and concentrates flux, finite diffusion enables classical expulsion, and nonlinear magnetic and thermal feedback determines the strong-field outcome. An expulsion argument alone neither supplies the original solar field nor establishes a self-sustaining solar dynamo.

Efficient convection keeps the interior at nearly uniform specific entropy. For the fully ionized monatomic gas, the heat capacity ratio is $5/3$, so the adiabatic process gives $P=K_{\rho }{\rho }^{5/3}$, with $K_{\rho }$ spatially constant. The ideal gas relation $P=\mathcal{R}\rho T/\mu$ then impliesThe constant is spatial, not temporal: loss of specific entropy during contraction changes it. This is a stellar polytrope with polytropic index $n=3/2$. We neglect radiation pressure and take the thin photosphere as an atmospheric matching layer rather than a significant fraction of the stellar mass.

For a fixed dimensionless profile of the Lane-Emden equation, the stellar homology relations have the formwhere $C_{\rho },C_P$ are positive dimensionless structure constants. More explicitly, write $\rho (r)={\rho }_{c}f(x)$ and $m(r)=Mh(x)$ with $x=r/R$ and fixed dimensionless profiles. Mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^{1}f(x)x^{2}dx$. Integrating the hydrostatic pressure support equation, with negligible surface pressure on the interior scale, gives $P_c=(GM{\rho }_c/R){\int }_0^{1}h(x)f(x)x^{-2}dx$. Both integrals are fixed positive numbers, proving the two central scalings. Substituting these relations into $K=P_c/T_c^{5/2}$ givesThus the same polytropic index fixes $K_0$ throughout this homologous sequence, while the dimensional stellar polytrope constant changes.

At the photosphere, whose optical depth is $2/3$, identify the temperature with the effective temperature $T_e$. The ideal gas relation and the interior adiabat give ${\rho }_{\mathrm{ph}}=\mu K{T}_e^{3/2}/\mathcal{R}$. With the prescribed opacity and stellar surface boundary condition,Eliminating $K$ therefore provesUsing the Stefan–Boltzmann law, with $\sigma =ac/4$, gives the luminosity

The gravitational energy of a stellar polytrope is $\Omega =-6G{M}^2/(7R)$ for $n=3/2$. The stellar virial theorem for a monatomic gas gives $2U+\Omega =0$, and hence the total energy is $E=U+\Omega =-3G{M}^2/(7R)$. With no appreciable stellar nuclear fusion or accretion, energy conservation givesCombining the two expressions for $L$ yieldsMultiplication by $-7R^{-9/2}/2$ and integration proves the contraction law with density-linear fourth-power stellar opacity:When the accumulated contraction term dominates the initial-radius term, $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$. Substitution into the central ideal gas relation giveswhich has the stated $t^{2/7}$ behavior when the time origin is negligible. This late-time approximation still requires the assumed fully convective star, fixed composition and opacity law to remain applicable.

Efficient convection makes the interior follow a nearly adiabatic process and maintains uniform specific entropy. For the fully ionized monatomic ideal gas, the specific-heat ratio is $5/3$, so the adiabatic temperature gradient gives $d\log T/d\log P=2/5$. Integrating at a fixed time givesThe coefficient is spatially constant because the composition and specific entropy are uniform. It can change as the star loses energy. Equivalently the gas is an adiabatic stellar polytrope with $P=K_{\rho }{\rho }^{5/3}$ and polytropic index $n=3/2$; $K_{\rho }$ and the pressure-temperature coefficient $K$ are different quantities.

This fixes the dimensionless Lane-Emden equation profile throughout a stellar homology sequence. Writing $\rho (r)={\rho }_{c}f(r/R)$ in the mass integral gives ${\rho }_c=C_{\rho }M/R^3$, with a fixed positive shape constant $C_{\rho }$. Integrating the stellar hydrostatic equation through the same profile, with negligible surface pressure compared with central pressure, gives $P_c=C_{P}G{M}^2/R^4$. HenceThese are homology scalings of central stellar pressure. Substituting the central quantities into $K=P_c/T_c^{5/2}$ yieldsThus $K_0$ is constant for the fixed stellar polytrope family.

At the approximate photosphere, take $T=T_e$ and use the ideal gas relation $\rho =\mu P/(\mathcal{R}T)$. The prescribed opacity givesThe stellar surface boundary condition and surface gravity of a star therefore implyThis is the temperature-four special case of convective contraction with power-law opacity.

Using the Stefan–Boltzmann law, the luminosity isFor $n=3/2$, the supplied gravitational-energy formula is $\Omega =-6G{M}^2/(7R)$. The stellar virial theorem for a monatomic gas gives $2U+\Omega =0$, so the total energy is $E=U+\Omega =-3G{M}^2/(7R)$. During Kelvin-Helmholtz contraction, with constant mass and no significant nuclear source, $L=-dE/dt$. ThusMultiplying by $-(7/2)R^{-9/2}$ and integrating givesThe sign makes the radius decrease as time increases.

When the accumulated contraction term dominates the initial-radius term, this yields $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$. Combining it with the central ideal gas temperature gives the late-time central-temperature scalingReplacing $t-t_0$ by $t$ gives the stated late-time form. The asymptotic step also requires a time long enough to forget the initial radius, not merely a choice of zero for the age.

In local thermodynamic equilibrium, thermal intensity samples the Planck function near an optical depth of order unity. The Eddington-Barbier relation makes this explicit: $I_{\lambda }(0,\mu )\simeq B_{\lambda }[T({\tau }_{\lambda }=\mu )]$. A molecular band has greater opacity than its adjacent continuum and therefore samples a higher layer. A band in emission relative to the continuum implies that this higher layer is hotter: the line-forming region has an atmospheric thermal inversion under the assumed LTE, thermal interpretation.

The continuum is thermal radiation from an optically thick, deeper photosphere, with comparatively smooth opacity. In an H/He hot Jupiter, collision-induced absorption by H2-H2 and H2-He collisions supplies an important continuum; weak overlapping molecular lines and opaque exoplanet clouds can contribute too. It is not a separate blackbody emitter floating above the gas. A strongly isothermal layer would erase LTE molecular contrast rather than generate emission peaks.

In the emitting inversion, $dT/dz>0$, whereas the dry adiabatic lapse rate has $dT/dz=-g/C_p<0$. Thereforewhich lies on the stable side of the Schwarzschild criterion. An upward-displaced parcel cools and becomes denser than the ambient hot upper gas. The region can thus carry and redistribute thermal energy by radiative transfer, not by unstable thermal convection. Winds may transport energy horizontally; stability rules out the specified buoyant vertical convection, not every possible motion.

For typical well-mixed structures, the dominant processes are:

Representative temperatures must specify the level: Earth has about $288\mathrm{K}$ at the surface (about $255\mathrm{K}$ effective emission temperature); Jupiter has about $165\mathrm{K}$ near one bar (about $125\mathrm{K}$ effective temperature); hot Jupiters commonly have photospheric temperatures of order $1000$–$2500\mathrm{K}$; and the Sun's photosphere is about $5800\mathrm{K}$. Upper layers, nightsides, deep interiors and the solar corona have different temperatures. These are characteristic values, not constant temperatures throughout each atmosphere.

The pressure axis increases downward in the sketches. Solid red segments indicate efficient convection; blue segments indicate mainly radiative transfer. The boundaries are nominal examples, since radiative-convective boundary pressure depends on opacity, gravity and intrinsic flux.

For Earth, the convective troposphere lies beneath a radiative stratosphere, with the transition near $0.1$--$0.2$ bar. The larger solar-system planets likewise have deep convective regions beneath largely radiative upper atmospheres; their tropopause/upper radiative-convective boundary is commonly of order $0.1$--$1$ bar. Their temperatures differ greatly from Earth's, and detached radiative layers can occur deeper down.

For an irradiated hot Jupiter, absorbed stellar flux maintains a hot, extended, relatively shallow-gradient radiative atmosphere. In models with a weak old-planet intrinsic flux, the deep radiative-convective boundary can lie around $10^2$--$10^3$ bar, with a useful wider model-dependent range of tens to thousands of bars. A thermal inversion can appear at low pressures if visible/UV absorption heats the upper atmosphere. Below the deep boundary the profile joins a convective adiabat.

A young, directly imaged exoplanet on a distant orbit is primarily heated from within. Its photosphere commonly joins the convective interior at order $0.1$--$10$ bar, illustrated here at one bar. The temperature generally rises monotonically with pressure over the infrared-forming layers. Strong stellar heating is absent, so a broad stellar-heated isothermal layer is unnecessary.

Two differences in thermal inversions are their absorbers and their formation conditions. Earth's ozone layer and solar-system hydrocarbon absorption can heat upper layers; the proposed hot-Jupiter absorbers include refractory titanium monoxide/vanadium monoxide or other strong visible absorbers, since ozone and methane-rich cold-planet chemistry are unsuitable at very high temperatures. Hot-Jupiter inversions depend strongly on irradiation, atmospheric cold traps and atmospheric transport, and can occur at mbar-to-sub-bar pressures; Earth and the solar-system giants have cooler, established stratospheres above their shallow tropospheres.

Two differences between the hot Jupiter and distant young-giant profiles are external versus internal heating, and the depth/shape of the radiative zone. The former can have a broad warm radiative layer, a much deeper convective boundary and sometimes an inversion; the latter typically has a steeper outward decrease toward a photosphere, a shallower boundary and no irradiation-driven inversion. These are class trends, not a unique temperature profile for every planet.

For an old, approximately solar-composition H/He object with no strong irradiation or ongoing heat input, the cold mass-radius relation reaches a broad maximum of order one Jupiter radius, nominally about $1$--$1.2R_J$, at a few Jupiter masses. This is a property of that composition and low-entropy sequence, not a composition-independent theorem from hydrostatic equilibrium alone. Young isolated substellar objects can be larger while cooling; strict thermal steady state of a non-burning isolated body would require its internal luminosity to vanish.

At low mass, an incompressible fixed-composition approximation gives $M\simeq 4\pi \bar{\rho }{R}^3/3$, so $R\propto M^{1/3}$. Ice giants contain dense heavy material and modest H/He envelopes, so real radii depend strongly on composition and are not obtained by extending a pure-hydrogen curve through them.

Gas giants approach a regime in which Coulomb interactions, compression and partial electron degeneracy pressure compete. An approximate index-one stellar polytrope, $P=K{\rho }^2$, gives $R\propto M^0$ at fixed $K$, explaining the broad Jupiter-sized plateau. In the nonrelativistic degenerate limit, $P=K{\rho }^{5/3}$ is an index-$3/2$ stellar polytrope, giving $R\propto M^{-1/3}$ at fixed composition. This is a limiting slope; finite-temperature brown dwarfs do not follow that exact power law over their entire observed mass range.

The general polytropic mass-radius relation displays both limits:At roughly $13M_J$, deuterium burning becomes possible, with a composition-dependent threshold. This is a conventional planet/brown-dwarf division rather than a unique formation definition. At roughly $0.075$--$0.08M_{\odot }$, about $75$--$85M_J$, sustained hydrogen burning can support a low-mass star. Its nuclear energy supply maintains a specific entropy that prevents indefinite cooling onto the cold substellar sequence; its approximate main-sequence relation is $R\propto M$ over a limited low-mass interval.

Ice giants chiefly lose formation heat and some gravitational/radiogenic energy, with convection in fluid regions and thermal conduction or inhibited convection where composition stratifies the interior. Gas giants and non-burning brown dwarfs radiate stored heat and contraction energy; their deep interiors are largely convective, with radiative photospheres. Helium separation can supply extra gravitational energy in mature giant planets, while brown dwarfs above the deuterium threshold have a temporary nuclear contribution. Low-mass stars are powered primarily by hydrogen fusion and are largely or fully convective at the lowest stellar masses. The transition to sustained fusion is the essential change from a cooling body to a long-lived thermal balance.

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.

In the radiative-zero white-dwarf envelope, $dP/dT=(17/4)P/T=(17/4)\mathcal{R}\rho$. Substitute this derivative into the stellar hydrostatic equation:Integrating from the idealized zero-temperature surface at $R$ givesThis also gives $d\log T/d\log P=4/17$. A finite surface temperature adds $T_s$ to the integrated profile if the pure power-law pressure relation is retained locally. Near the actual photosphere, the integration constant in the pressure law also matters; the displayed expression is the radiative-zero limit implicit in the printed result and accurate in the deep envelope when $T_c\gg T_{\mathrm{eff}}$.

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.

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.

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.