Astronomy & Celestial Mechanics
the sky, computed
Generator

The dispersion-support generator

Mass from a velocity dispersion, across nine decades
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.

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.

7 essays call dispersion-support. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Every figure listed here is the same construction drawn at different numbers, so a correction to one is a correction to all of them.

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. Gravitation

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.

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

A cluster that boils itself away

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

 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. Galaxies

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.

 Why an Einstein radius is a mass. The geometry, drawn at an angle some ten thousand times larger than the real one so that anything is visible at all. Light from a source directly behind a lens reaches the observer along every path that passes the lens at the same distance, so the image is a ring rather than a point. The ring's angular radius is θ_E = √(4GM/c² · D_ls/D_l D_s), which for a lens of 1.0×10¹² M☉ at these distances is 2.52 arcseconds and encloses 10 kpc at the lens. Rearranged, it is a mass in terms of an angle and three distances — and the mass so obtained is inside a cylinder rather than a sphere, and assumes nothing whatever about the lens being in equilibrium, which is the assumption every other weighing in this collection makes. Galaxies

Weighed by the light that bends past it

Every other mass in this collection is measured from something orbiting, which requires the system to have settled down. A gravitational lens weighs whatever is in the way with no such assumption — the light does not care whether the mass is in equilibrium.

 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. Galaxies

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.

 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. Galaxies

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.

 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. Galaxies

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.

The whole library · All essays