Astronomy & Celestial Mechanics
the sky, computed
Generator

The galaxy-profile generator

The mass of NGC 3198 inside each radius, weighed two ways
The mass of NGC 3198 inside each radius, weighed two ways. The mass inside radius r from the rotation curve, M = v²r/G, against the mass inside r from the light. The luminous curve flattens by about 10 kpc because there is no more disc to enclose; the dynamical curve rises almost linearly, because a flat rotation curve means exactly M ∝ r. At 30 kpc the two differ by a factor of 6.8. The straightness of the upper line is the finding: it is not that a galaxy has more mass than expected, it is that the extra mass keeps arriving at the same rate however far out one looks.

The mass inside radius r from the rotation curve, M = v²r/G, against the mass inside r from the light. The luminous curve flattens by about 10 kpc because there is no more disc to enclose; the dynamical curve rises almost linearly, because a flat rotation curve means exactly M ∝ r. At 30 kpc the two differ by a factor of 6.8. The straightness of the upper line is the finding: it is not that a galaxy has more mass than expected, it is that the extra mass keeps arriving at the same rate however far out one looks.

5 essays call galaxy-profile. 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 mass-to-light ratio of NGC 3198, radius by radius. The dynamical mass inside each radius divided by the light inside it, in solar units. Inside two disc scale lengths it is nearly flat at about 2.2 — a galaxy made of stars, weighing what stars weigh. Outside the disc the light stops and the mass does not, so the ratio climbs to 10 by 30 kpc with no sign of turning over. The curve is not fitted: it is the rotation curve's enclosed mass divided by the light profile's enclosed light, both drawn elsewhere in this collection. Galaxies

The mass that is not the light

A galaxy's mass-to-light ratio is not a number, it is a curve, and it rises without turning over. What stellar populations can plausibly weigh sets a ceiling; the dynamics sit far above it, and the gap is a shape rather than a discrepancy.

 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.

 A cusp and a core in the same halo: 39 times apart in density, 2.34 in rotation. Two dark-matter density profiles for a halo of 10¹⁰ solar masses, drawn on logarithmic axes so a power law is a straight line and its exponent is a slope. The upper curve is the profile cold dark matter simulations produce: density rising inward without limit, with a logarithmic slope tending to -1.01 — a cusp. It is not a fit to observations; it emerges from simulations run by many groups with different codes, and its robustness is what makes it a prediction worth testing. The lower curve is a cored profile with the same virial mass, flat inside a kiloparsec, slope -0.01. What the rotation curves of gas-rich dwarf galaxies prefer is the second. The disagreement is stark where it is drawn — a factor of 39 in density at 0.10 kiloparsecs — and much less stark in what is actually measured, because a rotation speed depends on the mass enclosed rather than on the local density, and integrating a cusp over a small radius does not accumulate very much. At 1 kiloparsecs the two halos differ by a factor of 2.34 in circular speed — a real difference, and a far smaller one than the 39 in density that produced it. That compression is the whole difficulty of the problem. The observable is an integral of the quantity in dispute; integrating a cusp over a small radius does not accumulate much mass, so the sharpest disagreement lives where the instrument is bluntest. And the innermost points of a rotation curve are also the ones most affected by non-circular motions, by beam smearing, and by the inclination assumed — so the measurement is hardest exactly where it matters most. Whether the resolution is astrophysical — supernova feedback moving gas repeatedly and dragging the dark matter outward with it — or a statement about what dark matter is remains open, and the figure deliberately shows only the alternatives rather than choosing. Galaxies

An argument about the innermost kiloparsec

Simulations of cold dark matter have produced the same halo for thirty years — a density rising inward without limit. The rotation curves of the smallest galaxies say the density flattens off. The disagreement is confined to a region a thousandth of the halo's size, and it has not been settled.

 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.

The whole library · All essays