Astronomy & Celestial Mechanics
the sky, computed
Concept

Limb darkening — where it appears

The fall in surface brightness from the centre of a stellar disc to its edge, which measures the temperature gradient rather than any property of a surface. It must be modelled before a transit depth becomes a planetary radius, and the coefficients used are among the larger systematic uncertainties in an exoplanet radius.

Named by 12 essays across 4 fields — each of them below, with the objects they name alongside it.

A shadow edge 10.3 metres wide, and a stellar diameter read off how blurred it is. A star disappearing behind the Moon, drawn as intensity against position across the shadow. The horizontal axis is in Fresnel scales of √(λD/2) = 10.3 metres at 550 nm and 3.844e+5 km, which is the only length the problem has; at a limb speed of 0.62 km s⁻¹ one of them takes 16.6 milliseconds to pass, so the whole event is over in a tenth of a second and needs photometry at a kilohertz. The Moon has no atmosphere and its limb is a knife edge, and a knife edge does not cast a shadow with an edge: the intensity at the geometric boundary is 0.250, a quarter rather than a half, and outside it the light overshoots to 1.37 before ringing down. Every one of those numbers is a property of the wave and of nothing else. What the star contributes is the blurring. Each point of the stellar disc casts its own copy of the pattern, displaced by its own position, so the observed trace is the pattern convolved with the star's projected disc — 22.4 metres wide for the 12 milliarcsecond curve, against a 10.3-metre fringe. The contrast falls from 0.28 to 0.04 across the four curves drawn, and inverting that fall is how several hundred stellar diameters were measured with a single telescope, no interferometer, and no resolution at all. The picture cannot show the limitation that ended the technique's dominance: the Moon goes where it goes, so only stars within a few degrees of the ecliptic are ever occulted, and each is occulted at whatever position angle the geometry happens to offer.

The edge of a shadow is a wave

An asteroid's shadow has an edge because the asteroid is large. The Moon's does not — at visible wavelengths and lunar distance the edge of a shadow is ten metres wide, so a lunar occultation is a diffraction pattern sweeping past at half a kilometre a second, and how blurred its fringes are is the star's own diameter.

sky · Occultations Fringe visibility for a 47 mas disc at 575 nm. Fringe visibility against the separation of the two apertures, for a disc 47 milliarcseconds across seen at 575 nm. The solid curve is a uniform disc, |2J₁(x)/x| with x = πθB/λ; it is exactly one at zero baseline, where both apertures see the same wavefront, and falls to zero at 3.08 m — read off the drawn samples, and equal to 1.21967 λ/θ to better than one part in a million. That is the measurement: not a brightness, a baseline. The dashed curve is the same disc with linear limb darkening u = 0.4, whose null is 4.9% further out at 3.23 m — so the same observed null implies 47 mas as a uniform disc and 49.3 mas limb-darkened, and a diameter quoted without its model is a number without a unit. At the 2.54 m aperture of the telescope this was done on, the visibility is still 0.17: one mirror cannot reach the null, which is the same statement as saying it cannot resolve the star.

An angle of five hundredths of an arcsecond

No telescope has ever resolved a star other than the Sun, and stellar diameters are measured anyway — by finding the separation of two apertures at which the star's interference fringes vanish. What that returns is an angle; the radius arrives only when a distance is brought in, and the distance is the worse-known half.

starlight · Angular diameter Three profiles of equal equivalent width, 28.6 mÅ. Three absorption profiles with the same equivalent width — 28.6 milliångström, 1.72 km/s at 500 nm, matched to better than 0.1% by root-finding over the quadrature — differing in nothing but shape. Left: the cores, on a common velocity axis. Right: the same three normalised to their own half widths, on a logarithmic depth scale. The thermal profile is a Gaussian set by Fe's mass at 6000 K, 1.34 km/s; the collisional one a Lorentzian of γ = 2.28e-3 nm; the rotational one the classical kernel of a disc turning at v sin i = 1.73 km/s, which is exactly zero beyond 1.28 half widths and is the only one of the three with an edge. Matching the areas does not match the widths: the half widths are 1.11 km/s (thermal), 1.35 km/s (rotational), 0.68 km/s (collisional), a factor of 1.98 between the widest and the narrowest. At three half widths the collisional wing is 49 times the thermal one and at five it is 1.27·10⁶ times — six decades, which is why a line's shape stays diagnostic long after its width has stopped being so. 11% of the Lorentzian's own equivalent width lies beyond the right-hand panel's edge and is not drawn anywhere.

The same width for three different reasons

Thermal motion, rotation and collisions each widen an absorption line, and they can be tuned to areas that agree to a part in a thousand. What separates them is the shape, and the shape carries a rotation speed from one profile and a surface gravity from another.

