No Arabic abstract
We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.
We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions.
We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass sigma/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in sigma/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with sigma/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged)
We study the evolution of cosmological perturbations in dark-matter models with elastic and velocity-independent self interactions. Such interactions are imprinted in the matter-power spectrum as dark acoustic oscillations, which can be experimentally explored to determine the strength of the self scatterings. Models with self interactions have similarities to warm dark matter, as they lead to suppression of power on small scales when the dark-matter velocity dispersion is sizable. Nonetheless, both the physical origin and the extent of the suppression differ for self-interacting dark matter from conventional warm dark matter, with a dark sound horizon controlling the reduction of power in the former case, and a free-streaming length in the latter. We thoroughly analyze these differences by performing computations of the linear power spectrum using a newly developed Boltzmann code. We find that while current Lyman-$alpha$ data disfavor conventional warm dark matter with a mass less than 5.3 keV, when self interactions are included at their maximal value consistent with bounds from the Bullet Cluster, the limits are relaxed to 4.4 keV. Finally, we make use of our analysis to set novel bounds on light scalar singlet dark matter.
We present the first cosmological simulations of dwarf galaxies, which include dark matter self-interactions and baryons. We study two dwarf galaxies within cold dark matter, and four different elastic self-interacting scenarios with constant and velocity-dependent cross sections, motivated by a new force in the hidden dark matter sector. Our highest resolution simulation has a baryonic mass resolution of $1.8times 10^2,{rm M}_odot$ and a gravitational softening length of $34,{rm pc}$ at $z=0$. In this first study we focus on the regime of mostly isolated dwarf galaxies with halo masses $sim10^{10},{rm M}_odot$ where dark matter dynamically dominates even at sub-kpc scales. We find that while the global properties of galaxies of this scale are minimally affected by allowed self-interactions, their internal structures change significantly if the cross section is large enough within the inner sub-kpc region. In these dark-matter-dominated systems, self-scattering ties the shape of the stellar distribution to that of the dark matter distribution. In particular, we find that the stellar core radius is closely related to the dark matter core radius generated by self-interactions. Dark matter collisions lead to dwarf galaxies with larger stellar cores and smaller stellar central densities compared to the cold dark matter case. The central metallicity within $1,{rm kpc}$ is also larger by up to $sim 15%$ in the former case. We conclude that the mass distribution, and characteristics of the central stars in dwarf galaxies can potentially be used to probe the self-interacting nature of dark matter.
Observations show that supermassive black holes (SMBHs) with a mass of $sim10^9 M_odot$ exist when the Universe is just $6%$ of its current age. We propose a scenario where a self-interacting dark matter halo experiences gravothermal instability and its central region collapses into a seed black hole. The presence of baryons in protogalaxies could significantly accelerate the gravothermal evolution of the halo and shorten collapse timescales. The central halo could dissipate its angular momentum remnant via viscosity induced by the self-interactions. The host halo must be on high tails of density fluctuations, implying that high-$z$ SMBHs are expected to be rare in this scenario. We further derive conditions for triggering general relativistic instability of the collapsed region. Our results indicate that self-interacting dark matter can provide a unified explanation for diverse dark matter distributions in galaxies today and the origin of SMBHs at redshifts $zsim6-7$.