Warm Dark Matter Models Research Articles

Differentiating warm dark matter models through 21-cm line intensity mapping: A convolutional neural network approach

Differentiating warm dark matter models through 21-cm line intensity mapping: A convolutional neural network approach

New Galaxy UV Luminosity Constraints on Warm Dark Matter from JWST

We exploit the recent James Webb Space Telescope (JWST) determination of galaxy UV luminosity functions over the redshift range z = 9–14.5 to derive constraints on warm dark matter (WDM) models. The delayed structure formations in WDM universes make high-redshift observations a powerful probe to set limits on the particle mass m x of WDM candidates. By integrating these observations with blank-field surveys conducted by the Hubble Space Telescope (HST) at redshifts z = 4–8, we impose constraints on both astrophysical parameters (β, γ, ϵ N, and M c for a double-power-law star formation efficiency, and σMUV for a Gaussian magnitude–halo mass relation) and the WDM parameter (dark matter particle mass m x) simultaneously. We find a new limit of m x ≥ 3.2 keV for the mass of thermal relic WDM particles at a 95% confidence level. This bound is tighter than the most stringent result derived using HST data before. Future JWST observations could further reduce the observation uncertainties and improve this constraint.

Velocity Acoustic Oscillations on Cosmic Dawn 21 cm Power Spectrum as a Probe of Small-scale Density Fluctuations

We investigate the feasibility of using the velocity acoustic oscillations (VAO) features on the Cosmic Dawn 21 cm power spectrum to probe small-scale density fluctuations. In the standard cold dark matter (CDM) model, Population III stars form in minihalos and affect the 21 cm signal through Lyα and X-ray radiation. Such a process is modulated by the relative motion between dark matter and baryons, generating the VAO wiggles on the 21 cm power spectrum. In the fuzzy or warm dark matter models for which the number of minihalos is reduced, the VAO wiggles are weaker or even fully invisible. We investigate the wiggle features in the CDM with different astrophysical models and in different dark matter models. We find that (1) in the CDM model the relative streaming velocities can generate the VAO wiggles for broad ranges of parameters f *, ζ X , and f esc,LW ζ LW, though for different parameters the wiggles would appear at different redshifts and have different amplitudes. (2) For the axion model with m a ≲ 10−19 eV, the VAO wiggles are negligible. In the mixed model, the VAO signal is sensitive to the axion fraction. For example, the wiggles almost disappear when f a ≳ 10% for m a = 10−21 eV. Therefore, the VAO signal can be an effective indicator for small-scale density fluctuations and a useful probe of the nature of dark matter. The Square Kilometre Array-low with ∼2000 hr observation time has the ability to detect the VAO signal and constrain dark matter models.

Unveiling dark matter free streaming at the smallest scales with the high redshift Lyman-alpha forest

This study introduces novel constraints on the free streaming of thermal relic warm dark matter (WDM) from Lyman-α forest flux power spectra. Our analysis utilizes a high resolution, high redshift sample of quasar spectra observed using the HIRES and UVES spectrographs (z=4.2–5.0). We employ a Bayesian inference framework and a simulation-based likelihood that encompasses various parameters including the free streaming of dark matter, cosmological parameters, the thermal history of the intergalactic medium, and inhomogeneous reionization to establish lower limits on the mass of a thermal relic WDM particle of 5.7 keV (at 95% CL). This result surpasses previous limits from the Lyman-α forest through reduction of the measured uncertainties due to a larger statistical sample and by measuring clustering to smaller scales (kmax=0.2 km−1 s). The approximately two-fold improvement due to the expanded statistical sample suggests that the effectiveness of Lyman-α forest constraints on WDM models at high redshifts are limited by the availability of high quality quasar spectra. Restricting the analysis to comparable scales and thermal history priors as in prior studies (kmax<0.1 km−1 s) lowers the bound on the WDM mass to 4.1 keV. As the precision of the measurements increases, it becomes crucial to examine the instrumental and modeling systematics. On the modeling front, we argue that the impact of the thermal history uncertainty on the WDM particle mass constraint has diminished due to improved independent observations. At the smallest scales, the primary source of modeling systematic arises from the structure in the peculiar velocity of the intergalactic medium and inhomogeneous reionization. Published by the American Physical Society 2024

Optimal 1D Ly α forest power spectrum estimation – III. DESI early data

ABSTRACT The 1D power spectrum P1D of the Ly α forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of P1D with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of 54 600 quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with 2D image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly α emission line. In a companion paper, we present a similar analysis based on the Fast Fourier Transform estimate of the power spectrum. We conclude with a comparison of these two approaches and discuss the key sources of systematic error that we need to address with the upcoming DESI Year 1 analysis.

