Inner Dark Matter Density Profiles

• Inner dark matter density profiles define how dark matter density changes in halo centers, quantified by the logarithmic slope α at small radii.
• Simulations and observations demonstrate that baryonic processes, initial conditions, and dark matter microphysics produce a range of profile shapes from cuspy to cored.
• Understanding these profiles is essential for resolving the cusp–core issue and testing dark matter models like SIDM and WDM across various mass scales.

Inner dark matter density profiles characterize the variation of dark matter density as a function of radius in the central regions of halos, typically parameterized by a logarithmic slope $\alpha = d\log\rho/d\log r$ evaluated at small radii. The properties of these inner profiles encode the interplay between initial conditions, dark matter microphysics, baryonic processes, and halo assembly history, and underpin persistent theoretical and observational issues including the cusp–core problem, halo contraction, and the apparent non-universality of density profiles across mass scales.

1. Definitions and Parametric Forms

The canonical description of inner halo structure is the power-law slope of the spherically averaged density profile: $\alpha(r) = \frac{d \log \rho}{d \log r}$ A variety of models describe structural diversity:

• Navarro–Frenk–White (NFW): $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$, with $\alpha\to -1$ as $r\to 0$ (cusp) (Dehghani et al., 2020).
• Generalized NFW (gNFW): $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$, where $\gamma$ is the asymptotic inner slope (Kollmann et al., 20 Oct 2025, Wasserman et al., 2017).
• Einasto profile: $$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$, with $\alpha_E$ controlling curvature, and $\alpha \to 0$ at $\alpha(r) = \frac{d \log \rho}{d \log r}$0 (core) (Salvador-Solé et al., 2011, Despali et al., 17 Dec 2025).
• Cored models: Burkert and similar profiles approach $\alpha(r) = \frac{d \log \rho}{d \log r}$1 const as $\alpha(r) = \frac{d \log \rho}{d \log r}$2.
• SIDM-specific forms: Explicitly cored, with core size set by cross section and mass (Despali et al., 17 Dec 2025, Kahlhoefer et al., 2019).
• Double power laws and entropy-based forms: Used to describe toy models or simulations with explicit initial entropy, e.g. $\alpha(r) = \frac{d \log \rho}{d \log r}$3 (Pilipenko et al., 2012).

Best-fit slopes and core sizes are typically quoted for $\alpha(r) = \frac{d \log \rho}{d \log r}$4 kpc scales in galaxies, $\alpha(r) = \frac{d \log \rho}{d \log r}$5 tens of parsecs in dwarfs, and $\alpha(r) = \frac{d \log \rho}{d \log r}$6 tens-hundreds of kpc in clusters, with the inner profile defined locally or as an average over a resolved interval (Relatores et al., 2019).

2. Theoretical Predictions and Simulation Results

Dissipationless CDM

Collisionless ΛCDM simulations robustly produce cuspy profiles with $\alpha(r) = \frac{d \log \rho}{d \log r}$7 (NFW) at $\alpha(r) = \frac{d \log \rho}{d \log r}$8 across a wide dynamic range in mass and redshift (Brun et al., 2017, Muni et al., 2023, Despali et al., 17 Dec 2025). Inner slopes measured at $\alpha(r) = \frac{d \log \rho}{d \log r}$9 in massive cluster simulations are $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$0, stable from $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$1 to $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$2 with low scatter (0.15 dex) (Brun et al., 2017). Dynamical, orbit-averaged estimators confirm that cuspy behavior persists down to the softening scale in well-relaxed halos, with uncertainties reduced by an order of magnitude relative to binned counts (Muni et al., 2023).

Baryonic Feedback and Halo Contraction

Hydrodynamical simulations and semi-empirical models demonstrate that baryonic physics strongly modifies the inner slope $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$3:

• Adiabatic contraction from dissipative gas inflow and star formation steepens profiles, yielding $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$4 in massive early-type galaxies, compared to NFW ($\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$5) (Chae et al., 2012, Grillo, 2012). Observed lensing and kinematics demand these contracted inner profiles, with enhanced densities within the effective radius by factors 3–4 at $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$6 (Chae et al., 2012).
• Supernova-driven outflows and bursty star formation inject energy into the DM core, flattening the profile ($\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$7) for a range of halo masses and star formation histories (Tollet et al., 2015, Sarrato-Alós et al., 23 Jan 2026). The efficiency of this process peaks for $\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$8–$\rho(r) = \rho_s / [(r/r_s)(1 + r/r_s)^2]$9 (Tollet et al., 2015, Sarrato-Alós et al., 23 Jan 2026).
• Temporal evolution: The same halo may cycle through cored and cuspy states. Halos of MW mass can form large cores at $\alpha\to -1$0 and recontract to steep cusps at late times if stellar mass inflow dominates (Tollet et al., 2015).
• AGN feedback can drive core formation and subsequent regeneration of cusps at late times in massive halos, with a three-phase evolution: baryonic contraction ($\alpha\to -1$1; steep slopes), quasar-mode core formation ($\alpha\to -1$2; flattened slopes), and cusp regeneration via radio-mode feedback or further contraction ($\alpha\to -1$3; re-steepening) (Peirani et al., 2016).

