Dark-Matter Halos

Weigh a spiral galaxy by clocking how fast its gas orbits, and something refuses to add up. Circular balance says $M(<r)=v_c^2 r/G$ , so once you run out of stars the orbital speed should fall away Keplerian, as $r^{-1/2}$ . It does not. Across dozens of kiloparsecs the rotation speed stays stubbornly flat, which means mass keeps accumulating, $M(<r)\propto r$ , where there is essentially no light. The same failure appears at every scale: cluster velocity dispersions and X-ray temperatures demand ten times more gravitating mass than the baryons supply, and gravitational lensing draws mass maps that dwarf the visible galaxies threading them. The name for the missing, dominant, diffuse mass is a dark-matter halo. Dynamically a halo is a collisionless self-gravitating system. Its two-body relaxation time exceeds the age of the Universe by roughly ten orders of magnitude, so individual particles never scatter; the matter is a smooth phase-space fluid $f(\mathbf{x},\mathbf{v},t)$ obeying the Collisionless Boltzmann Equation, its moments the Jeans equations, its energy budget fixed by the Virial theorem $2T+W=0$ rather than by pressure. Cosmologically, halos are the nonlinear end products of gravitational instability: tiny fluctuations grow, collapse, and — in the cold-dark-matter picture — merge hierarchically, small halos first, while violent relaxation drives the survivors toward a nearly universal shape, the Navarro–Frenk–White profile. This article develops the observational diagnostics, the canonical analytic profiles, the virial and concentration bookkeeping, and the equilibrium machinery tying density to velocity structure.

Key takeaways

A dark-matter halo is a quasi-equilibrium, self-gravitating cloud of collisionless dark matter inside which a galaxy condenses; it dominates the mass on kiloparsec-to-megaparsec scales.
The evidence is kinematic: flat rotation curves, cluster velocity dispersions, X-ray hydrostatics, and gravitational lensing all demand far more mass than the light.
Everything follows from one relation — $v_c^2(r)=GM(<r)/r$ — so a flat curve forces $M(<r)\propto r$ and hence $\rho\propto r^{-2}$ locally.
In CDM simulations relaxed halos follow the near-universal NFW profile: an inner cusp $\rho\propto r^{-1}$ steepening to an outer tail $\rho\propto r^{-3}$ .
A halo is summarized by three numbers — virial mass $M_{200}$ , virial radius $r_{200}$ , and concentration $c=r_{200}/r_s$ — with $c$ encoding formation epoch.
The workhorse tools are collisionless: the Boltzmann and Jeans equations and the Virial theorem, not gas pressure, hold a halo up.
Key idea: a dark-matter halo is a collisionless self-gravitating equilibrium whose density profile — universally close to NFW in CDM — is read off through $v_c^2(r)=GM(<r)/r$ , so a flat rotation curve directly demands $M(<r)\propto r$ , i.e. $\rho\propto r^{-2}$ in the outskirts.

The core insight

A dark-matter halo is a collisionless self-gravitating equilibrium whose density profile — universally close to NFW in CDM — is probed observationally through the enclosed-mass relation $v_c^2(r)=GM(<r)/r$ ; flat rotation curves demand $M(<r)\propto r$ , i.e. $\rho\propto r^{-2}$ in the outer regions.

Where this lives today

Halo structure is a live frontier at both ends of the mass scale. In galaxy surveys, weak gravitational lensing — the coherent one-percent distortion of background-galaxy shapes by foreground mass — is now measured for millions of lenses by DES, KiDS, HSC, and soon Euclid and Rubin/LSST, giving stacked mass profiles that test NFW and pin the concentration–mass relation. At cluster scales, lensing, X-ray hydrostatics from Chandra and XMM, and the Sunyaev–Zel'dovich effect together weigh $10^{14}$ – $10^{15}\,M_\odot$ halos, whose abundance versus mass and redshift is a sharp handle on structure growth and neutrino mass. The most contested arena is the smallest: dwarf-galaxy rotation curves often prefer cored inner profiles ( $\rho\to$ const) over the NFW cusp — the "core–cusp problem" — and the count of luminous satellites sits below the CDM subhalo count. Whether these are resolved by baryonic feedback flattening cusps, or by new physics — warm, self-interacting, or ultralight dark matter — is being tested with stellar streams, strong-lens substructure, and Milky-Way satellite censuses from Gaia and DES.

