Cusp-Core Transition Model

Updated 27 January 2026

Cusp-core transition model is a framework describing how dark matter halos evolve from steep, cuspy NFW profiles to flat, core-like structures through energy redistribution mechanisms.
It quantitatively maps the transformation to Burkert-like cores while conserving mass and density, reproducing key observational scaling relations.
Key mechanisms include baryonic feedback (e.g., supernova and AGN), dynamical heating by massive perturbers, and non-standard dark matter physics.

The cusp-core transition model addresses the longstanding discrepancy between the steep, cuspy central dark matter profiles predicted by cold dark matter (CDM)-only simulations—typically exhibiting $\rho \propto r^{-1}$ central slopes (NFW)—and the much shallower, often nearly constant-density cores inferred from galaxy rotation curves and stellar kinematics, particularly in dwarfs and low-surface-brightness systems. The model describes the physical processes and structural mapping by which an initially cuspy halo is transformed into a flat-core profile. Multiple distinct mechanisms are proposed and quantitatively explored across the literature, including energetic baryonic feedback (from supernovae or active galactic nuclei), dynamical heating (e.g. by massive perturbers), nonstandard dark matter properties (e.g. self-interactions, ultra-light scalar fields), and sharp dark-sector phase transitions. The model both encodes the physics of the transformation and provides a quantitative pathway for connecting simulated and observed structural scaling relations across the dark-halo mass spectrum.

1. From Cuspy to Cored Profiles: Structural Mapping

The initial state of a dark matter halo in the standard CDM paradigm is well described by the Navarro–Frenk–White (NFW) profile: $\rho_\text{NFW}(r) = \frac{\rho_s}{(r/r_s)(1 + r/r_s)^2}$ with two free parameters: the scale density ( $\rho_s$ ) and the scale radius ( $r_s$ ). The central logarithmic slope is $d\ln\rho/d\ln r\rightarrow -1$ as $r\to 0$ , i.e., a cusp.

In the cusp-core transition, some physical process redistributes dark matter in the central region—leaving the outer halo largely intact—so that the inner profile is shallower ( $d\ln\rho/d\ln r\to 0$ ). The transformed density is often modeled by the Burkert profile: $\rho_\text{B}(r) = \frac{\rho_0\,r_0^3}{(r + r_0)(r^2 + r_0^2)}$ where $r_0$ is the core radius and $\rho_0$ is the core density. The mapping is physically constrained by continuity of the mass and density at the virial radius and by conservation of halo mass, up to baryonic mass loss.

For a full analytic mapping, the dimensionless ratio $\eta = r_\text{core}/r_s$ is found by solving

$g(\eta) \equiv \frac{\left(c_{200}+\eta\right)\left(c_{200}^2+\eta^2\right)} {c_{200}(1+c_{200})^2} - \frac{F_N(c_{200})}{F_B(c_{200}/\eta)} = 0$

where $F_N$ and $F_B$ are the mass integrals for NFW and Burkert, respectively (Shinozaki et al., 20 Jan 2026, Kaneda et al., 2024).

2. Physical Mechanisms Driving the Transition

2.1 Supernova and Stellar Feedback

In galaxies with sufficient stellar content, repeated cycles of star formation and supernovae (SN) feedback expel and re-accrete gas in the central regions. The resulting time-varying gravitational potential irreversibly transfers energy to collisionless dark matter (and stars) via the impulsive approximation: $\Delta E \simeq -\langle \Delta\Phi \rangle$ per burst. Recurring episodes operate as a random walk in phase-space and heat the cusp into a core (Teyssier et al., 2012). The fundamental core-formation criterion is resonance: the oscillation period $T_\mathrm{osc}$ of the baryon potential should match the local dark matter dynamical time $t_\mathrm{dyn}(r_\text{core})$ (Ogiya et al., 2012).

The energetics can be encapsulated as: $\Delta E = \epsilon\, E_\mathrm{SNII}$ where $\epsilon$ is the energy-conversion efficiency, the fraction of SN energy that heats the dark matter. Empirical calibration yields $\epsilon \sim 0.01$ for SPARC galaxies (Shinozaki et al., 20 Jan 2026).

2.2 AGN Feedback

