Hubbry Logo
Cuspy halo problemCuspy halo problemMain
Open search
Cuspy halo problem
Community hub
Cuspy halo problem
logo
7 pages, 0 posts
0 subscribers
Be the first to start a discussion here.
Be the first to start a discussion here.
Cuspy halo problem
Cuspy halo problem
Cuspy halo problem
View on Wikipedia
from Wikipedia

The cuspy halo problem (also known as the core-cusp problem) is a discrepancy between the inferred dark matter density profiles of low-mass galaxies and the density profiles predicted by cosmological N-body simulations. Nearly all simulations form dark matter halos which have "cuspy" dark matter distributions, with density increasing steeply at small radii, while the rotation curves of most observed dwarf galaxies suggest that they have flat central dark matter density profiles ("cores").[1][2]

Several possible solutions to the core-cusp problem have been proposed. Many recent studies have shown that including baryonic feedback (particularly feedback from supernovae and active galactic nuclei) can "flatten out" the core of a galaxy's dark matter profile, since feedback-driven gas outflows produce a time-varying gravitational potential that transfers energy to the orbits of the collisionless dark matter particles.[3][4] Other works have shown that the core-cusp problem can be solved outside of the most widely accepted Cold Dark Matter (CDM) paradigm: simulations with warm or self-interacting dark matter also produce dark matter cores in low-mass galaxies.[5][6] It is also possible that the distribution of dark matter that minimizes the system energy has a flat central dark matter density profile.[7]

Simulation results

[edit]

According to W.J.G. de Blok "The presence of a cusp in the centers of CDM halos is one of the earliest and strongest results derived from N-body cosmological simulations."[8] Numerical simulations for CDM structure formation predict some structure properties that conflict with astronomical observations.

Observations

[edit]

The discrepancies range from galaxies to clusters of galaxies. "The main one that has attracted a lot of attention is the cuspy halo problem, namely that CDM models predict halos that have a high density core or have an inner profile that is too steep compared to observations."[9]

Potential solutions

[edit]

The conflict between numerical simulations and astronomical observations creates numerical constraints related to the core/cusp problem. Observational constraints on halo concentrations imply the existence of theoretical constraints on cosmological parameters. According to McGaugh, Barker, and de Blok,[10] there might be 3 basic possibilities for interpreting the halo concentration limits stated by them or anyone else:

  1. "CDM halos must have cusps, so the stated limits hold and provide new constraints on cosmological parameters."[11]
  2. "Something (e.g. feedback, modifications of the nature of dark matter) eliminates cusps and thus the constraints on cosmology."[12]
  3. "The picture of halo formation suggested by CDM simulations is wrong."

One approach to solving the cusp-core problem in galactic halos is to consider models that modify the nature of dark matter; theorists have considered warm, fuzzy, self-interacting, and meta-cold dark matter, among other possibilities.[13] One straightforward solution could be that the distribution of dark matter that minimizes the system energy has a flat central dark matter density profile.[7]

See also

[edit]

References

[edit]
Revisions and contributorsEdit on WikipediaRead on Wikipedia

Cuspy halo problem

