Galactic Dark Matter Halos and Density Profiles

A dark matter halo is the quasi-equilibrium, self-gravitating cloud of collisionless matter inside which a luminous galaxy sits like a small bright ember in a much larger, invisible bonfire. Because dark matter feels only gravity, its distribution is shaped by gravitational collapse and violent relaxation rather than by radiative cooling — it cannot shed energy and shrink the way gas does, so it settles into an extended, velocity-supported structure whose density falls off slowly with radius. This makes halos both simpler in principle than galaxies and astonishingly regular: dissipationless simulations find that halos spanning fifteen orders of magnitude in mass — from dwarf-galaxy hosts to galaxy clusters — collapse to nearly the same density shape once rescaled by a characteristic radius and density. Three profiles organize everything that follows. The singular isothermal sphere, $\rho\propto r^{-2}$ , is the historical starting point: it reproduces the flat 1970s rotation curves exactly, but its mass diverges and its centre is singular. The Navarro–Frenk–White (NFW) profile is the modern standard from cold-dark-matter simulations — a two-parameter family whose slope steepens continuously from $-1$ at the centre to $-3$ in the outskirts. And the Einasto profile, with a slope curving smoothly as a power of radius, fits the highest-resolution simulations even better and cures NFW's central divergence. These are not idle fitting formulae — they encode the cosmological assembly history of structure, and the mismatch between their cuspy predictions and the cored centres many dwarf galaxies show (the core–cusp problem) is one of the sharpest tests we have of what dark matter is.

Key takeaways

A dark matter halo is the extended, gravitationally dominant mass reservoir in which a galaxy forms; its radial density $\rho(r)$ is captured by a small family of parametric profiles.
The singular isothermal sphere $\rho=\sigma_v^2/2\pi G r^2$ gives a perfectly flat rotation curve $v_c=\sqrt2\,\sigma_v$ , but has infinite mass and a central singularity.
Cold-dark-matter simulations produce the near-universal NFW profile — inner $r^{-1}$ cusp, outer $r^{-3}$ falloff — set by a scale radius $r_s$ and concentration $c=r_{200}/r_s$ .
The Einasto profile, with slope $2(r/r_{-2})^\alpha$ , fits high-resolution halos best and has a finite central density.
Halo mass and size are set by an overdensity threshold ( $M_{200}$ , $r_{200}$ at $200\rho_{\rm crit}$ ); concentration correlates with mass and formation epoch, acting as a fossil clock.
Key idea: every halo profile is characterized by its local logarithmic slope $\gamma(r)=-\mathrm{d}\ln\rho/\mathrm{d}\ln r$ — constant at $2$ for the isothermal sphere, running $1\to3$ for NFW, rising as $2(r/r_{-2})^\alpha$ for Einasto — and the rotation curve $v_c^2(r)=GM(<r)/r$ follows by integrating each $\rho(r)$ .

The core insight

Every halo profile is characterized by its local logarithmic slope $\gamma(r)=-\mathrm{d}\ln\rho/\mathrm{d}\ln r$ : constant ( $=2$ ) for the isothermal sphere, running from $1$ to $3$ for NFW, and rising as $2(r/r_{-2})^\alpha$ for Einasto. The rotation curve $v_c^2(r)=GM(<r)/r$ follows directly by integrating each $\rho(r)$ .

Where this lives today

Halo profiles are the working currency of near-field cosmology, and their sharpest edge is the core–cusp problem. Dissipationless CDM predicts a cuspy $\gamma\to1$ centre, yet the rotation curves of gas-rich dwarf and low-surface-brightness galaxies — where dark matter dominates — often prefer nearly flat cores. Surveys like THINGS and SPARC quantify the discrepancy across hundreds of decomposed curves, and the proposed resolutions each point somewhere profound: supernova feedback that flattens the cusp, self-interacting dark matter that thermalizes a core, or warm dark matter that suppresses small-scale structure. Gravitational lensing supplies a complementary, equilibrium-free probe — the analytic NFW surface density (Wright & Brainerd) lets stacked weak lensing confirm the NFW shape and mass–concentration relation across decades in mass. And Gaia's six-dimensional stellar phase space now constrains the Milky Way's own profile and its local dark-matter density — the very number direct-detection experiments need to turn a null result into a bound on the dark-matter particle.

