Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dark matter distribution”

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

The DREAMS Project: Disentangling the Impact of Halo-to-Halo Variance and Baryonic Feedback on Milky Way Dark Matter Speed Distributions

Dark matter direct detection experiments require information about the local dark matter speed distribution to produce constraints on dark matter candidates, or infer their properties in the event of a discovery. In this paper, we analyze how the uncertainty in the dark matter speed distribution near the Sun is affected by baryonic feedback, halo-to-halo variance, and halo mass. To do so, we harness the statistical power of the new DREAMS cold dark matter simulation suite, which is comprised of 1024 zoom-in Milky Way–mass halos with varied initial conditions as well as cosmological and astrophysical parameters. Applying a normalizing flows emulator to these simulations, we find that the uncertainty in the local dark matter speed distribution is dominated by halo-to-halo variance and, to a lesser extent, uncertainty in host halo mass. Uncertainties in supernova and black hole feedback (from the IllustrisTNG model in this case) are negligible in comparison. Using the DREAMS suite, we present a state-of-the-art prediction for the dark matter speed distribution in the Milky Way. Although the standard halo model is contained within the uncertainty of this prediction, individual galaxies may have distributions that differ from it. Lastly, we apply our DREAMS results to the XENON1T experiment and demonstrate that the astrophysical uncertainties are comparable to the experimental ones, solidifying previous results in the literature obtained with a smaller sample of simulated Milky Way–mass halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dark Matter Velocity Distributions for Direct Detection: Astrophysical Uncertainties Are Smaller Than They Appear

The sensitivity of direct detection experiments depends on the phase-space distribution of dark matter near the Sun, which can be modeled theoretically using cosmological hydrodynamical simulations of Milky Way–like galaxies. However, capturing the halo-to-halo variation in the local dark matter speeds—a necessary step for quantifying the astrophysical uncertainties that feed into experimental results—requires a sufficiently large sample of simulated galaxies, which has been a challenge. In this Letter, we quantify this variation with nearly 100 Milky Way–like galaxies from the tng50 simulation, the largest sample to date at this resolution. Moreover, we introduce a novel phase-space scaling procedure that endows every system with a reference frame that accurately reproduces the local standard-of-rest speed of our Galaxy, providing a principled way of extrapolating the simulation results to real-world data. The ensemble of predicted speed distributions is well characterized by the standard halo model, a Maxwell-Boltzmann distribution truncated at the escape speed, though the individual distributions can deviate from it, especially at high speeds. The dark matter–nucleon cross section limits placed by these speed distributions vary by ∼ 60% about the median. This places the 1⁢𝜎 astrophysical uncertainty at or below the level of the systematic uncertainty of current ton-scale detectors, even down to the energy threshold. The predicted uncertainty remains unchanged when subselecting on those TNG 50 galaxies with merger histories similar to the Milky Way. Tabulated speed distributions, as well as Maxwell-Boltzmann fits, are provided for use in computing direct detection bounds or projecting sensitivities.

Milky Way↗
DOE OSTI2022

Evaluating Lyman- α constraints for general dark-matter velocity distributions: Multiple scales and cautionary tales

One of the most important observational constraints on possible models of dark-matter physics exploits the Lyman-α absorption spectrum associated with photons traversing the intergalactic medium. Because these data allow us to probe the linear matter power spectrum with great accuracy down to relatively small distance scales, finding ways of accurately evaluating such Lyman-α constraints across large classes of candidate models of dark-matter physics is of paramount importance. While such Lyman-α constraints have been evaluated for dark-matter models that give rise to relatively simple dark-matter velocity distributions, more complex models—particularly those whose dark-matter velocity distributions stretch across multiple scales—have recently been receiving increasing attention. In this paper, we undertake a study of the Lyman-α constraints associated with general dark-matter velocity distributions. Although such Lyman-α constraints are difficult to evaluate in principle, in practice there currently exist two classes of methods in the literature through which such constraints can be recast into forms which are easier to evaluate and which therefore allow a more rapid determination of whether a given dark-matter model is ruled in or out. Accordingly, we utilize both of these recasts in order to determine the Lyman-α bounds on different classes of dark-matter velocity distributions. Here we also develop a general method by which the results of these different recasts can be compared. For relatively simple dark-matter velocity distributions, we demonstrate that these two classes of recasts tend to align and give similar results. However, we find that the situation is far more complex for distributions involving multiple velocity scales: while these two classes of recasts continue to yield similar results within certain regions of parameter space, they nevertheless yield dramatically different results within precisely those regions of parameter space which are likely to be phenomenologically relevant. This, then, serves as a cautionary tale regarding the use of such recasts for complex dark-matter velocity distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