In massive halos, active galactic nuclei (AGN) feedback drives cycles of gas expulsion and fallback, generating large potential fluctuations on timescales comparable to the local dynamical time. Each cycle unbinds a fraction of the central dark matter: $\Delta E_\mathrm{AGN} \gtrsim \Delta M_\mathrm{gas} \Phi_\mathrm{bind}$ Repeated cycles cumulatively flatten the inner slope from $\gamma\sim -1$ to $\gamma \sim -0.5 ... -0.7$ , producing a core of size $r_\text{core} \sim 5-10$ kpc in clusters (Martizzi et al., 2012).

2.3 Dynamical Heating by Massive Perturbers

The passage of massive baryonic clumps, globular clusters, or primordial black holes through the halo center transfers orbital energy to the dark matter via dynamical friction or two-body relaxation, leading to a reduction of the central density: $t_\text{core} \sim \kappa\, T_r(r_c)$ where $T_r$ is the local relaxation time and $\kappa \sim 100$ –$300$ (Boldrini et al., 2019, Boldrini et al., 2019). The resulting core size and duration depend on perturber mass, orbit, and frequency of central passages.

2.4 Non-Standard Dark Sector Mechanisms

Bound Dark Matter (BDM): Phase transition at a critical energy density ( $E_c$ ) transforms dark matter from massive (CDM) to massless (HDM) above a threshold, enforcing a cored profile with $r_c \propto 1/E_c^4$ (Mastache et al., 2011).
Late-Time Annihilation via Dark Sector Oscillations: Reactivation of annihilation in asymmetric dark matter models with a small DM-number violating mass term erases cusps post-structure formation (Cline et al., 2020).
Ultra-light dark matter: Produces soliton (cored) centers, but the correlation between core density and radius ( $\rho_c \propto 1/R_c^\beta$ ) does not match observed $\beta \sim 1$ , instead predicting $\beta=4$ or higher (Deng et al., 2018, Kendall et al., 2019).
Self-interacting DM: Heat exchange flattens the cusp inside a characteristic radius; the relevant cross section is $\sigma/m \sim 0.1$ – $1\,\mathrm{cm}^2/\mathrm{g}$ (Popolo et al., 2022).

3. Scaling Relations and Mass-Dependence

The cusp-core transition model predicts scale-dependent outcomes for core formation. Baryonic feedback transforms cusps to cores only in halos with sufficient potential energy to capture SN-driven energy (typically $10^8\,M_\odot \lesssim M_{200} \lesssim 10^{11}\,M_\odot$ for $\epsilon\sim0.01$ ), with both ultra-faint dwarfs (insufficient star formation) and clusters (too deep potentials) remaining cuspy (Shinozaki et al., 20 Jan 2026, Hayashi et al., 29 Jul 2025).

The mapping from NFW to Burkert or pseudo-isothermal core produces scaling relations: $\Sigma_\text{core} \propto M_{200}^{0.30},\quad V_\mathrm{max}^\text{(core)} \propto M_{200}^{0.32},\quad r_\text{core}\propto M_{200}^{0.34}$ with near-universality satisfied up to logarithmic corrections (Kaneda et al., 2024).

A key empirical finding is the central surface density universality ( $\mu_{0D} \equiv \rho_0 r_0 \simeq 140\,M_\odot\,\mathrm{pc}^{-2}$ across a wide range of galaxies), naturally reproduced within the transition model as a consequence of mass-conserving inner rearrangement and the weak mass-dependence of halo concentration parameters at fixed redshift (Ogiya et al., 2013).

4. Observational Tests and Diagnostics

The mass-dependent outcome of the cusp-core transition is revealed by rotation-curve and stellar-kinematic fits. Bayesian analyses using flexible six-parameter DM profiles on SPARC data show a spectrum of inner slopes $\gamma$ from near-core ( $\sim 0$ ) to steep cusp ( $\sim 2$ ), with a pronounced transition around Milky Way–mass spirals ( $V_\mathrm{max}\sim100\,\mathrm{km/s}$ ) and diverging toward cuspy NFW values in ultra-faint dwarfs and clusters (Hayashi et al., 29 Jul 2025). The locus of galaxies in the $\Sigma_\text{DM}( <0.01\,r_{V_\mathrm{max}})$ – $V_\mathrm{max}$ plane shows a central surface-density dip at intermediate masses, consistent with simulated baryonic feedback predictions.

