NFW and Density Profiles of Halos

Every galaxy sits at the bottom of a much larger, invisible gravitational well — a dark-matter halo whose mass outweighs the stars and gas by an order of magnitude and whose structure sets the rotation curve, the lensing signal, and the orbits of satellites. Dark matter is collisionless: it feels only gravity, so its equilibrium is not set by pressure balance but by the collisionless Boltzmann equation, and for realistic cosmological initial conditions that equation has no closed-form solution. One might therefore expect halos to be a zoo of shapes, each remembering the accidents of its own collapse. The great surprise of 1990s computational cosmology was the opposite. When Navarro, Frenk & White followed the gravitational collapse and merging of cold dark matter in large N-body simulations, they found that the spherically-averaged density profiles of the resulting halos were nearly universal. Objects born from wildly different fluctuations, and spanning from dwarf galaxies to the richest clusters — more than five decades in mass — all settle onto one two-parameter curve. That curve, the NFW profile, has an inner cusp where density rises as $r^{-1}$ , an outer skirt falling as $r^{-3}$ , and a single scale radius $r_s$ marking the transition where the logarithmic slope passes through $-2$ . Universality is what makes the profile powerful: a halo is pinned down by just its virial mass $M_\Delta$ and its concentration $c = r_\Delta/r_s$ . This article develops the NFW form from the collapse physics that produces it, derives its mass, potential, and rotation curve, and confronts its two live tensions — the cuspy centre that clashes with dwarf-galaxy rotation curves, and the Einasto profile that fits the simulations slightly better.

Key takeaways

Cosmological N-body simulations show dark-matter halos across five decades in mass converge on one nearly universal spherically-averaged form — the Navarro–Frenk–White (NFW) profile.
The profile $\rho(r)=\rho_s/[(r/r_s)(1+r/r_s)^2]$ has an inner cusp $\rho\propto r^{-1}$ , an outer decline $\rho\propto r^{-3}$ , and a scale radius $r_s$ where the logarithmic slope is exactly $-2$ .
A halo needs only two numbers: its virial mass $M_\Delta$ and its concentration $c=r_\Delta/r_s$ ; concentration tracks formation epoch, so low-mass halos that assembled early are more concentrated.
The enclosed mass grows only as $\ln(1+r/r_s)-(r/r_s)/(1+r/r_s)$ , so the total mass diverges logarithmically and the profile must be truncated at a virial radius.
The rotation curve peaks at a fixed $r_{\max}\approx 2.163\,r_s$ , giving $V_{\max}$ as a concentration-sensitive observable, and the potential is finite at the centre despite the density cusp.
The $r^{-1}$ cusp clashes with cored dwarf rotation curves (the core–cusp problem), and the Einasto profile with a gently varying slope fits high-resolution halos better.
Key idea: collisionless gravitational collapse drives halos toward a universal two-parameter profile with an $r^{-1}$ inner cusp and $r^{-3}$ outer fall-off; the scale radius $r_s$ and concentration $c$ capture the entire structure.

The core insight

Collisionless gravitational collapse drives dark-matter halos toward a universal two-parameter profile with an $r^{-1}$ inner cusp and an $r^{-3}$ outer fall-off; a single scale radius $r_s$ and a concentration $c=r_\Delta/r_s$ fix the whole structure.

Where this lives today

The NFW profile is the standard model of the dark sector on halo scales, and it is doing daily work across observational cosmology. Weak-gravitational-lensing surveys — the Dark Energy Survey, KiDS, HSC, and soon Euclid and the Rubin Observatory's LSST — stack the shear signal around clusters and galaxies and fit NFW (or Einasto) profiles to recover $M_{200}$ and $c$ purely from gravity, independent of any baryonic tracer; the resulting mass–concentration relation is a precision test of $\Lambda$ CDM. In galaxy dynamics, NFW is the dark-halo term fitted jointly with stellar and gas components to rotation curves (via $v_c^2=GM/r$ ) and to pressure-supported dwarfs (via the Jeans equations), and the tension it generates — many dwarf and low-surface-brightness galaxies prefer a cored inner profile over the $r^{-1}$ cusp — is one of the sharpest small-scale challenges to cold dark matter, driving both baryonic-feedback models (supernova-driven potential fluctuations that flatten cusps) and alternative microphysics such as self-interacting dark matter. Closer to home, the Milky Way's halo is modelled as NFW to predict satellite orbits, stellar-stream tracks, and the local dark-matter density that sets the expected signal in direct-detection experiments like LZ and XENONnT. And every large cosmological simulation — IllustrisTNG, EAGLE, the FLAMINGO and Millennium suites — is validated in part by checking that its halos reproduce the NFW/Einasto form and its concentration statistics.

