Limb Darkening - Calculation of Limb Darkening

Calculation of Limb Darkening

In the figure on the right, 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 \, \textrm{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 ak=0 except a0=1. As another example, for the sun at 550 nm, the limb darkening is well expressed by N=2 and

(See Cox, 2000). Note - 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 which must sum to unity.

We can convert from ψ to θ using the relationship:


\cos(\psi) =
\frac{\sqrt{\cos^2(\theta)-\cos^2(\Omega)}}{\sin(\Omega)}

where Ω is the angle from the observer to the limb of the star.

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 Im is the integral of the intensity over the disk of the star divided by the solid angle subtended by the disk:

where dω=sin(θ)dθdφ is a solid angle element and the integrals are over the disk: 0≤φ≤2π and 0≤θ≤Ω. Although this equation can be solved analytically, it is rather cumbersome. However, for an observer at infinite distance from the star, the above equation simplifies to:

Read more about this topic:  Limb Darkening

Famous quotes containing the words calculation of, calculation, limb and/or darkening:

    “To my thinking” boomed the Professor, begging the question as usual, “the greatest triumph of the human mind was the calculation of Neptune from the observed vagaries of the orbit of Uranus.”
    “And yours,” said the P.B.
    Samuel Beckett (1906–1989)

    “To my thinking” boomed the Professor, begging the question as usual, “the greatest triumph of the human mind was the calculation of Neptune from the observed vagaries of the orbit of Uranus.”
    “And yours,” said the P.B.
    Samuel Beckett (1906–1989)

    You really were a panther, a wild-cat,
    who tore me limb from limb;
    my thanks for that.
    Hilda Doolittle (1886–1961)

    The only asylum
    Was the poorhouse, and those who could afford,
    Rather than send their folks to such a place,
    Kept them at home; and it does seem more human.
    But it’s not so: the place is the asylum.
    There they have every means proper to do with,
    And you aren’t darkening other people’s lives....
    Robert Frost (1874–1963)