Time-resolved simulations predict specific signatures:

Bursty, duty-cycled star formation ( $\Delta t_\text{burst} \sim t_\text{dyn}$ , $\text{SFR}_\text{peak} / \text{trough} \sim 5$ –10) (Teyssier et al., 2012).
Hot, thick stellar disks with $v/\sigma \sim 1$ .
Extended globular clusters ( $R_\mathrm{eff} > 10$ pc) and reduced tidal debris in systems with established cores (Orkney et al., 2019).
Cusp-core cycles and partial regeneration driven by perturber orbital timescales (Boldrini et al., 2019).

5. Universality, Limitations, and Controversies

The cusp-core transition model, under reasonable physical assumptions (mass and density conservation outside the core, rapid central rearrangement, neglecting outer mass loss in the simplest mappings), explains:

The empirical central surface-density $\mu_{0D}$ and Strigari $M(<300\,\mathrm{pc})$ relations for dwarf galaxies (Ogiya et al., 2013).
The lack of universal core formation despite universal NFW initial conditions: only galaxies with sufficient star formation and feedback efficiency experience transition (Shinozaki et al., 20 Jan 2026).
The observed diversity in central inner slopes at fixed $V_\mathrm{max}$ reflects stochasticity in SN coupling, SFR, and halo assembly.

However, the precise physical mechanism remains debated. While baryonic feedback is strongly favored in the intermediate-mass regime (Hayashi et al., 29 Jul 2025), purely dark-matter–based scenarios (self-interactions, PBHs, phase transitions, late-time annihilations) are not excluded and may dominate at the low- and high-mass extremes (Boldrini et al., 2019, Mastache et al., 2011, Cline et al., 2020). Ultra-light dark matter models' inability to reproduce the observed $\rho_c - R_c$ ( $\beta\approx 1$ ) scaling with stability further motivates baryonic or alternative dark-sector origins for the transition (Deng et al., 2018).

6. Future Directions and Observational Prospects

Next-generation facilities (SKA, ngVLA, TMT, JWST, ELT) will enable kinematic and lensing measurements at $\sim 10$ –100 pc scales, potentially distinguishing subtle structural distinctions (e.g., core radius, profile shape) and time-dependent features. Hydrodynamical cosmological simulations are refining estimates of feedback coupling efficiency $\epsilon$ and the stochasticity of core formation. Combining high-quality spatially resolved rotation curves, strong-lensing mass profiles, and chemodynamical data will be critical to further constrain the landscape of permissible transition mechanisms (Hayashi et al., 29 Jul 2025, Kaneda et al., 2024).

7. Summary Table: Mapping Mechanisms and Predictive Regimes

Mechanism	Core Formation Scale	Key Diagnostic
Supernova feedback	$10^{8-11}\,M_\odot$ halos	Bursty SFH, $v/\sigma \sim 1$
AGN feedback	Clusters ( $\gtrsim 10^{13}\,M_\odot$ )	Core $r_c \sim$ 5-10 kpc, rapid central gas cycles
Dynamical heating (GCs, PBH)	Dwarfs ( $10^7-9\,M_\odot$ )	Transient core, GC size ( $R_\text{eff} > 10$ pc)
Bound dark matter	Model-dependent, core set by $E_c$	$r_c \propto 1/E_c^4$
Ultra-light DM	Inconsistent with observed scaling	Soliton core, but $\beta \neq 1$
DM annihilation (oscillations)	Centrally peaked halos	Flat core, suppressed inner dispersion

Each of these mechanisms is supported by analytic models and/or numerical simulations providing quantitative core scalings and energetics. Concordance between simulation predictions and observed scaling relations is achieved for baryon-driven transition models in the mass regime of star-forming spirals, whereas alternative scenarios gain traction in the least and most massive halos, or in the presence of exotic dark matter properties (Kaneda et al., 2024, Mastache et al., 2011, Cline et al., 2020, Orkney et al., 2019, Boldrini et al., 2019).

References: