Self-interaction Cross Section Research Articles

Brightest Cluster Galaxy Offsets in Cold Dark Matter

The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above , we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of . We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using the the “light traces mass” approach, in this context meaning gas is used as a tracer for the potential, can overestimate offsets by factors of ∼10 and ∼30, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section.

Self-Interacting Dark Matter Solves the Final Parsec Problem of Supermassive Black Hole Mergers.

Evidence for a stochastic gravitational wave (GW) background, plausibly originating from the merger of supermassive black holes (SMBHs), is accumulating with observations from pulsar timing arrays. An outstanding question is how inspiraling SMBHs get past the "final parsec" of separation, where they have a tendency to stall before GW emission alone can make the binary coalesce. We argue that dynamical friction from the dark matter (DM) spike surrounding the black holes is sufficient to resolve this puzzle, if the DM has a self-interaction cross section of order cm^{2}/g. The same effect leads to a softening of the GW spectrum at low frequencies as suggested by the current data. For collisionless cold DM, the friction deposits so much energy that the spike is disrupted and cannot bridge the final parsec, while for self-interacting DM, the isothermal core of the halo can act as a reservoir for the energy liberated from the SMBH orbits. A realistic velocity dependence, such as generated by the exchange of a massive mediator like a dark photon, is favored to give a good fit to the GW spectrum while providing a large enough core. A similar velocity dependence has been advocated for solving the small-scale structure problems of cold DM.

HSTPROMO Internal Proper-motion Kinematics of Dwarf Spheroidal Galaxies. I. Velocity Anisotropy and Dark Matter Cusp Slope of Draco

