No Arabic abstract
We present N-body simulations of a new class of self-interacting dark matter models, which do not violate any astrophysical constraints due to a non-power-law velocity dependence of the transfer cross section which is motivated by a Yukawa-like new gauge boson interaction. Specifically, we focus on the formation of a Milky Way-like dark matter halo taken from the Aquarius project and re-simulate it for a couple of representative cases in the allowed parameter space of this new model. We find that for these cases, the main halo only develops a small core (~1 kpc) followed by a density profile identical to that of the standard cold dark matter scenario outside of that radius. Neither the subhalo mass function nor the radial number density of subhaloes are altered in these models but there is a significant change in the inner density structure of subhaloes resulting in the formation of a large density core. As a consequence, the inner circular velocity profiles of the most massive subhaloes differ significantly from the cold dark matter predictions and we demonstrate that they are compatible with the observational data of the brightest Milky Way dSphs in such a velocity-dependent self-interacting dark matter scenario. Specifically, and contrary to the cold dark matter case, there are no subhaloes that are more concentrated than what is inferred from the kinematics of the Milky Way dSphs. We conclude that these models offer an interesting alternative to the cold dark matter model that can reduce the recently reported tension between the brightest Milky Way satellites and the dense subhaloes found in cold dark matter simulations.
The development of methods and algorithms to solve the $N$-body problem for classical, collisionless, non-relativistic particles has made it possible to follow the growth and evolution of cosmic dark matter structures over most of the Universes history. In the best studied case $-$ the cold dark matter or CDM model $-$ the dark matter is assumed to consist of elementary particles that had negligible thermal velocities at early times. Progress over the past three decades has led to a nearly complete description of the assembly, structure and spatial distribution of dark matter haloes, and their substructure in this model, over almost the entire mass range of astronomical objects. On scales of galaxies and above, predictions from this standard CDM model have been shown to provide a remarkably good match to a wide variety of astronomical data over a large range of epochs, from the temperature structure of the cosmic background radiation to the large-scale distribution of galaxies. The frontier in this field has shifted to the relatively unexplored subgalactic scales, the domain of the central regions of massive haloes, and that of low-mass haloes and subhaloes, where potentially fundamental questions remain. Answering them may require: (i) the effect of known but uncertain baryonic processes (involving gas and stars), and/or (ii) alternative models with new dark matter physics. Here we present a review of the field, focusing on our current understanding of dark matter structure from $N$-body simulations and on the challenges ahead.
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 self-interacting dark matter (vdSIDM) in which the low velocity self-scattering cross-section, $sigma_{T}/m_{x}$, reaches 100 cm$^{2}$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, $V_{max}$, are significantly higher. We determine that these high $V_{max}$ subhaloes were objects in the mass range [$5times10^{6}$, $1times10^{8}$] $M_odot$ 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.
We investigate cosmological implications of an energy density contribution arising by elastic dark matter self-interactions. Its scaling behaviour shows that it can be the dominant energy contribution in the early universe. Constraints from primordial nucleosynthesis give an upper limit on the self-interaction strength which allows for the same strength as standard model strong interactions. Furthermore we explore the cosmological consequences of an early self-interaction dominated universe. Chemical dark matter decoupling requires that self-interacting dark matter particles are rather light (keV range) but we find that super-weak inelastic interactions are predicted by strong elastic dark matter self-interactions. Assuming a second, collisionless cold dark matter component, its natural decoupling scale exceeds the weak scale and is in accord with the electron and positron excess observed by PAMELA and Fermi-LAT. Structure formation analysis reveals a linear growing solution during self-interaction domination, enhancing structures up to ~ 10^(-3) solar masses long before the formation of the first stars.
We explore the phenomenology of having a second epoch of dark matter annihilation into dark radiation long after the standard thermal freeze-out. Such a hidden reannihilation process could affect visible sectors only gravitationally. As a concrete realization we consider self-interacting dark matter (SIDM) with a light force mediator coupled to dark radiation. We demonstrate how resonantly Sommerfeld enhanced cross sections emerge to induce the reannihilation epoch. The effect is a temporally local modification of the Hubble expansion rate and we show that the Cosmic Microwave Background (CMB) measurements -- as well as other observations -- have a high sensitivity to observe this phenomenon. Special attention is given to the model region where late kinetic decoupling and strong self-interactions can alleviate several small-scale problems in the cold dark matter paradigm at the same time. Interestingly, we find that reannihilation might here also simultaneously lower the tension between CMB and low-redshift astronomical observations of $H_{0}$ and $sigma_{8}$. Moreover, we identify reannihilation as a clear signature to discriminate between the phenomenologically otherwise almost identical vector and scalar mediator realizations of SIDM.
Self-interacting dark matter offers an interesting alternative to collisionless dark matter because of its ability to preserve the large-scale success of the cold dark matter model, while seemingly solving its challenges on small scales. We present here the first study of the expected dark matter detection signal taking into account different self-scattering models. We demonstrate that models with constant and velocity dependent cross sections, which are consistent with observational constraints, lead to distinct signatures in the velocity distribution, because non-thermalised features found in the cold dark matter distribution are thermalised through particle scattering. Depending on the model, self-interaction can lead to a 10% reduction of the recoil rates at high energies, corresponding to a minimum speed that can cause recoil larger than 300 km/s, compared to the cold dark matter case. At lower energies these differences are smaller than 5% for all models. The amplitude of the annual modulation signal can increase by up to 25%, and the day of maximum amplitude can shift by about two weeks with respect to the cold dark matter expectation. Furthermore, the exact day of phase reversal of the modulation signal can also differ by about a week between the different models. In general, models with velocity dependent cross sections peaking at the typical velocities of dwarf galaxies lead only to minor changes in the detection signals, whereas allowed constant cross section models lead to significant changes. We conclude that different self-interacting dark matter scenarios might be distinguished from each other through the details of direct detection signals. Furthermore, detailed constraints on the intrinsic properties of dark matter based on null detections, should take into account the possibility of self-scattering and the resulting effects on the detector signal.