Astronomy & Celestial Mechanics
the sky, computed
Concept

The point-spread function — where it appears

The image an optical system makes of a point source, into which every real image is convolved. Its width sets the resolution; its asymmetry is the more dangerous property, because an anisotropic point-spread function stretches every object the same way and imitates a gravitational shear.

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

Resolution stops improving at 10 cm of aperture. Angular resolution against aperture at 500 nm, both logarithmic. The falling line is diffraction alone, 1.03 λ/D, which is what a telescope in vacuum delivers and has no floor. The curve is the same telescope under an atmosphere of Fried parameter r₀ = 10 cm, combining diffraction and seeing in quadrature: it follows the diffraction line while D < r₀ and then bends onto a plateau at 0.98 λ/r₀ = 1.01″. An amateur's 100 mm at 0.1 m would resolve 1.062″ above the air and delivers 1.47″ through it; a metre at 1 m would resolve 0.106″ above the air and delivers 1.02″ through it; the VLT at 8.2 m would resolve 0.013″ above the air and delivers 1.01″ through it; the ELT at 39 m would resolve 0.003″ above the air and delivers 1.01″ through it. At 39 m the atmosphere is costing a factor of 371: the aperture is 390 coherence lengths across and every one of the 152,100 patches it collects arrives with a phase of its own. What the extra aperture still buys is photons and speckles, and those two are what adaptive optics and speckle interferometry respectively spend to get the falling line back.

A ten-metre mirror that resolves like a ten-centimetre one

The atmosphere delivers a wavefront in patches about ten centimetres across, and an aperture larger than a patch collects patches rather than detail. Resolution stops improving at that size — and what the extra aperture keeps buying is photons and speckles, which is why there are two entirely different ways out.

sky · Seeing What the error bar is made of, on a 1 m in 60 s. The four contributions to a photometric error, against the brightness of the star, for a 1-metre aperture, a 60-second exposure and a sky of 21 magnitudes per square arcsecond. The star's own photons give a line of slope exactly 0.2 — σ ∝ N^−1/2 and N ∝ 10^−0.4m, so a magnitude of extra faintness costs a fifth of a magnitude of precision, and no instrument changes that. The sky and the read noise are fixed counts, so their lines have slope 0.4, twice as steep, and they overtake the star at V = 18.25 — that crossing is the faint limit of the night, and it moves when the Moon rises rather than when the telescope changes. Scintillation is flat, because the atmosphere modulates a bright star and a faint one by the same fraction: at 4.09e-4 relative it is 0.44 millimagnitudes here and it is what caps the bright end, up to about V = 10.8. Below all of them is the systematic floor at 0.3 millimagnitudes, which is flat-fielding and colour terms and does not integrate down at all.

The faint star is measured against a brighter sky

For anything at the edge of detection the dominant source of noise is not the object. It is the sky in the same aperture, which is brighter than the star and is subtracted rather than measured — and once that is true, every rule of thumb about apertures, exposure times and image quality changes.

starlight · Photon noise A coherent one-per-cent distortion, invisible on every galaxy in the picture. 150 background galaxies behind a lens of Einstein radius 14″, each drawn at its own ellipticity: an intrinsic shape with a dispersion of 0.3 per component, plus the reduced shear the lens adds. The strongest shear on any galaxy here is 0.035, one part in 8 of the intrinsic scatter, so no object in this field is measurably distorted and the tangential alignment cannot be seen by eye at all. Averaged over these 150, the mean tangential ellipticity is 0.0088 ± 0.0245 against the 0.0130 the lens model predicts — consistent with the lens and equally consistent with nothing, because 150 galaxies buy a precision of 0.024 and the signal is 0.013. Detecting it at five sigma takes about 13,275 of them, which is not a picture anybody can draw. The cross component — every shape rotated by 45°, which gravitational lensing cannot produce — averages −0.0016 ± 0.0245, consistent with nothing, and that null is what separates a mass from a badly figured optic. The signal is not in any galaxy. It is in the sum, and the whole design of a lensing survey follows from that.

A one-per-cent distortion, and a million galaxies to see it

A galaxy's own shape is unknown and scatters with a dispersion of about 0.3, so a coherent one-per-cent shear is thirty times smaller than the noise on any single measurement. Nothing is ever measured about one object; the estimator is an average, and the whole design of a survey follows from 0.3 over the square root of N.

galaxies · Weak lensing A star drawn out into 3.0 arcseconds of spectrum. Atmospheric refraction relative to its value at 550 nanometres, against wavelength, at four zenith angles. The air's refractive index rises towards the blue, so the blue image of a star sits above the red one and the object is smeared into a short vertical spectrum. At 60 degrees from the zenith the separation across an optical band is 2.96 arcseconds — several times the size of the image at a good site, and comparable to the width of a spectrograph slit. Every curve here is the same curve multiplied by the tangent of the zenith angle, which is why one corrector with an adjustable strength works at every airmass. The practical consequences are three: a slit aligned other than vertically loses blue light or red light depending on where it was centred, a photometric aperture contains a different fraction of the light in each band, and an astrometric position depends on the colour of the star it is measured from.

The atmosphere is a prism as well as a lens

Refraction lifts a star towards the zenith, and everybody corrects for that. It lifts blue light further than red, and the difference is a short vertical spectrum a few arcseconds long — larger than the image, larger than a spectrograph's slit, and quietly present in every ground-based measurement not taken straight overhead.

