Astronomy & Celestial Mechanics
the sky, computed
Concept

Velocity dispersion — where it appears

The spread of speeds within a bound system of stars or galaxies, which is what supports it against its own gravity. It is measured from the width of a galaxy's absorption lines, and squaring it and multiplying by a radius gives a mass with no orbit solved anywhere.

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

An average that arrives, and a snapshot that never does. 2T/|U| against time for Burrau's problem — three masses released from rest, and never repeating, in two forms: the instantaneous ratio, and the ratio of the running time-averages. The instantaneous one runs between 0.03 and 1.97 — the system is a long way from equilibrium at almost every moment, because the bodies are alternately falling together and flying apart. The averaged one settles: after 50 time units it is 1.0142, against the exact 1 the virial theorem requires of any bound system. This system is not periodic and is not even permanently bound — Burrau's problem ejects its lightest body — so the average is drifting rather than converged, and that is the theorem's own condition made visible: it is a statement about bound systems and says nothing about anything that is leaving. What the figure cannot show is the error: a cluster observed once gives 2T/|U| with a scatter of this size, and what makes its mass believable is not the measurement but the relaxation.

An average that weighs what cannot be watched

A cluster's mass can be had from its speeds alone — no orbit followed, no period observed, no distance to any single star. The theorem that allows it is an average over time, which is exactly what a photograph is not.

gravitation · Virial theorem The wake, computed from the streamlines that make it. Left: 26 streamlines past a point mass, in the mass's own frame, each integrated from far upstream with the same speed and a different impact parameter, and mirrored about the axis. Nothing is drawn to converge — every track is the hyperbola its own impact parameter gives it, and they cross downstream because an attraction focuses. Right: the density that focusing produces at 3.2 focusing radii behind the mass, as (b/y)(db/dy) — the Jacobian of the map from starting radius to arrival radius — which peaks at 15.53 times the background at 0.03 radii off the axis. The overdensity is behind the mass, and that is the entire mechanism: the wake pulls backwards on the body that made it. What the figure cannot show is the steady state, because it has no time in it: a real wake is continuously replenished, and the drag is the sum over an infinite train of these encounters, which is where Chandrasekhar's logarithm comes from.

A drag with nothing to drag against

A massive body moving through a sea of light ones raises an overdensity behind itself and is pulled back by it. The force does not depend on the masses of the background bodies at all — and it is strongest at one particular speed.

gravitation · Dynamical friction Mass from a velocity dispersion, across nine decades. The virial estimate M = 5 σ²R_e/G, drawn for half-light radii from 0.2 to 20 kpc, with four systems placed on it: M32 at σ = 75 km/s and R_e = 0.1 kpc, giving 6.5e+8 M☉; the Milky Way's bulge at σ = 105 km/s and R_e = 1 kpc, giving 1.3e+10 M☉; M87 at σ = 320 km/s and R_e = 7 kpc, giving 8.3e+11 M☉; the Coma cluster at σ = 1000 km/s and R_e = 1500 kpc, giving 1.7e+15 M☉. The same expression spans from a dwarf elliptical to a cluster of galaxies — nine decades of mass — because the virial theorem does not care what the moving objects are. In the last case they are whole galaxies, and the estimate exceeds everything visible in the cluster by a factor of about fifty, which is where the problem was first noticed.

A galaxy held up by disorder

An elliptical galaxy barely rotates. What holds it up against its own gravity is the randomness of its stars' motions, and the virial theorem turns that randomness into a mass by exactly the reasoning that turns a rotation speed into one.

galaxies · Velocity dispersion A cluster weighed three ways, and its stars weighed once. A cluster of galaxies with a velocity dispersion of 1000 km/s, gas at 8 keV and a strong-lensing Einstein radius of 25 arcseconds, each turned into a mass inside 1.5 Mpc by its own relation and nothing else: 7.0, 8.9 and 15.2 × 10¹⁴ M☉. The three assume, respectively, that the galaxies are in equilibrium, that the gas is, and nothing whatever — so their agreement to within a factor of 2.2 is not three restatements of one assumption. The lensing bar is the loosest of the three and is drawn that way deliberately: it measures the mass inside a cylinder of radius 109 kpc, 1.11×10¹⁴ M☉, and carrying that out to 1.5 Mpc as though it were a sphere overstates it. The stars are 2.9 per cent of it.

A cluster weighed three ways

