Galaxy Rotation Curve

A single plot that broke the universe

If you take a disk galaxy — the Milky Way, Andromeda, NGC 6503 — and measure how fast its stars and gas are orbiting at every radius from the center, you get a function v(R). This function is a galaxy's rotation curve, and it is one of the simplest, most reproducible measurements in extragalactic astronomy. It is also one of the most consequential. The shape of v(R) in nearly every spiral galaxy disagrees with the shape predicted by Newton's law of gravity applied to the visible mass. The simplest fix is to add an enormous, invisible, roughly spherical mass distribution — a dark matter halo — containing about 90% of the total galactic mass. That conclusion, drawn from the spectra of starlight, is what made dark matter unavoidable.

The disagreement is large, persistent, and visible at a glance. A Newtonian disk galaxy should have its rotation speed peak at a few kpc — somewhere near the half-light radius — and then fall off as v ∝ R⁻¹/². The actual curves rise sharply in the central kpc, then settle into a flat plateau that continues out to whatever radius the observations can reach, sometimes three to five times the optical radius. The flat plateau is universal. It is the single empirical fact that motivates dark matter at galactic scales.

What Newton predicts

Newton's shell theorem says that outside a spherically symmetric mass distribution, the gravitational acceleration is the same as if all the mass were concentrated at the center. A test particle moving on a circular orbit at radius R then has

v²(R) = G M(R) / R

where M(R) is the mass enclosed within R. The galactic disk is not spherical, but at radii well outside the bulk of the luminous matter the asymmetry is a small correction. The dominant question becomes: how does M(R) grow with R?

For visible matter — stars in the bulge and disk, plus the thin layer of cold gas — M(R) rises sharply from R = 0 to R ~ R_optical, where most of the integrated mass is contained. Beyond R_optical, where there is little visible matter to add, M(R) plateaus. Inserting M(R) ≈ M_visible_total into v²(R) = G M(R) / R gives v(R) ∝ R⁻¹/². The orbital speed must fall off as the inverse square root of radius. This is the same Keplerian decline that governs the planets: Mercury at 0.39 AU moves at 47.4 km/s, Neptune at 30 AU moves at 5.4 km/s — a factor of √77 ≈ 8.8 slower. Apply the same logic to a galaxy and you predict the outer stars should be much slower than the inner ones.

What we actually measure

The first hint of trouble came in 1939 from Horace Babcock, whose Andromeda spectra showed unexpectedly high rotation in the outer disk; with limited radial coverage and noisy data the discrepancy was easy to dismiss. The decisive work began in 1970, when Vera Rubin and Kent Ford published their long-slit Hα spectroscopy of M31, extending the rotation curve out to about 24 kpc. The curve refused to bend down. Over the next decade Rubin, Ford, and collaborators measured 60+ spirals at Kitt Peak and Lowell; nearly every one showed the same flat outer curve. The trend was inescapable.

Albert Bosma's 1978 PhD thesis added 21 cm HI radio measurements that pushed v(R) out to 3–5 times R_optical, far beyond the reach of optical spectroscopy. The HI gas is rotationally cold, thin, and traces the full gravitational potential. Begeman's 1989 paper on NGC 6503 made the case canonical: a small spiral with an optical disk of about 6 kpc but a flat HI rotation curve at v ≈ 120 km/s out to 22 kpc. The visible mass beyond 6 kpc is negligible, but the rotation curve refuses to fall.

Worked example: the Milky Way

Our own galaxy provides the cleanest single test. From measurements compiled by Eilers et al. (2019) using Gaia DR2 RGB stars combined with APOGEE spectra, the Milky Way's rotation curve has the following form, accurate to ~5 km/s out to R = 25 kpc:

v(R = 5 kpc)  = 226 km/s
v(R = 8 kpc)  = 229 km/s   ← solar circle
v(R = 15 kpc) = 218 km/s
v(R = 25 kpc) = 205 km/s
v(R = 50 kpc) ≈ 190–210 km/s  (from BeSSeL masers and stellar streams)

Now compute what Newton would predict from the visible mass alone. The Milky Way's total stellar mass is about M_star ≈ 6 × 10¹⁰ M_⊙ (most concentrated within 10 kpc) and the gas mass is M_gas ≈ 8 × 10⁹ M_⊙. Take M_visible ≈ 7 × 10¹⁰ M_⊙ and plug in at R = 50 kpc:

v_predicted = √(G M_visible / R)
            = √[(6.67×10⁻¹¹) × (7×10¹⁰ × 2×10³⁰) / (50 × 3.086×10¹⁹)]
            = √[(9.3×10³⁰) / (1.54×10²¹)]
            = √(6.0×10⁹) m/s
            ≈ 77 km/s

The Newtonian prediction at 50 kpc is about 77 km/s. The measured value is around 200 km/s. The factor of ~2.6 in velocity corresponds to a factor of ~7 in enclosed mass, all of which must be invisible. Integrated to the virial radius R_vir ≈ 200 kpc, the missing mass becomes M_vir ≈ 10¹² M_⊙ — meaning the Milky Way is roughly 90% dark matter by mass. The same calculation in NGC 6503 gives a dark-to-luminous mass ratio of ~10 by the time HI runs out, and in dwarf spheroidals like Draco it reaches ~100.

The dark matter halo

The simplest model that produces a flat rotation curve is an isothermal sphere — a self-gravitating, pressure-supported distribution of particles with constant velocity dispersion. Its density falls as ρ(R) ∝ R⁻², which gives M(R) ∝ R and therefore v(R) = const. Real halos are not exactly isothermal at all radii; the form that emerges from cold-dark-matter cosmological N-body simulations is the Navarro-Frenk-White (NFW) profile:

