Astronomy & Celestial Mechanics
the sky, computed
The collection

Every essay — page 33

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 641–660 of 909.

The observed sky

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

A lucky frame stays lucky further from the star when the high air was quiet. The mean Strehl ratio of short exposures of a target at a separation from an on-axis reference star, against the shift of the high turbulent layer between the two lines of sight in units of the aperture, s = θh/D — 0.50 of an aperture is 32 arcseconds for a 2.5-metre telescope and a layer at 8 km. The atmosphere has 50 per cent of its turbulence at the ground and the rest high, with D/r₀ = 7, and 1600 simulated frames. Frames are selected by the reference star's sharpness. When the two layers' strengths vary from frame to frame (log-normal, width 0.6), the best 2 per cent start at a mean Strehl of 0.523 and keep 0.339 at the largest separation; when the strengths are steady, the best frames start at 0.321 and fall to 0.105. In the first case the Strehl never falls to six-tenths of its on-axis value within the range drawn; in the second it does so at s = 0.208, 13 arcseconds. Selecting frames by one star selects, preferentially, the moments when the high layer was weak — and a weak high layer is weak for every star in the field.

Luck that reaches as far as the quiet air

Lucky imaging keeps the rare short exposures in which a reference star happens to be sharp, and a sharp star says the air was flat along one line of sight. A target a few arcseconds away looks through different air. It turns out to be sharper in those frames than the geometry of a single turbulent layer allows, because choosing a lucky instant chooses, preferentially, an instant when the high turbulence was weak — and weakness, unlike the pattern of the turbulence, is shared across the field.

5 figures · Seeing One setting Sun drawn through a layer of warm air. The apparent shape of the Sun at five moments of its setting, seen from 20 m above the sea through an atmosphere with an inversion — 0.15 K per metre over 40 m starting 30 m up — computed by tracing rays from the observer out through the layer and mapping each apparent elevation back to the point of the solar disc it came from. The true elevation of the Sun's centre is given under each panel; the dotted outline is the Sun as a standard atmosphere would show it at the same moment, where it is above the sea horizon at all. The layer does not merely flatten the disc: it cuts it into two separate strips at some moments, stacks an inverted image above the upright one, and draws a thin neck between them, because each strip of apparent elevation maps to a strip of the true disc through a curve that folds. The shapes have been recorded for centuries — the Etruscan vase, the omega, the flattened stack — and each is a picture of the transfer curve for the air that night.

A sunset's shape is a sounding of the air

A layer of warm air over a cold sea bends the rays of the setting Sun into shapes — a flattened stack, an inverted copy balanced on the upright one, a vase with a thin neck — and every shape is a drawing of one curve: the map from where light appears to come from to where it really came from. Run backwards, the curve is a temperature profile. It can be read, but only for the lowest air, only from a known height, and only up to a trade between how strong the warm layer was and how thick.

5 figures · Refraction The zenith sky after sunset, with light scattered twice. The brightness of the zenith sky at 550 nm in magnitudes per square arcsecond, brighter upward, against the Sun's depression: light scattered once (dashed), light scattered twice (dash-dot), and their sum (solid), in a spherical atmosphere with US Standard densities, Rayleigh scattering and 300 Dobson units of ozone. The dashed horizontal line is the natural night sky, 21.9 magnitudes per square arcsecond. Once-scattered light reaches it at 9.2°; the sum, dominated beyond about seven degrees by light scattered twice, reaches it at 13.2°, and it falls to a hundredth of the night sky — five magnitudes fainter, about the smallest excess a careful measurement of the sky can detect — at 17.5°, where once-scattered light alone would have reached that level at 11.9°. The second order falls more slowly than the first because it can reach the zenith by two grazing hops — sunlight to the lit air low towards the Sun, and that air to the shadowed air overhead — each of which needs only part of the depression a single hop would. Orders beyond the second are left out, so the solid curve is a lower bound on the sky.

The night is late because the light bounces twice

