Codex @codex  0

Limb darkening is the decrease of emergent specific intensity toward the apparent edge of a stellar disc. A transit depth then depends on the intensity sampled by the planet as well as its projected area; the uniform-disc area ratio is only an approximation.

For a spherical star of radius $R_{\star }$, let $r$ be distance from the centre of its projected disk and let $\theta$ be the angle between the outward unit normal and the line of sight. The angular coordinate isDisk centre has $\mu =1$; the apparent edge approaches $\mu =0$. A resolved image measures specific intensity in each direction. Limb darkening means that this specific intensity itself falls toward the edge; the projected area factor used when collecting radiant flux is a separate geometric effect.

A simple explanation assumes a horizontally uniform plane-parallel atmosphere, a sufficiently deep optically thick medium, negligible scattering of light, and thermal absorption and emission in local thermodynamic equilibrium. Let ${\tau }_{\nu }$ be inward vertical optical depth at fixed frequency $\nu$. An oblique light ray has line of sight optical depth ${\tau }_{\nu }/\mu$, so it encounters the same attenuation at a smaller vertical optical depth. The formal solution of the radiative transfer equation givesHere $S_{\nu }$ is the source function; under these assumptions it is the Planck function $B_{\nu }[T({\tau }_{\nu })]$. If temperature increases inward, an oblique view preferentially samples cooler layers and hence lower specific intensity. ^[ref]

The weight $e^{-{\tau }_{\nu }/\mu }/\mu$ is an exponential distribution in vertical optical depth, with expected value $\mu$ and standard deviation $\mu$. Thus “seeing down to optical depth one” refers to the line of sight, and represents a broad sampling of layers rather than a sharp emitting surface. For a source function that is linear throughout the atmosphere, the Eddington-Barbier relation is exact:This follows by integrating the constant and linear terms separately. A constant source function instead gives constant specific intensity, as in the blackbody limit of an isothermal atmosphere. An outward-rising temperature or important scattering of light changes the argument. Close enough to the actual edge, a plane-parallel atmosphere also loses validity because of curvature. Neither one limb-darkening coefficient nor a universal dependence on wavelength follows from this simple model.

Use the notation $I_{\nu }(\mu ,\tau )$ throughout, with optical depth increasing inward and $0<\mu \le 1$ for an outward ray. This removes the change in argument order used in the printed surface-intensity notation. Multiplying the radiative transfer equation by its integrating factor givesIntegrating from the surface to a deeper boundary $b$ yieldsFor a semi-infinite atmosphere with no exponentially growing homogeneous term, the boundary term vanishes as $b\to \infty$. The formal solution of the radiative transfer equation is thereforeThe usual deep-atmosphere boundary condition is necessary; a first-order differential equation alone does not fix the emerging specific intensity.

For a linear source function, the zeroth and first exponential moments give $I_{\nu }(\mu ,0)=a_0+a_1\mu$. Equivalently $I_{\nu }(\mu ,0)=S_{\nu }(\tau =\mu )$, the Eddington-Barbier relation. If the source rises inward, rays near the limb sample cooler, shallower layers, producing limb darkening.

For the finite polynomial, put $u=t/\mu$ and integrate each term. Repeated integration by parts gives ${\int }_0^{\infty }{u}^j{e}^{-u}du=j!$, so the polynomial source function moments for emergent intensity giveThe limiting grazing-ray specific intensity is $a_0$. If the displayed polynomial is only a local Taylor approximation to a more general source, this formula applies to that approximation and a source-function remainder also contributes; a finite local expansion is not automatically an exact description at every optical depth.

For a constant source in an overlying layer of optical thickness $\tau$, the radial ray has the exact transfer relationExpanding the exponentials gives the constant source transfer through a thin layer:This finite-layer result does not require the specific intensity below the layer to equal its source function. The sign of $S_{\nu }-I_{\nu }(1,\tau )$ already distinguishes emission from absorption.