The speeds of a cluster's galaxies, the temperature of its gas and the bending of light behind it are three measurements with almost nothing in common. They agree within a factor of two, and all three exceed the mass of the stars by about a hundred.

galaxies · Clusters S2's orbit, and the mass it implies. The orbit of S2 about the centre of the Milky Way, drawn from its measured elements: a semi-major axis of 0.1251 arcseconds, which at 8.28 kpc is 1036 AU, an eccentricity of 0.8843, and a period of 16.05 years watched round twice. Kepler's third law in solar units — a³/P² — gives 4.31 million solar masses, by exactly the calculation that weighs a planetary system. At periapsis the star is 120 AU from the focus, 1407 Schwarzschild radii, and the mean density inside that radius is 5.3e+15 M☉ per cubic parsec — a million times the densest star cluster known. That density, not the mass, is the argument: there is no configuration of stars that would fit.

The mass at the centre that is not stars

One star at the centre of the Milky Way has been watched round a complete orbit. Its semi-major axis and its period give four million solar masses by Kepler's third law — and the volume that mass occupies is small enough to rule out every alternative to a black hole.

galaxies · Galactic nuclei A delay of 81 days, and a sheet nobody can see that moves H₀ to 82.4. Above: the arrival-time surface of a lensed source, along the line through the lens. The curve is the Fermat potential in days — the geometric cost of taking a longer path, minus the gravitational cost of climbing out of the potential — and the images sit at its stationary points, at -1.20″ and 2.04″, which for an isothermal sphere is β ± θ_E. Fermat's principle is doing all of the work here: light does not take the shortest path or the quickest one, it takes every stationary one, and the number of images is the number of stationary points. The vertical distance between the two is 81 days, and it is measurable — the source is a quasar, quasars vary, and the same wiggle appears in one image and then the other. That single number carries an absolute distance: the delay is D_Δt/c times a dimensionless function of the lens model, and D_Δt goes as 1/H₀, so a monitoring campaign gives the Hubble constant with no rung of any ladder beneath it. Below: the two light curves, shifted by exactly that delay. The dashed curve is the second image with the delay removed, and the agreement is the measurement. What the picture also shows is the reason the answer keeps moving. The second arrival-time curve is the same lens with a uniform sheet of convergence added and the source moved to compensate: every image sits at the same place, every flux ratio is the same, every image shape is the same, and the delay is λ = 0.85 times as long. A lens model fitted to positions alone cannot see the sheet, and inferring H₀ from the same delay under it gives 82.4 instead of 70 — a 15 per cent shift with no observable attached. Breaking it needs a mass measured some other way: the velocity dispersion of the deflector, or a count of everything else along the line of sight.

A distance measured with a stopwatch

The images of a lensed quasar sit at the stationary points of an arrival-time surface, and the height between two of them is a delay in days. That delay is proportional to a distance, and the distance is proportional to one over the Hubble constant — so a flickering quasar gives H₀ with no ladder under it.

galaxies · Strong lensing Fingers 3.1 times long, a large scale squashed to 0.92, and a test of gravity. Left: a 360-megaparsec slice of a clustered universe, as the galaxies actually sit. Right: the same galaxies as a redshift survey records them, with the line of sight up the page. Nothing has moved sideways, because an angle is an angle; every displacement is along the line of sight, because that coordinate came from a redshift and a redshift is the expansion plus whatever the galaxy is doing on its own. Two effects, opposite in sign and separated by scale. Inside a cluster the motions are virial and random, 720 kilometres a second of them, which at H₀ = 67.36 is 11 megaparsecs of smearing on an object a few across: the clusters become fingers 3.1 times longer than they are wide, all pointing at the observer, which is the one structure in cosmology that is definitely not real. On the scale of a supercluster the motions are coherent — everything is falling in, so the far side is approaching and the near side receding — and the structure is compressed rather than stretched, to 0.92 of its true extent here, measured on cluster centroids so the fingers have already averaged away. Below: why that compression is worth having. Its amplitude is the rate at which structure is currently growing, and the growth rate is Ωₘ(z) raised to a power that general relativity fixes at about 0.55. A theory of gravity that differs from general relativity on cosmological scales while matching every solar-system test changes that exponent and nothing else, and the two curves drawn — γ = 0.55 and γ = 0.68 — differ by only 4 per cent at redshift a half, against error bars of 12 per cent on the points beside them. The worst systematic in a redshift survey is the measurement — and it is a hard one, because a quarter of a change in the exponent that governs how gravity assembles structure moves the observable by less than the width of the curve it is drawn on. What the picture cannot show is the degeneracy that limits it: what is measured is fσ₈, a product, and separating the growth rate from the amplitude of clustering needs something else entirely.

