Core-cusp Problem Research Articles

Galaxy Formation in ΛCDM Cosmology

This is a golden age for galaxy formation: Existing and especially new telescopes are providing observations that challenge and illuminate rapidly improving theory and simulations. This review describes the formation of the cosmic web and the structure of the dark matter halos that provide the scaffolding of the Universe. It then summarizes how empirical models, semianalytic models, and hydrodynamic simulations attempt to account for key properties of the galaxy population, including the main sequence of star-forming galaxies, the inefficiency of star formation, the shape evolution and color bimodality of galaxies, and the phenomena that cause galaxies to quench their star formation. It concludes with a summary of observations that have challenged the cosmological constant cold dark matter (ΛCDM) paradigm of galaxy formation—including the Hubble and S 8 tensions, bright galaxies in the early Universe, an extragalactic background light mystery, missing satellite galaxies, the diversity of dwarf galaxies, the cusp–core problem, the too-big-to-fail problem, stellar clumps, planes of satellite galaxies, and galaxies without dark matter—and solutions that have been proposed.

Bigravity portal dark matter

We consider a model where the interaction between dark matter and the Standard Model particles is mediated by a ghost-free bigravity portal. The bigravity model invokes a massive spin-2 particle coupled to the usual massless graviton as well as generic bimetric matter couplings. The cross sections for dark matter direct detection are computed and confronted with the experimental bounds. The presence of the massive spin-2 mediator resolves the core-cusp problem, which in turn significantly constrains the dark matter coupling in such a bigravity theory. Yet, there remains a window of the parameter space where the model can be tested in the upcoming direct detection experiments such as XENONnT and PandaX-30T. The model also predicts a reheating temperature on the order of 106 GeV. Published by the American Physical Society 2024

Structure of the equivalent Newtonian systems in MONDN-body simulations

Aims.We investigate the core-cusp problem of the Λ cold dark matter (ΛCDM) scenario in the context of the modified Newtonian dynamics (MOND) paradigm while exploiting the concept of an equivalent Newtonian system (ENS).Methods.By means of particle-meshN-body simulations in MOND, we explored the processes of galaxy formation via cold dissipationless collapse and the merging of smaller substructures. From the end states of our simulations, we recovered the associated ENS and studied the properties of their dark matter halos. We compared the simulation results with simple analytical estimates with a family ofγ-models.Results.We find that the dark matter density of ENSs of most spherical cold collapses have a markedly cored structure, particularly for the lowest values of the initial virial ratios. End states of some simulations with initially clumpy conditions have more complex profiles, and some of their ENSs exhibit a moderate cusp, with the logarithmic density slope always shallower than one.Conclusions.In contrast to what one would expect from theoretical and numerical arguments in ΛCDM, these results seem to point towards the fact that the absence of a central DM cusp in most observed galaxies would be totally consistent in a MONDian description.

Aspherical ULDM collapse: variation in the core–halo mass relation

ABSTRACT Ultralight dark matter (ULDM) is an interesting alternative to the cold dark matter (CDM) paradigm. Due to the extremely low mass of the constituent particle (∼10−22 eV), ULDM can exhibit quantum effects up to kiloparsec scales. In particular, runaway collapse in the centres of ULDM haloes is prevented by quantum pressure, providing a possible resolution to the ‘core-cusp problem’ of CDM. However, the the detailed relationship between the ULDM core mass and that of the overall halo is poorly understood. We simulate the collapse of both spherical and aspherical isolated ULDM overdensities using AxioNyx, finding that the central cores of collapsed haloes undergo sustained oscillatory behaviour, which affects both their peak density and overall morphology. The asphericity of the post-collapse core fluctuates both quantitatively and qualitatively, with oblate initial overdensities generating cores, which fluctuate between prolate and oblate configurations, and more complicated dynamics arising in triaxial scenarios. Furthermore, the peak central densities are higher in spherical configurations. Consequently, astrophysically realistic haloes may exhibit substantial departures from theoretical core–halo profiles and we would expect a significant variance of the properties of haloes with the same mass.

Structure formation paradigm and axion quark nugget dark matter model

