Astronomy & Celestial Mechanics
the sky, computed
The collection

Every essay — page 30

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 581–600 of 909.

Gravitation

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

The delay a companion's gravity puts in a pulsar's pulses, over one orbit. The Shapiro delay of a pulsar's pulses against orbital phase, for a companion of 0.5 solar masses and orbital inclinations of 89.17°, 85°, 70°, from Δ = −2r ln(1 − sin i · sin Φ), where r is the companion's mass times G/c³ = 2.46 µs. Superior conjunction, when the pulsar is behind its companion, is at 90°. For the orbit nearly edge-on (89.17°, the inclination measured for PSR J1614−2230) the delay rises in a narrow spike 48.5 µs high as the pulses graze the companion; at 85° the spike is 30.8 µs and broad, at 70° 17.1 µs and nearly a sinusoid. The height is set by the companion's mass and the sharpness by the inclination, which is how one curve returns two numbers the orbit alone cannot give.

The part of a delay no orbit can imitate

A binary pulsar's pulses are delayed by up to fifty microseconds as they pass its companion, and most of that delay is invisible — a timing fit folds it into the orbit. What survives is the third orbital harmonic and above, whose size goes as the cube of the tangent of half the inclination, and that cube, not the delay, decides which neutron stars can be weighed.

5 figures · Relativistic orbits A cold sphere collapsing onto the virial theorem. The virial ratio 2K/|W| — twice the kinetic energy over the magnitude of the potential energy — counted over the bound shells, of 400 concentric shells released almost at rest (2K/|W| = 0.1) from a sphere whose density falls as one over the radius, against time in units of the sphere's free-fall time. It starts near zero; overshoots to 1.7 as the collapse reaches the centre, the inner shells first; rebounds and oscillates; and settles to 1.03, the virial value, staying within ten per cent of it from about 4.9 free-fall times on. The shells never collide with anything; each feels only the mass inside it. What damps the oscillation is phase mixing — shells of different energy oscillate with different periods and drift out of step — and what redistributes the energy is the potential changing while they move through it. A system of stars reaches equilibrium this way in a few dynamical times, where encounters between the stars would take far longer than the age of the universe.

An equilibrium reached without a single encounter

Every dynamical mass of a galaxy assumes its stars obey the virial theorem, and no two of those stars have ever exchanged energy. A collapse computed shell by shell shows how a system gets there anyway — in a few free-fall times, at half the size it fell from, by a potential that changes while its orbits are still moving through it.

5 figures · Virial theorem

Orbits

Kepler's three laws, and the family of curves a single force allows.

The eccentricity the planets force on an asteroid, and where it diverges. The forced eccentricity of a body in the inner Solar System against its semi-major axis, from the linear secular theory with all eight planets — its two giant-planet frequencies set to the measured 4.257 and 28.245 arcseconds a year — on a logarithmic axis: the largest value it reaches over the planetary cycles (solid) and its present value (dashed). Across most of the asteroid belt it is a few hundredths — 0.048 at 2.8 AU — the part of every asteroid's eccentricity that belongs to Jupiter and Saturn rather than to the asteroid. Near 2.04 AU it diverges: there the asteroid's own apsidal precession rate equals the frequency of the planetary mode dominated by Saturn, the denominator A − g₆ vanishes, and the theory's forced eccentricity goes to infinity. That is the ν₆ secular resonance. Real asteroids near it do not reach infinite eccentricity; they reach eccentricities large enough to cross the orbits of Mars and the Earth, and are removed. The location is for a body of low eccentricity and inclination; the full theory curves the resonance to larger distances at higher inclination, which is why the belt's inner edge is not a single number.

A denominator that emptied the inner belt

The eccentricity Jupiter and Saturn force on an asteroid is a fraction, and its denominator is the difference between the asteroid's own precession rate and one of the planets'. At two astronomical units that difference is zero. Where the zero sits depends on where Saturn is — and the belt's quiet orbits say Saturn did not take long to get there.

5 figures · Secular theory Perihelia of a young family, run backwards until they meet. The longitude of perihelion of each of 26 members of a synthetic young asteroid family at the Karin family's position, relative to the family's mean at each moment, against time before the present. Today the angles are spread round the whole circle: members only a few thousandths of an astronomical unit apart precess at rates that differ by a fraction of an arcsecond a year, which over millions of years is many full turns. Run backwards — with each member's own drift in semi-major axis included — they converge, and the spread is smallest at 5.75 million years ago, where they were made. The nodes turn the other way, and in the real family, whose rates come from a full theory in which the two frequencies differ, they give a second, independent convergence at the same instant. An age read this way is a measurement of a time in the way a light curve's period is, with an error of a few per cent, rather than an inference from a drift rate.

An age read from angles that meet

A family of asteroids a few million years old has had its members' orbits turned round and round by the planets, each at a slightly different rate, until their orientations look random. Turn them back and they come together at one instant — the collision — with an uncertainty of a few tens of thousands of years, but only if each member's own slow drift in sunlight is put back too.

5 figures · Asteroid families

The observed sky

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

How often Easter falls on each date, over the Gregorian computus's whole cycle. The fraction of years in which Gregorian Easter falls on each of its 35 possible dates, 22 March to 25 April, computed year by year over 5.7 million years — the length after which the computus repeats exactly, because the golden number, the weekday and the two century corrections only then return to the same combination. The distribution rises from 0.48 per cent on 22 March through a plateau near 3.3 per cent a date to a spike of 3.87 per cent on 19 April, then falls to 0.74 per cent on 25 April. The ramps at each end come from adding a Sunday to a full moon: the earliest dates need both the earliest full moon and the right weekday. The spike on 19 April has no astronomical meaning at all; it is made by a rule in the table that forbids the paschal full moon to fall on 19 April and moves it to the 18th, so the 18th is a full-moon date half as often again as its neighbours.

A full moon computed from the year number

The date of Easter is defined by the Moon and computed without looking at it — ten lines of integer arithmetic on the year number, a finite machine that repeats itself exactly every 5.7 million years. Its output has a spike on 19 April that no sky produces, a table that has kept pace with the real Moon for four centuries, and a Julian twin that is now four days behind it.

5 figures · Calendars How sharp a single short exposure is, for three sizes of telescope. The distribution of the Strehl ratio — the peak brightness of a star's image as a fraction of a perfect telescope's — for 600 independent short exposures through simulated Kolmogorov turbulence, with image motion removed, for telescopes 4, 7, 10 times the Fried parameter r₀ across. The screens' mean phase variance matches the analytic value for tilt-removed turbulence to within 2 per cent. At D/r₀ = 4 the exposures spread over the whole range and 31 per cent are better than 0.37 — a phase error under one radian, a nearly diffraction-limited image. At 7 the distribution has slid to low Strehl ratios and none of the 600 are that good; at 10, none. The best frames at the larger ratios are still several times sharper than the average, which is what practical selection keeps. The atmosphere is not uniformly bad: every so often, for a few milliseconds, the turbulence over a small enough aperture happens to be nearly flat, and those instants are a diffraction-limited telescope for free.

A sharp image kept by throwing the rest away

The atmosphere does not blur every instant equally. Now and then, for a few milliseconds, the air over a small telescope is almost flat, and a short exposure taken then is as sharp as the mirror allows. How often that happens falls as the exponential of the square of the telescope's size — which is why the method works on a two-metre mirror in the red and nowhere else.

5 figures · Seeing

Stars

The diagram that sorted them, and the one quantity that decides a star's whole life.

How long a star lives on the main sequence, down to the smallest stars. Main-sequence lifetime against mass, on logarithmic axes, as the hydrogen a star burns divided by its luminosity, scaled so that the Sun lives ten billion years. The dashed line assumes every star burns a tenth of its hydrogen, as the Sun does; the solid line lets stars below about a third of a solar mass, which are convective throughout, burn nearly all of it. The horizontal line is the age of the universe, 13.8 billion years: every star lighter than 0.90 solar masses that has ever formed is still on the main sequence. A star of a tenth of a solar mass lives about 8 trillion years — 597 times the present age of the universe — and its life so far is a few thousandths of one per cent of it.

A lifetime that is still a forecast

No star lighter than nine-tenths of the Sun has ever died, because none has had time. The smallest are convective to their centres, burn almost all their hydrogen rather than a tenth of it, and will shine for trillions of years. Most of the stars in the Galaxy therefore have a life history no observation has tested, and only one in thirteen stars ever born has finished.

5 figures · Stellar lifetimes 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.

5 figures · Variable stars

Spaceflight

Celestial mechanics used forwards — where to burn, and what it costs.

What it costs to reach Mars faster than the Hohmann transfer. The total velocity change to go from a 400 km orbit round the Earth to a 400 km circular orbit round Mars (solid), and the heliocentric part alone (dashed), against flight time, for transfers that leave the Earth's orbit tangentially and cross Mars's before reaching their far point; both planets are on circular, coplanar orbits. The right-hand end is the Hohmann transfer: 259 days and 5.65 km/s. Shortening it to 200 days costs 0.82 km/s more; to 150 days, 3.4 km/s more; below a hundred days the cost climbs past twenty. The curve leaves the Hohmann point with zero slope — the cheapest transfer is a minimum, so the first days saved cost almost nothing — and turns steeply upward only once the transfer ellipse has become long and thin. For the fastest transfers the orbit-to-orbit cost is below the heliocentric one, because a burn made deep in a planet's gravity well adds more to the speed at which the craft leaves it than the burn itself.

