---
title: Stellar Mass–Halo Mass Relation
url: https://www.emergentmind.com/topics/stellar-mass-halo-mass-relation
type: topic
---

# Stellar Mass–Halo Mass Relation

The stellar mass–halo mass relation, usually abbreviated SHMR or SMHM, is the empirical and theoretical relation linking a galaxy’s stellar mass \(M_\star\) to the mass of its host dark matter halo, \(M_{\rm halo}\). In its most generic form it is written as \(M_\star = f(M_{\rm halo})\), or as a stellar mass fraction \(M_\star/M_{\rm halo}\); some studies further normalize by the cosmic baryon fraction and define \(f_\star \equiv (M_\star/M_{\rm h})/(\Omega_b/\Omega_m)\) as a time-averaged star-formation efficiency [1311.5492, 2102.11282]. The relation is a central ingredient of abundance matching, halo occupation modeling, semi-analytic models, and hydrodynamical simulations, but it is not operationally unique: published work variously uses \(M_{200}\), \(M_{500}\), or \(M_{\rm 200c}^{\rm DMO}\), and the stellar component may refer to a central galaxy, a brightest cluster galaxy (BCG), or the total stellar mass of all member galaxies within \(r_{200}\) [2511.20165, 1501.01966, 1401.7329].

## 1. Definitions and parameterizations

A common phenomenological form is the double power law adopted in semi-analytic and empirical work,
\[
\frac{M_\star}{M_{\mathrm{H}}} = 2N \left[ \left(\frac{M_{\mathrm{H}}}{M_{1}}\right)^{-\beta} + \left(\frac{M_{\mathrm{H}}}{M_{1}}\right)^{\gamma} \right]^{-1},
\]
where \(N\) is the normalization, \(M_1\) is the characteristic halo mass, \(\beta\) governs the low-mass branch, and \(\gamma\) the high-mass branch [1510.08463]. In the GALFORM reference model, a single non-evolving fit over \(0<z<4\) gives \(N \approx 0.012\), \(\log_{10}(M_1/M_\odot)\approx 11.75\), \(\beta \approx 1.41\), and \(\gamma \approx 0.59\), implying \(M_\star \propto M_{\mathrm{H}}^{2.4}\) at low masses and \(M_\star \propto M_{\mathrm{H}}^{0.41}\) at high masses [1510.08463].

Other studies adopt regime-specific parameterizations. Local Group dwarfs over \(10^7 < M_\star/M_\odot \le 10^8\) were fitted with a single power law \(M_\star \propto M_{\rm halo}^{\alpha}\) with \(\alpha = 3.1\), steepening to \(\alpha = 3.5\) if the completeness limit is pushed to \(5\times10^6\,M_\odot\) [1311.5492]. At the cluster scale, BCG analyses often fit log-linear relations such as \(\log_{10} M_{*,\mathrm{BCG}} = m \log_{10} M_{200} + c\), emphasizing the sub-linear growth of the central galaxy relative to the halo [1908.01559].

The relation is also definition-dependent. In the COLIBRE analysis, the fiducial \(z=0\) SMHM is the median relation between the stellar mass of the central within a 50 physical kpc aperture and the halo mass \(M_{\rm 200c}^{\rm DMO}\) of the matched dark-matter-only halo [2511.20165]. In low-mass X-ray groups, by contrast, \(M_{\rm stars}\) denotes the aggregate stellar mass of all group members within \(r_{200}\), corrected to a common lower stellar-mass limit of \(10^9\,M_\odot\) [1501.01966]. At the high-mass end, some cluster studies include diffuse light smoothly connected to the BCG profile as part of \(M_{*,\rm BCG}\), while others explicitly exclude most intracluster light; this distinction materially affects the normalization of the inferred relation [1401.7329, 2602.10193].

## 2. Inference methods

The dominant empirical strategy is abundance matching: one rank-orders galaxies by stellar mass and halos by mass, then imposes a one-to-one monotonic mapping between the cumulative stellar mass function and the halo mass function [1311.5492, 1411.2597]. Parametric abundance matching was used in COSMOS-UltraVISTA to derive the SHMR from \(z\sim0.2\) to \(z\sim5\), exploiting a large area and highly complete stellar masses in the COSMOS field [1810.10557]. At the dwarf scale, Brook et al. applied the same logic inside a Local Group-sized volume, matching the observed Local Group stellar mass function within 1.8 Mpc to a simulated halo mass function that is well described by a single power law [1311.5492].