The physical picture

Put a test particle on a circular orbit at radius $r$ in a spherical mass distribution. The interior mass supplies the centripetal acceleration, $v_c^2/r = GM(<r)/r^2$ , so the orbital speed reads off the mass inside the orbit: $M(<r)=v_c^2 r/G$ . Now trace $v_c$ outward with 21-cm hydrogen or optical emission lines. Inside the disk the curve rises; past the last stars, if the light were all the mass, $v_c$ would peel off as $r^{-1/2}$ — the Keplerian decline of a point mass. Real spirals do the opposite: the curve flattens and stays flat, sometimes to ten disk scale-lengths. A flat $v_c$ means $M(<r)\propto r$ , so mass is still piling up where no light is emitted. That stubborn fact is the observational backbone of dark matter, and the first figure is the diagnostic an observer plots.

Figure

Flat rotation curve: the total circular speed stays near v_c ≈ 200 km/s far past the stellar disk, so the enclosed mass M(<r) = v_c² r / G grows linearly with radius instead of falling; the amber baryon-only curve drops off Keplerian while the green dark-matter halo term rises to hold the total flat

A halo is held up by orbits, not pressure. Because encounters are negligible, its constituents stream ballistically through the mean field; equilibrium means the phase-space density is stationary, not that particles have thermalized. This is why the governing statements are the Collisionless Boltzmann Equation, its velocity moments the Jeans equations, and the Virial theorem $2T+W=0$ . The density profile is a fossil of the collapse: cold, nearly pressureless matter relaxes violently in the rapidly varying potential of a forming halo, erasing memory of the initial conditions. What survives is a near-universal shape — a steep inner cusp, a shallower-then-steeper outer decline — captured by NFW, $\rho\propto r^{-1}(1+r/r_s)^{-2}$ . The second figure shows it on log–log axes, where its three regimes and the isothermal idealization become straight lines.

Figure

The NFW density profile on log–log axes, ρ = ρ_s / ((r/r_s)(1+r/r_s)²): an inner cusp ρ ∝ r⁻¹, logarithmic slope −2 at the scale radius r_s, and an outer tail ρ ∝ r⁻³. The amber straight line is the singular isothermal sphere ρ ∝ r⁻² that gives an exactly flat rotation curve; for concentration c = r₂₀₀/r_s = 10 the virial edge r₂₀₀ sits at r/r_s = 10

The most compact single descriptor of a halo is its concentration $c=r_{200}/r_s$ — how tightly mass is packed toward the center. Halos that assembled early collapsed when the mean cosmic density was higher, so they are more concentrated; massive halos, still assembling today, are less so. Concentration is thus a clock, which is why the mean $c(M,z)$ relation from simulations is a cosmological probe.

The core equations

The keystone is the circular-speed–mass relation. By the shell theorem only interior mass pulls, so the field at $r$ is $GM(<r)/r^2$ ; balancing it against the centripetal $v_c^2/r$ turns an orbital-speed measurement directly into an enclosed mass.

v_c^{2}(r)=\frac{G\,M(<r)}{r},\qquad M(<r)=\frac{v_c^{2}(r)\,r}{G} \qquad (1)

where: $v_c(r)$ — circular speed at radius $r$ ; $M(<r)$ — total gravitating mass within $r$ ; $G$ — Newton's constant. Reading it: a flat curve, $v_c=\text{const}$ , forces $M(<r)\propto r$ , and differentiating gives a local slope $\rho\propto r^{-2}$ . The shape $N$ -body simulations converge on is the Navarro–Frenk–White profile: fixed exponents and two scales — a density $\rho_s$ and a radius $r_s$ — with the logarithmic slope passing through $-2$ exactly at $r_s$ .

\rho_\mathrm{NFW}(r)=\frac{\rho_s}{(r/r_s)\,(1+r/r_s)^{2}} \qquad (2)