Light scattered once can reach the zenith only while the Earth's shadow over the observer is still low in the air, and by that reckoning the sky should reach the darkness of night about nine degrees after sunset. It does not until eighteen. The difference is light scattered twice — first by the sunlit air hundreds of kilometres towards the sunset, then again by the shadowed air overhead — and adding it to the calculation moves the end of twilight almost exactly to where observers have always put it.

5 figures · Twilight

Spaceflight

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

Each kilometre of aim error, turned into a velocity error. The error in the outgoing velocity per kilometre of error in where a spacecraft passes a moon, in metres per second per kilometre, against the flyby altitude on a logarithmic scale, for the four large moons of Jupiter at an encounter speed of 6 km/s. The flyby turns the velocity through an angle δ with tan(δ/2) = aₕ/b, where b is the aim point's distance from the moon's centre and aₕ = μ/v∞²; missing the aim by a kilometre along b changes δ, and missing it by a kilometre across b rotates the plane of the turn, and both move the outgoing velocity by about 2v∞aₕ/b² per kilometre. At 300 km above Ganymede that is 0.32 m/s per km along and 0.32 across. The closer and more massive the flyby, the more it turns — and the more it amplifies whatever error the approach carried.

A flyby that multiplies its own aiming error

A close flyby of a moon turns a spacecraft's path by ten degrees for nothing, and it turns every error in the aim with it. A kilometre of miss at three hundred kilometres above Ganymede becomes a third of a metre per second of velocity error, and a week later two hundred kilometres of miss at the next moon. A tour of flybys is an amplifier that must be re-zeroed after every pass, and whether a small spacecraft can use one in place of an engine burn depends on how many kilometres it can aim to.

5 figures · Gravity assist What it costs to follow a turning grapple point. The velocity a chaser must spend each minute to hold station on a grapple point carried round by a flat spin, ω²r × 60 s, against the spin rate in degrees per second, for points 1, 3, 10 m from the spin axis, both axes logarithmic. The dashed line is what holding station 10 m from a target in a 400 km orbit costs in the same units, from the orbital dynamics alone: 2.3 mm/s per minute. At 1° per second, a point 3 m out costs 5 cm/s per minute — 24 times the orbital term; at 5° per second, 1.4 m/s per minute, 595 times. The cost rises as the square of the rate, so a target spinning twice as fast is four times as expensive to follow, and above a few tenths of a degree per second the orbit — the whole subject of the approach up to this point — has become the small correction.

A grapple point that moves faster than the orbit

Every rendezvous problem up to the last few metres is orbital mechanics: two bodies in nearly the same orbit, drifting apart and together under gravity. A derelict satellite tumbling at a few degrees a second changes the problem entirely. Holding station on a point of its surface costs more each minute than a week of orbital station-keeping, the moment to grab it comes and goes in seconds, and a free tumble never repeats — so the orbit becomes a small correction to a problem in rotation.

5 figures · Rendezvous Heliocentric latitude bought three ways, for the same launch. The highest inclination to the ecliptic a spacecraft can reach, against the hyperbolic excess speed it leaves the Earth with. Launched directly into an inclined orbit at 1 AU (dashed), it pays 2v⊕ sin(i/2) — 15.4° for 8 km/s and 23.2° for 12. Launched on a transfer to Mars whose plane is tilted — the wall across the porkchop plot, where transfers near half a turn are forced steeply out of the plane — it reaches latitude i half-way, for 13.7° at 8 km/s and 21.5° at 12: the same list price, starting from the cost of getting to Mars at all. Launched to Jupiter, it can turn there: a flyby rotates its Jupiter-relative velocity, and the steepest heliocentric orbit a flyby at 6 Jupiter radii can leave it on climbs to 90° for a 12 km/s launch, reaching the 79° of Ulysses's orbit (dotted) from about 10.6 km/s — close to the launch speed that mission actually used. Inclination bought near the Earth costs the Earth's orbital speed; bought at Jupiter it costs Jupiter's, less than half as much, and the flyby pays.

The wall sells inclination at the dearest price