Implications of the Stellar Mass Density of High-z Massive Galaxies from JWST on Warm Dark Matter

A significant excess of the stellar mass density at high redshift has been discovered from the early data release of James Webb Space Telescope (JWST), and it may require a high star formation efficiency. However, this will lead to large number density of ionizing photons in the epoch of reionization (EoR), so that the reionization history will be changed, which can arise tension with the current EoR observations. Warm dark matter (WDM), via the free streaming effect, can suppress the formation of small-scale structure as well as low-mass galaxies. This provides an effective way to decrease the ionizing photons when considering a large star formation efficiency in high-z massive galaxies without altering the cosmic reionization history. On the other hand, the constraints on the properties of WDM can be derived from the JWST observations. In this work, we study WDM as a possible solution to reconcile the JWST stellar mass density of high-z massive galaxies and reionization history. We find that, the JWST high-z comoving cumulative stellar mass density alone has no significant preference for either CDM or WDM model. But using the observational data of other stellar mass density measurements and reionization history, we obtain that the WDM particle mass with keV and star formation efficiency parameter in 2σ confidence level can match both the JWST high-z comoving cumulative stellar mass density and the reionization history.

WIMP decay as a possible Warm Dark Matter model

Abstract The Weakly Interacting Massive Particles (WIMPs) have long been the favoured Cold Dark Matter (CDM) candidate in the standard ΛCDM model. However, owing to great improvement in the experimental sensitivity in the past decade, some parameter space of the Supersymmetric (SUSY)-based WIMP model is ruled out. In addition, a massive stable WIMP as the CDM particle is also at variance with other astrophysical observables at small scales. We consider a model that addresses both these issues. In the model, the WIMP decays into a massive particle and radiation. We study the background evolution and the first order perturbation theory (coupled Einstein-Boltzmann equations) for this model and show that the dynamics can be captured by a single parameter r = mL/q, which is the ratio of the lighter mass and the comoving momentum of the decay particle. We incorporate the relevant equations in the existing Boltzmann code CLASS to compute the matter power spectra and Cosmic Microwave Background (CMB) angular power spectra. The decaying WIMP model is akin to a non-thermal Warm Dark Matter (WDM) model and suppresses matter power at small scales, which could alleviate several issues that plague the CDM model at small scales. We compare the predictions of the model with CMB and galaxy clustering data. As the model deviates from the ΛCDM model at small scales, the evolution of the collapse fraction of matter in the universe is compared with the high-redshift Sloan Digital Sky Survey (SDSS) HI data. Both these data sets yield r ≳ 106, which can be translated into the bounds on other parameters. In particular, we obtain the following lower bounds on the thermally-averaged self-annihilation cross-section of WIMPs, ⟨σv⟩, and the lighter mass: ⟨σv⟩ ≳ 4.9 × 10-34 cm3 sec-1 and mL ≳ 2.4 keV. The lower limit on mL is comparable to constraints on the mass of thermally-produced WDM particle. The limit on the self-annihilation cross-section greatly expands the available parameter space as compared to the stable WIMP scenario.

JWST’s PEARLS: Mothra, a new kaiju star at z = 2.091 extremely magnified by MACS0416, and implications for dark matter models

We report the discovery of Mothra, an extremely magnified monster star, likely a binary system of two supergiant stars, in one of the strongly lensed galaxies behind the galaxy cluster MACS J0416.1−2403. Mothra is in a galaxy with spectroscopic redshift z = 2.091 in a portion of the galaxy that is parsecs away from the cluster caustic. The binary star is observed only on the side of the critical curve with negative parity but has been detectable for at least eight years, implying the presence of a small lensing perturber. Microlenses alone cannot explain the earlier observations of this object made with the Hubble Space Telescope. A larger perturber with a mass of at least 104 M⊙ offers a more satisfactory explanation. Based on the lack of perturbation on other nearby sources in the same arc, the maximum mass of the perturber is 2.5 × 106 M⊙, making this the smallest substructure constrained by lensing at z > 0.3. The existence of this millilens is fully consistent with expectations from standard cold dark matter cosmology. On the other hand, the existence of such a small substructure in a cluster environment has implications for other dark matter models. In particular, warm dark matter models with particle masses below 8.7 keV are excluded by our observations. Similarly, axion dark matter models are consistent with the observations only if the axion mass is in the range 0.5 × 10−22 eV < ma < 5 × 10−22 eV.

Inferring warm dark matter masses with deep learning

