Galactic Rotation Curves and Mass Modeling

A galactic rotation curve is the simplest possible experiment with a galaxy: watch how fast material orbits at each radius, and let gravity do the accounting. A star or a cloud of gas on a nearly circular orbit is in free fall, so its centripetal acceleration $v_c^2/R$ must be supplied by the gravitational pull of everything interior to it. That single balance, $v_c^2(R)=GM(<R)/R$ , turns a Doppler shift measured as a function of projected radius into a direct weight of the enclosed mass — without ever asking how bright the matter is or what its mass-to-light ratio might be. The rotation curve is the disk-galaxy counterpart of the pressure-support diagnostics used for ellipticals and clusters (the Jeans equations and the virial theorem), and it inherits its machinery from the two-body problem, galactic potential models, and orbits in spherical potentials. The result that made this technique famous is an anticlimax that became a revolution. Newtonian intuition, calibrated on the Solar System, says that once you are outside essentially all the mass the orbital speed should fall off Keplerian-fashion, $v_c\propto R^{-1/2}$ . Instead, spiral after spiral shows a curve that rises through the inner disk and then stays flat — not falling — out to and beyond the last glimmer of starlight, wherever cold hydrogen can be traced in the 21 cm line. A flat curve forces $M(<R)\propto R$ and $\rho\propto R^{-2}$ : mass keeps piling up exactly where the light has faded. Mass modeling makes this quantitative, decomposing the observed $v_c^2$ into the quadrature sum of bulge, stellar disk, gas, and a dark halo fit to close the gap.

Key takeaways

A rotation curve plots the circular speed $v_c(R)$ of stars and gas on near-circular orbits against galactocentric radius; for a spherical distribution $v_c^2(R)=GM(<R)/R$ , so kinematics alone weigh the enclosed mass.
The empirical shock is that $v_c(R)$ does not decline Keplerian-fashion past the starlight — it stays remarkably flat to the last measured point.
A flat curve requires $M(<R)\propto R$ and $\rho\propto R^{-2}$ , the signature of an extended dark-matter halo that dominates the mass budget at large radius.
Mass modeling decomposes $v_c^2=v_\mathrm{bulge}^2+v_\mathrm{disk}^2+v_\mathrm{gas}^2+v_\mathrm{halo}^2$ , because fields superpose and circular speeds add in quadrature, never linearly.
The measurement spans optical Hα and stellar Doppler shifts inside the disk and 21 cm neutral-hydrogen emission far beyond it — where the flatness, and thus dark matter, is clearest.
Key idea: for circular orbits $v_c^2(R)=GM(<R)/R$ ; the observed flatness of $v_c(R)$ at large radius forces $M(<R)\propto R$ and $\rho\propto R^{-2}$ , revealing an unseen dark-matter halo whose mass keeps growing where the light has faded.

Where this lives today

Rotation curves are still a frontier. The SPARC database (Lelli, McGaugh and Schombert, 2016) gathered homogeneous Spitzer photometry and HI/Hα curves for 175 disk galaxies, and from it emerged the radial acceleration relation — a startlingly tight link between the observed centripetal acceleration and that expected from the baryons alone, which dark-matter models must reproduce as an emergent regularity and MOND takes as fundamental. In the Milky Way, Gaia astrometry of luminous red giants (Eilers and collaborators, 2019) finds the circular speed gently declining beyond the Sun, roughly $-1.7\,\mathrm{km\,s^{-1}\,kpc^{-1}}$ , tracing the disk-to-halo transition in our own back yard. The cusp–core problem persists, now framed as a "diversity of rotation curves": dwarfs show a spread of inner slopes that neither NFW cusps nor a single cored law explains, driving work on baryonic feedback and self-interacting or warm dark matter. And ALMA and JWST now resolve rotation curves of massive disks only a few billion years after the Big Bang. The balance $v_c^2=GM(<R)/R$ remains the lever under it all.

The physical picture

