Self-interacting Dark Matter Models Research Articles

Motivations for a large self-interacting dark matter cross-section from Milky Way satellites

ABSTRACT We explore the properties of Milky Way (MW) subhaloes in self-interacting dark matter models for moderate cross-sections of 1–5 cm2 g−1 using high-resolution zoom-in N-body simulations. We include the gravitational potential of a baryonic disc and bulge matched to the MW, which is critical for getting accurate predictions. The predicted number and distribution of subhaloes within the host halo are similar for 1 and 5 cm2 g−1 models, and they agree with observations of MW satellite galaxies only if subhaloes with peak circular velocity over all time >7 km s−1 are able to form galaxies. We do not find distinctive signatures in the pericentre distribution of the subhaloes that could help distinguish the models. Using an analytical model to extend the simulation results, we are able to show that subhaloes in models with cross-sections between 1 and 5 cm2 g−1 are not dense enough to match the densest ultrafaint and classical dwarf spheroidal galaxies in the MW. This motivates exploring velocity-dependent cross-sections with values larger than 5 cm2 g−1 at the velocities relevant for the satellites such that core collapse would occur in some of the ultrafaint and classical dwarf spheroidals.

Modelling self-interacting dark matter substructures – I. Calibration with N-body simulations of a Milky-Way-sized halo and its satellite

ABSTRACT We study evolution of single subhaloes with their masses of ∼109 M⊙ in a Milky Way-sized host halo for self-interacting dark matter (SIDM) models. We perform dark-matter-only N-body simulations of dynamical evolution of individual subhaloes orbiting its host by varying self-scattering cross-sections (including a velocity-dependent scenario), subhalo orbits, and internal properties of the subhalo. We calibrate a gravothermal fluid model to predict time evolution in spherical mass density profiles of isolated SIDM haloes with the simulations. We find that tidal effects of SIDM subhaloes can be described with a framework developed for the case of collision-less cold dark matter (CDM), but a shorter typical time-scale for the mass loss due to tidal stripping is required to explain our SIDM simulation results. As long as the cross-section is less than $\sim \! 10\, \mathrm{cm}^2\,\mathrm{g}^{ -1}$ and initial states of subhaloes are set within a 2σ-level scatter at redshifts of ∼2 predicted by the standard ΛCDM cosmology, our simulations do not exhibit a prominent feature of gravothermal collapse in the subhalo central density for 10 Gyr. We develop a semi-analytic model of SIDM subhaloes in a time-evolving density core of the host with tidal stripping and self-scattering ram pressure effects. Our semi-analytic approach provides a simple, efficient, and physically intuitive prediction of SIDM subhaloes, but further improvements are needed to account for baryonic effects in the host and the gravothermal instability accelerated by tidal stripping effects.

Constraining SIDM with halo shapes: Revisited predictions from realistic simulations of early-type galaxies

ABSTRACT We study the effect of self-interacting dark matter (SIDM) and baryons on the shape of early-type galaxies (ETGs) and their dark matter haloes, comparing them to the predictions of the cold dark matter (CDM) scenario. We use five hydrodynamical zoom-in simulations of haloes hosting ETGs ($M_{\rm vir}\sim 10^{13}\, \mathrm{M}_{\odot }$ and $M_{*}\sim 10^{11}\, \mathrm{M}_{\odot }$), simulated in CDM and a SIDM model with constant cross-section of σT/mχ = 1 cm2g−1. We measure the 3D and projected shapes of the dark matter haloes and their baryonic content using the inertia tensor and compare our measurements to the results of three HST samples of gravitational lenses and Chandra and XMM–Newton X-ray observations. We find that the inclusion of baryons greatly reduces the differences between CDM and a SIDM, together with the ability to draw constraints based on shapes. Lensing measurements reject the predictions of CDM dark-matter-only simulations and prefer one of the hydro scenarios. When we consider the total sample of lenses, observational data prefer the CDM hydro scenario. The shapes of the X-ray emitting gas are compatible with observational results in both hydro runs, with CDM predicting higher elongations only in the very centre. Contrary to previous claims at the scale of elliptical galaxies, we conclude that both CDM and our SIDM model can still explain observed shapes once we include baryons in the simulations. Our results demonstrate that this is essential to derive realistic constraints and that new simulations are needed to confirm and extend our findings.