Transfers to Mars that go almost exactly half-way round the Sun are forced into planes tilted steeply out of the ecliptic, which is why they form a wall of expensive dates across every launch chart. A mission that wanted to leave the ecliptic could use the wall on purpose. It would be buying inclination at the Earth's orbital speed — the dearest place to buy it — and the routes that reach the Sun's poles have always gone somewhere slower to make the turn.

5 figures · Launch windows

Starlight

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

A synthetic spectrum fitted with a line list that is missing the weak lines. Four ångströms of a simulated metal-line spectrum near 5000 Å at resolving power 20,000 and a signal-to-noise ratio of 200 (thin, noisy line), and the best-fitting synthesis (heavy line) computed from a line list that omits every line with a central optical depth below 0.15 — 137 of the 160 lines in the ten-ångström window, carrying 32 per cent of its total absorption. The fit adjusts one abundance for all the listed lines, with the continuum held at its true level; it returns an abundance 0.188 dex too high: the missing lines' absorption has nowhere to go but into the parameters the synthesis has. The residuals (below, offset) show the missing lines as dips that no parameter of the synthesis can fill — which is how an incomplete list announces itself at high resolution, and why it does not at low.

The lines a synthesis keeps pay for the ones it lacks

Modern abundance work fits a computed spectrum to a whole stretch of the observed one rather than measuring lines one at a time, and it averages away the errors in each line's laboratory data as it does so. It cannot average away a line that is not in its list. The absorption of the missing lines has to be put somewhere, and the fit puts it into the parameters it has — the abundance of the lines it knows, or the level of the continuum — in a direction set by which of the two it was allowed to move.

5 figures · Line formation How much light a fixed mass of dust scatters, as its largest grains grow. The scattering cross-section per unit mass of dust, on a logarithmic scale relative to the largest value drawn, against the size of the largest grains in a power-law distribution n(a) ∝ a^−3.5 from 0.005 µm, at wavelengths of 1.25, 2.2, 3.6 µm. With the largest grains at a quarter of a micron — the size the interstellar extinction curve implies — the dust scatters 3.6 µm light 33 times less efficiently than it does once they reach a micron: grains much smaller than the wavelength scatter as the fourth power of their size over the wavelength and are all but invisible. At 3.6 µm the scattering per mass peaks when the largest grains are 3.5 µm, and falls beyond, as the mass moves into grains too large to scatter efficiently per gram. A dense core that glows at 3.6 µm in scattered starlight is therefore a core whose grains have grown to a micron or more.

A glow that weighs the grains reddening cannot see

The interstellar extinction curve reads the sizes of dust grains from how much more they dim blue light than red, and it goes blind at the top end: grains larger than about a micron dim every visible colour alike. Those are the grains that hold most of the mass, and the ones from which planets begin. Dense cloud cores glowing faintly at 3.6 microns in scattered starlight are the first direct sign of them — because only grains grown to a micron or more scatter that light at all.

5 figures · Extinction Amplitudes measured near the noise, against the amplitudes that were there. A simulated survey of 6000 stars, each with one planet whose true radial-velocity semi-amplitude K is drawn from a population rising towards small amplitudes (dN/d ln K ∝ K^−0.6, from 0.2 to 20 m/s), each measured with an amplitude uncertainty of 0.6 m/s; 1500 are drawn. Across: the true amplitude; up: the amplitude fitted at the planet's period; both logarithmic. The dashed diagonal is a perfect measurement; the horizontal line is the detection threshold, 4 times the uncertainty. Points above it are the detections (filled). Well above the threshold they sit on the diagonal; near it they sit above it — the detections within one unit of the threshold have measured amplitudes 1.33 times their true ones on average — and below it the measured amplitudes stop falling at all, because an amplitude is a length and noise alone gives it about 0.75 m/s.

A prior the survey measures through its own noise

Whether a weak periodic signal is a planet depends on how large planets are expected to be, and the best source for that expectation is the population the surveys themselves have measured. But the population's small end is built from the marginal detections whose reality depends on the prior. Breaking the circle means measuring the population and the planets together, with the noise and the threshold written into the model — and doing it that way changes the answers by factors of two where they matter most.