Everything begins with a single idealized orbit. A star at radius $R$ on a circle feels the gravity of all the mass interior to it — and, by Newton's shell theorem, of nothing exterior, provided the distribution is spherical. Balancing the inward field $g=GM(<R)/R^2$ against the centripetal requirement $v_c^2/R$ gives $v_c^2=GM(<R)/R$ directly. There is no photometry in that statement: measure the speed and the radius and you have weighed the interior. The first figure fixes this balance with a real number — the Sun's orbit around the Milky Way.

Figure

A star on a circular orbit at radius R feels only the mass interior to it, written M_enc. Balancing the inward gravity against the centripetal term v_c²/R gives v_c² = G·M_enc/R, so measuring v_c and R weighs everything inside. For the Sun, R = 8.2 kpc and v_c = 220 km/s give M_enc ≈ 9×10¹⁰ M⊙

Now walk outward through a real galaxy, tallying enclosed mass. If the mass stopped accumulating once you passed the visible disk — if all the matter were the matter that shines — gravity would weaken and orbital speeds would drop like $1/\sqrt{R}$ , as Neptune orbits the Sun more slowly than Mercury. Nature refuses: the speed stays constant, possible only if invisible mass keeps piling up in strict proportion to radius. Because fields superpose, the observed curve is the quadrature sum of each component's curve; the dark halo is the residual left after subtracting the best baryonic model, not a free-floating assumption. The second figure is the diagnostic itself — a decomposed curve in the style of the classic HI spiral NGC 3198.

Figure

Decomposed rotation curve of a spiral in the style of NGC 3198. The observed circular speed stays flat near 150 km/s while the stellar-disk contribution peaks near 2.2 disk scale lengths and declines; a dark halo rising toward 130 km/s fills the gap. Speeds add in quadrature, v_c² = v_disk² + v_gas² + v_halo². The baryons-only prediction (dashed) declines and falls far short

The gas curve pushes the measurement past the starlight: 21 cm emission from neutral hydrogen traces cold gas to radii where the stellar disk has faded into the sky background — where the dashed baryons-only prediction has peeled away below the flat observed curve, and dark matter becomes impossible to escape.

The core equations

The whole subject rests on one balance. For a spherically symmetric distribution only the interior mass contributes to the radial force, so setting the centripetal acceleration equal to the local gravity gives the master relation — exact for the bulge and halo, an excellent guide for the disk.

v_c^2(R)=\frac{GM(<R)}{R} \qquad (1)

where: $v_c$ — circular (rotation) speed at radius $R$ ; $M(<R)$ — mass enclosed within $R$ ; $G$ — Newton's gravitational constant; $R$ — galactocentric radius. A larger $v_c$ at fixed $R$ means more interior mass, so the flatness of $v_c$ is a statement about how mass accumulates. Inverting (1) is the entire point: kinematics deliver mass. Every dynamical weight of a galaxy flows from this rearrangement.

M(<R)=\frac{v_c^2\,R}{G}\;\approx\;2.33\times10^{5}\left(\frac{v_c}{\mathrm{km\,s^{-1}}}\right)^{2}\!\left(\frac{R}{\mathrm{kpc}}\right)M_\odot \qquad (2)

where: $M(<R)$ — dynamical mass within $R$ ; the numerical form packages $R/G$ into astronomer's units, so $200\,\mathrm{km\,s^{-1}}$ at $10\,\mathrm{kpc}$ encloses roughly $10^{11}M_\odot$ . Differentiating the enclosed mass gives the local density, so a rotation curve is also a direct probe of the density profile — how the inner slope that distinguishes cusps from cores is measured.

\rho(R)=\frac{1}{4\pi R^{2}}\frac{\mathrm{d}M(<R)}{\mathrm{d}R}=\frac{1}{4\pi G R^{2}}\frac{\mathrm{d}\!\left(v_c^{2}R\right)}{\mathrm{d}R} \qquad (3)

where: $\rho(R)$ — local mass density at radius $R$ ; $\mathrm{d}M/\mathrm{d}R$ — mass in a spherical shell per unit radius. A flat $v_c$ makes the derivative constant, forcing $\rho\propto R^{-2}$ . Because gravitational fields superpose, the total curve is the sum of component squares, not speeds. This quadrature rule is the backbone of every mass model — and why a declining disk and a rising halo can conspire into a flat total.