View on Grokipedia
from Grokipedia
The cuspy halo problem, also known as the core-cusp problem, refers to the discrepancy between the cuspy central density profiles of halos predicted by (CDM) simulations and the cored profiles inferred from observations of dwarf galaxies and low-surface-brightness galaxies. In CDM models, such as the Lambda (ΛCDM) paradigm, N-body simulations predict that dark matter density ρ increases steeply toward the halo center, following profiles like the Navarro-Frenk-White (NFW) form where the logarithmic slope approaches -1 (i.e., ρ ∝ r^{-1}). By contrast, observational data from rotation curves and reveal much flatter central slopes, typically between 0 and -0.5, indicating constant-density cores rather than cusps. This tension emerged prominently in the early 2000s through analyses of gas-rich dwarf galaxies, such as those in the THINGS survey, which showed central slopes around -0.2 ± 0.2, inconsistent with the steeper cusps from pure simulations. Subsequent studies, including those of low-mass galaxies, have reinforced the evidence for cores, with average observed slopes of approximately -0.22 ± 0.08, while simulations yield values closer to -0.8 to -1.4 depending on resolution and profile fits like NFW or Einasto. The problem is particularly acute in low-mass systems, where baryonic effects are expected to be minimal, highlighting potential shortcomings in the collisionless CDM framework. The cuspy halo problem represents one of the key small-scale challenges to the ΛCDM model, alongside issues like the missing satellites and too-big-to-fail problems, and has spurred extensive research into resolution mechanisms. Proposed solutions include baryonic processes, such as supernova feedback that drives gas outflows and flattens cusps, or from massive objects; alternative candidates like warm or self-interacting , which naturally produce cores; and modifications to . Recent high-resolution simulations and observations, including those from 2025, suggest that while some tensions persist—especially in ultra-faint dwarfs—advances in modeling baryonic physics may reconcile predictions with data in many cases.

Background and Definition

Core Concept

The cold dark matter (CDM) paradigm posits that the large-scale structure of the universe arises from the gravitational instability of a smooth, initially uniform distribution of non-relativistic, collisionless dark matter particles, leading to the hierarchical formation of dark matter halos with centrally cuspy density profiles. These profiles feature a density ρ\rho that diverges toward the galactic center, scaling asymptotically as ρ∝r−α\rho \propto r^{-\alpha} with α≈1\alpha \approx 1 in the inner regions. A cusp describes this steep density gradient, where the dark matter density rises sharply as the radius rr approaches zero, in contrast to a core, which is a central region with roughly constant density. The cuspy halo problem emerges from the apparent mismatch between these predicted cuspy profiles and evidence for cored central densities, a discrepancy most evident in low-mass galaxies where the inner dynamics are dominated by dark matter. This theoretical expectation is encapsulated in the Navarro-Frenk-White (NFW) profile, a universal form fitted to halo structures from N-body simulations of CDM collapse: ρ(r)=ρs(r/rs)(1+r/rs)2,\rho(r) = \frac{\rho_s}{(r/r_s)(1 + r/r_s)^2}, where ρs\rho_s is a characteristic density and rsr_s is a scale radius; the profile yields an inner cusp slope of −1-1.

Historical Development

Predictions of cuspy dark matter halo profiles emerged from numerical simulations of structure formation in the cold dark matter (CDM) paradigm in the early 1990s. However, the recognition of the cuspy halo problem, also known as the core-cusp problem, began in 1994 when Flores and Primack, along with Moore, highlighted the tension between these cuspy profiles and observations of cored density distributions in dwarf and low-surface-brightness galaxies. High-resolution N-body simulations conducted by Navarro, Frenk, and White in 1996 revealed that dark matter halos exhibit a universal density profile with a central cusp, where the density ρ scales as ρ ∝ r^{-1} at small radii, as encapsulated in the Navarro-Frenk-White (NFW) profile. This finding solidified the expectation that CDM halos should be cuspy rather than featuring constant-density cores, marking a key milestone in theoretical predictions for dark matter distribution. Initial observational tensions were further detailed in the mid-1990s and early 2000s through analyses of rotation curves in dwarf and low-surface-brightness galaxies, which favored cored profiles over the predicted cusps. For instance, de Blok et al. (2001) derived high-resolution rotation curves from HI observations of such systems, demonstrating that their central densities are better fit by pseudo-isothermal profiles with flat cores (ρ ≈ constant) rather than NFW cusps. These results highlighted a growing discrepancy between simulations and data, prompting further scrutiny of small-scale structure in the CDM model. The explicit naming of the "core-cusp problem" appeared in the literature around , framing the mismatch as a fundamental challenge to CDM predictions on galactic scales. Bosma (2004), in an observational overview, referred to the issue as the core/cusp problem, emphasizing its implications for in low-surface-brightness galaxies. This terminology encapsulated the ongoing debate over whether observations systematically indicate cores where theory demands cusps. Throughout the 2010s, the debate evolved with refined observations and simulations, intensifying scrutiny of the problem's implications for the ΛCDM model. Reviews such as de Blok (2010) synthesized evidence from gas-rich dwarfs showing persistent core preferences, while higher-resolution data from dwarf spheroidals added layers of complexity, revealing that the tension persists even after accounting for baryonic effects. Similarly, Bullock and Boylan-Kolchin (2017) highlighted how the core-cusp issue, alongside other small-scale challenges, tests the robustness of ΛCDM, spurring proposals for resolutions like self-interacting . Despite these advances, the problem remained unresolved, underscoring uncertainties in both theory and observation.