The physical picture

Picture the dark matter as a swarm of collisionless particles that fell together as an overdense patch of the early universe decoupled from the expansion, turned around, and collapsed. During collapse the orbits phase-mix and the distribution relaxes — Lynden-Bell's violent relaxation — into a smooth, time-independent cloud held up by the random spread of particle velocities. Because it cannot radiate energy away, the matter never condenses into a thin disk the way cooling gas does; it stays puffy and extended, and the luminous galaxy is a small baryonic afterthought at its centre. The most useful diagnostic is not the density at any one radius but the rate at which it falls — the local slope $\gamma(r)=-\mathrm{d}\ln\rho/\mathrm{d}\ln r$ , a dimensionless "steepness dial" that reads off each profile's character independent of normalization, and exactly what a rotation curve or lensing map tries to measure. The first figure plots $\gamma(r)$ for all three profiles; every one passes through $\gamma=2$ at its own scale radius, which is why that radius is special.

Figure

The logarithmic slope γ = −d ln ρ / d ln r for the three canonical halo profiles versus radius in scale-radius units. Isothermal is pinned at γ = 2; NFW climbs from 1 to 3 as γ = 1 + 2x/(1+x); Einasto rises smoothly as γ = 2x^α with α = 0.17. All three cross γ = 2 at the scale radius x = 1

A flat rotation curve is the fingerprint of an $r^{-2}$ interior: if $\rho\propto r^{-2}$ then $M(<r)\propto r$ and $v_c^2=GM/r$ is independent of radius. Real halos are not exactly isothermal — shallower ( $\gamma\to1$ ) at the centre and steeper ( $\gamma\to3$ ) outside — so their curves rise, round over to a broad maximum, and decline slowly beyond it. The concentration $c$ sets where that maximum sits; halos that assembled early, when the cosmos was denser, froze in that higher density and are more concentrated, so $c$ is a fossil record of assembly time. The second figure contrasts the flat isothermal curve with a Milky-Way-scale NFW halo.

Figure

Rotation curves for a Milky-Way-scale halo. The isothermal sphere gives a strictly flat v_c = √2 σᵥ = 220 km/s with σᵥ = 156 km/s. An NFW halo of M₂₀₀ = 10¹² M☉, c = 12 (so rₛ = 17 kpc, r₂₀₀ = 206 kpc) rises to a broad maximum near r = 37 kpc then declines; the dark-matter part at the Sun's r = 8 kpc is about 145 km/s

The core equations

The starting point is the singular isothermal sphere. Treat the halo as a self-gravitating "gas" whose pressure comes from a constant velocity dispersion $\sigma_v$ rather than temperature; hydrostatic balance then forces a pure inverse-square density law — the simplest profile giving a flat rotation curve, and the reference against which the simulation-based profiles are read.

\rho_{\rm iso}(r)=\frac{\sigma_v^{2}}{2\pi G\,r^{2}} \qquad (1)

where: $\sigma_v$ — one-dimensional velocity dispersion, constant with radius; $G$ — gravitational constant; $r$ — spherical radius. A hotter halo (larger $\sigma_v$ ) is denser at every radius, and the inverse-square falloff is exactly what makes the rotation curve flat. The modern standard is the Navarro–Frenk–White profile, the shape dissipationless cold-dark-matter halos relax to across all masses. It is a two-parameter form — a characteristic density $\rho_s$ and scale radius $r_s$ — with a mild $r^{-1}$ inner cusp and a steep $r^{-3}$ outer decline, the transition centred on $r_s$ .

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

where: $\rho_s$ — characteristic density; $r_s$ — scale radius, where the slope equals $-2$ ; $r$ — radius. With $x=r/r_s$ : for $x\ll1$ , $\rho\to\rho_s/x$ (cusp); for $x\gg1$ , $\rho\to\rho_s/x^3$ (outer halo). Higher-resolution simulations showed the true slope curves smoothly rather than following NFW's fixed rational function — captured by the three-parameter Einasto profile, whose slope is itself a power law in radius. Crucially, its central density is finite.

\ln\!\frac{\rho_{\rm Ein}(r)}{\rho_{-2}}=-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^{\alpha}-1\right] \qquad (3)