v_c^{2}(R)=v_\mathrm{bulge}^{2}(R)+v_\mathrm{disk}^{2}(R)+v_\mathrm{gas}^{2}(R)+v_\mathrm{halo}^{2}(R) \qquad (4)

where: $v_\mathrm{bulge}$ — central spheroid; $v_\mathrm{disk}$ — stellar disk, scaled by the stellar mass-to-light ratio $M/L$ ; $v_\mathrm{gas}$ — atomic plus molecular gas (scaled up ~1.4 for helium); $v_\mathrm{halo}$ — dark-matter halo. Each term is the curve that component would produce alone. The most-used cored halo is the pseudo-isothermal sphere: a constant-density core turning over into an $R^{-2}$ tail, giving a solid-body rise inside and a flat curve outside — the phenomenology of many dwarfs and low-surface-brightness disks.

\rho_\mathrm{iso}(R)=\frac{\rho_0 R_c^{2}}{R^{2}+R_c^{2}},\qquad v_c^{2}(R)=4\pi G\rho_0 R_c^{2}\left[1-\frac{R_c}{R}\arctan\!\frac{R}{R_c}\right] \qquad (5)

where: $\rho_0$ — central density; $R_c$ — core radius. For $R\gg R_c$ the bracket $\to1$ and $v_c\to\sqrt{4\pi G\rho_0R_c^2}$ ; for $R\ll R_c$ , $v_c\propto R$ . Cold-dark-matter simulations instead predict a cuspy profile — the Navarro–Frenk–White (NFW) form, with an $r^{-1}$ inner cusp and an $r^{-3}$ outer fall, and an analytic enclosed mass.

\rho_\mathrm{NFW}(r)=\frac{\rho_s}{(r/r_s)(1+r/r_s)^{2}},\qquad M(<r)=4\pi\rho_s r_s^{3}\left[\ln\!\Big(1+\frac{r}{r_s}\Big)-\frac{r/r_s}{1+r/r_s}\right] \qquad (6)

where: $\rho_s$ — characteristic density; $r_s$ — scale radius; concentration $c=r_{200}/r_s$ . The curve $v_c^2=GM(<r)/r$ rises, peaks near $\sim2.16\,r_s$ , and declines only slowly — flat over a broad range. Finally, the luminous disk is not spherical, so its true contribution comes from integrating the surface density. For the exponential disk $\Sigma(R)=\Sigma_0 e^{-R/R_d}$ , Freeman's closed form uses modified Bessel functions and, decisively, peaks and then falls.

v_c^{2}(R)=4\pi G\Sigma_0 R_d\,y^{2}\big[I_0(y)K_0(y)-I_1(y)K_1(y)\big],\qquad y\equiv\frac{R}{2R_d} \qquad (7)

where: $\Sigma_0$ — central surface density; $R_d$ — disk scale length; $I_n,K_n$ — modified Bessel functions of the first and second kind. The curve peaks near $R\approx2.2R_d$ and declines beyond — the disk alone can never make a flat curve.

Derivation AFlat rotation curve implies an inverse-square isothermal halo

A rotation curve that neither rises nor falls is an unusually strong constraint. If $v_c$ is independent of radius across a wide range, the enclosed mass and the density profile are not free — they are forced into a unique power-law form. The plan is to start from the spherical circular-speed relation, impose $v_c=\text{const}$ , and read off both $M(<R)$ and $\rho(R)$ , arriving at the singular isothermal sphere: the minimal halo reproducing a flat curve, and the reference profile against which every dark-matter model is compared.

Step 1. Rearrange the master relation (1) with $v_c$ held constant to get $M(<R)=v_c^2R/G$ : the enclosed mass grows linearly with radius, without bound, unless the halo is truncated at large $R$ . Why: a constant speed can only be held if the gravitating mass keeps pace with the growing radius, so the ratio $M/R$ stays fixed.

Step 2. Differentiate to find the mass in each spherical shell, $\mathrm{d}M/\mathrm{d}R=v_c^2/G$ , a constant — every shell of unit thickness adds the same mass. Why: equal mass per unit radius is the fingerprint of a flat curve, nothing like the exponentially declining light — the first hint the two do not trace the same stuff.