How to Build an Empirical Speed Distribution for Dark Matter in the Solar Neighborhood

The dark matter flux in a direct detection experiment depends on its local speed distribution. This distribution has been inferred from simulations of Milky Way–like galaxies, but such models serve only as proxies, given that no simulation directly captures the detailed evolution of our own Galaxy. This motivates alternative approaches that obtain this distribution directly from observations. In this work, we utilize 98 Milky Way analogues from the TNG50 simulation to develop and validate a procedure for inferring the dark matter speed distribution using the kinematics of nearby stars. We find that the dark matter that originated from old mergers, plus that from recent nonluminous accretions, is well described by a Maxwell–Boltzmann speed distribution centered at the local standard-of-rest velocity. Meanwhile, recently accreted dark matter from massive mergers has speeds that can be traced from the associated stellar debris of these events. The stellar populations systematically underestimate the velocity dispersion of their dark matter counterparts, but a simple kinematic boost brings the two into good alignment. Using the TNG50 host galaxies, we demonstrate that combining these two contributions provides an accurate reconstruction of the local dark matter speeds. As an application of the procedure to our own Galaxy, we utilize stellar kinematic data from Gaia to quantify how the dark matter remnants from the Milky Way’s last major merger impact its speed distribution in the solar neighborhood.

Shpigel, Tal [Princeton Univ., NJ (United States)]↗
DOE OSTI2023

Maximum entropy distributions of dark matter in ΛCDM cosmology

Context. Small-scale challenges to ΛCDM cosmology require a deeper understanding of dark matter physics. Aims. This paper aims to develop the maximum entropy distributions for dark matter particle velocity (denoted by X ), speed (denoted by Z ), and energy (denoted by E ) that are especially relevant on small scales where system approaches full virialization. Methods. For systems involving long-range interactions, a spectrum of halos of different sizes is required to form to maximize system entropy. While the velocity in halos can be Gaussian, the velocity distribution throughout the entire system, involving all halos of different sizes, is non-Gaussian. With the virial theorem for mechanical equilibrium, we applied the maximum entropy principle to the statistical equilibrium of entire system, such that the maximum entropy distribution of velocity (the X distribution) could be analytically derived. The halo mass function was not required in this formulation, but it did indeed result from the maximum entropy. Results. The predicted X distribution involves a shape parameter α and a velocity scale, v 0 . The shape parameter α reflects the nature of force ( α → 0 for long-range force or α → ∞ for short-range force). Therefore, the distribution approaches Laplacian with α → 0 and Gaussian with α → ∞. For an intermediate value of α , the distribution naturally exhibits a Gaussian core for v ≪ v 0 and exponential wings for v ≫ v 0 , as confirmed by N -body simulations. From this distribution, the mean particle energy of all dark matter particles with a given speed, v , follows a parabolic scaling for low speeds (∝ v 2 for v ≪ v 0 in halo core region, i.e., “Newtonian”) and a linear scaling for high speeds (∝ v for v ≫ v 0 in halo outskirt, i.e., exhibiting “non-Newtonian” behavior due to long-range gravity). We compared our results against N -body simulations and found a good agreement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Phase space distribution functions and energy distributions of dark matter particles in haloes

ABSTRACT For a spherical dark matter halo with isotropic velocity distribution, the phase space distribution function (DF), the energy distribution, and the density profile form a set of self-consistent description of its equilibrium state, and knowing one is sufficient to determine the other two. The Navarro–Frenk–White density profile (NFW profile) is known to be a good approximation to the spherically averaged density distribution in simulated haloes. The DARKexp energy distribution is also known to compare well with the simulated energy distribution. We present a quantitative assessment of the NFW and DARKexp fits to the simulated DF and energy distribution for a wide range of haloes in a dark-matter-only simulation from the IllustrisTNG Project. As expected, we find that the NFW fits work well except at low energy when the density at small radii deviates from the NFW profile. Further, the NFW and DARKexp fits have comparable accuracy in the region where both fit well, but the DARKexp fits are better at low energy because they require matching of the central gravitational potential. We also find an approximate relation between the energy scale parametrizing the DARKexp energy distribution and that defined by the characteristic density and radius of the NFW profile. This relation may be linked to the relaxation process during halo formation.

Astronomy & Astrophysics↗
DOE OSTI2026

ELG×LRG Distribution through Dark Matter Halo Dynamics

