This content will become publicly available on May 23, 2023

Core-collapse, evaporation, and tidal effects: the life story of a self-interacting dark matter subhalo
ABSTRACT Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass haloes, accelerated if haloes are subhaloes of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of ${\lesssim}10^8 \, \mathrm{M}_\odot$ subhaloes in lens-mass systems (${\sim}10^{13}\, \mathrm{M}_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhaloes with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhaloes (≲1/1000 host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We more »
Authors:
; ; ; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10348018
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
513
Issue:
4
Page Range or eLocation-ID:
4845 to 4868
ISSN:
0035-8711
1. 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 \ \rmmore » 2. 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,more » 3. 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 atmore » 4. ABSTRACT The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter haloes on sub-galactic scales. Using the image positions and flux ratios from eight quadruply imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhaloes and haloes along the line of sight, and account for warm dark matter free-streaming effects on the mass function and mass–concentration relation. By calibrating the scaling of the SHMF with host halo mass and redshift using a suite of simulated haloes, we infer a global normalization for the SHMF. We account for finite-size background sources, and marginalize over the mass profile of the main deflector. Parametrizing dark matter free-streaming through the half-mode mass mhm, we constrain the thermal relic particle mass mDM corresponding to mhm. At$95 \, {\rm per\, cent}$CI: mhm < 107.8 M⊙ ($m_{\rm {DM}} \gt 5.2 \ \rm {keV}$). We disfavour$m_{\rm {DM}} = 4.0 \,\rm {keV}$and$m_{\rm {DM}} = 3.0 \,\rm {keV}with likelihood ratios of 7:1 and 30:1, respectively, relative to the peak of the posterior distribution. Assuming cold dark matter, we constrainmore » 5. ABSTRACT We investigate the redshift evolution of the intrinsic alignments (IAs) of galaxies in the MassiveBlackII (MBII) simulation. We select galaxy samples above fixed subhalo mass cuts (M_h\gt 10^{11,12,13}\,\mathrm{M}_{\odot }\, h^{-1}) at z = 0.6 and trace their progenitors to z = 3 along their merger trees. Dark matter components of z = 0.6 galaxies are more spherical than their progenitors while stellar matter components tend to be less spherical than their progenitors. The distribution of the galaxy–subhalo misalignment angle peaks at ∼10 deg with a mild increase with time. The evolution of the ellipticity–direction (ED) correlation amplitude ω(r) of galaxies (which quantifies the tendency of galaxies to preferentially point towards surrounding matter overdensities) is governed by the evolution in the alignment of underlying dark matter (DM) subhaloes to the matter density of field, as well as the alignment between galaxies and their DM subhaloes. At scales\sim 1~\mathrm{Mpc}\, h^{-1}, the alignment between DM subhaloes and matter overdensity gets suppressed with time, whereas the alignment between galaxies and DM subhaloes is enhanced. These competing tendencies lead to a complex redshift evolution of ω(r) for galaxies at\sim 1~\mathrm{Mpc}\, h^{-1}$. At scales$\gt 1~\mathrm{Mpc}\, h^{-1}\$, alignment between DM subhaloes and matter overdensity does not evolve significantly; the evolutionmore »