Warm Dark Matter Halos Research Articles

Halo assembly in cold and warm dark matter during the JWST frontier epoch

ABSTRACT The JWST mission is in the process of probing the galaxy mass function at z > 10, when conceivably any delay in halo assembly due to the presence of a dwarf galaxy-scale power spectrum cutoff may drastically suppress the number of galaxies relative to the cold dark matter (CDM) expectation. We employ N-body simulations of CDM and warm dark matter (WDM) to explore how the difference in halo collapse time between these models scales with z = 0 descendant halo mass. We demonstrate that collapse begins first for the most massive haloes, and the delay in collapse time between CDM and WDM haloes correlates inversely with descendant mass. We thus infer that only present-day dwarf galaxies exhibit any difference in their assembly history between CDM and WDM at z = 10, and therefore support previous studies that have found JWST is unlikely to determine whether our Universe is better described by the CDM cosmology or the WDM cosmology.

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.

New Model of Density Distribution for Fermionic Dark Matter Halos

We formulate a new model of density distribution for halos made of warm dark matter (WDM) particles. The model is described by a single microphysical parameter – the mass (or, equivalently, the maximal value of the initial phase-space density distribution) of dark matter particles. Given the WDM particle mass and the parameters of a dark matter density profile at the halo periphery, this model predicts the inner density profile. In the case of initial Fermi–Dirac distribution, we successfully reproduce cored dark matter profiles from N-body simulations. We calculate also the core radii of warm dark matter halos of dwarf spheroidal galaxies for particle masses mFD = 100, 200, 300, and 400 eV.

The mass–concentration–redshift relation of cold and warm dark matter haloes

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.

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.

Fingerprints of the initial conditions on the density profiles of cold and warm dark matter haloes

We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radius. Virialization isotropizes the velocity dispersion in the inner regions of the halo but does not erase the memory of the initial conditions in phase space. The location of the observed deviations from CDM in the density profile and in phase space can be directly related to the ratio between the halo mass and the filtering mass and are most evident in small mass haloes, even for a 34 keV thermal relic WDM. The rearrangement of mass within the haloes supports analytic models of halo structure that include angular momentum. We also find evidence of a dependence of the slope of the inner density profile in CDM cosmologies on the halo mass with more massive haloes exhibiting steeper profiles, in agreement with the model predictions and with previous simulation results. Our work complements recent studies of microhaloes near the filtering scale in CDM and strongly argue against a universal shape for the density profile.

Hollow cores in warm dark matter halos from the Vlasov-Poisson equation

We report the results of extended high--resolution numerical integrations of the Vlasov--Poisson equation for the collapse of spherically symmetric WDM halos. For thermal relics with mass $m=1\,$keV/$c^2$, we find collapsed halos with cores of size from $0.1$ to $0.6\,$kpc. The typical core is hollow, with the mass density decreasing towards the core center by almost three orders of magnitude from its maximum near the core radius $r_{\rm c}$. The core is in equilibrium with the diffused part of the halo but far from virialization. These properties are rooted in the conservation of the squared angular momentum and in the original excess, proper of WDM initial conditions, of kinetic energy in the core region. In a sample of more than one hundred simulated collapses, the values of $r_{\rm c}$ and of the core density $\rho_{\rm c}$ are in the range typical of dwarf spheroids, while the maximal circular velocities $V_{\rm max}$ are proper of small disk galaxies. The product $\mu_{\rm c}=\rho_{\rm c}r_{\rm c}$ takes values between $116\,M_\odot/$pc$^2$ and $283\,M_\odot/$pc$^2$, while the surface density $\mu_0$, as determined from a Burkert fit, is roughly three times larger. From these data and data obtained at smaller values of $m$, we extrapolate for one particular halo $\mu_{\rm c}=263(308)\,M_\odot/$pc$^2$ and $\mu_0=754(855)\,M_\odot/$pc$^2$ at $m=2(3.3)\,$keV/$c^2$, to be compared with the observed value $140^{+83}_{-52} M_\odot$/pc$^2$. In view of the many improvements and enhancements available, we conclude that WDM is a viable solution for explaining the presence and the size of cores in low mass galaxies.

The first billion years of a warm dark matter universe

We present results of cosmological N-body hydrodynamic chemistry simulations of primordial structure growth and evolution in a scenario with warm dark matter (WDM) having a mass of 3 keV (thermal relic) and compare with a model consisting of standard cold dark matter (CDM). We focus on the high-redshift universe ($z>6$), where the structure formation process should better reflect the primordial (linear) differences in terms of matter power spectrum. We find that early epochs can be exceptional probes of the dark-matter nature. Non-linear WDM power spectra and mass functions are up to 2 dex lower than in CDM and show spreads of factor of a few persisting in the whole first Gyr. Runaway molecular cooling in WDM haloes results severely inhibited because of the damping of power at large $k$ modes and hence cosmic (population III and II-I) star formation rate (SFR) is usually suppressed with respect to CDM predictions. Luminous objects formed in a WDM background are very rare at $z>10$, due to the sparser and retarded evolution of early WDM mini-haloes during the dark ages and their lack can be fitted with a simple analytical formula depending only on magnitude and redshift. Future high-$z$ observations of faint galaxies have the potential to discriminate between CDM and WDM scenarios by means of cosmic stellar mass density (SMD) and specific SFR, as well. When compared to the effects of alternative cosmologies (e.g. non-Gaussian or dark-energy models) or of high-order corrections at large $z$ (e.g. primordial streaming motions or changes in the pristine IMF) the ones caused by WDM are definitely more dramatic.