X-ray morphology of cluster-mass haloes in self-interacting dark matter

ABSTRACT We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq (1\!-\!3)\times 10^{15}{\, \rm M_\odot }$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\!-\!0.2R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma /m)= 0.5\!-\!1\, {\rm cm}^2\, {\rm g}^{-1}$ than CDM and SIDM with $(\sigma /m)=0.1\, {\rm cm}^2\, {\rm g}^{-1}$. Our statistical analysis indicates that the latter two models are disfavoured at the $68{{\ \rm per\ cent}}$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.

A critical assessment of solutions to the galaxy diversity problem

Galactic rotation curves exhibit a diverse range of inner slopes. Observational data indicates that explaining this diversity may require a mechanism that correlates a galaxy's surface brightness with the central-most region of its dark matter halo. In this work, we compare several concrete models that capture the relevant physics required to explain the galaxy diversity problem. We focus specifically on a Self-Interacting Dark Matter (SIDM) model with an isothermal core and two Cold Dark Matter (CDM) models with/without baryonic feedback. In contrast to the CDM case, the SIDM model can lead to the formation of an isothermal core in the halo, and is also mostly insensitive to baryonic feedback processes, which act on longer time-scales. Using rotation curves from 90 galaxies in the Spitzer Photometry & Accurate Rotation Curves (SPARC) catalog, we perform a comprehensive model comparison that addresses issues of statistical methodology from prior works. The best-fit halo models that we recover are consistent with standard CDM concentration-mass and abundance matching relations. We find that both the SIDM and feedback-affected CDM models are better than a CDM model with no feedback in explaining the rotation curves of low and high surface brightness galaxies in the sample. However, when compared to each other, there is no strong statistical preference for either the SIDM or the feedback-affected CDM halo model as the source of galaxy diversity in the SPARC catalog.

The Signatures of Self-interacting Dark Matter and Subhalo Disruption on Cluster Substructure

The abundance, distribution, and inner structure of satellites of galaxy clusters can be sensitive probes of the properties of dark matter. We run 30 cosmological zoom-in simulations with self-interacting dark matter (SIDM), with a velocity-dependent cross section, to study the properties of subhalos within cluster-mass hosts. We find that the abundance of subhalos that survive in the SIDM simulations are suppressed relative to their cold dark matter (CDM) counterparts. Once the population of disrupted subhalos—which may host orphan galaxies—is taken into account, satellite galaxy populations in CDM and SIDM models can be reconciled. However, even in this case, the inner structures of subhalos are significantly different in the two dark matter models. We study the feasibility of using the weak-lensing signal from the subhalo density profiles to distinguish between the cold and self-interacting dark matter while accounting for the potential contribution of orphan galaxies. We find that the effects of self-interactions on the density profile of subhalos can appear degenerate with subhalo disruption in CDM, when orphans are accounted for. With current error bars from the Subaru Hyper Suprime-Cam Strategic Program, we find that subhalos in the outskirts of clusters (where disruption is less prevalent) can be used to constrain dark matter physics. In the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time will give precise measurements of the weak-lensing profile and can be used to constrain σ T /m at the ∼1 cm2 g−1 level at v ∼ 2000 km s−1.

Effective Self-interaction of Dark Matter from Gravitational Scattering

I show that gravitational scattering of dark matter objects of mass ∼104 M ⊙ and speeds of ∼10 km s−1, provides the cross section per unit mass required in self-interacting dark matter models that alleviate the small-scale structure challenges to the collisionless cold dark matter model. For primordial objects of mass 104 M 4 M ⊙, moving at the velocity dispersion characteristic of dwarf galaxies, 10v 1 km s−1, the cross section per unit mass for gravitational scattering is . The steep decline in interaction with increasing velocity explains why self-interaction is not evident in data on massive galaxies and clusters of galaxies.

Primordial black hole merger rate in self-interacting dark matter halo models