A map stretched by the thing it measures

A redshift survey plots galaxies at distances derived from their redshifts, and a galaxy's redshift contains its own motion as well as the expansion. So the map is systematically distorted — squashed on large scales, drawn out into radial spikes on small ones — and both distortions are caused by the gravity the survey exists to measure.

cosmology · Large-scale structure Dynamical friction is strongest at 0.97 times the dispersion, and weaker either side. The dynamical friction on a body moving through a Maxwellian sea of stars, against its own speed in units of the square root of two times the velocity dispersion. The dashed curve is the fraction of the field moving slower than the body, which is the only part of the field that contributes: a star overtaking from behind pulls the body forward exactly as often as one being overtaken pulls it back, so the fast tail cancels out entirely. That fraction rises from nothing to one. The drag is that fraction divided by the square of the speed, and the quotient of a saturating numerator and a growing denominator has a single maximum, here at 0.97. The consequences run in both directions. A body moving much faster than the stars around it is barely slowed at all, which is why a galaxy passing through a cluster at a thousand kilometres a second does not sink; and a body already at rest with respect to the field feels nothing either, because there is no wake to be behind it. Sinking is therefore fastest in the middle of the process rather than at the start of it.

A drag that is strongest in the middle

Dynamical friction is not like air resistance. Only the stars moving slower than the body contribute at all, so the drag rises from nothing, peaks when the body has slowed to about the speed of everything around it, and falls away again — which is why a fast satellite is never captured and a slow one sinks in a hurry.

gravitation · Dynamical friction Gas kept inside 3.6 to 6.8 kiloparsecs, depending only on the speed of the fall. The pressure holding a disc galaxy's gas down, against radius, with the ram pressure of the cluster gas it is falling through drawn as three levels. The restoring pressure is two pi times the gravitational constant times the product of the stellar and gaseous surface densities, both exponential with scale lengths of 3 and 5 kiloparsecs, so it falls steeply outward; the ram pressure is the intracluster density times the square of the infall speed and does not depend on radius at all. Where the level crosses the curve, the gas goes. A galaxy entering at 700 kilometres a second keeps everything inside 6.8 kiloparsecs; at 1600 it keeps only 3.6. Nothing in this touches the stars, which are not a fluid and feel no pressure at all, so the galaxy emerges with its stellar disc and its rotation curve intact and its star formation stopped from the outside in. That is the mechanism behind the most conspicuous fact about clusters: their spirals are red, gas-poor and still spiral-shaped, which no process acting on the stars could produce.

Red, gas-poor, and still spiral-shaped

A galaxy falling into a cluster meets a headwind of hot gas at two thousand kilometres a second. Where that ram pressure exceeds the disc's own grip on its gas, the gas goes; the stars, which are not a fluid and feel no pressure, do not. The result is a spiral with no fuel, and it is the commonest kind of galaxy in a cluster.

galaxies · Clusters Prograde and retrograde, the same encounter. The same encounter twice: a companion of 1 times the primary's mass passing at 1.5 disc radii, with the disc spinning the same way the companion orbits and then the other way. Nothing else differs. The prograde disc grows a bridge towards the companion and a tail away from it; the retrograde one is barely disturbed. The reason is a resonance rather than a force: in the prograde case the outer particles keep pace with the companion for a substantial part of the encounter and are pulled the whole time, and in the retrograde case they sweep past it and the impulses cancel.

Whether it merges is a ratio of two times

Two galaxies passing each other either merge or do not, and what decides it is not how close they come. It is whether the encounter lasts long enough for the internal motions to respond — a slow, prograde passage transfers orbital energy into stellar orbits and the pair is bound; a fast one leaves both galaxies heated and still moving.

