skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The onset of gravothermal core collapse in velocity-dependent self-interacting dark matter subhaloes
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
Award ID(s):
1814053 1814259
PAR ID:
10283695
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
505
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5327 to 5339
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. 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 https://github.com/JiangFangzhou/SIDM. 
    more » « less
  5. Abstract We study the evolution of isolated self-interacting dark matter (SIDM) halos that undergo gravothermal collapse and are driven deep into the short-mean-free-path regime.We assume spherical Navarro-Frenk-White (NFW) halos as initial conditions and allow for elastic dark matter self-interactions.We discuss the structure of the halo core deep in the core-collapsed regime and how it depends on the particle physics properties of dark matter, in particular, the velocity dependence of the self-interaction cross section. We find an approximate universality deep in this regime that allows us to connect the evolution in the short- and long-mean-free-path regimes, and approximately map the velocity-dependent self-interaction cross sections to constant ones for the full gravothermal evolution. We provide a semi-analytic prescription based on our numerical results for halo evolution deep in the core-collapsed regime.Our results are essential for estimating the masses of the black holes that are likely to be left in the core of SIDM halos. 
    more » « less