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 20 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 A sight line stops where its own optical depth is one. The source function of a grey atmosphere in radiative equilibrium — a straight line, S = (3/4)(τ + 2/3) in units of the flux — with the depth each sight line reads off it. The emergent intensity at an angle is the source function integrated along the ray with everything in front of it attenuating, and for a source function linear in optical depth that integral is exactly the source function evaluated at τ = μ. Not approximately: the Eddington–Barbier relation is an identity for a linear source function, and the numerical integral agrees with it here to a part in ten thousand, which is the quadrature's error and not the relation's. A vertical ray reads the source function at τ = 1 and a ray at 78 degrees reads it at τ = 0.2, which is higher in the atmosphere and therefore cooler — so the limb is dimmer than the centre by the amount the source function has risen between those two depths. Everything about limb darkening is the slope of this one line, and the famous 2/3 is the intercept: the average ray, over the whole disc, reads the source function at two thirds of a unit of optical depth, which is the depth a stellar "surface" actually means.

A sight line stops at two thirds

The whole of limb darkening is one sentence: what leaves a star at an angle is the source function evaluated where the slant optical depth is one. For a source function linear in depth that statement is exact rather than approximate, and it produces the grey atmosphere's (2+3μ)/5 with no fitting anywhere in it.

starlight · Limb darkening A Cepheid's velocity, radius and angular size through two pulsation cycles. Three curves for a model Cepheid with δ Cephei's period of 5.37 days and a mean radius of 44 solar radii, at 272 parsecs, through two cycles. Top: the radial velocity of its spectral lines, in km/s, oscillating about the star's own velocity of -16.8 km/s — negative is towards the observer, so the surface is expanding when the curve is lowest. Middle: the change in radius, from integrating that velocity over time and multiplying by the projection factor p = 1.3: a swing of 4.9 solar radii, 11 per cent of the star. Bottom: the angular diameter the same star presents, in milliarcseconds, swinging by 168 microarcseconds about 1.505. The middle and bottom curves have the same shape because they are the same motion — one in kilometres, one in angle — and their ratio is the distance.

A distance that rests on one unmeasured number

A pulsating star's swelling can be measured twice — in kilometres from the Doppler shift of its lines, and in milliarcseconds by an interferometer — and the ratio is its distance, with no calibration against any other star. Except one: the factor that turns a line's velocity into the speed of the surface, which geometry puts at 1.33 to 1.5, measurement puts at 1.2 to 1.35, and every distance carries in proportion.

stars · Variable stars A colour that predicts an angle. The surface-brightness parameter Fᵥ — a combination of the V magnitude and the angular diameter that depends only on how bright each unit of the star's surface is — against the colour V − K, the ratio of the star's light at 0.55 and 2.2 microns. The line is the relation log θ = 0.5175 + 0.2753 (V − K) − 0.2 V, fitted to interferometric diameters of dwarfs and subgiants: Fᵥ = 3.962 − 0.1376 (V − K). The dots are 40 calibrating stars drawn about it with a 1 per cent scatter in angular diameter, the precision of a good interferometric measurement; refitting them returns a slope of −0.1374 and a scatter of 0.0023 in Fᵥ, which is 1.1 per cent in angle. The dashed curve is what blackbodies of 3,500 to 9,000 K would give, pinned to the Sun: it departs from the stars' line by up to 0.035 over the calibrated range, because real atmospheres have molecular bands and limb darkening a blackbody lacks — which is why the relation is calibrated on measured stars rather than derived. Given a star's V and V − K, the line returns its angular diameter without resolving it.

An angle too small to measure, read from a colour

A giant star in the Large Magellanic Cloud subtends two microarcseconds, a hundred times below anything an interferometer resolves. Its angle is read instead from its colour, through a relation calibrated on nearby stars that can be resolved — and because dust dims a star and reddens it by amounts that nearly cancel in that relation, the distance it gives hardly depends on the dust at all.

starlight · Angular diameter The edge of the Sun, computed. The brightness of the Sun at 500 nanometres across its limb, as a fraction of the brightness at the centre of the disc, against height above the level where the photosphere's radial optical depth is one, from the radiative transfer along each line of sight through a grey model atmosphere. Inside the disc the brightness falls gently with limb darkening; at the edge it falls to nothing over about two hundred kilometres. The steepest point of the fall — the inflection point, which is what a solar radius conventionally means — is 335 km above the radial τ = 1 level, because a ray grazing the limb passes through 270 times more atmosphere than a vertical one and becomes opaque that much higher up, at about 336 km. The brightness falls to one per cent 555 km above the τ = 1 level and to a thousandth at 694.

The Sun is three hundred kilometres larger seen edge-on