where: $\alpha$ — shape parameter, $\approx0.16$ – $0.20$ for CDM halos; $r_{-2}$ — radius where the slope equals $-2$ ; $\rho_{-2}$ — density there. Smaller $\alpha$ curves the slope more slowly, giving a more NFW-like halo. The quantity unifying all three is the local logarithmic slope, the "steepness dial" of the first figure — what distinguishes a cusp from a core.

\gamma(r)\equiv-\frac{\mathrm{d}\ln\rho}{\mathrm{d}\ln r} \qquad (4)

where: $\gamma$ — local slope, the signature of each profile: $\gamma=2$ (isothermal), $\gamma=1+2x/(1+x)$ running $1\to3$ (NFW), $\gamma=2(r/r_{-2})^\alpha$ (Einasto). Integrating the NFW density in spherical shells gives a closed-form enclosed mass — the object fed into the rotation-curve and lensing formulae. It grows only logarithmically at large radius, so a finite halo mass needs a boundary.

M_{\rm 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 (5)

where: $M(<r)$ — mass within $r$ . The bracketed $m(x)=\ln(1+x)-x/(1+x)$ rises from $0$ (as $x^2/2$ ) and diverges logarithmically; it recurs throughout NFW bookkeeping. The dynamical payoff of any $\rho(r)$ is the circular speed of a test particle — what a rotation curve measures directly.

v_c^{2}(r)=\frac{GM(<r)}{r} \qquad (6)

where: $v_c$ — circular speed. With equation (1) it is constant; with equation (5) it rises, peaks near $r=2.16\,r_s$ , and declines. Because both masses diverge, a halo's mass and size are defined by an overdensity threshold: the radius $r_{200}$ inside which the mean density is $200\rho_{\rm crit}$ , and the mass $M_{200}$ it encloses.

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

where: $M_{200}$ — virial mass; $r_{200}$ — virial radius; $\rho_{\rm crit}$ — critical density; $H$ — Hubble parameter. The factor $200$ is a spherical-collapse convention; always state it, since $M_{200c}$ (vs. critical) and $M_{200m}$ (vs. mean) differ. Finally, the concentration ties the two NFW parameters together: requiring the mass inside $r_{200}$ to equal $M_{200}$ fixes $\rho_s$ in terms of $c$ alone.

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

where: $c$ — concentration, typically $5$ – $20$ (larger for earlier-forming halos). A halo is fully specified by two numbers, e.g. $(M_{200},c)$ ; everything else follows.

Derivation AFlat rotation curve from the singular isothermal sphere

The goal is to show that the inverse-square density law of equation (1) produces a strictly constant circular velocity — the property that, observed in real galaxies during the 1970s, forced astronomers to postulate extended dark halos. We derive the profile from hydrostatic equilibrium of an isothermal "gas", integrate it for the enclosed mass, then form the circular speed and watch every trace of radius cancel, leaving $v_c=\sqrt2\,\sigma_v$ .

Step 1. Model the halo as a fluid with "pressure" $P=\rho\sigma_v^2$ ( $\sigma_v$ constant) in hydrostatic balance against its own gravity: $\mathrm{d}P/\mathrm{d}r=-\rho\,GM(<r)/r^2$ , paired with $\mathrm{d}M/\mathrm{d}r=4\pi r^2\rho$ . Why: a self-gravitating isothermal system obeys the same equation as an isothermal gas — velocity dispersion plays the role of temperature, and the isotropic Jeans equation is formally this pressure balance.

Step 2. Try the power-law ansatz $\rho=A r^{-2}$ . Then $M(<r)=\int_0^r 4\pi r'^2(A r'^{-2})\,\mathrm{d}r'=4\pi A r$ , and the balance equation gives $-2\sigma_v^2 A r^{-3}=-4\pi G A^2 r^{-3}$ . Why: the inverse-square form makes both sides scale as $r^{-3}$ , so the equation collapses to an algebraic condition on $A$ with no leftover radial dependence.

