Cold Dark Matter Counterparts Research Articles

X-ray morphology of cluster-mass haloes in self-interacting dark matter

ABSTRACT We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq (1\!-\!3)\times 10^{15}{\, \rm M_\odot }$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\!-\!0.2R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma /m)= 0.5\!-\!1\, {\rm cm}^2\, {\rm g}^{-1}$ than CDM and SIDM with $(\sigma /m)=0.1\, {\rm cm}^2\, {\rm g}^{-1}$. Our statistical analysis indicates that the latter two models are disfavoured at the $68{{\ \rm per\ cent}}$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.

The feasibility of constraining DM interactions with high-redshift observations by JWST

ABSTRACT Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in cold dark matter (CDM) and alternative DM models which have varying relativistic coupling and self-interaction strengths. The interacting DM scenarios produce a cutoff in the linear power spectrum on small-scales, followed by a series of ‘dark acoustic oscillations’. We find that DM interactions suppress the abundance of galaxies below $M_\star \sim 10^8\, {\rm M}_\odot$ for the models considered. The cutoff in the power spectrum delays structure formation relative to CDM. Objects in ETHOS that end up at the same final masses as their CDM counterparts are characterized by a more vigorous phase of early star formation. While galaxies with $M_\star \lesssim 10^6\, {\rm M_\odot }$ make up more than 60 per cent of star formation in CDM at z ≈ 10, they contribute only about half the star formation density in ETHOS. These differences diminish with decreasing redshift. We find that the effects of DM self-interactions are negligible compared to effects of relativistic coupling (i.e. the effective initial conditions for galaxy formation) in all properties of the galaxy population we examine. Finally, we show that the clustering strength of galaxies at high redshifts depends sensitively on DM physics, although these differences are manifest on scales that may be too small to be measurable by JWST.

The central densities of Milky Way-mass galaxies in cold and self-interacting dark matter models

ABSTRACT We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the ‘Feedback In Realistic Environment’ (FIRE) project. The three simulated haloes have virial masses of $\sim 10^{12}\, \text{M}_\odot$ at z = 0, and we study velocity-independent self-interaction cross sections of 1 and 10 ${\rm cm^2 \, g^{-1}}$. We study star formation rates and the shape of dark matter density profiles of the parent haloes in both cold dark matter (CDM) and SIDM models. Galaxies formed in the SIDM haloes have higher star formation rates at z ≤ 1, resulting in more massive galaxies compared to the CDM simulations. While both CDM and SIDM simulations show diverse shape of the dark matter density profiles, the SIDM haloes can reach higher and more steep central densities within few kpcs compared to the CDM haloes. We identify a correlation between the build-up of the stars within the half-mass radii of the galaxies and the growth in the central dark matter densities. The thermalization process in the SIDM haloes is enhanced in the presence of a dense stellar component. Hence, SIDM haloes with highly concentrated baryonic profiles are predicted to have higher central dark matter densities than the CDM haloes. Overall, the SIDM haloes are more responsive to the presence of a massive baryonic distribution than their CDM counterparts.

The onset of gravothermal core collapse in velocity-dependent self-interacting dark matter subhaloes

ABSTRACT It has been proposed that gravothermal collapse due to dark matter self-interactions (i.e. self-interacting dark matter, SIDM) can explain the observed diversity of the Milky Way (MW) satellites’ central dynamical masses. We investigate the process behind this hypothesis using an N-body simulation of a MW-analogue halo with velocity-dependent SIDM (vdSIDM) in which the low-velocity self-scattering cross-section, $\sigma _{\rm T}/m_{\rm x}$, reaches 100 cm2 g−1; we dub this model the vd100 model. We compare the results of this simulation to simulations of the same halo that employ different dark models, including cold dark matter (CDM) and other, less extreme SIDM models. The masses of the vd100 haloes are very similar to their CDM counterparts, but the values of their maximum circular velocities, Vmax, are significantly higher. We determine that these high Vmax subhaloes were objects in the mass range [5 × 106, 1 × 108] M⊙ at z = 1 that undergo gravothermal core collapse. These collapsed haloes have density profiles that are described by single power laws down to the resolution limit of the simulation, and the inner slope of this density profile is approximately −3. Resolving the ever decreasing collapsed region is challenging, and tailored simulations will be required to model the runaway instability accurately at scales <1 kpc.

