Astronomy & Celestial Mechanics
the sky, computed
The collection

Every essay — page 4

One idea per essay, ordered so that the earlier ones set up the later ones — but nothing here depends on being read in sequence. Essays 61–80 of 514.

Gravitation

Two bodies pulling on each other, and everything that goes wrong at three.

7 decades of relaxation time, and a Hubble time between them. Crossing time and two-body relaxation time for 5 self-gravitating systems, against the number of bodies in each. The lower marks are the crossing time R/σ, which spans 3 decades; the upper ones are t_relax ≈ 0.1N/ln N times it, which spans 8. The horizontal rule is a Hubble time, and the same expression puts these systems on opposite sides of it. Two of them relax — an open cluster and a globular cluster — so their stars have exchanged energy, their heavy stars have sunk, and their present structure is not the one they were born with. The rest do not: a dwarf spheroidal, an elliptical galaxy, a cluster of galaxies, the elliptical missing it by a factor of 10⁶, which is why a galaxy is modelled as a smooth potential with no stars in it at all. The dependence is on N and hardly at all on anything else, because the Coulomb logarithm makes every decade of impact parameter contribute equally. A cluster that relaxes also evaporates: about a stellar mass in a hundred and forty leaves per relaxation time, so the globular cluster empties in 6.1·10¹¹ years.

How long a system takes to forget

A star crossing a galaxy is deflected by every other star, and the deflections add as a random walk. The time for that walk to change a star's energy by its own amount is a hundred million billion years for a galaxy and a billion for a globular cluster, and everything about how the two are modelled follows from which side of a Hubble time that falls on.

9 figures · Two-body relaxation Io's measured heat needs k₂/Q = 0.016, and Enceladus's needs 0.011. Tidal surface heat flux against orbital eccentricity, from Ė = (21/2)(k₂/Q)GM_p²R⁵ne²/a⁶ evaluated at each satellite's own orbit, drawn at a common k₂/Q of 0.015. Every line has slope 2 because the dissipation is quadratic in e and nothing else on this axis varies. The filled marks are each body at its actual eccentricity; the two ringed ones are the bodies with a measured surface heat flux, and they are the only points here that are observations. Solving each of those for k₂/Q gives 0.016 for Io and 0.011 for Enceladus — within a factor of 1.5 of one another, for a warm silicate body and a 500-kilometre ball of ice, which ought to be a coincidence and is instead the sharpest problem in the subject: nothing about Enceladus's ice can plausibly dissipate at 0.011, and the number is what the heat requires all the same. The dashed line is the Earth's measured surface heat flux, 0.087 W/m², which Io exceeds by a factor of 28 — the most volcanically active body in the solar system is the fourth largest moon of the fifth planet, and the reason is entirely in the orbit.

A moon heated by not being allowed to relax

Tidal dissipation goes as the square of an eccentricity that tides themselves destroy, so a moon radiating tidal heat is spending something it cannot have saved. Io's would be gone in a hundred and forty thousand years, and the resonance that keeps putting it back is the reason there are volcanoes.

9 figures · Tidal heating GW150914: 33 Hz to 250 Hz in 0.22 seconds. The strain of GW150914 — two black holes — through the last 0.22 seconds before merger, computed from the quadrupole sweep at the chirp mass its fit returned, 28.716 solar masses, and drawn at the luminosity distance it returned, 440 megaparsecs. Two things rise together and neither is free to rise on its own: the frequency goes from 33 Hz to 250 Hz, and the envelope — the outer curve — grows by a factor of 3.9, because the amplitude goes as f^2/3 and nothing else in it changes over so short a span. The vertical axis is in units of 10⁻²¹, so the peak here is a fractional length change of about 2.9·10⁻²¹: over the four kilometres of an interferometer arm that is 1.2·10⁻¹⁷ metres, a thousandth of the width of a proton. The chirp mass is not fitted to the amplitude at all — it comes from the spacing of these zero crossings, which is why it is the best-determined number in the whole event and why the distance, which does come from the amplitude, is the worst.