where: $\rho_s$ — characteristic density; $r_s$ — scale radius (slope $-2$ ); $r$ — radius from center. Reading it: $\rho\propto r^{-1}$ for $r\ll r_s$ (cusp) and $\rho\propto r^{-3}$ for $r\gg r_s$ (tail) — one formula interpolating the two power laws in the figure. Integrating in spherical shells gives a closed-form enclosed mass — what turns a fitted profile into a predicted rotation curve.

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

where: $M_\mathrm{NFW}(<r)$ — NFW mass within $r$ ; $\rho_s,\,r_s$ — the two NFW scales. Reading it: the bracket grows only logarithmically, so the total mass formally diverges and the halo needs a finite outer boundary. That boundary is a cosmological convention: the virial radius encloses a fixed overdensity — commonly $200$ times the critical density — and the virial mass is the mass inside it.

M_{200}=\frac{4}{3}\pi r_{200}^{3}\,\big(200\,\rho_\mathrm{crit}\big),\qquad \rho_\mathrm{crit}=\frac{3H^{2}}{8\pi G} \qquad (4)

where: $M_{200}$ — virial mass; $r_{200}$ — radius of mean interior density $200\rho_\mathrm{crit}$ ; $\rho_\mathrm{crit}$ — critical density; $H$ — Hubble parameter at the halo's epoch. Reading it: $r_{200}$ is a bookkeeping edge, not a physical surface, and a halo mass is ambiguous without its overdensity definition. Combining the virial definition with the NFW mass fixes the characteristic density in terms of one number — the concentration.

c\equiv\frac{r_{200}}{r_s},\qquad \rho_s=\frac{200}{3}\,\rho_\mathrm{crit}\,\frac{c^{3}}{\ln(1+c)-c/(1+c)} \qquad (5)

where: $c$ — concentration (typically $5$ – $20$ , higher for low-mass or early halos); $\rho_s$ — NFW density set by $c$ and cosmology. Reading it: $\rho_s\propto c^3/[\ln(1+c)-c/(1+c)]$ climbs steeply with $c$ , so a more concentrated halo is far denser at its scale radius. The bridge from density to observable velocities is the spherical Jeans equation — the radial moment of the Collisionless Boltzmann Equation for a steady spherical system.

\frac{1}{\rho}\frac{\mathrm{d}(\rho\,\sigma_r^{2})}{\mathrm{d}r}+\frac{2\beta\,\sigma_r^{2}}{r}=-\frac{\mathrm{d}\Phi}{\mathrm{d}r}=-\frac{G M(<r)}{r^{2}} \qquad (6)

where: $\sigma_r$ — radial velocity dispersion; $\beta=1-\sigma_t^2/(2\sigma_r^2)$ — velocity anisotropy; $\Phi$ — halo potential. Reading it: the dispersion (a "pressure") and its anisotropy support the halo against gravity — hydrostatic equilibrium for a collisionless fluid. Two closed-form limits anchor the scheme: the singular isothermal sphere, the scale-free equilibrium giving a perfectly flat curve, and the Einasto profile, the mild NFW refinement high-resolution simulations prefer.

\rho_\mathrm{iso}(r)=\frac{\sigma^{2}}{2\pi G\,r^{2}},\quad v_c=\sqrt{2}\,\sigma;\qquad \rho_\mathrm{Ein}(r)=\rho_{-2}\exp\!\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^{\alpha}-1\right]\right\} \qquad (7)

where: $\sigma$ — constant one-dimensional dispersion; $\alpha\approx0.16$ – $0.20$ — Einasto shape index; $r_{-2},\rho_{-2}$ — radius and density at slope $-2$ . Reading it: the isothermal sphere approximates any halo where the slope passes $-2$ , while Einasto rounds off the NFW cusp with a gently curving log-slope.

Derivation AFlat rotation curves force an isothermal halo

The goal is to run the observation backward: take the one measured fact that the circular speed is constant over the halo-dominated regime, and deduce what density profile is required to produce it, together with the velocity dispersion of the matter that carries it. We convert a constant $v_c$ into an enclosed mass growing linearly with radius, differentiate to recover the local density, and match to the singular isothermal sphere — the chain that makes "flat curve" and "inverse-square density" the same statement.

