---
title: Radial Galaxy Number Density Profile
url: https://www.emergentmind.com/topics/radial-galaxy-number-density-profile
type: topic
---

# Radial Galaxy Number Density Profile

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 [1207.2542][1302.2912][1201.5491].

## 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 [1201.5491][1201.1296][2105.05914].

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
$$
\rho_g(r) \propto r^2 \exp \left[-\beta (c r)^\alpha\right]
$$
with mass-dependent parameters $\alpha$ and $\beta$, providing better agreement with both the halo occupation distribution and the two-point correlation function on small scales [2308.03298].

Complex models, such as the Diemer & Kravtsov (DK14) form, introduce an explicit transition (“splashback”) feature:
$$
\rho(r) = \rho_{\rm inner}(r) f_{\rm trans}(r) + \rho_{\rm outer}(r)
$$
where $f_{\rm trans}(r)$ steepens the profile near the halo edge, typically at $r_{\rm sp} \sim 2\,R_{200}$ [2105.05914][2311.10854].

## 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:
$$
n(R) =
\begin{cases}
n_0 (R/R_0)^{-\alpha_1} & (R \leq r_b)\\
n_0 (r_b/R_0)^{-\alpha_1} (R/r_b)^{-\alpha_2} & (R > r_b)
\end{cases}
$$
with $n_0 = 0.35 \pm 0.18\,\mathrm{kpc}^{-3}$, break radius $r_b = 21 \pm 2$ kpc, inner slope $\alpha_1 = 2.8 \pm 0.4$, and outer slope $\alpha_2 = 4.8 \pm 0.4$ [1802.00570]. The thick disk is mapped in $(R_\mathrm{GC}, |Z|)$ bins without a global power-law fit.

In the innermost regions, the projected density of RRab stars follows three distinct power laws: $-0.94 \pm 0.051$ (cusp, $0 < R < 2.2^\circ$), $-1.50 \pm 0.019$ (bulge, $2.2^\circ < R < 8.0^\circ$), and $-2.43 \pm 0.043$ (halo transition, $8.0^\circ < R < 30.0^\circ$), with half-population radii demonstrating that RRab stars are more extended than red-clump giants but more centrally peaked within $\sim 0.3$ kpc [2010.06603].

## 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 $c_{\rm gal} \simeq 2.6$—nearly independent of mass and roughly half that of the dark matter profile $c_{\rm dm} \sim 4-6$ [1201.5491]. Both projected and 3D number density profiles are strikingly self-similar across a wide range of halo masses and redshifts; deviations occur at $R \lesssim 0.3 R_{500}$ due to tidal stripping and central merging.

The splashback radius $R_{\rm sp}$, a physical boundary set by the first apocenter of recently accreted galaxies, manifests as a steepening (minimum in $d\ln n / d\ln r$) at $R_{\rm sp}\sim(1-2) R_{200}$ [2311.10854][2105.05914][1509.07506]. 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 $R_{\rm sp}$, 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 $r_{\rm sp}\sim2.1h^{-1}$ Mpc and slopes of $-3.4$ near $r_{\rm sp}$, steeper than classical NFW ($-2.7$) [2105.05914].

## 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 $c$ increasing as a function of decreasing satellite luminosity and minimal dependence on host properties for typical masses [1201.1296]:
- Bright satellites: $c \sim 1-1.5$
- Faint satellites: $c \sim 4-5$

Deviations from NFW arise at small radii for the faintest satellites, with up to 20% excess at $r/r_{200}<0.1$, possibly due to unresolved ultra-compact companions or intrinsic profile steepening.

Mass segregation is weak: outside the central baryon-dominated region ($R > 25$ kpc), satellites of different luminosities have nearly identical radial profiles. Only for $R < 25$ kpc does a mild overabundance of bright satellites appear, likely reflecting dynamical friction and tidal survivorship effects [1201.5114].

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 $z=1.6$ show remarkably little evolution in the radial projected number density profile around massive galaxies. From $z=1.6$ to $z=0$, both the slope ($\alpha\simeq-1.8\pm0.3$) and normalization of the projected profile remain consistent, indicating a steady-state balance between merger-driven depletion and accretion-driven replenishment of satellites [1302.2912].

On cosmological scales, the radial galaxy number density as a function of redshift, $n(z)$, must account for the luminosity function, flux limit, survey geometry, and cosmological distance measures. In a flux-limited survey, observed profiles $n_i(z)$ deviate from the constant comoving density predicted for a complete, volume-limited sample, and distortions grow with redshift [1207.2542]. The relevant expressions are:
$$
n_i(z) = \phi^* \Gamma(\alpha+1,\,4\pi D_L^2(z) f_\mathrm{lim} / L^*)
$$
where $\phi^*$ and $L^*$ are Schechter function parameters, and $D_L(z)$ 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 $\sim0.3\,R_{500}$ [1201.5491].
- 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 [2311.10854][2105.05914][1509.07506].
- Baryonic dominance at small radii: The excess of satellites over NFW in central $\lesssim25$ kpc is consistent with the stellar mass profile and can be modeled with Sèrsic components, impacting the local dark-to-baryonic mass ratio [1201.5114].
- 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 [2010.06603].
- 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 [2308.03298][2311.10854]. Observational confirmation of steepening (splashback/density jumps) is now routine, with position and sharpness tied to accretion rate and cluster mass.

## References

| Phenomenon / Measurement    | Key Paper(s)        | arXiv ID      |
|----------------------------|---------------------|---------------|
| Broken power-law halo (MW) | Ablimit & Zhao      | 1802.00570    |
| NFW profile in clusters    | Budzynski et al.    | 1201.5491     |
| New exponential profile    | Qin et al.          | 2308.03298    |
| Satellite systematics      | Guo et al.          | 1201.1296     |
| Splashback feature         | Diemer & Kravtsov, Shin et al. | 2105.05914, 2311.10854 |
| RR Lyrae stellar bulge/halo| Navarro et al.      | 2010.06603    |
| Luminosity function/statistics | Iribarrem et al.| 1207.2542     |
| Profile redshift evolution | Tal et al., van Dokkum et al. | 1302.2912, 1201.5114 |

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.

Source: https://www.emergentmind.com/topics/radial-galaxy-number-density-profile