The first weeks off a Hohmann transfer are nearly free

The Hohmann transfer to Mars takes 259 days and is the cheapest route there is. Because it is a minimum, the cost of going faster starts flat — ten days saved costs twenty metres a second — and only climbs steeply once weeks have become months. For a crew, the time is itself a cost, paid in radiation at about two millisieverts a day.

5 figures · Orbital transfer How far a small error at departure moves the arrival at Mars. The distance by which a craft on the Hohmann transfer from the Earth's orbit to Mars's misses its aim point at the nominal arrival time, against an error in its departure speed, on logarithmic axes, integrated numerically for the 259-day flight. An error along the direction of motion (solid) moves the arrival by 76 km for every millimetre a second; one across it (dashed) by 23 km. The relation is linear over the whole range drawn: a metre a second at departure — the precision of a good launch — is 76 thousand kilometres at Mars, several times the planet's own radius and a thousand times what a landing needs. The ratio is the corresponding element of the state transition matrix, and it is the number every navigation plan is built around.

A millimetre a second at launch, seventy-six kilometres at Mars

A craft on its way to Mars arrives seventy-six kilometres from where it was aimed for every millimetre a second its departure speed was wrong. No launch is that good, so no trajectory is flown as designed; it is flown as corrected, by Newton's method run on the matrix that says how the arrival depends on the departure — which doubles its correct digits at every step and makes a two-body sketch into a real trajectory in two iterations.

5 figures · Patched conics Two clocks for a radiation-belt electron, at Kp 2. Two times against drift shell L, on a logarithmic scale, for electrons of about 1 MeV. The first is how long radial diffusion takes to move an electron by one Earth radius, the inverse of an empirical diffusion coefficient that grows as the tenth power of L — at Kp 2, 550 years at L = 2, 10 years at L = 3, 21 days at L = 5. The second is how long the electron survives before it is scattered into the loss cone. Inside the plasmapause, here at L = 4.7, that is set by plasmaspheric hiss: years at L = 1.3, 14 days at L = 2.5, 5.2 days at L = 3; outside it a single lifetime of 30 days stands for every other loss. Inside L = 4.8 the electron is lost before diffusion can carry it further in, which is the shaded region; outside, diffusion wins and the belt is fed. The two curves cross at a single place, and that crossing is the inner edge of the outer belt.

The gap between the belts is a race between two clocks

Electrons reach the inner magnetosphere by diffusing inward from the outer edge, and they are lost there to waves that scatter them into the atmosphere. The diffusion slows as the tenth power of distance and the loss does not, so the two clocks cross at one sharp place. Inside it lies the slot, and at megaelectronvolt energies the inner electron belt is not a steady population at all but what the last great storm left behind.

6 figures · Radiation belts

Starlight

Distance, brightness and colour: what can be measured when nothing can be visited.

The evidence for a period, one trial period at a time. The natural logarithm of the Bayes factor for a sinusoid at each trial period against no signal, for 120 radial velocities spread over 359 days with 2.0 m/s white noise per point and a planet of semi-amplitude 1.0 m/s at 11.3 days. The amplitude prior is Gaussian with scale 3.0 m/s in each component. At the planet's period the evidence is 14.4 — odds of about 1.8·10⁶ to one — and almost everywhere else it is near zero or negative, because a sinusoid that fits nothing is penalised for the amplitude range it wasted. The quantity that answers the question "is there a planet at some period" is not the peak but the average of the Bayes factor over the prior on the period, here uniform in the logarithm of the period from 1.5 to 1,000 days, and that average is ln B = 7.7: 6.7 lower than the peak, the Bayesian form of the look-elsewhere penalty, paid because the data had to find the period as well as the amplitude.

Odds for a planet that depend on a guess about its size

A false-alarm probability asks how often noise alone would produce a peak this tall. A Bayes factor asks something better — how much more probable the data are with a planet than without one — and pays for it with two numbers nobody measured. The width allowed for the period barely matters. The width allowed for the amplitude can turn strong evidence into none, from the same velocities.

6 figures · Periodograms 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.

