🔭 astrophysics

Universal Fitting Formulae for the Peak Concentration of Dark Matter Halos

By analyzing extensive N-body simulations across diverse cosmologies, the authors derive a universal, low-scatter fitting formula that predicts the most probable dark matter halo concentration based on a revised peak height parameter, effectively capturing the relationship's dependence on mass, redshift, and cosmological parameters.

Original authors: Dao-zhou Wang, Weipeng Lin, Tian-Cheng Luan

Published 2026-07-21

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Dao-zhou Wang, Weipeng Lin, Tian-Cheng Luan

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, invisible ocean made of dark matter. We can't see this ocean directly, but we know it's there because it acts like a cosmic glue, holding together the islands of normal matter we call galaxies. Just like how a whirlpool in a bathtub has a specific shape and density, these dark matter "islands" (called halos) have their own internal structure. Scientists have long been trying to figure out exactly how "squished" or "concentrated" these halos are. Think of concentration like the difference between a fluffy cloud and a dense marble; both are made of the same stuff, but one is spread out while the other is packed tight.

Why does this matter? Because the shape of these invisible halos dictates how galaxies form and evolve. If we can predict how concentrated a halo is, we can better understand why some galaxies are bright and busy while others are quiet and dark. For decades, scientists have tried to find a simple rule—a "universal formula"—that links a halo's mass and age to its concentration. It's like trying to find a single recipe that predicts the perfect crust for every pizza, no matter the size or the oven. But the universe is tricky; different cosmic environments and different times in history seemed to break the rules, leaving scientists with a messy pile of data that didn't quite fit together.

This paper, written by Dao-zhou Wang, Weipeng Lin, and Tian-Cheng Luan, dives into a massive collection of computer simulations to solve this puzzle. The authors didn't just look at one type of universe; they ran simulations for different cosmological models, including some where dark matter is "warm" and moves faster, and others where the universe expands differently. They treated these simulations like a giant cosmic laboratory, tracking billions of dark matter particles to see how they clump together.

The team's big breakthrough was realizing that the old way of measuring a halo's "peak height" (a fancy term for how rare and special a halo is in the cosmic crowd) was missing a crucial ingredient: the exact moment the halo was born. They introduced a new, "revised" peak height that accounts for the universe's growth rate at the time the halo formed. It's like realizing that to judge a runner's speed, you don't just look at their finish time, but also the wind speed and the track conditions at the moment they started running.

By using this new, more precise measurement, the authors found something amazing: a universal, tight relationship that holds true across almost all their simulations. They discovered that if you know the "peak height" of a halo's formation, you can predict its most likely concentration with incredible accuracy. The relationship follows a simple curve: as halos get "rarer" (higher peak height), their concentration drops, but it never goes below a certain "floor." This floor is a concentration value of about 2.54, meaning even the most spread-out halos have a minimum level of density.

The paper also addresses a few specific scenarios. For instance, they looked at halos in "Open CDM" universes (where the geometry of space is different) and found that these halos can be slightly more concentrated than the rule predicts, likely because the expanding universe stretches their outer edges. They also noted that for very massive halos or those formed very early, there might be a tiny uptick in concentration, but they suspect this is just a glitch from the computer resolution rather than a real physical law.

Ultimately, the authors provide a new, simple formula that acts like a cosmic translator. It takes the complex history of how a halo formed and translates it into a single, reliable number for its concentration. This formula works for different types of dark matter and different universe models, making it a powerful tool for future astronomers. They even made a software package available so anyone can use this new rule to predict how the invisible scaffolding of our universe is built. It's a step toward a unified theory of how the dark matter web holds the cosmos together, turning a chaotic scatter of data into a clear, predictable pattern.

Technical Summary: Universal Fitting Formulae for the Peak Concentration of Dark Matter Halos

Problem Statement
The prediction of structural properties of dark matter halos, specifically their concentration, is fundamental to modern cosmology and galaxy formation studies. While the Navarro-Frenk-White (NFW) profile is widely accepted, accurately predicting halo concentration remains challenging due to its dependence on dynamical state, environment, mass, redshift, formation history, and cosmological parameters. Previous attempts to link concentration to the halo peak height parameter ( $\nu$ ) or assembly history have yielded relations that are not fully universal across different cosmological models and redshift regimes. Existing models often lack a parameter that effectively incorporates the linear growth factor at the halo formation epoch, limiting their ability to track assembly history physically. Furthermore, the scatter in concentration relations often varies with simulation resolution and cosmology, hindering robust predictions.

