--- title: Cusp-to-Core Transition in Dark Matter Halos url: https://www.emergentmind.com/topics/cusp-to-core-transition type: topic --- # Cusp-to-Core Transition in Dark Matter Halos The cusp-to-core transition denotes the transformation of the inner density structure of a dark-matter halo from a centrally divergent cusp, typically of Navarro–Frenk–White type, to a shallow or approximately constant-density core. It is the dynamical content behind the broader core–cusp problem: collisionless $\Lambda$CDM simulations predict inner profiles close to $\rho \propto r^{-1}$, whereas many observed dwarf, low-surface-brightness, and late-type galaxy rotation curves favor much flatter central mass distributions. Recent cross-scale analyses further indicate that the phenomenon is not universal in the sense of all halos becoming cored; rather, halo structure appears to depend on mass scale, with late-type galaxies often showing shallower-than-NFW inner profiles, while lower-mass dwarf spheroidals and higher-mass groups and clusters tend to remain cuspy [0910.3538] [2507.22155]. ## 1. Formal definition and parametrization A cusp-to-core transition is usually expressed through the inner logarithmic density slope $$ \alpha \equiv \left.\frac{d\ln\rho}{d\ln r}\right|_{r\to 0}. $$ A cuspy profile has $\alpha \approx -1$, as in the NFW form, whereas a truly cored profile has $\alpha \approx 0$. In the generalized Hernquist fits used by Hayashi et al., the free parameter $\gamma$ plays exactly this role, with $\gamma \simeq -\alpha$; thus $\gamma \to 1$ denotes a steep cusp and $\gamma \to 0$ a core [2507.22155]. The canonical cuspy reference model is the NFW profile, $$ \rho_{\rm NFW}(r)=\frac{\rho_s}{(r/r_s)(1+r/r_s)^2}, $$ whose asymptotic inner behavior is $\rho \propto r^{-1}$. Two widely used cored alternatives are the pseudo-isothermal profile, $$ \rho_{\rm PI}(r)=\frac{\rho_0}{1+(r/r_c)^2}, $$ and the Burkert profile, $$ \rho_{\rm B}(r)=\frac{\rho_0}{(1+r/r_c)(1+(r/r_c)^2)}, $$ both of which approach an approximately constant central density in the inner region [2209.14151]. Two surface-density-like quantities are central in the literature. Ogiya et al. define the cored-halo central surface density as $$ \mu_{0D}\equiv \rho_0 r_0, $$ and show that it can remain nearly mass-independent under simple cusp-to-core mappings [1309.1646]. Hayashi et al. instead adopt a mass-scaled central quantity designed for cross-comparison from dwarfs to clusters, $$ \Sigma_{\rm DM}(<0.01r_{V_{\rm max}})\equiv \frac{M_{\rm DM}(<0.01r_{V_{\rm max}})}{\pi [0.01r_{V_{\rm max}}]^2}, $$ with $$ M_{\rm DM}(<0.01r_{V_{\rm max}})=\int_0^{0.01r_{V_{\rm max}}}4\pi r'^2\rho_{\rm DM}(r')\,dr'. $$ The choice of $0.01\,r_{V_{\rm max}}$ is a compromise: it lies well inside the region that rotation curves typically resolve across galaxy types, while still scaling with halo size so as not to bias comparisons of dwarfs, spirals, and clusters [2507.22155]. ## 2. Observational status across galaxy mass scales The empirical basis for the cusp-to-core transition was established from high-resolution studies of low-surface-brightness and gas-rich dwarf galaxies. HI synthesis data with beam $\sim 10$–$20''$ and linear resolution $\sim 1$–$2\,{\rm kpc}$, together with H$\alpha$ long-slit observations with seeing $\sim 1$–$2''$ and $\lesssim 200\,{\rm pc}$ resolution, showed slowly rising inner rotation curves and inferred inner density slopes $\alpha_{\rm obs}\approx 0$ to $-0.3$ in the inner $0.1$–$1\,{\rm kpc}$. Two-dimensional velocity fields found non-circular motions of typical amplitude $\lesssim 5\,{\rm km\,s^{-1}}$ and photometric–kinematic center offsets $\lesssim 200\,{\rm pc}$, both too small to hide a genuine cusp in all but contrived cases [0910.3538]. Hayashi et al. extended the problem to a uniform, mass-dependent comparison using 115 high-quality SPARC galaxies with inclination $>30^\circ$ and at least 9 data points. Their best-fit generalized-Hernquist inner slopes span $\gamma=0.01$ to $1.96$. Roughly half the sample has $\gamma \lesssim 0.5$, much shallower than NFW, and typical median values are $\gamma \simeq 0.3$–$0.6$, with cores $(\gamma<0.2)$ in many Milky Way-mass disks. In the same framework, Milky Way dwarf spheroidals, including classical dSphs and ultra-faint dwarfs, generally lie near the cuspy NFW prediction in the $\Sigma_{\rm DM}$–$V_{\rm max}$ plane and have $\gamma$ in the range $0.5$–$1.5$, while galaxy groups and clusters show $\Sigma_{\rm DM}(<0.01r_{V_{\rm max}})\sim 500$–$1\,000\,M_\odot\,{\rm pc}^{-2}$ at $V_{\rm max}\sim 500$–$2\,000\,{\rm km\,s^{-1}}$, again matching a cuspy NFW baseline. Late-type SPARC galaxies with $V_{\rm max}\sim 50$–$250\,{\rm km\,s^{-1}}$ scatter systematically below the NFW line and track FIRE-2 core predictions, whereas ultra-faint dwarfs with $V_{\rm max}<30\,{\rm km\,s^{-1}}$ and clusters with $V_{\rm max}>500\,{\rm km\,s^{-1}}$ revert to cuspy behavior. This yields a mass-dependent “cusp-to-core-to-cusp” pattern across $V_{\rm max}\approx 10 \to 2\,000\,{\rm km\,s^{-1}}$ [2507.22155]. That pattern is significant because it disfavors the idea of a single, mass-independent inner profile family. It also narrows the transition region to galaxy scales where baryonic processes are dynamically important but not overwhelmingly suppressed by either shallow star-formation budgets or deep cluster potentials. ## 3. Inference from rotation curves and mass-scaled halo diagnostics Modern analyses infer cusp-to-core behavior from rotation-curve decomposition rather than from direct density measurements. In the SPARC analysis, each observed rotation curve is modeled as $$ V_{\rm tot}^2=V_{\rm DM}^2+\Upsilon_{\rm disk}V_{\rm disk}^2+\Upsilon_{\rm bul}V_{\rm bul}^2+V_{\rm gas}^2. $$ The dark component is obtained by numerically integrating an axisymmetric generalized Hernquist density re-parameterized by $(V_{200},\,C_{200},\,Q,\,\alpha,\,\beta,\,\gamma)$, where $\gamma$ is the inner slope and $\alpha$ controls the sharpness of the transition to outer slope $\beta$ [2507.22155]. The priors are deliberately broad. For the dark halo, flat priors are adopted over $$ 10\,{\rm km\,s^{-1}}\le V_{200}\le 500\,{\rm km\,s^{-1}},\quad 0\le C_{200}\le 500,\quad 0.1\le Q\le 2, $$ $$ 0.5\le \alpha \le 3,\quad 3\le \beta \le 10,\quad 0\le \gamma \le 2. $$ The stellar mass-to-light ratios follow log-normal priors $\Upsilon_{\rm disk}=0.5\pm 0.1\,{\rm dex}$ and $\Upsilon_{\rm bul}=0.7\pm 0.1\,{\rm dex}$, while distance $D$ and inclination $i$ are also varied, entering as $V_k\propto \sqrt{D/D_0}$ and $V_{\rm obs}\propto \sin i_0/\sin i$, with Gaussian priors from SPARC catalog uncertainties. The posterior is sampled with `emcee` under $$ \mathcal{L}\propto \exp(-\tfrac12\chi^2),\qquad \chi^2=\sum_R\left[\frac{V_{\rm obs}-V_{\rm tot}}{\delta V_{\rm obs}}\right]^2. $$ This construction propagates uncertainties in $\Upsilon$, $D$, and $i$ into the posterior rather than treating them as fixed nuisance choices [2507.22155]. The importance of the mass-scaled surface density $\Sigma_{\rm DM}(<0.01r_{V_{\rm max}})$ lies in comparative inference. A fixed physical aperture would mix systems at different fractions of their characteristic halo radii; the chosen aperture instead normalizes to halo size and allows dwarfs, spirals, groups, and clusters to be placed on a common structural plane. Within that plane, departures below the NFW baseline can be interpreted as a compact empirical signature of core formation, while proximity to the baseline indicates retention of cuspy structure [2507.22155]. ## 4. Baryonic routes to core formation Within $\Lambda$CDM, the most developed physical explanation for the cusp-to-core transition is repeated baryonic feedback. In hydrodynamical simulations such as FIRE-2, repeated supernova outbursts generate rapid, oscillatory changes in the central potential. Dark-matter particles gain energy through resonant interactions, and the inner profile flattens toward $\gamma\to 0$. Hayashi et al. emphasize that core-formation efficiency peaks at stellar-to-halo mass ratio $M_\star/M_{\rm halo}\sim 10^{-2}$, corresponding to $V_{\rm max}\sim 50$–$150\,{\rm km\,s^{-1}}$, which coincides with the mass range where the SPARC galaxies most clearly depart from the cuspy baseline [2507.22155]. A controlled demonstration was provided by Teyssier et al. in an isolated adaptive-mesh-refinement experiment. They evolved an initially cuspy halo with $M_{200}=1.4\times 10^{10}\,M_\odot$, $V_{200}=35\,{\rm km\,s^{-1}}$, and concentration $c=10$, using RAMSES with quasi-Lagrangian refinement to maximum level $\ell_{\max}=14$, corresponding to $\Delta x=18\,{\rm pc}$. Their feedback model introduces a non-thermal energy reservoir to represent unresolved turbulence, cosmic rays, and magnetic fields, with dissipation time $t_{\rm diss}\simeq 10\,{\rm Myr}$ and cooling suppressed where $\sigma_{\rm turb}>10\,{\rm km\,s^{-1}}$. In the run with feedback, a clear core appears after $\sim 0.5\,{\rm Gyr}$, and by $2\,{\rm Gyr}$ the dark-matter profile is well fit by a pseudo-isothermal profile with core radius $r_c\approx 800\,{\rm pc}$ and central density $\rho_0\approx 0.1$–$0.2\,M_\odot\,{\rm pc}^{-3}$. Over the same interval the inner slope evolves from $\gamma\approx -1$ at $t=0$ to $\gamma\approx -0.3\dots -0.5$, while the star-formation history oscillates by factors of $5$–$10$ with duty cycle $\sim 50$–$100\,{\rm Myr}$, and the final stellar component becomes thick and hot with $v_{\rm rot}/\sigma_\theta\sim 1$ [1206.4895]. Ogiya and Mori formulated the same basic process as a resonance problem. If the baryonic potential oscillates with period $T_{\rm osc}$, halo particles are efficiently heated at the radius where the local dynamical time satisfies $$ t_{\rm dyn}(r_{\rm core})\approx T_{\rm osc}. $$ For an NFW halo, their analytic estimate gives $$ r_{\rm core}\simeq 25\, \left(\frac{M_{\rm vir}}{10^9\,M_\odot}\right)^{1/3} \left(\frac{T_{\rm osc}}{10^7\,{\rm yr}}\right)^2 {\rm pc} $$ for $c=10$ and $z=0$. Their $N$-body tests show that runs with $T\lesssim 3\tau$ form a clear core at the radius predicted by $t_{\rm dyn}(r)=T$, and that the oscillation period fixes the eventual core radius, while baryon mass and oscillation amplitude mainly determine the rate of energy transfer and the number of cycles required [1206.5412]. Taken together, these results define the baryonic cusp-to-core transition not as a single explosive event but as a cumulative, non-adiabatic reorganization of collisionless orbits by repeated potential fluctuations. ## 5. Alternative dynamical and non-CDM channels Baryonic supernova feedback is not the only proposed route. On cluster scales, Martizzi, Teyssier, and Moore showed in idealized simulations that repeated AGN-driven outflow/inflow cycles can flatten a cuspy NFW halo. In their setup, a halo with $M_{200}=1.4\times 10^{13}\,M_\odot$ and $c=10$ hosts a central black hole of initial mass $10^9\,M_\odot$ with radiative efficiency $\epsilon_r=0.1$ and coupling $\epsilon_c\approx 0.15$. After $\sim 15$–$20$ AGN cycles over $2$–$3\,{\rm Gyr}$, the dark mass within $r<5\,{\rm kpc}$ drops by $\sim 60\%$, the profile flattens over $r_{\rm core}\approx 5$–$10\,{\rm kpc}$, and the inner slope evolves from $\alpha\approx -1.0$ to $\alpha_{\rm final}\approx -0.7$. However, the observational comparison by Hayashi et al. places groups and clusters close to the cuspy NFW baseline, so cluster-scale core formation remains system-dependent rather than an empirically established generic outcome [1211.2648] [2507.22155]. In dwarf spheroidals, several mechanisms invoke dynamical heating by perturbers rather than gas feedback. Boldrini, Mohayaee, and Silk simulated Fornax with globular clusters embedded in dark-matter minihalos of mass $M_{\rm mh}=2\times 10^7\,M_\odot$. In their recent-accretion scenario, five clusters on eccentric orbits pass through the central $r\lesssim 300\,{\rm pc}$ region 3–4 times within $\simeq 3\,{\rm Gyr}$, with a typical inner orbital period of $\simeq 0.8\,{\rm Gyr}$. The fitted inner slope evolves from $\alpha_{\rm inner}=1.00$ at $t=0$ to $0.12$ at $t=3.5\,{\rm Gyr}$, while the fitted core radius grows from $0$ to $380\,{\rm pc}$ after four crossings. Between passages, phase mixing permits partial cusp regeneration, with $\alpha_{\rm inner}$ recovering by $+0.1$–$0.2$ and $r_c$ shrinking by $\simeq 20$–$30\,{\rm pc}$, yielding a sawtooth sequence of flattening and partial re-steepening [1909.07404]. Related $N$-body calculations of globular-cluster evolution in static and time-varying halo potentials show that clusters in a static cusp remain compact, with $R_{\rm eff}^M\simeq 3$–$6\,{\rm pc}$ by $14\,{\rm Gyr}$, whereas those in a static core expand to $R_{\rm eff}^M\gtrsim 10\,{\rm pc}$; time-varying cusp-to-core histories yield intermediate final sizes around $\simeq 8\,{\rm pc}$ but retain extra tidal debris from the earlier cusp phase [1906.04759]. Such channels are especially relevant for gas-poor dwarf spheroidals, where the stellar-feedback mechanism appears negligible in the classical-dwarf regime and core formation must be sought either in perturber-driven reconfiguration of the potential or in non-standard dark-matter physics [2201.01056]. The review literature groups the non-CDM alternatives into self-interacting dark matter, warm or wave-like dark matter, and other particle-physics mechanisms [2209.14151]. SIDM with cross section per unit mass $\sigma/m\sim 0.1$–$1\,{\rm cm^2\,g^{-1}}$ produces cored profiles with core radii of a few kiloparsecs in halos of $10^9$–$10^{12}\,M_\odot$, while larger $\sigma/m>10\,{\rm cm^2\,g^{-1}}$ is excluded by cluster lensing and ellipticity constraints. Warm dark matter with $m_x\approx 2$–$3\,{\rm keV}$ produces only $\lesssim 100\,{\rm pc}$ cores in dwarfs and does not solve cusp–core and too-big-to-fail simultaneously [2209.14151]. More specialized proposals include primordial-black-hole heating, oscillating asymmetric dark matter, and scalar-field dark matter. In two-component dwarf halos, PBHs of mass $25$–$100\,M_\odot$ comprising at least $1\%$ of the dark matter can heat the cold component through dynamical friction and two-body relaxation, producing robust $10$–$40\,{\rm pc}$ cores in $1$–$8\,{\rm Gyr}$ [1909.07395]. In late-time dark-matter oscillation models, a tiny dark-number-violating mass term $\delta m\sim 10^{-30}$–$10^{-28}\,{\rm eV}$ reactivates annihilation in halo centers; simulations then produce dwarf-galaxy cores of order $0.5$–$2\,{\rm kpc}$ and cluster cores of order $10$–$100\,{\rm kpc}$ for $m_\chi\sim 0.1$–$1\,{\rm GeV}$ [2010.12583]. In an exponential-potential scalar-field model, $$ V(\phi)=V_0 e^{-\alpha\phi}, $$ the self-pressure of the field yields a central slope $\gamma=0$ and an outer $\rho\propto r^{-2}$ envelope, so the halo interpolates analytically between a constant-density core and a flat-rotation-curve regime [1008.3717]. ## 6. Scaling relations, energetics, and unresolved issues A major development in the cusp-to-core literature is the connection between inner-profile transformation and halo scaling relations. Ogiya et al. studied a mapping from an initial NFW halo to a final Burkert halo under two assumptions: conservation of total virial mass $M_{200}$ and preservation of the outer density profile beyond $\simeq r_s$. These imply $$ \rho_s r_s^3=\rho_0 r_0^3. $$ Under this mapping the central surface density $$ \mu_{0D}\equiv \rho_0 r_0 $$ becomes nearly mass-independent, with the observed value $$ \mu_{0D}\simeq 140\,M_\odot\,{\rm pc}^{-2} $$ over more than 8 orders of magnitude in $M_{200}$. For dwarf-galaxy core radii $r_0\sim 0.5$–$1\,{\rm kpc}$, the same relation yields $$ M(<300\,{\rm pc})\approx 10^7\,M_\odot, $$ thereby recovering the Strigari relation as a consequence of constant $\mu_{0D}$ in that regime. Because $\rho_0\approx \rho_s$ and $\rho_s\propto \rho_{\rm crit}(z_{\rm form})\propto (1+z_{\rm form})^3$, the central density of a core also records halo formation redshift [1309.1646]. Kaneda, Mori, and Otaki recast this program in terms of the $\Lambda$CDM concentration–mass relation. They derive an analytic NFW-to-Burkert mapping by enforcing equality of density and enclosed mass at $r_{200}$, solving for $r_c=\eta r_s$ and then recomputing $V_{\rm max}$, $r_{\rm max}$, and characteristic surface densities for the cored analogue. In that framework, the predicted cored and cuspy $\bar\Sigma(