Astronomy & Celestial Mechanics
the sky, computed
Concept

General relativity — where it appears

The theory of gravity that adds one further term to the Newtonian orbit, turning the apsidal line and draining the orbit's energy at a computed rate. Its effects in the solar system are small and unambiguous: 43 arcseconds a century of Mercury's precession, and a binary pulsar's orbit shrinking by 76 microseconds a year.

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

Mercury's perihelion, term by term. The observed advance of Mercury's perihelion is 5600 arcseconds per century against the equinox. Almost all of it is the equinox: the coordinate frame itself turns, and removing it leaves 574. Subtracting the perturbations of the other planets — computed by Le Verrier in 1859 and refined many times since — leaves 42.98 arcseconds a century that nothing in Newtonian gravitation accounts for. The bars are logarithmic in nothing; they are the real proportions, which is why the residual is barely visible beside the frame term.

Forty-three arcseconds, after everything else

Mercury's perihelion moves through 5,600 arcseconds a century. Nearly all of that is the coordinate system turning, and almost all of the rest is the other planets. What was left over was 43 — under one per cent of the raw number, and the most consequential residual in the history of the subject.

gravitation · Relativistic orbits PSR B1913+16: 45 years of periastron arriving 88 seconds early. The cumulative shift of periastron passage for PSR B1913+16, Δt = ½(Ṗ_b/P_b)T², over 45 years from its discovery in 1975. The orbital period is 0.322997449 days and is measured to be shortening at 2.423e-12 seconds per second — a change in the twelfth decimal place, which over a career accumulates to 87.5 seconds. That is the whole reason the effect is measurable. Ṗ_b is taken here as a measured quantity and no radiated power is derived from it; the parabola is the general-relativistic prediction as published, and the points are that prediction scaled by the published ratio of observed to predicted decay, 0.997 ± 0.002 (Hulse & Taylor 1975; Weisberg & Huckins 2016), which is how the agreement is quoted. The error bars are drawn: at the last point the bar is 1.25 px tall against a dot 9 pixels across, so they are invisible, and their invisibility is the result. The campaign has run 45 years.

An orbit measured to be shrinking

A 7.75-hour orbital period that shortens by 68 nanoseconds each turn is beyond any single measurement, and unmissable after fifty thousand of them, because the shift accumulates as the square of the elapsed time. Forty-five years of pulse arrival times have made it 87.5 seconds, which is why the rate is a measured quantity rather than an inferred one.

gravitation · Relativistic orbits Relativity switches the cycle off, halving its reach at a ratio of 0.80. The greatest eccentricity a Kozai–Lidov cycle reaches, against the strength of the orbit's own relativistic pericentre precession, measured in units of the cycle's own precession rate at zero eccentricity. The horizontal axis is logarithmic and spans three decades. At the left the relativistic term is negligible and the cycle reaches 0.838, which is the closed-form value for a start at 65 degrees and is what the integration is checked against. At the right it is gone. The mechanism depends on the pericentre staying put while the outer body pulls on the same side of the orbit for a whole cycle, and the relativistic precession is a competing rotation of that same pericentre; when it is faster, the pull averages away. The threshold sits near one by construction and the transition is sharp rather than gradual, with the reach halved at 0.80. What makes it matter is where the relativistic term is largest: it grows as the pericentre falls, so it strengthens exactly as the cycle drives the orbit inward, and it therefore sets a floor on the pericentre distance that this mechanism can deliver a body to.

The precession that switches the cycle off

A distant companion can trade an orbit's inclination for its eccentricity, over and over, and drive a pericentre almost onto the central body. General relativity's own precession competes with the mechanism for the same pericentre, and when it wins the cycle stops — sharply, at a ratio of one, which puts a floor on how close anything can be delivered.