We study the merger rate of primordial black holes (PBHs) in self-interacting dark matter (SIDM) halo models. To explore a numerical description for the density profile of SIDM halo models, we use the result of a previously performed simulation for SIDM halo models with $\sigma/m=10~{\rm cm^{2}g^{-1}}$. We also propose a concentration-mass-time relation that can explain the evolution of the halo density profile related to SIDM models. Furthermore, we investigate the encounter condition of PBHs that may have been randomly distributed in the medium of dark matter halos. Under these assumptions, we calculate the merger rate of PBHs within each halo considering SIDM halo models and compare the results with that obtained for cold dark matter (CDM) halo models. To do this, we employ the definition of the time after halo virialization as a function of halo mass. We indicate that SIDM halo models for $f_{\rm PBH}>0.32$ can generate sufficient PBH mergers in such a way that those exceed the one resulted from CDM halo models. By considering the spherical-collapse halo mass function, we obtain similar results for the cumulative merger rate of PBHs. Moreover, we calculate the redshift evolution of the PBH total merger rate. To determine a constraint on the PBH abundance, we study the merger rate of PBHs in terms of their fraction and masses and compare those with the black hole merger rate estimated by the Advanced LIGO (aLIGO)-Advanced Virgo (aVirgo) detectors during the third observing run. The results demonstrate that within the context of SIDM halo models, the merger rate of $10~M_{\odot}-10~M_{\odot}$ events can potentially fall within the aLIGO-aVirgo window. We also estimate a relation between the fraction of PBHs and their masses, which is well consistent with our findings.

Self-interacting dark matter in cosmology: accurate numerical implementation and observational constraints

This paper presents a systematic and accurate treatment of the evolution of cosmological perturbations in self-interacting dark matter models, for particles which decoupled from the primordial plasma while relativistic. We provide a numerical implementation of the Boltzmann hierarchies developed in a previous paper [JCAP, 09 (2020) 041] in a publicly available Boltzmann code and show how it can be applied to realistic DM candidates such as sterile neutrinos either under resonant or non-resonant production mechanisms, and for different field mediators. At difference with traditional fluid approximations — also known as a c eff-c vis parametrizations — our approach follows the evolution of phase-space perturbations under elastic DM interactions for a wide range of interaction models, including the effects of late kinetic decoupling. Finally, we analyze the imprints left by different self interacting models on linear structure formation, which can be constrained using Lyman-α forest and satellite counts. We find new lower bounds on the particle mass that are less restrictive than previous constraints.

Constraining self-interacting dark matter with dwarf spheroidal galaxies and high-resolution cosmological N -body simulations

We study the density structures of dark matter subhalos for both cold dark matter and self-interacting dark matter models using high-resolution cosmological $N$-body simulations. We quantify subhalo's central density at 150 pc from the center of each subhalo at the classical dwarf spheroidal and ultrafaint dwarf scales found in Milky-Way sized halos. By comparing them with observations, we find that the self-interacting scattering cross section of $\ensuremath{\sigma}/m<3\text{ }\text{ }{\mathrm{cm}}^{2}\text{ }{\mathrm{g}}^{\ensuremath{-}1}$ is favored. Due to the combination of hosts' tide and self-interactions, the central density of subhalos with small pericenter shows a noticeable difference between the cold and the self-interacting models, indicating that Milky-Way satellites with small pericenter are ideal sites to further constrain the nature of dark matter by future large spectroscopic surveys.

Super interacting dark sector: An improvement on self-interacting dark matter via scaling relations of galaxy clusters

Self-interacting dark matter is known as one of the most appropriate candidates for dark matter. Due to its excellent success in removing many astrophysical problems, particularly in small scale structure, studying this model has taken on added significance. In this paper, we focus on the results of two previously performed simulations of cluster sized halos with self-interacting dark matter and introduce a new function for the density profile of galaxy clusters, which can perfectly describe the result of these simulations. This density profile helps to find a velocity dispersion profile and also a relation between cluster mass and concentration parameter. Using these relations, we investigate two scaling relations of galaxy clusters, namely mass-velocity dispersion and mass-temperature relations. The scaling relations reveal that in the self-interacting dark matter model, halos are more massive than what the standard non-interacting model predicts for any fixed temperature. We also study the mass-temperature relation for a hybrid interacting model, which is a combination of self-interacting dark matter idea with another model of the dark sector in which dark matter particle mass is determined according to its interaction with dark energy. This super interacting dark sector (SIDS) model can change the mass-temperature relation to a modified form that has the same result as a non-interacting model. Finally, we provide quantitative expressions which can describe the constants of this interacting model with the value of cross-section per unit mass of dark matter particles.