We advocate an idea that ``non-baryonic" dark matter in form of nuggets made of standard model quarks and gluons (similar to the old idea of the Witten's strangelets) could play a crucial role in structure formation. The corresponding macroscopically large objects, which are called the axion quark nuggets (AQN) behave as {\it chameleons}: they do not interact with the surrounding material in dilute environment, but they become strongly interacting objects in sufficiently dense environment. The AQN model was invented long ago with a single motivation to explain the observed similarity $\Omega_{\rm DM}\sim \Omega_{\rm visible}$ between visible and DM components. This relation represents a very generic feature of this framework, not sensitive to any parameters of the construction. We argue that the strong visible-DM interaction may dramatically modify the conventional structure formation pattern at small scales to contribute to a resolution of a variety of interconnected problems (such as Core-Cusp problem, etc) which have been a matter of active research and debates in recent years. We also argue that the same visible-DM interaction at small scales is always accompanied by a broad band diffuse radiation. We speculate that the recently observed excesses of the UV emission by JWST at high redshifts and by GALEX in our own galaxy might be a direct manifestation of this AQN-induced radiation. We also speculate that the very same source of energy injection could contribute to the resolution of another long standing problem related to the Extragalactic Background Light (EBL) with known discrepancies in many frequency bands (from UV to optical, IR light and radio emissions).

A global analysis of resonance-enhanced light scalar dark matter

We study a minimal model for a light scalar dark matter, requiring a light scalar mediator to address the core-cusp problem and interact with the standard model particles. We analyze the model comprehensively by focusing on the Breit-Wigner resonance for dark matter annihilation and self-scattering channels, considering the thermal relic abundance condition that includes the early kinetic decoupling effect, as well as the present and future constraints from collider, direct, and indirect dark matter detections. We found that the scalar dark matter with the mass of 0.3–2 GeV remains uncharted, which will be efficiently tested by the near future MeV gamma-ray observations.

Is the Core-cusp Problem a Matter of Perspective? Jeans Anisotropic Modeling against Numerical Simulations

Mock member stars for 28 dwarf galaxies are constructed from the cosmological auriga simulation, which reflects the dynamical status of realistic stellar tracers. Axisymmetric Jeans Anisotropic Multi-Gaussian Expansion (jam) modeling is applied to 6000 star particles for each system to recover the underlying matter distribution. The stellar or dark matter component individually is poorly recovered, but the total profile is constrained more reasonably. The mass within the half-mass radius of tracers is recovered the tightest, and the mass between 200 and 300 pc, M(200–300 pc), is an unbiasedly constrained ensemble, with a scatter of 0.167 dex. If using 2000 particles and only line-of-sight velocities with typical errors, the scatter in M(200–300 pc) is increased by ∼50%. Quiescent Saggitarius dSph–like systems and star-forming systems with strong outflows show distinct features, with M(200–300 pc) mostly underestimated for the former, and likely overestimated for the latter. The biases correlate with the dynamical status, which is a result of contraction motions due to tidal effects in quiescent systems or galactic winds in star-forming systems, driving them out of equilibrium. After including Gaia DR3 proper motion errors, we find proper motions can be as useful as line-of-sight velocities for nearby systems at < ∼60 kpc. By extrapolating the actual density profiles and the dynamical constraints down to scales below the resolution, we find the mass within 150 pc can be an unbiasedly constrained ensemble, with a scatter of ∼0.255 dex. We show that the contraction of member stars in nearby systems is detectable based on Gaia DR3 proper motion errors.

The quantum character of the Scalar Field Dark Matter

ABSTRACT The scalar field dark matter (SFDM) model, also called Fuzzy, Wave, Bose–Einstein, and Ultra-light Dark Matter, has received a lot of attention because it has been able to provide simpler and more natural explanations for various features of galaxies, such as the number of satellite galaxies and the cusp-core problem. We recently showed that this model is able to explain the vast polar orbits of satellite galaxies around their host, the so-called VPO, and to explain the X-ray and gamma-ray emissions in the vacuum regions of our galaxy, that is, the Fermi Bubbles. In all these phenomena, the quantum character of SFDM has been crucial. In this work, we study the quantum effects of SFDM at the cosmological level, to see these effects not only at the galactic scale, but also at the cosmological scale. Using a convenient ansatz, we were able to integrate the perturbed equations to show that the shape of the SFDM haloes resembling atoms is a generic result. The main conclusion of this work is that quantum mechanics, the successful microworld theory, could also explain the dark side of the Cosmos.

Constraints on dark matter self-interaction from the internal density profiles of X-COP galaxy clusters

The fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter α ≈ 0.18 in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if DM particles self-interact via additional forces, we expect the mass density profiles of DM halos to flatten in their central regions, thereby increasing the Einasto shape parameter. We measured the structural properties of 12 massive galaxy clusters from observations of their hot gaseous atmosphere, using the X-ray observatory XMM-Newton, and of the Sunyaev-Zel’dovich effect using the Planck all-sky survey. After removing morphologically disturbed systems, we measured Einasto shape parameters with mean ⟨α⟩=0.19 ± 0.03 and intrinsic scatter σα = 0.06, which is in close agreement with the prediction of the cold dark matter paradigm. We used cosmological hydrodynamical simulations of cluster formation with self-interacting DM (BAHAMAS-SIDM) to determine how the Einasto shape parameter depends on the self-interaction cross section. We used the fitted relation to turn our measurements of α into constraints on the self-interaction cross section, which imply σ/m &lt; 0.19 cm2 g−1 (95% confidence level) at collision velocity vDM − DM ∼ 1000 km s−1. This is lower than the interaction cross section required for DM self-interactions to solve the core-cusp problem in dwarf spheroidal galaxies, unless the cross section is a strong function of velocity.

Predictions from the logotropic model: The universal surface density of dark matter halos and the present proportions of dark matter and dark energy

The logotropic model (Chavanis, 2015) may be an interesting alternative to the ΛCDM model. It is able to account for the present accelerating expansion of the universe while solving at the same time the core–cusp problem of the CDM model. In the logotropic model, there is a single dark fluid. Its rest-mass plays the role of dark matter and its internal energy plays the role of dark energy. We highlight two remarkable predictions of the logotropic model. It yields cored dark matter halos with a universal surface density equal to Σ0th=0.01955cΛ/G=133M⊙/pc2, without free parameter, in very good agreement with the observational value Σ0obs=141−52+83M⊙/pc2. It also predicts the present ratio of dark energy and dark matter to be the Euler number Ωde,0th/Ωdm,0th=e=2.71828... in very good agreement with the observations giving Ωde,0obs/Ωdm,0obs=2.669±0.08. Using the measured present proportion of baryonic matter Ωb,0obs=0.0486±0.0010, we find that the values of the present proportions of dark matter and dark energy are Ωdm,0th=11+e(1−Ωb,0)=0.2559 and Ωde,0th=e1+e(1−Ωb,0)=0.6955 in very good agreement with the observational values Ωdm,0obs=0.2589±0.0057 and Ωde,0obs=0.6911±0.0062 within the error bars. These theoretical predictions are obtained by advocating a mysterious strong cosmic coincidence (dubbed “dark magic”) implying that our epoch plays a particular role in the history of the universe. We review the three types of logotropic models introduced in our previous papers depending on whether the equation of state is expressed in terms of the energy density, the rest-mass density, or the pseudo-rest mass density of a complex scalar field. We discuss the similarities and the differences between these models. Finally, we point out some intrinsic difficulties with the logotropic model similar to those encountered by the Chaplygin gas model and discuss possible solutions.

Incorporating baryon-driven contraction of dark matter halos in rotation curve fits

The condensation of baryons within a dark matter (DM) halo during galaxy formation should result in some contraction of the halo as the combined system settles into equilibrium. We quantify this effect on the cuspy primordial halos predicted by DM-only simulations for the baryon distributions observed in the galaxies of the SPARC database. We find that the DM halos of high surface brightness galaxies (with Σeff ≳ 100 L⊙ pc−2 at 3.6 μm) experience strong contraction. Halos become more cuspy as a result of compression: the inner DM density slope increases with the baryonic surface mass density. We iteratively fit rotation curves to find the balance between initial halo parameters (constrained by abundance matching), compression, and stellar mass-to-light ratio. The resulting fits often require lower stellar masses than expected for stellar populations, particularly in galaxies with bulges: stellar mass must be reduced to make room for the DM it compresses. This trade off between dark and luminous mass is reminiscent of the cusp-core problem in dwarf galaxies, but occurs in more massive systems: the present-epoch DM halos cannot follow from cuspy primordial halos unless (1) the stellar mass-to-light ratios are systematically smaller than expected from standard stellar population synthesis models, and/or (2) there is a net outward mass redistribution from the initial cusp, even in massive galaxies widely considered to be immune from such effects.

Predictive model of fermionic dark matter halos with a quantum core and an isothermal atmosphere

We develop a thermodynamical model of fermionic dark matter halos at finite temperature. Statistical equilibrium states may be justified by a process of violent collisionless relaxation in the sense of Lynden-Bell or from a collisional relaxation of nongravitational origin if the fermions are self-interacting. The most probable state (maximum entropy state) generically has a "core-halo" structure with a quantum core (fermion ball) surrounded by an isothermal atmosphere. The quantum core is equivalent to a polytrope of index $n=3/2$. The Pauli exclusion principle creates a quantum pressure that prevents gravitational collapse and solves the core-cusp problem of the cold dark matter model. The isothermal atmosphere (which is similar to the NFW profile of cold dark matter) accounts for the flat rotation curves of the galaxies at large distances. We numerically solve the equation of hydrostatic equilibrium with the Fermi-Dirac equation of state and determine the density profiles and rotation curves of fermionic dark matter halos.