Clustering-based methods infer halo masses statistically through halo occupation modeling. In VUDS, the projected two-point function \(w_p(r_p)\) for \(\sim 3000\) galaxies at \(2<z<5\) was fitted with a three-parameter HOD model,
\[
\langle N_g|M\rangle =
\begin{cases}
1 + \left(M/M_1\right)^\alpha, & M > M_{\min} \\
0, & \text{otherwise,}
\end{cases}
\]
so that \(M_{\min}\) provides the characteristic halo mass for a stellar-mass-threshold sample [1412.5584]. BOSS clustering has also been used to constrain the scatter of the high-mass SMHM relation by comparing the large-scale bias \(b(M_{\rm gal})\) of samples ranked by stellar mass to abundance-matching models with different \(\sigma_{\log M}\) [1607.04678].

A second class of methods infers halo masses dynamically for individual systems. For massive early-type galaxies, halo masses have been derived from the phase-space distribution of globular cluster systems using an action-angle distribution function embedded in an NFW potential, with a prior on concentration from the \(\Lambda\)CDM concentration–mass relation [2102.11282]. For isolated field dwarfs, HI rotation curves were modeled with the coreNFW profile, which accounts for feedback-driven cusp–core transformations and yields direct estimates of \(M_{200}\) over \(5\times10^5 \lesssim M_\star/M_\odot \lesssim 10^8\) [1607.03127].

At group and cluster masses, X-ray methods are especially important. Halo masses may be inferred from \(L_X\) scaling relations [1501.01966, 1908.01559], from X-ray temperatures [1401.7329], or, more recently, from stacked eROSITA spectra whose average gas temperatures are converted to \(M_{500}\) and then to \(M_{200}\) via the \(M\)–\(T_X\) relation [2602.10193]. These X-ray-based measurements are complementary to weak lensing, satellite kinematics, and group catalogs.

## 3. Mean shape and redshift evolution

Across many determinations, the stellar fraction rises with halo mass at low masses, peaks near a characteristic halo mass of order \(10^{12}\,M_\odot\), and declines again toward group and cluster scales. In COSMOS-UltraVISTA, the ratio of stellar-to-halo mass content peaks at \(M_{\rm h}=10^{12}M_\odot\) at \(z\sim0.2\), increases to \(M_{\rm h}=10^{12.5}M_\odot\) at \(z\sim2.3\), and then remains flat up to \(z=4\) [1810.10557]. For early-type galaxies inferred from globular-cluster dynamics, the baryon-fraction-normalized stellar fraction peaks at \(M_\star\sim5\times10^{10}M_\odot\) and \(M_{\rm h}\sim10^{12}M_\odot\), then declines at higher masses [2102.11282]. An X-ray-based SHMR for central galaxies across \(10^{12}\)–\(10^{15}\,M_\odot\) likewise finds the relation peaking near \(M_{\rm halo}\sim10^{12}\,M_\odot\), followed by a declining central stellar fraction at higher masses [2602.10193].

Direct high-redshift constraints broadly support the same picture. In VUDS, galaxies with \(M_\star=1.3\times10^9\) to \(7.4\times10^9\,M_\odot\) occupy halos with \(M_h=1.3\times10^{11}\) to \(3\times10^{11}\,M_\odot\), implying SMHR values from \(1\%\) to \(2.5\%\); the corresponding integrated star-formation efficiency is \(6\)–\(9\%\) for lower-mass galaxies and \(16\%\) for galaxies with median stellar mass \(\sim 7\times10^9\,M_\odot\) at \(z\sim3\) [1412.5584]. In the GALFORM semi-analytic model, the standard implementation of supernova feedback and gas reincorporation predicts only weak evolution of the median SHM relation over \(0<z<4\), with nearly constant \(N\), \(M_1\), and low-mass slope [1510.08463].

At the high-mass end, BCG-focused studies find similarly slow evolution. For X-ray-selected groups and clusters over \(0.1 \le z \le 0.65\), the central-galaxy stellar mass–halo mass relation is a simple power law in log space, and the low- and high-redshift fits agree within their uncertainties; no notable redshift evolution is detected since \(z\sim0.65\) [1908.01559]. A plausible implication is that, in this regime, central stellar mass growth broadly tracks halo growth, even though the stellar mass fraction declines with increasing halo mass.

## 4. Mass-regime dependence