ABSTRACT We present a new suite of over 1500 cosmological N-body simulations with varied warm dark matter (WDM) models ranging from 2.5 to 30 keV. We use these simulations to train Convolutional Neural Networks (CNNs) to infer WDM particle masses from images of DM field data. Our fiducial setup can make accurate predictions of the WDM particle mass up to 7.5 keV with an uncertainty of ±0.5 keV at a 95 per cent confidence level from (25 h−1Mpc)2 maps. We vary the image resolution, simulation resolution, redshift, and cosmology of our fiducial setup to better understand how our model is making predictions. Using these variations, we find that our models are most dependent on simulation resolution, minimally dependent on image resolution, not systematically dependent on redshift, and robust to varied cosmologies. We also find that an important feature to distinguish between WDM models is present with a linear size between 100 and 200 h−1 kpc. We compare our fiducial model to one trained on the power spectrum alone and find that our field-level model can make two times more precise predictions and can make accurate predictions to two times as massive WDM particle masses when used on the same data. Overall, we find that the field-level data can be used to accurately differentiate between WDM models and contain more information than is captured by the power spectrum. This technique can be extended to more complex DM models and opens up new opportunities to explore alternative DM models in a cosmological environment.

Extending the unified subhalo model to warm dark matter

ABSTRACT Using a set of high-resolution N-body simulations, we extend the unified distribution model of cold dark matter (CDM) subhaloes to the warm dark matter (WDM) case. The same model framework combining the unevolved mass function, unevolved radial distribution, and tidal stripping can predict the mass function and spatial distribution of subhaloes in both CDM and WDM simulations. The dependence of the model on the DM particle property is universally parametrized through the half-mode mass of the initial power spectrum. Compared with the CDM model, the WDM model differs most notably in two aspects. (1) In contrast to the power-law form in CDM, the unevolved subhalo mass function for WDM is scale-dependent at the low mass end due to the cut-off in the initial power spectrum. (2) WDM subhaloes are more vulnerable to tidal stripping and disruption due to their lower concentrations at accretion time. Their survival rate is also found to depend on the infall mass. Accounting for these differences, the model predicts a final WDM subhalo mass function that is also proportional to the unevolved subhalo mass function. The radial distribution of WDM subhaloes is predicted to be mass-dependent. For low mass subhaloes, the radial distribution is flatter in the inner halo and steeper in the outer halo compared to the CDM counterpart, due to the scale-dependent unevolved mass function and the enhanced tidal stripping. The code for sampling subhaloes according to our generalized model is available at https://github.com/fhtouma/subgen2.

Pushing the limits of detectability: mixed dark matter from strong gravitational lenses.

One of the frontiers for advancing what is known about dark matter lies in using strong gravitational lenses to characterize the population of the smallest dark matter haloes. There is a large volume of information in strong gravitational lens images - the question we seek to answer is to what extent we can refine this information. To this end, we forecast the detectability of a mixed warm and cold dark matter scenario using the anomalous flux ratio method from strong gravitational lensed images. The halo mass function of the mixed dark matter scenario is suppressed relative to cold dark matter but still predicts numerous low-mass dark matter haloes relative to warm dark matter. Since the strong lensing signal receives a contribution from a range of dark matter halo masses and since the signal is sensitive to the specific configuration of dark matter haloes, not just the halo mass function, degeneracies between different forms of suppression in the halo mass function, relative to cold dark matter, can arise. We find that, with a set of lenses with different configurations of the main deflector and hence different sensitivities to different mass ranges of the halo mass function, the different forms of suppression of the halo mass function between the warm dark matter model and the mixed dark matter model can be distinguished with 40 lenses with Bayesian odds of 30:1.

Implications of multi-axion dark matter on structure formation

Axions are candidates for dark matter in the universe.We develop an accurate Boltzmann code to calculate the linear growth of the plasma. As an interesting example, we investigate a mixed dark matter model consisting of cold dark matter (CDM) and two-axion dark matter. We analyze the growth of the structure numerically and analytically. We find that an effective single axion with an effective mass and an effective abundance is useful to characterize the two-axion cosmology. Moreover, we generalize the effective single axion description to multi-axion dark matter cosmology. We also compare the results with those of warm dark matter (WDM) model. Moreover, we calculate halo mass functions for the mixed model by using the Press-Schechter model and linear perturbations and then determine the mass function as a function of masses and axion abundance.

JWST high-redshift galaxy constraints on warm and cold dark matter models