We investigate the clustering and halo occupation distribution (HOD) of DESI Y1 emission-line (ELGs) and luminous red (LRGs) galaxies at 0.8 < z < 1.1, including their cross-correlation (ELG×LRG), using the A BACUS S UMMIT suite and a new Halo Occupation Model (H OME ) for galaxy multitracers. This integrates intrahalo dynamics, halo exclusion, and quenching, bridging insights from hydrodynamical, HOD, abundance-matching, and semianalytic studies. Leveraging full phase-space information from the Uchuu N-body simulation, and sampling satellites from dark-matter particle positions via physically motivated prescriptions, Home reproduces the anisotropic clustering down to s = 200 h −1 kpc with unprecedented accuracy. Model parameters are inferred solely from two-point statistics using a two-level Bayesian framework, yielding high-fidelity ELG, LRG, and cross-reference catalogs. We find that satellite ELGs behave as incoherent flows within their parent halos, dominating the clustering below 4 h −1 Mpc. The HOD from the best-fit Home has the following properties: (i) 90.50% (85.91%) of ELGs (LRGs) are central galaxies without satellites, residing in halos of M vir ∼ 6.6 × 10 11 (1.2 × 10 13 ) h −1 M ⊙ ; (ii) the ELG×LRG cross-correlation is governed by central-central pairs and shaped by halo exclusion on 2–5 h −1 Mpc scales; (iii) 9.50% (14.09%) of ELGs (LRGs) are satellites, of which 1.09% (3.52%) inhabit halos with a central galaxy of the same species in a maximally conformal configuration, 7.02% (0.005%) orbit complementary hosts in a minimally conformal state, and 0.58% (10.57%) are orphans. The high sensitivity of Home precisely captures the dynamics of satellites in different host environments, opening a promising avenue for understanding systematics and the dynamical nature of dark matter, potentially distinguishing gravity models.

