The NFW Dark Matter Profile

The surprising universality

Run a cosmological N-body simulation: seed a patch of universe with cold dark matter particles tracing the tiny density ripples left by inflation, then let gravity do its work for 13.8 billion years. Overdense regions pull in their neighbours, merge, and relax into bound clumps called halos. The remarkable discovery — made in the mid-1990s — is that no matter how massive the halo, no matter when it formed, and across a wide range of cosmological parameters, the spherically averaged density runs along the same curve. A dwarf-galaxy halo of 10⁹ M☉ and a galaxy-cluster halo of 10¹⁵ M☉ have density profiles that are scaled copies of one another.

That curve is the NFW profile. Its defining feature is a cusp: rather than flattening to a constant central density (a "core"), the density keeps climbing toward the centre as ρ ∝ r⁻¹. Move outward and the slope steepens continuously, passing through −2 at a characteristic radius r_s and tending toward −3 far out — the falling "skirt." A single transition scale separates cusp from skirt, and once you specify the halo's total mass and that one transition scale, the entire density curve is fixed. A two-parameter family describing structures spanning twenty decades in mass is exactly the kind of simplicity physicists hope for and rarely get.

The governing equation

The NFW density profile is written compactly as

ρ(r) = ρ_s / [ (r / r_s) (1 + r / r_s)² ]

Here r_s is the scale radius and ρ_s is a characteristic density (four times the density at r = r_s). Reading off the two limits:

r ≪ r_s :  ρ → ρ_s (r_s / r)        →  ρ ∝ r⁻¹   (the cusp)
r ≫ r_s :  ρ → ρ_s (r_s / r)³        →  ρ ∝ r⁻³   (the skirt)

The logarithmic slope d(ln ρ)/d(ln r) runs smoothly from −1 to −3, and equals exactly −2 at r = r_s — which is the operational definition of the scale radius. Integrating the density in spherical shells gives the enclosed mass

M(<r) = 4π ρ_s r_s³ [ ln(1 + r/r_s) − (r/r_s) / (1 + r/r_s) ]

This is the single most-used NFW formula in practice — it appears in every rotation-curve fit and lensing mass estimate. Note its behaviour at large r: the bracket grows like ln(r/r_s), so the total mass diverges logarithmically. The halo therefore must be truncated, conventionally at the virial radius r_vir.

The halo is normalised by a virial overdensity condition. The virial radius r_vir (sometimes r₂₀₀) encloses a mean density Δ times the critical density of the universe ρ_crit, with Δ ≈ 200 a common convention (or the redshift-dependent value from spherical collapse, Δ ≈ 178 in an Einstein–de Sitter universe). Then the characteristic density is tied to the concentration c = r_vir/r_s by

ρ_s = (Δ/3) ρ_crit · c³ / [ ln(1 + c) − c/(1 + c) ]

So the three derived quantities ρ_s, r_s and r_vir collapse onto just two free numbers: a mass (which sets r_vir via M_vir = (4/3)π Δ ρ_crit r_vir³) and a concentration c.

Concentration: the one shape knob

Because the NFW shape is fixed, the only thing that distinguishes a "fat" halo from a "centrally peaked" one is the concentration parameter

c = r_vir / r_s

A high-c halo packs more of its mass deep inside r_s and has a sharply rising rotation curve; a low-c halo is more spread out. Crucially, concentration anti-correlates with mass. Low-mass halos collapsed earlier in cosmic history, when the mean density of the universe was higher, so they "froze in" a denser core and ended up more concentrated. The empirical relation, calibrated on simulations, is roughly

c ≈ 9 (M_vir / 10¹² M☉)^(−0.1) · (1 + z)⁻¹   (approximate, z = redshift)

This mass-concentration relation, refined by groups such as Bullock et al. (2001) and Dutton & Macciò (2014), is gentle — only a factor of a few across the full mass range — which is itself part of the "universality." Typical values: galaxy clusters c ≈ 4–7, Milky-Way-scale halos c ≈ 10–15, dwarf-galaxy halos c ≈ 15–25.

The key numbers

For our own Galaxy, the dark matter halo has a virial mass of roughly M₂₀₀ ≈ 1.0–1.5 × 10¹² M☉, a virial radius r₂₀₀ ≈ 200–250 kpc, and a concentration c ≈ 12, putting the scale radius near r_s ≈ 16–20 kpc — comfortably outside the Sun's 8.2 kpc orbit. The local dark matter density at the Sun's position is measured to be about ρ_⊙ ≈ 0.3–0.5 GeV cm⁻³, or roughly 0.008–0.013 M☉ pc⁻³ — about one proton mass per cubic centimetre's worth of dark matter, the headline input for every direct-detection experiment.

