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.

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.

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.

At a molecular line, the large opacity moves the optical-depth-one surface to lower pressure and higher altitude than the neighboring continuum. The Eddington-Barbier relation makes the emergent intensity approximately the Planck function at that layer. In an ordinary outward-cooling atmosphere, the line-forming layer is cooler and the feature is in absorption. In an atmospheric thermal inversion, it is hotter, so the spectral-line emission from an atmospheric thermal inversion exceeds the continuum brightness and the feature appears in emission.

Factors that can create or suppress inversions include the abundance of high-altitude optical absorbers such as TiO, VO, or atomic metals; the host star's irradiation level and spectral energy distribution; and clouds, hazes, composition, and day-night circulation, all of which alter where stellar and thermal radiation are absorbed.

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.

Without scattering and in local thermodynamic equilibrium, the radiative transfer equation has the emergent solutionThe Eddington-Barbier relation gives the useful approximationA molecular band has larger opacity than its neighboring continuum and therefore reaches optical depth unity at lower pressure.

If temperature decreases outward, the band samples cooler gas and appears in absorption. If the atmosphere is isothermal, both levels have the same source function and the feature disappears. If an atmospheric thermal inversion makes the upper layer hotter, the band appears in emission. For a weak separation of formation pressures,which explicitly shows that feature sign and amplitude measure the vertical temperature gradient.