Favole, Ginevra [Universidad de La Laguna (Spain);↗
DOE OSTI2025

Neutrino and gamma-ray signatures of inelastic dark matter annihilating outside neutron stars

We present a new inelastic dark matter search: neutron stars in dark matter-rich environments capture inelastic dark matter which, for interstate mass splittings between about 45 - 285 MeV, will annihilate away before becoming fully trapped inside the object. This means a sizable fraction of the dark matter particles can annihilate while being outsidethe neutron star, producing neutron star-focused gamma-rays and neutrinos. We analyze this effect for the first time and target the neutron star population in the Galactic Center, where the large dark matter and neutron star content makes this signal most significant. Depending on the assumed neutron star and dark matter distributions, we set constraints on the dark matter-nucleon inelastic cross-section using existing H.E.S.S. observations. We also forecast the sensitivity of upcoming gamma-ray and neutrino telescopes to this signal, which can reach inelastic cross-sections as low as ∼ 2 × 10 -47 cm 2 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Search for interacting galaxy clusters from SDSS DR-17 employing optimized friends-of-friends algorithm and multimessenger tracers

ABSTRACT In the theoretical framework of hierarchical structure formation, galaxy clusters evolve through continuous accretion and mergers of substructures. Cosmological simulations have revealed the best picture of the universe as a 3D filamentary network of dark-matter distribution called the cosmic web. Galaxy clusters are found to form at the nodes of this network and are the regions of high merging activity. Such mergers being highly energetic, contain a wealth of information about the dynamical evolution of structures in the Universe. Observational validation of this scenario needs a colossal effort to identify numerous events from all-sky surveys. Therefore, such efforts are sparse in literature and tend to focus on individual systems. In this work, we present an improved search algorithm for identifying interacting galaxy clusters and have successfully produced a comprehensive list of systems from SDSS DR-17. By proposing a set of physically motivated criteria, we classified these interacting clusters into two broad classes, ‘merging’ and ‘pre-merging/postmerging’ systems. Interestingly, as predicted by simulations, we found that most cases show cluster interaction along the prominent cosmic filaments of galaxy distribution (i.e. the proxy for dark matter filaments), with the most violent ones at their nodes. Moreover, we traced the imprint of interactions through multiband signatures, such as diffuse cluster emissions in radio or X-rays. Although we could not find direct evidence of diffuse emission from connecting filaments and ridges; our catalogue of interacting clusters will ease locating such faintest emissions as data from sensitive telescopes such as eROSITA or SKA, becomes accessible.

Oak, Tejas↗
DOE OSTI2023

Uncovering dark matter density profiles in dwarf galaxies with graph neural networks

Dwarf galaxies are small, dark-matter-dominated galaxies, some of which are embedded within the Milky Way. Their lack of baryonic matter (e.g., stars and gas) makes them perfect test beds for probing the properties of dark matter—understanding the spatial dark matter distribution in these systems can be used to constrain microphysical dark matter interactions that influence the formation and evolution of structures in our Universe. Here, we introduce a new method that leverages simulation-based inference and graph-based machine learning in order to infer the dark matter density profiles of dwarf galaxies from observable kinematics of stars gravitationally bound to these systems. Our approach aims to address some of the limitations of established methods based on dynamical Jeans modeling. We show that this novel method can place stronger constraints on dark matter profiles and, consequently, has the potential to weigh in on some of the ongoing puzzles associated with the small-scale structure of dark matter halos, such as the core-cusp discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamical friction in dark matter superfluids: The evolution of black hole binaries

Abstract The theory of superfluid dark matter is characterized by self-interacting sub-eV particles that thermalize and condense to form a superfluid core in galaxies. Massive black holes at the center of galaxies, however, modify the dark matter distribution and result in a density enhancement in their vicinity known as dark matter spikes. The presence of these spikes affects the evolution of binary systems by modifying their gravitational wave emission and inducing dynamical friction effects on the orbiting bodies. In this work, we assess the role of dynamical friction for bodies moving through a superfluid core enhanced by a central massive black hole. As a first step, we compute the dynamical friction force experienced by bodies moving in a circular orbit. Then, we estimate the gravitational wave dephasing of the binary, showing that the effect of the superfluid drag force is beyond the reach of space-based experiments like LISA, contrarily to collisionless dark matter, therefore providing an opportunity to distinguish these dark matter models.

Astronomy & Astrophysics↗
DOE OSTI2024

Slant, Fan, and Narrow: The Response of Stellar Streams to a Tilting Galactic Disk

Stellar streams are sensitive tracers of the gravitational potential, which is typically assumed to be static in the inner Galaxy. However, massive mergers like Gaia-Sausage-Enceladus can impart torques on the stellar disk of the Milky Way that result in the disk tilting at rates of up to 10°–20° Gyr –1 . Here, we demonstrate the effects of disk tilting on the morphology and kinematics of stellar streams. Through a series of numerical experiments, we find that streams with nearby apocenters (r apo ≲ 20 kpc) are sensitive to disk tilting, with the primary effect being changes to the stream's on-sky track and width. Interestingly, disk tilting can produce both more diffuse streams and more narrow streams, depending on the orbital inclination of the progenitor and the direction in which the disk is tilting. Our model of Pal 5's tidal tails for a tilting rate of 15° Gyr –1 is in excellent agreement with the observed stream's track and width, and reproduces the extreme narrowing of the trailing tail. We also find that failure to account for a tilting disk can bias constraints on shape parameters of the Milky Way's local dark matter distribution at the level of 5%–10%, with the direction of the bias changing for different streams. Disk tilting could therefore explain discrepancies in the Milky Way's dark matter halo shape inferred using different streams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

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

Abstract 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 cm 2 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 R max – V max relation can be used to select deeply core-collapsed halos, where V max is a halo’s maximum circular velocity, and R max is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈10 8 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 ≈10 9 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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Enhancing direct detection of Higgsino dark matter

While much supersymmetric weakly interacting massive particle (WIMP) parameter space has been ruled out, one remaining important candidate is Higgsino dark matter. The Higgsino can naturally realize the “inelastic dark matter” scenario, where the scattering off a nucleus occurs between two nearly-degenerate states, making it invisible to WIMP direct detection experiments if the splitting is too large to be excited. It was realized that a “luminous dark matter” detection process, where the Higgsino upscatters in the Earth and subsequently decays into a photon in a large neutrino detector, offers the best sensitivity to such a scenario. We consider the possibility of adding a large volume of a heavy element, such as Pb or U, around the detector. We also consider the presence of U and Th in the Earth itself, and the effect of an enhanced high-velocity tail of the dark matter distribution due to the presence of the Large Magellanic Cloud. These effects can significantly improve the sensitivity of detectors such as JUNO, SNO + , KamLAND, and Borexino, potentially making it possible in the future to cover much of the remaining parameter space for this classic supersymmetric WIMP dark matter. Published by the American Physical Society 2025

Graham, Peter W. (ORCID:0000000216001601)↗
DOE OSTI2025

The Dark Matter Content of Milky Way Dwarf Spheroidal Galaxies: Draco, Sextans, and Ursa Minor

The Milky Way Survey of the Dark Energy Spectroscopic Instrument (DESI) has so far observed three classical dwarf spheroidal galaxies (dSphs): Draco, Sextans, and Ursa Minor. Based on the observed line-of-sight velocities and metallicities of their member stars, we apply the axisymmetric Jeans Anisotropic Multi-Gaussian Expansion modeling (JAM) approach to recover their inner dark matter distributions. In particular, both the traditional single-population Jeans model and the multiple population chemodynamical model are adopted. With the chemodynamical model, we divide member stars of each dSph into metal-rich and metal-poor populations. The metal-rich populations are more centrally concentrated and dynamically colder, featuring lower velocity dispersion profiles than the metal-poor populations. We find a diversity of the inner density slopes γ of dark matter halos, with the best constraints by the single-population or chemodynamical models consistent with each other. The inner density slopes are $0.7{1}_{-0.35}^{+0.34}$, $0.2{6}_{-0.12}^{+0.22}$, and $0.3{3}_{-0.16}^{+0.20}$ for Draco, Sextans, and Ursa Minor, respectively. We also present the measured astrophysical J and D factors of the three dSphs. Our results indicate that the study of the dark matter content of dSphs through stellar kinematics is still subject to uncertainties behind both the methodology and the observed data, through comparisons with previous measurements and datasets.

Dark matter↗
DOE OSTI2025

Gravitational waves and galaxies cross-correlations: a forecast on GW biases for future detectors

ABSTRACT Gravitational waves (GWs) have rapidly become important cosmological probes since their first detection in 2015. As the number of detected events continues to rise, upcoming instruments like Einstein Telescope (ET) and Cosmic Explorer (CE) will observe millions of compact binary mergers. These detections, coupled with galaxy surveys by instruments such as the Dark Spectroscopic Energy Instrument (DESI), Euclid, and the Vera Rubin Observatory, will provide unique information on the large-scale structure of the universe by cross-correlating GWs with the distribution of galaxies hosting them. In this paper, we focus on how cross-correlations constrain the clustering bias of GWs emitted by the coalescence of binary black holes (BBHs). This parameter links BBHs to the underlying dark matter distribution, hence informing us how they populate galaxies. Using a multitracer approach, we forecast the precision of these measurements under different survey combinations. Our results indicate that current GW detectors will have limited precision, with measurement errors as high as $\displaystyle \sim 50~{{\ \rm per\ cent}}$. However, third-generation detectors like ET, when cross-correlated with Legacy Survey of Space and Time (LSST) data, can improve clustering bias measurements to within 2.5 per cent. Furthermore, we demonstrate that these cross-correlations can enable a per cent-level measurement of the magnification lensing effect on GWs. Despite this, there is a degeneracy between magnification and evolution biases, which hinders the precision of both. This degeneracy is most effectively addressed by assuming knowledge of one bias or targeting an optimal redshift range of $\displaystyle 1 \lt z \lt 2.5$. Our analysis opens new avenues for studying the distribution of BBHs and testing the nature of gravity through large-scale structure.

Zazzera, Stefano (ORCID:0000000158979221)↗
DOE OSTI2024

A detailed chemical study of the extreme velocity stars in the galaxy

ABSTRACT Two decades on, the study of hypervelocity stars is still in its infancy. These stars can provide novel constraints on the total mass of the Galaxy and its dark matter distribution. However how these stars are accelerated to such high velocities is unclear. Various proposed production mechanisms for these stars can be distinguished using chemo-dynamic tagging. The advent of Gaia and other large surveys have provided hundreds of candidate hyper velocity objects to target for ground-based high-resolution follow-up observations. We conduct high-resolution spectroscopic follow-up observations of 16 candidate late-type hyper velocity stars using the Apache Point Observatory and the McDonald Observatory. We derive atmospheric parameters and chemical abundances for these stars. We measure up to 22 elements, including the following nucleosynthetic families: $\alpha$ (Mg, Si, Ca, and Ti), light/odd-Z (Na, Al, V, Cu, and Sc), Fe-peak (Fe, Cr, Mn, Co, Ni, and Zn), and neutron capture (Sr, Y, Zr, Ba, La, Nd, and Eu). Our kinematic analysis shows one candidate is unbound, two are marginally bound, and the remainder are bound to the Galaxy. Finally, for the three unbound or marginally bound stars, we perform orbit integration to locate possible globular cluster or dwarf galaxy progenitors. We do not find any likely candidate systems for these stars and conclude that the unbound stars are likely from the the stellar halo, in agreement with the chemical results. The remaining bound stars are all chemically consistent with the stellar halo as well.

79 ASTRONOMY AND ASTROPHYSICS↗