Orbits and adiabatic contraction in scalar-field dark matter halos: revisiting the cusp-core problem in dwarf galaxies

ABSTRACT Bose–Einstein-condensed dark matter, also called scalar field dark matter (SFDM), has become a popular alternative to cold dark matter (CDM), because it predicts galactic cores, in contrast to the cusps of CDM halos (‘cusp-core problem’). We continue the study of SFDM with a strong, repulsive self-interaction; the Thomas–Fermi (TF) regime of SFDM (SFDM-TF). In this model, structure formation is suppressed below a scale related to the TF radius RTF, which is close to the radius of central cores in these halos. We investigate for the first time the impact of baryons onto realistic galactic SFDM-TF halo profiles by studying the process of adiabatic contraction (AC) in such halos. In doing so, we first analyse the underlying quantum Hamilton–Jacobi framework appropriate for SFDM and calculate dark matter orbits, in order to verify the validity of the assumptions usually required for AC. Then, we calculate the impact of AC onto SFDM-TF halos of mass $\sim 10^{11}\, {\rm M}_{\odot }$, with various baryon fractions and core radii, RTF ∼ (0.1–4) kpc, and compare our results with observational velocity data of dwarf galaxies. We find that AC-modified SFDM-TF halos with kpc-size core radii reproduce the data well, suggesting stellar feedback may not be required. On the other hand, halos with sub-kpc core radii face the same issue than CDM, in that they are not in accordance with galaxy data in the central halo parts.

Cuspy dark matter density profiles in massive dwarf galaxies

ABSTRACT Rotation curves of galaxies probe their total mass distributions, including dark matter. Dwarf galaxies are excellent systems to investigate the dark matter density distribution, as they tend to have larger fractions of dark matter compared to higher mass systems. The core-cusp problem describes the discrepancy found in the slope of the dark matter density profile in the centres of galaxies (β*) between observations of dwarf galaxies (shallower cores) and dark matter-only simulations (steeper cusps). We investigate β* in six nearby spiral dwarf galaxies for which high-resolution CO J = 1–0 data were obtained with ALMA (Atacama Large Millimeter/submillimeter Array). We derive rotation curves and decompose the mass profile of the dark matter using our CO rotation curves as a tracer of the total potential and 4.5 $\mu$m photometry to define the stellar mass distribution. We find 〈β*〉 = 0.6 with a standard deviation of ±0.1 among the galaxies in this sample, in agreement with previous measurements in this mass range. The galaxies studied are on the high stellar mass end of dwarf galaxies and have cuspier profiles than lower mass dwarfs, in agreement with other observations. When the same definition of the slope is used, we observe steeper slopes than predicted by the FIRE and NIHAO simulations. This may signal that these relatively massive dwarfs underwent stronger gas inflows towards their centres than predicted by these simulations, that these simulations overpredict the frequency of accretion or feedback events, or that a combination of these or other effects are at work.

A heuristic wave equation parameterizing BEC dark matter halos with a quantum core and an isothermal atmosphere

The Gross–Pitaevskii–Poisson equations that govern the evolution of self-gravitating Bose–Einstein condensates, possibly representing dark matter halos, experience a process of gravitational cooling and violent relaxation. We propose a heuristic parametrization of this complicated process in the spirit of Lynden-Bell’s theory of violent relaxation for collisionless stellar systems. We derive a generalized wave equation that was introduced phenomenologically in Chavanis (Eur Phys J Plus 132:248, 2017) involving a logarithmic nonlinearity associated with an effective temperature $$T_\mathrm{eff}$$ and a damping term associated with a friction $$\xi $$ . These terms can be obtained from a maximum entropy production principle and are linked by a form of Einstein relation expressing the fluctuation-dissipation theorem. The wave equation satisfies an H-theorem for the Lynden-Bell entropy and relaxes towards a stable equilibrium state which is a maximum of entropy at fixed mass and energy. This equilibrium state represents the most probable state of a Bose–Einstein condensate dark matter halo. It generically has a core-halo structure. The quantum core prevents gravitational collapse and may solve the core-cusp problem. The isothermal halo leads to flat rotation curves in agreement with the observations. These results are consistent with the phenomenology of dark matter halos. Furthermore, as shown in a previous paper (Chavanis in Phys Rev D 100:123506, 2019), the maximization of entropy with respect to the core mass at fixed total mass and total energy determines a core mass–halo mass relation which agrees with the relation obtained in direct numerical simulations. We stress the importance of using a microcanonical description instead of a canonical one. We also explain how our formalism can be applied to the case of fermionic dark matter halos.

