Core Galactic Dark Matter Profiles

• Core galactic dark matter profiles describe the radial density of dark matter in galaxy halos, featuring both cuspy (e.g., NFW) and cored (e.g., Burkert, pISO) structures.
• Baryonic feedback and self-interacting dark matter transform initial cusps into cores, with simulations indicating core sizes of 1–5 kpc at optimal stellar-to-halo mass ratios.
• Empirical constraints from rotation curves and kinematic studies, alongside advanced simulations, provide insights to distinguish core versus cusp models and refine dark matter physics.

Core galactic dark matter profiles describe the radial structure of dark matter (DM) density in galactic halos, especially the central regions (≲few kpc), where key physical processes imprint signatures accessible to both observations and simulations. The form and diversity of these profiles play a critical role in ΛCDM cosmology, baryonic feedback studies, and in probing possible departures from cold, collisionless dark matter physics.

1. Parametric Forms of Core and Cusp Dark Matter Profiles

The density profile ρ(r) serves as the primary descriptor of the spatial distribution of dark matter in galactic halos. Canonical forms include:

Profile	Formula	Central Slope
NFW	$\rho(r) = \rho_s/(r/r_s)[1+(r/r_s)]^2$	$r^{-1}$ (cusp)
Einasto	$$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$	$\to$ flat for small $\alpha$; quasi-core for $\alpha\ll1$
Burkert	$\rho(r)=\rho_0/[(1+r/r_c)(1+(r/r_c)^2)]$	constant (core)
Pseudo-Isothermal	$\rho(r)=\rho_0/(1+(r/r_c)^2)$	constant (core)
Generalized/“Zhao”	$$\rho(r) = \rho_0(r/b_{\rm halo})^{-\gamma}[1+(r/b_{\rm halo})^\alpha]^{-(\beta-\gamma)/\alpha}$$	$r^{-\gamma}$

Cuspy profiles feature $r^{-1}$0 (NFW), while cored profiles exhibit $r^{-1}$1 as $r^{-1}$2 (Burkert, pISO). The Einasto and generalized double-power-law forms interpolate between these behaviors, allowing for arbitrary inner slopes, transition sharpness, and outer slopes (Sarkar et al., 4 Feb 2026, Hayashi et al., 2020, Lazar et al., 2020).

2. Physical Mechanisms for Core Formation

Theoretical and simulation work has established several baryonic and DM-intrinsic mechanisms capable of transforming initial cusps into (quasi-)cores.

2.1 Baryonic Feedback

Energetic feedback from repeated, impulsive star-formation-driven outflows (e.g., supernovae, radiation pressure) can generate “potential fluctuations” that heat the central DM phase-space, leading to the flattening of cusps. This process is maximally effective at $r^{-1}$3 (or $r^{-1}$4), producing cores of size $r^{-1}$5 kpc (core-Einasto, DC14 profiles) (Lazar et al., 2020, Brook et al., 2015, Tollet et al., 2015, Hayashi et al., 29 Jul 2025). Above this mass scale, deep potentials inhibit outflows; below it, feedback is energetically insufficient.

2.2 Dark Matter Self-Interaction (SIDM)

Self-interacting dark matter induces heat conduction into the inner halo, isotropizing velocities and converting initial cusps into constant-density cores. Simulations and analytic expectations yield core radii scaling with self-interaction cross-section: for $r^{-1}$6 cm² g⁻¹, $r^{-1}$7 kpc at $r^{-1}$8 (Ray et al., 2022). SIDM cores can be erased or modified by baryonic contraction in massive galaxies (Despali et al., 17 Dec 2025).

2.3 Alternative Exotic Physics

Scenarios such as annihilation-induced “late-time core formation” (reactivated in asymmetric DM after oscillations), ultra-light scalar DM (e.g., chameleon models, fuzzy DM), and quantum-pressure-supported fermionic or bosonic cores have been invoked to explain observed dark matter flattening in low-mass systems (Cline et al., 2020, Chanda et al., 2017, Siutsou et al., 2014). These models often generate cores on small scales (∼100 pc–kpc) within otherwise NFW-like outer profiles.

3. Empirical Constraints from Milky Way and Local Galaxies

Rotation curve and stellar kinematic data robustly constrain the form of central DM profiles:

• Classical Dwarf Spheroidals (dSphs): Non-spherical Jeans modeling finds a range of $r^{-1}$9 (inner slope) values: Draco ($$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$0), UMi ($$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$1), Leo I ($$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$2), Leo II ($$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$3) are strongly cusped; Carina, Sextans, Sculptor, Fornax allow for cores ($$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$4–0.77) within uncertainties (Hayashi et al., 2020).
• Milky Way and M31: Global rotation curve fits favor mildly cusped NFW or Hernquist inner profiles ($$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$5–1.2), with the Einasto index $$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$6–0.43 in the allowed range; large, constant-density Einasto cores are disfavored by the data (Kumar et al., 23 Jan 2025, Pato et al., 2015).
• Nearby Spirals and LSBs: Mass modeling yields acceptable fits with both cuspy and cored profiles; non-parametric inversion, however, closely aligns with NFW-like cusps in the inner kiloparsec, with cored models underpredicting central densities (Sarkar et al., 4 Feb 2026). Several LSB and dwarf galaxies retain evidence for cored profiles, especially at lower mass or stellar-to-halo mass ratios (Hague et al., 2014, Hayashi et al., 29 Jul 2025).