Self-interacting dark matter and the delay of supermassive black hole growth

ABSTRACT Using cosmological hydrodynamic simulations with physically motivated models of supermassive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass (Mvir ≈ 7 × 1011 M⊙ at z = 0) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm2 g−1. We find that SMBH formation is suppressed in the early Universe due to SIDM interactions. SMBH–SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyr of their evolution. Nonetheless, by z = 0.8 the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the MBH–M* relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.

Self-Interacting Dark Matter Subhalos in the Milky Way’s Tides

We study evolution of self-interacting dark matter subhalos in the MilkyWay tidal field. The interaction between the subhalos and the MilkyWay's tides lead to more diverse dark matter distributions in the inner region, compared to their cold dark matter counterparts. We test this scenario with two MilkyWay satellite galaxies, Draco and Fornax, opposite extremes in the inner dark matter content, and find that they can be accommodated within the self-interacting dark matter model proposed to explain the diverse rotation curves of spiral galaxies in the field.

ETHOS – an effective theory of structure formation: formation of the first haloes and their stars

Abstract A cut-off in the linear matter power spectrum at dwarf galaxy scales has been shown to affect the abundance, formation mechanism and age of dwarf haloes, and their galaxies at high and low redshifts. We use hydrodynamical simulations of galaxy formation within the ETHOS framework in a benchmark model that has such a cut-off and that has been shown to be an alternative to the cold dark matter (CDM) model that alleviates its dwarf-scale challenges. We show how galaxies in this model form differently to CDM, on a halo-by-halo basis, at redshifts z ≥ 6. We show that when CDM haloes with masses around the ETHOS half-mode mass scale are resimulated with the ETHOS matter power spectrum, they form with 50 per cent less mass than their CDM counterparts due to their later formation times, yet they retain more of their gas reservoir due to the different behaviour of gas and dark matter during the monolithic collapse of the first haloes in models with a galactic-scale cut-off. As a result, galaxies in ETHOS haloes near the cut-off scale grow rapidly between z = 10 and 6 and by z = 6 end up having very similar stellar masses, higher gas fractions and higher star formation rates relative to their CDM counterparts. We highlight these differences by making predictions for how the number of galaxies with old stellar populations is suppressed in ETHOS for both z = 6 galaxies and for gas-poor Local Group fossil galaxies. Interestingly, we find an age gradient in ETHOS between galaxies that form in high- and low-density environments.

The edge of galaxy formation III: the effects of warm dark matter on Milky Way satellites and field dwarfs

In this third paper of the series, we investigate the effects of warm dark matter with a particle mass of $m_\mathrm{WDM}=3\,\mathrm{keV}$ on the smallest galaxies in our Universe. We present a sample of 21 hydrodynamical cosmological simulations of dwarf galaxies and 20 simulations of satellite-host galaxy interaction that we performed both in a Cold Dark Matter (CDM) and Warm Dark Matter (WDM) scenario. In the WDM simulations, we observe a higher critical mass for the onset of star formation. Structure growth is delayed in WDM, as a result WDM haloes have a stellar population on average two Gyrs younger than their CDM counterparts. Nevertheless, despite this delayed star formation, CDM and WDM galaxies are both able to reproduce the observed scaling relations for velocity dispersion, stellar mass, size, and metallicity at $z=0$. WDM satellite haloes in a Milky Way mass host are more susceptible to tidal stripping due to their lower concentrations, but their galaxies can even survive longer than the CDM counterparts if they live in a dark matter halo with a steeper central slope. In agreement with our previous CDM satellite study we observe a steepening of the WDM satellites' central dark matter density slope due to stripping. The difference in the average stellar age for satellite galaxies, between CDM and WDM, could be used in the future for disentangling these two models.