A stringent upper limit on dark matter self-interaction cross-section from cluster strong lensing

ABSTRACT We analyse strongly lensed images in eight galaxy clusters to measure their dark matter density profiles in the radial region between 10 kpc and 150 kpc, and use this to constrain the self-interaction cross-section of dark matter (DM) particles. We infer the mass profiles of the central DM haloes, bright central galaxies, key member galaxies, and DM subhaloes for the member galaxies for all eight clusters using the qlens code. The inferred DM halo surface densities are fit to a self-interacting dark matter model, which allows us to constrain the self-interaction cross-section over mass σ/m. When our full method is applied to mock data generated from two clusters in the Illustris-TNG simulation, we find results consistent with no dark matter self-interactions as expected. For the eight observed clusters with average relative velocities of $1458_{-81}^{+80}$ km s−1, we infer $\sigma /m = 0.082_{-0.021}^{+0.027} \rm cm^2\, g^{ -1}$ and $\sigma /m \lt 0.13~ \rm cm^2\, g^{ -1}$ at the 95 per cent confidence level.

Testing extensions to ΛCDM on small scales with forthcoming cosmic shear surveys

ABSTRACT We investigate the constraining power of forthcoming Stage-IV weak lensing surveys (Euclid, lsst, and NGRST) for extensions to the Lambda cold dark matter model on small scales, via their impact on the cosmic shear power spectrum. We use high-resolution cosmological simulations to calculate how warm dark matter (WDM), self-interacting dark matter (SIDM), and a running of the spectral index affect the non-linear matter power spectrum, P(k), as a function of scale and redshift. We evaluate the cosmological constraining power using synthetic weak lensing observations derived from these power spectra and that take into account the anticipated source densities, shape noise, and cosmic variance errors of upcoming surveys. We show that upcoming Stage-IV surveys will be able to place useful, independent constraints on both WDM models (ruling out models with a particle mass of ≲0.5 keV) and SIDM models (ruling out models with a velocity-independent cross-section of ≳10 cm2 g−1) through their effects on the small-scale cosmic shear power spectrum. Similarly, they will be able to strongly constrain cosmologies with a running spectral index. Finally, we explore the error associated with the cosmic shear cross-spectrum between tomographic bins, finding that it can be significantly affected by Poisson noise (the standard assumption is that the Poisson noise cancels between tomographic bins). We provide a new analytic form for the error on the cross-spectrum that accurately captures this effect.

Constraining self-interacting dark matter with the full dataset of PandaX-II

Self-interacting Dark Matter (SIDM) is a leading candidate proposed to solve discrepancies between predictions of the prevailing cold dark matter theory and observations of galaxies. Many SIDM models predict the existence of a light force carrier that mediate strong dark matter self-interactions. If the mediator couples to the standard model particles, it could produce characteristic signals in dark matter direct detection experiments. We report searches for SIDM models with a light mediator using the full dataset of the PandaX-II experiment, based on a total exposure of 132 tonne-days. No significant excess over background is found, and our likelihood analysis leads to a strong upper limit on the dark matter-nucleon coupling strength. We further combine the PandaX-II constraints and those from observations of the light element abundances in the early universe, and show that direct detection and cosmological probes can provide complementary constraints on dark matter models with a light mediator.

The central densities of Milky Way-mass galaxies in cold and self-interacting dark matter models

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.

X-ray shapes of elliptical galaxies and implications for self-interacting dark matter

Several proposed models for dark matter posit the existence of self-interaction processes that can impact the shape of dark matter halos, making them more spherical than the ellipsoidal halos of collisionless dark matter. One method of probing the halo shapes, and thus the strength of the dark matter self-interaction, is by measuring the shape of the X-ray gas that traces the gravitational potential in relaxed elliptical galaxies. In this work we identify a sample of 11 relaxed, isolated elliptical galaxies and measure the ellipticity of the gravitating matter using X-ray images from the XMM-Newton and Chandra telescopes. We explore a variety of different mass configurations and find that the dark matter halos of these galaxies have ellipticities around ϵ ≃ 0.2-0.5. While we find non-negligible scatter in the ellipticity distribution, our results are consistent with some degree of self-interaction at the scale of σ/m ∼ 1 cm2/g, yet they also remain compatible with a cold dark matter scenario. We additionally demonstrate how our results can be used to directly constrain specific dark matter models and discuss implications for current and future simulations of self-interacting dark matter models.