The dwarf regime remains the sharpest test of the low-mass branch. For Local Group galaxies with \(10^7 \lesssim M_\star/M_\odot \lesssim 10^8\), Brook et al. find a steep power-law SHMR with slope \(3.1\), significantly steeper than most abundance-matching extrapolations based on large surveys [1311.5492]. For isolated field dwarfs, rotation-curve modeling with coreNFW yields a monotonic \(M_\star\)–\(M_{200}\) relation with little scatter, and abundance matching based on the SDSS field stellar mass function agrees with the directly inferred dwarf relation down to \(M_{200}\sim 5\times10^9\,M_\odot\), or to \(5\times10^8\,M_\odot\) if the SDSS stellar mass function is extrapolated as a power law below \(M_\star\sim10^7\,M_\odot\) [1607.03127].

In low-mass X-ray groups at \(0.5<z<1\), the relevant observable is often the total stellar mass in member galaxies rather than the central galaxy alone. Patel et al. measure
\[
\log M_{\rm stars} = (0.47 \pm 1.33) + (0.84 \pm 0.10)\,\log M_{200}
\]
for groups, with an observed scatter of \(0.25\) dex, and find that stars comprise \(\sim 3\)–\(4\%\) of the total mass for halos with \(10^{12.8}<M_{200}/M_\odot<10^{13.5}\) [1501.01966]. This sub-unity slope means that \(M_{\rm stars}/M_{200}\) decreases toward higher halo masses.

At cluster scales, the central-galaxy relation is distinctly shallow. For BCGs, the relation \(\log_{10}(M_{*,\mathrm{BCG}}/M_\odot) = (0.41\pm0.04)\,x + (5.59\pm0.54)\) at \(0.1\le z\le0.3\) and \((0.31\pm0.02)\,x + (7.00\pm0.37)\) at \(0.3<z\le0.65\), with \(x=\log_{10}(M_{200}/M_\odot)\), shows that BCG stellar mass increases sub-linearly with halo mass [1908.01559]. A deeper photometric re-analysis of nearby clusters finds \(M_{*,\rm BCG}\propto M_{500}^{0.33\pm0.11}\), \(M_{*,\rm sat}\propto M_{500}^{0.75\pm0.085}\), and \(M_{*,\rm tot}\propto M_{500}^{0.59\pm0.08}\), again implying a declining stellar fraction with increasing cluster mass [1401.7329]. This regime is also particularly sensitive to aperture and surface-brightness systematics: older abundance-matching calibrations based on standard SDSS photometry were argued to underestimate BCG stellar masses by factors of \(\sim 2\)–4 [1401.7329].

## 5. Scatter and secondary dependencies

A recurrent result is that the scatter in stellar mass at fixed halo mass is small, typically \(\lesssim 0.2\) dex, but the physical origin of that scatter is not unique. From BOSS clustering of galaxies with \(\log M_\star \gtrsim 11.4\), the preferred scatter is \(\sigma_{\log M}=0.18^{+0.01}_{-0.02}\), including measurement error, and repeated spectra imply an upper limit of \(0.16\) dex on the intrinsic scatter [1607.04678]. In a merger-tree analysis anchored to the Behroozi et al. relation, hierarchical assembly alone produces \(\approx 0.16\) dex scatter at \(M_{\rm vir}>10^{14}\,M_\odot\), while adding an intrinsic in-situ scatter of \(0.2\) dex yields an approximately constant \(\sim 0.2\) dex scatter from \(10^{12}\) to \(10^{14.75}\,M_\odot\) at \(z=0\); this suggests that the observed flatness of the scatter with halo mass is largely a coincidence of two different growth channels [1602.01099].

Secondary halo properties measurably structure that scatter. In COLIBRE, at fixed \(M_{\rm 200c}^{\rm DMO}\), stellar mass correlates positively with halo concentration, with Spearman \(\mathcal R_{\rm s}\approx 0.50\) near \(10^{11}\,M_\odot\) and \(\approx 0.46\) near \(10^{11.5}\,M_\odot\); the correlation weakens but remains positive above \(10^{12}\,M_\odot\) [2511.20165]. Concentration also correlates with stellar age, but stellar age itself correlates only weakly with stellar mass at fixed halo mass, whereas stellar metallicity correlates strongly with both concentration and stellar mass; moreover, the concentration–metallicity correlation persists at fixed stellar mass and halo mass with mean \(\mathcal R_{\rm s}\approx 0.34\). This supports a potential-depth interpretation in which deeper halos suppress feedback-driven outflows and retain more baryons and metals [2511.20165].