Step 3. Convert shell mass to density with $\mathrm{d}M=4\pi R^2\rho\,\mathrm{d}R$ , giving $\rho(R)=v_c^2/(4\pi G R^2)$ , falling exactly as $R^{-2}$ . Why: this is the singular isothermal sphere ( $v_c^2=2\sigma^2$ ); the baryons fall much faster than $R^{-2}$ , so the outer density must be dark.

Step 4. Check it the other way: integrating $\rho$ over a shell returns $\mathrm{d}M/\mathrm{d}R=v_c^2/G$ , and feeding $M(<R)$ back into (1) returns a constant $v_c$ . The loop closes. Why: consistency guarantees no algebra was lost, and the profile matches Binney & Tremaine §2.2. *The round-trip is self-consistent: $M(<R)=v_c^2R/G$ and $\rho=v_c^2/(4\pi GR^2)$ feed back into (1) and return a constant $v_c$ . Numerically, $v_c=220\,\mathrm{km\,s^{-1}}$ at $8.2\,\mathrm{kpc}$ gives $M\approx9\times10^{10}M_\odot$ and $\rho\approx0.013\,M_\odot\,\mathrm{pc^{-3}}$ .*

Go deeper · Derivation B · 1 minThe dark-matter fraction from a decomposed curve

Derivation BThe dark-matter fraction from a decomposed curve

Because gravitational fields superpose, the observed rotation curve is the quadrature sum of the curves each component would produce alone. This lets us turn a measured curve plus a photometric baryon model into a radius-by-radius accounting of how much enclosed mass is dark. The goal is to derive the dark fraction $f_\mathrm{dark}(R)$ from squared velocities alone, and show it climbs toward unity wherever the baryonic contribution sinks below the flat total — the quantitative form of "dark-matter dominated in the outskirts."

Step 1. At fixed radius every spherical component obeys $M_i(<R)=v_{c,i}^2R/G$ , and the total obeys $M(<R)=v_c^2R/G$ . Why: the shell theorem applies component-by-component, so each mass shares the common factor $R/G$ .

Step 2. Take the ratio $M_i(<R)/M(<R)=v_{c,i}^2/v_c^2$ ; the $R/G$ cancels, so mass fractions equal squared-velocity fractions. Why: this makes the decomposition measurable — you need only the component speeds, never $R$ or $G$ .

Step 3. Group the luminous contributions into $v_\mathrm{bar}^2=v_\mathrm{bulge}^2+v_\mathrm{disk}^2+v_\mathrm{gas}^2$ , so the halo is the remainder $v_\mathrm{halo}^2=v_c^2-v_\mathrm{bar}^2$ . Why: the halo is defined operationally as whatever the curve requires beyond the counted baryons — a residual, not an input.

Step 4. The dark fraction is $f_\mathrm{dark}(R)=v_\mathrm{halo}^2/v_c^2=1-v_\mathrm{bar}^2/v_c^2$ . As $R$ grows, $v_\mathrm{bar}$ declines while $v_c$ stays flat, so $f_\mathrm{dark}\to1$ . Why: it equals $0$ where baryons account for all rotation and $1$ where they contribute nothing, matching best-fit decompositions in Binney & Tremaine Ch. 6. *As a check: $f_\mathrm{dark}=1-(v_\mathrm{bar}/v_c)^2$ ; with $v_\mathrm{bar}=120$ and $v_c=200\,\mathrm{km\,s^{-1}}$ this gives $0.64$ , and with $v_\mathrm{bar}=90$ it gives $0.80$ — confirming that the dark fraction rises outward.*

Go deeper · Derivation C · 1 minTwo halos, one flat curve: the isothermal core versus the NFW cusp

Derivation CTwo halos, one flat curve: the isothermal core versus the NFW cusp

The two standard halo laws are nearly indistinguishable where the data are best and disagree most where they are worst — the very center. Deriving each rotation curve from its profile shows the pseudo-isothermal sphere rising as a solid body, $v_c\propto R$ , out of a constant-density core, while the NFW halo rises only as $v_c\propto\sqrt{r}$ out of its $\rho\propto r^{-1}$ cusp. The gap between these inner slopes is the entire content of the cusp–core problem, one of the sharpest confrontations between simulations and observed dwarfs.

