Pith. Sign in

REVIEW

The Halo Occupation Distribution: Towards an Empirical Determination of the Relation Between Galaxies and Mass

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Andreas A. Berlind, David H. Weinberg (The Ohio State University)

Astrophysics · Submitted 1 Sep 2001 · 102 references

Pith's verdict: A solid, durable paper that maps which clustering statistic probes which occupation parameter — but it demonstrates sensitivity, not identifiability, and the abstract's cosmology claim runs ahead of the evidence.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/0109001 (version 2)
Pith IDTTDVQ46S, Pith's permanent record of this paper and its review
CitationAndreas A. Berlind, David H. Weinberg (The Ohio State University). "The Halo Occupation Distribution: Towards an Empirical Determination of the Relation Between Galaxies and Mass." arXiv:astro-ph/0109001 (2001). https://pith.science/paper/TTDVQ46S
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Galaxy clustering can be read backwards to reveal what sits inside dark matter halos.

The reading

The paper asks a deceptively practical question: given only maps of where galaxies are, can we work out how galaxies are distributed among dark matter halos, rather than assuming an answer from a theory of galaxy formation? Its language for that distribution is the halo occupation distribution, or HOD: for a halo of virial mass $M$, the probability $P(N|M)$ that it holds $N$ galaxies, together with rules for where inside the halo those galaxies sit and how fast they move. The authors populate the halos of a $256^3$-particle $N$-body simulation of a flat $0.3$-$0.7$ cold dark matter universe with many such prescriptions and measure how each clustering statistic shifts when one feature of the occupation is changed. They find that the statistics divide the labour: the large-scale amplitude of galaxy clustering depends only on the mean occupation, while the small-scale shape depends on the second factorial moment of the counts, so the observed power-law correlation function is a balance between two regimes rather than a generic outcome. If the argument holds, redshift surveys can yield an empirical occupation distribution for each class of galaxy, which is a target for formation theory and, if cosmological and bias effects really are separable, a sharper test of cosmology.

What carries the argument

The load-bearing device is the one-halo / two-halo decomposition of the correlation function, used as an interpretive scalpel rather than a proof. Equation (8) gives the large-scale bias as an integral of the mean occupation against the halo bias factor; equation (11) gives the one-halo term as an integral of $N(N-1)_M$ against the distribution of pair separations inside halos, $F(r/2R_{vir})$. Because the two integrals weight different moments of $P(N|M)$, changing the width of the count distribution moves the small-scale amplitude while leaving the large-scale amplitude untouched, and changing $M_{min}$ or the slope moves both by different amounts. The authors verify numerically that $F(r/2R_{vir})$ is nearly independent of halo mass and only weakly sensitive to the galaxy profile, which is what licenses reading the small-scale correlation function as a measurement of $N(N-1)_M$ in a localized mass range. The second piece of machinery is the empirical inversion of Section 6: match the cumulative group multiplicity function to the cumulative halo mass function to recover $N_{avg}(M)$ at high mass, invert the one-halo term via equation (30) to recover $N(N-1)_M$ at low mass, and use the mean density, the large-scale bias factor, the void probability function, and the pairwise velocity dispersion as cross-checks.

Extended reading notes

Core claim

The central claim is that galaxy bias, described at the level of virialized halos, is measurable rather than merely parametrizable, and that different clustering statistics carry complementary slices of it. The mechanism the authors isolate is a clean split by scale. Above the virial diameter of the largest halos, all pairs lie in separate halos, and the galaxy correlation function is $b^2$ times the mass correlation function with $b = {n_g}^{-1} dM (dn/dM) N_{avg}(M) b_h(M)$, a weighted average of halo bias that sees only the mean occupation $N_{avg}(M)$. Below that scale the correlation function is dominated by pairs inside a single halo, and its amplitude runs through the second factorial moment $N(N-1)_M$, which is larger for wider count distributions at fixed mean. From this follows the paper's most concrete result: a power-law galaxy correlation function is not a generic prediction of the framework. Reproducing the observed form requires a large fraction of galaxies in low-multiplicity halos, which in practice means sub-Poisson count fluctuations, a minimum halo mass below roughly $3 10^{11} h^{-1} M_$, and a shallow mass slope near $= 0.4-0.6$; Poisson models need implausibly low minimum masses or a broken power-law occupation. The authors argue on this basis that the success of semi-analytic and hydrodynamic models in producing a power-law correlation function is non-trivial, and they sketch a procedure by which the occupation could be extracted from survey data instead of assumed.