The physical picture

Picture dark matter raining into a growing potential well through a hierarchy of mergers: small clumps collapse first, then assemble into larger ones. Material that turns around and collapses early does so from a denser background universe, and it ends up deep inside the halo as a tightly bound core; material accreting later falls in from a diluter cosmos and settles at large radius with lower density. The halo is therefore not a single power law but a graded structure — shallower than isothermal in the middle, steeper than isothermal outside — and the crossover happens at one characteristic radius. The first figure is that structure drawn in the only honest coordinates for a scale-free-looking object, log density against log radius, so the two asymptotic slopes and their pivot become straight lines.

Figure

Log-log NFW density profile ρ = ρ_s divided by x(1+x)² with x = r/r_s. The logarithmic slope dlnρ/dlnr = −1 − 2x/(1+x) runs from −1 (inner cusp) through exactly −2 at r = r_s to −3 (outer envelope), so the curve bends smoothly between two straight power laws

The scale radius $r_s$ is the pivot: inside it the halo is cusp-dominated and the rotation curve is still rising; outside it the halo is envelope-dominated and the enclosed mass creeps up only logarithmically. Because the slope at $r_s$ is exactly $-2$ — the slope of a singular isothermal sphere — one can read $r_s$ as the radius at which the halo momentarily "looks isothermal" before steepening. Concentration $c=r_\Delta/r_s$ is just the ratio of the outer truncation radius to this pivot: a big number means the truncation radius sits far out on the steep envelope, so most of the mass is packed deep inside. Two halos of identical virial mass but different assembly histories differ only in $c$ : the early-forming one has a small $r_s$ , high $c$ , and a deeper central potential; the late-forming one is fluffier. The second figure turns this into the observable that a dynamicist actually measures — the circular-velocity curve of a Milky-Way-scale halo.

Figure

Circular velocity of a Milky-Way-scale NFW halo with M₂₀₀ = 10¹² M_☉ and concentration c = 12, giving r_s = 17.6 kpc and r₂₀₀ = 212 kpc. The curve rises, peaks at V_max = 179 km/s at r_max = 2.163 r_s = 38 kpc, then declines — V₂₀₀ = 143 km/s at the virial radius

Two features of that curve encode the whole profile. The peak at $r_{\max}\approx 2.16\,r_s$ is a fixed multiple of the scale radius, so measuring where the rotation curve turns over locates $r_s$ and hence $c$ without ever touching the normalization. And the fact that $v_c$ stays large and only gently falling out to hundreds of kiloparsecs — never plunging like a Keplerian $r^{-1/2}$ — is the dynamical fingerprint of the extended $r^{-3}$ envelope that first announced dark matter through flat galactic rotation curves.

The core equations

The profile itself is the starting point. NFW write the density as a broken power law with one normalization and one length, and everything else — mass, potential, rotation curve, slope — follows by calculus.

\rho(r)=\frac{\rho_s}{\dfrac{r}{r_s}\left(1+\dfrac{r}{r_s}\right)^{2}} \qquad (1)

where: $\rho(r)$ — spherically-averaged mass density at radius $r$ ; $\rho_s$ — characteristic density normalization, four times the density at $r_s$ ; $r_s$ — scale radius, where the logarithmic slope is $-2$ ; $r$ — radius from the halo centre. Reading it: near the centre ( $r\ll r_s$ ) the bracket is $\approx 1$ and $\rho\propto r^{-1}$ ; far out ( $r\gg r_s$ ) the bracket dominates and $\rho\propto r^{-3}$ . Integrating the density over spherical shells gives the enclosed mass. The integral is elementary and produces the signature NFW combination of a logarithm and a rational term.