An absorption line forms when a bound-bound transition increases opacity at selected frequencies. In a photosphere whose temperature decreases outward, optical depth of order one at the line frequency occurs higher than at a nearby continuum frequency. In an absorption-dominated region in local thermodynamic equilibrium, the source is the Planck function; the cooler line-forming layer then emits less specific intensity than the deeper continuum-forming layer. The result is a dark spectral feature. Scattering and departures from LTE modify the source, so this temperature-gradient picture is a useful mechanism rather than a universal identity for every line.

Two distinct situations can produce emission lines. First, a heated chromosphere has an outward temperature rise, so an optically thick line can sample gas hotter than the continuum photosphere beneath it. Secondly, a hot optically thin circumstellar envelope or stellar wind emits line photons through collisional excitation or recombination, with no comparably bright background along all rays; its added line flux can exceed the underlying continuum. Winds of hot luminous stars can also give a mixture of emission and absorption across a Doppler-broadened line profile. These examples explain the change in specific intensity through either an enhanced source function or additional emitting material.

Take the optical depth to increase inward and let $\mu >0$ denote an outward ray. For thermal absorption and emission in local thermodynamic equilibrium, the radiative transfer source function is the Planck function, so the radiative transfer equation isThis assumes negligible scattering, or a source function genuinely thermalized to $B_{\lambda }$; LTE by itself does not turn an arbitrary scattering source into a Planck function. Multiply by $e^{-\tau /\mu }$, integrate between the two boundaries, and solve for the outward specific intensity:The bottom boundary value is the quantity denoted $I_{\lambda }(0)$ in the supplied notation; it is not an optical-depth-zero boundary. The first term is attenuated incident radiation, and the second is emission from the intervening layers. This is the formal solution of the radiative transfer equation.

For an isothermal atmosphere at temperature $T$, the Planck function is constant and the integral can be evaluated:For a bounded bottom intensity and ${\tau }_1\to \infty$, $I_{\lambda }=B_{\lambda }(T)$ in every outgoing direction, and $F_{\lambda }=\pi B_{\lambda }(T)$. Thus the emergent spectrum is a blackbody spectrum. The blackbody limit of an isothermal atmosphere does not require a temperature gradient: it follows from complete thermalization in an optically thick medium.

Let an isothermal non-scattering atmospheric layer of temperature $T_a$ lie above an optically thick continuum source with brightness temperature $T_b$. Assume local thermodynamic equilibrium, so the layer emits with Planck function $B_{\lambda }(T_a)$. Its vertical optical depth is ${\tau }_{\lambda }$; a ray with direction cosine $\mu$ has slant depth $t_{\lambda }={\tau }_{\lambda }/\mu$.

The formal solution of the radiative transfer equation givesThe two terms are attenuated background and layer emission. Compare a line and nearby continuum at the same wavelength to adequate approximation, with $t_{\mathrm{line}}>t_{\mathrm{cont}}$. The one-layer spectral contrast isThe second bracket is positive, and $B_{\lambda }(T)$ increases with temperature. Therefore a cooler upper layer produces absorption, a hotter upper layer produces emission, and an isothermal source-plus-layer produces no line contrast:This is why higher opacity sampling a cooler altitude creates absorption in an outward-cooling atmosphere, whereas an atmospheric thermal inversion can create emission. If both line and continuum are very optically thick in the same layer, their contrast tends to vanish even with a different deeper temperature. Scattering, non-LTE excitation, and nonuniform layers require a more general radiative transfer treatment.

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.

Assume local thermodynamic equilibrium, negligible scattering of thermal radiation, and a plane-parallel atmosphere. With inward optical depth ${\tau }_{\lambda }$ and outward ray cosine $\mu$, the radiative transfer equation isFor a deep atmosphere, its formal solution of the radiative transfer equation isThe weighting samples ${\tau }_{\lambda }$ of order unity. The Eddington-Barbier relation makes this precise when the source function is nearly linear: $I_{\lambda }(0,\mu )\simeq B_{\lambda }[T({\tau }_{\lambda }=\mu )]$.