galaxies · Galaxy interactions An arm that moves angular momentum outwards, and is spent doing it. Where a 2-armed spiral pattern of 25 km/s per kiloparsec takes angular momentum from the disc and where it gives it back, against galactocentric radius, for the Milky Way's rotation curve. The resonances are found on that curve rather than placed: the inner Lindblad resonance at 1.91 kiloparsecs, corotation at 8.67, the outer Lindblad resonance at 13.80. The lower curve is the angular momentum removed per unit radius and the upper one what the wave delivers; between them runs the flux the wave carries, which is flat across the whole region where nothing is resonant, because a wave that is not interacting with anything simply travels. The two deposits are required to cancel to a part in a million — the pattern is a conveyor and keeps nothing — so the disc's total angular momentum is unchanged while its distribution is not. That is the sense in which a spiral is a machine: it moves mass inwards by moving angular momentum outwards, which is the same operation an accretion disc performs and the same one a protostellar disc must perform to make a star. The cost is paid in random motion. Every exchange heats the stellar disc, a hotter disc supports a weaker wave, and the pattern that did the work is the thing the work destroys. What the figure cannot show is the pattern's own lifetime, which is the unsettled part of the subject.

An arm that is undone by the work it does

A spiral pattern takes angular momentum from the inner disc and delivers it to the outer, which lets mass move inwards without violating anything. It is paid for in random motion, and a disc with too much random motion cannot carry a wave — so the pattern destroys the conditions it needs to exist.

galaxies · Spiral structure A virial mass inflated 23.2-fold by orbits nobody resolved. The factor by which a virial mass is overestimated when the velocity dispersion is measured from single-epoch spectra, against the system's true dispersion, for four numbers of observing epochs. Every star in a binary carries its own orbital velocity, which adds in quadrature to the system's own, and the orbital velocities of ordinary binaries are of order a kilometre a second. A system whose real dispersion is 0.3 km/s therefore measures 1.45, and since a virial mass goes as the square of the dispersion the mass comes out 23.2 times too large. Repeat epochs fix it: the orbital velocities are uncorrelated between visits while the system's own are not, so the binary variance falls as one over the number of epochs and the correction is measured rather than modelled.

A dispersion inflated by orbits nobody resolved

A velocity dispersion measured from single spectra is not the dispersion of the system's centres of mass. Every star in a binary carries its own orbital velocity of a kilometre or so, and for a dwarf galaxy whose real dispersion is smaller than that, the measured value — and the dark-matter content computed from its square — is mostly binaries.

galaxies · Velocity dispersion Three mass models, one rotation curve. Rotation curves for three decompositions of the same galaxy, from a disc contributing 30 per cent of the outer rotation to one contributing 95. Each is the quadrature sum of a stellar disc, whose shape is fixed by the light distribution and whose amplitude is an unknown mass-to-light ratio, and a dark halo with parameters of its own. All three reproduce the same flat outer rotation, because the halo's amplitude is adjusted to make up whatever the disc does not supply. They differ in the inner few kiloparsecs, by 27 kilometres a second here — which is more than the measurement error and less than the uncertainty in the disc's own contribution, since that depends on a mass-to-light ratio nobody measures directly. The maximum-disc assumption picks the largest disc consistent with the data, and it is a convention rather than a result.

Three mass models that fit the same curve

A rotation curve is one function of radius and it is fitted with three components, only one of which is known. The stellar disc's contribution scales with a mass-to-light ratio nobody measures, and whatever the disc does not supply the dark halo does — so a family of models reproduces the same curve exactly, and choosing among them is a convention rather than a measurement.

galaxies · Dark matter A single-epoch mass, good to a factor of 2.6. The error budget of a black-hole mass obtained from one spectrum, as probability densities in the logarithm of the ratio of the estimate to the truth. Four independent widths go in. The radius comes from the luminosity through a relation with 0.19 dex of scatter about it. The line width has to be measured on a profile that is blended with narrow lines and with iron, and it enters the mass squared, so 0.12 dex. The luminosity was not measured at the epoch the relation was calibrated for and the object has varied since, which is 0.053 dex. And the term this figure exists for is orientation: the virial factor is one number fitted to a sample, every individual object has its own according to how its broad-line region is inclined, and the spread of that is 0.35 dex — the largest of the four by some margin. In quadrature the total is 0.42 dex — a factor of 2.6 — and the drawn total's half-width is checked against that combination rather than assumed. The bracket at the foot is a different kind of error and is not in the quadrature: the absolute calibration of the virial factor, uncertain by about a factor of 3, which moves every quasar mass ever published by the same 0.24 dex and widens nothing.

The factor that multiplies every quasar mass

A reverberation lag and a line width give a length times a velocity squared, which is a mass multiplied by an unknown number of order one. That number is fixed by making a few dozen nearby active galaxies lie on a correlation measured in quiescent ones — so every black-hole mass at redshift two rests on a fit performed at redshift zero.