orbits · Kozai–Lidov Central pressure bracketed without a model: 6 bodies, 23 decades apart. What can be said about the middle of a body from its mass and its radius alone. The lower end of each bar is GM²/8πR⁴, which follows from hydrostatic equilibrium and nothing else — no equation of state, no composition, no temperature, no assumption whatever about how the density is arranged inside. The upper end costs one more assumption, that the density does not increase outward, and it needs a central density, which is a model output rather than an observation and is why that edge is drawn as the softer one. The dot is what a full structural model gives. For the first five bodies every dot lies inside its bar, and what is worth noticing is how wide the bar is: Sun's rigorous floor is 4.48e+13 pascals against a modelled 2.34e+16, a factor of 522. The bound is true and nearly useless there, because most of a centrally condensed body's pressure comes from the concentration and the derivation deliberately knows nothing about it. The relativistic entry is the exception, and the reason to draw the figure at all. neutron star's modelled central pressure is 8.5 times the Newtonian ceiling — a body no Newtonian arrangement of matter with density falling outward can produce. The floor still holds, and holds for a statable reason: relativity makes the pressure gradient steeper than Newtonian gravity does, so the true central pressure can only exceed what the Newtonian derivation demands. The bracket therefore does more than constrain an interior. Applied at a small enough radius it breaks, and where it breaks is where Newtonian hydrostatics has stopped being the right equation.

A floor under the centre that assumes nothing

There is a lower bound on the pressure at the centre of any body in hydrostatic equilibrium, and it needs no equation of state, no composition and no temperature — only a mass and a radius. For the Sun it is nearly useless. For a neutron star it says which theory of gravity the interior needs.

stars · Hydrostatic equilibrium 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.

gravitation · Relativistic orbits One hot spot, seen through more and more gravity. The X-ray flux from a small hot spot at colatitude 60° on a rotating neutron star seen from inclination 60°, over one rotation, each curve divided by its own mean, for no bending, 2GM/Rc² = 0.20, 2GM/Rc² = 0.34, 2GM/Rc² = 0.48 — the last two are a 1.4 M☉ star of 12 km and a 2.0 M☉ star of 12.4 km. Light bending is computed exactly; Doppler effects are left out. With no bending the spot is behind the star for 39 per cent of the turn and the pulsed fraction is 1.00. As the compactness rises the star's gravity bends light from the far side into view: the dark interval shrinks to 27 per cent, then vanishes, and the pulsed fraction falls to 1.00, 0.99, 0.64. The more compact the star, the more of its own back it shows and the smoother its pulse.

The pulse that shows the back of the star

A hot spot on a spinning neutron star flashes once a turn, and the flash is weaker than it would be on an ordinary star because the star's gravity bends light from its far side into view. How much weaker measures the ratio of mass to radius. But on its own the flash cannot tell a compact star from a tilted one — the two trade exactly — and what separates them is the lean the spin puts into the pulse, which measures the radius directly.

stars · Degeneracy The Sun bends light from every direction in the sky. The gravitational deflection of starlight by the Sun, as seen from the Earth, against the angle between the star and the Sun, as general relativity predicts it (solid) and as a Newtonian particle of light would be bent (dashed), exactly half. At the Sun's limb the deflection is 1.75 arcseconds, the value measured at the eclipse of 1919. It falls only as the cotangent of half the angle, so it is 4.07 milliarcseconds at right angles to the Sun and 9.8 at 45°, the closest a space astrometry mission looks; it reaches zero only directly opposite the Sun. Every star in a modern catalogue has been displaced by an amount far larger than the precision of its position, and the shading marks the angles such a mission actually observes.

Every star in the catalogue is bent by the Sun

The bending of starlight by the Sun is remembered as an eclipse measurement of 1.75 arcseconds at the solar limb. The same bending reaches every direction in the sky, falling only slowly away from the Sun, and at right angles to it is still four milliarcseconds — two hundred times the precision of a modern star position. Every catalogue corrects every star for it, and the correction has become the most precise test of how much space a mass curves.