Theoretical Framework

N-Body Simulations

N-body simulations are computational methods that model the gravitational evolution of a large number of particles representing in the , treating them as a collisionless fluid under Newtonian . These simulations discretize the distribution into N discrete particles, each with mass m = ρ L³ / N, where ρ is the mean cosmic and L is the simulation box size, and evolve their positions and velocities over cosmic time using the equations of motion derived from for the . By solving these equations numerically, the simulations track the growth of structure from tiny initial fluctuations to the formation of galaxies and clusters. Key techniques in cosmological N-body simulations include particle-mesh (PM) methods, which compute the on a fixed Cartesian grid using fast Fourier transforms (FFT) to solve efficiently for long-range forces, achieving O(N log N) scaling but limited resolution due to grid discretization. Tree-based algorithms, such as the Barnes-Hut method introduced in 1986, approximate distant particle interactions by grouping them into a hierarchical structure, enabling higher resolution for short-range forces at O(N log N) cost while reducing computational expense compared to direct particle-particle summation. Hybrid approaches like tree-particle-mesh (TPM) combine PM for large-scale modes with tree or particle-particle calculations for small scales, balancing accuracy and efficiency in large-volume runs. In these simulations, halos form through the hierarchical merging of smaller structures originating from primordial fluctuations in the early , seeded by quantum fluctuations during and amplified by gravitational instability as the expands. Initial conditions are generated from linear theory power spectra consistent with observations, with particles placed according to a , and the simulation evolves forward from high (e.g., z ≈ 100) to the present, allowing overdensities to collapse into bound halos via nonlinear gravitational dynamics. High-resolution examples include the Millennium Simulation of 2005, which used the code to follow 10 billion particles in a 500 h⁻¹ Mpc volume under a ΛCDM cosmology, resolving halo structures down to scales. Recent high-resolution simulations, such as those from the IllustrisTNG project (as of 2018 and later refinements), continue to support the predicted cuspy profiles with inner slopes around -1, with improved resolution minimizing artificial relaxation effects. A primary limitation of pure N-body simulations is their assumption of collisionless dark matter dynamics, neglecting the effects of baryonic matter such as gas cooling, , and feedback processes that can alter halo structures in realistic galaxy formation scenarios. Early simulations ignored baryons entirely to focus on clustering, though subsequent hydrodynamical extensions incorporate them to address these shortcomings. These simulations typically predict central density profiles for halos resembling the Navarro-Frenk-White (NFW) form.

Predicted Density Profiles