A distance with no ladder under it

The frequency sweep of an inspiral fixes the chirp mass with no distance in it, and the amplitude then gives the luminosity distance directly, because one expression fixes both. That is a distance measured with nothing calibrated beneath it — and its error budget is one angle.

8 figures · Gravitational waves Nothing visible in any pulsar, and a quadrupole in the angle between them. Above: 4 millisecond pulsars' timing residuals over 15 years, at the few hundred nanoseconds a good one reaches. Each wanders, and none of them shows anything a reader could call a signal; a gravitational-wave background of amplitude 2.4·10⁻¹⁵ at one cycle per year contributes a common part to all of them that is smaller than each pulsar's own red noise. Below: the correlation between pairs, against the angle on the sky between them. 2211 pairs out of 67 pulsars, binned into 15 angles, against three curves with no free parameters between them. A quadrupolar background gives the Hellings–Downs shape — positive for nearby pulsars, negative near 83°, and back up to exactly half its zero-separation value at 180° because a background looks the same in opposite directions. An error in the observatory clock would give a flat line, because it shifts every pulsar identically. An error in the solar-system ephemeris would give a cosine, because it moves the barycentre in one direction. The drawn points prefer the quadrupole over the flat line by Δχ² = 358. That is the detection: not a waveform, not an event, not a moment — a shape in an angle, accumulated over fifteen years, on data taken for another purpose entirely.

A detector the size of the galaxy

At a nanohertz no instrument can be built, so the clocks already in the sky are used instead. The signal is in no single pulsar's data — it is in the correlation between pairs as a function of the angle between them, and that curve has no free parameters at all.

8 figures · Gravitational waves The adaptive step starts 5.7× more accurate and ends 7.5× worse. The envelope of the relative energy error against elapsed revolutions — the worst error within each plotted interval rather than the error at one sample in it — for the same second-order symplectic integrator run two ways on a Kepler orbit at e = 0.5. The flat band is a fixed step of one 200th of a period: its error oscillates once per revolution and the envelope does not grow, because a symplectic method at constant h is the exact solution of a nearby Hamiltonian and is conserving that one. The rising curve is the same method with the step refined where the orbit is fast — h ∝ r^3/2, the free-fall time, varying by a factor of 5 around the orbit — which is the first thing anybody reaches for at a close encounter and which is strictly more accurate step for step: over its first three revolutions it stays a factor of 5.7 below the fixed run. By 1500 revolutions it is a factor of 7.5 above it and still climbing. Changing the step changes which Hamiltonian is being conserved, the errors from successive steps stop cancelling, and what is left is a random walk with no bound at all. The practical consequence is that a solar-system integration cannot adapt its step: it either keeps a step short enough for the closest encounter it will ever meet, or it detects the encounter and hands that piece of the trajectory to an entirely different, non-symplectic method for the duration — which is what every long-term integration of the planets actually does.

A step that must not be adapted

A symplectic integrator's bounded energy error is a property of a fixed step. It is conserving a Hamiltonian a step-size away from the intended one, and changing the step changes which Hamiltonian — so refining the step at a close encounter, which is the first thing anybody does, destroys the only property the method was chosen for.

8 figures · Numerical integration A wave whose crests count out 35 kilograms per square metre of ring. Optical depth across a spiral density wave in Saturn's A ring, drawn outwards from the 5:3 inner Lindblad resonance with Mimas at 131,988 km. The wave is launched at the resonance on the left and damps away to the right. Its wavelength is not constant: it starts at about 5.6 km and has shortened to 1.13 km by the last crest drawn, because the wavenumber grows in proportion to the distance from resonance and the ring's own self-gravity is the only restoring force in the dispersion relation. That makes the pattern a chirp with exactly one unknown in it. Fitting the 24 crest positions actually drawn here — the square of each one's distance from resonance against its number, which is a straight line — returns a surface density of 35.0 kilograms per square metre against the 35 the profile was built from. The rings have been weighed this way rather than by anything touching them: the mass per unit area follows from counting bright bands in a light curve as a star sets behind the ring.