Step 1. From equation (1) with $v_c(r)=v_0$ constant, the enclosed mass is $M(<r)=v_0^2 r/G$ , growing linearly with radius. Why: a flat curve is not "no more mass" — it is mass accumulating at exactly the rate needed to keep $GM/r$ fixed as $r$ increases.

Step 2. Differentiate the shell relation $\mathrm{d}M=4\pi r^2\rho\,\mathrm{d}r$ , giving $4\pi r^2\rho(r)=\mathrm{d}M/\mathrm{d}r=v_0^2/G$ , a constant. Why: the mass in each spherical shell must be independent of radius, which pins down how fast the density falls.

Step 3. Solve for the density, $\rho(r)=v_0^2/(4\pi G r^2)$ — an inverse-square profile, the density of a singular isothermal sphere. Why: only $\rho\propto r^{-2}$ keeps shell mass constant, so flat curve and isothermal density are equivalent over the observed range.

Step 4. Compare with $\rho_\mathrm{iso}=\sigma^2/(2\pi G r^2)$ ; the two agree when $v_0^2=2\sigma^2$ , i.e. $v_c=\sqrt{2}\,\sigma$ . Why: this ties the orbital speed to the velocity spread, so the flat-curve amplitude measures the halo's dispersion.

Go deeper · Derivation B · 1 minIntegrating the NFW profile and fixing its density

Derivation BIntegrating the NFW profile and fixing its density

The goal is twofold: to evaluate the mass of an NFW halo in closed form, and then to remove one of its two normalizations by demanding consistency with the virial definition, so the whole profile is set by a single number, the concentration. We perform the shell integral by nondimensionalizing on $r_s$ , evaluate the NFW mass at the virial radius, equate it to the overdensity definition of $M_{200}$ , and solve for $\rho_s$ — turning NFW from a two-parameter fit into a one-parameter, cosmology-anchored family.

Step 1. Insert the NFW density into $M(<r)=4\pi\int_0^r r'^2\rho\,\mathrm{d}r'$ and substitute $x\equiv r'/r_s$ ; one power of $x$ cancels, leaving $M=4\pi\rho_s r_s^3\int_0^{r/r_s} x\,(1+x)^{-2}\,\mathrm{d}x$ . Why: scaling out $r_s$ makes the mass depend on radius only through the ratio $r/r_s$ .

Step 2. Evaluate by writing $x=(1+x)-1$ : $\int_0^{y} x(1+x)^{-2}\mathrm{d}x=\ln(1+y)-y/(1+y)$ , giving equation (3). Why: the logarithm is the fingerprint of the outer $r^{-3}$ tail — mass grows without bound, if only logarithmically.

Step 3. Evaluate at $r=r_{200}=c\,r_s$ , so $M_{200}=4\pi\rho_s r_s^3\,m(c)$ with $m(c)\equiv\ln(1+c)-c/(1+c)$ . Why: this is the NFW mass out to the adopted edge, one handle on $M_{200}$ .

Step 4. Set it equal to the overdensity definition $M_{200}=\tfrac43\pi c^3 r_s^3(200\rho_\mathrm{crit})$ ; the $r_s^3$ cancels and $\rho_s=\tfrac{200}{3}\rho_\mathrm{crit}\,c^3/m(c)$ , equation (5). Why: the cancellation makes $\rho_s$ purely a function of $c$ and cosmology, so equal- $c$ halos share it regardless of mass — the origin of NFW universality. *The characteristic overdensity is $\delta_c=\rho_s/\rho_\mathrm{crit}=(200/3)c^3/m(c)$ ; for $c=10$ this evaluates to $\delta_c\approx4.5\times10^4$ .*

Go deeper · Derivation C · 1 minThe isothermal sphere from the isotropic Jeans equation

Derivation CThe isothermal sphere from the isotropic Jeans equation

The goal is to show that the inverse-square halo is not merely consistent with a flat rotation curve but is an exact, self-supporting equilibrium of the collisionless dynamics — a genuine steady solution of the spherical Jeans equation. We specialize the Jeans equation to isotropic orbits and constant dispersion, guess the scale-free power law a flat curve suggests, compute its enclosed mass, and verify that its pressure gradient exactly balances its own gravity, fixing the remaining constant.