ρ_NFW(R) = ρ_s / [ (R/R_s) (1 + R/R_s)² ]

This is a universal two-parameter form with characteristic density ρ_s and scale radius R_s, transitioning smoothly from an inner ρ ∝ R⁻¹ cusp to an outer ρ ∝ R⁻³ decline. For a Milky Way-like halo R_s ≈ 20 kpc and the virial mass M_200 ≈ 10¹² M_⊙. Around R_s, the NFW circular velocity is approximately flat, perfectly mimicking what is observed.

The NFW form is not derived from first principles — it is a fitting function from simulations of structure formation under pure CDM gravity. But its near-universality across decades of cosmological N-body work, combined with its excellent match to observed rotation curves, makes it the standard. Variations include the Einasto profile, the Burkert profile (with a constant-density core instead of an inner cusp), and modifications driven by baryonic feedback.

Alternatives: MOND and beyond

In 1983 Mordehai Milgrom proposed Modified Newtonian Dynamics (MOND): below a critical acceleration a_0 ≈ 1.2 × 10⁻¹⁰ m/s², gravitational force is no longer GM/R² but √(GM a_0)/R — meaning v(R) = (GM a_0)^(1/4) at low acceleration, which is a constant. With one parameter (a_0) and one fit parameter per galaxy (the stellar mass-to-light ratio), MOND reproduces the rotation curves of the SPARC sample (175 galaxies, McGaugh, Lelli & Schombert 2016) remarkably well. As a near-derivation it predicts the empirical Baryonic Tully-Fisher Relation v_flat⁴ ∝ M_baryon to within observed scatter.

But MOND has serious problems at larger scales. In galaxy clusters it still needs about half the mass to be missing — even with MOND, dark matter has not been eliminated. On cosmological scales, MOND alone cannot reproduce the cosmic microwave background acoustic peak structure or the matter power spectrum at the level dark matter does effortlessly. Relativistic extensions like Bekenstein's TeVeS (2004) and Skordis & Złośnik's RMOND (2020) have been proposed, but each new test of structure formation pushes them back. The mainstream consensus remains: dark matter exists, and MOND describes some effective phenomenology that emerges from how that dark matter interacts with baryons via feedback. The story is not closed.

Famous rotation curves

Galaxy	Type	v_flat (km/s)	R_last (kpc)	M_total / M_baryon
Milky Way	SBb	~ 220	~ 50 (stellar streams)	~ 15
M31 (Andromeda)	Sb	~ 240	~ 38	~ 12
NGC 3198	SBc	150	30	~ 9
NGC 6503	SAcd	120	22	~ 6
UGC 2885	Sc	300	80	~ 20
NGC 1560	Sd (small)	80	10	~ 10
DDO 154	IB (dwarf)	50	8	~ 30
NGC 1052-DF2	UDG (anomaly)	~ 0 (no rotation)	—	~ 1

The exception at the bottom is dramatic: van Dokkum et al. (2018) reported that the ultra-diffuse galaxy NGC 1052-DF2 has stellar kinematics consistent with no dark matter at all. The result is contentious — the distance, and therefore the inferred dynamical mass, is disputed — but if confirmed it is a strong argument for dark matter as a distinct constituent that can be separated from baryons, rather than a manifestation of modified gravity that should apply everywhere.

Where rotation curves matter

• Baryonic Tully-Fisher Relation. v_flat⁴ ∝ M_baryon (McGaugh 2000, 2005) is the tightest scaling relation in extragalactic astronomy, with ~ 0.1 dex scatter across 5 decades in mass. It connects rotation kinematics to total baryonic content and remains a benchmark for galaxy formation models.
• Galactic cosmography. Combining rotation curves with stellar kinematics from Gaia gives the Milky Way potential out to 100 kpc, constraining the local dark matter density ρ_DM(R_⊙) ≈ 0.4 GeV/cm³ — a critical input to direct-detection experiments.
• Galaxy formation models. The mass-concentration relation M_vir ↔ c (where c = R_vir / R_s) calibrated from observed rotation curves anchors the comparison between ΛCDM N-body simulations and reality.
• Dwarf galaxy puzzles. The core-cusp problem (observed cores in dwarf rotation curves vs. cusps predicted by CDM), the missing-satellites problem, and the "too-big-to-fail" problem are all rotation-curve-based tensions between ΛCDM and observation that drive current research.
• Modified gravity tests. Future surveys (Euclid, LSST, SKA HI) will map rotation curves of millions of galaxies, distinguishing dark-matter halos from MOND-type universal acceleration laws with statistical leverage no individual galaxy can provide.

Common pitfalls

• Confusing rotation curve flatness with solid-body rotation. In the inner kpc, v(R) ∝ R (solid-body), reflecting the rising mass inside small R. The plateau is at large R, after the rise.
• Forgetting the inclination correction. Observed line-of-sight velocity is v_rot sin(i) where i is the disk inclination. A poorly constrained inclination scales the entire curve by csc(i) and can hide or amplify the discrepancy.
• Treating the dark halo as flat. The dark matter distribution is roughly spherical, not disk-like. The disk-plane rotation curve probes the spherically averaged enclosed mass, but the halo extends far above and below the plane.
• Quoting "Vera Rubin discovered dark matter" without nuance. Zwicky inferred missing cluster mass in 1933; Babcock saw the hint in Andromeda in 1939; Rubin's contribution was making the case overwhelming across the spiral galaxy population — the discovery is a 40-year arc with many contributors.
• Confusing the Tully-Fisher and Baryonic Tully-Fisher relations. Tully-Fisher (1977) links v_flat to total luminosity. Baryonic Tully-Fisher links v_flat to total baryonic mass (stars + gas), with a steeper slope of ~4 instead of ~3 and dramatically smaller scatter.