We analyze four epochs of Hubble Space Telescope imaging over 18 yr for the Draco dwarf spheroidal galaxy. We measure precise proper motions for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass–anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally averaged value βB¯=−0.20−0.53+0.28 , (where 1 – βB≡〈vtan2〉/〈vrad2〉 in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, Γdark, is −0.83−0.37+0.32 , consistent with the inner cusp predicted in ΛCDM cosmology. As expected given Draco’s low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, and 942 pc at 1σ, 2σ, and 3σ confidence, respectively, thus imposing important constraints on the self-interacting DM cross section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is 24.19−2.97+6.31kms−1 . We infer a dynamical distance of 75.37−4.00+4.73 kpc and show that Draco has a modest LOS rotation, with v/σ=0.22±0.09 . Our results provide a new stringent test of the so-called “cusp–core” problem that can be readily extended to other dwarfs.

On the late-time evolution of velocity-dependent self-interacting dark matter halos

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.

Simulations of galaxy cluster mergers with velocity-dependent, rare, and frequent self-interactions

ABSTRACT Self-interacting dark matter (SIDM) has been proposed to solve small-scale problems in $\rm {\Lambda CDM}$ cosmology. In previous work, constraints on the self-interaction cross-section of dark matter have been derived assuming that the self-interaction cross-section is independent of velocity. However, a velocity-dependent cross-section is more natural in most theories of SIDM. Using idealized N-body simulations without baryons, we study merging clusters with velocity-dependent SIDM. In addition to the usual rare scattering in the isotropic limit, we also simulate these systems with anisotropic, small-angle (frequent) scatterings. We find that the collisionless brightest cluster galaxy (BCG) has an offset from the DM peak that grows at later stages. Finally, we also extend the existing upper bounds on the velocity-independent, isotropic self-interaction cross-section to the parameter space of rare and frequent velocity-dependent self-interactions by studying the central densities of dark matter-only isolated haloes. For these upper-bound parameters, the DM-BCG offsets just after the first pericentre in the dark matter-only simulations are found to be ≲10 kpc. On the other hand, because of BCG oscillations, we speculate that the distribution of BCG offsets in a relaxed cluster is a statistically viable probe. Therefore, this motivates further studies of BCG off-centring in hydrodynamic cosmological simulations.

R-modes as a new probe of dark matter in neutron stars

In this work, we perform the first systematic investigation of effects of the presence of dark matter on r-mode oscillations in neutron stars (NSs). Using a self-interacting dark matter (DM) model based on the neutron decay anomaly and a hadronic model obtained from the posterior distribution of a recent Bayesian analysis, we impose constraints on the DM self-interaction strength using recent multimessenger astrophysical observations. We also put new constraints on the DM fraction for this model of DM. The constrained DM interaction strength is then used to estimate DM self-interaction cross section and shear viscosity resulting from DM, which is found to be several orders of magnitude smaller than shear viscosity due to hadronic matter. Assuming chemical equilibrium among DM fermions and neutrons, we estimate the bulk viscosity resulting from the dark decay of neutrons considering different scenarios for the temperature dependence of the reaction rate and investigate the effect on the r-mode instability window. We conclude that DM shear and bulk viscosity may significantly modify the r-mode instability window compared with the minimal hadronic viscosities, depending on the temperature dependence of the process. We also found that for the window to be compatible with the X-ray and pulsar observational data, the rate for the dark decay process must be fast.

Constraining self-interacting fermionic dark matter in admixed neutron stars using multimessenger astronomy

ABSTRACT We investigate the structure of admixed neutron stars with a regular hadronic component and a fraction of fermionic self-interacting dark matter. Using two limiting equations of state for the dense baryonic interior, constructed from piecewise generalized polytropes, and an asymmetric self-interacting fermionic dark component, we analyse different scenarios of admixed neutron stars depending on the mass of dark fermions mχ, interaction mediators mϕ, and self-interacting strengths g. We find that the contribution of dark matter to the masses and radii of neutron stars leads to tension with mass estimates of the pulsar J0453+1559, the least massive neutron star, and with the constraints coming from the GW170817 event. We discuss the possibilities of constraining dark matter model parameters g and y ≡ mχ/mϕ, using current existing knowledge on neutron star estimations of mass, radius, and tidal deformability, along with the accepted cosmological dark matter freeze-out values and self-interaction cross-section to mass ratio, σSI/mχ, fitted to explain Bullet, Abell, and dwarf galaxy cluster dynamics. By assuming the most restrictive upper limit, σSI/mχ < 0.1 cm2 g−1, along with dark matter freeze-out range values, the allowed g–y region is 0.01 ≲ g ≲ 0.1, with 0.5 ≲ y ≲ 200. For the first time, the combination of updated complementary restrictions is used to set constraints on self-interacting dark matter.

Anisotropic strong lensing as a probe of dark matter self-interactions

ABSTRACT Galaxy-scale strongly lensed systems have been shown to provide a unique technique for exploring the underlying physics of dark matter at sub-galactic scales. In the past, much attention was given to detecting and studying individual haloes in a strong lens system. In addition to the subhaloes, line-of-sight haloes contribute significantly to the small perturbations in lensed images. In prior work, we demonstrated that these line-of-sight haloes imprint a distinctive anisotropic signature and hence give rise to a detectable non-zero parity-even quadrupole moment in the effective convergence field’s two-point correlation function. In this study, we show that these line-of-sight haloes also produce a non-zero curl component of the effective deflection field with a parity-odd quadrupole moment of the two-point function. These multipole moments have the ability to statistically separate line-of-sight haloes from dark matter substructure. In this paper, we examine how these multipole moments evolve in the presence of warm dark matter and self-interacting dark matter in terms of central density evolution and dark matter halo abundance. Importantly, we show that these different multipole moments display exquisite sensitivity to both the amplitude and the velocity dependence of the dark matter self-interaction cross-section. Our approach opens the door for strong lensing observations to probe dark matter self-interaction over a broad range of relative velocities.

An excursion into the core of the cluster lens Abell 1689

ABSTRACT Abell 1689 is a well-studied cluster of galaxies and one of the largest gravitational lens systems ever observed. We have obtained a reconstruction of the cluster Abell 1689 using grale, a free-form lens inversion method that relies exclusively on the multiple image data. Non-inclusion of any data related to cluster member galaxies ensures an unbiased measure of the mass distribution, which is the most notable feature of free-form methods like grale. We used two different sets of multiple image systems from the available strong lensing data – one containing only the secure systems (107 images) and the other containing all available systems, only excluding some very non-secure systems (151 images). For the very well-constrained central ∼100 kpc region of the cluster, we made a detailed comparison of the grale reconstructed lensing mass and stellar mass retrieved by the Spectral Energy Distribution (SED) fitting software fast++. We found a light-unaccompanied mass peak in this region, whose existence, while tentative, is favoured by the distribution of nearby images that are local maxima in the Fermat potential. However, further tests, using different methodologies are needed to confirm the reality of this feature. If it is shown to be real, this light-unaccompanied mass peak is consistent with dark matter self-interaction cross-section σ ≲ 1 cm2 g−1, while being in tension with larger cross-sections.

Revisiting puffy dark matter with novel insights: partial wave analysis

We present a comprehensive study on the self-interaction cross-section of puffy dark matter (DM) particles, which have a significant intrinsic size compared to their Compton wavelength. For such puffy DM self-interaction cross-section in the resonant and classical regimes, our study demonstrates the significance of the Yukawa potential and the necessity of partial wave analysis: (i) Due to the finite-size effect of puffy DM particles, the new Yukawa potential of puffy DM is found to enlarge the Born-effective regime for the self-interaction cross-section, compared with the point-like DM; (ii) Our partial wave analysis shows that depending on the value of the ratio between Rχ (radius of a puffy DM particle) and 1/mϕ (force range), the three regimes (Born-effective, resonant and classical) for puffy DM self-interaction cross-section can be very different from the point-like DM; (iii) We find that to solve the small-scale anomalies via self-interacting puffy DM, the Born-effective and the resonant regimes exist for dwarf galaxies, while for the cluster and Milky Way galaxy the non-Born regime is necessary.

Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way

We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)–like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches ≈100 cm2 g−1 at relative velocities of ≈10 km s−1, motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos’ gravothermal evolution, demonstrating that deviations from the collisionless – relation can be used to select deeply core-collapsed halos, where is a halo’s maximum circular velocity, and is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈108 M ⊙ is deeply core collapsed in our SIDM model. Core-collapsed systems form ≈10% of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above ≈109 M ⊙ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies.