Step 1. In equation (6) set the anisotropy $\beta=0$ and take $\sigma_r=\sigma$ constant; the support term collapses to $\sigma^2\,\mathrm{d}\ln\rho/\mathrm{d}r=-GM(<r)/r^2$ . Why: isotropy kills the $2\beta\sigma_r^2/r$ term and constancy pulls $\sigma^2$ out, leaving hydrostatic balance for an isothermal collisionless gas.

Step 2. Try the scale-free ansatz $\rho=\rho_0(r/r_0)^{-2}$ , for which $\mathrm{d}\ln\rho/\mathrm{d}r=-2/r$ . Why: we test whether the flat-curve slope $\rho\propto r^{-2}$ is a true self-consistent solution, not just a local slope.

Step 3. Its enclosed mass is $M(<r)=4\pi\int_0^r\rho_0(r'/r_0)^{-2}r'^2\,\mathrm{d}r'=4\pi\rho_0 r_0^2\,r$ , linear in $r$ . Why: the $r^{-2}$ density keeps shell mass constant, so the running mass is linear — now from dynamics, not kinematics.

Step 4. Substitute into the balance: $\sigma^2(-2/r)=-4\pi G\rho_0 r_0^2/r$ , an identity provided $\rho_0 r_0^2=\sigma^2/(2\pi G)$ , i.e. $\rho(r)=\sigma^2/(2\pi G r^2)$ . Why: the radius dependence cancels, so balance holds at every radius — an exact scale-free equilibrium, recovering the isothermal sphere with $v_c=\sqrt2\,\sigma$ .

Limiting cases

Case	Result	Meaning
Inner NFW cusp	$\rho\;\xrightarrow{r\ll r_s}\;\rho_s r_s/r\propto r^{-1}$	Divergent but integrable central density; $M(<r)\propto r^2\to0$ , so $v_c\to0$ at the very center.
Outer NFW tail	$\rho\;\xrightarrow{r\gg r_s}\;\rho_s r_s^3/r^3\propto r^{-3}$	Mass diverges only logarithmically; total mass is formally infinite, forcing a finite $r_{200}$ .
Singular isothermal sphere	$\rho=\sigma^2/(2\pi G r^2),\;v_c=\sqrt2\,\sigma=\text{const}$	Exactly flat rotation curve; the intermediate-radius approximation where the NFW slope passes $-2$ .
Keplerian (no halo)	$v_c\;\xrightarrow{r>R_\star}\;\sqrt{GM_\star/r}\propto r^{-1/2}$	What one would see beyond the baryons if no halo existed; its absence is the primary evidence for dark matter.
High concentration	$\rho_s\propto c^3/\ln c$ grows steeply	Early-forming, low-mass halos are far denser at $r_s$ ; concentration is a clock for formation epoch.

Common trapA flat rotation curve means the halo is isothermal everywhere.

It does not. The flat portion of a curve constrains only the enclosed mass over the observed radii — it says $M(<r)\propto r$ there, hence $\rho\propto r^{-2}$ locally, and nothing about the profile inside or outside that range. Simulated NFW halos have $\rho\propto r^{-1}$ cusps and $\rho\propto r^{-3}$ tails yet still produce approximately flat curves across a decade in radius. A second trap: the "virial radius" $r_{200}$ is a convention (mean interior density $=200\rho_\mathrm{crit}$ ), not a physical edge. Quoting a halo mass without stating the overdensity — $200\rho_\mathrm{crit}$ , $200\bar\rho_m$ , or $\Delta_\mathrm{vir}\approx100$ – $200$ — is ambiguous and can shift the mass by tens of percent.

Worked examples

1 · Weighing a spiral from its flat curve (confidence builder)

A spiral's rotation curve is flat at $v_c=200\ \mathrm{km\,s^{-1}}$ out to $r=20\ \mathrm{kpc}$ . Find the enclosed mass, the implied mass-to-light ratio if the enclosed light is $L\approx2\times10^{10}\,L_\odot$ , and the local density slope.

Step 1. Use equation (1) in galaxy-friendly units, $G=4.301\times10^{-6}\ \mathrm{kpc\,(km\,s^{-1})^2}\,M_\odot^{-1}$ : $M(<r)=v_c^2 r/G=(200)^2(20)/(4.301\times10^{-6})=8.0\times10^{5}/4.301\times10^{-6}\ M_\odot$ . $M(<20\ \mathrm{kpc})\approx1.9\times10^{11}\,M_\odot$ .

Step 2. The mass-to-light ratio inside $20\ \mathrm{kpc}$ is $M/L\approx1.9\times10^{11}/2\times10^{10}\approx\mathbf{9\,M_\odot/L_\odot}$ , rising with radius because $M\propto r$ while the light has stopped. A stellar population has $M/L\sim1$ – $5$ ; a value climbing past $10$ is the dark halo announcing itself.

Step 3. A flat curve means $M(<r)\propto r$ , so $\mathrm{d}M/\mathrm{d}r=4\pi r^2\rho$ is constant and $\rho\propto r^{-2}$ locally — the isothermal slope. The corresponding one-dimensional dispersion is $\sigma=v_c/\sqrt2\approx\mathbf{141\ km\,s^{-1}}$ .

2 · Building a Milky-Way-scale NFW halo (fully worked)

Take a halo of virial mass $M_{200}=10^{12}\,M_\odot$ and concentration $c=10$ , with $H_0=70\ \mathrm{km\,s^{-1}\,Mpc^{-1}}$ . Find $\rho_\mathrm{crit}$ , the virial radius $r_{200}$ , the scale radius $r_s$ , the characteristic density $\rho_s$ , and the circular speed at $r_s$ .

Step 1. Critical density from equation (4): $\rho_\mathrm{crit}=3H_0^2/8\pi G=9.2\times10^{-27}\ \mathrm{kg\,m^{-3}}=\mathbf{136\,M_\odot\,kpc^{-3}}$ . The cosmic yardstick against which the $\times200$ overdensity is set.

Step 2. Virial radius from $M_{200}=\tfrac43\pi r_{200}^3(200\rho_\mathrm{crit})$ : $r_{200}=\big[3M_{200}/(4\pi\cdot200\rho_\mathrm{crit})\big]^{1/3}\approx\mathbf{206\ kpc}$ . Then $r_s=r_{200}/c=\mathbf{20.6\ kpc}$ — the scale radius sits around the edge of the visible disk.

Step 3. Characteristic density from equation (5) with $m(10)=\ln11-10/11=1.489$ : $\rho_s=(200/3)\rho_\mathrm{crit}c^3/m(c)=66.7\times136\times1000/1.489\approx\mathbf{6.1\times10^{6}\,M_\odot\,kpc^{-3}}$ , i.e. $\delta_c=\rho_s/\rho_\mathrm{crit}\approx\mathbf{4.5\times10^{4}}$ . The halo core is tens of thousands of times denser than the mean Universe.

Step 4. Mass inside $r_s$ : $M(<r_s)=M_{200}[\ln2-\tfrac12]/m(c)=10^{12}\times0.193/1.489\approx1.3\times10^{11}\,M_\odot$ . Circular speed there, $v_c(r_s)=\sqrt{G M(<r_s)/r_s}\approx\mathbf{164\ km\,s^{-1}}$ ; the curve peaks near $174\ \mathrm{km\,s^{-1}}$ at $r\approx2.16\,r_s$ , and the virial velocity is $V_{200}=\sqrt{GM_{200}/r_{200}}\approx144\ \mathrm{km\,s^{-1}}$ — all in the Milky-Way range.

3 · A cluster the Zwicky way (order of magnitude)

The Coma-like cluster has a line-of-sight velocity dispersion $\sigma\approx1000\ \mathrm{km\,s^{-1}}$ over a radius $R\approx1\ \mathrm{Mpc}$ . Estimate its gravitating mass.

Step 1. For a virialized system $2T+W=0$ gives $M\sim\sigma^2 R/G$ . With $G=4.301\times10^{-9}\ \mathrm{Mpc\,(km\,s^{-1})^2}\,M_\odot^{-1}$ : $M\approx(1000)^2(1)/(4.301\times10^{-9})\approx\mathbf{2\times10^{14}\,M_\odot}$ .

Step 2. The full virial estimator carries a geometric prefactor of order $3$ – $5$ , pushing the mass toward $\mathbf{10^{15}\,M_\odot}$ . Against the roughly $10^{13}\,M_\odot$ in stars and hot gas, this is the "missing mass" Zwicky flagged in 1933 — the first detection of dark matter, decades before rotation curves confirmed it in galaxies.

Physics of the portrait

The story of dark-matter halos is the story of a shape: for fifty years a list of anomalies, it became physics when simulations showed that gravity, given cold collisionless matter, builds the same density profile every time.

Origins · 1 minNavarro–Frenk–White profile and how this was discovered

Navarro–Frenk–White profile

the universal CDM halo density ρ_s/((r/r_s)(1+r/r_s)²) · cold-dark-matter structure formation

The idea that galaxies float in far more mass than they show is older than any theory of what that mass is. In 1933 Fritz Zwicky applied the virial theorem to the Coma cluster and found its galaxies moving too fast for the visible mass to bind them — he named the shortfall dunkle Materie, dark matter, and was largely ignored. The galactic case matured slowly, but in the 1970s Vera Rubin and Kent Ford turned flat rotation curves from curiosity into an inescapable pattern across spiral after spiral, reinforced by Albert Bosma's 21-cm maps reaching far beyond the starlight. By the 1980s the missing mass was undeniable; what it was and how it was arranged stayed open. The answer came from computation. As cold-dark-matter $N$ -body simulations grew through the 1980s and early 1990s, Julio Navarro, Carlos Frenk and Simon White noticed that relaxed halos — regardless of mass or cosmology — settled into one two-parameter shape with an $r^{-1}$ inner cusp and an $r^{-3}$ outer fall, published in 1996–97 as the profile that now bears their initials. It was not derived from first principles but discovered in the machine, and it remains the organizing fact of halo physics. Jaan Einasto's gently curving profile (1965) fits the highest-resolution halos slightly better, but NFW is still the reference against which every rotation curve, lensing map, and cluster is measured.

History

1933

— Zwicky infers dark matter from the Coma cluster velocity dispersion

1970

— Rubin and Ford measure flat rotation curves in Andromeda

1978

— Bosma establishes flat 21-cm rotation curves as generic in spirals

1985

— cold-dark-matter N-body simulations reproduce hierarchical clustering

1996

— Navarro, Frenk and White propose the universal NFW density profile

Going further

Texts. Binney & Tremaine, Galactic Dynamics, 2nd ed., §2.2 (spherical systems), §4.3 (Jeans equations), §4.8 (violent relaxation) · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 7 (gravitational collapse), Ch. 11 (halo structure) · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §24.3 (galaxy masses and dark matter) · Dodelson, Modern Cosmology, 2nd ed., Ch. 12 (nonlinear structure and halos) · Ryden, Introduction to Cosmology, 2nd ed., Ch. 6–8. The equilibrium framework — the Collisionless Boltzmann Equation, the Jeans equations, and the Virial theorem — is applied throughout. The key relations used are: (1) $v_c^2=GM/r$ and the flat-curve $\Rightarrow\rho\propto r^{-2}$ argument; (2) the NFW integral $\int_0^y x(1+x)^{-2}\mathrm{d}x=\ln(1+y)-y/(1+y)$ and $\rho_s=\tfrac{200}{3}\rho_\mathrm{crit}c^3/m(c)$ ; (3) $v_c=\sqrt2\,\sigma$ , $\rho_\mathrm{iso}=\sigma^2/(2\pi G r^2)$ ; (4) $\rho_\mathrm{crit}=3H^2/8\pi G$ and the $\Delta=200$ overdensity definition. The worked-example numbers use $G=4.301\times10^{-6}\ \mathrm{kpc\,(km/s)^2}M_\odot^{-1}$ and $H_0=70\ \mathrm{km\,s^{-1}\,Mpc^{-1}}$ , giving $M(<20\,\mathrm{kpc})=1.9\times10^{11}M_\odot$ , $r_{200}=206\,\mathrm{kpc}$ , $\rho_s=6.1\times10^{6}M_\odot\,\mathrm{kpc}^{-3}$ , $v_c(r_s)=164\,\mathrm{km/s}$ , $V_{200}=144\,\mathrm{km/s}$ , and cluster $M\sim2\times10^{14}M_\odot$ .

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

Ask the tutor: Ask me to fit an NFW profile to a real galaxy's rotation curve, or to show why a cored dark-matter halo would bend the inner curve differently from a cusp.

Why is dark matter needed at all — could the rotation curves be modified gravity instead?

They could, and that is a serious minority program: Modified Newtonian Dynamics (MOND) reproduces flat rotation curves and the baryonic Tully–Fisher relation with a single acceleration scale and no dark matter. But halos are favored because the same invisible mass, arranged the same way, simultaneously explains cluster dynamics, X-ray gas, lensing (including the mass–gas offset in the Bullet Cluster), the CMB acoustic peaks, and structure growth — which modified gravity struggles to fit all at once. Dark matter's strength is not any single curve but the consistency of one hypothesis across a dozen independent scales.

What actually is the dark matter that makes up a halo?

Dynamically, a halo only cares that its matter is gravitating, cold (slow-moving when structure forms), and nearly collisionless — enough to fix the NFW shape. The particle identity is unknown. Leading candidates are weakly interacting massive particles (WIMPs), axions or other ultralight bosons, and sterile neutrinos, with primordial black holes a constrained possibility. Essentially ruled out is ordinary baryonic matter hiding as gas, dust, or faint stars: Big Bang nucleosynthesis and the CMB both cap the baryon density far below the total matter density. So the halo is real and its gravity mapped in detail, while the microphysics of its constituents remains a central open problem in physics.

Why is the NFW profile "universal" — where does it come from?

No one has derived it cleanly from first principles; it was found empirically in simulations and is best understood as an attractor of dissipationless collapse. As cold matter merges hierarchically, violent relaxation — the scattering of orbits off the rapidly changing potential of each merger — erases memory of the initial conditions and drives the system toward a preferred phase-space arrangement. The result is insensitive to halo mass, formation history, and even the initial power spectrum, which is why one two-parameter form fits dwarf galaxies to clusters, with mild residual deviations captured by the extra shape parameter of the Einasto profile.

What is the concentration parameter really telling me?

Concentration $c=r_{200}/r_s$ measures how centrally packed a halo is, and it is essentially a clock. A halo's central density reflects the density of the Universe when its core assembled: cores that formed early, when the cosmic mean density was high, are dense and therefore concentrated. Because small halos collapse first, low-mass halos are on average more concentrated than massive ones, and concentration at fixed mass falls toward the present as low-density material accretes and pushes $r_{200}$ outward faster than $r_s$ . This is why the mean concentration–mass–redshift relation, calibrated from simulations, both predicts a halo's full profile from its mass and tests the cosmological model.

Does the baryonic galaxy change the dark-matter halo, or are they independent?

They interact, and ignoring it is a known pitfall. As baryons cool and sink to form the disk, their added gravity pulls dark matter inward — adiabatic contraction — steepening the inner halo above the dark-matter-only prediction. Feedback runs the other way: repeated bursts of supernova- or AGN-driven outflow rapidly change the central potential, and this heating can transfer energy to the dark matter and flatten a cusp into a core, especially in dwarfs. The inner profile of a real halo is thus negotiated between collisionless dark matter and messy baryonic astrophysics — which is why the core–cusp problem is hard to adjudicate and why hydrodynamic, not gravity-only, simulations are needed for dwarf rotation curves.

How do we weigh a halo we cannot see — what are the independent methods?

Every method is gravitational, and their agreement is the point. Rotation curves and stellar kinematics use orbiting tracers through $v_c^2=GM(<r)/r$ and the Jeans equations. Cluster velocity dispersions use the virial theorem, $M\sim\sigma^2 R/G$ . X-ray observations assume the hot gas is in hydrostatic equilibrium and invert its temperature and density for the binding mass. Gravitational lensing measures the projected mass directly from light-bending, with no equilibrium assumption, and the Sunyaev–Zel'dovich effect weighs the gas via its CMB imprint. That these independent probes, with different systematics, return consistent masses is what makes the halo a measured object, not a modelling convenience.