Absorption line 2026-10-07
An absorption line is a deficit of observed spectral specific intensity relative to the surrounding continuum. In a photosphere with temperature falling outward, increased line opacity makes photons emerge from cooler layers than the continuum. For an absorption-dominated LTE source, the Planck function then supplies a smaller line specific intensity. A foreground layer with source below the incident continuum gives another direct transfer description.
Blackbody limit of an isothermal atmosphere Created 2026-10-06 Updated 2026-10-07
For a thermalized radiative transfer source function constant along a ray, the specific intensity after optical-depth interval is for . As , it tends to the Planck function, independently of a bounded boundary input. A temperature gradient is unnecessary for blackbody emission. In a scattering medium, local thermodynamic equilibrium alone does not guarantee ; thermalization must be justified.
Consider a hypothetical star whose emitted specific intensity is exactly that of a uniform blackbody. Assume transparent propagation with negligible Doppler shift and interstellar extinction, so geometric dilution multiplies every spectral interval by the same factor. The stipulated temperature gives
Using the converted frequency peak as though it maximized would return , a substantial error caused solely by mixing spectral conventions.
Now stipulate a photodetector with constant quantum efficiency over the relevant wavelengths, collecting equal narrow wavelength bins. For fixed accepted projected area times solid angle , with dilution included in , its expected photon count rate density is
This follows by dividing the received power by the photon energy . It peaks where , at , giving for the same model star. Constant quantum efficiency is not constant energy spectral responsivity: the latter includes the wavelength-dependent photon energy conversion. For example, , and a bin centered at give approximately detected photons per second in the narrow-bin approximation. These are calculated model values, not observations of a named star.
A more general linear detector channel has output
where is the background reading and includes spectral selection, optical transmission and detector spectral responsivity. A raw channel maximum need not locate either Planck function maximum. Recovering a temperature requires an appropriate calibrated response and an emission model. For a real stellar atmosphere, opacity and radiative transfer can reshape the spectrum; interstellar extinction can reshape it again. A single peak does not establish a unique physical temperature without those assumptions.
Figure 1.
One model star, three peak conventions, and a spectral measurement chain
. All curves use the same Planck function and are normalized separately. The curve retains density per frequency while displaying converted wavelength. The photon curve assumes constant quantum efficiency in equal narrow wavelength bins. The lower diagram shows a generic measurement sequence, not a reconstruction of a particular historical apparatus. Details.
At fixed , substitute into the Planck function:
Here is the Planck constant, the Boltzmann constant, and the speed of light. Since , an interior maximum in wavelength is also a stationary point in . A logarithmic derivative gives
The required positive root is . The rearranged equation also has , but that endpoint is not a peak. For , changes sign once, at . Thus first increases and then decreases to negative values, crossing zero exactly once for . Since tends to zero at both endpoints, this stationary point is the unique global maximum. Solving the definition of for wavelength gives
The inverse-temperature dependence follows before any numerical root-finding. More generally for one fixed function , so spectra at all temperatures have the same rescaled shape. This is stronger than a relation between two peak positions. Replacing the denominator by in the Wien approximation would instead give and a peak wavelength about too short.
A brightness temperature equates a measured intensity with a Planck function; it does not determine its contribution pressure. In hydrostatic equilibrium, . Estimating therefore needs gravity and opacity, as well as a model of the contribution function. Different pressures can produce the same band brightness.
At a fixed temperature, the frequency maximizing and the wavelength maximizing in Planck's law are different ways of describing the same black-body radiation. They are maxima of different densities. The inverse relation between frequency and wavelength does not preserve those maxima because its interval-conversion factor varies with wavelength.
For a short extremum calculation, substitute into , leaving a positive constant times . Substitute into , leaving a positive constant times . Setting each derivative to zero gives
Select the positive roots: and . Although solves the rearranged equations, it is an endpoint, not either peak. For or , the function initially increases, its derivative crosses zero once, and it eventually becomes negative. Hence it has exactly one positive root. Each density tends to zero at both ends, so this root gives its unique maximum. The bisection method brackets these roots in and .
Consequently Wien's displacement law gives
These constants agree with NIST's tabulation. [ref] The wavelength peaks at , , and are , , and . Doubling temperature halves this peak wavelength.
At , the frequency peak is , whose converted wavelength is . The peak is instead , corresponding to about . These numbers are compatible: equal-width frequency bins and equal-width wavelength bins partition the same radiation differently. Specifying which density was maximized is essential when using a peak to infer temperature.
Figure 1.
One spectrum, two spectral-coordinate peaks
. The left panel plots the Planck function per terahertz and marks its frequency peak. The right panel plots the same black-body radiation per micrometre and marks its wavelength peak. Dashed lines locate the counterpart of the other panel's peak. The converted frequency peak at is not the peak at : the change of variables changes density heights through . Details.
Limb darkening Created 2026-10-06 Updated 2026-10-14
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 , let be distance from the centre of its projected disk and let be the angle between the outward unit normal and the line of sight. The angular coordinate is
Disk centre has ; the apparent edge approaches . 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 be inward vertical optical depth at fixed frequency . An oblique light ray has line of sight optical depth , so it encounters the same attenuation at a smaller vertical optical depth. The formal solution of the radiative transfer equation gives
Here is the source function; under these assumptions it is the Planck function . If temperature increases inward, an oblique view preferentially samples cooler layers and hence lower specific intensity. [ref]
The weight is an exponential distribution in vertical optical depth, with expected value and standard deviation . 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.
When radiative rates compete with or dominate collisions, level populations need statistical equilibrium coupled to radiative transfer, rather than thermal population laws at the local temperature. Radiation from other depths can make the local radiation field differ from the Planck function. This changes ionization and spectral-line strengths, so a geometrically thin atmosphere may still require non-LTE abundance modelling.
A non-scattering isothermal atmospheric layer in local thermodynamic equilibrium, of temperature and slant optical depth , lies above a blackbody background at . The specific intensity is
Thus an increase in opacity shifts the specific intensity toward the layer's Planck function. If a line is more opaque than its continuum, it appears in absorption for , in emission for , and has no contrast for . This is a local explanation, not a claim that all features diagnose one uniform layer.
The Rosseland mean opacity is the transport opacity appropriate to an optically thick region close to local thermodynamic equilibrium. If is monochromatic transport opacity and the Planck function, then
It is a weighted harmonic mean: low-opacity frequency windows carry much of the diffusive flux. One must combine the monochromatic processes before averaging, rather than adding their separately averaged opacities.
In hot fully ionized material, electron-scattering opacity supplies an almost density- and temperature-independent contribution at fixed composition. In the nonrelativistic Thomson scattering limit, . Free-free opacity, the absorption counterpart of thermal bremsstrahlung, comes from photon absorption during encounters between free electrons and ions. Its approximate mean obeys Kramers' opacity law, , with composition and quantum correction factors.
When ions retain bound electrons, bound-free absorption is photoionization; edges, occupations and ionization fractions strongly affect its frequency dependence. Bound-bound absorption produces absorption lines. Numerous metal lines can close low-opacity windows and produce important opacity features in partially ionized layers. Line populations depend on excitation as well as ionization, while Doppler broadening, pressure broadening and overlapping lines change the transport windows. Stimulated inverse processes must be included consistently with detailed balance and the local thermal radiation field.
In cooler hydrogen-rich atmospheres, negative hydrogen ion opacity from both bound-free and free-free processes is often important. Its strong temperature dependence helps determine convective surface relations. In still cooler layers, molecular bands, and at sufficiently low temperatures dust absorption/scattering, may be important. These processes cannot be represented by one Kramers power law over the whole density-temperature plane.
An opacity calculation therefore needs the chemical abundances, the equation of state, electron density and bound-state populations. Saha ionization equation and thermal excitation populations give useful dilute LTE estimates; high-density pressure ionization, nonideal effects and altered level populations require corrections. Composition matters even when one writes only : hydrogen, helium and the metal abundance pattern are additional fixed inputs.
In dense degenerate cores, thermal conduction by electrons can dominate energy transport. An equivalent conductive opacity may be introduced so that , because the two diffusive fluxes add; conduction is not itself a photon absorption process. Very hot regimes also require departures from the simple nonrelativistic scattering approximation. Near the optically thin photosphere, a Rosseland diffusion description alone is inadequate, and frequency-dependent radiative transfer with appropriate boundary conditions is required. These regime changes are why realistic stellar opacity tables, rather than one global power law, enter stellar structure.
Write , with at linear order. Since , the background Planck function depends only on comoving momentum. Expanding its temperature-shifted argument gives
Using rather than in the perturbation denominator changes only higher-order terms.
The photon brightness perturbation is the energy-weighted angular perturbation
Integration by parts gives ; the boundary term vanishes at both endpoints for a Planck spectrum. Therefore . Multiplying the transport equation by and integrating yields the synchronous photon brightness equation
The PDF's definition is ; the TeX's denominator is a transcription error.
Expand in Legendre moments in . The monopole is the density/temperature perturbation, and the dipole is the bulk velocity. Before decoupling, rapid Thomson scattering isotropizes photons in the electron rest frame. In the tight-coupling approximation, the scattering rate is much larger than and , so the quadrupole and higher moments are suppressed by the mean-free-path parameter, with the quadrupole typically of order times the dipole. To leading order only monopole and dipole are needed. This suppression uses the collision term restored before decoupling; it does not follow from the collisionless equation by itself. Finite quadrupole corrections generate polarization and diffusion damping.
Use the notation throughout, with optical depth increasing inward and 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 gives
Integrating from the surface to a deeper boundary yields
For a semi-infinite atmosphere with no exponentially growing homogeneous term, the boundary term vanishes as . The formal solution of the radiative transfer equation is therefore
The 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 . Equivalently , 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 and integrate each term. Repeated integration by parts gives , so the polynomial source function moments for emergent intensity give
The limiting grazing-ray specific intensity is . 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 , the radial ray has the exact transfer relation
Expanding 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 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.
The high luminosity-to-mass ratio makes radiative acceleration and atmospheric departures from the simple approximations important. For an illustrative fully ionized hydrogen-rich composition with electron-scattering opacity , the stellar electron-scattering Eddington factor is
Thus electron scattering alone removes roughly half of the effective surface gravity. It is below the Eddington luminosity, so electron scattering alone does not prove that a hydrostatic atmosphere is impossible. Line opacity can provide substantially more radiative acceleration and can drive a stellar wind.
For a plane-parallel atmosphere, the geometrical condition is that the line-forming region be thin compared with the radius. In a simple isothermal photospheric estimate, the atmospheric scale height is
This is plane-parallel validity from pressure scale height. The mass and luminosity alone do not specify or the atmospheric temperature, so they cannot decide the geometry uniquely. For example, if and , the luminosity relation gives and . Plane-parallel geometry could then be adequate for deep weak photospheric lines. It would still fail for lines formed over an extended wind, and a cooler, more extended supergiant needs a separate geometrical assessment. This numerical example is conditional, not an additional datum about the stated star.
LTE requires collisions and local thermal processes to establish the relevant excitation and ionization populations at the local gas temperature. In the dilute outer layers of a luminous hot star, radiative transition rates can dominate collisions, and the radiation field arriving from other depths is not the local Planck function. A non-LTE stellar atmosphere should therefore use statistical equilibrium with radiative and collisional rates coupled to transfer. Deep layers may approach LTE, but surface ionization balances and abundance-sensitive spectral lines are particularly susceptible to departures. An LTE abundance can be systematically wrong even when a plane-parallel approximation is geometrically good.
A static atmosphere requires velocities to be negligible where the diagnostic lines form. Massive luminous stars commonly have line-driven outflows, so static models cannot represent wind-formed profiles or their mass density and velocity gradients. A hydrostatic approximation can still be useful below the sonic region if the chosen lines form there, provided radiative acceleration is included in the effective gravity. Large outflows call for spherical moving-atmosphere models rather than merely a changed abundance in a static model.
The three assumptions cannot be accepted on the quoted mass and luminosity alone. LTE is especially doubtful in the outer layers; wind-sensitive abundance diagnostics require non-LTE moving models, and geometry must be checked through atmospheric extension and line-formation depth. Reliable stellar chemical abundances should be supported by observed wind signatures, temperature and gravity constraints, and agreement between several suitable ionization stages or weak photospheric lines.
Take the optical depth to increase inward and let 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 is
This assumes negligible scattering, or a source function genuinely thermalized to ; LTE by itself does not turn an arbitrary scattering source into a Planck function. Multiply by , integrate between the two boundaries, and solve for the outward specific intensity:
The bottom boundary value is the quantity denoted 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 , the Planck function is constant and the integral can be evaluated:
For a bounded bottom intensity and , in every outgoing direction, and . 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.
In local thermodynamic equilibrium, thermal intensity samples the Planck function near an optical depth of order unity. The Eddington-Barbier relation makes this explicit: . 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, , whereas the dry adiabatic lapse rate has . Therefore
which 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 , , uniform disk-averaged thermal emission, a nearly blackbody stellar continuum, negligible reflected light at , and a flat response over the narrow band. If the eclipse depth is normalized to the out-of-eclipse star-plus-planet flux, the planet/star flux ratio is ; using changes the estimate by less than one percent. The exoplanet secondary eclipse then gives
Using the Planck law at gives , or about at the precision justified by the radius assumption. The Rayleigh-Jeans law instead gives , a useful rough check. This is a brightness temperature, not the bolometric planetary equilibrium temperature. With zero Bond albedo and global redistribution, the latter would be only about at ; a thermal window can sample hotter, deeper gas.
A continuum photon normally emerges from an optical depth of order unity. In hydrostatic equilibrium,
The emission-pressure degeneracy of a brightness temperature means that pressure cannot be uniquely recovered from the measured brightness temperature without and an opacity model. For an explicit nominal estimate, take and an effective continuum mass absorption coefficient . Then , or approximately one bar. In a hydrogen-rich atmosphere, collision-induced absorption by transient pairs supplies continuum opacity even without ordinary molecular lines; its density dependence makes the constant-opacity calculation an illustrative estimate. With a normal deep atmospheric pressure-temperature profile, temperature increases at larger pressures, eventually following an interior adiabatic temperature gradient.
Two mechanisms for are wavelength-dependent sampling of a vertically varying temperature, and wavelength-dependent weighting of a horizontally nonuniform dayside. In the first, carbon monoxide-rich opacity near can sample cooler gas around bar while the continuum samples hotter gas around one bar. In the second, hot and cool areas contribute different proportions of the total flux in different bands because the Planck function depends nonlinearly on temperature. The vertical mechanism alone produces the illustrated pair of temperatures; a single isothermal, horizontally uniform blackbody would have the same brightness temperature in both bands.
Figure 1.
Illustrative pressure-temperature profile with different infrared brightness-temperature contribution depths
. Details.
If the J-band observation is thermal emission, its low-opacity window normally samples at least as deeply as the continuum at . On the illustrated non-inverted profile, . The ordering is conditional: molecular-line absence alone does not specify the relative continuum opacities, and an inversion or reflected-light contamination can change the inference. The exoplanet emission spectrum measures weighted contribution regions rather than three exact thermometers at known pressures.
A natural interpretation of the short-wavelength peak is reflected starlight, while the longer-wavelength peak is planetary thermal radiation. A reflected spectrum approximately follows the stellar spectrum multiplied by the wavelength-dependent geometric albedo; thermal radiative flux approximately follows the planet's Planck function, modulated by molecular opacity. Thus two peaks need not represent two planetary surface temperatures.
For a quantitative estimate assume both are broad peaks, reflection has a slowly varying geometric albedo, and star and planet have approximately blackbody spectral envelopes. Wien's displacement law then gives
The stellar estimate is compatible at order of magnitude with an old, relatively small main-sequence star; equal age does not mean equal temperature to the Sun. Assume the far-infrared signal is thermal and sufficiently long-wavelength that the Rayleigh-Jeans law applies to both bodies. The planet-star radius estimate in the Rayleigh-Jeans limit gives
and therefore
This is a small volatile-rich-planet-sized estimate, not a Jupiter-sized one. It is conditional: peaks caused by molecular windows, strongly chromatic reflection, or a spectrum expressed as rather than do not support those two Wien estimates. Without , the far-infrared ratio only fixes .
For a close-in hot planet around a small star, transit-based atmospheric observations are the natural route if the orbital geometry allows them. Exoplanet transmission spectra gain from the small stellar radius, with a limb signal scaling as ; they probe composition at the day-night terminator. The stated thermal contrast also makes exoplanet secondary eclipse measurements a particularly useful route to the dayside exoplanet emission spectrum. Resolving such a close-in small planet by exoplanet direct imaging is much harder. Transit and secondary-eclipse spectroscopy are favoured for a transiting close-in interpretation; the supplied spectrum alone does not determine the orbital geometry or a unique best technique.