Small-Scale Perturbations and Non-Universality

Toy N-body models show that adding sufficiently strong small-scale initial perturbations can flatten the inner profile, raising the protohalo's entropy and yielding $\alpha\to -1$4 for specific ranges ($\alpha\to -1$5, early collapse) (Pilipenko et al., 2012). Only under such initial conditions are genuine, long-lived cores expected in collisionless simulations.

Dependence on Halo Mass and Assembly

Both semi-analytic and numerical approaches predict systematic variation of inner slope with mass and redshift:

• In models including baryons, $\alpha\to -1$6 varies from $\alpha\to -1$7 (core) in dwarfs, $\alpha\to -1$8–$\alpha\to -1$9 in MW-sized halos, to $r\to 0$0 in clusters (Popolo, 2010, Chae et al., 2012). Non-universality is a generic outcome when baryon effects or initial conditions are included (Popolo, 2010).

3. Observational Constraints across Mass Scales

Disk Galaxies and Dwarfs

Rotation curve modeling using high-quality samples (e.g., SPARC, coadded rotation curves) consistently disfavors pure NFW inner profiles. The best-fit inner slopes in late-type disks and dwarfs span $r\to 0$1, with many systems forming constant-density cores ($r\to 0$2), particularly at $r\to 0$3–$r\to 0$4 km/s (Dehghani et al., 2020, Relatores et al., 2019, Hayashi et al., 29 Jul 2025). Bayesian frameworks fitting flexible gNFW or Zhao profiles show a pronounced diversity, with $r\to 0$5 (inner slope) ranging from $r\to 0$6 (strong cores) to $r\to 0$7 (steep cusps) within the same population (Hayashi et al., 29 Jul 2025).

Early-Type Galaxies and Clusters

Strong lensing, stellar kinematics, and dynamical modeling in massive early-types and clusters indicate systematically steeper than NFW inner profiles. Projected logarithmic slopes inferred from ensembles of lens galaxies lie near $r\to 0$8 (Chabrier IMF), corresponding to $r\to 0$9 (Grillo, 2012). Direct Jeans analyses yield even steeper $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$0 with significant halo-to-halo scatter (Chae et al., 2012). In clusters, simulations and X-ray/strong-lensing-based studies converge on NFW-like cusps with no convincing evidence for cores at $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$1 (Brun et al., 2017).

Dwarf Spheroidals and Milky Way Satellites

Comprehensive Jeans modeling (including axisymmetry) for classical MW dSphs finds that most favor cuspy (NFW-like) or mildly cuspy profiles (median $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$2), with robust cores allowed only in a minority (e.g., Fornax, Sculptor) (Hayashi et al., 2020). The correlation between $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$3 and $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$4 aligns with core formation scenarios from baryonic feedback, but a genuine diversity in observed slopes is confirmed. Ultra-faint dwarfs exhibit high densities at sub-kpc scales—populating both the upper and lower envelope of predicted core/cusp possibilities (Zavala et al., 2019, Kahlhoefer et al., 2019).

4. Impact of Dark Matter Microphysics: SIDM and Alternatives

Self-interacting dark matter (SIDM) and, to a lesser extent, warm dark matter (WDM) fundamentally alter the physics of core formation:

• SIDM (constant cross section): DMO SIDM simulations generically produce large, constant-density cores (inner slope $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$5), with core radii and densities set by the self-interaction rate and mass (Despali et al., 17 Dec 2025, Kahlhoefer et al., 2019). At large enough cross sections ($$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$6), high-concentration halos undergo gravothermal core collapse, steepening the inner slope to $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$7 (Kahlhoefer et al., 2019).
• SIDM + baryons: Baryon-driven contraction in central regions partially counters self-interaction, erasing cores in more massive halos and yielding a broader possible distribution of inner slopes—from cored to “over-cuspy” profiles ($$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$8) in compact galaxies (Despali et al., 17 Dec 2025).
• Velocity-dependent SIDM: Bimodal populations—cuspy collapsed vs. cored—can emerge in satellite halos, naturally explaining the extreme density diversity of MW ultra-faints (Zavala et al., 2019, Kahlhoefer et al., 2019).
• WDM: Thermal relic WDM models yield mild inner flattening of profiles in low-mass halos (inner $$\rho(r) = \rho_0 (r/r_s)^{-\gamma} (1 + r/r_s)^{\gamma-3}$$9 decreases from $\gamma$0 to $\gamma$1 at $\gamma$2), but lack significant core formation (Despali et al., 17 Dec 2025).

