Limb darkening

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.

The linear limb-darkening law parametrizes angular specific intensity by a stipulated limb-darkening coefficient $u$:Here $I$ denotes specific intensity at one frequency, or radiance integrated over a specified interval of frequency or wavelength; it does not denote radiant intensity. Within this stipulated range, $u=0$ gives a uniform disk and $I(0)/I_c=1-u$. The linear limb-darkening law is linear in $\mu$, while its radial profile $I(r)=I_c[1-u+u\sqrt{1-(r/R_{\star }{)}^2}]$ is curved. In the linear source function example above, $u=b_{\nu }/(a_{\nu }+b_{\nu })$ when $a_{\nu },b_{\nu }\ge 0$ and $a_{\nu }+b_{\nu }>0$.

For the illustrative choice $u=0.6$, the specific intensity at $r/R_{\star }=0$, $0.8$, and the limiting edge is respectively $I_c$, $0.76I_c$, and $0.4I_c$. These are synthetic values, not measurements of the Sun. Real optical spectra and stellar atmospheres need not obey a linear limb-darkening law; more flexible profiles can be necessary for accurate exoplanet transit photometry. ^[ref]

The parameter $u$ in the linear limb-darkening law specifies its contrast between disk centre and the limiting edge. Its interpretation requires the chosen law and interval of frequency or wavelength. A limb-darkening coefficient fitted to an exoplanet transit may trade off against the transit impact parameter and planet radius; it is not a direct measurement of one layer's temperature.

For a circular stellar disk whose specific intensity is independent of azimuth, its mean specific intensity over projected area isThe polar coordinates element is $d{A}_{\mathrm{proj}}=2\pi rdr$ after angular integration, and $rdr=-R_{\star }^2\mu d\mu$. Substitution of the linear limb-darkening law givesFor a distant observer at distance $D\gg R_{\star }$ receiving light normally to the detecting plane, the observed radiative flux isIt is spectral irradiance, with units $\mathrm{W}{\mathrm{m}}^{-2}\mathrm{H}{\mathrm{z}}^{-1}$ when $I_{\nu }$ is per hertz. The local outward radiative flux per actual surface area is instead $F_{\nu }=2\pi {\int }_0^1{I}_{\nu }(\mu )\mu d\mu =\pi {\overline{I}}_{\nu }$, a spectral radiant exitance. In a spherical model with the same stellar atmosphere everywhere, the luminosity is $L_{\nu }=4\pi R_{\star }^2{F}_{\nu }$ and $f_{\nu }=L_{\nu }/(4\pi D^2)$, consistent with the inverse-square law. Compare disks at fixed $I_c$ or fixed total observed radiative flux explicitly: increasing $u$ lowers the total in the first comparison, while the second requires a larger $I_c$.

A dark, opaque circular planet of radius $R_p$ blocks specific intensity over its overlap with the projected star. With $k=R_p/R_{\star }$ and no transit dilution, the instantaneous exoplanet transit depth is exactlyWhen $k\ll 1$, the entire planetary disk is inside the stellar disk, and specific intensity varies little across that disk, small-planet transit weighting reduces toHere $z$ is the projected centre separation divided by $R_{\star }$. The numerator samples the local specific intensity; the denominator normalizes by the whole star's mean specific intensity. The small-size argument and the distinction between full overlap and ingress are explicit in the original analytic treatment. ^[ref]

For a short approximately straight crossing at constant transverse velocity $v_{\perp }$, write $x=v_{\perp }t/R_{\star }$ with $t=0$ at mid-transit, and let $b$ be the transit impact parameter. Then $z^2=x^2+b^2$ and ${\mu }_p=\sqrt{1-b^2-x^2}$. With stipulated $u=0.6$ and $k=0.05$, a central crossing has mid-transit exoplanet transit depth approximately $0.003125$, whereas $b=0.8$ gives $0.002375$. A uniform disk would give $k^2=0.0025$ for either fully overlapping crossing. The brighter central chord produces a deeper dip and the off-centre chord a shorter crossing at the same $v_{\perp }$.

At fixed $z<1$ away from the edge, a two-dimensional Taylor's theorem expansion over the symmetric planetary disk cancels the first derivatives; the first correction to the blocked radiant flux is of order $k^4$ if the second derivatives remain bounded. This error estimate is not uniform at the edge, where the radial derivative of the linear limb-darkening law diverges. Ingress, egress, a grazing exoplanet transit, and a crossing of a starspot require the overlap integral or a suitable more complete model. Inferring $k=\sqrt{\delta }$ without the local-to-mean specific intensity ratio can bias the estimated planet radius.

A motivating measurement used the McMath solar optical telescope at Kitt Peak. The Sun's image drifted across a $10\mathrm{mm}$-high, $0.2\mathrm{mm}$-wide optical slit; a double-pass spectrograph selected wavelength intervals from $0.0015\mathrm{nm}$ wide at $303\mathrm{nm}$ to $0.007\mathrm{nm}$ at $1099\mathrm{nm}$. An EMI 9750 QB photodetector or silicon photodiode supplied a digitized signal: $1\mathrm{kHz}$ samples were summed into 4096 positions per crossing, including dark and off-limb sky readings. ^[ref]

Controls included daily realignment of optical telescope mirrors, wavelength settings between spectral lines insensitive to the rotational Doppler shift, selection of detectors without detectable hysteresis, and measured electronic nonlinearity corrections below $0.1\%$ of maximum signal. Limb-position calibration combined centres of gravity of the edge derivative profiles with almanac semidiameters; the preliminary edge shifts were consistent with astronomical seeing of approximately $5$ arcseconds full width at half maximum. ^[ref]

The decreasing centre-to-limb readout gives a specific intensity profile, constraining stratification of the photosphere under a radiative transfer model. Astronomical seeing, instrumental scattering of light, and solar granulation limit this inference; brightness alone does not uniquely determine temperature versus geometric height. These observations need not follow the stipulated linear limb-darkening law.