Step 3. Cancel the common $r^{-3}$ and solve $2\sigma_v^2 A=4\pi G A^2$ , giving $A=\sigma_v^2/2\pi G$ , hence $\rho_{\rm iso}=\sigma_v^2/2\pi G r^2$ . The enclosed mass is then $M(<r)=4\pi A r=2\sigma_v^2 r/G$ , strictly linear in radius. Why: a mass rising in exact proportion to radius is precisely the condition that keeps $GM/r$ fixed.

Step 4. Form the circular speed: $v_c^2=GM(<r)/r=(G/r)(2\sigma_v^2 r/G)=2\sigma_v^2$ , so $v_c=\sqrt2\,\sigma_v$ , independent of $r$ ; the Milky Way's $220$ km/s implies $\sigma_v\approx156$ km/s. Why: the cancellation of $r$ is the content of a flat rotation curve, and inverting it reads a measured rotation speed directly as a halo dispersion. *Verified: Binney & Tremaine §4.3.3; $M(<r)=2\sigma_v^2 r/G$ and $v_c=\sqrt2\,\sigma_v$ exact, $220/\sqrt2=155.6$ km/s.*

Go deeper · Derivation B · 1 minThe NFW logarithmic slope runs from one to three

Derivation BThe NFW logarithmic slope runs from one to three

The goal is to compute the NFW local slope $\gamma(r)$ explicitly and show it climbs continuously from $1$ at the centre to $3$ in the outskirts, passing through $2$ exactly at the scale radius. This one calculation explains why NFW is "universal but not a power law", why $r_s$ is the natural pivot, and why the deep centre is a mild integrable cusp rather than a flat core — the crux of the core–cusp debate.

Step 1. Write $x=r/r_s$ and take the log of equation (2): $\ln\rho=\ln\rho_s-\ln x-2\ln(1+x)$ . Since $x$ and $r$ differ only by the constant $r_s$ , $\mathrm{d}\ln x=\mathrm{d}\ln r$ . Why: taking the log turns a product of powers into a sum, so the slope becomes a sum of simple terms — the cleanest route to $\mathrm{d}\ln\rho/\mathrm{d}\ln r$ .

Step 2. Differentiate term by term: $-\ln x$ contributes $-1$ , and $-2\ln(1+x)$ contributes $-2x/(1+x)$ via $\mathrm{d}\ln(1+x)/\mathrm{d}\ln x=x/(1+x)$ . Hence $\gamma=-\mathrm{d}\ln\rho/\mathrm{d}\ln r=1+2x/(1+x)$ . Why: the two terms carry the two limiting slopes — the $-1$ is the inner cusp, the $-2x/(1+x)$ is the extra steepening that switches on only once $x\gtrsim1$ .

Step 3. Read off the limits. As $x\to0$ , $\gamma\to1$ (the divergent-but-integrable cusp); as $x\to\infty$ , $\gamma\to3$ ; at $x=1$ , $\gamma=2$ , the isothermal value — which is why $r_s$ is the "slope $=-2$ " radius and why all three profiles meet there in the first figure. Why: a slope of $1$ inside and $3$ outside is exactly what makes a rotation curve rise, peak near $2.16\,r_s$ , and slowly fall — the shape simulations robustly predict. *Verified: matches Navarro, Frenk & White (1996, 1997); $\gamma(1)=2$ by construction, $\gamma\to1$ and $\gamma\to3$ at the two extremes.*

Go deeper · Derivation C · 1 minConcentration fixes the characteristic density

Derivation CConcentration fixes the characteristic density

The goal is to derive the NFW normalization $\rho_s$ in terms of the concentration $c$ by demanding self-consistency: the mass the profile encloses inside $r_{200}$ must equal $M_{200}$ , the mass the overdensity definition assigns there. This turns a two-parameter fit into a physical statement — because $c$ correlates with formation epoch, the resulting $\rho_s\propto c^3/\ln c$ quantifies how much denser early-forming halos are at their centres.

Step 1. Evaluate equation (5) at $r=r_{200}$ (where $r_{200}/r_s=c$ ) and set it equal to equation (7): $4\pi\rho_s r_s^3\big[\ln(1+c)-c/(1+c)\big]=\tfrac43\pi r_{200}^3(200\rho_{\rm crit})$ . Why: both sides are "the mass inside $r_{200}$ " — one from the density profile, one from the virial boundary — and equating them is the only consistency condition available, enough to fix the normalization.