M(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 (2)

where: $M(r)$ — mass within radius $r$ . Reading it: the bracket, call it $g(x)$ with $x=r/r_s$ , grows like $x^2/2$ at small $x$ (so $M\propto r^2$ , and the cusp holds negligible mass) and like $\ln x$ at large $x$ (so $M$ diverges — the halo must be truncated). A halo is normally quoted by its virial mass and concentration rather than by $\rho_s$ and $r_s$ . The virial radius $r_\Delta$ encloses a mean overdensity $\Delta$ times a reference cosmic density, and the concentration is its ratio to the scale radius.

M_\Delta=\frac{4}{3}\pi r_\Delta^{3}\,\Delta\,\rho_{\rm ref},\qquad c\equiv\frac{r_\Delta}{r_s} \qquad (3)

where: $M_\Delta$ — mass within $r_\Delta$ ; $\Delta$ — overdensity threshold (commonly 200); $\rho_{\rm ref}$ — reference density, either the critical density $\rho_{\rm crit}$ or the mean matter density $\bar\rho_m$ ; $c$ — concentration. Reading it: choosing $\Delta$ and $\rho_{\rm ref}$ turns a mathematically infinite object into a finite halo of definite mass and size. Matching $M(r_\Delta)$ from (2) to the virial mass in (3) fixes the characteristic density purely in terms of the concentration — the single most useful bookkeeping relation for NFW halos.

\rho_s=\frac{\Delta}{3}\,\rho_{\rm ref}\,\frac{c^{3}}{\ln(1+c)-\dfrac{c}{1+c}} \qquad (4)

where: all symbols as above. Reading it: $\rho_s/\rho_{\rm ref}$ is a pure function of $c$ (often written $\delta_c$ ); a more concentrated halo has a much larger characteristic density — for $\Delta=200$ , $c=5$ gives $\delta_c\approx 8.7\times10^{3}$ while $c=15$ gives $\delta_c\approx 1.2\times10^{5}$ . Because the mass is analytic, so is the potential from the spherical Poisson equation. The remarkable feature is that it is finite at the centre even though the density diverges there.

\Phi(r)=-\frac{4\pi G\rho_s r_s^{3}}{r}\,\ln\!\left(1+\frac{r}{r_s}\right),\qquad \Phi(0)=-4\pi G\rho_s r_s^{2} \qquad (5)

where: $\Phi(r)$ — gravitational potential; $G$ — Newton's constant. Reading it: the $r^{-1}$ cusp is weak enough that $M(r)\propto r^2\to 0$ faster than $1/r$ diverges, so orbits near the centre feel a finite well, not an infinite spike. Setting $v_c^2=GM/r$ gives the rotation curve, whose shape is universal in $x=r/r_s$ ; only the amplitude depends on the halo's mass and concentration.

v_c^{2}(r)=\frac{GM(r)}{r}=\frac{4\pi G\rho_s r_s^{3}}{r}\left[\ln(1+x)-\frac{x}{1+x}\right],\qquad x=\frac{r}{r_s} \qquad (6)

where: $v_c(r)$ — circular speed at radius $r$ . Reading it: $v_c$ rises inside $r_s$ , peaks at $r_{\max}=2.163\,r_s$ , and declines slowly outside — the gentle, non-Keplerian fall-off that is the signature of dark halos. The logarithmic slope makes the "runs from $-1$ to $-3$ " statement precise and pins the meaning of $r_s$ .

\frac{\mathrm{d}\ln\rho}{\mathrm{d}\ln r}=-1-\frac{2x}{1+x},\qquad x=\frac{r}{r_s} \qquad (7)

where: $x$ — dimensionless radius. Reading it: at $x\to0$ the slope is $-1$ ; at $x=1$ it is exactly $-2$ ; at $x\to\infty$ it approaches $-3$ — a single smooth function with no free parameters. High-resolution simulations show a small but systematic improvement if the slope is allowed to vary as a power of radius rather than through the fixed NFW form; this is the Einasto profile, now the standard refinement.

\rho_{\rm E}(r)=\rho_{-2}\,\exp\!\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^{\!\alpha}-1\right]\right\} \qquad (8)

where: $\alpha$ — shape index (typically $0.15$ – $0.25$ , mildly mass-dependent); $r_{-2}$ — radius where the slope is $-2$ ; $\rho_{-2}$ — density there. Reading it: unlike NFW the Einasto density is finite at the centre (no cusp singularity) and its slope $-2(r/r_{-2})^\alpha$ varies continuously, which is what lets it track simulated halos to sub-percent accuracy.

Derivation AThe enclosed mass and why the halo must be truncated

The goal is to integrate the NFW density over spherical shells to obtain the enclosed mass $M(r)$ , and in doing so to expose the two behaviours that dominate NFW life: the cusp holds almost no mass, while the envelope holds so much that the total diverges. We non-dimensionalize with $x=r/r_s$ , reduce the density integral to a standard rational integral, and read off the small- $x$ and large- $x$ limits to see why a virial truncation is not optional but forced.

Step 1. Write the mass in shells, $M(r)=\int_0^r 4\pi r'^2\,\rho(r')\,\mathrm{d}r'$ , valid for any spherically symmetric distribution because each shell of radius $r'$ contributes $4\pi r'^2\rho\,\mathrm{d}r'$ . Why: spherical symmetry lets a three-dimensional integral collapse to one dimension, the only reason NFW admits closed forms at all.