A ring weighed by the wave crossing it

Saturn's rings are a few tens of metres thick and spread over an area larger than the Earth, made of pieces nobody can resolve, and nothing has ever landed on them. Their mass per unit area is nevertheless known to a few per cent — from the rate at which the crests of a wave crowd together as it travels outwards.

9 figures · Planetary rings One measured number, and every pair of masses that produces it. The plane of the two component masses of an inspiralling binary, with three curves of constant chirp mass across it. The middle one is GW150914's value of 28.7 solar masses, and the chirp mass recovered from the coordinates of the drawn curve varies along its whole length by 7.4e-14 per cent — which is the point: every binary on that line radiates the same frequency sweep at leading order, so the early inspiral cannot tell them apart. Two of them are marked. An equal pair of 33.0 and 33.0 solar masses and a lopsided pair of 63.6 and 18.4 sit on the same contour, and their total masses differ by a factor of 1.24. What separates them is the mass ratio, which enters the phasing only at the first post-Newtonian order, suppressed by the square of the orbital speed in units of the speed of light — small through the hundreds of cycles that carry most of the signal, and appreciable only in the last few, where that speed approaches a third of c. So the chirp mass is a measurement and the individual masses are an inference from the end of the signal, which is exactly the part a detector's high-frequency noise eats first.

One number where two masses were

The hundreds of orbits an inspiralling binary completes inside a detector's band depend on its two masses only through one combination of them. Every pair on that contour radiates an identical sweep, so the early signal — which carries nearly all the signal-to-noise — cannot say which pair it was.

8 figures · Gravitational waves 7 decades of relaxation time, and a Hubble time between them. Crossing time and two-body relaxation time for 5 self-gravitating systems, against the number of bodies in each. The lower marks are the crossing time R/σ, which spans 3 decades; the upper ones are t_relax ≈ 0.1N/ln N times it, which spans 8. The horizontal rule is a Hubble time, and the same expression puts these systems on opposite sides of it. Two of them relax — an open cluster and a globular cluster — so their stars have exchanged energy, their heavy stars have sunk, and their present structure is not the one they were born with. The rest do not: a dwarf spheroidal, an elliptical galaxy, a cluster of galaxies, the elliptical missing it by a factor of 10⁶, which is why a galaxy is modelled as a smooth potential with no stars in it at all. The dependence is on N and hardly at all on anything else, because the Coulomb logarithm makes every decade of impact parameter contribute equally. A cluster that relaxes also evaporates: about a stellar mass in a hundred and forty leaves per relaxation time, so the globular cluster empties in 6.1·10¹¹ years.

A cluster that boils itself away

A star cluster has no thermostat. Encounters between its members push a few of them above escape speed, the cluster loses them, and losing them makes it contract — which makes it hotter, which makes more of them escape. A self-gravitating system heats up as it loses energy, and the process ends by destroying the system.

8 figures · Two-body relaxation

The observed sky

Coordinates, seasons, phases and shadows — geometry seen from inside it.

The sky from latitude 52°. The celestial sphere seen from latitude 52 degrees. The pole stands 52 degrees above the horizon, the celestial equator meets the horizon due east and west, and a star at declination 20 degrees traces the drawn circle once a day. Everything below the horizon is drawn faint.

The sphere that is not there, and why it is still the right model

The stars are at wildly different distances and the celestial sphere is a fiction. It is also the most useful fiction in observational astronomy, because for pointing at things, distance is exactly the information to throw away.

