Radial Galaxy Number Density Profile

• Radial galaxy number density profiles describe how galaxy counts change with distance from a central reference, crucial for interpreting galaxy evolution and dark matter distributions.
• Key analytic models like the NFW, exponential, and DK14 profiles capture essential features such as inner depletion and splashback steepening, enhancing our understanding of halo dynamics.
• Robust observational methods addressing survey completeness, flux limits, and background subtraction are vital for accurately mapping profiles in clusters, groups, and the Milky Way.

A radial galaxy number density profile describes how the number density of galaxies or specific galaxy tracers varies as a function of distance from a reference point, such as the center of a galaxy, galaxy group, cluster, or the Milky Way itself. These profiles are foundational both for interpreting galaxy formation and evolution, and for linking the observed distribution of galaxies to the three-dimensional structure and dynamics of dark matter halos.

1. Definitions and Observational Framework

The galaxy number density profile can be defined in both three-dimensional (3D) and two-dimensional (projected) forms. The 3D (spherical) number density, $n(r)$, describes the number of galaxies per unit volume at distance $r$ from the center:

$n(r) = \frac{dN}{dV}$

The projected or surface density, $\Sigma(R)$, gives the number per projected area at projected distance $R$:

$$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$

Survey completeness, flux limits, background/foreground subtraction, and the choice of radial bins (logarithmic, linear) all impact the estimation of these profiles.

In cosmology, the profile is often measured with respect to the center of a halo (for clusters), central/satellite galaxies, or the Galactic center (for the Milky Way). The choice of distance measure (projected vs. 3D, and particular cosmological distances) and the normalization—such as per primary galaxy or per cluster—must be made explicit (Iribarrem et al., 2012, Tal et al., 2013, Budzynski et al., 2012).

2. Canonical Parametric Models

The most widely-adopted analytic form for the radial galaxy number density profile in halos is the Navarro-Frenk-White (NFW) profile, originally developed for dark matter:

$n(r) = \frac{n_s}{(r/r_s)(1 + r/r_s)^2}$

where $r_s$ is the scale radius and $n_s$ is a normalization constant. The concentration $c = R_{200}/r_s$ describes the relative sizes of the halo and its inner region. This profile fits both three-dimensional and properly projected observations in groups and clusters (Budzynski et al., 2012, Guo et al., 2012, Shin et al., 2021).

Recent developments show that the NFW profile systematically overpredicts galaxy counts near the core and underpredicts counts in the outskirts. The new two-parameter exponential profile introduced by Qin et al. takes the form

$r$0

with mass-dependent parameters $r$1 and $r$2, providing better agreement with both the halo occupation distribution and the two-point correlation function on small scales (Qin et al., 2023).

Complex models, such as the Diemer & Kravtsov (DK14) form, introduce an explicit transition (“splashback”) feature:

$r$3

where $r$4 steepens the profile near the halo edge, typically at $r$5 (Shin et al., 2021, Pizzardo et al., 2023).

3. Milky Way and Local Group

For the Milky Way, the radial number density profile of RR Lyrae stars has been extensively used to trace the spatial distribution of old stars in both the disk and halo. Ablimit & Zhao find that the stellar halo is best described by a broken power-law:

$r$6

with $r$7, break radius $r$8 kpc, inner slope $r$9, and outer slope $n(r) = \frac{dN}{dV}$0 (Ablimit et al., 2018). The thick disk is mapped in $n(r) = \frac{dN}{dV}$1 bins without a global power-law fit.

In the innermost regions, the projected density of RRab stars follows three distinct power laws: $n(r) = \frac{dN}{dV}$2 (cusp, $n(r) = \frac{dN}{dV}$3), $n(r) = \frac{dN}{dV}$4 (bulge, $n(r) = \frac{dN}{dV}$5), and $n(r) = \frac{dN}{dV}$6 (halo transition, $n(r) = \frac{dN}{dV}$7), with half-population radii demonstrating that RRab stars are more extended than red-clump giants but more centrally peaked within $n(r) = \frac{dN}{dV}$8 kpc (Navarro et al., 2020).

4. Galaxy Groups and Clusters: NFW and the Splashback Feature

In groups and clusters, the NFW projected profile remains empirically successful. Budzynski et al. demonstrate that galaxy concentration is $n(r) = \frac{dN}{dV}$9—nearly independent of mass and roughly half that of the dark matter profile $\Sigma(R)$0 (Budzynski et al., 2012). Both projected and 3D number density profiles are strikingly self-similar across a wide range of halo masses and redshifts; deviations occur at $\Sigma(R)$1 due to tidal stripping and central merging.

The splashback radius $\Sigma(R)$2, a physical boundary set by the first apocenter of recently accreted galaxies, manifests as a steepening (minimum in $\Sigma(R)$3) at $\Sigma(R)$4 (Pizzardo et al., 2023, Shin et al., 2021, Patej et al., 2015). In simulations and in high-S/N weak lensing and galaxy density measurements, this location matches the boundary where radial infall ceases and the inner halo is virialized.

