Limb darkening

Limb darkening is an optical effect seen in stars (including the Sun) and planets, where the central part of the disk appears brighter than the edge, or limb. Its understanding offered early solar astronomers an opportunity to construct models with such gradients. This encouraged the development of the theory of radiative transfer.

A filtered image of the Sun in visible light, showing the limb-darkening effect as a dimmer luminosity towards the edge or limb of the solar disk. The image was taken during the 2012 transit of Venus (seen here as the dark spot at the upper right).

Basic theory

An idealized case of limb darkening. The outer boundary is the radius at which photons emitted from the star are no longer absorbed. L is a distance for which the optical depth is unity. High-temperature photons emitted at A will just barely escape from the star, as will the low-temperature photons emitted at B. This drawing is not to scale. For example, for the Sun, L would be only a few hundred km.

Optical depth, a measure of the opacity of an object or part of an object, combines with effective temperature gradients inside the star to produce limb darkening. The light seen is approximately the integral of all emission along the line of sight modulated by the optical depth to the viewer (i.e. 1/e times the emission at 1 optical depth, 1/e² times the emission at 2 optical depths, etc.). Near the center of the star, optical depth is effectively infinite, causing approximately constant brightness. However, the effective optical depth decreases with increasing radius due to lower gas density and a shorter line of sight distance through the star, producing a gradual dimming, until it becomes zero at the apparent edge of the star.

The effective temperature of the photosphere also decreases with increasing distance from the center of the star. The radiation emitted from a gas is approximately black-body radiation, the intensity of which is proportional to the fourth power of the temperature. Therefore, even in line of sight directions where the optical depth is effectively infinite, the emitted energy comes from cooler parts of the photosphere, resulting in less total energy reaching the viewer.

The temperature in the atmosphere of a star does not always decrease with increasing height. For certain spectral lines, the optical depth is greatest in regions of increasing temperature. In this scenario, the phenomenon of "limb brightening" is seen instead. In the Sun, the existence of a temperature minimum region means that limb brightening should start to dominate at far-infrared or radio wavelengths. Above the lower atmosphere, and well above the temperature-minimum region, the Sun is surrounded by the million-kelvin solar corona. For most wavelengths this region is optically thin, i.e. has small optical depth, and must, therefore, be limb-brightened if it is spherically symmetric.

Calculation of limb darkening

Limb darkening geometry. The star is centered at O  and has radius R . The observer is at point P  a distance r  from the center of the star, and is looking at point S  on the surface of the star. From the point of view of the observer, S  is at an angle θ from a line through the center of the star, and the edge or limb of the star is at angle Ω.

In the figure shown here, as long as the observer at point P is outside the stellar atmosphere, the intensity seen in the direction θ will be a function only of the angle of incidence ψ. This is most conveniently approximated as a polynomial in cos ψ: $${\frac {I(\psi )}{I(0)}}=\sum _{k=0}^{N}a_{k}\cos ^{k}\psi ,$$ where I(ψ) is the intensity seen at P along a line of sight forming angle ψ with respect to the stellar radius, and I(0) is the central intensity. In order that the ratio be unity for ψ = 0, we must have $$\sum _{k=0}^{N}a_{k}=1.$$

For example, for a Lambertian radiator (no limb darkening) we will have all a_k = 0 except a₁ = 1. As another example, for the Sun at 550 nanometres (5.5×10⁻⁷ m), the limb darkening is well expressed by N = 2 and $${\begin{aligned}a_{0}&=1-a_{1}-a_{2}=0.3,\\a_{1}&=0.93,\\a_{2}&=-0.23\end{aligned}}$$

The equation for limb darkening is sometimes more conveniently written as $${\frac {I(\psi )}{I(0)}}=1+\sum _{k=1}^{N}A_{k}(1-\cos \psi )^{k},$$ which now has N independent coefficients rather than N + 1 coefficients that must sum to unity.