The radius of the Sun measured from images and eclipses and the radius measured from its oscillations disagree by about three hundred kilometres, and neither is wrong. A line of sight grazing the limb passes through hundreds of times more atmosphere than a vertical one and becomes opaque higher up — so the edge a telescope sees sits above the surface the waves feel, and every other definition of the edge, from the detector's threshold to the colour of the light, moves it again.

sky · Phases and eclipses A star brighter at its poles, and three paths a planet might take across it. A star rotating at 0.6 of its critical angular velocity, seen with its spin axis tipped 30° towards the observer from the plane of the sky, drawn with that axis vertical. Its outline is the Roche surface, 1.05 times as wide as it is tall in projection. The contours are the surface brightness, from von Zeipel's law with exponent 0.2 and a linear limb darkening: gravity at the equator is weaker, so the equator is cooler and its intensity there is 81 per cent of the pole's before limb darkening. Three transit chords at the same impact parameter, 0.3 polar radii, cross at projected spin–orbit angles of 0°, 45°, 90°. Only the first is symmetric about the star's axis; the others run from bright regions to dark ones or the reverse, and the light curve records which.

A transit that runs differently backwards

A star spinning near break-up is brighter at its poles than its equator, so a planet crossing it at an angle covers bright surface on the way in and dim surface on the way out. The light curve stops being symmetric in time, and the part that differs forwards and backwards measures the true spin–orbit angle from photometry alone — on exactly the stars where the line-profile method struggles.

exoplanets · Spin–orbit alignment A feature read back from the limb. A bump in the source function at optical depth 0.1, a quarter of a decade wide (black dashed), and what an inversion of the limb profile recovers on average when the profile is measured at 24 angles with errors of 1, 0.1, 0.01, 0.001 per cent, the inversion assuming only that the departure from a grey atmosphere is smooth. At one per cent the recovered bump has 56 per cent of its height and is spread over about twice its width; at a thousandth of a per cent, 82 per cent, and still wider than the truth. A thousandfold gain in precision buys a modest gain in resolution — each factor of ten a step, not a factor.

A profile that averages a decade of depth

Each point on a star's disc reports the temperature at a different depth, so a limb-darkening profile looks like a depth sounding of the atmosphere waiting to be read. It is one, but a blurred one. Every angle averages a whole decade of optical depth, the averages overlap, and recovering the temperature structure from them is the inversion of a Laplace transform — a problem whose resolution grows only as the logarithm of the data's precision.

starlight · Limb darkening The same two lines, seen as frequencies. The amplitude of the Fourier transform of the two line profiles of the previous comparison, on a logarithmic scale, against frequency in cycles per kilometre per second. At low frequency — the overall width — they agree, as the profiles did. At higher frequency they separate completely: the first star's transform falls to zero at 0.110 cycles per km/s and the second's at 0.132, because the rotation kernel's transform vanishes at 0.660/v sin i and no Gaussian-like broadening has a zero anywhere. The horizontal dashed line marks a noise level of 0.0005 in the transform, which a high signal-to-noise spectrum averaged over several lines reaches; the dip and the side lobe beyond it rise above it for both stars, so the rotation of each is measured by where its transform dips, independent of the macroturbulence. In a noisier spectrum the second star's weaker side lobe would be lost first.

A rotation found where the transform goes dark

A star's rotation and the turbulent flows on its surface both widen its spectral lines, and a profile fitted for one can be refitted for the other with almost no change. Convolution becomes multiplication in Fourier space, though, and the transform of a rotating disc has zeros — frequencies where it is exactly nothing — while the transforms of turbulence and of the spectrograph have none. Where the transform of a line goes dark is a measurement of rotation that turbulence cannot fake.

starlight · Line formation A star polarised only at its edge. The degree of linear polarisation of light from an atmosphere where electron scattering is the only opacity, against projected distance from the disc's centre in stellar radii (solid), and the relative brightness of the disc (dashed, scaled so the centre is 12). The polarisation is zero at the centre, 0.2 per cent halfway out, 1.9 per cent at nine-tenths of the radius and 11.7 per cent at the limb, always parallel to the limb. The brightness falls to a third towards the edge, so the most polarised light is also the faintest.

A polarisation that a round star cancels

In 1946 Chandrasekhar worked out that light leaving a hot star's atmosphere, where scattering by free electrons dominates, is polarised parallel to the star's edge — not at all at the centre of the disc, 11.7 per cent at the very limb. Over a round, unresolved star the pattern adds up to exactly nothing. It took seventy years to see it in the light of a whole star, because something has to break the circle first: a companion crossing the disc, a planet, or a star spun so fast that its outline is no longer round.

starlight · Polarimetry

Named alongside it

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

Angular diameterEclipsing binaryEffective temperatureInterferometryOptical depthRadiative transferStellar radiusEddington barbierGravity darkeningImpact parameterInclinationLine-broadening

All concepts