Step 2. Substitute $r_{200}=c\,r_s$ so $r_{200}^3=c^3 r_s^3$ , and abbreviate $m(c)=\ln(1+c)-c/(1+c)$ : $4\pi\rho_s r_s^3\,m(c)=\tfrac43\pi c^3 r_s^3(200\rho_{\rm crit})$ , with $4\pi r_s^3$ on both sides. Why: expressing $r_{200}$ through $c$ and $r_s$ makes $r_s$ cancel, leaving a relation between $\rho_s/\rho_{\rm crit}$ and the pure number $c$ — the profile shape is scale-free.

Step 3. Cancel $4\pi r_s^3$ and solve: $\rho_s=\dfrac{200}{3}\rho_{\rm crit}\dfrac{c^3}{\ln(1+c)-c/(1+c)}$ . Since the denominator grows only logarithmically while $c^3$ grows as a cube, $\rho_s$ climbs steeply with concentration — a $c=20$ halo is far denser at its scale radius than a $c=5$ halo of the same mass. Why: this is the fossil-clock statement — $\rho_s\propto c^3/\ln c$ says an early-forming, high-concentration halo carries the higher background density of its birth epoch frozen into its core. *Verified: standard NFW normalization (Navarro, Frenk & White 1997); for $c=12$ , $m(c)=1.642$ and $\rho_s\approx7.0\times10^4\,\rho_{\rm crit}$ .*

Limiting cases

Each profile is valid in a regime, and the way their small- and large-radius behaviour differs is where the physics lives:

Case	Result	Meaning
Singular isothermal, all $r$	$\rho=\dfrac{\sigma_v^2}{2\pi G\,r^{2}},\quad v_c=\sqrt2\,\sigma_v$	Perfectly flat rotation curve; but $M(<r)=2\sigma_v^2 r/G$ diverges linearly, so the profile must be truncated at some outer radius.
NFW inner cusp, $r\ll r_s$	$\rho\to\dfrac{\rho_s r_s}{r},\quad\gamma\to1$	Density diverges but the enclosed mass stays finite ( $M\propto r^2$ ); this mild cusp is the CDM prediction at the heart of the core–cusp problem.
NFW outer halo, $r\gg r_s$	$\rho\to\dfrac{\rho_s r_s^{3}}{r^{3}},\quad\gamma\to3$	Mass diverges only logarithmically, so a finite halo mass requires an overdensity boundary such as $r_{200}$ .
Einasto centre, $r\to0$	$\rho\to\rho_{-2}\,e^{2/\alpha}$ (finite), $\gamma\to0$	Finite central density and zero slope — a genuine core, unlike NFW; for $\alpha=0.17$ the centre is $e^{2/\alpha}\approx1.3\times10^5$ times $\rho_{-2}$ .
Rotation curve, any $\rho$	$v_c^2(r)=GM(<r)/r$	A running integral of the density; its shape reports enclosed mass, not local density.

Common trapthe misconception that the scale radius is the halo's edge

It is tempting to read $r_s$ as "the size of the halo", but $r_s$ is only the interior pivot where the NFW slope passes through $-2$ ; the actual size is the virial radius $r_{200}=c\,r_s$ , five to twenty times larger. Two traps ride along. First, quoting a halo mass without the overdensity convention is ambiguous: $M_{200c}$ (relative to critical), $M_{200m}$ (relative to mean), and the virial $M_{\rm vir}$ differ by tens of percent for the same halo. Second, an NFW "cusp" is not an infinite central mass — the density diverges as $r^{-1}$ , but $\int r^2\rho\,\mathrm{d}r\propto\int r\,\mathrm{d}r$ converges, so the enclosed mass goes smoothly to zero. Cuspy density, finite mass: the core–cusp debate is about the slope, never a mass singularity.

Worked examples

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

Take a halo of virial mass $M_{200}=1.0\times10^{12}\,M_\odot$ and concentration $c=12$ , with $H_0=70$ km/s/Mpc so that $\rho_{\rm crit}=1.36\times10^{2}\,M_\odot\,\mathrm{kpc}^{-3}$ . Find $r_{200}$ , $r_s$ , and $\rho_s$ , then the dark-matter circular speed at the Sun's radius $r=8$ kpc.