The a_k constants can be related to the A_k constants. For N = 2, $${\begin{aligned}A_{1}&=-(a_{1}+2a_{2}),\\A_{2}&=a_{2}.\end{aligned}}$$

For the Sun at 550 nm, we then have $${\begin{aligned}A_{1}&=-0.47,\\A_{2}&=-0.23.\end{aligned}}$$

This model gives an intensity at the edge of the Sun's disk of only 30% of the intensity at the center of the disk.

We can convert these formulas to functions of θ by using the substitution $$\cos \psi ={\frac {\sqrt {\cos ^{2}\theta -\cos ^{2}\Omega }}{\sin \Omega }}={\sqrt {1-\left({\frac {\sin \theta }{\sin \Omega }}\right)^{2}}},$$ where Ω is the angle from the observer to the limb of the star. For small θ we have $$\cos \psi \approx {\sqrt {1-\left({\frac {\theta }{\sin \Omega }}\right)^{2}}}.$$

We see that the derivative of cos ψ is infinite at the edge.

The above approximation can be used to derive an analytic expression for the ratio of the mean intensity to the central intensity. The mean intensity I_m is the integral of the intensity over the disk of the star divided by the solid angle subtended by the disk: $$I_{m}={\frac {\int I(\psi )\,d\omega }{\int d\omega }},$$

where dω = sin θ dθ dφ is a solid angle element, and the integrals are over the disk: 0 ≤ φ ≤ 2π and 0 ≤ θ ≤ Ω. We may rewrite this as $$I_{m}={\frac {\int _{\cos \Omega }^{1}I(\psi )\,d\cos \theta }{\int _{\cos \Omega }^{1}d\cos \theta }}={\frac {\int _{\cos \Omega }^{1}I(\psi )\,d\cos \theta }{1-\cos \Omega }}.$$

Although this equation can be solved analytically, it is rather cumbersome. However, for an observer at infinite distance from the star, $d\cos \theta$ can be replaced by $\sin ^{2}\Omega \cos \psi \,d\cos \psi$, so we have $$I_{m}={\frac {\int _{0}^{1}I(\psi )\cos \psi \,d\cos \psi }{\int _{0}^{1}\cos \psi \,d\cos \psi }}=2\int _{0}^{1}I(\psi )\cos \psi \,d\cos \psi ,$$ which gives $${\frac {I_{m}}{I(0)}}=2\sum _{k=0}^{N}{\frac {a_{k}}{k+2}}.$$

For the Sun at 550 nm, this says that the average intensity is 80.5% of the intensity at the center.

References

• Billings, Donald E. (1966). A Guide to the Solar Corona. Academic Press, New York.
• Milne, E.A. (1921). "Radiative Equilibrium in the Outer Layers of a Star: the Temperature Distribution and the Law of Darkening". MNRAS. 81 (5): 361–375. Bibcode:1921MNRAS..81..361M. doi:10.1093/mnras/81.5.361 – via Zenodo.
• Minnaert, M. (1930). "On the Continuous Spectrum of the Corona and its Polarisation". Zeitschrift für Astrophysik. 1: 209. Bibcode:1930ZA......1..209M – via NASA/ADS.
• Neckel, H.; Labs, D. (1994). "Solar Limb Darkening 1986-1990". Solar Physics. 153 (1–2): 91–114. Bibcode:1994SoPh..153...91N. doi:10.1007/BF00712494. S2CID 119704650 – via NASA/ADS.
• van de Hulst; H. C. (1950). "The Electron Density of the Solar Corona". Bulletin of the Astronomical Institutes of the Netherlands. 11 (410): 135. Bibcode:1950BAN....11..135V.
• Mariska, John (1993). The Solar Transition Region. Cambridge University Press, Cambridge. ISBN 0521382610.
• Steiner, O. (2007). "Photospheric processes and magnetic flux tubes". AIP Conference Proceedings. 919: 74–121. arXiv:0709.0081. Bibcode:2007AIPC..919...74S. doi:10.1063/1.2756784. S2CID 16932214.