Self-interacting Dark Matter Research Articles

Freezing-in cannibal dark sectors

Self-Interacting Dark Matter models can successfully explain dark matter (DM) production through interactions confined within the dark sector. However, they often lack measurable experimental signals due to their secluded nature. Including a feeble interaction with the visible sector through a Higgs portal leads not only to potential detection avenues and richer thermal production dynamics, but also to a possible explanation of the initial dark sector population through the freeze-in mechanism. In this work we study, by solving the full system of coupled Boltzmann equations for the number densities and temperatures of all the involved states, three scenarios of this type where the DM is: a real scalar with broken ℤ2, a complex scalar with unbroken ℤ3, and a ℤ3 scalar with an additional scalar mediator. All of these models have viable dark matter candidates in a cannibal phase while having different detection profiles. We show that cosmological bounds can be either exacerbated or evaded by changing the dark sector interactions, leading to potential signatures in long-lived particle and indirect detection experiments.

Stealth dark matter spectrum using Laplacian Heaviside smearing and irreducible representations

We present nonperturbative lattice calculations in the quenched approximation of the low-lying meson and baryon spectrum of the SU(4) gauge theory with fundamental fermion constituents. This theory is one instance of stealth dark matter, a class of strongly coupled theories, where the lowest mass stable baryon is the dark matter candidate. This work constitutes the first milestone in the program to study stealth dark matter self-interactions. Here, we focus on reducing excited state contamination in the single-baryon channel by applying the Laplacian Heaviside method, as well as projecting our baryon operators onto the irreducible representations of the octahedral group. We compare our resulting spectrum to previous work involving Gaussian smeared nonprojected operators and find good agreement with reduced statistical uncertainties. We also present the spectrum of the low-lying odd-parity baryons for the first time. Published by the American Physical Society 2024

Dynamical friction in rotating ultralight dark matter galactic cores

Abstract Dynamical friction and stellar orbital motion in spiral galaxies with dark matter composed of ultralight bosons in the state of rotating Bose-Einstein condensate (BEC) are studied. It is found that the dynamical friction force is significantly affected by the topological charge of the vortex structure of the BEC core with the strongest effect at distances near the galactic center. It is also shown that the ultralight dark matter self-interaction plays an important role in studying the dynamical friction.

Constraints on dark matter from dynamical heating of stars in ultrafaint dwarfs. II. Substructure and the primordial power spectrum

There is a large and growing interest in observations of small-scale structure in dark matter. We propose a new way to probe dark matter structures in the ∼10–108M⊙ range. This allows us to constrain the primordial power spectrum over shorter distances scales than possible with direct observations from the CMB. For k in the range ∼10–1000 Mpc−1 our constraints on the power spectrum are orders of magnitude stronger than previous bounds. We also set some of the strongest constraints on dark matter isocurvature perturbations. Our method relies on the heating effect such dark matter substructures would have on the distribution of stars in an ultrafaint dwarf galaxy. Many models of inflation produce enhanced power at these short distance scales and can thus be constrained by our observation. Further, many dark matter models such as axion dark matter, self-interacting dark matter and dissipative dark matter, produce dense structures which could be constrained this way. Published by the American Physical Society 2024

Self-interacting scalar dark matter around binary black holes

Gravitational waves can provide crucial insights about the environments in which black holes live. In this work, we use numerical relativity simulations to study the behavior of self-interacting scalar (wavelike) dark matter clouds accreting onto isolated and binary black holes. We find that repulsive self-interactions smoothen the “spike” of an isolated black hole and saturate the density. Attractive self-interactions enhance the growth and result in more cuspy profiles, but can become unstable and undergo explosions akin to the superradiant bosenova that reduce the local cloud density. We quantify the impact of self-interactions on an equal-mass black hole merger by computing the dephasing of the gravitational-wave signal for a range of couplings. We find that repulsive self-interactions saturate the density of the cloud, thereby reducing the dephasing. For attractive self-interactions, the dephasing may be larger, but if these interactions dominate prior to the merger, the dark matter can undergo bosenova during the inspiral phase, disrupting the cloud and subsequently reducing the dephasing. Published by the American Physical Society 2024

The Imprint of Dark Matter on the Galactic Acceleration Field

Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, pulsar timing, and eclipsing binary stars in the solar neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density and infer the amount of dark matter (DM) by accounting for the mass in stars, gas, and dust. Leveraging FIRE-2 simulations of Milky Way–mass galaxies we compare vertical acceleration profiles between cold DM (CDM) and self-interacting DM (SIDM) with a constant cross section of 1 cm2 g−1 across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos manifest a more oblate shape and consistently exhibit higher local stellar and DM densities and steeper vertical acceleration gradients, up to 10%–30% steeper near the solar neighborhood. However, similar magnitude changes can arise from azimuthal variations in the baryonic components at a fixed radius and external influences like mergers, making it difficult to distinguish between CDM and SIDM using acceleration measurements in a single galaxy.