Step 1. Virial radius from equation (7): $r_{200}^3=\dfrac{3M_{200}}{4\pi(200\rho_{\rm crit})}=\dfrac{3\times10^{12}}{3.42\times10^5}=8.78\times10^{6}\,\mathrm{kpc}^3$ , so $r_{200}=\mathbf{206\ kpc}$ . The scale radius follows from the concentration: $r_s=r_{200}/c=206/12=\mathbf{17.2\ kpc}$ .

Step 2. Characteristic density from equation (8). The mass function $m(c)=\ln13-12/13=2.565-0.923=\mathbf{1.642}$ , so $\rho_s=\dfrac{200}{3}\rho_{\rm crit}\dfrac{c^3}{m(c)}=66.7\cdot(1.36\times10^2)\cdot\dfrac{1728}{1.642}=\mathbf{9.5\times10^{6}\ M_\odot\,kpc^{-3}}$ — about $7\times10^4$ times critical, so even an ordinary halo core is enormously overdense.

Step 3. Enclosed dark mass at the Sun. With $x=8/17.2=0.465$ , equation (5) gives $M(<8)=4\pi\rho_s r_s^3\,m(x)$ with $m(0.465)=\ln1.465-0.465/1.465=0.0645$ . The prefactor $4\pi\rho_s r_s^3=6.1\times10^{11}\,M_\odot$ , so $M(<8)=6.1\times10^{11}\times0.0645=\mathbf{3.9\times10^{10}\ M_\odot}$ .

Step 4. Circular speed from equation (6), with $G=4.30\times10^{-6}\,\mathrm{kpc\,(km/s)^2}\,M_\odot^{-1}$ : $v_c^2=GM(<8)/8=4.30\times10^{-6}\cdot3.9\times10^{10}/8=2.1\times10^{4}\,(\mathrm{km/s})^2$ , so $v_c\approx\mathbf{145\ km/s}$ . This is the dark-matter-only contribution at the solar circle; the observed $\approx220$ km/s is completed by the disk and bulge — disentangling them is the rotation-curve decomposition problem. Two numbers $(M_{200},c)=(10^{12}M_\odot,12)$ fix everything: $r_{200}=206$ kpc, $r_s=17.2$ kpc, $\rho_s=9.5\times10^6\,M_\odot\,\mathrm{kpc}^{-3}$ , DM $v_c(8\,\mathrm{kpc})\approx145$ km/s.

2 · The isothermal sphere and why it cannot be the whole story

The Milky Way's rotation curve is roughly flat at $v_c\approx220$ km/s. Model it as a singular isothermal sphere: find the velocity dispersion, the mass enclosed within 8 kpc, and explain why the model must break down.

Step 1. From the flat-curve result $v_c=\sqrt2\,\sigma_v$ , the dispersion is $\sigma_v=v_c/\sqrt2=220/1.414=\mathbf{156\ km/s}$ . A single number sets the entire profile — an isothermal halo is a one-parameter object.

Step 2. Enclosed mass from $M(<r)=2\sigma_v^2 r/G$ . At $r=8$ kpc: $M(<8)=2(156)^2(8)/(4.30\times10^{-6})=\mathbf{9.0\times10^{10}\ M_\odot}$ . Because $M\propto r$ , doubling the radius doubles the mass — no convergence.

Step 3. That linear growth is the fatal flaw: $M(<r)=2\sigma_v^2 r/G\to\infty$ , so a real halo must be truncated by hand. The central singularity $\rho\to\infty$ as $r\to0$ is the second flaw; the cored "pseudo-isothermal" variant $\rho=\rho_0/(1+r^2/r_c^2)$ patches it with a core radius $r_c$ .

Step 4. As a sanity check, the implied local density $\rho_{\rm iso}(8)=\sigma_v^2/(2\pi G r^2)\approx0.53$ GeV/cm³ overshoots the favoured $0.3$ – $0.4$ GeV/cm³ — because the model attributes all of the $220$ km/s to a dark $r^{-2}$ halo, ignoring the disk. $\sigma_v=156$ km/s, $M(<8\,\mathrm{kpc})=9.0\times10^{10}\,M_\odot$ ; the mass diverges and the centre is singular, so the model needs truncation and a core.