Step 1. Expand the pseudo-isothermal mass (5) at small radius using $\arctan x\approx x-x^3/3$ : $M(<R)\to\tfrac{4\pi}{3}\rho_0R^3$ . This is a uniform sphere, so $v_c^2=GM/R=\tfrac{4\pi}{3}G\rho_0R^2$ and $v_c\propto R$ . Why: a constant-density core has no central concentration, so the inner rotation is solid-body — the smoking-gun shape seen in many dwarfs.

Step 2. Take the same halo at large radius, where $\rho\to\rho_0R_c^2/R^2$ and the bracket in (5) tends to $1$ : $v_c\to\sqrt{4\pi G\rho_0R_c^2}=\text{const}$ . Why: the cored halo recovers the isothermal flat curve far out, so core and cusp look identical in the outskirts.

Step 3. Now expand the NFW mass (6) at small $x=r/r_s$ with $\ln(1+x)-x/(1+x)\approx x^2/2$ : $M(<r)\to2\pi\rho_s r_s\,r^2$ , so $v_c^2=GM/r\propto r$ and $v_c\propto\sqrt{r}$ . Why: the $\rho\propto r^{-1}$ cusp packs more mass at the very center than a core, so the NFW curve rises faster than solid-body — a $\sqrt{r}$ rise, not a straight $\propto R$ .

Step 4. Both profiles flatten across the observed disk; NFW peaks near $2.16\,r_s$ and declines gently outside. The discriminant is entirely the innermost slope of $v_c(R)$ , which is why high-resolution inner curves — not the flat outskirts — decide the debate. Why: the outer curve is degenerate between halo and disk $M/L$ , so profile information lives where the curve is hardest to measure. *The results are internally consistent: the pseudo-isothermal profile gives $v_c\propto R$ in the core and $v_c\to\mathrm{const}$ outside; the NFW profile gives $v_c\propto\sqrt{r}$ in the cusp with a peak near $2.16\,r_s$ ; and the series expansion $\ln(1+x)-x/(1+x)=x^2/2-\tfrac{2}{3}x^3+\cdots$ underpins the small- $x$ limits used above.*

Limiting cases

The circular-speed relation collapses to a handful of readable regimes, each a shape you can spot by eye on a measured curve.

Case	Result	Meaning
Uniform-density core (solid body)	$v_c(R)\propto R$	Inner regions with $\rho\approx\text{const}$ give $M(<R)\propto R^3$ ; rotation rises linearly, as in the innermost parts of many disks.
All mass enclosed (Keplerian)	$v_c(R)=\sqrt{GM/R}\propto R^{-1/2}$	Outside a point mass or fully enclosed distribution — the decline expected if a galaxy ended at its starlight.
Flat rotation curve	$v_c=\text{const}\Rightarrow M(<R)=\dfrac{v_c^2R}{G},\ \rho=\dfrac{v_c^2}{4\pi GR^2}$	The observed large-radius regime: mass grows linearly, density falls as $R^{-2}$ — the isothermal signature.
Thin exponential disk alone	$v_c^2=4\pi G\Sigma_0R_d\,y^2[I_0K_0-I_1K_1],\ y=\tfrac{R}{2R_d}$	Peaks near $R\approx2.2R_d$ , then declines; cannot by itself make a flat curve — the Freeman result.

Common trapA flat curve can be read straight off a spherical formula for a disk

Two related traps hide here. First, a real stellar disk is flattened, not spherical, so the shell theorem does not apply exactly and $v_c^2=GM(<R)/R$ is only a guide for the luminous disk — its honest contribution is the Freeman/Bessel result (7), which peaks near $2.2$ scale lengths and then declines. That decline is not a measurement failure; it is precisely why the observed flatness cannot come from the disk alone and genuinely requires extra mass. Second, do not add component speeds linearly. Fields superpose, so the squares add: $v_c^2=\sum_i v_{c,i}^2$ . Writing $v_c=v_\mathrm{disk}+v_\mathrm{halo}$ instead of $v_c^2=v_\mathrm{disk}^2+v_\mathrm{halo}^2$ overstates the total and mis-apportions the dark mass at every radius.

