skip to main content

Title: Inferring warm dark matter masses with deep learning

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
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Medium: X Size: p. 739-755
["p. 739-755"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)–mass cosmological zoom-in simulations and an empirical galaxy–halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy–halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses ofM50=108Mat the 1σlevel; however, a tail toward lowM50prevents a 2σmeasurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties onM50at the 1σlevel, but the 2σtail toward lowM50persists. We project that observations of one (two) complete satellite populations can constrain warm DM models withmWDM≈ 10 keV (20 keV). Subhalo mass function (SHMF) suppression can be constrained to ≈70%, 60%, and 50% that in cold dark matter (CDM) at peak virial masses of 108, 109, and 1010M, respectively; SHMF enhancement constraints are weaker (≈20, 4, and 2 times that in CDM, respectively) due to galaxy–halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics.

    more » « less
  2. 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

    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
  4. 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
  5. 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