The surprising accuracy of isothermal Jeans modelling of self-interacting dark matter density profiles

ABSTRACT Recent claims of observational evidence for self-interacting dark matter (SIDM) have relied on a semi-analytic method for predicting the density profiles of galaxies and galaxy clusters containing SIDM. We present a thorough description of this method, known as isothermal Jeans modelling, and then test it with a large ensemble of haloes taken from cosmological simulations. Our simulations were run with cold and collisionless dark matter (CDM) as well as two different SIDM models, all with dark matter only variants as well as versions including baryons and relevant galaxy formation physics. Using a mix of different box sizes and resolutions, we study haloes with masses ranging from 3 × 1010 to $3 \times 10^{15} \mathrm{\, M_\odot }$. Overall, we find that the isothermal Jeans model provides as accurate a description of simulated SIDM density profiles as the Navarro–Frenk–White profile does of CDM haloes. We can use the model predictions, compared with the simulated density profiles, to determine the input DM–DM scattering cross-sections used to run the simulations. This works especially well for large cross-sections, while with CDM our results tend to favour non-zero (albeit fairly small) cross-sections, driven by a bias against small cross-sections inherent to our adopted method of sampling the model parameter space. The model works across the whole halo mass range we study, although including baryons leads to DM profiles of intermediate-mass ($10^{12} - 10^{13} \mathrm{\, M_\odot }$) haloes that do not depend strongly on the SIDM cross-section. The tightest constraints will therefore come from lower and higher mass haloes: dwarf galaxies and galaxy clusters.

Constraint on Hybrid Stars with Gravitational Wave Events

Motivated by the recent discoveries of compact objects from LIGO/Virgo observations, we study the possibility of identifying some of these objects as compact stars made of dark matter called dark stars, or the mix of dark and nuclear matters called hybrid stars. In particular, in GW190814, a new compact object with 2.6 M⊙ is reported. This could be the lightest black hole, the heaviest neutron star, and a dark or hybrid star. In this work, we extend the discussion on the interpretations of the recent LIGO/Virgo events as hybrid stars made of various self-interacting dark matter (SIDM) in the isotropic limit. We pay particular attention to the saddle instability of the hybrid stars which will constrain the possible SIDM models.

Self-interacting dark matter and the delay of supermassive black hole growth

ABSTRACT Using cosmological hydrodynamic simulations with physically motivated models of supermassive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass (Mvir ≈ 7 × 1011 M⊙ at z = 0) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm2 g−1. We find that SMBH formation is suppressed in the early Universe due to SIDM interactions. SMBH–SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyr of their evolution. Nonetheless, by z = 0.8 the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the MBH–M* relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.

The impact of inelastic self-interacting dark matter on the dark matter structure of a Milky Way halo

Abstract We study the effects of inelastic dark matter self-interactions on the internal structure of a simulated Milky Way (MW)-size halo. Self-interacting dark matter (SIDM) is an alternative to collisionless cold dark matter (CDM) which offers a unique solution to the problems encountered with CDM on sub-galactic scales. Although previous SIDM simulations have mainly considered elastic collisions, theoretical considerations motivate the existence of multi-state dark matter where transitions from the excited to the ground state are exothermic. In this work, we consider a self-interacting, two-state dark matter model with inelastic collisions, implemented in the Arepo code. We find that energy injection from inelastic self-interactions reduces the central density of the MW halo in a shorter timescale relative to the elastic scale, resulting in a larger core size. Inelastic collisions also isotropize the orbits, resulting in an overall lower velocity anisotropy for the inelastic MW halo. In the inner halo, the inelastic SIDM case (minor-to-major axis ratio s ≡ c/a ≈ 0.65) is more spherical than the CDM (s ≈ 0.4), but less spherical than the elastic SIDM case (s ≈ 0.75). The speed distribution f(v) of dark matter particles at the location of the Sun in the inelastic SIDM model shows a significant departure from the CDM model, with f(v) falling more steeply at high speeds. In addition, the velocity kicks imparted during inelastic collisions produce unbound high-speed particles with velocities up to 500 km s−1 throughout the halo. This implies that inelastic SIDM can potentially leave distinct signatures in direct detection experiments, relative to elastic SIDM and CDM.