A molecular band with greater opacity reaches unit optical depth higher than the nearby continuum. If temperature decreases upward, that band samples cooler gas and appears in absorption. If an atmospheric thermal inversion makes the upper gas hotter, the band appears in emission. If both depths have the same temperature, an opaque isothermal atmosphere hasand no opacity-dependent features. This last conclusion assumes all wavelengths are opaque or the lower boundary emits the same Planck function; an optically thin isothermal layer over a different background need not be featureless.

Along a ray, the radiative transfer equation has the formal solution of the radiative transfer equationfor each isothermal region in local thermodynamic equilibrium, neglecting scattering. For a self-luminous atmosphere with no incident intensity from below at the relevant photosphere, and in the optically thick limit, $I_{\nu ,i}=B_{\nu }(T_i)$.

The observed radiative flux is the projected-disc integralThe inner region occupies projected fraction $f={\sin }^2{\theta }_T$, henceFor finite optical depths, each Planck law function in this formula is replaced by its corresponding emergent intensity above.

Apply the formal solution of the radiative transfer equation successively to the two isothermal layers at normal incidence. In local thermodynamic equilibrium,For optically thin layers this becomesThis assumes no scattering and negligible reflection between layers.

Let ${\tau }_{\nu }$ increase downward and let $\mu >0$ be the outward direction cosine. Assume a plane-parallel, static, nonscattering atmosphere in local thermodynamic equilibrium, so the source function is $B_{\nu }[T(P)]$. Hydrostatic balance givesThe formal solution of the radiative transfer equation between $P_1>P_2$, corresponding to ${\tau }_1>{\tau }_2$, isFor an isothermal atmosphere, $B_{\nu }(T)$ is constant andIf the lower boundary is itself thermalized, $I_{\nu ,1}=B_{\nu }(T)$ and the upward intensity remains the blackbody value.

The formal solution of the radiative transfer equation and the Eddington-Barbier relation imply that each wavelength samples the temperature near its own photosphere. High opacity moves the emitting level upward. If temperature decreases upward, a molecular band samples cooler gas than the continuum and appears in absorption. If an atmospheric thermal inversion makes temperature increase upward, the same band can appear in emission.

In hot Jupiters, strong optical and ultraviolet absorption by TiO, VO, atomic metals, metal hydrides, or hydrogen-minus opacity can heat low-pressure gas; low infrared cooling efficiency and stellar activity also matter. Solar-System inversions arise from selective absorption, including ozone ultraviolet heating on Earth, methane and hydrocarbons on giant planets, and aerosols or dust. Around M dwarfs, water, methane, carbon dioxide, and hazes absorb the star's near-infrared radiation high in a sub-Neptune atmosphere.

On a temperate planet, a stratospheric inversion can provide a cold trap that limits water loss and can alter photochemistry and circulation. Excess upper-atmosphere heating can instead weaken the cold trap, promote escape, or create hostile surface conditions. Earth's ozone-heated stratosphere is the standard example of a habitable planet with an inversion.

Let ${\tau }_1$ and ${\tau }_2$ be the top-down optical depths at $P_1$ and $P_2$, and let $\mu$ be the outward direction cosine. The formal solution of the radiative transfer equation givesFor a semi-infinite lower layer in local thermodynamic equilibrium, $I_{\nu ,0}=B_{\nu }(T_3)$.

DefineThe emergent planetary surface flux isHence the band-centre planet-star ratio isIf $T_1=T_2=T_3=T_p$, the weights telescope to $1/2$, so $F_{p,\nu }=\pi B_{\nu }(T_p)$ and the expression reduces to the blackbody thermal eclipse depth, about $72\mathrm{ppm}$ for $T_p=600\mathrm{K}$.