Worked examples

1 · Weighing the Milky Way inside the solar circle (fully worked)

The Sun orbits at $R_0=8.2\,\mathrm{kpc}$ with a local circular speed $v_c=220\,\mathrm{km\,s^{-1}}$ . How much mass lies interior to the Sun, and what does the flat outer curve imply for the Galaxy's total mass?

Step 1. Use the enclosed-mass form (2). The compact astronomer's version, $M(<R)\approx2.33\times10^{5}(v_c/\mathrm{km\,s^{-1}})^2(R/\mathrm{kpc})\,M_\odot$ , already folds in $R/G$ , so no SI conversion is needed. Every ingredient is a pure number here.

Step 2. Substitute: $M(<R_0)=2.33\times10^{5}\times(220)^2\times8.2$ . With $48\,400\times8.2=3.97\times10^5$ , this is $2.33\times10^{5}\times3.97\times10^5=\mathbf{9.2\times10^{10}\,M_\odot}$ . A full SI calculation with $G=6.674\times10^{-11}$ gives $9.23\times10^{10}M_\odot$ — the two agree.

Step 3. Interpret. The Galaxy's stellar mass is roughly $5$ – $6\times10^{10}M_\odot$ , so at the Sun the dynamical and stellar masses are comparable — baryons still dominate here. But the curve stays flat to tracers tens of kpc out (globular clusters, halo stars, satellites), and $M\propto R$ drives the total toward $\sim10^{12}M_\odot$ by $\sim200\,\mathrm{kpc}$ . Interior to the Sun, $M\approx9\times10^{10}M_\odot$ ; the full halo is roughly ten times larger and overwhelmingly dark.

2 · The dark-matter fraction, and how it grows outward

At an inner radius the observed circular speed is $v_c=200\,\mathrm{km\,s^{-1}}$ and the baryonic mass model predicts $v_\mathrm{bar}=120\,\mathrm{km\,s^{-1}}$ . Farther out the curve is still flat at $200\,\mathrm{km\,s^{-1}}$ but the baryons have dropped to $v_\mathrm{bar}=90\,\mathrm{km\,s^{-1}}$ . Find the dark fraction at both radii.

Step 1. Apply $f_\mathrm{dark}=1-(v_\mathrm{bar}/v_c)^2$ from the decomposition derivation. Inner: $(120/200)^2=0.36$ , so $f_\mathrm{dark}=1-0.36=\mathbf{0.64}$ . About two-thirds of the enclosed mass is already dark.

Step 2. Outer: $(90/200)^2=0.2025$ , so $f_\mathrm{dark}=1-0.2025\approx\mathbf{0.80}$ . As $v_\mathrm{bar}$ falls while $v_c$ holds flat, the halo's share climbs. The galaxy goes from baryon-supported to dark-matter dominated outward — $f_\mathrm{dark}$ marches from $0.64$ toward $1$ .

3 · The local dark-matter density from a singular isothermal sphere

If the Sun's entire $220\,\mathrm{km\,s^{-1}}$ rotation at $8.2\,\mathrm{kpc}$ were produced by a singular isothermal halo, what local dark-matter density would that imply?

Step 1. Use the flat-curve density law $\rho=v_c^2/(4\pi GR^2)$ from Derivation A, with $v_c=2.20\times10^5\,\mathrm{m\,s^{-1}}$ and $R=8.2\,\mathrm{kpc}=2.53\times10^{20}\,\mathrm{m}$ . This is an upper bound: in reality much of the inner rotation is baryonic, not halo.

Step 2. Compute $\rho=(2.20\times10^5)^2/[4\pi(6.674\times10^{-11})(2.53\times10^{20})^2]=9.0\times10^{-22}\,\mathrm{kg\,m^{-3}}$ . Converting with $1\,\mathrm{pc}^3=2.94\times10^{49}\,\mathrm{m}^3$ gives $\rho\approx\mathbf{0.013\,M_\odot\,pc^{-3}}$ . This overestimate lands within a factor of a few of the measured local dark-matter density, $\sim0.01\,M_\odot\,\mathrm{pc^{-3}}$ ( $\approx0.4\,\mathrm{GeV\,cm^{-3}}$ ) — reassuring for a one-line model.