Context. Warm dark matter is a possible alternative to cold dark matter to explain cosmological structure formation. Aims. We study the implications of the latest JWST data on the nature of dark matter. Methods. We compare properties of high-redshift galaxies observed by JWST with hydrodynamical simulations, in the standard cold dark matter model and in warm dark matter models with a suppressed linear matter power spectrum Results. We find that current data are neither in tension with cold dark matter nor with warm dark matter models with mWDM > 2 keV, since they probe bright and rare objects whose physical properties are similar in the different scenarios. Conclusions. We also show how two observables, the galaxy luminosity functions and the galaxy correlation function at small scales of faint objects, can be promising tools for discriminating between the different dark-matter scenarios. Further hints may come from early stellar-mass statistics and galaxy CO emission.

Massive prompt cusps: a new signature of warm dark matter

ABSTRACT Every dark matter halo and subhalo is expected to have a prompt ρ ∝ r−1.5 central density cusp, which is a relic of its condensation out of the smooth mass distribution of the early Universe. The sizes of these prompt cusps are linked to the scales of the peaks in the initial density field from which they formed. In warm dark matter (WDM) models, the smoothing scale set by free streaming of the dark matter can result in prompt cusps with masses of order 107 M⊙. We show that WDM models with particle masses ranging from 2 to 6 keV predict prompt cusps that could detectably alter the observed kinematics of Local Group dwarf galaxies. Thus, prompt cusps present a viable new probe of WDM. A prompt cusp’s properties are highly sensitive to when it formed, so prospects can be improved with a better understanding of when the haloes of the Local Group dwarfs originally formed. Tidal stripping can also affect prompt cusps, so constraints on satellite galaxy orbits can further tighten WDM inferences.

Sensitivity of strong lensing observations to dark matter substructure: a case study with Euclid

ABSTRACTWe introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a vmax–rmax relation. We then estimate the dark matter subhalo sensitivity in 16 000 simulated strong lens observations with depth and resolution resembling Euclid VIS images. We find that with a 3σ detection threshold, 2.35 per cent of pixels inside twice the Einstein radius are sensitive to subhaloes with a mass Mmax ≤ 1010 M⊙, 0.03 per cent are sensitive to Mmax ≤ 109 M⊙, and the limit of sensitivity is found to be Mmax = 108.8 ± 0.2 M⊙. Using our sensitivity maps and assuming CDM, we estimate that Euclid-like lenses will yield $1.43^{+0.14}_{-0.11}[f_\mathrm{sub}^{-1}]$ detectable subhaloes per lens in the entire sample, but this increases to $35.6^{+0.9}_{-0.9}[f_\mathrm{sub}^{-1}]$ per lens in the most sensitive lenses. Estimates are given in units of the inverse of the substructure mass fraction $f_\mathrm{sub}^{-1}$. Assuming fsub = 0.01, one in every 70 lenses in general should yield a detection, or one in every ∼ three lenses in the most sensitive sample. From 170 000 new strong lenses detected by Euclid, we expect ∼2500 new subhalo detections. We find that the expected number of detectable subhaloes in warm dark matter models only changes relative to cold dark matter for models which have already been ruled out, i.e. those with half-mode masses Mhm > 108 M⊙.

Probing dark matter interactions with 21cm observations

Similarly to warm dark matter which features a cut-off in the matter power spectrum due to free-streaming, many interacting dark matter models predict a suppression of the matter power spectrum on small length scales through collisional damping. Forecasts for 21cm line intensity mapping have shown that an instrument like the SKA will be able to probe a suppression of power in warm dark matter scenarios in a statistically significant way. Here we investigate the implications of these findings on interacting dark matter scenarios, particularly dark matter-neutrino interactions, which we use as an example. Using a suite of cosmological N-body simulations, we demonstrate that interacting scenarios show a suppression of the non-linear power spectrum similar to warm dark matter models. This implies that 21cm line intensity mapping will be able to set the strongest limits yet on dark matter-neutrino scattering, improving the constraints by two orders of magnitude over current Lyman-α bounds, and by four orders of magnitude over cosmic microwave background and baryon acoustic oscillations limits. However, to distinguish between warm dark matter and interacting scenarios, our simulations show that percent-level precision measurements of the matter power spectrum at redshifts z ≳ 15 are necessary, as the key features of interacting scenarios are washed out by non-linear evolution at later times.

Small-scale clumping of dark matter and the mean free path of ionizing photons at z = 6