Physics of the portrait

The dark matter halo is where cosmology and galactic dynamics meet: inferred from missing mass, measured through flat rotation curves, and finally given a universal shape by supercomputers. No single person owns it — it is the collective achievement of the observers who found the mass and the simulators who explained its distribution.

Origins · 1 mindark matter halo and how this was discovered

dark matter halo

the collisionless mass reservoir with a near-universal density profile ρ(r) · cosmological structure formation and galactic dynamics

The idea that galaxies are embedded in far more mass than they show began with Fritz Zwicky, who in 1933 applied the virial theorem to the Coma cluster and found galaxies moving so fast that unseen "dunkle Materie" must bind them — ignored for decades. The evidence became inescapable when Vera Rubin and Kent Ford, measuring optical rotation curves of Andromeda and dozens of spirals through the 1970s, found circular speeds staying flat far beyond the visible disk instead of falling off in Keplerian fashion; radio astronomers led by Albert Bosma traced the same flatness in neutral-hydrogen curves to even larger radii. A flat curve demands $M(<r)\propto r$ , hence an extended $r^{-2}$ dark halo. The theoretical shape came from the computer: building on the cold-dark-matter program of the 1980s, Julio Navarro, Carlos Frenk and Simon White ran dissipationless N-body simulations in the mid-1990s and found that halos of every mass relax to one two-parameter form — the NFW profile — distinguished only by concentration. Higher resolution later revealed a gently curving slope, reviving a 1965 formula of the Estonian astronomer Jaan Einasto. What began as Zwicky's anomaly is now a precision object whose profile tests the microphysics of dark matter itself.

History

1933

— Zwicky infers unseen mass binding the Coma cluster from galaxy velocities

1970

— Rubin and Ford measure the flat optical rotation curve of Andromeda

1978

— Bosma establishes flat neutral-hydrogen rotation curves as generic in spirals

1996

— Navarro, Frenk and White identify the universal two-parameter halo profile

2010

— high-resolution simulations favour the curved-slope Einasto form

Going further

Texts. Binney & Tremaine, Galactic Dynamics, 2nd ed., §2.2 (spherical systems and the isothermal sphere) and §4.3 · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 24 (galactic dynamics and dark matter) · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 7 (halo structure and the NFW profile) · Dodelson, Modern Cosmology, 2nd ed., Ch. 12 (structure formation) · Ryden, Introduction to Cosmology, 2nd ed., Ch. 8.

Provenance. A verifier should confirm: (1) $\rho_{\rm iso}=\sigma_v^2/2\pi G r^2$ gives $v_c=\sqrt2\,\sigma_v$ , so $220/\sqrt2=156$ km/s (Binney & Tremaine §4.3); (2) the NFW mass integral $\int x/(1+x)^2\,\mathrm{d}x=\ln(1+x)+1/(1+x)$ and slope $\gamma=1+2x/(1+x)$ with $\gamma(1)=2$ (Navarro, Frenk & White 1996, 1997); (3) the Einasto slope $\gamma=2(r/r_{-2})^\alpha$ and finite centre $\rho(0)=\rho_{-2}e^{2/\alpha}$ ; (4) $\rho_s=(200/3)\rho_{\rm crit}c^3/[\ln(1+c)-c/(1+c)]$ with $m(12)=1.642$ ; (5) the worked halo $M_{200}=10^{12}M_\odot$ , $c=12$ , $H_0=70$ giving $r_{200}=206$ kpc, $r_s=17.2$ kpc, $\rho_s=9.5\times10^6\,M_\odot\,\mathrm{kpc}^{-3}$ , DM $v_c(8\,\mathrm{kpc})\approx145$ km/s. Standard values ( $v_c\approx220$ km/s, local DM density $0.3$ – $0.4$ GeV/cm³, $\alpha\approx0.16$ – $0.20$ , $c\sim5$ – $20$ ) are textbook.

Prerequisites. Galactic Rotation Curves · Evidence for Dark Matter in Galaxies · Structure of the Milky Way. Leads to. Galactic Halo and Globular Cluster System · Neutral Hydrogen 21-cm Line and Galactic Kinematics · Elliptical Galaxies and the de Vaucouleurs Profile.

