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
Motivations for a large self-interacting dark matter cross-section from Milky Way satellites
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
-
-
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.
-
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
-
more » « less
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.
-
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
