Astronomy & Celestial Mechanics
the sky, computed
Generator

The rotation-curve generator

The rotation curve of NGC 3198, decomposed
The rotation curve of NGC 3198, decomposed. Circular speed against radius for a three-component model of NGC 3198: a Hernquist bulge of 1.0×10⁹ M☉, an exponential disc of 2.20×10¹⁰ M☉ with a scale length of 2.6 kpc, and a pseudo-isothermal halo whose asymptotic speed is not chosen but solved for — it is whatever brings the total to the measured 150 km/s at 30 kpc, and comes out at 171 km/s. The components add in quadrature because accelerations add. The disc alone peaks at 5.7 kpc and falls away; the total does not.

Circular speed against radius for a three-component model of NGC 3198: a Hernquist bulge of 1.0×10⁹ M☉, an exponential disc of 2.20×10¹⁰ M☉ with a scale length of 2.6 kpc, and a pseudo-isothermal halo whose asymptotic speed is not chosen but solved for — it is whatever brings the total to the measured 150 km/s at 30 kpc, and comes out at 171 km/s. The components add in quadrature because accelerations add. The disc alone peaks at 5.7 kpc and falls away; the total does not.

12 essays call rotation-curve. 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.

The rotation curve of NGC 3198, decomposed. Circular speed against radius for a three-component model of NGC 3198: a Hernquist bulge of 1.0×10⁹ M☉, an exponential disc of 2.20×10¹⁰ M☉ with a scale length of 2.6 kpc, and a pseudo-isothermal halo whose asymptotic speed is not chosen but solved for — it is whatever brings the total to the measured 150 km/s at 30 kpc, and comes out at 171 km/s. The components add in quadrature because accelerations add. The disc alone peaks at 5.7 kpc and falls away; the total does not. Galaxies

A rotation curve that refuses to fall

Beyond the edge of a galaxy's light there is nothing left to enclose, so the orbital speed should fall away as the inverse square root of radius. It does not fall at all, and the shape of that refusal says the missing mass arrives at a constant rate for as far out as anyone can measure.

 the Milky Way measured from inside it — the tangent-point construction. Left: the disc in plan, with the Sun at 8.2 kpc from the centre and four lines of sight at galactic longitudes 20°, 40°, 60°, 80°. Each one grazes a circle of radius R₀ sin l, and at that tangent point the whole circular velocity lies along the line of sight, so the largest velocity in the spectrum belongs to a radius the geometry fixes and no distance has to be measured. Right: the four radii and speeds that yields, on the model curve they are read against. The construction reaches only radii inside the Sun's, which is why the outer curve — the half of it that carries the argument — needs distances after all. Galaxies

A galaxy measured from inside it

The Milky Way is the one galaxy nobody can photograph, and the only one whose rotation can be measured without knowing a single distance. A line of sight at a chosen longitude grazes a circle of known radius, and the fastest gas along it is orbiting there.

 Three mass models, 2.3× apart in the disc, agreeing to 3.1 km/s everywhere. NGC 3198's rotation curve, decomposed three ways. The stellar disc has been scaled by 0.2, 0.65, 1.1 times its photometric mass, and for each scaling the halo's asymptotic speed and core radius have been fitted — not chosen — to reproduce the same total. The heavy curve and the marked points are that total: the three models agree with it to 3.08 km/s at every radius, well inside a measurement error of 4.5. The light curves below are the disc's own contribution, and at 17 kpc they differ by a factor of 2.3 — from 36 to 85 km/s — with the halo taking up exactly the slack, 126 down to 97. The curve is one function and the decomposition asks for two. The free parameter is the mass-to-light ratio of the stars, which the kinematics never measures, and it is why a "maximum disc" fit and a halo-dominated fit are both published for the same galaxy. The degeneracy does have one hard edge: scaling the disc to 1.35 times its photometric mass cannot be fitted by any halo in the family — the best leaves 5.6 km/s rms — because past the maximum-disc solution the stars alone already overshoot the curve and a halo cannot have negative mass. That is the one thing a rotation curve says about M/L on its own, and it is an upper limit. What separates the rest has to come from somewhere else: the vertical velocity dispersion of the disc, which weighs the stars alone; gas-rich dwarfs where there is scarcely a disc to argue about; or the baryonic Tully–Fisher relation, which ties the halo's speed to the baryons and would be a coincidence if the two were independent. Galaxies

