skip to main content


Title: Warm dark matter chills out: constraints on the halo mass function and the free-streaming length of dark matter with eight quadruple-image strong gravitational lenses
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
Award ID(s):
1714953
NSF-PAR ID:
10175463
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
491
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
6077 to 6101
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $0.5{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation. 
    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

    We analyse strongly lensed images in eight galaxy clusters to measure their dark matter density profiles in the radial region between 10 kpc and 150 kpc, and use this to constrain the self-interaction cross-section of dark matter (DM) particles. We infer the mass profiles of the central DM haloes, bright central galaxies, key member galaxies, and DM subhaloes for the member galaxies for all eight clusters using the qlens code. The inferred DM halo surface densities are fit to a self-interacting dark matter model, which allows us to constrain the self-interaction cross-section over mass σ/m. When our full method is applied to mock data generated from two clusters in the Illustris-TNG simulation, we find results consistent with no dark matter self-interactions as expected. For the eight observed clusters with average relative velocities of $1458_{-81}^{+80}$ km s−1, we infer $\sigma /m = 0.082_{-0.021}^{+0.027} \rm cm^2\, g^{ -1}$ and $\sigma /m \lt 0.13~ \rm cm^2\, g^{ -1}$ at the 95 per cent confidence level.

     
    more » « less
  4. ABSTRACT

    We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple massive black hole (MBH) systems at z = 2–3. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find 4 per cent of the $M_{\rm BH}\gt 10^7\, M_\odot$ are in tuples with $\Delta r_{\rm max} \lt 200\, {\rm kpc}$. The tuple systems span a range of separations with the majority of the observable AGN systems at Δr ∼ 50–100 kpc. They include some of the most massive BHs (up to $10^{10} \, M_\odot$) but with at least one of the components of $M_{\rm BH} \sim 10^7 \, {\rm M}_{\odot }$. Tuples’ host galaxies are typically massive with $M_* \sim 10^{10-11} \, M_\odot$. We find that $\gt 10~{{\ \rm per\ cent}}$ massive haloes with Mhalo > 1013 M⊙ host MBH tuples. Following the subsequent interactions between MBHs in tuples, we found that in $\sim 5~{{\ \rm per\ cent}}$ of the triplets all three MBHs merge within a Gyr, and 15 per cent go through one merger. As a by-product of the complex multigalaxy interaction of these systems, we also find that up to $\sim 5~{{\ \rm per\ cent}}$ of tuples lead to runaway MBHs. In ASTRID, virtually all of the ultramassive black holes ($\gt 10^{10} \, M_\odot$) have undergone a triple quasar phase, while for BHs with $M_{\rm BH} \sim 10^9 \, M_\odot$, this fraction drops to 50 per cent.

     
    more » « less
  5. ABSTRACT

    We use GRUMPY, a simple regulator-type model for dwarf galaxy formation and evolution, to forward model the dwarf galaxy satellite population of the Milky Way (MW) using the Caterpillar zoom-in simulation suite. We show that luminosity and distance distributions of the model satellites are consistent with the distributions measured in the DES, PS1, and SDSS surveys, even without including a model for the orphan galaxies. We also show that our model for dwarf galaxy sizes can simultaneously reproduce the observed distribution of stellar half-mass radii, r1/2, of the MW satellites and the overall r1/2–M⋆ relation exhibited by observed dwarf galaxies. The model predicts that some of the observed faint stellar systems with r1/2 < 10 pc are ultra-faint dwarf galaxies. Scaling of the stellar mass M⋆ and peak halo mass Mpeak for the model satellites is not described by a power law, but has a clear flattening of M⋆–Mpeak scaling at $M_{\rm peak}\lt 10^8\, \, M_{\odot }$ imprinted by reionization. As a result, the fraction of low mass haloes ($M_{\rm peak}\lt 10^8 \, M_{\odot }$) hosting galaxies with MV < 0 is predicted to be 50 per cent at $M_{\rm peak}\sim 3.6 \times 10^7\, \, M_{\odot }$. We find that such high fraction at that halo mass helps to reproduce the number of dwarf galaxies discovered recently in the HSC-SSP survey. Using the model we forecast that there should be the total of $440^{+201}_{-147}$ (68 per cent confidence interval) MW satellites with MV < 0 and r1/2 > 10 pc within 300 kpc and make specific predictions for the HSC-SSP, DELVE-WIDE, and LSST surveys.

     
    more » « less