Warm FIRE: simulating galaxy formation with resonant sterile neutrino dark matter

We study the impact of a warm dark matter (WDM) cosmology on dwarf galaxy formation through a suite of cosmological hydrodynamical zoom-in simulations of $M_{\rm halo} \approx10^{10}\,M_{\odot}$ dark matter halos as part of the Feedback in Realistic Environments (FIRE) project. A main focus of this paper is to evaluate the combined effects of dark matter physics and stellar feedback on the well-known small-scale issues found in cold dark matter (CDM) models. We find that the $z=0$ stellar mass of a galaxy is strongly correlated with the central density of its host dark matter halo at the time of formation, $z_{\rm f}$, in both CDM and WDM models. WDM halos follow the same $M_{\star}(z=0)-V_{\rm max}(z_{\rm f})$ relation as in CDM, but they form later, are less centrally dense, and therefore contain galaxies that are less massive than their CDM counterparts. As a result, the impact of baryonic effects on the central gravitational potential is typically diminished relative to CDM. However, the combination of delayed formation in WDM and energy input from stellar feedback results in dark matter profiles with lower overall densities. The WDM galaxies studied here have a wider diversity of star formation histories (SFHs) than the same systems simulated in CDM, and the two lowest $M_{\star}$ WDM galaxies form all of their stars at late times. The discovery of young ultra-faint dwarf galaxies with no ancient star formation -- which do not exist in our CDM simulations -- would therefore provide evidence in support of WDM.

Scalar field dark matter: helping or hurting small-scale problems in cosmology?

Building upon results of cosmological simulations of ultra-light scalar field dark matter (SFDM), we present a comprehensive model for the density profiles of SFDM haloes as a function of halo virial mass $M_{\rm h}$ and scalar field mass $m$. The central regions of SFDM haloes are dominated by solitons with characteristic densities that increase with increasing halo mass and asymptote to CDM-like profiles at large radii. For scalar field masses $m \sim 10^{-22}$ eV, consistent with large-scale structure observations, $M_{\rm h} \sim 10^{10} \,M_\odot$ haloes have lower core densities than their Cold Dark Matter (CDM) counterparts and this alleviates the Too Big to Fail problem (TBTF) in a regime where feedback is less effective. However, higher-mass SFDM haloes with $M_{\rm h} \sim 10^{11} \,M_\odot$ are denser than their CDM counterparts at small, observationally relevant radii. We use rotation curves of $V \sim 100$ km s$^{-1}$ galaxies from the SPARC database to show that SFDM exacerbates the cusp/core and central density problems seen in CDM at this scale. We conclude that if the conventional cosmological SFDM scaling relations are correct, then baryonic feedback is required to lower densities in SFDM haloes even more so than in CDM. This motivates cosmological and self-consistent hydrodynamic simulations of SFDM to determine whether central soliton structure can be altered by realistic feedback implementations.

ETHOS - an effective theory of structure formation: predictions for the high-redshift Universe - abundance of galaxies and reionization.

We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36 Mpc3 volume with gas cell mass of ∼ 105 M⊙ and minimum gas softening of ∼ 180 pc) within ETHOS to date - plus a CDM counterpart - to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM τ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep James Webb Space Telescope surveys of strongly lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.

Halo abundance and assembly history with extreme-axion wave dark matter at z ≥ 4

