Dark Matter Halos and the Halo Mass Function
A dark matter halo is a gravitationally bound clump of cold, collisionless matter that has decoupled from the cosmic expansion, turned around, collapsed, and settled into approximate virial equilibrium. Every galaxy — from a faint dwarf to the brightest cluster central — is thought to sit inside one, its potential well the reservoir into which gas falls, cools, and lights up. So the census of halos is, to first approximation, the census of the places where galaxies can form, and the central object of that census is the halo mass function dn/dM: the comoving number density of halos per unit mass. The theory that predicts it bridges two regimes that look nothing alike. On one side is a smooth, nearly Gaussian primordial density field, fixed by inflation and the transfer function and known statistically through its variance σ²(M); on the other is a discrete population of lumpy, nonlinear, virialized objects. The bridge is a single threshold: follow one overdense sphere through the spherical-collapse model and it collapses precisely when its linearly-extrapolated contrast reaches δ_c ≈ 1.686, leaving a halo of mean interior overdensity Δ_vir ≈ 178. Press and Schechter's insight was that you can then count halos by asking what fraction of the Gaussian field sits above that barrier. The pay-off is exponential leverage. Because the most massive clusters are the rarest high peaks of the field, their abundance is exponentially sensitive to the amplitude σ₈ and to the growth history set by dark energy — which is why counting clusters rivals the CMB and baryon acoustic oscillations as a cosmological probe. Beyond counting, halos share a universal Navarro–Frenk–White internal structure, so the mass function, halo bias, and density profile together form the "halo model" of modern weak lensing and galaxy clustering.
- Dark matter halos are quasi-equilibrium, gravitationally bound overdensities formed by nonlinear collapse of the primordial field — the sites where gas cools and galaxies form.
- Spherical collapse gives a universal barrier: a region collapses and virializes once its linearly-extrapolated contrast reaches δ_c ≈ 1.686, producing a halo of mean overdensity Δ_vir = 18π² ≈ 178 in an Einstein–de Sitter universe.
- The mass function dn/dM comes from counting Gaussian field points above that barrier — the Press–Schechter formalism — and depends on mass only through the peak height ν = δ_c/σ(M,z).
- Rare massive halos have ν ≫ 1 and are exponentially suppressed, e^(−ν²/2); low-mass halos are abundant as a near power law — the content of hierarchical structure formation.
- N-body calibrations (Sheth–Tormen, Tinker) correct PS for ellipsoidal collapse, adding more rare massive halos and matching simulations to ≈ 10 %; cluster abundance is a sharp probe of σ₈, Ω_m, and dark energy.
- Key idea: map the continuous Gaussian primordial density field onto a discrete population of bound objects — a region of mass M collapses when its smoothed linear contrast exceeds δ_c ≈ 1.686, and counting such regions in a field of variance σ²(M) gives dn/dM ∝ (ρ̄/M²) ν f(ν) |dlnσ/dlnM| with peak height ν = δ_c/σ.
Map the continuous, Gaussian primordial density field onto a discrete population of bound objects: a region of mass M forms a halo when its linear density contrast smoothed on scale M exceeds the collapse threshold δ_c ≈ 1.686, and counting such regions in a Gaussian field of variance σ²(M) gives the mass function dn/dM ∝ (ρ̄/M²) ν f(ν) |dlnσ/dlnM| with peak height ν = δ_c/σ.
Where this lives today
Counting halos is now a front-line cosmological measurement. X-ray surveys with eROSITA and Sunyaev–Zel'dovich surveys with the South Pole Telescope and the Atacama Cosmology Telescope detect thousands of clusters and fit their abundance dn/dM(z) against Tinker-type mass functions to constrain the σ₈–Ω_m plane — the same combination the CMB constrains, so the comparison tests the growth of structure since recombination and is one arena of the mild "S8 tension." The mass calibration comes from stacked weak-lensing shear, which measures the projected NFW profile directly and without dynamical assumptions; the Dark Energy Survey, the Euclid satellite (launched 2023), and the Vera Rubin Observatory's LSST are pushing this to millions of galaxies. At the other extreme, the low-mass tail sets the abundance of the first minihalos that host Population III stars and drive reionization — a regime where the dark-matter free-streaming cutoff (cold vs. warm) leaves a signature that JWST's earliest-galaxy counts and 21 cm experiments are beginning to probe.
The physical picture
Picture the primordial density field, smoothed on mass scale M, as a gently rolling landscape whose typical height above the mean is one standard deviation, σ(M). Gravity is a ratchet on the hills: any patch whose smoothed contrast climbs above a critical height δ_c eventually stops expanding, turns around, and collapses into a bound halo. Bigger halos require averaging over larger patches, and the field is smoother on large scales, so σ falls with M — massive halos are therefore rare, tall peaks and are exponentially scarce. The whole problem reduces to one ratio, the peak height ν = δ_c/σ(M), the number of standard deviations a collapsing region of mass M sits above the mean.
As the universe ages, structure grows: D(z) increases, σ rises, ν = δ_c/σ falls, and progressively rarer, more massive peaks are pulled across the barrier — that single sentence is hierarchical structure formation, small halos first, big ones assembled late by merging. One famous subtlety is the "factor of 2": naive Gaussian counting captures only half the collapsed mass, because an underdense patch can be swallowed inside a larger overdense region (the cloud-in-cloud problem), which the excursion-set reformulation recovers exactly. Stitched across mass scales, this gives the two-part shape of the mass function: a shallow power-law bank at low mass and a cliff-like exponential cutoff at high mass.
Once formed, halos are not featureless blobs. Violent relaxation drives them to a nearly universal Navarro–Frenk–White profile with an inner ρ ∝ r⁻¹ cusp and an outer ρ ∝ r⁻³ falloff, and the concentration — how centrally packed the halo is — records when it assembled, since an early-forming halo inherited a denser universe. The statistical population (mass function and bias) and the internal structure (profile and concentration) are the two halves of the halo model.
The core equations
Everything begins with the variance of the smoothed field, which turns the continuous power spectrum P(k) into a single number attached to a mass scale. We smooth the linear density field with a top-hat window of comoving radius R enclosing mass M, and integrate P(k) against the window in Fourier space.
where: σ²(M) — variance of the linear density contrast smoothed on mass scale M; P(k) — linear matter power spectrum; W̃(kR) — Fourier transform of the top-hat window of radius R; ρ̄ — mean comoving matter density; R — smoothing radius enclosing mass M. Because larger M means larger R and the window kills small-scale power, σ falls monotonically with M — the reason massive halos are rare. The single most important derived quantity is the peak height ν: it converts a mass and redshift into "how many σ above the mean," the only combination the multiplicity depends on. All redshift evolution enters through the growth factor D(z), which rescales σ uniformly.
where: ν — peak height, the number of standard deviations needed to reach the barrier; δ_c — linear collapse threshold, ≈ 1.686; σ(M,z) — rms fluctuation at redshift z; D(z) — linear growth factor. As z decreases D grows, σ rises, and ν falls, so the exponential cutoff sweeps toward higher mass — halos grow. The general mass function is a statement of mass conservation dressed as a change of variables: every mass element belongs to some halo, so the fraction of mass in halos of a given ν is a universal multiplicity ν f(ν), and the number density follows by the chain rule.
where: dn/dM — comoving number density of halos per unit mass; ν f(ν) — multiplicity function, the mass fraction in halos per unit ln ν; |dlnσ/dlnM| — logarithmic slope of the variance (always negative), the Jacobian from σ to M. Read it as: abundance = (density / mass²) × (collapse odds) × (how fast σ changes with mass). The Press–Schechter multiplicity is what you get from raw Gaussian statistics plus the factor-of-2 fix — the first closed form, and it already captures the exponential tail correctly.
where: ν f(ν) — the PS multiplicity, normalized so that ∫ f dν = 1 once the factor of 2 is included. The e^(−ν²/2) is the Gaussian tail; the leading ν is the Jacobian remnant. It overpredicts low-mass halos and underpredicts the rarest clusters. The Sheth–Tormen multiplicity generalizes PS to a "moving barrier" for triaxial (ellipsoidal) collapse: the same exponential character with a rescaled threshold and a low-ν upturn, matching simulations far better — crucially, it adds more rare massive halos than PS.
where: a ≈ 0.707 — lowers the effective barrier and broadens the distribution; p ≈ 0.3 — controls the low-mass (low-ν) enhancement; A ≈ 0.322 — normalization enforcing that all mass is in halos. Setting a = 1, p = 0, A = ½ recovers Press–Schechter exactly. The barrier δ_c is not a fitted number — it drops out of the spherical-collapse model as an exact analytic constant in an Einstein–de Sitter universe, and is remarkably insensitive to cosmology.
where: δ_c — the linearly-extrapolated density contrast a spherical perturbation would have at the instant it actually collapses. Real collapse reaches a formally infinite nonlinear contrast; δ_c is its linear-theory shadow — precisely the quantity Gaussian statistics of the initial field can supply. The companion constant is the virial overdensity, the mean density interior to a freshly virialized halo relative to the background — the number that fixes the edge of a halo and hence its mass.
where: Δ_vir — mean halo overdensity relative to background matter at virialization; ρ̄_halo — mean density inside the virial radius; ρ̄_m — background matter density. In ΛCDM (Ω_m ≈ 0.3) the Bryan–Norman fit raises this to ≈ 337 relative to mean matter — why halo mass is definition-dependent. Halos share a universal internal profile, the two-parameter NFW form, whose single shape parameter is the concentration c — the ratio of the outer radius to the scale radius where the slope passes through −2.
where: ρ_s — characteristic density; r_s — scale radius dividing the inner r⁻¹ cusp from the outer r⁻³ decline; c — concentration (higher c = earlier formation); r₂₀₀ — radius enclosing 200 times the critical density. This profile is what stacked weak lensing and rotation curves actually measure. Finally, halos are biased tracers: their overdensity is an amplified copy of the matter overdensity, rising steeply for rare peaks. The peak-background split gives the Press–Schechter linear bias in closed form.
where: b(ν) — linear halo bias relating halo to matter overdensity, δ_h = b δ_m. At ν = 1 (an M∗ halo) the bias is exactly 1; the rarest clusters (ν ≈ 3–4) have b ≈ 6–7, clustering far more strongly than the mass itself.
The goal of this derivation is to obtain the two constants on which the entire mass function rests — the collapse barrier δ_c ≈ 1.686 and the virial overdensity Δ_vir ≈ 178 — from nothing but Newtonian gravity applied to a single uniform overdense sphere in an Einstein–de Sitter background. We treat the sphere as its own closed sub-universe, solve its Friedmann equation parametrically, find when it turns around and collapses, and ask what value linear theory would have assigned to its contrast at that nonlinear moment — which is exactly what δ_c is.
Step 1. A closed dust sphere obeys the same cycloid solution as a closed Friedmann universe, r = A(1 − cos θ), t = B(θ − sin θ), with A³ = GMB² relating the amplitudes to the enclosed mass M and development angle θ. Why: by Birkhoff's theorem the sphere feels only its interior mass, so it evolves as a closed cosmology.
Step 2. Maximum radius — turnaround — is at θ = π, where r_ta = 2A and t_ta = πB; the sphere then symmetrically recollapses to r = 0 at θ = 2π, so t_coll = 2 t_ta. Why: cycloid symmetry pins the collapse time to twice the turnaround time — the fact that later fixes the barrier.
Step 3. Comparing the sphere to the EdS background (which dilutes as ρ̄ ∝ t⁻²) at turnaround gives the true nonlinear overdensity 1 + δ_ta = 9π²/16 ≈ 5.55. Why: collapse is genuinely nonlinear — the sphere is already 5.5× the background by turnaround.
Step 4. Expand the cycloid for small θ and propagate with the EdS linear law δ ∝ t^(2/3); extrapolating to the collapse time gives δ_lin(t_coll) = (3/20)(12π)^(2/3). Why: the statistics know only the linear field, so we need its value at collapse, not the real contrast.
Step 5. Evaluate: (3/20)(12π)^(2/3) = 0.15 × (37.70)^(2/3) = 0.15 × 11.24 = 1.686, the barrier δ_c. Why: this is the height a smoothed linear fluctuation must clear to have collapsed by that redshift.
Step 6. Energy conservation and the virial theorem give r_vir = ½ r_ta; folding in the density growth to virialization yields Δ_vir = 18π² ≈ 178. Why: this sets where a halo ends, and hence its mass. *Note: (3/20)(12π)^(2/3) = 1.6865 and 18π² = 177.65 to the quoted precision; the cycloid symmetry t_coll = 2 t_ta and 1 + δ_ta = 9π²/16 = 5.552 are exact EdS results.*
Go deeper · Derivation B · 1 minCounting peaks: the Press–Schechter mass function
The aim here is to derive the closed-form mass function dn/dM by literally counting how many points of a Gaussian random field sit above the collapse barrier, and to confront honestly the notorious factor of 2 that this counting is missing. We assume the smoothed contrast δ_M is Gaussian with variance σ²(M), identify the collapsed mass fraction with a barrier-crossing probability, patch in the factor of 2, and differentiate to turn a mass fraction into a number density — a result that depends on M only through ν = δ_c/σ.
Step 1. For a Gaussian field the fraction of points with δ_M > δ_c is P(>δ_c) = ½ erfc(δ_c / √2 σ), interpreted as the mass fraction bound in halos of mass ≥ M. Why: a point above the barrier at scale M is taken to lie inside a collapsed region of at least that mass — the leap from field statistics to object counts.
Step 2. As M → 0, σ → ∞ and P → ½: only half the mass is ever accounted for. Why: the cloud-in-cloud problem — an underdense point can sit inside a larger collapsed region, and thresholding misses exactly those.
Step 3. Multiply by 2 by hand, so F(>M) = erfc(δ_c / √2 σ) now tends to 1 as M → 0 — all mass is in halos. Why: the excursion-set random walk against an absorbing barrier (sharp-k filter) later shows this factor is exact, not a fudge.
Step 4. Convert to number density: each halo carries mass M, so dn/dM = (ρ̄/M) |dF/dM|, and the chain rule through σ(M) produces the Jacobian |dlnσ/dlnM|. Why: dividing the mass fraction by the mass per halo and differentiating gives a differential number count.
Step 5. Assembling gives dn/dM = √(2/π) (ρ̄/M²)(δ_c/σ) |dlnσ/dlnM| exp(−δ_c²/2σ²), i.e. ν f(ν) = √(2/π) ν e^(−ν²/2). Why: the dependence collapses onto ν alone, giving one near-universal curve across twenty decades of mass. *Note: the erfc normalization and the factor of 2 are the historically canonical forms, and the multiplicity integrates to unity by construction.*
Go deeper · Derivation C · 1 minThe peak-background split and why clusters cluster
This derivation explains why halos are biased — why the rarest clusters trace the cosmic web far more strongly than the dark matter itself — and produces the closed-form Press–Schechter bias b(ν) = 1 + (ν² − 1)/δ_c. The idea is to split the field into a long-wavelength background and short-wavelength fluctuations, recognize that the background shifts the local collapse threshold, and read the halo response off the mass function's sensitivity to that shift.
Step 1. Write δ = δ_ℓ + δ_s, a slowly varying long-wavelength background δ_ℓ plus short-wavelength peaks δ_s that actually collapse. Why: separating scales lets us treat the large-scale environment as a locally modified cosmology.
Step 2. A positive background δ_ℓ means a peak needs less of its own amplitude to reach δ_c, so the local barrier is lowered to δ_c − δ_ℓ. Why: collapse cares only about total contrast, so overdense environments cross the threshold preferentially.
Step 3. The local abundance therefore responds to δ_ℓ; expanding to first order, the fractional overabundance is b × δ_ℓ, and for PS the derivative gives b(ν) = 1 + (ν² − 1)/δ_c. Why: differentiating e^(−ν²/2) at the shifted barrier brings down a factor ∝ ν², so rare peaks respond most violently.
Step 4. Check regimes: ν = 1 gives b = 1 (M∗ halos trace mass), while ν ≈ 3.4 (a 10¹⁵ M⊙ cluster) gives b ≈ 7. Why: this is why cluster surveys see a filamentary sky, and why bias must be modeled to read cosmology from clustering. *Note: b(ν = 1) = 1 exactly, and b(ν = 3.37) = 1 + (3.37² − 1)/1.686 = 7.14, as expected from the peak-background split of Press–Schechter.*
Limiting cases
The mass function and its ingredients simplify cleanly in the regimes that matter observationally:
| Case | Result | Meaning | ||
|---|---|---|---|---|
| Rare massive halos (ν ≫ 1) | dn/dlnM ∝ exp(−δ_c²/2σ²) | Gaussian-tail cutoff; cluster abundance is exponentially sensitive to σ₈ and D(z) — the sharpest lever in cluster cosmology. | ||
| Low-mass halos (ν ≪ 1) | dn/dM ∝ M⁻² \ | dlnσ/dlnM\ | Near power-law regime; abundance rises toward low mass, regulated by the free-streaming cutoff and baryonic feedback. | |
| Einstein–de Sitter (Ω_m = 1) | δ_c = (3/20)(12π)^(2/3) ≈ 1.686, Δ_vir = 18π² ≈ 178 | Exact analytic collapse constants; the cosmology dependence of δ_c is negligible at the percent level. | ||
| ΛCDM at z = 0 | Δ_vir ≈ 337 (rel. to ρ̄_m), δ_c ≈ 1.675 | Bryan–Norman fit: with Ω_m ≈ 0.3 the virial overdensity relative to mean matter rises above 178; δ_c shifts only slightly. |
Worked examples
1 · The collapse constants from pure geometry (confidence builder)
Derive the numerical values of the collapse barrier δ_c, the turnaround overdensity, and the virial overdensity Δ_vir for an Einstein–de Sitter universe, using only the spherical-collapse results.
Step 1. The barrier is δ_c = (3/20)(12π)^(2/3). Evaluate the inside: 12π = 37.699, and 37.699^(2/3) = (3.3522)² = 11.237, so δ_c = 0.15 × 11.237 = 1.686 — the height a smoothed linear fluctuation must reach to have collapsed.
Step 2. The nonlinear overdensity at turnaround is 1 + δ_ta = 9π²/16 = 88.826/16 = 5.55: a proto-halo is already 5.5× the background when it stops expanding — collapse is fully nonlinear, and δ_c ≈ 1.686 is only its linear-theory shadow.
Step 3. The virial overdensity is Δ_vir = 18π² = 177.7 ≈ 178. A virialized EdS halo is 178× the mean matter density; the ΛCDM Bryan–Norman fit lifts this to ≈ 337 with Ω_m ≈ 0.3. δ_c = 1.686, 1 + δ_ta = 5.55, Δ_vir = 178 — three exact numbers, no free parameters.
2 · Peak height, rarity, and redshift evolution (fully worked)
Take σ(M,0) ≈ 2.0 for a galaxy halo of M = 10¹² M⊙ and σ(M,0) ≈ 0.5 for a cluster of M = 10¹⁵ M⊙. Find each peak height, compare their Gaussian suppressions, and see what happens to the cluster at z = 1 (use the Einstein–de Sitter growth D ∝ a, so σ(z) = σ(0)/(1+z)).
Step 1. Peak heights: for the galaxy, ν = δ_c/σ = 1.686/2.0 = 0.84 — a sub-σ fluctuation, so galaxy-mass halos are common. For the cluster, ν = 1.686/0.5 = 3.37 — a rare 3.4σ peak.
Step 2. Gaussian suppression of the tail: e^(−ν²/2). Galaxy: e^(−0.35) = 0.70. Cluster: e^(−5.68) = 3.4 × 10⁻³. The ratio of exponential factors alone is ≈ 200, so before any prefactors the cluster tail is already hundreds of times more suppressed — the cliff in the second figure.
Step 3. Now push the cluster to z = 1. With D ∝ a, σ halves to 0.25, so ν = 1.686/0.25 = 6.75 and the suppression becomes e^(−22.7) ≈ 1.3 × 10⁻¹⁰. A 10¹⁵ M⊙ cluster is essentially absent at z = 1 in this estimate — massive halos are assembled late, the quantitative statement of hierarchical growth.
Step 4. Check against Sheth–Tormen and bias at the cluster scale: at ν = 3.37, ST/PS ≈ 2.2 (ellipsoidal collapse adds rare halos), and b = 1 + (3.37² − 1)/1.686 = 7.1 — clusters are ≈ 7× more contrasted than the dark matter. Galaxy ν ≈ 0.84 (common), cluster ν ≈ 3.37 rising to 6.75 at z = 1, bias ≈ 7.
Physics of the portrait
The halo mass function is less the invention of one person than the convergence of a generation's ideas about how a smooth beginning becomes a lumpy universe — measured, modeled, and finally calibrated on the largest computers of each decade.
Origins · 1 minhalo mass function and how this was discovered
The halo mass function is one of those results that arrived before the tools to justify it. In 1972 James Gunn and J. Richard Gott followed a single overdense sphere through turnaround and collapse, giving the spherical-collapse picture and the constants δ_c and Δ_vir that still anchor the subject. Two years later William Press and Paul Schechter took the audacious step of counting halos as points of a Gaussian field above that barrier — and, confronted with a factor of 2 that would not go away, simply multiplied by it, guessing rightly that the missing mass hid in the cloud-in-cloud problem. For seventeen years that fudge was an embarrassment, until Bond, Cole, Efstathiou, and Kaiser recast collapse in 1991 as an excursion-set random walk against an absorbing barrier and showed the factor was exact for a sharp-k filter. Navarro, Frenk, and White then found in 1996–97 that halos share a universal two-parameter density profile regardless of mass or cosmology; Sheth and Tormen extended the barrier to ellipsoidal collapse in 1999; and by 2008 Jeremy Tinker and collaborators had calibrated the mass function on N-body simulations to percent accuracy — turning an inspired guess into a precision tool that today weighs the universe by counting its clusters.
Going further
Texts. Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 7 (gravitational collapse, mass function, bias) · Dodelson, Modern Cosmology, 2nd ed., Ch. 12 (nonlinear structure; Press–Schechter) · Ryden, Introduction to Cosmology, 2nd ed., Ch. 12 (structure formation) · Binney & Tremaine, Galactic Dynamics, 2nd ed., §9 (violent relaxation, virial equilibrium) · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., cosmology chapters. The collapse constants δ_c = (3/20)(12π)^(2/3) ≈ 1.686 and Δ_vir = 18π² ≈ 178 are exact EdS spherical-collapse results, reproducible from the cycloid with t_coll = 2 t_ta and 1 + δ_ta = 9π²/16 = 5.55. The PS multiplicity is √(2/π) ν e^(−ν²/2); the Sheth–Tormen parameters are a ≈ 0.707, p ≈ 0.3, A ≈ 0.322; and Δ_vir ≈ 337 (rel. mean matter) comes from the Bryan & Norman fit. The peak-height and multiplicity numbers used throughout — ν = 0.84 at 10¹² M⊙, ν = 3.37 at 10¹⁵ M⊙, ST/PS ≈ 2.2, b ≈ 7.1 at the cluster scale — are computed from these relations; the σ(M) values are representative ΛCDM figures.
Prerequisites. Density Perturbations and the Density Contrast · Linear Growth of Structure · Jeans Instability and Cosmological Perturbation Regimes · Transfer Function and the Matter Power Spectrum · The Primordial Power Spectrum · Cold Dark Matter. Leads to. The ΛCDM Concordance Model · Dark Energy and the Cosmological Constant · Baryon Acoustic Oscillations.
Ask the tutor: Ask me to walk through the spherical-collapse cycloid that gives δ_c = 1.686, or to show how a change in σ₈ moves the cluster counts that weigh the universe.
Why is δ_c ≈ 1.686 a linear quantity if collapse is violently nonlinear?
Because it is deliberately a linear-theory bookmark, not a real density. A collapsing sphere reaches a formally infinite contrast at zero radius, and even at turnaround it is already 5.5× the background — nothing linear about it. But the statistics we have describe only the initial, linear field. So the trick is to ask: if we had evolved the initial contrast with linear theory (δ ∝ growth factor) and read it off at the collapse epoch, what value would we get? That number is (3/20)(12π)^(2/3) ≈ 1.686 — a fictitious linear label attached to a nonlinear event, and exactly what lets us convert Gaussian field statistics into halo counts.
What exactly is the peak height ν, and why does everything depend on it?
The peak height is ν = δ_c/σ(M,z): the ratio of the fixed collapse barrier to the rms amplitude of the smoothed field on mass scale M at redshift z — how many standard deviations a collapsing region of that mass sits above the mean. It matters because the mass function's dependence on mass, redshift, and cosmology all funnel through this one number: the multiplicity ν f(ν) is nearly a universal function of ν alone. Large ν (small σ, big M) means a rare high peak and exponential suppression; ν ≈ 1 marks the characteristic mass M∗ dividing the power-law halos from the exponential ones; ν ≪ 1 is the crowded low-mass regime.
Why must every quoted halo mass function specify a mass definition?
Because "the edge of a halo" is a convention, not a physical surface. Common choices are M₂₀₀c (within the radius enclosing 200× the critical density), M₂₀₀m (200× the mean matter density), M_vir (using Δ_vir from spherical collapse), and friends-of-friends masses. Because NFW halos have extended outskirts, these can differ by tens of percent for the same object — and Δ_vir itself is 178 in EdS but ≈ 337 relative to mean matter in ΛCDM. The multiplicity f(ν) is only approximately universal within a fixed definition, so calibrated fits like Tinker's are tabulated per overdensity. Quoting an abundance without its mass definition is like quoting a length without a unit.
How does the mass function actually probe dark energy?
Through growth. Dark energy governs the expansion history and hence the linear growth factor D(z), which sets σ(M,z) = σ(M,0) D(z)/D(0). Because cluster abundance depends on σ through the exponential e^(−δ_c²/2σ²), a small change in how fast structure grows produces a large change in how many clusters exist at a given redshift. By counting clusters as a function of mass and redshift, surveys measure the growth history directly and can detect whether dark energy is a cosmological constant or something evolving — complementing geometric probes like baryon acoustic oscillations that measure distances rather than growth.
If small halos are so abundant, why doesn't the low-mass end diverge?
In pure cold dark matter theory the abundance really does keep rising toward low mass as a near power law dn/dM ∝ M⁻², with no intrinsic cutoff in the idealized calculation. Reality supplies two. First, the dark-matter particle has a free-streaming scale: if it was ever relativistic (warm dark matter) it erased power below some mass; for cold dark matter this cutoff can be as low as Earth-mass, but it exists. Second, baryonic physics — photoheating during reionization, supernova feedback — prevents gas from cooling into the shallowest wells, so the tiniest halos stay dark even if they form. The low-mass slope is therefore a sensitive probe of the dark-matter particle's nature.
Where does the universal NFW profile come from, and does the theory predict it?
The NFW profile — ρ ∝ 1/[(r/r_s)(1 + r/r_s)²], with an inner r⁻¹ cusp and an outer r⁻³ tail — was an empirical discovery from N-body simulations: halos across an enormous range of mass and cosmology relax to nearly the same shape, differing only in one concentration parameter. It is not derived from the mass-function machinery, which is purely statistical and says nothing about internal structure. The physics behind it is violent relaxation plus the hierarchical accretion history, which is why concentration correlates with formation time: early-forming halos, built when the universe was denser, end up more centrally concentrated. A full first-principles derivation of the exact NFW shape remains an active research question.