The local dark matter distribution in self-interacting dark matter halos

We study the effects of dark matter self-interactions on the local dark matter distribution in selected Milky Way-like galaxies in the EAGLE hydrodynamical simulations. The simulations were run with two different self-interacting dark matter models, a constant and velocity-dependent self-interaction cross-section. We find that the local dark matter velocity distribution of the Milky Way-like halos in the simulations with dark matter self-interactions and baryons are generally similar to those extracted from cold collisionless dark matter simulations with baryons. In both cases, the local dark matter speed distributions agree well with their best fit Maxwellian distributions. Including baryons in the simulations with or without dark matter self-interactions increases the local dark matter density and shifts the dark matter speed distributions to higher speeds. To study the implications for direct detection, we compute the dark matter halo integrals obtained directly from the simulations and compare them to those obtained from the best fit Maxwellian velocity distribution. We find that a Maxwellian distribution provides a good fit to the halo integrals of most halos, without any significant difference between the results of different dark matter self-interaction models.

Exploring the self interacting dark matter properties from low redshift observations

We calculate the shear (eta ) and bulk viscosity (zeta ) of Self Interacting Dark Matter (SIDM) fluid using the kinetic theory formalism. Using the astrophysical constraints on dark matter self-interaction cross section over mass sigma /m , we demonstrate that viscous SIDM fluid violates the lower bound on the ratio of shear viscosity to its entropy density, eta /mathfrak {s}=frac{1}{4pi }. Then, considering the eta /mathfrak {s} bound as universal, we derive a theoretical upper limit on the ratio of velocity average dark matter self interaction cross-section to its mass and also estimate an upper limit on SIDM mass. We report that mass of the SIDM particle should be sub-GeV scale. Furthermore, with the assumption of a power-law form of eta and zeta , we study its evolution in the light of low redshift observations. We show that at the large redshift, the SIDM viscosity is small, but at the small redshift, it becomes sufficiently large and contributes significantly to cosmic dissipation. As a consequence, viscous SIDM can explain the low redshift observations and also consistent with the standard cosmological prediction.