sky · Aberration Five clocks around the black hole at the centre of the Galaxy. Timescales against distance from a black hole of 4.3·10⁶ solar masses at the centre of a cluster whose stellar mass inside r is 10⁶ (r/1 pc)^1.5 solar masses, made of stars of 1 solar mass. Ordinary two-body relaxation (orange) is slow everywhere, 8.4·10¹⁰ years at 0.01 pc. Orbits that close are nearly Keplerian and keep their orientation until they precess — through the stars' own mass farther out, through relativity close in — and over that coherence time the torques between them add rather than random-walk. That makes scalar resonant relaxation (blue) faster than two-body relaxation inside about 0.57 pc, and vector resonant relaxation (green), which turns orbital planes and is not stopped by precession within a plane, faster still: 3.5·10⁷ years at 0.1 pc, against the 6·10⁶-year age of the young stars there. Coefficients of order one are set to one; the slopes are the content.

Torques that add because the orbits hold still

Near a massive black hole every star's orbit is almost a fixed ellipse, repeated thousands of times before it precesses. The torques the other stars exert on it therefore point the same way orbit after orbit, and they add in proportion to time instead of cancelling in a random walk — so angular momentum near the hole relaxes hundreds of times faster than ordinary encounters allow, and orbital planes faster still.

gravitation · Two-body relaxation Where relativity overtakes the margin: 306 Schwarzschild radii for 10⁵ solar masses. Two things against the size of a contracting star, measured in its own Schwarzschild radius, on logarithmic axes. The flat lines are the margin by which the star's pressure-weighted adiabatic exponent exceeds 4/3 — the Newtonian condition for stability — for stars of 1000, 10⁴, 10⁵, 10⁶, 10⁷, 10⁸ solar masses with mean molecular weight 0.61; it does not change as the star shrinks, because it is fixed by the gas fraction of the pressure, and it falls with mass from 0.036 to 1.2·10⁻⁴. The diagonal is what general relativity adds to the requirement, 1.1245 × 2GM/(Rc²), which in these units is the same line for every mass. Contraction moves each star leftward along its flat line, and where the line meets the diagonal the star becomes dynamically unstable: at 31, 97, 306, 966, 3,052, 9,649 Schwarzschild radii respectively. Every one of those crossings lies between 31 and 9,649 Schwarzschild radii out, where the star is not remotely compact — the relativistic correction is tiny, and it decides the matter only because the margin it is compared with is tinier.

A star tipped over at a thousand times its horizon

A star of a million suns is held up almost entirely by its own light, which leaves it stable by about a thousandth. General relativity raises the bar by an amount that is tiny until the star is compact, and a thousandth is tiny too — so the star collapses at nearly a thousand times the radius of the black hole it becomes, at a mean density below that of air, and with no fuel yet burning.

stars · Polytropes What DI Herculis's orbit should do, and what it does. The rate at which the periastron of the eclipsing binary DI Herculis advances, in degrees per century. Top: the prediction for stars whose spins are aligned with the orbit — 2.34 from general relativity, computed from the masses (5.17 and 4.52 solar masses), the 10.550164-day period and the eccentricity of 0.489; and 1.93 from the stars' tidal and rotational distortion, as published stellar models give it, of which the rotational flattening is 1.50 and the tides 0.43. Total 4.27. Bottom: the observed rate, between about 0.6 and 1.4 degrees per century in the published timing analyses — a quarter to a seventh of the prediction. For a quarter of a century after the discrepancy was established it was cited as possible evidence against general relativity, which contributes more than half the predicted total by itself.

A precession slowed by stars lying on their sides

The eclipsing binary DI Herculis turns the long axis of its orbit at a quarter to a seventh of the rate relativity and the stars' own distortion predict, and for a quarter of a century that gap was offered as evidence against general relativity. The answer was in the stars: both spin nearly on their sides, and a flattened star tipped past 54.7 degrees drags the orbit's periastron backwards instead of forwards. Explained, the slow precession becomes a measurement of something no spectrum shows — how far the spins lean towards the Earth.

stars · Binary stars

Named alongside it

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

Apsidal precessionNeutron starBinary pulsarEccentricityEquation of stateGravitational lensingMass functionPolytropePulsar timingResidualThe Shapiro delayAberration

All concepts