6 figures · Angular diameter Two lobes, one halo, and two different paths out. A double radio source drawn in the plane containing the line of sight, the observer far to the right. The host galaxy sits at the centre of a hot halo — the rings mark 20, 40, 60, 100 kpc, in a β-model of core radius 40 kpc — and its two lobes lie 60 kpc out along an axis 45° from the line of sight. Light from the near lobe leaves the halo through its outskirts; light from the far lobe has to cross the halo's middle. With a tangled field of 3 µG at the centre in cells of 5 kpc, the rotation measure in front of the near lobe varies across the beam with a dispersion of 9 rad m⁻² and in front of the far lobe 28 rad m⁻². On the sky the two lobes are the same distance from the core, 42 kpc, and nothing in a total-intensity image says which is which.

The fainter polarisation belongs to the lobe further away

The two lobes of a radio galaxy lie at the same distance from the core on the sky, and nothing in a picture of them says which is nearer. Their polarisation does. Light from the far lobe crosses the magnetised hot halo of the host on its way out and light from the near lobe mostly does not, so at long wavelengths the far lobe loses its polarisation first — and the one-sided jet, it turned out, points at the lobe that keeps it.

6 figures · Faraday rotation

Gravitation

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

The pull of a ring, a shell and a point, in the plane. Radial force in the plane of a thin ring of radius 1, against distance from its centre, in units where the same mass at a point would pull with one over the distance squared; positive is outward. Inside a spherical shell of the same radius and mass the force is exactly zero. Inside the ring it is not: it points outward, toward the nearer part of the ring, at 0.34 at half the radius and 2.75 at nine tenths, growing without limit toward the rim. Outside, the ring pulls inward harder than a point would, by 25 per cent at twice its radius and 9 per cent at three times, because the near side of the ring is closer than the centre. The shell is indistinguishable from the point everywhere outside it.

Inside a ring the pull points outward

A hollow sphere exerts no force on anything inside it, and the proof is a pair of cones whose areas grow exactly as fast as gravity weakens. Cut the sphere down to a ring and the arcs grow only as the distance, the nearer side wins, and a planet set in the middle of a rigid ring falls into it within a handful of orbits — which is why Saturn's rings cannot be solid, and why a flat galaxy spins faster than its own enclosed mass allows.

6 figures · Shell theorem The same energy given to a cluster two ways. The distribution of the specific energies of 20,000 stars in a Plummer cluster (units G = M = a = 1), before and after an impulsive compressive shock of strength k = 0.06, which changes each star's vertical velocity by k times its height. Stars whose energy ends above the potential at the tidal radius, −0.196, marked by the vertical line, are lost. The thin curve is the cluster before. The dashed curve adds only the half-square term, the heating the textbook integral computes: every star gains a little, and 0.8 per cent of the stars cross the line. The solid curve adds the whole change, including the cross term between each star's own velocity and its kick, which averages to zero — the mean energy input is 1.59e-3 against 1.52e-3 — but broadens the distribution: 1.6 per cent of the stars cross. The same energy, 2.2 times as many escapers.

Stars lost to a kick that adds no energy

A tidal shock changes each star's energy by two terms. One is the heating every textbook integral computes, always positive and second order in the shock. The other is the product of the star's own velocity and its kick, first order in the shock and exactly zero on average — and for the weak shocks a cluster meets most often it decides who leaves, doubling the losses without adding any energy at all.

6 figures · Impulse approximation

Galaxies

Where the unknown stops being a number and becomes a profile — and most of it is not light.

Polarised dust, synchrotron and the signal between them. The amplitude of B-mode polarisation at angular multipole 80, in microkelvin of CMB temperature on a logarithmic scale, against frequency from 20 to 450 GHz, for a patch of sky as clean as the cleanest large field observed. Dust is a modified blackbody with emissivity index 1.53 at 19.6 K, normalised to 4.6 µK² at 353 GHz; synchrotron falls as frequency to the power −3.1 in antenna temperature, normalised to an illustrative 1 µK² at 23 GHz. The horizontal lines are a primordial signal with tensor-to-scalar ratio 0.01, about 0.048 µK, and the lensing B-modes, about 0.071 µK. The foregrounds together are faintest at 81 GHz, at 0.038 µK, 0.8 times the signal. At 150 GHz the dust alone is 0.105 µK, the amplitude a tensor ratio of 0.046 would produce; at 353 GHz it is 2.23 µK, 47 times the signal. No frequency has the signal above the foreground.

A foreground drawn by the field it traces