5 figures · Periodograms

Gravitation

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

A flattened shell that pulls nothing inside it. The strength of the gravitational field of a thin homoeoid — the shell between two similar concentric ellipsoids, here with axis ratio 0.5 — along its long axis (solid) and its short axis (dashed), in units of GM/a² against distance from the centre in units of the long semi-axis a, from 20,000 equal masses placed on it. The dotted curve is a point mass. Inside the shell the field is zero along both axes to within 0.0001 — the numerical residue of a finite sampling — however flattened the shell is: this is Newton's theorem for a homoeoid, the flattened generalisation of the spherical shell's empty interior. Outside, the shell does not pull like a point: at 1.2 a along the long axis the field is 1.003 against the point's 0.694, and at 0.60 a along the short axis, just outside the shell there, 0.901. The exterior field remembers the shape; the interior does not know there is one.

The shell can be flattened if it thins at the poles

A spherical shell exerts no force anywhere inside it, and a ring exerts a force that points outward. Between them lies a shape Newton found that keeps the sphere's empty interior however far it is flattened: the homoeoid, a shell bounded by two similar ellipsoids, thick at the equator and thin at the poles in exactly the right proportion. It is why the gravity inside any flattened body can be computed at all — and why the planes of stellar streams measure the shapes of dark haloes.

5 figures · Shell theorem

Galaxies

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

The orbit that turns across the bar. The radial displacement with which a nearly circular orbit answers a weak rotating bar, against radius in kiloparsecs, in linear theory, scaled to its largest value: positive means the orbit's long axis lies along the bar, negative across it. The disc's rotation curve rises within 0.35 kpc to 220 km/s, the bar ends at 3 kpc and turns with corotation at 3.6 kpc, and the disc then has two inner Lindblad resonances, at 0.37 and 0.78 kpc (vertical lines). Without dissipation (dashed) the response diverges at each resonance and changes sign across it: orbits between corotation and the outer resonance are elongated along the bar — the x1 family that builds the bar itself — and orbits between the two resonances are elongated across it — the x2 family. Gas, which dissipates (solid, damping 8 km/s/kpc), cannot follow the flip abruptly; its orbits turn through ninety degrees over a band round the resonance, where streams of gas on neighbouring orbits cross, collide and shock. That band is where a nuclear ring forms. Near corotation the linear theory itself fails, and the curves stop short of it.

The bar leaves its gas where the orbits turn

A rotating bar in a disc galaxy drives gas inwards from kiloparsecs out, and the gas does not fall to the centre. It stops in a ring a few hundred parsecs across and forms stars there at rates the rest of the galaxy cannot match. The ring sits at the bar's inner Lindblad resonance, where the orbits that support the bar give way to orbits turned across it — so its existence says there is mass concentrated at the centre, and its size says how fast the bar turns.

5 figures · Spiral structure

Exoplanets

Planets nobody has seen, weighed and measured from a dip, a wobble and a delay.

Four transiting planets, one star, and a density that must agree. The stellar density implied by each of four transiting planets of one simulated star, each fitted as though its orbit were circular, relative to the star's true density, against orbital period, both logarithmic, with one-sigma errors of 8, 10, 14, 20 per cent from the transit fits. Three planets have eccentricities of a few hundredths and agree with one another and with the truth (dashed line: their weighted mean, 0.950). The planet at 15.9 days has e = 0.2 with its periastron 70° from the line of sight, which makes it cross the star faster than a circular orbit would, and its implied density is 2.10 times the star's — 5.2 standard deviations from the others. A single planet's density can be reconciled with any star by choosing its eccentricity; several planets must share one star, and the odd one out is caught.

Every planet of a star must agree on its density

A single transiting planet gives the density of its star only if its orbit is known to be circular, and an eccentric orbit can disguise almost any star as almost any other. A star with four transiting planets has four such measurements of one number. They must agree — and the way they disagree measures how eccentric the orbits are, catches the candidate that orbits some other star, and pins down the star better than its own spectrum does.

5 figures · Transits

Stars

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

