skip to main content

Explore Research Products in the NSF-PAR

Title: The Milky Way’s halo and subhaloes in self-interacting dark matter
ABSTRACT

We perform high-resolution simulations of an MW-like galaxy in a self-interacting cold dark matter model with elastic cross-section over mass of $1~\rm cm^2\, g^{-1}$ (SIDM) and compare to a model without self-interactions (CDM). We run our simulations with and without a time-dependent embedded potential to capture effects of the baryonic disc and bulge contributions. The CDM and SIDM simulations with the embedded baryonic potential exhibit remarkably similar host halo profiles, subhalo abundances, and radial distributions within the virial radius. The SIDM host halo is denser in the centre than the CDM host and has no discernible core, in sharp contrast to the case without the baryonic potential (core size ${\sim}7 \, \rm kpc$). The most massive subhaloes (with $V_{\mathrm{peak}}\gt 20 \, \rm km\, s^{-1}$) in our SIDM simulations, expected to host the classical satellite galaxies, have density profiles that are less dense than their CDM analogues at radii less than 500 pc but the deviation diminishes for less massive subhaloes. With the baryonic potential included in the CDM and SIDM simulations, the most massive subhaloes do not display the too-big-to-fail problem. However, the least dense among the massive subhaloes in both these simulations tend to have the smallest pericenter values, a trend that is not apparent among the bright MW satellite galaxies.

  more » « less
NSF-PAR ID:
10121811
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
490
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2117-2123
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT It has been proposed that gravothermal collapse due to dark matter self-interactions (i.e. self-interacting dark matter, SIDM) can explain the observed diversity of the Milky Way (MW) satellites’ central dynamical masses. We investigate the process behind this hypothesis using an N-body simulation of a MW-analogue halo with velocity-dependent SIDM (vdSIDM) in which the low-velocity self-scattering cross-section, $\sigma _{\rm T}/m_{\rm x}$, reaches 100 cm2 g−1; we dub this model the vd100 model. We compare the results of this simulation to simulations of the same halo that employ different dark models, including cold dark matter (CDM) and other, less extreme SIDM models. The masses of the vd100 haloes are very similar to their CDM counterparts, but the values of their maximum circular velocities, Vmax, are significantly higher. We determine that these high Vmax subhaloes were objects in the mass range [5 × 106, 1 × 108] M⊙ at z = 1 that undergo gravothermal core collapse. These collapsed haloes have density profiles that are described by single power laws down to the resolution limit of the simulation, and the inner slope of this density profile is approximately −3. Resolving the ever decreasing collapsed region is challenging, and tailored simulations will be required to model the runaway instability accurately at scales <1 kpc. 
    more » « less
  2. ; ; ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the ‘Feedback In Realistic Environment’ (FIRE) project. The three simulated haloes have virial masses of $\sim 10^{12}\, \text{M}_\odot$ at z = 0, and we study velocity-independent self-interaction cross sections of 1 and 10 ${\rm cm^2 \, g^{-1}}$. We study star formation rates and the shape of dark matter density profiles of the parent haloes in both cold dark matter (CDM) and SIDM models. Galaxies formed in the SIDM haloes have higher star formation rates at z ≤ 1, resulting in more massive galaxies compared to the CDM simulations. While both CDM and SIDM simulations show diverse shape of the dark matter density profiles, the SIDM haloes can reach higher and more steep central densities within few kpcs compared to the CDM haloes. We identify a correlation between the build-up of the stars within the half-mass radii of the galaxies and the growth in the central dark matter densities. The thermalization process in the SIDM haloes is enhanced in the presence of a dense stellar component. Hence, SIDM haloes with highly concentrated baryonic profiles are predicted to have higher central dark matter densities than the CDM haloes. Overall, the SIDM haloes are more responsive to the presence of a massive baryonic distribution than their CDM counterparts. 
    more » « less
  3. ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    We explore the properties of Milky Way (MW) subhaloes in self-interacting dark matter models for moderate cross-sections of 1–5 cm2 g−1 using high-resolution zoom-in N-body simulations. We include the gravitational potential of a baryonic disc and bulge matched to the MW, which is critical for getting accurate predictions. The predicted number and distribution of subhaloes within the host halo are similar for 1 and 5 cm2 g−1 models, and they agree with observations of MW satellite galaxies only if subhaloes with peak circular velocity over all time >7 km s−1 are able to form galaxies. We do not find distinctive signatures in the pericentre distribution of the subhaloes that could help distinguish the models. Using an analytical model to extend the simulation results, we are able to show that subhaloes in models with cross-sections between 1 and 5 cm2 g−1 are not dense enough to match the densest ultrafaint and classical dwarf spheroidal galaxies in the MW. This motivates exploring velocity-dependent cross-sections with values larger than 5 cm2 g−1 at the velocities relevant for the satellites such that core collapse would occur in some of the ultrafaint and classical dwarf spheroidals.

     
    more » « less
  4. ; ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT We derive a new mass estimator that relies on internal proper motion measurements of dispersion-supported stellar systems, one that is distinct and complementary to existing estimators for line-of-sight velocities. Starting with the spherical Jeans equation, we show that there exists a radius where the mass enclosed depends only on the projected tangential velocity dispersion, assuming that the anisotropy profile slowly varies. This is well-approximated at the radius where the log-slope of the stellar tracer profile is −2: r−2. The associated mass is $M(r_{-2}) = 2 G^{-1} \langle \sigma _{\mathcal {T}}^{2}\rangle ^{*} r_{-2}$ and the circular velocity is $V^{2}({r_{-2}}) = 2\langle \sigma _{\mathcal {T}}^{2}\rangle ^{*}$. For a Plummer profile r−2 ≃ 4Re/5. Importantly, r−2 is smaller than the characteristic radius for line-of-sight velocities derived by Wolf et al. Together, the two estimators can constrain the mass profiles of dispersion-supported galaxies. We illustrate its applicability using published proper motion measurements of dwarf galaxies Draco and Sculptor, and find that they are consistent with inhabiting cuspy NFW subhaloes of the kind predicted in CDM but we cannot rule out a core. We test our combined mass estimators against previously published, non-spherical cosmological dwarf galaxy simulations done in both cold dark matter (CDM; naturally cuspy profile) and self-interacting dark matter (SIDM; cored profile). For CDM, the estimates for the dynamic rotation curves are found to be accurate to $10\rm { per\, cent}$ while SIDM are accurate to $15\rm { per\, cent}$. Unfortunately, this level of accuracy is not good enough to measure slopes at the level required to distinguish between cusps and cores of the type predicted in viable SIDM models without stronger priors. However, we find that this provides good enough accuracy to distinguish between the normalization differences predicted at small radii (r ≃ r−2 < rcore) for interesting SIDM models. As the number of galaxies with internal proper motions increases, mass estimators of this kind will enable valuable constraints on SIDM and CDM models. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_{\rm halo} \approx 10^{10}\, {\rm M}_{\odot }$ at z = 0, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r1/2 < 500 pc) have lower σ⋆/Vmax ratios, reinforcing the idea that smaller dwarfs may reside in haloes that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in ΛCDM. The V1/2–r1/2 relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Overall, the inclusion of baryons substantially reduces any distinct signatures of dark matter physics in the observable properties of dwarf galaxies. Spatially resolved rotation curves in the central regions (<400 pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r1/2 values can only originate from dark matter self-interactions. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email