Load-bearing premise

The whole framework rests on the assumption that a halo's galaxy content depends only on its virial mass and not on whether the halo sits in a dense or empty region, an assumption the authors defend with simulation evidence but explicitly say deserves more detailed theoretical investigation.

Editorial extensions

If this is right

  • The occupation distribution becomes an observable: group multiplicity matching plus the one-halo inversion recovers $N_{avg}(M)$ and $N(N-1)_M$ separately, so the mean and the width of the halo counts can be measured rather than fitted.
  • Sub-Poisson fluctuations in low-multiplicity halos are required for the observed power-law correlation function, which turns a qualitative feature of semi-analytic models into a quantitative constraint on cooling, feedback, and merger physics.
  • Galaxy-galaxy lensing measures the galaxy-mass cross-correlation function, which depends on $N_{avg}(M)$ even where the galaxy autocorrelation depends on $N(N-1)_M$, so the two together separate the mean occupation from the width of the count distribution.
  • The void probability function constrains the low-mass cutoff $M_{min}$ sharply, while the group multiplicity function constrains the high-mass slope, giving the two ends of the occupation relation independent handles.
  • If cosmological and occupation effects are non-degenerate as argued, an incorrect cosmological model should fail to reproduce the full suite of clustering statistics for any choice of occupation, making galaxy clustering a test of cosmology at nonlinear scales.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The inversion strategy is only as good as the group finder: the paper works in real space with a single friends-of-friends linking length, and the mapping from groups to halos degrades for poor groups, so a practical pipeline will inherit a systematic error budget that the idealized test does not show.
  • Because the one-halo term localizes $N(N-1)_M$ in mass through $R_{vir} r$, the method has a natural ceiling around $5 10^{12} h^{-1} M_$ where the two-halo term takes over; pushing it higher requires either subtracting a modelled two-halo contribution or using a numerical inversion of equation (11), both of which reintroduce cosmological dependence.
  • The claim that cosmology and bias are non-degenerate is argued qualitatively rather than demonstrated; a useful extension would be to run the same occupation-fitting procedure on a simulation with a different $m$ or $8$ and check whether any occupation choice reproduces the fiducial clustering down to the precision of a survey like SDSS.
  • The environmental-independence assumption is the natural weak point to test directly: stacking halos of fixed mass in voids versus clusters and comparing their galaxy content would either validate the framework's foundation or show that the one-halo and two-halo terms cannot be separated as cleanly as the formalism requires.
Original abstract

We investigate galaxy bias in the framework of the ``Halo Occupation Distribution'' (HOD), which defines the bias of a population of galaxies by the conditional probability P(N|M) that a dark matter halo of virial mass M contains N galaxies, together with prescriptions that specify the relative spatial and velocity distributions of galaxies and dark matter within halos. By populating the halos of a cosmological N-body simulation using a variety of HOD models, we examine the sensitivity of different galaxy clustering statistics to properties of the HOD. The galaxy correlation function responds to different aspects of P(N|M) on different scales. Obtaining the observed power-law form of xi(r) requires rather specific combinations of HOD parameters, implying a strong constraint on the physics of galaxy formation; the success of numerical and semi-analytic models in reproducing this form is entirely non-trivial. Other clustering statistics such as the galaxy-mass correlation function, the bispectrum, the void probability function, the pairwise velocity dispersion, and the group multiplicity function are sensitive to different combinations of HOD parameters and thus provide complementary information about galaxy bias. We outline a strategy for determining the HOD empirically from redshift survey data. This method starts from an assumed cosmological model, but we argue that cosmological and HOD parameters will have non-degenerate effects on galaxy clustering, so that a substantially incorrect cosmological model will not reproduce the observations for any choice of HOD. Empirical determinations of the HOD as a function of galaxy type from the 2dF and SDSS redshift surveys will provide a detailed target for theories of galaxy formation, insight into the origin of galaxy properties, and sharper tests of cosmological models.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,609 input tokens.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.42 · human participants 0.01 · controlled comparison 0.04 · randomized assignment 0.01 · states quantitative result 0.08 · compares with prior results 0.22 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