Constraints on self-interacting Bose-Einstein condensate dark matter using large-scale observables

Constraints on the cosmic history of self-interacting Bose-Einstein condensed (SIBEC) dark matter (DM) are obtained using the cosmic microwave background (CMB), baryonic acoustic oscillations (BAO), growth factor measurements, and type Ia supernovae (SNIa) distances. Four scenarios are considered, one with purely SIBEC-DM, and three in which SIBEC-DM is the final product of some transition from different initial states, which are either cold, warm, or has a constant equation of state. Using a fluid approximation for the self-interacting scalar field it is found that in the simplest scenario of purely SIBEC-DM the self-interaction necessary for solving the cusp-core problem, with core-radii of low-mass halos of order R c ≳ 1kpc, is excluded at 2.4σ, or 98.5% confidence. Introducing a transition, however, relaxes this constraint, but the transitions are preferred to be after matter-radiation equality, and the initial phase to be cold.

The Cusp–Core Problem in Gas-Poor Dwarf Spheroidal Galaxies

This review deals with the inconsistency of inner dark matter density profiles in dwarf galaxies, known as the cusp–core problem. In particular, we aim to focus on gas-poor dwarf galaxies. One of the most promising solutions to this cold dark matter small-scale issue is the stellar feedback, but it seems to be only designed for gas-rich dwarfs. However, in the regime of classical dwarfs, this core mechanism becomes negligible. Therefore, it is required to find solutions without invoking these baryonic processes as dark matter cores tend to persist even for these dwarfs, which are rather dark-matter-dominated. Here, we have presented two categories of solutions. One consists of creating dark matter cores from cusps within cold dark matter by altering the dark matter potential via perturbers. The second category gathers solutions that depict the natural emergence of dark matter cores in alternative theories. Given the wide variety of solutions, it becomes necessary to identify which mechanism dominates in the central region of galaxies by finding observational signatures left by them in order to highlight the true nature of dark matter.

Review of Solutions to the Cusp-Core Problem of the ΛCDM Model

This review aims at proposing to the field an overview of the Cusp-core problem, including a discussion of its advocated solutions, assessing how each can satisfactorily provide a description of central densities. Whether the Cusp-core problem reflects our insufficient grasp on the nature of dark matter, of gravity, on the impact of baryonic interactions with dark matter at those scales, as included in semi-analytical models or fully numerical codes, the solutions to it can point either to the need for a paradigm change in cosmology, or to to our lack of success in ironing out the finer details of the ΛCDM paradigm.

Reshaping our understanding on structure formation with the quantum nature of the dark matter

We study the nonlinear structure formation in cosmology accounting for the quantum nature of the dark matter (DM) particles in the initial conditions at decoupling, as well as in the relaxation and stability of the DM halos. Different from cosmological N-body simulations, we use a thermodynamic approach for collisionless systems of self-gravitating fermions in general relativity, in which the halos reach the steady state by maximizing a coarse-grained entropy. We show the ability of this approach to provide answers to crucial open problems in cosmology, among others: the mass and nature of the DM particle, the formation and nature of supermassive black holes in the early Universe, the nature of the intermediate mass black holes in small halos, and the core-cusp problem.

Thermal relic of self-interacting dark matter with retarded decay of mediator

The existence of a light mediator is beneficial to some phenomena in astroparticle physics, such as the core-cusp problem and diversity problem. It can decouple from Standard Model to avoid direct detection constraints, generally realized by retard decay of the mediator. Their out-of-equilibrium decay process changes the dark matter (DM) freeze-out via temperature discrepancy. This type of hidden sector (HS) typically requires a precision calculation of the freeze-out process considering HS temperature evolution and the thermal average of the cross-section. If the mediator is light sufficiently, we can not ignore the s-wave radiative bound state formation process from the perspective of CMB ionization and Sommerfeld enhancement. We put large mass splitting between DM and mediator, different temperature evolution on the same theoretical footing, discussing the implication for DM relic density in this HS. We study this model and illustrate its property by considering the general Higgs-portal dark matter scenario, which includes all the relevant constraints and signals. It shows that the combination of BBN and CMB constraint favors the not-too-hot HS, rinf< 102, for the positive cubic interaction of mediator scenario. On the other hand, the negative cubic interaction is ruled out except for our proposed blind spot scenario.