The same curve, two galaxies

A rotation curve is one function of radius. Decomposing it asks for two — how much of the speed is the stars and how much is the halo — and the mass-to-light ratio that trades one against the other is not measured by anything in the kinematics. A maximum disc and a halo-dominated fit run through identical points.

 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.

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

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.

 One velocity, two distances — 4.7 and 9.5 kiloparsecs. The radial velocity of gas along a line of sight at galactic longitude 30 degrees, against distance from the Sun, for a flat rotation curve. The velocity rises to a maximum at the tangent point — where the line of sight is tangent to a circle of radius R₀ sin l — and falls again beyond it, so every velocity below the maximum corresponds to two distances. A cloud observed at 84 kilometres a second is at either 4.7 or 9.5 kiloparsecs, and nothing about its velocity says which. The two possibilities differ by a factor in distance and by its square in luminosity and mass, so the ambiguity is not a refinement — it decides whether a star-forming region is an ordinary one nearby or a monster on the far side of the Galaxy. Galaxies

One velocity and two distances

Inside the Sun's orbit a line of sight crosses each galactocentric radius twice, and the two crossings have identical radial velocities. So a cloud's velocity gives two candidate distances, near and far, differing by a factor — and nothing about the velocity says which, though the choice decides whether the object is ordinary or extraordinary.

 The ratio that is √2 only for a point. Escape speed divided by circular speed, against distance in units of the Sun's, for a point mass and for flat-rotation-curve haloes truncated at four different radii. The point-mass value is √2 at every distance and it is a floor rather than a rule. Mass outside an orbit contributes nothing to the circular speed — a spherical shell exerts no force inside itself — and everything to the depth of the well, so any extended distribution has a ratio above √2, by an amount that measures how much lies beyond. For a flat rotation curve the escape speed is not even defined without an outer edge: the potential grows as the logarithm of radius and diverges, so the ratio is √(2 ln(R/r)) with R the truncation radius, and the answer depends on a boundary nobody can see. The measured value near the Sun is about 2.27, which on this curve puts the edge at 13 times the Sun's distance — and the strength of that inference is exactly the weakness of the method, because a ratio of 2.50 would put it at 23. Spaceflight

The ratio that is not √2

Escape speed is exactly √2 times circular speed at every distance from a point mass, and the Sun's own escape speed from the Galaxy is 2.3 times its orbital speed. The excess is not an error in either number. It is a measurement of the mass that lies outside the Sun's orbit, which pulls on nothing and holds everything.

 The observed acceleration against the acceleration the visible matter explains. The acceleration measured from a rotation curve, v²/r, against the acceleration the stars and gas alone would produce at the same radius, both on logarithmic axes in metres per second squared. The solid curve is the fitting function to the measured relation for 153 disc galaxies, g(obs) = g(bar) / (1 − e^(−√(g(bar)/a₀))) with a₀ = 1.2 × 10⁻¹⁰ m s⁻², and the band is its observed scatter of ±0.13 dex, most of which is measurement error. Above a₀ the relation follows the dashed line of equality, where the baryons account for everything; below it, it bends onto the dotted line g(obs) = √(g(bar)·a₀). At ten times a₀ the observed acceleration exceeds the baryonic by a factor of 1.06; at a₀, by 1.58; at a hundredth of a₀, by 10.5. The discrepancy the dark matter is invoked to explain is a function of one variable, and the variable is an acceleration. Galaxies

One acceleration, and every galaxy's missing mass

Plot the acceleration measured at each radius of a rotation curve against the acceleration its stars and gas would produce there, for every galaxy at once, and the points fall on a single curve. Where the discrepancy begins is set not by a radius or a mass but by an acceleration.

 A disc weighed by how far its stars rise. The share of the circular speed at 2.2 disc scale lengths supplied by the disc, against the vertical velocity dispersion of its stars measured at 1.5 scale lengths, for a disc of scale length 3 kpc whose rotation curve reaches 200 km/s, with scale heights of 0.30, 0.41, 0.60 kpc. The share grows linearly with the dispersion, because the surface density goes as its square and the disc's speed as the root of the surface density. A maximal disc, supplying about 85 per cent of the speed (dashed), needs a dispersion of 44 km/s at the middle thickness. A measured 24 ± 3 km/s (shaded) gives 0.41–0.52: a sub-maximal disc, with the halo supplying most of the speed even at the disc's own peak. Galaxies

