Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SELF-INTERACTING DARK MATTER”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Mergers Matter: Gravothermal Collapse in Dwarf Halos with Self-Interacting Dark Matter

Self-Interacting Dark Matter (SIDM) models with large cross sections at relative velocities below $\sim100\,{\rm km \, s}^{-1}$ can be tested with dwarf galaxy observations. We analyze six dark-matter-only zoom-in $\sim10^{10}\,{\rm M}_\odot$ halos with diverse assembly histories, adopting a cross section over mass of $σ/m = 70\,cm^2 \, g^{-1}$. We find that mergers inject orbital kinetic energy into the halo, altering the heat transport and the gravothermal evolution of the core. Three of the six halos -- those with the most quiescent merger histories -- show clear signs of core collapse in these simulations. Halos with sustained mergers do not collapse. Furthermore, merger-induced heat transport drives two non-collapsing halos to central densities well below the predictions of the gravothermal fluid model. These findings suggest a novel mechanism for producing dark-matter-deficient galaxies and expanding the diversity of rotation curves beyond what halo concentration alone predicts. Merger histories are thus essential for understanding central density distributions of dwarf galaxy halos in SIDM.

Silverman, Maya [UC, Irvine; Bohr Inst.] (ORCID:00↗
DOE OSTI2025

SIDM Concerto: Compilation and Data Release of Self-interacting Dark Matter Zoom-in Simulations

We present SIDM Concerto: 14 cosmological zoom-in simulations in cold dark matter (CDM) and self-interacting dark matter (SIDM) models based on the Symphony and Milky Way-est suites. SIDM Concerto includes one Large Magellanic Cloud– (LMC-) mass system (host mass ∼10 11 M ⊙ ), two Milky Way (MW) analogs (∼10 12 M ⊙ ), two group-mass hosts (∼10 13 M ⊙ ), and one low-mass cluster (∼10 14 M ⊙ ). Each host contains ≈2 × 10 7 particles and is run in CDM and one or more strong, velocity-dependent SIDM models. Our analysis of SIDM (sub)halo populations over seven subhalo mass decades reveals that (1) the fraction of core-collapsed isolated halos and subhalos peaks at a maximum circular velocity corresponding to the transition of the SIDM cross section from a v −4 to v 0 scaling; (2) SIDM subhalo mass functions are suppressed by ≈50% relative to CDM in LMC, MW, and group-mass hosts but are consistent with CDM in the low-mass cluster host; (3) subhalos’ inner density profile slopes, which are more diverse in SIDM than in CDM, are sensitive to both the amplitude and shape of the SIDM cross section. Our simulations provide a benchmark for testing SIDM predictions with astrophysical observations of field and satellite galaxies, strong lensing systems, and stellar streams. Data products are publicly available at doi:10.5281/zenodo.14933624.

dark matter↗
DOE OSTI2025

Novel challenges in tracking self-interacting dark matter subhalos

Cosmological N-body simulations are among the primary tools for studying structure formation in the Universe. Analyses of these simulations critically depend on accurately identifying and tracking dark matter subhalos over time. In recent years, several new algorithms have been developed to improve the accuracy and consistency of subhalo tracking in cold dark matter (CDM) simulations. These algorithms should be revisited in the context of new physics beyond gravity, which can modify the evolution and final properties of subhalo populations. In this work, we apply the particle-tracking-based subhalo finder Symfind to velocity-dependent self-interacting dark matter (SIDM) simulations with large cross section amplitudes to assess the performance of particle-tracking methods beyond the CDM paradigm. We find that the core-particle-tracking technique, which is key to the success of these algorithms in CDM, does not always yield accurate results in SIDM. In particular, the interplay between dark matter self-interactions and tidal stripping can cause the diffusion of core particles to larger radii, leading particle-tracking-based algorithms to prematurely lose track of SIDM subhalos. For massive core-expansion subhalos and core-collapse subhalos that experience close or repeated pericentric passages, a significant fraction of core particles can be lost, and particle-tracking-based finders such as Symfind offer no clear advantage over traditional methods that rely on identifying phase-space overdensities. On the other hand, for subhalos with large pericentric distances or fewer, more distant passages, Symfind tends to outperform. These differences depend sensitively on the cross section amplitude and turnover velocity of the underlying SIDM model. We therefore recommend a hybrid approach that leverages the strengths of both techniques to produce complete and robust catalogs of core-expansion and core-collapse SIDM subhalos.

cosmological simulations↗
DOE OSTI2025

Effects of Dark Matter Self-interactions on Sagittarius and its Stream

This work explores how assumptions regarding the particle-physics nature of dark matter can alter the evolution of the Sagittarius (Sgr) dwarf spheroidal galaxy and its expansive stellar stream. We run a large suite of N-body simulations to model the infall of a spherically symmetric Sgr-like dwarf, exploring how the presence of dark matter self-interactions impacts its evolution. For a scattering cross section of σ/m χ = 30 cm 2 g −1 (at orbital velocity scales), these interactions result in significantly less stellar mass and little to no dark matter bound to the progenitor at the present day. To isolate the cause of this mass loss, we introduce a novel technique for controlling which pairs of dark matter simulation particles can interact. This enables us to identify ram-pressure evaporation—the scattering of satellite and host dark matter particles—as the primary source of the enhanced mass loss. The rapid disintegration of the Sgr progenitor when self-interactions are allowed alters some key properties of the resulting stellar stream, most dramatically suppressing the presence of a “spur” on the apocenter of the trailing stream arm that correlates with the mass of the satellite at last pericenter. We demonstrate how the effects on the Sgr system scale with the particular choice of self-interaction cross section, which affects the degree of ram-pressure evaporation. These findings generalize beyond the Sgr system, underscoring that dwarf stellar streams and dwarf galaxies with close passages may serve as sensitive probes for dark matter self-interactions.

Milky Way dark matter halo↗
DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗
DOE OSTI2025

COZMIC. III. Cosmological Zoom-in Simulations of Self-interacting Dark Matter with Suppressed Initial Conditions

We present eight cosmological dark matter (DM)-only zoom-in simulations of a Milky Way–like system that include suppression of the linear matter power spectrum P(k), and/or velocity-dependent DM self-interactions, as the third installment of the COZMIC suite. We consider a model featuring a massive dark photon that mediates DM self-interactions and decays into massless dark fermions. The dark photon and dark fermions suppress linear matter perturbations, resulting in dark acoustic oscillations in P(k), which ultimately affect dwarf galaxy scales. The model also features a velocity-dependent elastic self-interaction between DM particles (SIDM), with a cross section that can alleviate small-scale structure anomalies. For the first time, our simulations test the impact of P(k) suppression on gravothermal evolution in an SIDM scenario that leads to core collapse in (sub)halos with present-day virial masses below ≈10 9 M ⊙ . In simulations with P(k) suppression and self-interactions, the lack of low-mass (sub)halos and the delayed growth of structure reduce the fraction of core-collapsed systems relative to SIDM simulations without P(k) suppression. In particular, P(k) suppression that saturates current warm DM constraints almost entirely erases core collapse in isolated halos. Models with less extreme P(k) suppression produce core collapse in ≈20% of subhalos and ≈5% of isolated halos above 10 8 M ⊙ , and also increase the abundance of extremely low-concentration isolated low-mass halos relative to SIDM. These results reveal a complex interplay between early and late-Universe DM physics, revealing new discovery scenarios in the context of upcoming small-scale structure measurements.

dark matter↗
DOE OSTI2025

Gravothermal collapse and the diversity of galactic rotation curves

The rotation curves of spiral galaxies exhibit a great diversity that challenges our understanding of galaxy formation and the nature of dark matter. Previous studies showed that in self-interacting dark matter (SIDM) models with a cross section per unit mass of σ/m ≈ $\mathscr{O}$(1) cm 2 /g, the predicted dark matter central densities are a good match to the observed densities in galaxies. Here, in this work, we explore a regime with a larger cross section of σ/m ≈ 20−40 cm 2 /g in dwarf galactic halos. We will show that such strong dark matter self-interactions can further amplify the diversity of halo densities inherited from their assembly history. High concentration halos can enter the gravothermal collapse phase within 10 Gyr, resulting in a high density, while low concentration ones remain in the expansion phase and have a low density. We fit the rotation curves of 14 representative low surface brightness galaxies and demonstrate how the large range of observed central densities is naturally accommodated in the strong SIDM regime of σ/m ≈ 20 – 40 cm 2 /g. Galaxies that are outliers in the previous studies, due to their high halo central densities, are no longer outliers in this SIDM regime as their halos would be in the collapse phase. For galaxies with a low density, the SIDM fits are robust to the variation of the cross section. Our findings open up a new window for testing gravothermal collapse, the unique signature of strong dark matter self-interactions, and exploring a broader SIDM model space. As an example, we illustrate how the larger cross sections favored by our fits, together with upper limits from strong lensing observations in clusters, pick out the preferred SIDM model space for a dark matter particle coupled to a light gauge boson in the Born regime.

dark matter↗
DOE OSTI2026

Using strong lensing to detect subhaloes with steep inner density profiles

ABSTRACT The inner region of a subhalo’s density distribution is particularly sensitive to dark matter microphysics, with alternative dark matter models leading to both cored and steeply-rising inner density profiles. This work investigates how the lensing signature and detectability of dark matter subhaloes in mock HST-, Euclid-, and JWST-like strong lensing observations depend on the subhalo’s radial density profile, especially with regards to the inner power-law slope, $\beta$. We demonstrate that the minimum subhalo mass detectable along the Einstein ring of a system is strongly dependent on $\beta$. In particular, we show that subhaloes with $\beta = 2.2$ can be detected down to masses over an order-of-magnitude lower than their Navarro–Frenk–White (NFW) counterparts with $\beta = 1$. Importantly, we find that the detectability of subhaloes with steep inner profiles is minimally affected by increasing the complexity of the main lens galaxy’s mass model. This is a notable characteristic of these subhaloes, as those with NFW or shallower profiles become essentially undetectable when multipole perturbations are added to the lens model. The results of this work highlight how the underlying dark matter physics can significantly impact the expected number of subhalo detections from strong gravitational lensing observations. This is important for testing Cold Dark Matter against alternative models, such as Self-Interacting Dark Matter, that predict a diverse range of subhalo inner density profiles.

dark matter↗
DOE OSTI2026

Massive black holes seeded by dark matter: Implications for little red dots and gravitational wave signatures

Observations of supermassive black holes (SMBHs) at high redshifts challenge standard seeding scenarios. We examine a dissipative self-interacting dark matter (dSIDM) model in which gravothermal collapse leads to the formation of massive BH seeds ab initio. We utilize a semianalytical framework to predict properties of the dSIDM-seeded SMBH population. Billion solar mass quasars are reproduced along with low-mass faint active galactic nuclei (known as little red dots) with SMBH-to-galaxy stellar mass ratios consistent with recent James Webb Space Telescope observations. To match the abundance of the observed bright quasars, a subpercent-level duty cycle is suggested, implying a large population of dormant SMBHs. The gravitational wave (GW) signals from mergers of these massive SMBHs can be detected by LISA while remaining within the NANOGrav constraints on the GW background. These results provide testable signatures of DM-driven SMBH formation, offering a pathway to probe hidden-sector physics through SMBH and GW observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Spin-dependent scattering of sub-GeV dark matter: Models and constraints

We calculate the scattering rate of sub-GeV dark matter in solid-state targets for spin-dependent dark matter–nucleon interactions. For dark matter particles with mass below 100 MeV, the scattering occurs predominantly through incoherent phonon production. For dark matter heavier than 100 MeV, we match onto the nuclear recoil calculation. To compare the sensitivity of future direct detection experiments with existing constraints, we consider three models with interactions that are mediated by spin-0 or spin-1 particles. This allows us to derive bounds on the cross section from searches for the mediating particle, including bounds from stellar cooling, beam dump experiments, meson factories, and dark matter self-interactions. The existing bounds are very stringent, though for 𝑚 𝜒 ≳100 MeV there is parameter space, which may be accessible with direct detection, depending on the exposure and background rates.

Gori, Stefania [University of California, Santa Cr↗
DOE OSTI2026

SEARCHING FOR MESONIC DARK MATTER WITH THE HEAVY PHOTON SEARCH EXPERIMENT

Several highly-sensitive astrophysical experiments over the past couple of decades have demonstrated that the current abundance of visible Standard Model matter cannot explain galactic rotation curves, the expansion history of the Universe, or the apparent warping of light in empty space. Instead, one finds strong agreement with this body of experimental results upon positing the existence of an invisible particulate field, dark matter. Namely, a cold, weakly interacting dark matter component can explain all these phenomena. A number of accelerator-based experiments have been developed to search for the weak couplings/interactions of these particles, many of them concentrating on particle models with masses of tens to thousands of GeV. A relatively new, well-motivated model is a dark sector coupled to the Standard Model via a dark photon. The current abundance of dark matter can be obtained if one assumes that dark matter is coupled to light by a MeV to GeV particle with a U(1) symmetry. The parameter space of these models remains largely unexplored because they are difficult to probe experimentally. In this thesis, I analyze data from the Heavy Photon Search (HPS) detector, whose two detector halves closely surround the electron beam, providing acceptance to far-forward boosted interactions. This forward acceptance to highly boosted particles yields unprecedented sensitivity to MeV-scale invariant masses. I exhaustively optimize the offline reconstruction of the HPS detector. Each reconstruction object, from Silicon Vertex Tracker hits to tracks, is studied to maximize acceptance of dark matter events. I then use the 2021 run data to search for one model of dark-photon-mediated matter, the Strongly Interacting Massive Particle (SIMP). SIMP models provide self-interacting dark matter candidates that can form bound states resembling dark mesons. HPS can detect SIMPs through the decay of a dark vector boson (either a dark ¿ or ¿) into e+e- pairs. I obtain exclusion contours for SIMPs using both an optimized cuts-based selection and a machine-learning-based selection, advancing our knowledge of the nature of dark matter.

O'Dwyer, Rory [Stanford Univ., CA (United States).↗
DOE OSTI2026

Self-interacting sub-GeV dark matter with strong MeV gamma-ray

Sub-GeV dark matter (DM) with s-channel resonant self-scattering provides a promising framework for addressing small-scale structure problems. However, models that also account for the observed relic abundance through the same resonance are strongly constrained by current γ-ray observations, since the associated signals are significantly enhanced. To overcome this limitation, we propose a framework in which the relic abundance and self-scattering are governed independently by two distinct mediators. As a concrete realization, we present a singlet scalar DM model in which self-scattering is mediated by a vector boson associated with a gauged baryon number, while the relic density is determined by forbidden annihilation into dark Higgs bosons that generate the gauge boson mass. By imposing cosmological, experimental, and theoretical constraints, We identify viable parameter regions that reproduce the observed relic density, alleviate small-scale problems, and remain consistent with current bounds. Notably, the model predicts multiple distinctive MeV γ-ray signals, a significant fraction of which will be testable with next-generation MeV γ-ray telescopes, including the Compton Spectrometer and Imager (COSI).

Models for Dark Matter↗
DOE OSTI2026

Naturally resonant dark matter from extra dimensions

We explore the mass-resonance structure that naturally arises in extradimensional models. The resonance can enhance both dark matter annihilation and self-interactions. We demonstrate such a resonance structure by considering the fermionic dark matter and dark photon models on an S 1 / ( Z 2 × Z 2 ′ ) orbifold. We also note that this model embeds a dark matter axial-vector coupling to the dark photon, thereby opening up the viable dark matter parameter space. We then present the unique predictions for direct-detection experiments and accelerator searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Plausible constraints and inflationary production for dark photons

Generic constraints on dark photons are generally presented assuming they have Stueckelberg masses. These constraints are strengthened if instead the mass is due to the Higgs mechanism and the dark Higgs is light. First, we show that under reasonable assumptions on the origin of kinetic mixing 𝜀 and perturbativity, the strengthened constraints on 𝜀 cannot be arbitrarily relaxed by making the Higgs heavy. Second, we demonstrate a simple mechanism for generating dark photon dark matter after inflation, where fluctuations of a dark Higgs by stochastic misalignment can produce stable dark photons through ℎ →𝐴′⁢𝐴′ decay. Third, we point out new lower bounds on 𝜀 in the case where the dark photon mediates thermal freeze-out of light dark matter by 𝑠-channel exchange, taking account of generic expectations for the size of 𝜀, and astrophysical upper bounds on the self-interaction cross section.

inflation↗
DOE OSTI2026

Topological freeze-out by semi-annihilation

We point out that a QCD-like dark sector can be coupled to the Standard Model by gauging the topological Skyrme current, which measures the dark baryon number in the infrared, to give a technically natural model for dark matter. This coupling allows for a semi-annihilation process χχ → χX μ , where X μ is the gauge boson mediator and χ a dark pion field, which plays the dominant role in setting the dark matter relic abundance. The topological interaction is purely p-wave and so free from indirect detection constraints. We show that the dark matter pion mass needs to be in the range 10 MeV ≲ m χ ≲ 1 TeV; towards the lighter end of this range, there can moreover be significant self-interactions. We discuss prospects for probing this scenario at collider experiments, ranging from the LHC to low-energy e + e − colliders, future Higgs factories, and beam-dump experiments.

Models for Dark Matter↗
DOE OSTI2025

Non-Gaussianity from explicit U(1)-breaking interactions

We investigate primordial non-Gaussianity (NG) arising from the explicit U(1) symmetry-breaking interactions during inflation involving a nearly massless axial component of a complex scalar field P. We analyze the induced NG parameter f NL under scenarios where the axial field functions as either a curvaton or cold dark matter (CDM). In the curvaton framework, there is a conventional contribution to the local NG of f NL ≃ -O(1). Additional positive local NG can result from either the self-interactions of axial field fluctuations, their interactions with a light radial partner, or kinetic mixing with the inflaton via U(1) symmetry-breaking terms. We identify parameter regions where the interactions lead to cancellations, suppressing the overall local NG to |f loc NL | ≲ O(0.1). In the CDM scenario, these interactions enhance the NG in the isocurvature fluctuations. Moreover, interactions between the axial field and another light scalar, such as a curvaton, can generate O(1) curvature NG signals and significant mixed curvature-isocurvature NGs that are within the reach of future experiments with σ(f loc NL ) ∼ 1. We also explore the role of a heavy radial field in generating oscillating correlation signals, noting that such signals can dominate the shape of the mixed adiabatic-isocurvature bispectrum. In certain cases, an oscillatory isocurvature bispectrum signal may be observable in the future, aiding in distinguishing between certain types of the U(1)-breaking self-interactions of the axial field.

axions↗
DOE OSTI2025

Nonperturbative and perturbative dynamics of a light QCD axion: Dark matter and the strong 𝐶⁢𝑃 problem

Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶⁢𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶⁢𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.

Axions↗
DOE OSTI2025

The Atacama Cosmology Telescope: DR6 constraints on extended cosmological models

We use new cosmic microwave background (CMB) primary temperature and polarization anisotropy measurements from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) to test foundational assumptions of the standard cosmological model, ΛCDM, and set constraints on extensions to it. We derive constraints from the ACT DR6 power spectra alone, as well as in combination with legacy data from the Planck mission. To break geometric degeneracies, we include ACT and Planck CMB lensing data and baryon acoustic oscillation data from DESI Year-1. To test the dependence of our results on non-ACT data, we also explore combinations replacing Planck with WMAP and DESI with BOSS, and further add supernovae measurements from Pantheon+ for models that affect the late-time expansion history. We verify the near-scale-invariance (running of the spectral index dn s /d ln k = 0.0062 ± 0.0052) and adiabaticity of the primordial perturbations. Neutrino properties are consistent with Standard Model predictions: we find no evidence for new light, relativistic species that are free-streaming (N eff = 2.86 ± 0.13, which combined with astrophysical measurements of primordial helium and deuterium abundances becomes N eff = 2.89 ± 0.11), for non-zero neutrino masses (∑m ν < 0.089 eV at 95% CL), or for neutrino self-interactions. We also find no evidence for self-interacting dark radiation (N idr < 0.134), or for early-universe variation of fundamental constants, including the fine-structure constant (α EM /α EM,0 = 1.0043 ± 0.0017) and the electron mass (m e /m e,0 = 1.0063 ± 0.0056). Our data are consistent with standard big bang nucleosynthesis (we find Y p = 0.2312 ± 0.0092), the COBE/FIRAS-inferred CMB temperature (we find T CMB = 2.698 ± 0.016 K), a dark matter component that is collisionless and with only a small fraction allowed as axion-like particles, a cosmological constant (w = -0.986 ± 0.025), and the late-time growth rate predicted by general relativity (γ = 0.663 ± 0.052). We find no statistically significant preference for a departure from the baseline ΛCDM model. In fits to models invoking early dark energy, primordial magnetic fields, or an arbitrary modified recombination history, we find H 0 = 69.9 +0.8 -1.5 , 69.1 ± 0.5, or 69.6 ± 1.0 km/s/Mpc, respectively; using BOSS instead of DESI BAO data reduces the central values of these constraints by 1–1.5 km/s/Mpc while only slightly increasing the error bars. In general, models introduced to increase the Hubble constant or to decrease the amplitude of density fluctuations inferred from the primary CMB are not favored over ΛCDM by our data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