skip to main content


Title: 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.

 
more » « less
NSF-PAR ID:
10473076
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
527
Issue:
1
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 739-755
Size(s):
["p. 739-755"]
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Ly α forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of λhm = 0.089 Mpc h−1 at the 95 per cent confidence level, which we show to be stable for a broad range of prior choices. Assuming a Planck cosmology and that WDM particles are thermal relics, this corresponds to an upper limit on the half-mode mass of Mhm < 3 × 107 M⊙ h−1, and a lower limit on the particle mass of mth > 6.048 keV, both at the 95 per cent confidence level. We find that models with λhm > 0.223 Mpc h−1 (corresponding to mth > 2.552 keV and Mhm < 4.8 × 108 M⊙ h−1) are ruled out with respect to the maximum likelihood model by a factor ≤1/20. For lepton asymmetries L6 > 10, we rule out the 7.1 keV sterile neutrino dark matter model, which presents a possible explanation to the unidentified 3.55 keV line in the Milky Way and clusters of galaxies. The inferred 95 percentiles suggest that we further rule out the ETHOS-4 model of self-interacting DM. Our results highlight the importance of extending the current constraints to lower half-mode scales. We address important sources of systematic errors and provide prospects for how the constraints of these probes can be improved upon in the future. 
    more » « less
  2. ABSTRACT

    Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, bv – defined as the ratio of galaxy and DM velocity dispersions – as a function of redshift, host halo mass, and galaxy stellar mass threshold ($M_{\rm \star , sat}$), for massive haloes ($M_{\rm 200c}\gt 10^{13.5} \, {\rm M}_\odot$) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of bv increasing with M200c and decreasing with redshift and $M_{\rm \star , sat}$. The ensemble-estimated theoretical uncertainty in bv is 2–3 per cent, but becomes percent-level when considering only the three highest resolution simulations. We update the mass–richness normalization for an SDSS redMaPPer cluster sample, and find our improved bv estimates reduce the normalization uncertainty from 22 to 8 per cent, demonstrating that dynamical mass estimation is competitive with weak lensing mass estimation. We discuss necessary steps for further improving this precision. Our estimates for $b_v(M_{\rm 200c}, M_{\rm \star , sat}, z)$ are made publicly available.

     
    more » « less
  3. ABSTRACT The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter haloes on sub-galactic scales. Using the image positions and flux ratios from eight quadruply imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhaloes and haloes along the line of sight, and account for warm dark matter free-streaming effects on the mass function and mass–concentration relation. By calibrating the scaling of the SHMF with host halo mass and redshift using a suite of simulated haloes, we infer a global normalization for the SHMF. We account for finite-size background sources, and marginalize over the mass profile of the main deflector. Parametrizing dark matter free-streaming through the half-mode mass mhm, we constrain the thermal relic particle mass mDM corresponding to mhm. At $95 \, {\rm per\, cent}$ CI: mhm < 107.8 M⊙ ($m_{\rm {DM}} \gt 5.2 \ \rm {keV}$). We disfavour $m_{\rm {DM}} = 4.0 \,\rm {keV}$ and $m_{\rm {DM}} = 3.0 \,\rm {keV}$ with likelihood ratios of 7:1 and 30:1, respectively, relative to the peak of the posterior distribution. Assuming cold dark matter, we constrain the projected mass in substructure between 106 and 109 M⊙ near lensed images. At $68 \, {\rm per\, cent}$ CI, we infer $2.0{-}6.1 \times 10^{7}\, {{\rm M}_{\odot }}\,\rm {kpc^{-2}}$, corresponding to mean projected mass fraction $\bar{f}_{\rm {sub}} = 0.035_{-0.017}^{+0.021}$. At $95 \, {\rm per\, cent}$ CI, we obtain a lower bound on the projected mass of $0.6 \times 10^{7} \,{{\rm M}_{\odot }}\,\rm {kpc^{-2}}$, corresponding to $\bar{f}_{\rm {sub}} \gt 0.005$. These results agree with the predictions of cold dark matter. 
    more » « less
  4. ABSTRACT Reticulum II (Ret II) is a satellite galaxy of the Milky Way (MW) and presents a prime target to investigate the nature of dark matter (DM) because of its high mass-to-light ratio. We evaluate a dedicated INTEGRAL observation campaign data set to obtain γ-ray fluxes from Ret II and compare those with expectations from DM. Ret II is not detected in the γ-ray band 25–8000 keV, and we derive a flux limit of ${\lesssim}10^{-8}\, \mathrm{erg\, cm^{-2}\, s^{-1}}$. The previously reported 511 keV line is not seen, and we find a flux limit of ${\lesssim}1.7 \times 10^{-4}\, \mathrm{ph\, cm^{-2}\, s^{-1}}$. We construct spectral models for primordial black hole (PBH) evaporation and annihilation/decay of particle DM, and subsequent annihilation of e+s produced in these processes. We exclude that the totality of DM in Ret II is made of a monochromatic distribution of PBHs of masses ${\lesssim}8 \times 10^{15}\, \mathrm{g}$. Our limits on the velocity-averaged DM annihilation cross section into e+e− are $\langle \sigma v \rangle \lesssim 5 \times 10^{-28} \left(m_{\rm DM} / \mathrm{MeV} \right)^{2.5}\, \mathrm{cm^3\, s^{-1}}$. We conclude that analysing isolated targets in the MeV γ-ray band can set strong bounds on DM properties without multi-year data sets of the entire MW, and encourage follow-up observations of Ret II and other dwarf galaxies. 
    more » « less
  5. Cosmological simulations of galaxy formation are limited by finite computational resources. We draw from the ongoing rapid advances in artificial intelligence (AI; specifically deep learning) to address this problem. Neural networks have been developed to learn from high-resolution (HR) image data and then make accurate superresolution (SR) versions of different low-resolution (LR) images. We apply such techniques to LR cosmological N-body simulations, generating SR versions. Specifically, we are able to enhance the simulation resolution by generating 512 times more particles and predicting their displacements from the initial positions. Therefore, our results can be viewed as simulation realizations themselves, rather than projections, e.g., to their density fields. Furthermore, the generation process is stochastic, enabling us to sample the small-scale modes conditioning on the large-scale environment. Our model learns from only 16 pairs of small-volume LR-HR simulations and is then able to generate SR simulations that successfully reproduce the HR matter power spectrum to percent level up to16h1Mpcand the HR halo mass function to within10%down to1011M. We successfully deploy the model in a box 1,000 times larger than the training simulation box, showing that high-resolution mock surveys can be generated rapidly. We conclude that AI assistance has the potential to revolutionize modeling of small-scale galaxy-formation physics in large cosmological volumes.

     
    more » « less