The least a star can average, whatever it is made of. The lower bound on the mass-averaged temperature of a main-sequence star of ideal gas, μmₕGM/6kR, against its mass, both logarithmic, with the radius taken from the main-sequence mass–radius relation, for mean molecular weights of 0.5 and 0.62 — 0.5 is fully ionised hydrogen, the lightest possible gas and so the lowest floor; 0.62 is the solar mixture. For a star of the Sun's mass and radius the floor is 1.9 million kelvin whatever the star is made of. The dashed curve is the average a centrally condensed star actually has, for the index-three polytrope that approximates the Sun: three times the floor, 7.2 million kelvin for the solar mixture. The floor rises slowly with mass, because radius grows almost as fast. Nothing in it refers to nuclear reactions: mechanics alone requires a star to be hot enough inside that hydrogen is ionised throughout and fusion is within reach.

Mechanics sets a star's temperature before fusion does

A ball of gas with the Sun's mass and radius must average at least two million kelvin inside, whatever it is made of and whatever powers it, because nothing else can hold it up. The bound comes from the virial theorem and one inequality about gravitational binding. Where it is met, a body can be a star; where it demands temperatures a body plainly does not have — the Earth, Jupiter, a white dwarf — it proves that something other than hot gas is holding that body up.

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

7 figures · Degeneracy

Cosmology

One object, seen once, from inside — and every number in it the output of a model.

The polarisation peaks sit in the temperature troughs. A tight-coupling sketch of the two acoustic sources against multipole, each scaled to its own maximum: the temperature, from the compression of the photon–baryon fluid at last scattering plus its Doppler term, and the E-mode polarisation, which is driven by the fluid's velocity gradient. The sound horizon and the distance to last scattering put the acoustic scale at ℓ = 302, the measured phase of the peaks is 0.267 and the damping scale is ℓ = 1400; the baryons' weight, which changes the heights of the temperature peaks and not their places, is left out. The temperature peaks fall at ℓ ≈ 218, 515, 813, 1110, 1408; the polarisation peaks at ℓ ≈ 130, 391, 681, 975, 1272, between them, because a compression is at its extreme when the fluid is momentarily still, and the velocity is largest a quarter of a cycle later. The model has no projection, no decaying potentials and no thickness, so it places the peaks and says nothing about their shapes.

Polarised only where the glow was lopsided

An electron bathed in light from every side scatters it unpolarised, and so does one bathed in light brighter on one side. Only a glow brighter along one axis than across it leaves a polarisation behind, and in the microwave background that glow existed for a short time, in a thin part of a thin shell, because the plasma was moving.

5 figures · Microwave background The spectrum ahead is the spectrum behind, shifted. The difference between the local power spectrum of the microwave background measured in the direction of the observer's motion and in the opposite direction, as a percentage of the average, against multipole, at 369.8 km/s. Ahead, the pattern is compressed and every feature moves to a higher multipole by a factor 1 + β; behind, to a lower one by 1 − β, with β = 1.234e-3. The difference is largest where the spectrum is steepest — 2.30 per cent at most in this range — and passes through zero at every peak and trough, where a small shift changes nothing. The spectrum is an interpolation through the measured peaks and troughs, so its steepness between them is known only to that resolution. A per cent of a quantity estimated from millions of independent modes is a measurable number.

A speed read from the finest spots in the map

The dipole in the microwave sky is taken to be the observer's motion and subtracted. The same motion also squeezes the small-scale pattern towards the direction of travel by a part in a thousand, and that squeeze measures the speed a second time without using the dipole at all — which is the only way to ask how much of the dipole is the universe's own.

5 figures · Microwave background

Exoplanets

Planets nobody has seen, weighed and measured from a dip, a wobble and a delay.