Physics of the portrait

The rotation curve carries a name because one astronomer, more than anyone, made the flat curve undeniable — through a relentless, galaxy-by-galaxy campaign that left no room for doubt.

Origins · 1 minVera Rubin and how this was discovered

Vera Rubin

1928–2016 · flat rotation curves and the dynamical case for dark matter

Vera Rubin spent a career being right before it was fashionable. Denied access to Palomar until 1965, she teamed with instrument-builder Kent Ford, whose sensitive image-tube spectrograph could measure Doppler shifts of faint Hα emission far out in a galaxy's disk. Beginning with Andromeda (M31) around 1970 and sweeping through dozens of Sc and Sb spirals, Rubin and Ford found the same thing every time: the rotation speed climbed through the inner disk and then refused to fall, staying flat to the outermost point they could reach. Their 1980 synthesis of 21 bright spirals made the result impossible to dismiss as a quirk of one object. Rubin herself was cautious — she reported curves, not conclusions — but the pattern demanded either enormous unseen mass or a change to gravity itself. The dark-matter reading, foreshadowed by Fritz Zwicky's 1933 cluster masses, sharpened by Kenneth Freeman's 1970 note on HI disks, and cemented by Albert Bosma's 21 cm curves and the Ostriker–Peebles argument that disks need massive halos to stay stable, became the consensus. Her optical curves and the radio HI curves reinforced each other, turning dark matter from a cluster-scale curiosity into a property of ordinary galaxies. She never won the Nobel her data arguably earned; the flat curve is her monument.

History

1933

— Zwicky infers large unseen mass in the Coma cluster from galaxy velocities

1939

— Babcock measures a rising, non-Keplerian rotation curve for Andromeda

1970

— Freeman notes HI curves of M33 and NGC 300 imply undetected mass; Rubin and Ford publish the M31 optical curve

1978

— Bosma's 21 cm HI curves stay flat well beyond the starlight

1980

— Rubin, Ford and Thonnard establish flat curves across 21 Sc spirals

1996

— Navarro, Frenk and White propose the universal cold-dark-matter halo profile

Going further

Texts. Binney & Tremaine, Galactic Dynamics, 2nd ed., §2.2 (spherical potentials), §2.6 (disk potentials and the Freeman formula), Ch. 6 (mass modeling) · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 24 (Milky Way rotation, dark matter) · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 11 (halo profiles, NFW) · Sparke & Gallagher, Galaxies in the Universe, 2nd ed., Ch. 5 · Ryden, Introduction to Cosmology, 2nd ed., Ch. 8 · Dodelson, Modern Cosmology, Ch. 1. The key relations used throughout — $v_c^2=GM(<R)/R$ , the isothermal density law $\rho=v_c^2/(4\pi GR^2)$ , the pseudo-isothermal curve $v_c^2=4\pi G\rho_0R_c^2[1-(R_c/R)\arctan(R/R_c)]$ , the NFW enclosed mass $M(<r)$ , and the Freeman Bessel-function disk curve — are standard results. The Milky Way numbers adopted here ( $R_0\approx8.2\,\mathrm{kpc}$ , $v_c\approx220$ – $240\,\mathrm{km\,s^{-1}}$ , $M(<R_0)\approx9\times10^{10}M_\odot$ , total halo $\sim10^{12}M_\odot$ ) are standard textbook ranges. The unit shortcut $M\approx2.33\times10^5\,v^2R\,M_\odot$ and the local isothermal density $0.013\,M_\odot\,\mathrm{pc^{-3}}$ follow directly from the defining constants. The dark-fraction example uses only its stated inputs.

Prerequisites. Gravitational Two-Body Problem · Galactic Potential Models · Orbits in Spherical Potentials. Leads to. Jeans Equations · Virial Theorem in Stellar Systems · Distribution Function and Phase Space.