Abstract Wave dark matter (ψDM) composed of extremely light bosons (mψ ∼ 10 − 22 eV), with quantum pressure suppressing structures below a kpc-scale de Broglie wavelength, has become a viable dark matter candidate. Compared to the conventional free-particle ψDM (FPψDM), the extreme-axion ψDM model (EAψDM) proposed by Zhang & Chiueh features a larger cut-off wavenumber and a broad spectral bump in the matter transfer function. Here, we conduct cosmological simulations to compare the halo abundances and assembly histories at z = 4–11 between three different scenarios: FPψDM, EAψDM and cold dark matter (CDM). We show that EAψDM produces significantly more abundant low-mass haloes than FPψDM with the same mψ, and therefore could alleviate the tension in mψ required by the Lyα forest data and by the kpc-scale dwarf galaxy cores. We also find that, compared to the CDM counterparts, massive EAψDM haloes are, on average, 3–4 times more massive at z = 10–11 due to their earlier formation, undergo a slower mass accretion at 7 ≲ z ≲ 11, and then show a rapidly rising major merger rate exceeding CDM by ∼ 50 per cent at 4 ≲ z ≲ 7. This fact suggests that EAψDM haloes may exhibit more prominent starbursts at z ≲ 7.

Substructure and galaxy formation in the Copernicus Complexio warm dark matter simulations

We use the Copernicus Complexio (coco) high-resolution N-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle with a rest mass of 3.3 keV and a sterile neutrino with mass 7 keV and leptogenesis parameter L6 = 8.7. The latter corresponds to the ‘coldest’ model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass ≲ 109 h−1 M⊙ although they predict similar, nearly universal, normalized subhalo radial density distributions. Haloes and subhaloes in both models have cuspy Navarro-Frenk-White profiles, but WDM subhaloes below the cut-off scale in the power spectrum (corresponding to maximum circular velocities Vmaxz = 0 ≤ 50 kms− 1) are less concentrated than their CDM counterparts. We make predictions for observable properties using the galform semi-analytic model of Galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift z > 7 than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.

THE INNER STRUCTURE OF DWARF-SIZED HALOS IN WARM AND COLD DARK MATTER COSMOLOGIES

ABSTRACT By means of N-body + hydrodynamic zoom-in simulations we study the evolution of the inner dark matter and stellar mass distributions of central dwarf galaxies formed in halos of virial masses M v = (2–3) × 1010 h −1 M ⊙ at z = 0, both in a warm dark matter (WDM) and cold dark matter (CDM) cosmology. The half-mode mass in the WDM power spectrum of our simulations is M f = 2 × 1010 h −1 M ⊙. In the dark matter (DM) only simulations halo density profiles are well described by the Navarro–Frenk–White parametric fit in both cosmologies, though the WDM halos have concentrations lower by factors of 1.5–2.0 than their CDM counterparts. In the hydrodynamic simulations, the effects of baryons significantly flatten the inner density, velocity dispersion, and pseudo phase space density profiles of the WDM halos but not of the CDM ones. The density slope, measured at ≈0.02R v, α 0.02, becomes shallow in periods of 2–5 Gyr in the WDM runs. We explore whether this flattening process correlates with the global star formation (SF), M s/M v ratio, gas outflow, and internal specific angular momentum histories. We do not find any clear trends, but when α 0.02 is shallower than −0.5, M s/M v is always between 0.25% and 1%. We conclude that the main reason for the formation of the shallow core is the presence of strong gas mass fluctuations inside the inner halo, which are a consequence of the feedback driven by a very bursty and sustained SF history in shallow gravitational potentials. Our WDM halos, which assemble late and are less concentrated than the CDM ones, obey these conditions. There are also (rare) CDM systems with extended mass assembly histories that obey these conditions and form shallow cores. The dynamical heating and expansion processes behind the DM core flattening apply also to the stars in such a way that the stellar age and metallicity gradients of the dwarfs are softened, their stellar half-mass radii strongly grow with time, and their central surface densities decrease.

The Copernicus Complexio: statistical properties of warm dark matter haloes

The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological $N$-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have $\sim13.4$bn particles, each of mass $\sim 10^5\, h^{-1} M_\odot$, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales ($M_{200}\, \leq\, 2 \times 10^9\,h^{-1}\,M_\odot$) begins slightly later in COCO-Warm than in COCO-Cold. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass $\sim 2 \times 10^{9}\,h^{-1}\,M_\odot$ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass $M_{200}= 10^{8}\,h^{-1}\,M_\odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-Warm rotate slightly more slowly than their CDM counterparts.