The Diemer & Kravtsov (DK14) and follow-up models include this sharp transition:

• Inner profile: Einasto or NFW-like
• Transition function: controls steepening at $\Sigma(R)$5, dependent on accretion rate
• Outer profile: accounts for correlated large-scale structure

These models reproduce the observed splashback steepening in both galaxy and mass profiles, with best-fit splashback radii $\Sigma(R)$6 Mpc and slopes of $\Sigma(R)$7 near $\Sigma(R)$8, steeper than classical NFW ($\Sigma(R)$9) (Shin et al., 2021).

5. Satellite Systems and Mass Segregation

Statistical studies of satellite systems around isolated galaxies (primaries) show that the projected number density of satellites is well fit by the NFW form, with concentration $R$0 increasing as a function of decreasing satellite luminosity and minimal dependence on host properties for typical masses (Guo et al., 2012):

• Bright satellites: $R$1
• Faint satellites: $R$2

Deviations from NFW arise at small radii for the faintest satellites, with up to 20% excess at $R$3, possibly due to unresolved ultra-compact companions or intrinsic profile steepening.

Mass segregation is weak: outside the central baryon-dominated region ($R$4 kpc), satellites of different luminosities have nearly identical radial profiles. Only for $R$5 kpc does a mild overabundance of bright satellites appear, likely reflecting dynamical friction and tidal survivorship effects (Tal et al., 2012).

Color and morphology of both satellites and primaries modulate the amplitude and concentration, with red/early-type centrals hosting more and more concentrated satellites compared to blue/late types.

6. Redshift Evolution and Cosmological Considerations

Measurements out to $R$6 show remarkably little evolution in the radial projected number density profile around massive galaxies. From $R$7 to $R$8, both the slope ($R$9) and normalization of the projected profile remain consistent, indicating a steady-state balance between merger-driven depletion and accretion-driven replenishment of satellites (Tal et al., 2013).

On cosmological scales, the radial galaxy number density as a function of redshift, $$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$0, must account for the luminosity function, flux limit, survey geometry, and cosmological distance measures. In a flux-limited survey, observed profiles $$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$1 deviate from the constant comoving density predicted for a complete, volume-limited sample, and distortions grow with redshift (Iribarrem et al., 2012). The relevant expressions are:

$$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$2

where $$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$3 and $$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$4 are Schechter function parameters, and $$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$5 is the luminosity distance. Completeness corrections and distance reciprocity must be carefully included.

7. Physical Interpretation and Astrophysical Implications

The radial galaxy number density profile encodes significant information about formation histories, dynamical processes, and the interplay between baryonic and dark matter physics. Key insights include:

• Inner flattening (relative to NFW): Tidal stripping, dynamical friction, and mergers with the central galaxy deplete the satellite population within $$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$6 (Budzynski et al., 2012).
• Outer steepening (splashback): The density jump or splashback feature demarcates the boundary between infalling and virialized material, observed as a sharp steepening in both galaxy and mass profiles (Pizzardo et al., 2023, Shin et al., 2021, Patej et al., 2015).
• Baryonic dominance at small radii: The excess of satellites over NFW in central $$\Sigma(R) = \int_{-\infty}^{+\infty} n \left( \sqrt{R^2 + z^2} \right) dz$$7 kpc is consistent with the stellar mass profile and can be modeled with Sèrsic components, impacting the local dark-to-baryonic mass ratio (Tal et al., 2012).
• Population differences: Old, metal-poor stellar populations—traced by RR Lyrae and globular clusters—have distinct, more extended inner cusps and shallower outer halos compared to more metal-rich populations (Navarro et al., 2020).
• Profile universality: High self-similarity across halo mass and epoch supports universality in halo assembly under ΛCDM.

Recent models, validated with SDSS and simulations (DARKSAGE, IllustrisTNG), enable percent-level precision mocks for next-generation galaxy surveys, essential for constraining clustering, lensing, and cosmological parameters (Qin et al., 2023, Pizzardo et al., 2023). Observational confirmation of steepening (splashback/density jumps) is now routine, with position and sharpness tied to accretion rate and cluster mass.

Phenomenon / Measurement	Key Paper(s)	arXiv ID
Broken power-law halo (MW)	Ablimit & Zhao	(Ablimit et al., 2018)
NFW profile in clusters	Budzynski et al.	(Budzynski et al., 2012)
New exponential profile	Qin et al.	(Qin et al., 2023)
Satellite systematics	Guo et al.	(Guo et al., 2012)
Splashback feature	Diemer & Kravtsov, Shin et al.	(Shin et al., 2021, Pizzardo et al., 2023)
RR Lyrae stellar bulge/halo	Navarro et al.	(Navarro et al., 2020)
Luminosity function/statistics	Iribarrem et al.	(Iribarrem et al., 2012)
Profile redshift evolution	Tal et al., van Dokkum et al.	(Tal et al., 2013, Tal et al., 2012)

These studies collectively establish the radial galaxy number density profile as an essential diagnostic for both baryonic and dark matter structure, with robust methodologies for measuring and interpreting its form across all cosmic environments and epochs.