Dynamical friction from self-interacting dark matter

Context. Merging compact objects such as binary black holes provide a promising probe for the physics of dark matter (DM). The gravitational waves emitted during inspiral potentially allow one to detect DM spikes around black holes. This is because the dynamical friction force experienced by the inspiralling black hole alters the orbital period and thus the gravitational wave signal. Aims. The dynamical friction arising from DM can potentially differ from the collisionless case when DM is subject to self-interactions. This paper aims to understand how self-interactions impact dynamical friction. Methods. To study the dynamical friction force, we use idealised N-body simulations, where we include self-interacting dark matter. Results. We find that the dynamical friction force for inspiralling black holes would be typically enhanced by DM self-interactions compared to a collisionless medium (ignoring differences in the DM density). At lower velocities below the sound speed, we find that the dynamical friction force can be reduced by the presence of self-interactions. Conclusions. DM self-interactions have a significant effect on the dynamical friction for black hole mergers. Assuming the Chandrasekhar formula may underpredict the deceleration due to dynamical friction.

Tachyonic production of dark relics: classical lattice vs. quantum 2PI in Hartree truncation

We study the out-of-equilibrium production of non-minimally coupled self-interacting scalar dark matter during reheating using classical lattice simulations. The outcomes of the classical simulations are in qualitative agreement with the previous results obtained using the quantum 2PI approach in the Hartree truncation. In particular, the novel non-linear resonance found in the 2PI Hartee study is present also in the classical lattice simulations and can dominate the final dark matter yield. For the parameters considered, the difference in final value of the scalar two-point function between the two approaches is a factor of O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(1).

Effects of dark matter on the spontaneous scalarization in neutron stars

Dark matter, an important portion of compact objects, can influence different phenomena in neutron stars. The spontaneous scalarization in the scalar-tensor gravity has been proposed for neutron stars. Here, we investigate the spontaneous scalarization in dark matter admixed neutron stars. Applying the dark matter equations of state, we calculate the structure of scalarized neutron stars containing dark matter. The dark matter equations of state are based on observational data from the rotational curves of galaxies and the fermionic self-interacting dark matter. Our results verify that the spontaneous scalarization is affected by the dark matter pressure in neutron stars. Depending on the central density of scalarized dark matter admixed neutron stars, the dark matter pressure alters the central scalar field. The increase of dark matter pressure in low-density scalarized stars amplifies the central scalar field. However, the pressure of dark matter in high-density scalarized stars suppresses the central scalar field. Our calculations confirm that the stars in the merger event GW170817 and in the low-mass X-ray binary 4U 1820-30 can be scalarized dark matter admixed neutron stars.

A deep-learning algorithm to disentangle self-interacting dark matter and AGN feedback models

A deep-learning algorithm to disentangle self-interacting dark matter and AGN feedback models

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature.To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm2/g).When simulating these halos, we also vary the massresolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses.The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (106 particles per halo), we find that variations of order 10% in collapse time are still present.The results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

Gravothermal catastrophe in resonant self-interacting dark matter models

Gravothermal catastrophe in resonant self-interacting dark matter models

A quantitative comparison between velocity dependent SIDM cross-sections constrained by the gravothermal and isothermal models

ABSTRACT One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross-section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest – the effects of these processes can be highly degenerate with those of dark matter self-interactions at small scales. In this work, we select 9 isolated galaxies and brightest cluster galaxies (BCGs) with baryonic components small enough such that the baryonic gravitational potentials do not significantly influence the halo gravothermal evolution processes. We then constrain the parameters of Rutherford and Møller scattering cross-section models with the measured rotation curves and stellar kinematics through the gravothermal fluid formalism and isothermal method. Cross-sections constrained by the two methods are consistent at $1\sigma$ confidence level, but the isothermal method prefers cross-sections greater than the gravothermal approach constraints by a factor of $\sim 3$.

Asymmetric self-interacting dark matter with a canonical seesaw model