Halo	M_vir	r_vir	Concentration c	Scale radius r_s
Ultra-faint dwarf	~10⁹ M☉	~30 kpc	~20	~1.5 kpc
Classical dwarf (Fornax)	~10¹⁰ M☉	~60 kpc	~16	~4 kpc
Milky Way	~1.2 × 10¹² M☉	~220 kpc	~12	~18 kpc
Galaxy group	~10¹³ M☉	~500 kpc	~8	~60 kpc
Galaxy cluster (Coma-like)	~10¹⁵ M☉	~2 Mpc	~5	~400 kpc

Notice how r_s scales up with halo size while c slides down — exactly the mass-concentration trend. In every row the cusp lives well inside r_s and the r⁻³ skirt dominates the outer halo.

How NFW is measured in real halos

NFW was born from simulations, but it is tested against several independent observational probes, each sensitive to a different radial range:

• Galaxy rotation curves. The circular velocity v_c(r) = √[GM(<r)/r] is reconstructed from H I and Hα Doppler shifts. An NFW halo predicts a velocity that rises, peaks near ~2.16 r_s, and slowly declines. This is the cleanest probe of the inner-to-mid halo and the arena of the cusp-core debate.
• Gravitational lensing. Both strong lensing (multiple images, Einstein rings) and weak lensing (statistical galaxy-shape distortion) map the projected mass. Stacking weak-lensing signals around millions of galaxies — as in the Dark Energy Survey and KiDS — recovers the NFW shape out to and beyond r_vir, where rotation curves cannot reach.
• X-ray hydrostatic mass. In clusters, the temperature and density of the hot intracluster gas give M(<r) under hydrostatic equilibrium. Chandra and XMM-Newton clusters are routinely fit with NFW, yielding the low concentrations expected for the most massive halos.
• Stellar kinematics and streams. The velocity dispersions of dwarf-galaxy stars, and the shapes of disrupted stellar streams in the Milky Way halo (now mapped exquisitely by ESA's Gaia mission), constrain the enclosed mass at specific radii.

Across these probes, NFW (or its close cousin Einasto) fits to roughly the ~10 percent level over the radii that matter — an impressive achievement for a two-parameter formula.

Worked example: the Milky Way's dark mass inside the Sun

How much dark matter lies inside the Sun's orbit at R₀ = 8.2 kpc? Take an NFW halo with c = 12, r_vir = 220 kpc, so r_s = r_vir/c = 18.3 kpc, and a virial mass M_vir = 1.2 × 10¹² M☉. First we need ρ_s r_s³ = M_vir / {4π [ln(1+c) − c/(1+c)]}. The bracket is

μ(c) = ln(1 + 12) − 12/13 = ln(13) − 0.923 = 2.565 − 0.923 = 1.642

So 4π ρ_s r_s³ = M_vir / μ(c) = 1.2 × 10¹² / 1.642 = 7.31 × 10¹¹ M☉. Now evaluate the enclosed mass at r = R₀ = 8.2 kpc, where x = r/r_s = 8.2 / 18.3 = 0.448:

M(<R₀) = 4π ρ_s r_s³ [ ln(1 + x) − x/(1 + x) ]
       = 7.31 × 10¹¹ × [ ln(1.448) − 0.448/1.448 ]
       = 7.31 × 10¹¹ × [ 0.370 − 0.309 ]
       = 7.31 × 10¹¹ × 0.0610
       ≈ 4.5 × 10¹⁰ M☉

So about 4.5 × 10¹⁰ solar masses of dark matter sit inside the Sun's orbit — comparable to the stellar disk mass of ~5 × 10¹⁰ M☉, which is exactly why the rotation curve stays flat rather than falling Keplerian beyond the visible disk. The corresponding dark-matter contribution to the circular speed is v_c = √(G M / R₀) ≈ √(4.30 × 10⁻⁶ × 4.5 × 10¹⁰ / 8.2) kpc-based ≈ 154 km/s, which adds in quadrature with the baryonic disk and bulge to give the observed ~230 km/s.

Discovery: Navarro, Frenk and White

The profile is named for Julio Navarro, Carlos Frenk and Simon White, who announced it in two landmark papers. The first, "The Structure of Cold Dark Matter Halos" (Navarro, Frenk & White 1996, ApJ 462, 563), examined a handful of high-resolution cluster and galaxy halos and noticed they shared a common shape. The second and definitive paper, "A Universal Density Profile from Hierarchical Clustering" (1997, ApJ 490, 493), simulated halos across four decades of mass in several cosmologies and showed the form was genuinely universal — and that all the diversity collapsed onto the single mass-concentration relation.

The work built on the 1980s development of cold dark matter cosmology (Peebles; Blumenthal, Faber, Primack & Rees 1984) and on White & Rees's 1978 picture of hierarchical galaxy formation inside dark halos. It was enabled by the explosion in N-body computing power: the GADGET code (Volker Springel) and successors let later teams push to billions of particles in the Millennium (2005), Aquarius (2008), and Via Lactea II (2008) simulations, which refined the profile and exposed the subtle inner deviations now better described by the Einasto form. Frenk and White — alongside Marc Davis and George Efstathiou, their collaborators on the foundational 1985 cold-dark-matter simulations — went on to share the 2011 Gruber Cosmology Prize for that body of computational work.

The cusp-core controversy

NFW's sharpest confrontation with data is in the centres of small galaxies. The profile demands a cusp — density rising without bound as r⁻¹ toward r = 0. But the rotation curves of many dwarf irregulars and low-surface-brightness galaxies rise gently in their inner kiloparsec, implying a roughly constant-density core, better fit by phenomenological forms like the Burkert or pseudo-isothermal profile. This is the "cusp-core problem," and it has resisted easy resolution for over two decades.

Profile	Inner slope	Centre behaviour	Origin	Best for
NFW	ρ ∝ r⁻¹	Cusp	Pure CDM N-body	Clusters, mid-halo, lensing
Einasto	slope → 0 gradually	Mild cusp	High-res CDM N-body	Modern high-resolution halos
Burkert	ρ ∝ const	Core	Empirical fit to dwarfs	LSB / dwarf rotation curves
Pseudo-isothermal	ρ ∝ const	Core	Empirical / kinematic	Inner galaxy kinematics
Cored NFW (coreNFW)	tunable 0 to −1	Feedback-flattened cusp	Hydro sims (FIRE, NIHAO)	Reconciling CDM with cores
SIDM	ρ ∝ const inner	Thermalised core	Self-interacting DM	Diverse-rotation-curve solutions

Two broad classes of resolution are debated. The first keeps cold dark matter but invokes baryonic feedback: repeated bursts of supernova-driven gas outflow rapidly change the central gravitational potential, gravitationally "heating" the dark matter and turning a cusp into a core. Hydrodynamic simulations (FIRE, NIHAO) reproduce this for dwarfs in a specific stellar-mass window. The second class modifies the dark matter, most prominently self-interacting dark matter (SIDM), in which dark-matter particles scatter off one another, conduct heat inward, and thermalise the centre into a core. Distinguishing feedback from new physics is an active observational frontier.

Variants and related forms

• Einasto profile. ρ(r) ∝ exp{−(2/α)[(r/r_−2)^α − 1]} with α ≈ 0.16–0.18. Has no fixed inner cusp; the slope flattens gradually toward the centre. Fits modern high-resolution halos slightly better than NFW at the cost of one extra parameter.
• Generalised NFW (gNFW). ρ ∝ (r/r_s)^(−γ) (1 + r/r_s)^(γ−3), with the inner slope γ now a free parameter (γ = 1 recovers NFW). Used when the data want to constrain the central slope directly.
• Hernquist profile. ρ ∝ 1 / [(r/a)(1 + r/a)³], identical r⁻¹ cusp but a steeper r⁻⁴ skirt, so its total mass is finite — convenient for analytic galaxy models and stellar bulges.
• Burkert / pseudo-isothermal. Cored alternatives with ρ → const at the centre, motivated empirically by dwarf-galaxy rotation curves.
• Adiabatically contracted NFW. When baryons cool and settle into the centre, their added gravity pulls dark matter inward (Blumenthal et al. 1986), steepening the inner profile — the opposite of feedback-driven coring.

Common misconceptions and subtleties

• "NFW has a constant central density." No — it diverges. The density rises without bound as r → 0 (the cusp). What stays finite is the enclosed mass, because the volume shrinks faster than the density rises.
• "r_s is where the density is highest." The density is monotonically decreasing everywhere; r_s is simply where the logarithmic slope equals −2, the inflection between cusp and skirt. The rotation-curve peak, by contrast, sits at about r ≈ 2.16 r_s.
• "The NFW total mass is well defined." Only with a truncation radius. Integrated to infinity the mass diverges logarithmically because of the r⁻³ skirt; you must adopt r_vir (or r₂₀₀) to get a finite number.
• "Concentration measures how massive a halo is." Concentration is a shape parameter, not a mass. It happens to anti-correlate with mass — small halos are more concentrated because they formed earlier — but two halos of identical mass can have different concentrations reflecting different formation histories.
• "NFW is exact." It is a fitting formula accurate to ~10 percent, not a derived law. Real halos are triaxial, lumpy with substructure, and show the gentler Einasto inner curvature. The spherically averaged NFW form is a remarkably good approximation, not the truth.