The properties of warm dark matter haloes

Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmological N-body simulations of galactic haloes of mass similar to the Milky Way's. On scales larger than the free-streaming cutoff, the initial conditions have the same power spectrum and phases as one of the cold dark matter (CDM) haloes previously simulated by Springel et al as part of the Virgo consortium Aquarius project. We have simulated four haloes with WDM particle masses in the range 1.4-2.3keV and, for one case, we have carried out further simulations at varying resolution. N-body simulations in which the power spectrum cutoff is resolved are known to undergo artificial fragmentation in filaments producing spurious clumps which, for small masses (<10^7Msun in our case) outnumber genuine haloes. We have developed a robust algorithm to identify these spurious objects and remove them from our halo catalogues. We find that the WDM subhalo mass function is suppressed by well over an order magnitude relative to the CDM case for masses <10^9Msun. Requiring that there should be at least as many subhaloes as there are observed satellites in the Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density distributions that are well described by NFW or Einasto profiles. Their central densities are lower for lower WDM particle masses and none of the models we have considered suffer from the "too big to fail" problem recently highlighted by Boylan-Kolchin et al.

Typical density profile for warm dark matter haloes

Abstract Using the model for (bottom-up) hierarchical halo growth recently developed by Salvador-Solé et al., we derive the typical spherically averaged density profile for haloes with several relevant masses in the concordant warm dark matter (ΛWDM) cosmology with non-thermal sterile neutrinos of two different masses. The predicted density profiles become flat at small radii, as expected from the effects of the spectrum cut-off. The core cannot be resolved, however, because the non-null particle velocity yields the fragmentation of minimum mass protohaloes in small nodes, which invalidates the model at the corresponding radii.

Cores and cusps in warm dark matter halos

The apparent presence of large core radii in Low Surface Brightness galaxies has been claimed as evidence in favor of warm dark matter. Here we show that WDM halos do not have cores that are large fractions of the halo size: typically, rcore/r200≲10−3. This suggests an astrophysical origin for the large cores observed in these galaxies, as has been argued by other authors.

Tidal debris of dwarf spheroidals as a probe of structure formation models

Recent observations suggest that Carina and other nearby dwarf spheroidal galaxies (dSphs) are surrounded by unbound stars tidally stripped by the Milky Way. We have run high-resolution N-body simulations of dwarf galaxies orbiting within the Milky Way halo to determine if such observations can be explained with dark matter potentials like those implied by current structure formation models. We show that tidal forces acting on dwarfs with constant density cores or with cuspy profiles having a low concentration parameter c≲ 7) lead to flat outer stellar density profiles like that of Carina for a variety of orbital configurations. On the contrary, it is more difficult to remove stars from cuspy dark matter haloes with concentrations as high as predicted by cold dark matter (CDM) models at the mass scale of dwarf galaxies (c > 10), and the data can only be reproduced assuming nearly radial orbits. Our simulations show that Carina is losing mass at a fractional rate <0.1 Gyr−1 and its mass-to-light ratio could be inflated by at most a factor of 2 due to unbound stars projected along the line of sight. We follow the evolution of the tidal debris within a triaxial clumpy cold dark matter Milky Way halo which causes differential precession and small-scale heating of the stellar streams. This renders them useless as a dynamical tracer of the Galactic potential in CDM cosmogonies. Models with warm dark matter (WDM) or collisional fluid dark matter (FDM) produce dwarf haloes with lower central densities than CDM and would be consistent with the observed tidal tails even for orbits with eccentricities as low as indicated by current data on nearby dwarf spheroidals. Galactic haloes in models with FDM are expected to be smooth and spherical and would be favoured by the detection of cold coherent streams such as that associated with the Sagittarius dwarf spheroidal. In this respect the distribution of stellar streams on the plane of the sky can provide a novel test of the nature of the dark matter. On the other end, the streams could be similarly dispersed in lumpy CDM haloes and in smoother, triaxial WDM haloes.

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.

Formation and Structure of Halos in a Warm Dark Matter Cosmology

Using high-resolution cosmological N-body simulations, we study how the density profiles of dark matter halos are affected by the filtering of the density power spectrum below a given scale length and by the introduction of a thermal velocity dispersion. In the warm dark matter (WDM) scenario, both the free-streaming scale, Rf, and the velocity dispersion, v, are determined by the mass, mW, of the WDM particle. We found that v is too small to affect the density profiles of WDM halos. Down to the resolution attained in our simulations (~0.01 virial radii), there is not any significant difference in the density profiles and concentrations of halos obtained in simulations with and without the inclusion of v. Resolved soft cores appear only when we artificially increase the thermal velocity dispersion to a value that is much higher than v. We show that the size of soft cores in a monolithic collapse is related to the tangential velocity dispersion. The density profiles of the studied halos with masses down to ~0.01 the filtering mass Mf can be described by the Navarro-Frenk-White shape; soft cores are not formed. Nevertheless, the concentrations of these halos are lower than those of the CDM counterparts and are approximately independent of mass. The cosmogony of halos with masses Mf is not hierarchical: they form through monolithic collapse and by fragmentation of larger structures. The formation epoch of these halos is slightly later than that of halos with masses ≈Mf. The lower concentrations of WDM halos with respect to their CDM counterparts can be accounted for by their late formation epoch. Overall, our results point to a series of advantages of a WDM model over the CDM one. In addition to solving the substructure problem, a WDM model with Rf ~ 0.16 Mpc (mW ≈ 0.75 keV; flat cosmology with ΩΛ = h = 0.7) also predicts concentrations, a Tully-Fisher relation, and formation epochs for small halos, which seems to be in better agreement with observations than CDM predictions.