galaxies · Quasars Where the stars that are torn apart come from, round a 10⁶ solar-mass hole. The rate at which stars are delivered onto orbits reaching the tidal radius, per logarithmic interval of their orbital radius, for a hole of 10⁶ solar masses at the centre of a nucleus that is isothermal outside the hole's influence and a Bahcall–Wolf cusp within it, with a velocity dispersion of 54 km/s set by the black hole mass–dispersion relation. The radius is in units of the hole's influence radius, GM/σ² = 1.47 pc. Close in, a star's orbit is so short and its angular momentum changes so slowly that the loss cone empties every orbit, and the rate is limited by relaxation diffusing stars into it; far out, the angular momentum wanders across the whole cone in one orbit, the cone is full, and the rate is limited by how many stars there are and how long they take to come in. The flux peaks at 0.27 influence radii, near where the two regimes meet (0.19), and the total is 7.4·10⁻⁵ disruptions a year. The nucleus model is the simplest there is and real nuclei differ from it by factors of several; the shape — a narrow band of radii supplying most of the flares — is what survives.

The stars a black hole eats come from a narrow band

A star is torn apart only if its orbit happens to point almost exactly at the hole — within a millionth of the possible directions of its angular momentum. Stars on those orbits are gone within one orbit, and the only way new ones arrive is by the slow random walk of encounters with other stars. Close to the hole that walk is too slow; far from it the orbits take too long; and most of the stars a hole destroys come from a narrow band between.

galaxies · Tidal disruption Five speeds from the same galaxies, and the relation each one gives. The exponent and the scatter of the baryonic Tully–Fisher relation fitted to the same 90 model galaxies with five different speeds, fitted as speed on mass. outer speed: exponent 3.98, scatter 0.080 dex in mass; peak speed: exponent 3.86, scatter 0.110 dex in mass; at 2.2 scale lengths: exponent 3.29, scatter 0.167 dex in mass; W50 ÷ 2: exponent 3.73, scatter 0.119 dex in mass; W20 ÷ 2: exponent 4.05, scatter 0.118 dex in mass. The galaxies were built with the relation in their outer halo speed, and the other four speeds each lose some of it: the ones read from the inner curve inherit how concentrated each galaxy's stars are, and the line widths add the turbulence of the gas, which matters most in the slowest galaxies.

The speed a line width stands in for

A galaxy does not have a rotation speed. It has a rotation curve, rising in the smallest galaxies and peaking early in the largest, and the Tully–Fisher relation is fitted to whichever single number is read off it. Build galaxies with the relation placed in their outer speed and read four other speeds from the same curves, and each gives a shallower or looser relation — which is why the choice of speed is a statement about where the relation lives.

galaxies · Tully–Fisher How the scatter depends on how much dispersion is added to rotation. The scatter of the baryonic Tully–Fisher relation, in dex of mass, when the same 80 turbulent model discs are measured with S = √(K Vᵣₒₜ² + σ²), against the weight K given to the rotation, for rotation measured at 1 and 2.2 disc scale lengths. Measured at 1 scale length the scatter is smallest, 0.076 dex, at K = 0.55, and is 0.076 dex at K = 0.5; measured at 2.2 scale lengths the scatter is smallest, 0.077 dex, at K = 0.23, and is 0.134 dex at K = 0.5. For an exponential disc with constant dispersion the pressure correction is exactly K = one over twice the radius in scale lengths, so the best weight depends on where the rotation is measured. The widely used K = 0.5 is the value that makes a pure rotator and a pure isothermal sphere of the same mass agree; it is also the pressure correction for rotation measured at one scale length, and not at any other.

A disc that turns slower than its mass requires

The star-forming discs of ten billion years ago were not the thin, cold, orderly discs of today. Their gas moved randomly at tens of kilometres a second, and that motion is a pressure that holds up part of each disc, so it turns more slowly than its mass alone would require. Read those rotation speeds as if rotation did all the work and the Tully–Fisher relation tilts and scatters as if galaxies had evolved, when what differs is how much of their weight is carried by disorder.

galaxies · Tully–Fisher

Named alongside it

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

Dark matterVirial theoremThe baryonic Tully–Fisher relationCoulomb logarithmDynamical frictionImpulse approximationM sigma relationMass-to-light ratioRotation curveDynamical massGalaxy clusterGalaxy merger

All concepts