Step 2. Substitute $\rho=\rho_s/[x'(1+x')^2]$ and $r'=r_s x'$ so $\mathrm{d}r'=r_s\,\mathrm{d}x'$ , giving $M=4\pi\rho_s r_s^3\int_0^{x}\dfrac{x'^2}{x'(1+x')^2}\,\mathrm{d}x'$ . The factor $r_s^3$ emerges because volume scales as length cubed. Why: pulling out $r_s^3$ isolates the pure shape integral, which depends only on the dimensionless upper limit $x$ .

Step 3. Cancel one power of $x'$ to reach $\int_0^x \dfrac{x'}{(1+x')^2}\,\mathrm{d}x'$ . Split the integrand as $\dfrac{x'}{(1+x')^2}=\dfrac{1}{1+x'}-\dfrac{1}{(1+x')^2}$ , which integrates term by term to $\ln(1+x)+\dfrac{1}{1+x}-1=\ln(1+x)-\dfrac{x}{1+x}$ . Why: partial fractions turn an unfamiliar integral into two textbook ones, and the constant of integration is fixed by $M(0)=0$ .

Step 4. Assemble $M(r)=4\pi\rho_s r_s^3\left[\ln(1+x)-\dfrac{x}{1+x}\right]$ . Expand for small $x$ : $\ln(1+x)-x/(1+x)\approx x^2/2$ , so $M\propto r^2$ and the cusp is mass-light. Expand for large $x$ : the bracket $\to\ln x$ , so $M\to\infty$ logarithmically. Why: the divergence is the whole reason a "halo mass" is meaningless without a radius — one must cut the integral at $r_\Delta$ and quote $M_\Delta$ . *The partial-fraction integral does indeed yield $\ln(1+x)-x/(1+x)$ ; the small- $x$ expansion gives $x^2/2$ , confirming $M\propto r^2$ , and the large- $x$ behaviour grows as $\ln x$ , confirming the logarithmic divergence.*

Go deeper · Derivation B · 1 minThe logarithmic slope and the meaning of the scale radius

Derivation BThe logarithmic slope and the meaning of the scale radius

The goal is to show that the single NFW form contains a slope that slides continuously from $-1$ at the centre to $-3$ far out, passing through exactly $-2$ at one special radius, and that this special radius is precisely $r_s$ . This is what justifies calling $r_s$ "the" scale radius and what lets NFW interpolate between a cusp and a steep envelope without any additional shape parameter. We take the logarithmic derivative of the density and evaluate its limits and its midpoint.

Step 1. Take logs of $\rho=\rho_s/[x(1+x)^2]$ to get $\ln\rho=\ln\rho_s-\ln x-2\ln(1+x)$ , converting a product into a sum. Why: the logarithmic slope $\mathrm{d}\ln\rho/\mathrm{d}\ln r$ is the natural diagnostic for near-power-law density fields, and logs make the derivative trivial.

Step 2. Differentiate with respect to $\ln x$ (equal to $\ln r$ up to the constant $\ln r_s$ ): $\dfrac{\mathrm{d}\ln\rho}{\mathrm{d}\ln x}=x\dfrac{\mathrm{d}\ln\rho}{\mathrm{d}x}=x\left(-\dfrac{1}{x}-\dfrac{2}{1+x}\right)$ . Why: multiplying by $x$ converts an ordinary derivative into a logarithmic one, giving a dimensionless slope directly comparable across radii.

Step 3. Simplify to $\dfrac{\mathrm{d}\ln\rho}{\mathrm{d}\ln r}=-1-\dfrac{2x}{1+x}$ . As $x\to0$ the second term vanishes and the slope $\to-1$ (the cusp); as $x\to\infty$ it $\to-1-2=-3$ (the envelope). Why: a single rational expression encodes both asymptotes, which is exactly the universality NFW discovered — no piecewise fitting needed.

Step 4. Evaluate at $x=1$ , i.e. $r=r_s$ : $-1-\dfrac{2}{2}=-2$ , the slope of a singular isothermal sphere. Why: this fixes the operational meaning of the scale radius — it is the pivot where the halo transitions from cusp-dominated to envelope-dominated, and where it momentarily mimics the flat-rotation-curve isothermal model. *The slope formula evaluates to $-1.02$ at $x=0.01$ , exactly $-2$ at $x=1$ , and $-2.98$ at $x=100$ , tracing the full $-1\to-2\to-3$ progression.*

Go deeper · Derivation C · 1 minConcentration, characteristic density and the virial condition

Derivation CConcentration, characteristic density and the virial condition

The goal is to eliminate the abstract normalization $\rho_s$ in favour of the observable concentration $c$ , by demanding that the mean density inside the virial radius equal $\Delta$ times the reference density. This turns the two mathematical parameters $(\rho_s,r_s)$ into the two physical ones $(M_\Delta,c)$ that simulations and observations actually report, and it makes explicit that $\rho_s$ is not a central density but a concentration-dependent bookkeeping constant. We equate the enclosed-mass formula at $r_\Delta$ to the definition of virial mass and solve.

Step 1. Evaluate the enclosed mass at the virial radius using $r_\Delta=c\,r_s$ : $M(r_\Delta)=4\pi\rho_s r_s^3\left[\ln(1+c)-\dfrac{c}{1+c}\right]$ . Why: the concentration is precisely the upper limit of the shape integral, so all the profile's mass content is now carried by the single function $g(c)=\ln(1+c)-c/(1+c)$ .

Step 2. Set this equal to the virial-mass definition $M_\Delta=\dfrac{4}{3}\pi r_\Delta^3\,\Delta\,\rho_{\rm ref}=\dfrac{4}{3}\pi c^3 r_s^3\,\Delta\,\rho_{\rm ref}$ . Why: the virial condition is what closes the system — it says "the halo ends where the mean interior density drops to $\Delta\rho_{\rm ref}$ ", supplying the one extra equation needed to fix $\rho_s$ .

Step 3. Divide both sides by $4\pi r_s^3$ and cancel: $\dfrac{1}{3}c^3\Delta\rho_{\rm ref}=\rho_s\,g(c)$ . Why: the shared geometric factor $4\pi r_s^3$ drops out, showing the result cannot depend on the absolute size of the halo — only on the dimensionless $c$ .

Step 4. Solve for the characteristic density: $\rho_s=\dfrac{\Delta}{3}\rho_{\rm ref}\,\dfrac{c^3}{g(c)}$ . Because $g(c)$ grows only logarithmically while $c^3$ grows fast, $\rho_s$ climbs steeply with concentration — a $c=15$ halo has $\sim14\times$ the characteristic density of a $c=5$ halo. Why: this quantifies "early-forming halos are denser inside"; concentration is a fossil record of the background density at collapse. *The formula $\delta_c=\rho_s/\rho_{\rm crit}=(200/3)c^3/g(c)$ evaluates to $8.7\times10^3$ at $c=5$ , $7.0\times10^4$ at $c=12$ , and $1.2\times10^5$ at $c=15$ .*

Limiting cases

The NFW profile is an interpolation between two power laws, so its character is clearest in its asymptotes and in contrast with the models it replaced.

Case	Result	Meaning
Inner cusp ( $r\ll r_s$ )	$\rho\to \rho_s r_s/r\propto r^{-1}$	Density diverges but $M\propto r^2\to0$ : the cusp is integrable and holds negligible mass.
Outer envelope ( $r\gg r_s$ )	$\rho\to \rho_s r_s^3/r^3\propto r^{-3}$	Enclosed mass diverges only as $\ln r$ , forcing truncation at a virial radius $r_\Delta$ .
Pivot at scale radius	$\left.\mathrm{d}\ln\rho/\mathrm{d}\ln r\right\|_{r_s}=-2$	The halo momentarily matches a singular isothermal sphere; marks cusp-to-envelope crossover.
Rotation-curve peak	$r_{\max}=2.163\,r_s$	Fixed multiple of $r_s$ ; measuring $V_{\max}$ and its radius fixes the concentration.
Singular isothermal sphere (contrast)	$\rho_{\rm SIS}=\dfrac{\sigma^2}{2\pi G r^2}\propto r^{-2}$	Single power law with a flat rotation curve; NFW is shallower inside and steeper outside.
Einasto (refinement)	slope $=-2(r/r_{-2})^\alpha$ , finite centre	No central singularity; a continuously varying slope fits simulations slightly better.

Common trapThe NFW total mass is finite, and ρ_s is the central density

Both halves of that sentence are false, and both are common traps. The NFW total mass is infinite: because $\rho\propto r^{-3}$ in the outskirts, $M(r)$ grows without bound as $\ln(1+r/r_s)$ , so any statement of "the halo mass" is meaningless until you name a truncation radius and overdensity — quote $M_{200}$ (mean density $200\rho_{\rm crit}$ ), $M_{500}$ , or $M_{\rm vir}$ , and never assume two authors mean the same $\Delta$ . And $\rho_s$ is not the central density, which diverges; it is a normalization equal to four times the density at $r_s$ , tied to concentration through $\rho_s=(\Delta/3)\rho_{\rm ref}\,c^3/g(c)$ . A corollary: concentration $c$ and even the halo mass $M_\Delta$ shift when you change $\Delta$ , so " $c=10$ " is only meaningful once the overdensity convention is stated.

Worked examples

1 · A Milky-Way-scale halo from mass and concentration (fully worked)

A halo has $M_{200}=10^{12}\,M_\odot$ and concentration $c=12$ (overdensity $\Delta=200$ relative to $\rho_{\rm crit}$ ). Take $\rho_{\rm crit}=126\,M_\odot\,\mathrm{kpc}^{-3}$ (for $h=0.674$ ) and $G=4.30\times10^{-6}\,\mathrm{kpc}\,(\mathrm{km/s})^2 M_\odot^{-1}$ . Find $r_{200}$ , $r_s$ , the characteristic density $\rho_s$ , and the peak circular velocity $V_{\max}$ .

Step 1. Virial radius from $M_{200}=\frac{4}{3}\pi r_{200}^3(200\rho_{\rm crit})$ . Solving, $r_{200}=\left[\dfrac{3M_{200}}{4\pi\cdot200\rho_{\rm crit}}\right]^{1/3}=\left[\dfrac{3\times10^{12}}{4\pi\times200\times126}\right]^{1/3}\mathrm{kpc}=\mathbf{212\ kpc}$ .

Step 2. Scale radius from the concentration: $r_s=r_{200}/c=212/12=\mathbf{17.6\ kpc}$ .

Step 3. Characteristic density from (4). First $g(c)=\ln(1+12)-12/13=2.565-0.923=1.642$ . Then $\rho_s=\dfrac{200}{3}\rho_{\rm crit}\dfrac{c^3}{g(c)}=66.7\times126\times\dfrac{1728}{1.642}\,M_\odot\,\mathrm{kpc}^{-3}=\mathbf{8.8\times10^{6}\ M_\odot\,kpc^{-3}}$ , i.e. $\delta_c=\rho_s/\rho_{\rm crit}\approx\mathbf{7.0\times10^{4}}$ .

Step 4. Peak circular velocity. The rotation curve peaks at $r_{\max}=2.163\,r_s=38.1\,$ kpc, where the enclosed mass is $M(r_{\max})=4\pi\rho_s r_s^3\,g(2.163)$ with $g(2.163)=\ln(3.163)-2.163/3.163=1.152-0.684=0.468$ . Numerically $4\pi\rho_s r_s^3=4\pi(8.8\times10^6)(17.6)^3=6.0\times10^{11}\,M_\odot$ , so $M(r_{\max})=6.0\times10^{11}\times0.468=2.8\times10^{11}\,M_\odot$ . Then $V_{\max}=\sqrt{GM(r_{\max})/r_{\max}}=\sqrt{4.30\times10^{-6}\times2.8\times10^{11}/38.1}=\mathbf{179\ km/s}$ . For comparison the velocity at the virial radius is $V_{200}=\sqrt{GM_{200}/r_{200}}=\sqrt{4.30\times10^{-6}\times10^{12}/212}=143\,$ km/s, so $V_{\max}/V_{200}=1.26$ — a ratio set entirely by concentration. A $10^{12}\,M_\odot$ , $c=12$ halo has $r_{200}=212$ kpc, $r_s=17.6$ kpc, $\rho_s=8.8\times10^{6}\,M_\odot\,\mathrm{kpc}^{-3}$ , and $V_{\max}=179$ km/s at $38$ kpc — matching the Milky Way's dark halo.

2 · Why the rotation curve peaks at 2.163 r_s

Show that the NFW circular velocity is maximized at a fixed dimensionless radius, independent of halo mass or normalization, and locate it.

Step 1. From (6), $v_c^2\propto g(x)/x$ with $g(x)=\ln(1+x)-x/(1+x)$ ; the prefactor $4\pi G\rho_s r_s^3/r_s$ is a constant that cannot move the maximum. So the peak location depends on the shape function alone — it is universal.

Step 2. Set $\dfrac{\mathrm{d}}{\mathrm{d}x}\!\left[\dfrac{g(x)}{x}\right]=0$ , which gives $x\,g'(x)=g(x)$ . Differentiating $g$ , one finds the tidy result $g'(x)=\dfrac{x}{(1+x)^2}$ .

Step 3. Substitute to get the transcendental condition $\dfrac{x^2}{(1+x)^2}=\ln(1+x)-\dfrac{x}{1+x}$ . This has no closed-form root; solving numerically gives $x_{\max}=\mathbf{2.163}$ .

Step 4. Interpret: because $r_{\max}=2.163\,r_s$ is fixed, the observed radius of the rotation-curve peak directly measures $r_s$ , and combined with $r_{200}$ (from the total mass) it yields the concentration $c=r_{200}/r_s$ . This is why $V_{\max}$ and $r_{\max}$ are the workhorse observables for characterizing simulated and real halos alike. The NFW rotation curve always peaks at $r_{\max}=2.163\,r_s$ , making $V_{\max}$ a normalization-free concentration diagnostic.

Physics of the portrait

The NFW profile is unusual among the "named equations" of physics: it was not derived from a theory but discovered in a computer, an emergent regularity of nonlinear gravitational collapse that no one had predicted and that still lacks a complete first-principles explanation.

Origins · 1 minnavarro-frenk-white profile and how this was discovered

navarro-frenk-white profile

the universal two-parameter halo form found in cold-dark-matter simulations · cosmological structure formation and dark-matter halos

The profile is named for Julio Navarro, Carlos Frenk, and Simon White, who introduced it in two papers (1996, 1997) while studying the outcome of hierarchical structure formation in cold dark matter. Their route was computational rather than analytic: they evolved gravitational N-body simulations of halos forming through the collapse and merging of dark-matter clumps, then measured the spherically-averaged density profiles of the relaxed objects. The surprise was universality — halos spanning from dwarf-galaxy to galaxy-cluster masses, born from very different initial conditions, all fell onto the same two-parameter curve once rescaled by $r_s$ and $\rho_s$ . This built directly on the framework Simon White and Martin Rees had laid out in 1978, in which luminous galaxies condense inside a hierarchy of dark-matter halos, itself resting on the dynamical evidence for dark matter accumulated since Fritz Zwicky's 1933 Coma-cluster measurement and Vera Rubin's flat rotation curves of the 1970s. The NFW form quickly became the lingua franca of the field: it is how simulators report their halos, how observers fit lensing and rotation data, and how theorists model the Galactic potential. Its two enduring puzzles — why the inner slope is a cusp rather than a core, and why the Einasto form fits marginally better — remain active research, a reminder that a formula can dominate a field for thirty years while its origin is still being worked out.

History

1933

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

1970

— Rubin and Ford measure flat rotation curves in Andromeda, evidence for dark halos

1978

— White and Rees place galaxy formation inside a hierarchy of growing dark halos

1996

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

2004

— high-resolution simulations show the Einasto form fits relaxed halos slightly better

Going further

Texts. Binney & Tremaine, Galactic Dynamics, 2nd ed., §2.2 (spherical potentials) and §4.3 (equilibria) · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 7 (structure formation) and Ch. 11 (halo structure) · Carroll & Ostlie, An Introduction to Modern Astrophysics, galaxies and dark-matter chapters · Dodelson, Modern Cosmology, nonlinear structure and halos · Navarro, Frenk & White, ApJ 462 (1996) and 490 (1997) — the original papers. The key analytic results are: the enclosed-mass integral $\int_0^x u/(1+u)^2\,\mathrm{d}u=\ln(1+x)-x/(1+x)$ ; the potential at the centre $\Phi(0)=-4\pi G\rho_s r_s^2$ ; the slope running $-1\to-2\to-3$ with its $-2$ pivot at $x=1$ ; the rotation-curve peak condition $x^2/(1+x)^2=g(x)$ yielding $x_{\max}\approx2.163$ ; and the worked halo: $M_{200}=10^{12}M_\odot$ , $c=12$ gives $r_{200}=212$ kpc, $r_s=17.6$ kpc, $\delta_c\approx7.0\times10^4$ , $V_{\max}\approx179$ km/s. Concentration ranges, the mass–concentration relation, the core–cusp problem, and the Einasto improvement are qualitative summaries from the simulation literature; no non-standard numbers are asserted.

Prerequisites. Gravitational N-Body Problem · Collisionless Boltzmann Equation · Virial Theorem in Stellar Systems · Jeans Equations. Leads to. Galactic Potential Models · Orbits in Spherical Potentials · Jeans Theorem.

Ask the tutor: Ask me to fit an NFW profile to a rotation curve you have, convert between $M_{200c}$ and $M_{500c}$ , or work through why baryonic feedback might turn a cusp into a core.

Is the NFW profile derived from theory or just fitted to simulations?

Fitted — and that is the remarkable part. NFW is an empirical regularity: Navarro, Frenk & White ran N-body simulations of hierarchical collapse and found that the relaxed halos all followed the same two-parameter curve, across more than five decades in mass. There is still no complete first-principles derivation of why gravitational collapse produces exactly an $r^{-1}$ inner cusp and an $r^{-3}$ outer slope; partial arguments invoke violent relaxation, the shape of the initial power spectrum, and the mixing of material accreted at different epochs. This is unusual for a "textbook" formula, and it is why the profile is quoted as a fitting function with acknowledged residuals rather than as a law — and why the Einasto refinement, also empirical, is now often preferred.

What is concentration, and why do small halos have more of it?

Concentration is $c=r_\Delta/r_s$ , the ratio of the virial radius to the scale radius — a measure of how centrally packed a halo is. Simulations reveal a mass–concentration relation: lower-mass halos have higher $c$ (typically $c\sim10$ –15 for galaxies, falling to $c\sim4$ –5 for massive clusters). The physical reason is timing. Low-mass halos collapse earlier, when the mean density of the universe was higher, so their inner regions inherit that higher density and end up with a small $r_s$ relative to $r_\Delta$ . Massive clusters are still assembling today, from a diluter universe, so they are fluffier. Concentration is thus a fossil of a halo's formation epoch, and the relation is a genuine prediction of $\Lambda$ CDM that lensing and dynamics test directly.

What is the core–cusp problem?

NFW predicts an inner density cusp, $\rho\propto r^{-1}$ , so the rotation curve should rise steeply near the centre. But many dwarf and low-surface-brightness galaxies show rotation curves that rise gently, consistent with a constant-density core instead. This mismatch is the core–cusp problem, one of the main small-scale tensions of cold dark matter. There are two broad resolutions. Baryonic feedback: repeated gas inflows and supernova-driven outflows make the central potential fluctuate, dynamically heating the dark matter and flattening the cusp into a core — this works well in simulations of gas-rich dwarfs. Or new physics: self-interacting dark matter would thermalize the centre and naturally produce a core. Distinguishing them is an active observational program, since the very faintest dwarfs, with little gas, should retain cusps if feedback is the answer.

Why does the density diverge at the centre but the mass and potential stay finite?

Because the cusp is gentle. The density goes as $\rho\propto r^{-1}$ , which does blow up at $r=0$ , but the mass in a thin shell scales as $r^2\rho\propto r$ , so the enclosed mass $M(r)\propto r^2$ vanishes at the centre — the singular point contains no mass. The potential $\Phi(r)=-GM(r)/r+\ldots$ then approaches the finite value $\Phi(0)=-4\pi G\rho_s r_s^2$ . Contrast this with steeper cusps: a slope of $-2$ (isothermal) still gives finite central potential but a divergent central density, while a slope steeper than $-3$ would make even the enclosed mass diverge at the centre. NFW's $-1$ is comfortably in the integrable regime, so orbits through the centre are perfectly well-behaved.

When should I use Einasto instead of NFW?

Use Einasto when you need sub-percent accuracy against high-resolution simulations, or when the inner structure matters and you want to avoid the hard $r^{-1}$ cusp. The Einasto profile replaces NFW's fixed slopes with a slope that varies as a power of radius, $-2(r/r_{-2})^\alpha$ , adding one shape parameter $\alpha$ (typically $0.15$ – $0.25$ , mildly increasing with halo mass). It has a finite central density (no singularity) and fits relaxed simulated halos noticeably better, especially in the inner regions and the transition zone. For most analytic work, order-of-magnitude estimates, and lensing/dynamics fits where two parameters suffice, NFW remains the default — it has closed-form mass, potential, and rotation curve, whereas Einasto's mass involves the incomplete gamma function. In short: NFW for tractability, Einasto for fidelity.

Does the "200" in M₂₀₀ and r₂₀₀ mean anything physical, or is it a convention?

It is a convention with a loose physical motivation. The virial radius is defined as the radius inside which the mean density equals $\Delta$ times a reference cosmic density, and $\Delta=200$ is a common choice because the spherical-collapse model predicts a virialized overdensity of roughly $18\pi^2\approx178$ (often rounded to 200) relative to the critical density in an Einstein–de Sitter universe. But the choice is not unique: people also use $\Delta=500$ (better matched to X-ray-observable regions of clusters), $\Delta=200\bar\rho_m$ relative to the mean matter density, or a redshift-dependent $\Delta_{\rm vir}(z)$ from the full collapse calculation. Because $M_\Delta$ , $r_\Delta$ , and even $c$ all shift with the choice, you must always state which $\Delta$ and which reference density you mean — comparing an $M_{200c}$ to an $M_{500c}$ without conversion is a common and avoidable error.