Recently, the mean free path of ionizing photons in the z = 6 intergalactic medium (IGM) was measured to be very short, presenting a challenge to existing reionization models. At face value, the measurement can be interpreted as evidence that the IGM clumps on scales M ≲ 108 M ⊙, a key but largely untested prediction of the cold dark matter (CDM) paradigm. Motivated by this possibility, we study the role that the underlying dark matter cosmology plays in setting the z > 5 mean free path. We use two classes of models to contrast against the standard CDM prediction: (1) thermal relic warm dark matter (WDM), representing models with suppressed small-scale power; (2) an ultralight axion exhibiting a white noise-like power enhancement. Differences in the mean free path between the WDM and CDM models are subdued by pressure smoothing and the possible contribution of neutral islands to the IGM opacity. For example, comparing late reionization scenarios with a fixed volume-weighted mean neutral fraction of 20% at z = 6, the mean free path is 19 (45)% longer in a WDM model with mx = 3 (1) keV. The enhanced power in the axion-like model produces better agreement with the short mean free path measured at z = 6. However, drawing robust conclusions about cosmology is hampered by large uncertainties in the reionization process, extragalactic ionizing background, and thermal history of the Universe. This work highlights some key open questions about the IGM opacity during reionization.

Warm dark matter constraints using Milky Way satellite observations and subhalo evolution modeling

Warm dark matter (WDM) can potentially explain small-scale observations that currently challenge the cold dark matter (CDM) model, as warm particles suppress structure formation due to free-streaming effects. Observing small-scale matter distribution provides a valuable way to distinguish between CDM and WDM. In this work, we use observations from the Dark Energy Survey and PanSTARRS1, which observe 270 Milky-Way satellites after completeness corrections. We test WDM models by comparing the number of satellites in the Milky Way with predictions derived from the Semi-Analytical SubHalo Inference ModelIng (SASHIMI) code, which we develop based on the extended Press-Schechter formalism and subhalos' tidal evolution prescription. We robustly rule out WDM with masses lighter than 4.4 keV at 95% confidence level for the Milky-Way halo mass of $10^{12} M_\odot$. The limits are a weak function of the (yet uncertain) Milky-Way halo mass, and vary as $m_{\rm WDM}>3.6$-$5.1$ keV for $(0.6$-$2.0) \times 10^{12} M_\odot$. For the sterile neutrinos that form a subclass of WDM, we obtain the constraints of $m_{\nu_s}>11.6$ keV for the Milky-Way halo mass of $10^{12} M_{\odot}$. These results based on SASHIMI do not rely on any assumptions of galaxy formation physics or are not limited by numerical resolution. The models, therefore, offer a robust and fast way to constrain the WDM models. By applying a satellite forming condition, however, we can rule out the WDM mass lighter than 9.0 keV for the Milky-Way halo mass of $10^{12} M_\odot$.

Lensed radio arcs at milli-arcsecond resolution: Methods, science results, and current status

AbstractStrong gravitational lensing by galaxies provides us with a powerful laboratory for testing dark matter models. Various particle models for dark matter give rise to different small-scale distributions of mass in the lens galaxy, which can be differentiated with sensitive observations. Th. The sensitivity of a gravitational lens observation to the presence (or absence) of low-mass dark structures in the lens galaxy is determined mainly by the angular resolution of the instrument and the spatial structure of the source. Here, I discuss results from the analysis of a global VLBI observation of a gravitationally lensed radio jet. With an angular resolution better than 5 milli-arcseconds and a highly extended, spatially resolved source, we are able to place competitive constraints on the particle mass in fuzzy dark matter models using this single observation. I also discuss preliminary results from our analysis of warm dark matter models using this lens system.

Molecular Chemistry for Dark Matter. II. Recombination, Molecule Formation, and Halo Mass Function in Atomic Dark Matter

Dissipative dark matter predicts rich observable phenomena that can be tested with future large-scale structure surveys. As a specific example, we study atomic dark matter, consisting of a heavy particle and a light particle charged under a dark electromagnetism. In particular, we calculate the cosmological evolution of atomic dark matter focusing on dark recombination and dark molecule formation. We have obtained the relevant interaction rate coefficients by rescaling the rates for normal hydrogen, and evolved the abundances for ionized, atomic, and molecular states using a modified version of Recfast++ (which we have released publicly at a a https://github.com/jamesgurian/RecfastJulia ). We also provide an analytical approximation for the final abundances. We then calculate the effects of atomic dark matter on the linear power spectrum, which enter through a dark photon diffusion and dark acoustic oscillations. At formation time, the atomic dark matter model suppresses halo abundances on scales smaller than the diffusion scale, just as warm dark matter models suppress the abundance below the free-streaming scale. The subsequent evolution with radiative cooling, however, will alter the halo mass function further.