Strong lensing predictions are highly sensitive to the inner slope: subhalos with $\gamma$3 (steep cusps from SIDM collapse) are detectable down to $\gamma$4-fold lower masses relative to standard NFW cases and maintain detectability even against macro-model degeneracies, making measurement of the inner profile slope a direct probe of dark sector physics (Kollmann et al., 20 Oct 2025).

5. Role of Star Formation Histories and Baryonic Feedback

Detailed analysis of simulated galaxy populations demonstrates that the burstiness and temporal duration of star formation critically shape the range of realized inner slopes:

• Bursty and extended SFHs maximize core formation, as strong, rapidly varying central potentials transfer energy to the dark matter (Sarrato-Alós et al., 23 Jan 2026).
• Early, smooth SFHs tend to preserve or regenerate cusps, particularly in high-mass or quenched systems.
• General relations expressing $\gamma$5 as a function of $\gamma$6 can be improved by incorporating SFH metrics, reducing the scatter in predicted slopes (Sarrato-Alós et al., 23 Jan 2026).
• Both NIHAO and FIRE-2 simulation suites, when analyzed homogeneously, exhibit minimized tensions and consistent locus/dispersion in the $\gamma$7–$\gamma$8 plane.

The mass-dependent core formation efficiency, peaking in galaxies with $\gamma$9–$$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$0, underpins the non-universality of profiles and matches the mass regime where observations reveal most pronounced departure from NFW (Tollet et al., 2015, Hayashi et al., 29 Jul 2025).

6. Cusp–Core Demographics and the Limits of Universality

The assembled body of simulation and observational evidence reveals:

• Clusters: Universally cuspy, NFW/Einasto profiles (Brun et al., 2017, Despali et al., 17 Dec 2025).
• Groups/Early-Type Galaxies: Steepened inner profiles due to baryonic contraction, $$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$1–$$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$2 (Chae et al., 2012, Grillo, 2012).
• Milky Way–mass Disks: Strong diversity; cores and cusps coexist, with the mean slope shallower than NFW ($$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$3–$$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$4) (Dehghani et al., 2020, Hayashi et al., 29 Jul 2025, Relatores et al., 2019).
• Dwarfs/Satellites: Genuine spread, from core-collapsed ($$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$5), to cored ($$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$6), to NFW-like ($$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$7), controlled by SFH, self-interaction, and dynamical state (Hayashi et al., 2020, Zavala et al., 2019, Kahlhoefer et al., 2019).
• Low-mass dwarfs ($$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$8): Observed slopes shallower still ($$\rho(r) = \rho_{-2} \exp\left\{ -\frac{2}{\alpha_E} [(r/r_{-2})^{\alpha_E} - 1] \right\}$$9–$\alpha_E$0), demanding robust core-formation by feedback or DM microphysics.

Comparisons across samples show that universality of the inner profile—long assumed for CDM—breaks down when baryonic effects, assembly bias, or non-trivial microphysics are incorporated (Popolo, 2010, Tollet et al., 2015, Despali et al., 17 Dec 2025, Hayashi et al., 29 Jul 2025).

7. Interpretational Challenges and Future Directions

Observational uncertainties—including pressure support, non-circular gas motions, projection/deprojection effects, and limited spatial resolution—complicate precise measurement of the innermost slopes (Hayashi et al., 29 Jul 2025, Relatores et al., 2019). Key open questions include:

• Can hydrodynamical simulations reproduce both the observed median and scatter in slopes across mass scales?
• What is the dominant baryonic process shaping cores, especially at high masses?
• Can the diversity of observed inner slopes be reproduced naturally within standard or alternative DM models, absent fine-tuning?
• How do core–cusp transitions relate to measurable properties (SFH, compactness, pericenter history, dynamical times)?
• What are the limits on DM self-interaction cross section or power-spectrum cutoff imposed by the inner structure in dwarfs and satellites?

Future progress hinges on high-fidelity kinematic data (Subaru/PFS, Roman, Gaia), more precise SFH measurements, and statistical lensing surveys (Euclid, JWST) targeting population statistics of subhalo inner slopes. The mapping of density profile shape to observable quantities remains central for distinguishing ΛCDM from its alternatives and for using galactic centers as laboratories for fundamental physics and galaxy formation (Hayashi et al., 29 Jul 2025, Kollmann et al., 20 Oct 2025, Despali et al., 17 Dec 2025).