SIMULATIONS OF GALAXIES FORMED IN WARM DARK MATTER HALOS OF MASSES AT THE FILTERING SCALE

We present zoom-in N-body + hydrodynamic simulations of dwarf central galaxies formed in warm dark matter (WDM) halos with present-day masses of 2–4 × M⊙. Two different cases are considered: the first one when halo masses are close to the corresponding half-mode filtering scale, Mf ( = 1.2 keV), and the second when they are 20 to 30 times the corresponding Mf ( = 3.0 keV). The WDM simulations are compared with the respective cold dark matter (CDM) simulations. The dwarfs formed in halos of masses (20–30)Mf have roughly similar properties and evolution to their CDM counterparts; on the contrary, those formed in halos of masses around Mf, are systematically different from their CDM counterparts. As compared to the CDM dwarfs, they assemble the dark and stellar masses later, having mass-weighted stellar ages 1.4–4.8 Gyr younger; their circular velocity profiles are shallower, with maximal velocities 20%–60% lower; their stellar distributions are much less centrally concentrated and with larger effective radii, by factors of 1.3–3. The WDM dwarfs at the filtering scale ( = 1.2 keV) have disk-like structures, and end in most cases with higher gas fractions and lower stellar-to-total mass ratios than their CDM counterparts. The late halo assembly, low halo concentrations, and the absence of satellites of the former with respect to the latter are at the basis of the differences.

Massive Milky Way satellites in cold and warm dark matter: dependence on cosmology

We investigate the claim that the largest subhaloes in high resolution dissipationless cold dark matter (CDM) simulations of the Milky Way are dynamically inconsistent with observations of its most luminous satellites. We find that the inconsistency is largely attributable to the large values of \sigma_8 and n_s adopted in the discrepant simulations producing satellites that form too early and therefore are too dense. We find the tension between observations and simulations adopting parameters consistent with WMAP9 is greatly diminished making the satellites a sensitive test of CDM. We find the Via Lactea II halo to be atypical for haloes in a WMAP3 cosmology, a discrepancy that we attribute to its earlier formation epoch than the mean for its mass. We also explore warm dark matter (WDM) cosmologies for 1--4 keV thermal relics. In 1 keV cosmologies subhaloes have circular velocities at kpc scales ~ 60% lower than their CDM counterparts, but are reduced by only 10% in 4 keV cosmologies. Since relic masses < 2-3 keV are ruled out by constraints from the number of Milky Way satellites and Lyman-\alpha\ forest, WDM has a minor effect in reducing the densities of massive satellites. Given the uncertainties on the mass and formation epoch of the Milky Way, the need for reducing the satellite densities with baryonic effects or WDM is alleviated.

Angular Momentum Profiles of Warm Dark Matter Halos

We compare the specific angular momentum profiles of virialized dark halos in cold dark matter (CDM) and warm dark matter (WDM) models, using high-resolution dissipationless simulations. The simulations were initialized using the same set of modes, except on small scales, where the power was suppressed in WDM below the filtering length. Remarkably, WDM as well as CDM halos are well described by the two-parameter angular momentum profile of Bullock and coworkers, even though the halo masses are below the filtering scale of the WDM. Although the best-fit shape parameters change quantitatively for individual halos in the two simulations, we find no systematic variation in profile shapes as a function of the dark matter type. The scatter in shape parameters is significantly smaller for the WDM halos, suggesting that substructure and/or merging history plays a role in producing scatter about the mean angular momentum distribution, but that the average angular momentum profiles of halos originate from larger scale phenomena or a mechanism associated with the virialization process. The known mismatch between the angular momentum distributions of dark halos and disk galaxies is, therefore, present in WDM as well as CDM models. Our WDM halos tend to have a less coherent (more misaligned) angular momentum structure and smaller spin parameters than do their CDM counterparts, although we caution that this result is based on a small number of halos.