We investigate galaxy bias in the framework of the ``Halo Occupation Distribution'' (HOD), which defines the bias of a population of galaxies by the conditional probability P(N|M) that a dark matter halo of virial mass M contains N galaxies, together with prescriptions that specify the relative spatial and velocity distributions of galaxies and dark matter within halos. By populating the halos of a cosmological N-body simulation using a variety of HOD models, we examine the sensitivity of different galaxy clustering statistics to properties of the HOD. The galaxy correlation function responds to different aspects of P(N|M) on different scales. Obtaining the observed power-law form of xi(r) requires rather specific combinations of HOD parameters, implying a strong constraint on the physics of galaxy formation; the success of numerical and semi-analytic models in reproducing this form is entirely non-trivial. Other clustering statistics such as the galaxy-mass correlation function, the bispectrum, the void probability function, the pairwise velocity dispersion, and the group multiplicity function are sensitive to different combinations of HOD parameters and thus provide complementary information about galaxy bias. We outline a strategy for determining the HOD empirically from redshift survey data. This method starts from an assumed cosmological model, but we argue that cosmological and HOD parameters will have non-degenerate effects on galaxy clustering, so that a substantially incorrect cosmological model will not reproduce the observations for any choice of HOD. Empirical determinations of the HOD as a function of galaxy type from the 2dF and SDSS redshift surveys will provide a detailed target for theories of galaxy formation, insight into the origin of galaxy properties, and sharper tests of cosmological models.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We investigate galaxy bias in the framework of the ``Halo Occupation Distribution'' (HOD), which defines the bias of a population of galaxies by the conditional probability P(N|M) that a dark matter halo of virial mass M contains N galaxies, together with prescriptions that specify the relative spatial and velocity distributions of galaxies and dark matter within halos. background 0.79 0.14 0.05 asserted 0.94
s1 By populating the halos of a cosmological N-body simulation using a variety of HOD models, we examine the sensitivity of different galaxy clustering statistics to properties of the HOD. method 0.91 0.95 0.04 asserted 0.90
s2 The galaxy correlation function responds to different aspects of P(N|M) on different scales. result 0.98 0.68 0.09 asserted 1.00
s3 Obtaining the observed power-law form of xi(r) requires rather specific combinations of HOD parameters, implying a strong constraint on the physics of galaxy formation; the success of numerical and semi-analytic models in reproducing this form is entirely non-trivial. interpretation 0.75 0.66 0.25 asserted 0.99
s4 Other clustering statistics such as the galaxy-mass correlation function, the bispectrum, the void probability function, the pairwise velocity dispersion, and the group multiplicity function are sensitive to different combinations of HOD parameters and thus provide complementary information about galaxy bias. result 0.71 0.57 0.11 asserted 0.99
s5 We outline a strategy for determining the HOD empirically from redshift survey data. contribution 0.51 0.60 0.08 asserted 0.90
s6 This method starts from an assumed cosmological model, but we argue that cosmological and HOD parameters will have non-degenerate effects on galaxy clustering, so that a substantially incorrect cosmological model will not reproduce the observations for any choice of HOD. method 0.52 0.27 0.17 hedged possibility 0.48
s7 Empirical determinations of the HOD as a function of galaxy type from the 2dF and SDSS redshift surveys will provide a detailed target for theories of galaxy formation, insight into the origin of galaxy properties, and sharper tests of cosmological models. future work 0.65 0.10 0.32 asserted 0.50