The grains that polarise starlight along the Galaxy's magnetic field glow, at millimetre wavelengths, with light polarised across it — and that glow is brighter than the primordial signal cosmologists are looking for at every frequency above eighty gigahertz and in every patch of sky. Its pattern on the sky is set by how loosely the interstellar filaments follow the field, and removing it rests on a bet that the clouds along each line of sight share one temperature.

6 figures · Dust polarisation The matter power spectrum with galaxies' feedback taken out. The ratio of the matter power spectrum with baryonic feedback to the spectrum without it, against wavenumber in h per megaparsec on a logarithmic scale, for three strengths of feedback that span the range the large hydrodynamic simulations disagree over. Gas blown out of haloes by supernovae and active nuclei removes mass from scales below a few megaparsecs, so the ratio falls below one from about 0.5 h per megaparsec and is deepest near 4 h per megaparsec, at 9, 17 and 26 per cent; at still smaller scales gas that has cooled into galaxy centres adds mass back. The marks show the wavenumber that dominates a lensing measurement at multipoles 100, 1000, 5000, taken where the lensing kernel for sources at redshift 1 peaks: 0.09, 0.87, 4.35 h per megaparsec. The shape is a parametric stand-in placed where simulations put the effect, not a simulation.

A measurement of dark matter that needs a model of galaxies

Weak lensing weighs all the matter, dark and ordinary, by how it bends light. But the ordinary matter does not sit where gravity alone would put it — galaxies' supernovae and black holes blow a large part of it out of their haloes — and the scales where that happens are the scales that carry most of a lensing survey's statistical power. Ignored, it reads as a less clumpy universe, in the direction of the disagreement lensing has with the microwave background.

6 figures · Weak lensing The separation of the Milky Way and Andromeda, run forward. The distance between the Milky Way and Andromeda against time from now, in billions of years, integrated from today's separation of 770 kpc and radial approach of 110 km/s, for four values of Andromeda's sideways (tangential) velocity: 0, 17, 57, 100 km/s. The haloes are extended Hernquist spheres of 1.15 and 1.5 trillion solar masses, and dynamical friction drains the orbit at each passage. At 0 km/s the first approach is so nearly head on that the centres pass within 10 kpc and the pair merges at 4.0 Gyr. At 17 km/s the first approach is so nearly head on that the centres pass within 10 kpc and the pair merges at 4.8 Gyr. At 57 km/s the first pericentre is 117 kpc at 4.5 Gyr, and the pair merges at 9.2 Gyr. At 100 km/s the first pericentre is 321 kpc at 4.8 Gyr, and the pair has not merged by 14 Gyr. The measured sideways speed runs from about 17 km/s (Hubble) to about 57 km/s (Gaia), and the outcomes across that range differ by billions of years.

A collision decided by a speed across the sky

Andromeda is coming toward the Milky Way at a hundred and ten kilometres a second, and that part is measured to a few per cent. Whether the two galaxies merge in the next ten billion years depends on how fast Andromeda is also moving sideways — a drift of fifteen millionths of an arcsecond a year, measured twice, with answers of seventeen and fifty-seven kilometres a second. One is a near head-on collision in under five billion years, the other a coin toss.

6 figures · Galaxy interactions

Stars

The diagram that sorted them, and the one quantity that decides a star's whole life.

Plasma beta from the corona's base outward. The ratio β of gas pressure to magnetic pressure against distance from the Sun's centre, one to ten solar radii on logarithmic scales, for a static isothermal corona in a dipole field. quiet Sun, 5 G at the base, with 10⁹ electrons per cubic centimetre at 1.5 million kelvin: β starts at 0.41 and crosses one again at 3.46 solar radii; quiet Sun, 10 G at the base, with 10⁹ electrons per cubic centimetre at 1.5 million kelvin: β starts at 0.1 and crosses one again at 5.00 solar radii; active region, 30 G at the base, with 3·10⁹ electrons per cubic centimetre at 2.5 million kelvin: β starts at 0.058 and crosses one again at 2.97 solar radii. The gas pressure falls steeply at first and then levels off, because gravity weakens with distance; the field's pressure, falling as the sixth power, never levels. Every corona therefore has a second place where the gas takes over, and it sits at two to four solar radii — where the closed loops of the corona end and the streamers' tips begin.

A corona held by its field and opened by its gas

A thousand kilometres above the Sun's surface the gas pressure has fallen so far that the magnetic field takes charge, and the corona's loops are the shape the field gives it. But the gas pressure of a hot atmosphere stops falling a few solar radii out, while the field's keeps falling as the sixth power of distance. So there is a second place where the gas takes over — and it is where the loops end, the streamers begin, and the solar wind is let out.

6 figures · Plasma beta

Every sequence · Every named object · Search