9 figures · Celestial sphere The Sun's altitude through the day at latitude 52°. Solar altitude against the hour of the day, at one latitude, for the solstices and the equinox. Where a curve crosses zero is sunrise or sunset, and the width between the crossings is the length of the day.

The Sun's path, and the tilt that makes the seasons

Summer is not when the Earth is closest to the Sun — that happens in January. It is when the Sun climbs higher and stays up longer, and both come from a 23.4° tilt.

8 figures · Seasons Phases are a viewing angle, not a shadow. A satellite at eight points of its orbit. Exactly half of it is lit at every one of them; what changes is how much of the lit half faces the centre. Nothing is in shadow except at an eclipse.

Phases are not shadows, and eclipses are

Half the Moon is lit at every instant of every month. The phase is which part of the lit half faces the Earth — and confusing that with a shadow is the commonest error in astronomy.

9 figures · Phases and eclipses The path of the celestial pole over 25,772 years. The circle the Earth's rotation axis traces among the stars, at a radius equal to the obliquity, with the bright stars that fall near it and the years at which each is closest. Polaris is the pole star for a few centuries either side of now, and nothing else on the circle is nearly as close.

The pole star has a shelf life, and the sky has a slow hand

The Earth's axis traces a circle among the stars once every 25,772 years. Polaris is at the pole now, was not four thousand years ago, and will not be in two thousand more.

9 figures · Precession The equation of time, and its two causes. The difference between a sundial and a clock over the year, in minutes, computed from Kepler's equation and the tilt. The eccentricity term has one cycle a year and the obliquity term has two; their sum runs from −14.2 to 16.4 minutes.

The Sun is a bad clock, by up to sixteen minutes

Solar noon and twelve o'clock are not the same instant, and the discrepancy runs through a fixed annual cycle. It has two causes, one from the shape of the orbit and one from the tilt of the axis.

9 figures · Equation of time One star, two coordinate systems, at latitude 52°. The equatorial grid and the horizon grid drawn on the same sphere for an observer at latitude 52°. The star marked has declination 20° and hour angle -40° in the first, and altitude 45.5° and azimuth 239.4° in the second. The two frames differ by a single rotation through the co-latitude 38°, which is why the celestial pole stands 52° above the northern horizon.

Where a star is depends on who is asking

The sky needs two coordinate systems because two different things stay still in it — the observer's horizon and the stars themselves. One rotation converts between them, and the angle of that rotation is the time.

8 figures · Celestial sphere Three periods that nearly share a multiple. How far the draconic and anomalistic months are from a whole number, after a whole number of synodic months, in hours. At 223 synodic months — 6585.321 days — both are within an hour of closing, which is what makes an eclipse repeat. The draconic residual is 0.87 hours and the anomalistic 5.19 hours.

The eclipse that repeats a third of a world away

Three lunar periods nearly share a multiple after 6,585 days. The word "nearly" is what makes eclipses predictable, and the leftover third of a day is what moves each repeat a third of the way round the Earth.

10 figures · Phases and eclipses Refraction against apparent altitude. How far the air lifts an object, in arcminutes, against where the object appears to be. The lift is 34.5′ at the horizon — larger than the Sun's own diameter — and falls below one arcminute above 45°. At the horizon it is also the least reliable number in positional astronomy, because it depends on the temperature profile of the air the sight line crosses.

The Sun sets before it sets

At the moment the Sun's lower edge appears to touch the horizon, the whole of it is already below. Refraction lifts it by more than its own diameter, squashes it while it is there, and makes every sunrise and sunset time in every almanac a statement about the air rather than about the sky.

10 figures · Refraction The sub-Earth point over 400 days. Where on the Moon the Earth stands overhead, in selenographic longitude and latitude, sampled over 400 days. The excursion in longitude is the equation of centre — the Moon's rotation is uniform and its orbital motion is not, so the face runs alternately ahead and behind by up to 6.3° — and the excursion in latitude is the tilt of its equator, up to 6.7°. The two run on months of different length, so the track never repeats.