We study the possibility of generating dark matter (DM) and baryon asymmetry of the Universe (BAU) simultaneously in an asymmetric DM framework, which also alleviates the small-scale structure issues of cold DM. While the thermal relic of such self-interacting DM remains underabundant due to efficient annihilation into light mediators, a nonzero asymmetry in the dark sector can lead to the survival of the required DM in the Universe. The existence of a light mediator leads to the required self-interactions of DM at small scales while keeping DM properties similar to cold DM at large scales. It also ensures that the symmetric DM component annihilates away, leaving the asymmetric part in the spirit of cogenesis. The particle physics implementation is done in canonical seesaw models of light neutrino mass, connecting it to the origin of DM and BAU. In particular, we consider type-I and type-III seesaw origin of neutrino mass for simplicity and minimality of the field content. We show that the desired self-interactions and relic of DM together with BAU while satisfying relevant constraints lead to strict limits on DM mass O(GeV)≲MDM≲460 GeV. In spite of being a high-scale seesaw, the models remain verifiable in different experiments, including direct and indirect DM searches as well as colliders. Published by the American Physical Society 2024

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.

Dianoga SIDM: Galaxy cluster self-interacting dark matter simulations

Context. Self-interacting dark matter (SIDM) can tackle or alleviate small-scale issues within the cosmological standard model ΛCDM, and diverse flavours of SIDM can produce unique astrophysical predictions, resulting in different possible signatures which can be used to test these models with dedicated observations of galaxy clusters. Aims. This work aims to assess the impact of dark matter self-interactions on the properties of galaxy clusters. In particular, the goal is to study the angular dependence of the cross section by testing rare (large angle scattering) and frequent (small angle scattering) SIDM models with velocity-dependent cross sections. Methods. We re-simulated six galaxy cluster zoom-in initial conditions with a dark matter-only run and with full-physics set-up simulations that include a self-consistent treatment of baryon physics. We tested the dark matter-only setup and the full physics setup with the collisionless cold dark matter, rare self-interacting dark matter, and frequent self-interacting dark matter models. We then studied their matter density profiles as well as their subhalo population. Results. Our dark matter-only SIDM simulations agree with theoretical models, and when baryons are included in simulations, our SIDM models substantially increase the central density of galaxy cluster cores compared to full-physics simulations using collisionless dark matter. SIDM subhalo suppression in full-physics simulations is milder compared to the one found in the dark matter-only simulations because of the cuspier baryonic potential that prevents subhalo disruption. Moreover, SIDM with small-angle scattering significantly suppresses a larger number of subhaloes compared to large-angle scattering SIDM models. Additionally, SIDM models generate a broader range of subhalo concentration values, including a tail of more diffuse subhaloes in the outskirts of galaxy clusters and a population of more compact subhaloes in the cluster cores.

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.

Self-interacting dark matter implied by nano-Hertz gravitational waves

The self-interacting dark matter (SIDM) paradigm offers a potential solution to the small-scale structure problems faced by collisionless cold dark matter. This framework incorporates self-interactions among dark matter particles, typically mediated by a particle with a MeV-scale mass. Recent evidences of nano-Hertz gravitational waves from pulsar timing arrays (PTAs) such as NANOGrav, CPTA, EPTA, and PPTA suggest the occurrence of a first-order phase transition (FOPT) at a MeV-scale temperature. Considering the close proximity between these two scales, we propose that the mediator mass in the SIDM model originates from the spontaneous breaking of a U(1)′ symmetry, which is driven by the FOPT indicated by PTA data. Consequently, the alignment of these two scales is believed to be deeply connected by the same underlying physics. By extensively exploring the parameter space, remarkably, we find that the parameter space favored by SIDM just provides an explanation for the PTA data. Published by the American Physical Society 2024

Exploring the nature of dark matter with the extreme galaxy AGC 114905

i ) morphology, remains somewhat uncertain. We present unmatched ultra-deep optical imaging of AGC 114905 reaching surface brightness limits $ r,lim 32$ mag/arcsec$^2$ ($3 10 arcsec times 10 arcsec) obtained with the 10.4 m Gran Telescopio Canarias. With the new imaging, we characterise the galaxy's optical morphology, surface brightness, colours, and stellar mass profiles in great detail. The stellar disc has a similar extent to the i disc, presents spiral arms-like features, and shows a well-defined edge. Stars and gas have a similar morphology, and crucially, we find an inclination of $31 in agreement with the previous determinations. We revisit the rotation curve decomposition of the galaxy, and we explore different mass models in the context of CDM, self-interacting dark matter (SIDM), fuzzy dark matter (FDM) or Modified Newtonian dynamics (MOND). We find that the last does not fit the circular speed of the galaxy, while CDM only does so with dark halo parameters rarely seen in cosmological simulations. Within the uncertainties, SIDM and FDM remain feasible candidates to explain the observed kinematics of AGC 114905.