Methodology
The authors utilized a comprehensive suite of N-body simulations spanning diverse cosmologies and box sizes to derive a universal prescription for halo concentration. The dataset includes:

Public Simulations: Data from the TNG project (TNG50, TNG100, TNG300) and the MultiDark-Planck project (VSMDPL, SMDPL, MDPL2, BigMDPL).
Custom Simulations: A suite of simulations covering Standard CDM (SCDM), Open CDM (OCDM), a scale-free model (SF-n2), $w$ CDM models (with $w = -0.816$ and $w = -1.294$ ), and Warm Dark Matter (WDM) models with varying thermal relic masses (0.5, 0.8, 1.1, 1.5 keV).
Resolution and Selection: Simulations range from $30$ to $2500 \, h^{-1}\text{Mpc}$ in box size with particle masses from $3.65 \times 10^5$ to $2.36 \times 10^{10} \, h^{-1}M_\odot$ . To ensure reliability, only halos resolved by at least 2000 dark matter particles were selected. NFW profile fitting was performed with inner radii set to the convergence radius ( $r_{\text{conv}}$ ) or $3\epsilon$ (whichever was larger) to minimize collisional relaxation effects.

Key Methodological Innovations

Revised Peak Height ( $\nu_f$ ): The authors modified the standard halo peak height definition to incorporate the linear growth factor at the halo formation epoch ( $D(z_f)$ ). The formation peak height is defined as $\nu_f(M, z, z_f) = \nu(M, z) / D(z_f)$ , where $z_f$ is the redshift at which the main progenitor reached half the final halo mass. This revision aims to better capture the dependence of concentration on assembly history.
Peak Value (PV) Analysis: Recognizing that both halo concentration and formation time distributions follow lognormal profiles, the authors focused on the "peak value" (mode) of these distributions within narrow mass bins, rather than the mean or median. They argue the mode is more representative for lognormal distributions.
Effective Peak Height ( $\nu_{\text{eff}}$ ): To account for the suppression of low-mass halo formation in WDM cosmologies, an effective peak height was defined: $\nu_{\text{eff}} = \nu_{f,\text{peak}} (1 + a M_{\text{hm}}/M)$ , where $M_{\text{hm}}$ is the half-mode mass.

Key Results

Universal Relation: The authors derived a universal, tight relation between the peak concentration ( $c_{\text{peak}}$ ) and the effective peak height ( $\nu_{\text{eff}}$ ). This relation is invariant to redshift, box size, initial power spectrum, and cosmology, exhibiting remarkably small scatter.
Fitting Formula: The best-fit relation for virial overdensity ( $\Delta_{\text{vir}}$ ) is given by:
$c_{\text{peak}} = 12.03 \, \nu_{\text{eff}}^{-1.09} + 2.54$
For the $\Delta_{200c}$ overdensity, the relation is:
$c_{\text{peak}} = 8.81 \, \nu_{\text{eff}}^{-0.95} + 2.31$
Concentration Floor: The relation flattens systematically with increasing revised peak height, approaching an asymptotic value (concentration floor) of $2.54$ (for $\Delta_{\text{vir}}$ ), consistent with previous reports.
Cosmological Dependence: The formula successfully captures dependencies on redshift, halo mass, and cosmological parameters (including WDM suppression scales) through the variables $D(z)$ , $\sigma(M)$ , and $M_{\text{hm}}$ .
Resolution Effects: Comparisons across different resolution runs (e.g., within the MultiDark-Planck suite) show that concentration measurements are broadly consistent in overlapping mass ranges, with only modest residual offsets. A slight upturn in concentration for halos with very high formation peak heights is observed but is attributed to potential residual resolution effects rather than a robust physical signal.

Significance and Claims
The paper claims to provide a robust framework for predicting the most probable (peak) halo concentration using theoretical or semi-analytical mass assembly histories. By incorporating the linear growth factor at the formation epoch and utilizing the peak value of lognormal distributions, the proposed fitting formula offers a universal description that minimizes scatter across diverse cosmological models. The authors emphasize that this approach quantifies how concentration relations depend on redshift and cosmological parameters without requiring complex, non-universal corrections. A software package for calculating these peak concentrations is made publicly available to facilitate these predictions. The work does not claim to resolve all physical mechanisms driving concentration but offers a precise empirical tool for cosmological modeling.

More like this