Two models through the same points, one with a planet. 60 simulated velocities over 200 days of a star rotating every 25 days, with spot activity drawn from a quasi-periodic process of 3 m/s amplitude and a 40-day decay, white noise of 1 m/s, and a planet of 2.5 m/s semi-amplitude on a 27.0-day orbit — 1.08 times the rotation period. Two models are fitted, each with its activity amplitude and decay time chosen by maximum likelihood: activity plus the planet (solid), whose planet comes out at 3.3 m/s, and activity alone (dashed), which has no planet in it. Both pass through the points; the likelihoods differ by 3.6 in the logarithm, a ratio of 34.9, which is not decisive. Near the rotation period a flexible model of the star can carry the planet's signal by itself.

A model of the star that can carry the planet

Starspots are modelled as correlated noise with the star's rotation period built in, and the model is good enough to remove signals that would otherwise be announced as planets. It is also good enough to remove planets. Near the rotation period no amount of data in one season separates the two, and away from it a flexible noise model quietly takes most of the evidence a planet would have had.

5 figures · Reflex velocity The red velocity against the blue one. Simulated velocities of two stars, each measured simultaneously in a blue channel at 0.45 micrometres (across) and a red one at 0.8 (up), with 0.6 m/s of noise in each. One star (filled) has a 4 m/s planet on a 9.3-day orbit: its points lie along the diagonal, slope 0.97, because a Doppler shift is the same velocity in every colour. The other (open) has no planet and a spot 900 K cooler than its 5100 K photosphere, rotating every 12 days and evolving: its points lie along a shallower line, slope 0.64, which is the ratio of the spot's contrasts in the two channels. Neither star's blue velocities alone say which it is. The flux-deficit spot model is the simplest case; a spot whose signal is mostly suppressed blueshift lies closer to the diagonal.

A planet has one velocity in every colour

A Doppler shift moves every line in a spectrum by the same velocity, because it is the whole star that moves. A starspot does not move anything; it removes light, and how much light it removes depends on the colour it is seen in. Measuring the velocity in the blue and the red at once separates the two in a single night — until the star's magnetic field, which grows more visible towards the red, takes part of the separation back.

5 figures · Reflex velocity Curvature turns a curve into an island. The companions — separation across, mass up, on edge-on circular orbits — that fit simulated velocities of a star with a 12 Jupiter-mass companion at 22 AU (period 103 years), within a Δχ² of 6.18 of the best fit, for observing spans of 6, 12, 24 years with 60 velocities of 2 m/s noise each. Over 6 years the data show a slope and the allowed region is a long band along M ∝ a², spanning 3.3 decades in mass. As the span grows the curvature of the orbit appears, and the band contracts — to 1.0 decades at 24 years — towards the true companion. Inclined orbits, which the fit leaves out, would extend every region to higher masses.

A slope that is a whole curve of companions

A companion on an orbit longer than the observations shows in a star's velocity only as a steady drift. A drift is one number, an acceleration, and an acceleration is a mass divided by the square of a distance — so a planet at ten astronomical units and a brown dwarf at fifty produce the same straight line. Turning the line into a companion needs the orbit to bend, a picture, or the other two components of the acceleration.

5 figures · Reflex velocity

Galaxies

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

How often a merged black hole is kicked faster than a given speed. The fraction of black-hole mergers whose remnant recoils faster than the speed on the horizontal axis, for mass ratios spread evenly in logarithm between 0.1 and 1 and spins of 0.9 of the maximum, spins within 10° of the orbit; spins within 30°; spins in random directions, from a fitting formula to numerical relativity over 4000 random mergers each. The median kicks are 113, 234, 534 km/s. Faster than 100 km/s: 62%, 85%, 96%. Faster than 300 km/s: 5%, 40%, 74%. Faster than 1000 km/s: 0%, 2%, 25%. Gas accreted before the merger is expected to drag both spins into line with the orbit, which is why a gas-rich merger keeps its black hole far more often than a gas-poor one.

A black hole thrown out by its own merger

When two black holes merge, the gravitational waves they radiate carry momentum as well as energy, and unless the two are identical and spinning in line the waves leave more in one direction than the other. The remnant recoils, at anything from a few tens of kilometres a second to several thousand — and whether a galaxy keeps its central black hole comes down to how its spins were pointing in the last orbit before the merger.

5 figures · Galactic nuclei

Every sequence · Every named object · Search