Ask the tutor: Want to fit a real decomposed curve, derive the Freeman disk peak, or see how baryonic feedback turns an NFW cusp into a core? Ask and we will work it through.

Why do the component speeds add in quadrature instead of just adding up?

Because gravity, not velocity, is the thing that superposes. Each mass component sources its own field, and fields add as vectors: the total radial acceleration is the sum of the individual ones. Since a circular orbit ties acceleration to speed through $v_c^2/R=g$ , each component contributes $v_{c,i}^2/R$ , and summing those and multiplying by $R$ gives $v_c^2=\sum_i v_{c,i}^2$ . The speeds themselves never add — only their squares do. Adding speeds linearly double-counts the acceleration and badly overstates the total, so the quadrature rule is the difference between a correct mass model and a nonsensical one.

Does a flat rotation curve mean the galaxy has infinite mass?

Taken literally, an exactly flat curve extending to infinity would imply $M(<R)\propto R$ growing without bound — an infinite mass. Real halos are not infinite: they are truncated near the virial radius, where the density steepens (the NFW profile falls as $r^{-3}$ ) and the curve turns down. The flatness is an excellent description across the observable disk and the near halo, but it does not continue forever. What it robustly establishes is that far more mass is enclosed at large radius than the light accounts for — not that the mass is literally unbounded. The truncation is where the halo's total mass and virial radius are defined.

What is the cusp–core problem, and is it settled?

Cold-dark-matter simulations of collisionless particles predict halos with a central density cusp, $\rho\propto r^{-1}$ (the NFW form), making the inner rotation curve rise as $v_c\propto\sqrt{r}$ . Many observed dwarf and low-surface-brightness galaxies instead show constant-density cores, with a solid-body rise $v_c\propto R$ . That mismatch is the cusp–core problem, and it is not fully settled. Proposed resolutions include baryonic feedback — repeated gas inflows and supernova-driven outflows that dynamically heat the central dark matter and flatten a cusp into a core — as well as genuinely different dark matter (self-interacting or warm). The modern framing emphasizes the diversity of inner slopes at fixed mass, which any successful theory must reproduce.

Is modified gravity (MOND) a viable alternative to a dark halo here?

It is the most serious alternative, and rotation curves are its home turf. Modified Newtonian dynamics posits that below a characteristic acceleration $a_0\approx1.2\times10^{-10}\,\mathrm{m\,s^{-2}}$ the effective gravity strengthens, reproducing flat curves and the radial acceleration relation from the baryons alone, with no halo. Within individual galaxies MOND is impressively economical. Its difficulties appear elsewhere: galaxy clusters still need extra unseen mass even in MOND, and the cosmic microwave background and large-scale structure are described extraordinarily well by cold dark matter and poorly by MOND without added fields. Most cosmologists therefore read the flat curve as dark matter while treating the radial acceleration relation as a real regularity the dark picture must explain.

Why use 21 cm hydrogen gas rather than stars to measure the outer curve?

Stars in a spiral disk fade exponentially with radius, so optical kinematics run out of signal not far past the bright disk. Neutral atomic hydrogen, by contrast, extends in a thin gaseous disk to two or three times the optical radius and emits the 21 cm hyperfine line whether or not any stars are nearby. That line gives a clean line-of-sight velocity, so radio interferometers can trace $v_c(R)$ far into the region where the dark halo dominates — exactly where the curve stays flat and the baryons have all but vanished. Historically it was combining optical Hα curves inside with 21 cm curves in the outskirts that made the flatness, and thus the dark matter, undeniable.

What is the maximum-disk degeneracy, and how do you break it?

In the inner galaxy the stellar disk and the halo trade off almost freely: raise the stellar mass-to-light ratio and the disk explains more of the curve, needing a lighter halo; lower it and the halo must do more. Because both fit the same total curve, the decomposition is degenerate, and the "maximum-disk" solution deliberately pushes the stellar contribution as high as the curve allows. Breaking the degeneracy needs information the rotation curve alone lacks: stellar-population synthesis or colours to pin $M/L$ , vertical stellar velocity dispersions to weigh the disk surface density directly, or the outer curve where the halo dominates unambiguously. Without such a handle, the inner disk/halo split stays uncertain.