The face that is not quite fixed

The Moon keeps one face turned toward the Earth, and the sentence is exactly true only of a fictitious Moon on a circular orbit. The real one rocks by a few degrees each month, in two directions and for two unrelated reasons, and the rocking has shown 59 per cent of its surface to people who never left the ground.

9 figures · Libration Twilight at latitude 52°. Solar altitude through the second half of the day at 52°, on June solstice, equinox, December solstice, with the four thresholds that define twilight marked. Sunset is at −0.833° rather than at 0° because the Sun's own semi-diameter and the atmosphere's refraction together lift it by that much when it is geometrically already down. Civil twilight ends at −6°, when the brightest stars appear and outdoor work stops; nautical at −12°, when the horizon can no longer be seen against the sky and a sextant becomes useless; astronomical at −18°, when the sky stops contributing to a photometric measurement. Their durations at equinox here are 34 min, 40 min, 42 min — twilight is not a fixed length, it is the reciprocal of how steeply the Sun descends, and the Sun descends at an angle of about 90° − φ to the horizon. At this latitude the curve for June solstice never reaches −18° at all: astronomical twilight does not end, and there is no astronomically dark night from about 48.6° upward.

Three definitions of night

Twilight ends at three different depression angles, and each threshold is a statement about what can no longer be done. Its length is not a duration but a rate — how fast the Sun goes down — and above one latitude the deepest of the three never arrives at all.

8 figures · Twilight Mars through 8.3 months of sky. The geocentric ecliptic longitude and latitude of Mars over 252 days — 8.3 months — centred on opposition, computed as the direction of P − E with both orbits taken as circles: the Earth's of radius 1 AU, the planet's of 1.5237 AU inclined 1.850°. The motion reverses for 72.7 days and backs up 15.94° of longitude, and that interval is centred on opposition to within 0.0000° of longitude. The track closes on itself: over 149 days the planet visits the same point of the sky twice, and the loop it encloses is 15.9° long and 2.94° tall. Longitude and latitude are at different scales: a degree of latitude is drawn 4.6 times a degree of longitude, because the sweep of 51° in longitude and 4.5° in latitude will not share a scale on one page. Circular orbits mean one loop for every apparition, and that is the real cost of this figure: Mars is 0.525 AU away at the opposition drawn here, its true distance at opposition varies with its eccentricity, and the real loops differ in size from one apparition to the next because of it. The marks are 20 days apart, and they crowd where the motion stops.

The loop a planet does not make

Mars stops in the sky, backs up for ten weeks, and goes on. Neither orbit reverses anywhere — the loop belongs to the difference of two position vectors seen from one of them, and its width and duration are fixed by the ratio of the two radii and by nothing else.

8 figures · Apparent motion Why the solar day is 3m 55.91s longer than the sidereal one. Two positions of the Earth one day apart, with the direction of a fixed star and the direction of the Sun marked at each. After one sidereal day — 23h 56m 4.0905s — the Earth has turned through exactly 360° and the star is back on the meridian; the Sun is not, because the Earth has moved 0.9856° along its orbit, and the further 0.9856° of turning takes 3m 55.91s. That is the whole of the difference: 366.2422 turns against the stars in the 365.2422 solar days of a year, one more turn than the 365.2422 against the Sun. Nothing here is to scale. The orbital arc is drawn at 40°, which exaggerates the real 0.9856° by 41 times, and the Earth's disc is about 4073 times too large for its orbit. The two star sight lines are drawn parallel because they are: a star's distance cannot be put on the same page as an orbit.

The day that is four minutes short

A star crosses the meridian 3 minutes 55.91 seconds earlier each night, and the same star slips 3 minutes 56.56 seconds a day against a civil clock. Those are two different numbers, and the gap between them is the one extra turn the Earth makes against the stars in every year.

8 figures · Sidereal time

Every sequence · Every named object · Search