sky · Refraction A planet that subtracts 73 per cent of itself at the inner working angle. The fraction of a planet's flux that survives an angular differential imaging subtraction, against its separation from the star in resolution elements, for a sequence covering 25 degrees of field rotation. The reference image is built from the target's own frames, so a planet that has not moved far between them is present in the reference and is removed along with the speckles. How far it moves is the arc length, which is proportional to the separation — so the self-subtraction is severe close in and negligible far out, and the half-throughput point is at 2.1 resolution elements here. The consequence for any contrast curve is that it is a statement about an algorithm as well as about an instrument: the depth reached has to be measured by injecting fake planets into the data and recovering them, because no calculation predicts what fraction of a real one survives.

A star subtracted using the star

Imaging a planet means removing a halo of scattered starlight a hundred million times brighter than the planet, and no model of that halo is good enough to subtract. So it is built from the star's own exposures — and since the planet is in those exposures too, it subtracts part of itself.

exoplanets · Direct imaging A bright star wants a wide aperture and a faint one wants 0.68 of the seeing. Signal-to-noise of simple aperture photometry against the aperture radius, in units of the seeing's full width at half maximum (1″), each divided by what optimal pixel weighting achieves for the same star, for stars of V = 12, 17, 20, 23 observed for 60 s through a 1 m telescope under a sky of 21 mag/arcsec². A small aperture loses starlight; a large one admits sky, and the balance depends on which dominates. For a bright star its own photons are most of the noise, so a wider aperture keeps gaining light almost for free and the best radius is large — 1.63 FWHM at V = 12, reaching 100.0 per cent of the optimum. For a star fainter than its sky the best radius shrinks to 0.680 FWHM and the best aperture reaches only 90.5 per cent of what weighting each pixel by its share of starlight divided by its variance achieves. That residual is exact in the background-limited limit: the best aperture captures 71.5 per cent of the light and 0.902 of the optimal signal-to-noise, so optimal weighting is worth 11 per cent in signal-to-noise, or 23 per cent in exposure time, and no more. The image is taken to be Gaussian; a real point-spread function has broader wings, which makes a fixed aperture a little worse and the optimal weights harder to know.

The best aperture throws away a tenth

Aperture photometry counts every pixel inside a circle equally and every pixel outside it not at all. For a faint star against its sky the best circle is two-thirds of the seeing wide, catches 71.5 per cent of the light, and reaches 90.2 per cent of the signal-to-noise that weighting each pixel by what it is worth achieves — a loss of 23 per cent in exposure time that no algorithm can beat by more.

starlight · Photon noise A sky wrong by a thousandth, and the disc that ends or does not. The surface-brightness profile of an exponential disc with a central surface brightness of 21.65 magnitudes per square arcsecond, against radius in units of its scale length, drawn straight — a pure exponential is a straight line in magnitudes — and then as it would be measured after subtracting a sky of 21.5 magnitudes per square arcsecond wrong by +0.1%, −0.1%, +0.3%, −0.3%. A sky left in bends the profile up into a shallower outer part, the shape called an anti-truncation; a sky taken out too hard bends it down into a break and then drops it out of existence where the residual goes negative. The +0.1% error departs by 0.2 magnitudes at 5.2 scale lengths, 27.3 magnitudes per square arcsecond; The −0.1% error departs by 0.2 magnitudes at 5.0 scale lengths, 27.1 magnitudes per square arcsecond; The +0.3% error departs by 0.2 magnitudes at 4.1 scale lengths, 26.1 magnitudes per square arcsecond; The −0.3% error departs by 0.2 magnitudes at 3.9 scale lengths, 25.9 magnitudes per square arcsecond. Real discs are observed to break both ways, at 24 to 27 magnitudes; every one of those breaks lies in the range where a sky error of a few tenths of a per cent would make one.

Where the sky is the measurement

The outskirts of a galaxy at 29 or 30 magnitudes per square arcsecond are a thousand times fainter than the night sky they are seen through. Photon noise averages down with exposure; an error in the sky subtracted from under them does not. A sky wrong by a tenth of a per cent draws a break into a disc that has none, or a halo around a galaxy that has none, and the telescope's own scattered light draws another.

galaxies · Surface brightness A disc that bends up into a halo only the telescope made. The surface brightness of a face-on exponential disc against radius in units of its scale length, in magnitudes per square arcsecond with the centre at 20, as it is (black, a straight line) and as a telescope sees it if 2% or 5% of each point's light is spread into faint wings falling as the inverse square of distance out to ten arcminutes. Scattered light from the bright centre lands in the outskirts, and once the disc has faded below it the profile stops falling and bends up into a shallow component that looks like a stellar halo: the observed profile rises three-quarters of a magnitude above the disc beyond 10.3 and 9.3 scale lengths, near 29.2 magnitudes per square arcsecond. At ten scale lengths the disc itself is at 30.9 and the image at 30.2 and 29.6.

A halo the telescope paints on every galaxy

A telescope puts almost all of a star's light into a dot and a few per cent of it into faint wings that reach arcminutes away. Around a galaxy those wings add up. The scattered light from the bright centre lands in the outskirts, and below about twenty-nine magnitudes per square arcsecond it outshines the galaxy's own — bending the profile up into a halo, or an edge-on disc out into a thick one, at the same depth around every galaxy whatever its size.

starlight · Photometric systems

Named alongside it

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

Sky backgroundAperture photometryBackground limitedExponential discPhoton noiseRead noiseScattered lightSeeingStellar haloSurface brightnessSystematic errorAdaptive optics

All concepts