skip to main content

Title: Endothermic self-interacting dark matter in Milky Way-like dark matter haloes

Self-interacting dark matter (SIDM) offers the potential to mitigate some of the discrepancies between simulated cold dark matter (CDM) and observed galactic properties. We introduce a physically motivated SIDM model to understand the effects of self interactions on the properties of Milky Way and dwarf galaxy sized haloes. This model consists of dark matter with a nearly degenerate excited state, which allows for both elastic and inelastic scattering. In particular, the model includes a significant probability for particles to up-scatter from the ground state to the excited state. We simulate a suite of zoom-in Milky Way-sized N-body haloes with six models with different scattering cross sections to study the effects of up-scattering in SIDM models. We find that the up-scattering reaction greatly increases the central densities of the main halo through the loss of kinetic energy. However, the physical model still results in significant coring due to the presence of elastic scattering and down-scattering. These effects are not as apparent in the subhalo population compared to the main halo, but the number of subhaloes is reduced compared to CDM.

more » « less
Award ID(s):
2007355 2107724 2019786
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 288-306
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    Self-interacting dark matter (SIDM) is increasingly studied as a potential solution to small-scale discrepancies between simulations of cold dark matter (CDM) and observations. We examine a physically motivated two-state SIDM model with both elastic and inelastic scatterings. In particular, endothermic, exothermic, and elastic scattering have equal transfer cross-sections at high relative velocities ($v_{\rm rel}\gtrsim 400~{\rm km\, s}^{-1})$. In a suite of cosmological zoom-in simulation of Milky Way-size haloes, we vary the primordial state fractions to understand the impact of inelastic dark matter self-interactions on halo structure and evolution. In particular, we test how the initial conditions impact the present-day properties of dark matter haloes. Depending on the primordial state fraction, scattering reactions will be dominated by either exothermic or endothermic effects for high and low initial excited state fractions, respectively. We find that increasing the initial excited fraction reduces the mass of the main halo, as well as the number of subhaloes on all mass scales. The main haloes are cored, with lower inner densities and higher outer densities compared with CDM. Additionally, we find that the shape of the main halo becomes more spherical the higher the initial excited state fraction is. Finally, we show that the number of satellites steadily decreases with initial excited state fraction across all satellite masses.

    more » « less
  2. 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
  3. ABSTRACT We combine the isothermal Jeans model and the model of adiabatic halo contraction into a semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. The model agrees well with cosmological SIDM simulations over the entire core-forming stage up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable feature in the context of the structural diversity of bright dwarfs. The fast speed of the method facilitates analyses that would be challenging for numerical simulations – notably, we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, Mb/Mvir, and galaxy compactness, r1/2/Rvir; we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and we delineate the regime of core-collapse in the Mb/Mvir − r1/2/Rvir space, for a given cross section and concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible improvements of the fluid model. We have made our model publicly available at 
    more » « less

    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.

    more » « less
  5. ABSTRACT Core formation and runaway core collapse in models with self-interacting dark matter (SIDM) significantly alter the central density profiles of collapsed haloes. Using a forward modelling inference framework with simulated data-sets, we demonstrate that flux ratios in quadruple image strong gravitational lenses can detect the unique structural properties of SIDM haloes, and statistically constrain the amplitude and velocity dependence of the interaction cross-section in haloes with masses between 106 and 1010 M⊙. Measurements on these scales probe self-interactions at velocities below $30 \ \rm {km} \ \rm {s^{-1}}$, a relatively unexplored regime of parameter space, complimenting constraints at higher velocities from galaxies and clusters. We cast constraints on the amplitude and velocity dependence of the interaction cross-section in terms of σ20, the cross-section amplitude at $20 \ \rm {km} \ \rm {s^{-1}}$. With 50 lenses, a sample size available in the near future, and flux ratios measured from spatially compact mid-IR emission around the background quasar, we forecast $\sigma _{20} \lt 11\rm {\small {--}}23 \ \rm {cm^2} \rm {g^{-1}}$ at $95 {{\ \rm per\ cent}}$ CI, depending on the amplitude of the subhalo mass function, and assuming cold dark matter (CDM). Alternatively, if $\sigma _{20} = 19.2 \ \rm {cm^2}\rm {g^{-1}}$ we can rule out CDM with a likelihood ratio of 20:1, assuming an amplitude of the subhalo mass function that results from doubly efficient tidal disruption in the Milky Way relative to massive elliptical galaxies. These results demonstrate that strong lensing of compact, unresolved sources can constrain SIDM structure on sub-galactic scales across cosmological distances, and the evolution of SIDM density profiles over several Gyr of cosmic time. 
    more » « less