An asymmetric SIMP dark matter model

In this paper, we construct the first asymmetric strongly interacting massive particles (SIMP) dark matter (DM) model, where a new vector-like fermion and a new complex scalar both having nonzero chemical potentials can be asymmetric DM particles. After the spontaneous breaking of a U(1)D dark gauge symmetry, these two particles can have accidental ℤ4 charges making them stable. By adding one more complex scalar as a mediator between the SIMP DM, the relic density of DM is determined by 3 → 2 and two-loop induced 2 → 2 annihilations in this model. On the other hand, the SIMP DM can maintain kinetic equilibrium with the thermal bath until the DM freeze-out temperature via the new gauge interaction. Interestingly, this model can have a bouncing effect on DM, whereby the DM number density rises after the chemical freeze-out of DM. With this effect, the prediction of the DM self-interacting cross section in this model can be consistent with astrophysical observations, and the ratio of the DM energy density to the baryonic matter energy density can be explained by primordial asymmetries. We also predict the DM-electron elastic scattering cross section that can be used to test this model in future projected experiments.

Galactic satellite systems in CDM, WDM and SIDM

ABSTRACT We investigate the population of bright satellites ($M_{*} \ge 10^{5} \, \mathrm{M}_{\odot }$) of haloes of mass comparable to that of the Milky Way in cosmological simulations in which the dark matter (DM) is either cold, warm, or self-interacting (CDM, WDM, and SIDM, respectively). The nature of the DM gives rise to differences in the abundance and structural properties of field haloes. In WDM, the main feature is a reduction in the total number of galaxies that form, reflecting a suppression of low-mass DM haloes and lower galaxy formation efficiency compared to CDM. For SIDM, the changes are structural, restricted to the central regions of haloes and dependent on the assumed self-interaction cross-section. We also consider different baryonic subgrid physics models for galaxy formation, in which supernova gas blowouts can or cannot induce the formation of a core in dwarf galaxies. Overall, the inclusion of baryons lessen the differences in the halo properties in the different DM models compared to DM-only simulations. This affects the satellite properties at infall and therefore their subsequent tidal stripping and survival rates. None the less, we find slightly less concentrated satellite radial distributions as the SIDM cross-section increases. Unfortunately, we also find that the satellite populations in simulations with baryon-induced cores in CDM and WDM can mimic the results found in SIDM, making the satellite stellar mass and maximum circular velocity functions heavily degenerate on the assumed nature of the DM and the adopted subgrid modelling. These degeneracies preclude using the brightest satellites of the Milky Way to constrain the nature of DM.

New pathways to the relic abundance of vector-portal dark matter

We fully explore the thermal freezeout histories of a vector-portal dark matter model, in the region of parameter space in which the ratio of masses of the dark photon $A^{\prime}$ and dark matter $\chi$ is in the range $1 \lesssim m_{A^{\prime}}/m_{\chi} \lesssim 2$. In this region $2 \rightarrow 2$ and $3 \rightarrow 2$ annihilation processes within the dark sector, as well as processes that transfer energy between the dark sector and the Standard Model, play important roles in controlling the thermal freezeout of the dark matter. We carefully track the temperatures of all species, relaxing the assumption of previous studies that the dark and Standard Model sectors remain in thermal equilibrium throughout dark matter freezeout. Our calculations reveal a rich set of novel pathways which lead to the observed relic density of dark matter, and we develop a simple analytic understanding of these different regimes. The viable parameter space in our model provides a target for future experiments searching for light (MeV-GeV) dark matter, and includes regions where the dark matter self-interaction cross section is large enough to affect the small-scale structure of galaxies.

Self-interacting dark baryons