4. Simulation Insights and Mass-Dependent Diversity

Numerical hydrodynamical simulations (FIRE-2, NIHAO, AIDA-TNG) have mapped the diversity of DM cores/cusps across mass and redshift:

• The inner DM profile slope ($$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$7 or $$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$8) is set primarily by $$\rho(r) = \rho_{-2}\,\exp[-(2/\alpha)\{(r/r_{-2})^\alpha-1\}]$$9, evolving as halos undergo star formation and feedback. There is a prominent cusp-to-core-to-cusp sequence: low baryon fraction systems are cuspy, intermediate fractions develop and maintain cores, and highly massive halos recontract to cusps that may be steeper than NFW (Lazar et al., 2020, Tollet et al., 2015, Despali et al., 17 Dec 2025).
• Core-Einasto and related parameterizations capture the observed and simulated diversity, with empirical fitting functions for $\to$0 and concentration provided (Lazar et al., 2020). The core size $\to$1 peaks at $\to$2–$\to$3 kpc for $\to$4, declining at both higher and lower mass.
• Warm DM (WDM) produces only mild inner flattening, well fit by larger Einasto shape indices (α ≳ 0.22) but showing no explicit large cores; SIDM, in contrast, generates explicit cores absent baryons, which are significantly contracted in massive halos (Despali et al., 17 Dec 2025, Ray et al., 2022).

5. Observational Diagnostics and Degeneracies

Interpreting the presence of cores is complicated by degeneracies and projection effects:

• Photometry vs. Kinematics: Purely photometric (surface density) measurements cannot decisively distinguish 3D strong cores from weak cores or projected cusps; both produce flat inner profiles in projection. Only a joint analysis of velocity-dispersion data, anisotropy, and distribution-function consistency can robustly discriminate cored vs. cuspy DM halos (Valenciano et al., 17 Dec 2025).
• Circular Velocity Profiles: Cored halos produce flat, slowly rising $\to$5 in the central region; cuspy halos show a steeper inner rise. The circular velocity at the half-light radius is often insensitive to cusp/core structure, but the maximum circular velocity and detailed shape encode information about $\to$6 and potential flattening (Hayashi et al., 2020).
• Evolution: High-redshift ($\to$7) star-forming galaxies have smaller, denser cores than local ($\to$8) disks, consistent with ongoing “core expansion” driven by baryonic feedback over cosmic time (Sharma et al., 2021). The core scaling relations with $\to$9 and $\alpha$0 exhibit systematic shifts between epochs.

6. Constraints from Indirect Detection and Microphysical Models

Astrophysical quantities, such as the $\alpha$1-factor for DM annihilation and $\alpha$2-factor for decay, are sensitive to the shape of the central DM profile. Cuspy profiles (e.g., Draco dSph) yield higher $\alpha$3-factors; the spread in profile slope introduces large uncertainties in expected indirect detection signals (Hayashi et al., 2020).

The existence and properties of cores thus inform bounds on DM microphysics: strong (≳kpc-scale) cores across the galaxy mass function would favor SIDM or certain baryon–DM coupling scenarios; the observed diversity and mass dependence argue for a dominant baryonic-feedback origin, but leave open avenues for distinguishing among models via precision rotation curves and stellar kinematic studies (Hjorth et al., 2015, Cline et al., 2020, Chanda et al., 2017).

7. Synthesis and Outlook

Core galactic dark matter profiles are non-universal and encode crucial information about both dark matter microphysics and baryonic evolution. The preponderance of evidence indicates:

• The inner structure of galactic halos is diverse, spanning sharp cusps ($\alpha$4–1.5) in low- and high-mass systems, extended kpc-scale cores at intermediate mass, and a mass-dependent progression controlled by $\alpha$5.
• Baryonic feedback is the dominant driver of core formation, as substantiated by hydrodynamical simulations (FIRE-2, NIHAO) and the observed assembly of Local Group dwarfs and spirals (Hayashi et al., 29 Jul 2025, Brook et al., 2015, Tollet et al., 2015).
• Certain alternative DM physics (SIDM, annihilation/oscillation models, ultralight or fermionic DM) can produce cored profiles, but current data favor their effects only within stringent cross-section or mass bounds (Ray et al., 2022, Cline et al., 2020, Siutsou et al., 2014).
• The core–cusp problem is best regarded as a diagnostic probe of baryonic and DM physics on sub-kpc scales; refined kinematic analyses and systematic velocity-dispersion studies, particularly for classical dwarfs and low-mass spirals, remain essential to resolve residual degeneracies and further constrain the nature of dark matter.

Ongoing and next-generation wide-field spectroscopic campaigns, high-precision proper motion datasets, and sophisticated non-parametric modeling will progressively sharpen constraints, enabling robust discrimination among core-formation mechanisms and anchoring the astrophysical landscape of dark matter in galactic environments (Hayashi et al., 2020).