What kind of paper the classifier reports this is: simulation study 0.49 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. — P(N|M) relation for two particular HOD models. Each point represents the number of galaxies that occupy a single halo in the N-body simulation. Points for halos that contain no galaxies are arbitrarily placed at logN = −0.5. The HOD prescriptions shown have a power-law Navg(M) with α = 0.5 and Mmin = 2.8 × 1011h −1M⊙ (see definitions in § 2), as indicated by the solid lines. The P(N|Navg) distribution are Average (top panel) and Poisson (bottom panel)
Figure 2. — The cumulative fraction of galaxy pairs within halo
Figure 3. — Influence of Mmin on the galaxy correlation function. Curves show galaxy correlation functions for HOD models with a power-law Navg(M), α = 0.5, Average P(N|Navg), and different values of Mmin, which are listed in the legend in units of 1011h −1M⊙. The solid curve shows the mass correlation function, and points show the correlation function measured from the APM galaxy survey (Baugh 1996)
Figure 4. — Influence of α on the galaxy correlation function. Curves show galaxy correlation functions for HOD models with a power-law Navg(M), Mmin = 2.8 × 1011h −1M⊙, Average P(N|Navg), and different values of α, which are listed in the legend
Figure 5. — Influence of P(N|Navg) on the galaxy correlation function. Curves show galaxy correlation functions for HOD models with a power-law Navg(M), Mmin = 2.8 × 1011h −1M⊙, α = 0.5, and different forms of P(N|Navg), which are listed in the legend
Figure 6. — Influence of the galaxy profiles within halos on the ga
Figure 7. — Contributions to the galaxy correlation function fr
Figure 8. — Power-law fits to a sample of HOD model galaxy correla
Figure 9. — Real-space correlation functions of galaxy distri
Figure 10. — Galaxy-mass cross-correlation functions for diffe
Figure 11. — Effective bias functions defined by ratios of correla
Figure 12. — Reduced bispectrum Qeq(k) for equilateral configurations as a function of scale, for different HOD prescriptions. In panels (a)-(c), the dotted curve shows Qeq(k) for a model with power-law Navg(M), Mmin = 2.8 × 1011h −1M⊙, α = 0.5, and Average P(N|Navg). Other curves show Qeq(k) for models with different Mmin (panel a), different α (panel b), or different P(N|Navg) (panel c). Panel (d) sgows models with different P(N|Navg) for α = 0.8. Thick solid curves show Qeq(k) for the unbiased mass distribution
Figure 13. — The reduced bispectrum Q(θ) as a function of triangle shape, where k1 = 2k2 = 1hMpc−1 (k = 2π/λ). The four panels show the same HOD models as in
Figure 14. — Void probability functions (VPF) for different HOD p
Figure 15. — Mean pairwise radial velocity dispersions
Figure 16. — Ratios of the redshift to real space power spectra
Figure 17. — Cumulative group multiplicity functions for differ
Figure 18. — Richness vs. galaxy group virial mass estimates. E
Figure 19. — Same as Fig. 18, except that large points represent
Figure 20. — Recovery of Navg(M) from the group multiplicity function. Panel (a) shows the halo mass function nh(M) derived from the GIF N-body simulation, where halos were identified with a FoF algorithm. Panel (b) shows group multiplicity functions ngrp(> N) for three HOD models, where galaxy groups were identified with the same FoF algorithm. Panel (c) shows the resulting Navg(M) relations from matching nh(M) to ngrp(> N) (points). Also shown, for comparison, are the true Navg(M) functions that were used as inputs for the HOD models (lines)
Figure 21. — Recovery of hN(N − 1)iM from ξ1h(r) using eq. (30). The solid and open circles show the recovered hN(N − 1)iM for Poisson and Average P(N|Navg) HOD models, respectively. Also shown for comparison, are the “true” hN(N − 1)iM relations that were used as inputs for the HOD models (lines)

Figures reproduced from arXiv:astro-ph/0109001 by the authors.

References
How later papers used it

302 later papers cite it

Evidence: substantial. 285 of the citing papers share no author with it, and 50 of those carry a confident label from Pith's citation classifier.

Pith read 39 of the 359 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 12 of 52 papers with this label · 3 by the same authors

By the same authors

Uses its data or values 1 paper
Cites it as support 1 of 2 papers with this label
How the literature uses it

78 arXiv papers cite this in 90 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 41 sentences without a confident label

Uses its method or tool 18 papers
Uses its data or values 1 paper
Cites it as support 1 paper
Cites it as background 25 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.