Morphology, color, and size also act as second parameters. Dynamical modeling of early types and rotation-curve modeling of late types indicate that the SHMR is not universal across morphology: early types follow the familiar peaked relation, while late types show a monotonically rising \(f_\star\) over the mass range probed, with the two branches differing by a factor of \(\sim 7\) in \(f_\star\) at \(M_\star\sim10^{11}M_\odot\) [2102.11282]. At fixed very high stellar mass \(10^{11.3}<M_\star/M_\odot<10^{11.7}\) and \(z_s\sim0.6\), PAC measurements show that the most compact galaxies with Sérsic index \(n>6\) have halo masses around \(5.5\) times larger than disk-like systems with \(n<2\); red galaxies inhabit halos \(2.6\) times more massive than blue galaxies, and large galaxies halos \(2.3\) times more massive than small galaxies [2110.05760]. Environment matters as well: in hydrodynamical simulations, central galaxies in large-scale overdense regions have larger \(M_\star/M_{\rm halo}\) than counterparts in underdense regions at the same halo mass, and this persists even when the local density within 2 Mpc at \(z=0\) is matched [1509.05039].

## 6. Physical interpretation, controversies, and open directions

The low-mass branch is widely interpreted as feedback- and reionization-regulated. In field dwarfs, the decline of \(M_\star/M_{200}\) toward low halo mass is consistent with efficient supernova-driven outflows, with sufficient supernova energy to create dark-matter cores down to \(M_{200}\sim 5\times10^8\,M_\odot\) if star formation has persisted for long enough [1607.03127]. In fully cosmological SPH simulations with metal-line cooling, dust and self-shielding, H\(_2\)-based star formation, and supernova-driven outflows, the present-day SHM relation agrees well with abundance matching over \(M_\star \sim 2.2\times10^3\) to \(4.5\times10^{10}\,M_\odot\); the agreement improves once stellar masses are measured from synthetic photometry and halo masses are compared to dark-matter-only runs, because photometric stellar masses can underestimate true stellar masses and dark-matter-only runs overestimate halo masses when baryon loss is ignored [1209.1389].

The high-mass decline is usually attributed to AGN feedback, long cooling times, and ex-situ assembly. In the eROSITA-based SHMR, the decline of \(M_{\star,\mathrm{BCG}}/M_{200}\) at group and cluster scales is interpreted as the result of AGN feedback, reduced cooling efficiency, and the increasing dominance of ex-situ assembly while halos continue to grow through mergers and accretion [2602.10193]. GALFORM makes this logic explicit: weak evolution of the median SHM relation over \(0<z<4\) follows if gas reincorporation efficiency is effectively constant and the average efficiency of supernova feedback is approximately constant at fixed halo mass, while AGN feedback introduces a break whose location evolves only modestly [1510.08463]. This suggests that the near-stationarity of the median relation is a constraint on gas cycling as much as on star formation itself.

Several controversies remain methodological rather than purely physical. The first is the treatment of extended light in massive galaxies: improved surface-brightness photometry and Sérsic-based fits can raise BCG stellar masses enough to steepen the inferred high-mass SMHM slope and reduce the apparent tension between direct cluster measurements and abundance matching [1401.7329, 1411.2597]. The second is environmental non-universality: abundance matching based on the field stellar mass function produces a dwarf-scale SHMR consistent with isolated field rotation curves, whereas using the stellar mass function of nearby groups recovers the familiar missing-satellites and too-big-to-fail problems; this indicates that those problems are confined to group environments and arise from galaxy-formation physics rather than exotic cosmology [1607.03127]. The third is that morphology-blind or environment-blind parameterizations can absorb real secondary structure into an apparently larger intrinsic scatter; one estimate suggests that ignoring the late-type branch could inflate the inferred intrinsic scatter by \(\approx 0.08\) dex at \(M_\star=10^{11}M_\odot\) [2102.11282].

Future progress in the literature represented here is correspondingly twofold. Deeper Local Group surveys are expected to push completeness to \(10^{6.5}M_\odot\) or below and tighten the low-mass slope [1311.5492]; larger samples of massive galaxies are needed to test the redshift dependence of the characteristic halo mass and galaxy-formation efficiency at \(z\gtrsim2\) [1810.10557]. On the group and cluster side, deeper X-ray data and larger stacked samples should enable SHMR measurements conditioned on morphology, star-formation rate, concentration, or luminosity gap, directly probing the second-parameter structure that is already evident in simulations and high-mass observations [2602.10193].

Source: https://www.emergentmind.com/topics/stellar-mass-halo-mass-relation