Using results from lattice QCD, it is possible to quantitatively design models of dark baryons leading to velocity-dependent self-interaction cross sections that match the values needed for solving small-scale structure problems of standard cold dark matter. However it is not obvious that the main dark matter component in such models will be nucleons rather than large nuclei, or dark pions or atoms, whose scattering properties would be different. We first identify the parameters of a dark SU(3) sector analogous to QCD---the confinement scale $\mathrm{\ensuremath{\Lambda}}$ and pion mass ${m}_{\ensuremath{\pi}}$---needed to reproduce desired self-interaction cross sections. Then we show that these values can generically be compatible with the absence of a sufficiently stable deuteron bound state, and hence leading to no heavier nuclei, thus establishing the consistency of the scenario for self-interacting dark nucleons. The range of dark photon masses needed to avoid dominant pion or atomic dark matter is determined, as well as allowed values for the kinetic mixing parameter. The dark proton might be detected directly in future searches, by dark photon exchange.

New thermal relic targets for inelastic vector-portal dark matter

We examine the vector-portal inelastic dark matter (DM) model with DM mass $m_\chi$ and dark photon mass $m_{A'}$, in the `forbidden dark matter' regime where $1 \lesssim m_{A'}/m_\chi \lesssim 2$, carefully tracking the dark sector temperature throughout freezeout. The inelastic nature of the dark sector relaxes the stringent cosmic microwave background (CMB) and self-interaction constraints compared to symmetric DM models. We determine the CMB limits on both annihilations involving excited states and annihilation into $e^+e^-$ through initial-state-radiation of an $A'$, as well as limits on the DM self-scattering, which proceeds at the one-loop level. The unconstrained parameter space serves as an ideal target for accelerator $A'$ searches, and provides a DM self-interaction cross section that is large enough to observably impact small-scale structure.

Constraints on dark matter self-interaction from the internal density profiles of X-COP galaxy clusters

The fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter α ≈ 0.18 in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if DM particles self-interact via additional forces, we expect the mass density profiles of DM halos to flatten in their central regions, thereby increasing the Einasto shape parameter. We measured the structural properties of 12 massive galaxy clusters from observations of their hot gaseous atmosphere, using the X-ray observatory XMM-Newton, and of the Sunyaev-Zel’dovich effect using the Planck all-sky survey. After removing morphologically disturbed systems, we measured Einasto shape parameters with mean ⟨α⟩=0.19 ± 0.03 and intrinsic scatter σα = 0.06, which is in close agreement with the prediction of the cold dark matter paradigm. We used cosmological hydrodynamical simulations of cluster formation with self-interacting DM (BAHAMAS-SIDM) to determine how the Einasto shape parameter depends on the self-interaction cross section. We used the fitted relation to turn our measurements of α into constraints on the self-interaction cross section, which imply σ/m < 0.19 cm2 g−1 (95% confidence level) at collision velocity vDM − DM ∼ 1000 km s−1. This is lower than the interaction cross section required for DM self-interactions to solve the core-cusp problem in dwarf spheroidal galaxies, unless the cross section is a strong function of velocity.

Self-interacting freeze-in dark matter in a singlet doublet scenario

We examine the non-thermal production of dark matter in a scalar extended singlet doublet fermion model where the lightest admixture of the fermions constitutes a suitable dark matter candidate. The dark sector is non-minimal with the MeV scale singlet scalar, which is stable in the Universe lifetime and can mediate the self-interaction for the multi-GeV fermion dark matter mitigating the small scale structure anomalies of the Universe. If the dark sector is strongly coupled to yield a velocity dependent large self-interaction cross section, it undergoes internal dark thermal equilibrium after freeze-in production. We essentially end up with suppressed relic abundance for the fermion dark matter in a traditional radiation dominated Universe. In contrast, the presence of a modified cosmological phase in the early era drives the fermion dark matter to satisfy nearly the whole amount of observed relic. It also turns out that the assumption of an unconventional cosmological history can allow the GeV scale dark matter to be probed at LHC from displaced vertex signature with improved sensitivity.