A disc weighed by how far its stars climb

The rotation curve cannot say how much of a galaxy's speed belongs to its disc, because any split between disc and halo fits it equally well. The stars' motions perpendicular to the disc can, with no halo in the argument at all — and when that was done for thirty face-on spirals the discs came out light. But the vertical weighing needs the disc's thickness, which a face-on galaxy hides, and it needs the dispersion of the stars that hold the disc up, which the light measured may not belong to.

 Arms drawn into a disc that has none. Left: 1400 gas clouds spread uniformly over a disc between 3 and 13 kpc from the Galactic centre (cross), seen from above with the Sun (ring) 8.2 kpc out; the two thin curves mark where a two-armed spiral density wave of pitch 12° has its arms. The clouds have no arms in their positions, only in their motions: each streams by 10 km/s in a pattern locked to the arms. Right: the same clouds placed at their kinematic distances, read from their velocities with a circular rotation curve and on the correct side of every tangent point. The streaming has moved them by up to a few kiloparsecs, along lines of sight, and into bands. Lines of sight within 12° of the Galactic centre and 18° of the anticentre, where the velocities carry no distance, are left out. Counted in square kiloparsecs between 4 and 12 kpc, the densest fifth of the map holds 30 per cent of the clouds in reality and 45 per cent in the kinematic map. Galaxies

A spiral drawn by velocities alone

A map of the Milky Way's gas is built by turning each cloud's velocity into a distance with a circular rotation curve. Gas in a spiral arm does not move on circles: it streams by ten or fifteen kilometres a second in a pattern locked to the arms, so the distance errors are locked to the arms too. Read that way, a disc with no arms in it at all comes out of the calculation with arms.

 Two places one tidal dwarf could sit. The radial acceleration relation — observed acceleration against the acceleration the baryons alone produce, logarithmic, in m/s² — with the fitted relation for ordinary disc galaxies (solid) and the line of equality (dashed). The marks are a tidal dwarf galaxy of 10⁹ solar masses of gas and stars in an exponential disc of scale length 1.5 kpc, sampled every half kiloparsec to 6 kpc. If it carries no dark matter, as debris pulled from a disc should not, its points lie on the dashed line (open). If the relation is a law that every galaxy obeys, they lie on the curve (filled), and at the outermost radius the acceleration is 5.6 times the baryons' — a rotation speed 2.37 times higher. Every point is below a₀, where the two predictions are furthest apart. Galaxies

A galaxy born without a halo

When two galaxies collide, some of the gas flung out along their tidal tails condenses into new dwarf galaxies. Debris launched from a disc cannot take the dark halo with it, so these dwarfs should rotate at the speed their own gas and stars allow and no faster. The relation that predicts every ordinary galaxy's rotation from its baryons says they should rotate twice as fast. One rotation curve decides between a halo and a law — except that the inclination of a lumpy young disc can move the answer by more than the difference.

 The clusters scale with the halo and the stars do not. The total mass in a galaxy's globular clusters (lower, with its ±0.3 dex scatter shaded) and its stellar mass (upper), against the mass of its dark-matter halo, all logarithmic. The clusters follow a straight line of slope one, 3.9×10⁻⁵ of the halo, across five and a half decades. The stars follow a curve that bends at both ends round a peak near 5.8×10¹¹ M☉, because star formation is suppressed in small halos and in large ones. Stars outweigh clusters by 12 times in a halo of 10¹⁰ M☉, 879 times at 10¹², 31 times at 3.2×10¹⁴. Galaxies

A halo weighed by counting its oldest clusters

The globular clusters round a galaxy add up to about four parts in a hundred thousand of the mass of its dark-matter halo — in a dwarf, in the Milky Way, in a giant elliptical and in a cluster of galaxies, across five decades of halo mass. The galaxy's own stars follow no such rule. So a count of a few dozen old star clusters weighs a halo that its stars cannot, and the strangest thing about the rule may be that it is less strange than it looks.

The whole library · All essays