Simulations of (CDM) halos in a ΛCDM cosmology predict a universal density profile characterized by a cuspy inner region and a steeper outer decline. The Navarro-Frenk-White (NFW) profile, derived from high-resolution N-body simulations, describes this structure with an inner density slope of approximately -1 (where ρ∝r−α\rho \propto r^{-\alpha} and α≈1\alpha \approx 1) transitioning to an outer slope of approximately -3. This form arises from hierarchical , where smaller subhalos merge to build larger halos, leading to a central concentration of dark matter. While the NFW profile provides a simple two-parameter fit (scale radius rsr_s and characteristic density ρs\rho_s), variations in halo properties introduce deviations. An alternative, the Einasto profile, offers a better fit across a wider range of radii, given by ρ(r)=ρsexp⁡[−2α((rrs)α−1)],\rho(r) = \rho_s \exp\left[-\frac{2}{\alpha} \left( \left(\frac{r}{r_s}\right)^\alpha - 1 \right)\right], with the shape parameter α≈0.17\alpha \approx 0.17 for Milky Way-mass halos. This profile captures a smoother transition without sharp breaks, reflecting the cumulative effects of multiple accretion episodes in halo assembly. Simulations predict an inner slope close to -1 that is nearly independent of halo mass and formation redshift, consistent with the universal NFW profile. While the concentration parameter shows dependence on mass and accretion history, the central cusp slope remains around α ≈ 1 across a wide range of halo masses. Numerical convergence remains a challenge for resolving the innermost regions. Early simulations suggested steeper inner slopes (α≈1.5\alpha \approx 1.5), but higher-resolution runs with more particles reveal a convergence toward α≈1\alpha \approx 1 beyond a minimum resolved radius proportional to the particle mass. This resolution dependence highlights the need for simulations with at least 10610^6 particles per halo to reliably predict cuspy profiles, as artificial two-body relaxation can flatten or steepen the core artificially.

Observational Constraints

Evidence from Dwarf Galaxies

Dwarf spheroidal galaxies, such as Draco and Sculptor, serve as ideal probes of distributions due to their high mass-to-light ratios, often exceeding 100 in solar units, indicating dominance by over baryonic components throughout their extent. Kinematic studies of these systems rely on measurements of stellar dispersions, which reveal flat or slowly varying profiles in the central regions, consistent with constant densities rather than steeply rising ones. In the seminal analysis by Walker et al. (2009), application of the Jeans equation to data from eight classical dwarf spheroidals, including Draco and Sculptor, yielded mass profiles consistent with cored dark matter halos with scale radii around 150 pc, as well as cuspy Navarro–Frenk–White (NFW) profiles with larger scale radii (~800 pc). Subsequent analyses have reinforced these findings, estimating core radii in the range of 280 pc to 1.3 kpc for several satellites by modeling velocity dispersion profiles with isothermal or similar distributions. For instance, in Draco and Sculptor, the inferred central densities remain roughly constant out to these scales, supporting the presence of sizable cores. More recent studies using proper motions, such as those from 2022, have tightened constraints, finding evidence for cuspy profiles in massive dwarf spheroidals while cores remain favored in lower-mass systems. Complementary evidence comes from neutral hydrogen (HI) observations of gas-rich dwarf irregular galaxies, where rotation curves rise slowly in the inner regions, implying flat dark matter density profiles with core sizes on the order of hundreds of parsecs. These kinematic signatures from both stellar and gas tracers contrast with the cuspy central densities predicted by cold dark matter simulations.

Rotation Curve Measurements

Rotation curves, which measure the orbital velocities v(r)v(r) of gas or stars as a function of galactocentric radius rr, provide key observational constraints on the inner dark matter density profiles of galaxies. These velocities are typically derived from the 21 cm emission line of neutral hydrogen (H I) using radio interferometry, allowing mapping of the velocity fields in the gaseous disks of spiral, low-surface-brightness (LSB), and dwarf irregular galaxies. In the context of the cuspy halo problem, rotation curves reveal discrepancies between predicted cuspy profiles from cold dark matter (CDM) simulations and observed behaviors suggesting cored distributions. High-resolution observations, such as those from The H I Nearby Galaxy Survey (THINGS) conducted with the , have targeted LSB and dwarf irregular galaxies to probe the central regions where dominance is expected. The THINGS survey, analyzing 19 such galaxies with spatial resolutions down to sub-kpc scales and velocity resolutions of about 5 km/s, finds that inner rotation curves often exhibit a shallow, nearly linear rise in v(r)v(r) at small radii, rather than the steeper r\sqrt{r}
Add your contribution
Related Hubs
User Avatar
No comments yet.