starlight · Line formation The limb is 40% as bright as the centre, and that is a temperature gradient. Left: a stellar disc shaded by the grey-atmosphere law I(μ)/I(1) = (2 + 3μ)/5, in 26 steps, with μ = cos θ read from the centre outwards. Right: that law against μ, with linear laws at the measured solar coefficients from 400 to 1600 nm. The grey law's coefficient is exactly 3/5 and its very limb is exactly 2/5 of the central brightness — both read off the drawn curve rather than quoted — because the Eddington–Barbier relation makes the emergent intensity at angle μ the source function at optical depth τ = μ, and in radiative equilibrium that source function is linear in τ. The limb is not a cooler part of the star. A sight line entering at the edge reaches unit optical depth higher up, where the gas is cooler, so what the darkening measures is the run of temperature with depth; a star with an isothermal atmosphere would show a uniform disc, and one with a steeper gradient a darker limb. The measured coefficients fall from 0.9 at 400 nm to 0.35 at 1600 — the same gradient seen through a less steep Planck function — which is why a radius measured from a transit or a fringe null has to say which colour it was measured in.

The light that is missing from the edge

The Sun's limb is forty per cent as bright as its centre, and the reason is not that the edge is cooler. A sight line entering at the edge stops higher up, so what the darkening measures is the temperature gradient — and it is worth seventeen per cent on Betelgeuse's radius.

starlight · Limb darkening Two radial-velocity curves, and one mass ratio. The line-of-sight velocity of each star through one orbit of AI Phoenicis. Both curves are computed from the two masses and the period; what a spectrograph delivers is the reverse. The ratio of the amplitudes is the inverse ratio of the masses — 48.2 to 50.3 kilometres a second, so the heavier star moves more slowly — and the sum of the amplitudes with the period gives the mass sum, 2.437 solar masses, once the inclination is known from the eclipses.

The only stars whose masses are known

A star's mass cannot be measured by looking at it. It can be measured by watching two stars pull on each other, and if the pair also eclipses, the same observations give both radii as well — with no stellar model anywhere in the chain. A few hundred such systems calibrate everything else.

stars · Binary stars The light curve of AI Phoenicis, computed from its elements. Total light against orbital phase, computed by overlapping two discs of radius 1.805 and 2.9303 solar radii at an inclination of 88.5°, each weighted by its own surface brightness. The two eclipses hide the same area of sky and have different depths — 48.0 and 19.1 per cent — because what is lost is the light of whichever star is behind, and the ratio of the depths is therefore the ratio of the two surface brightnesses. Two things are left out and both matter to a real solution: this is the bolometric light rather than the light in a filter, and the discs are uniform, where a real one is limb-darkened and so has a deeper, rounder eclipse than the flat-bottomed one drawn here.

Two radii, from a light curve alone

The radius of a star is not measured. It is inferred, from a temperature and a luminosity, through a model. There is one exception — a pair of stars that eclipse each other, whose light curve and velocity curves between them give both radii, both masses and the ratio of temperatures with no model of a stellar interior anywhere in the chain.

stars · Binary stars A transit of a planet 0.103 of its star's radius. The star's brightness through one transit, computed by integrating the uniform stellar disc over the region the planet covers. The depth is 1.05%, which is exactly (Rp/R⋆)² = 0.01055. The four contact points are where the two discs are externally and internally tangent, at separations 1 ± 0.103 stellar radii.

A planet measured by the light it removes

A transit gives a depth, and the depth is a ratio of two radii rather than a size. Everything a transit says about a planet is said in units of a star nobody has visited either.

exoplanets · Transits A transit of a planet 0.103 of its star's radius. The star's brightness through one transit, computed by integrating the limb-darkened stellar disc over the region the planet covers. The depth is 1.26%, deeper than (Rp/R⋆)² = 0.01055 because the planet crosses a limb-darkened disc whose centre is brighter than its average. The four contact points are where the two discs are externally and internally tangent, at separations 1 ± 0.103 stellar radii.

Four contact points, and what they fix

The depth of a transit gives a radius ratio. The shape gives the impact parameter, and then — through nothing but Kepler's third law — the mean density of the star being crossed.

exoplanets · Transits A 27.5 m/s velocity the star does not have. Left: a rotating stellar disc, approaching on one side and receding on the other, with the chord a planet of 0.1 stellar radii takes across it at impact parameter 0.5 and a sky-projected obliquity of 0°. Right: the apparent radial velocity that results, computed by covering the disc cell by cell — the flux hidden at each phase and its mean line-of-sight velocity — rather than from a fitted formula. The star's centre of mass does not move at any point in this: the anomaly is entirely a statement about which parts of the line profile are missing. Its amplitude is 27.5 m/s at v sin i = 4.5 km/s, and the two numbers are related by the depth of the transit, since blocking a fraction f of light of mean velocity v shifts a flux-weighted centroid by f·v. The curve is antisymmetric about mid-transit to 0.00% of its own amplitude, which is what an aligned transit gives: equal time on the blue half and the red. Limb darkening is included at u = 0.6, and it matters: it weights the hidden light towards the centre of the disc, where the rotation velocity is smallest.

A velocity measured from a shape

A transiting planet hides part of a rotating disc, so the star's line profile loses a slice at one velocity and its fitted centroid moves. The star has not moved at all — and the lopsidedness of that motion is the whole measurement of whether the orbit lies in the star's own equatorial plane.

exoplanets · Spin–orbit alignment A transit depth of 1.200 per cent for a planet of area 1.000 per cent. Three transits of the same planet across the same star, differing only in how the star's brightness falls toward its edge. A planet of radius ratio 0.1 covers 1.000 per cent of the stellar disc's area, and if the disc were uniformly bright that would be the depth. It is not uniformly bright: a sight line near the limb leaves the photosphere at a shallow angle and therefore from a cooler layer, so the edge is dimmer than the centre, and a planet crossing near the middle blocks light that is brighter than average. The transit drawn with realistic coefficients is 1.200 per cent deep — 20 per cent deeper than the area — and it is also rounder, because the covered brightness changes through the crossing instead of staying flat. The consequence is stated in the numbers beside the curves. Each is a least-squares fit of the radius ratio to the realistic curve, performed with a different assumed limb-darkening law, and the recovered radius moves by up to 3.6 per cent depending on which law is assumed. Fitting with the law the curve was made from returns the input to five figures, which is the control: the bias is the mis-specification and not the fitter. Since the coefficients come from a model atmosphere rather than from the light curve, every published planetary radius carries a systematic from stellar physics that no amount of photometric precision removes — and it is the dominant one for the best-measured planets. The picture holds the impact parameter fixed; a grazing transit is worse, because it samples only the limb, where the disagreement between laws is largest.

The depth is not the area

A planet covering one per cent of its star's disc does not make a transit one per cent deep. The star is brighter in the middle, so a planet crossing the middle blocks more than its share — and the correction depends on coefficients that come from a stellar atmosphere model rather than from the light curve.

starlight · Limb darkening A star 1.24 times wider than it is tall, and 19 per cent brighter pole-on. Left, the meridional section of a star rotating at ω = 0.9257 of its critical angular velocity, computed from the Roche potential rather than sketched: the equator sits at 1.245 polar radii, and at the critical rate that ratio is exactly 1.5 whatever the star is made of. The same rotation expressed as a fraction of the critical equatorial speed is 0.768, and the two conventions differ by the distortion itself — a figure that prints one under the other's name is wrong by an amount that looks like rounding. Effective gravity at the equator is 0.329 of its polar value, so von Zeipel's flux law makes the pole hotter than the equator by a factor 1.320 at the theoretical exponent 0.25 and 1.232 at the 0.188 that interferometric imaging actually fits. Right, the apparent bolometric brightness against viewing inclination, integrated over the visible gravity-darkened surface: pole-on the star is 1.19 times brighter than edge-on, and the apparent temperature falls with it. The consequence is that a rapid rotator's place on the Hertzsprung–Russell diagram is partly a statement about the observer's position, which no spectrum taken alone can undo.

A temperature that depends on where the observer stands

A star turning near its break-up rate is half again as wide as it is tall, and its equator is thousands of degrees cooler than its poles. Neither of those is a small correction to a spectrum — the effective temperature and the luminosity such a star appears to have are partly statements about which way its axis happens to point.

starlight · Gravity darkening A few per cent of diameter, hidden in the second lobe. Visibility against baseline in the natural variable πθB/λ, for stars with linear limb-darkening coefficients of 0, 0.3, 0.6, 0.9, each rescaled to the uniform disc that best fits its own first lobe. In the first lobe the four curves are within 0.72 per cent of one another; past the first null they differ by up to 4.2 per cent. That is the whole difficulty of measuring a stellar diameter. A limb-darkened star has a faint edge, so a uniform-disc fit returns a diameter 9.7 per cent too small at u = 0.9 and 2.5 per cent too small at u = 0.3 — and the information needed to tell which is in a region where the visibility is under five per cent and the calibration errors of a real interferometer are comparable to the signal. The correction from what is measured to what is wanted is taken from a model atmosphere, because the observation that would supply it is the hardest one there is. The fit here is done by least squares on the drawn curves rather than read from a conversion table, and the zero-coefficient case returns the uniform disc to 0.00 per cent, which is the fit checking itself.

A diameter that depends on a model atmosphere

An interferometer measures fringe visibilities and somebody fits a disc. A uniform disc and a limb-darkened one agree to within a per cent across the whole of the first lobe and differ by five in the second — where the visibility is under five per cent and the calibration errors are the same size.

starlight · Angular diameter

Named alongside it

The objects these essays reach for when they reach for this one.

Angular diameterStellar radiusEclipsing binaryEffective temperatureImpact parameterInclinationRadius ratioTransitDiffraction limitInterferometryLine-broadeningOccultation

All concepts