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- Award ID(s):
- 1915005
- NSF-PAR ID:
- 10381881
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 518
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2418-2435
- Size(s):
- p. 2418-2435
- Sponsoring Org:
- National Science Foundation
More Like this
-
ABSTRACT Within lambda cold dark matter ($\Lambda$CDM), dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhaloes, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$$\gt 10^5$ M$_{\odot }$; Carina and Fornax) as well as several ultrafaint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$$\sim$6$\times 10^ {3}$ M$_{\odot }$) substructure population of isolated LMC-mass hosts ($M_{\text{200m}}$ = 1–3$\times 10^ {11}$ M$_{\odot }$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhaloes in simulations with and without baryons, we find that, within 0.2 $r_{\text{200m}}$, LMC-mass hosts deplete $\sim$30 per cent of their substructure, significantly lower than the $\sim$70 per cent of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot }$), $\sim 5\!-\!10$ subhaloes form galaxies with $M_\ast$$\ge 10^{4}$ M$_{\odot }$ , in agreement with the seven observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultrafaint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association.more » « less
-
ABSTRACT Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section σ/m. We compare the shape of the main dark matter (DM) halo at redshift z = 0 predicted by SIDM simulations (at σ/m = 0.1, 1, and 10 cm2 g−1) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on σ/m. This transition could potentially be used to set limits on the SIDM cross-section in the MW.
-
null (Ed.)ABSTRACT Recent progress in constraining the massive accretions (>1:10) experienced by the Milky Way (MW) and the Andromeda galaxy (M31) offers an opportunity to understand the dwarf galaxy population of the Local Group. Using zoom-in dark matter-only simulations of MW-mass haloes and concentrating on subhaloes that are thought to be capable of hosting dwarf galaxies, we demonstrate that the infall of a massive progenitor is accompanied with the accretion and destruction of a large number of subhaloes. Massive accreted progenitors do not increase the total number of infalling subhaloes on to a MW-mass host, but instead focus surrounding subhaloes on to the host causing a clustering in the infall time of subhaloes. This leads to a temporary elevation in the number of subhaloes as well as changes in their cumulative radial profile within the virial radius of the host. Surviving subhaloes associated with a massive progenitor have a large diversity in their orbits. We find that the star formation quenching times of Local Group dwarf spheroidal galaxies ($10^{5} \mathrm{\, M_{\odot }} \lesssim \mathrm{\mathit{ M}}_{*} \lesssim 10^{7} \mathrm{\, M_{\odot }}$) are clustered around the times of the most massive accretions suffered by the MW and M31. Our results imply that (a) the quenching time of dwarf spheroidals is a good proxy of their infall time and b) the absence of recently quenched satellites around M31 suggests that M33 is not on its first infall and was accreted much earlier.more » « less
-
The onset of gravothermal core collapse in velocity-dependent self-interacting dark matter subhaloesnull (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
-
Abstract We present the results of dark matter (DM) searches in a sample of 31 dwarf irregular (dIrr) galaxies within the field of view of the HAWC Observatory. dIrr galaxies are DM-dominated objects in which astrophysical gamma-ray emission is estimated to be negligible with respect to the secondary gamma-ray flux expected by annihilation or decay of weakly interacting massive particles (WIMPs). While we do not see any statistically significant DM signal in dIrr galaxies, we present the exclusion limits (95% C.L.) for annihilation cross section and decay lifetime for WIMP candidates with masses between 1 and 100 TeV. Exclusion limits from dIrr galaxies are relevant and complementary to benchmark dwarf Spheroidal (dSph) galaxies. In fact, dIrr galaxies are targets kinematically different from benchmark dSph, preserving the footprints of different evolution histories. We compare the limits from dIrr galaxies to those from ultrafaint and classical dSph galaxies previously observed with HAWC. We find that the constraints are comparable to the limits from classical dSph galaxies and ∼2 orders of magnitude weaker than the ultrafaint dSph limits.more » « less