Ask the tutor: Ask me to build a halo from any $(M_{200},c)$ you pick and trace its rotation curve, or to walk through why dwarf galaxies turn the core–cusp problem into a test of dark-matter microphysics.

Is a dark matter halo really spherical, or is that just an idealization?

It is an idealization, and a decent one. Simulated CDM halos are actually triaxial — mildly flattened ellipsoids with axis ratios around $0.6$ – $0.8$ and orientations that twist with radius, because they accrete anisotropically along cosmic filaments. Spherical profiles like NFW and Einasto describe the angle-averaged density, which is what a rotation curve or a stacked lensing measurement mostly senses, and central baryons make the inner halo rounder anyway. So "spherical" is shorthand for the radial trend; when you need the shape itself — for stellar streams or satellite tides — you restore the triaxiality explicitly.

What is the difference between the scale radius rₛ and the virial radius r₂₀₀?

They answer different questions. The virial radius $r_{200}$ is the halo's size — the edge of the collapsed object, inside which the mean density is $200\rho_{\rm crit}$ . The scale radius $r_s$ is an interior landmark where the NFW slope passes through $-2$ , marking the transition from inner cusp to outer skirt and sitting near the rotation-curve peak. Their ratio is the concentration $c=r_{200}/r_s$ , typically $5$ – $20$ , so the halo extends well beyond $r_s$ . Confusing the two is the most common error in reading a halo: $r_s$ is a feature of the profile, $r_{200}$ is the boundary.

Why do simulations predict cusps but many galaxies show cores?

This is the core–cusp problem, and it is genuinely unresolved. Dissipationless CDM simulations robustly produce a central $\rho\propto r^{-1}$ cusp, yet the innermost rotation curves of many dwarf and low-surface-brightness galaxies are better fit by a constant-density core. Three camps compete: one says the simulations are incomplete, because bursts of star formation drive gas in and out and the fluctuating potential heats the dark matter and erodes the cusp; a second says dark matter is weakly self-interacting, so particle scatterings thermalize a core; a third says small-scale power is suppressed, as in warm dark matter. Distinguishing them from real galaxy data is a frontier of the field.

How can an NFW cusp have infinite central density but finite mass?

Because mass weights density by volume, and a thin shell's volume shrinks faster than the density grows. Near the centre $\rho\propto r^{-1}$ , but the shell mass is $\mathrm{d}M=4\pi r^2\rho\,\mathrm{d}r\propto r\,\mathrm{d}r$ , which integrates to $M(<r)\propto r^2\to0$ smoothly. So the "cusp" is about density, not mass — no central point mass, no divergent gravity. The divergence would only bite if the inner slope were steeper than $-3$ ; NFW's mild $-1$ is comfortably integrable, which is exactly why the profile is physically usable despite its formally infinite central density.

Why 200 times the critical density — is that number physical?

Partly. It comes from the spherical-collapse model, in which a just-virialized overdensity has a mean interior density of about $18\pi^2\approx178$ times the background of an Einstein–de Sitter universe — conventionally rounded to $200$ . But it is a convention, and several coexist: $200$ relative to critical ( $M_{200c}$ ), $200$ relative to mean matter ( $M_{200m}$ ), or a redshift-dependent virial $\Delta_c(z)$ (about $103$ today). These differ by tens of percent for the same halo, and concentration is quoted with respect to whichever radius is chosen. The lesson is operational: a halo mass is meaningless until you state the overdensity definition attached to it.

Do baryons change the dark-matter profile, or just add to it?

They change it, because the coupling runs both ways. When baryons cool and sink to the centre they deepen the potential and the dark matter contracts inward — adiabatic contraction steepens the inner halo. Pulling the other way, bursty feedback from supernovae and active galactic nuclei moves large gas masses on short timescales, transferring energy to the dark matter and puffing the cusp into a core. Which wins depends on galaxy mass and star-formation history: contraction in massive quiescent systems, feedback-driven flattening in gas-rich dwarfs. So a real galaxy's dark-matter profile is a baryon-processed version of the pristine simulation prediction — which is why the cleanest tests use the most dark-matter-dominated dwarfs.