Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dark matter simulations”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Strong Lensing Cosmology with Population-level Calibrated Neural Ratio Estimation

Strong gravitational lensing contains key information about cosmic acceleration. Modern and next-generation galaxy imaging surveys are expected to provide high-quality data on $\mathcal{O}(10^5)$ galaxy-galaxy lensing systems. The plethora and complexity of the data are likely to present computational challenges for parameter inference methods for fitting high-dimensional likelihoods, which are often analytically intractable. Neural Ratio Estimation (NRE) efficiently computes individual likelihood ratios that can be combined into population-level posteriors. We use simulations to study the capacity of NRE to jointly predict the dark energy equation-of-state parameter $w$ and the total matter density $Ω_{m}$ from lensing images and companion spectroscopic information. We also introduce a post hoc posterior coverage calibration procedure that mitigates the model overconfidence that is typically found in neural density estimation applications. Our experiments show that the errors on both parameters decrease with increasing inference population sizes. In particular, for 100 lenses in a standard $Λ$CDM Universe, our calibrated NRE model achieves median fractional uncertainty of $22.8\%$ in $w$ and $2.9\%$ in $Ω_{m}$. This proof of concept demonstrates a potentially scalable approach for efficient cosmological parameter inference with large populations of galaxy-scale lenses observed in future surveys.

Jarugula, Sreevani [Fermilab] (ORCID:0000000253867↗
DOE OSTI2024

Advances in ArborX to support exascale applications

ArborX is a performance portable geometric search library developed as part of the Exascale Computing Project (ECP). In this paper, we explore a collaboration between ArborX and a cosmological simulation code HACC. Large cosmological simulations on exascale platforms encounter a bottleneck due to the in-situ analysis requirements of halo finding, a problem of identifying dense clusters of dark matter (halos). This problem is solved by using a density-based DBSCAN clustering algorithm. With each MPI rank handling hundreds of millions of particles, it is imperative for the DBSCAN implementation to be efficient. In addition, the requirement to support exascale supercomputers from different vendors necessitates performance portability of the algorithm. We describe how this challenge problem guided ArborX development, and enhanced the performance and the scope of the library. We explore the improvements in the basic algorithms for the underlying search index to improve the performance, and describe several implementations of DBSCAN in ArborX. Further, we report the history of the changes in ArborX and their effect on the time to solve a representative benchmark problem, as well as demonstrate the real world impact on production end-to-end cosmology simulations.

97 MATHEMATICS AND COMPUTING↗
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

The GD-1 Stellar Stream Perturber as a Core-collapsed Self-interacting Dark Matter Halo

The GD-1 stellar stream exhibits spur and gap structures that may result from a close encounter with a dense substructure. When interpreted as a dark matter subhalo, the perturber is denser than predicted in the standard cold dark matter (CDM) model. In self-interacting dark matter (SIDM), however, a halo could evolve into a phase of gravothermal collapse, resulting in a higher central density than its CDM counterpart. We conduct high-resolution controlled N-body simulations to show that a collapsed SIDM halo could account for the GD-1 perturber's high density. We model a progenitor halo with a mass of 3 × 10 8 M ⊙ , motivated by a cosmological simulation of a Milky Way analog, and evolve it in the Milky Way's tidal field. For a cross section per mass of σ/m ≈ 30–100 cm 2 g −1 at ${V}_{{\rm{\max }}}\unicode{x0007E}10\,{\rm{km}}\,{{\rm{s}}}^{-1}$, the enclosed mass of the SIDM halo within the inner 10 pc can be increased by more than 1 order of magnitude compared to its CDM counterpart, leading to a good agreement with the properties of the GD-1 perturber. Our findings indicate that stellar streams provide a novel probe into the self-interacting nature of dark matter.

dark matter↗
DOE OSTI2024

CTA and SWGO can discover Higgsino dark matter annihilation

Thermal Higgsino dark matter (DM), with a mass near 1.1 TeV, is one of the most well-motivated and untested DM candidates. Leveraging recent hydrodynamic cosmological simulations that give DM density profiles in Milky Way analog galaxies we show that the linelike gamma-ray signal predicted from Higgsino annihilation in the Galactic Center could be detected at high significance with the upcoming Cherenkov Telescope Array (CTA) and Southern Wide-field Gamma-ray Observatory (SWGO) for all but the most pessimistic DM profiles. We perform the most sensitive search to-date for the linelike signal using 15 years of data from the Fermi Large Area Telescope, coming within an order one factor of the necessary sensitivity to detect the Higgsino for some Milky Way analog DM density profiles. We show that H.E.S.S. has subleading sensitivity relative to Fermi for the Higgsino at present. In contrast, we analyze H.E.S.S. inner Galaxy data for the thermal wino model with a mass near 2.8 TeV; we find no evidence for a DM signal and exclude the wino by over a factor of two in cross section for all DM profiles considered. In the process, we identify and attempt to correct what appears to be an inconsistency in previous H.E.S.S. inner Galaxy analyses for DM annihilation related to the analysis effective area, which may weaken the DM cross-section sensitivity claimed in those works by around an order of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Halo formation from Yukawa forces in the very early Universe

If long-range attractive forces exist and are stronger than gravity then cosmic halo formation can begin in the radiation-dominated era. We study a simple realization of this effect in a system where dark matter fermions have Yukawa interactions mediated by scalar particles, analogous to the Higgs boson in the standard model. We develop a self-consistent description of the system including exact background dynamics of the scalar field, and precise modeling of the fermion density fluctuations. For the latter, we provide accurate approximations for the linear growth as well as quantitative modeling of the nonlinear evolution using N-body simulations. We find that halo formation occurs exponentially fast and on scales substantially larger than simple estimates predict. Here, the final fate of these halos remains uncertain, but could be annihilation, dark stars, primordial black holes, or even the existence of galaxy-sized halos at matter-radiation equality. More generally, our results demonstrate the importance of mapping scalar-mediated interactions onto structure formation outcomes and constraints for beyond the standard model theories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A fully Lagrangian, non-parametric bias model for dark matter halos

Here, we present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at z=0.5. Unlike the perturbative halo bias model that has been widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity δ 1 in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of δ 1 and other quantities. Especially for massive halos with 6×10 12 h -1 M ⊙ , the halo-to-mass ratio starts soaring up at δ 1 >0, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider (M>3×10 11 h -1 M ⊙ ) a non-parametric halo-to-mass ratio as a function of δ 1 and its local derivative ∇^2δ 1 can recover the halo power spectra to sub-percent level accuracy for wavenumbers k=0.01-0.1 h Mpc -1 given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At k<0.01 h Mpc -1 and for massive halos with M>6x10 12 h -1 M ⊙ , our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Milky Way Satellite Census -- II. Galaxy-Halo Connection Constraints Including the Impact of the Large Magellanic Cloud

The population of Milky Way (MW) satellites contains the faintest known galaxies, and thus provides essential insight into galaxy formation and dark matter microphysics. Here, we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high-Galactic-latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6.5\pm 1.5$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from $\textit{Gaia}$ proper motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $2.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The phase-space structure of dark matter haloes

ABSTRACT The phase space structure of dark matter haloes can be used to measure the mass of the halo, infer mass accretion rates, and probe the effects of modified gravity. Previous studies showed that the splashback radius can be measured in position space using a sharp drop in the density profile. Using N-body simulations, we model the distribution of the kinematically distinct infalling and orbiting populations of subhaloes and haloes. We show that the two are mixed spatially all the way to redge, which extends past the splashback radius defined by the drop in the spherically averaged density profile. This edge radius can be interpreted as a radius that contains a fixed fraction of the apocentres of dark matter particles. Our results highlight the possibility of measuring the outer boundary of a dark matter halo using its phase space structure and provide a firm theoretical foundation to the satellite galaxy model adopted in the companion paper, where we analysed the phase space distribution of Sloan Digital Sky Survey redMaPPer clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Nearest neighbour distributions: New statistical measures for cosmological clustering

ABSTRACT The use of summary statistics beyond the two-point correlation function to analyse the non-Gaussian clustering on small scales, and thereby, increasing the sensitivity to the underlying cosmological parameters, is an active field of research in cosmology. In this paper, we explore a set of new summary statistics – the k-Nearest Neighbour Cumulative Distribution Functions (kNN-CDF). This is the empirical cumulative distribution function of distances from a set of volume-filling, Poisson distributed random points to the k-nearest data points, and is sensitive to all connected N-point correlations in the data. The kNN-CDF can be used to measure counts in cell, void probability distributions, and higher N-point correlation functions, all using the same formalism exploiting fast searches with spatial tree data structures. We demonstrate how it can be computed efficiently from various data sets – both discrete points, and the generalization for continuous fields. We use data from a large suite of N-body simulations to explore the sensitivity of this new statistic to various cosmological parameters, compared to the two-point correlation function, while using the same range of scales. We demonstrate that the use of kNN-CDF improves the constraints on the cosmological parameters by more than a factor of 2 when applied to the clustering of dark matter in the range of scales between 10 and $40\, h^{-1}\, {\rm Mpc}$. We also show that relative improvement is even greater when applied on the same scales to the clustering of haloes in the simulations at a fixed number density, both in real space, as well as in redshift space. Since the kNN-CDF are sensitive to all higher order connected correlation functions in the data, the gains over traditional two-point analyses are expected to grow as progressively smaller scales are included in the analysis of cosmological data, provided the higher order correlation functions are sensitive to cosmology on the scales of interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Synthetic light-cone catalogues of modern redshift and weak lensing surveys with abacussummit

ABSTRACT The joint analysis of different cosmological probes, such as galaxy clustering and weak lensing, can potentially yield invaluable insights into the nature of the primordial Universe, dark energy, and dark matter. However, the development of high-fidelity theoretical models is a necessary stepping stone. Here, we present public high-resolution weak lensing maps on the light-cone, generated using the N-body simulation suite abacussummit, and accompanying weak lensing mock catalogues, tuned to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument. Available in this release are maps of the cosmic shear, deflection angle, and convergence fields at source redshifts ranging from z = 0.15 to 2.45 as well as cosmic microwave background convergence maps for each of the 25 base-resolution simulations ($L_{\rm box} = 2000\, h^{-1}\, {\rm Mpc}$ and Npart = 69123) as well as for the two huge simulations ($L_{\rm box} = 7500\, h^{-1}\, {\rm Mpc}$ and Npart = 86403) at the fiducial abacussummit cosmology. The pixel resolution of each map is 0.21 arcmin, corresponding to a healpix Nside of 16 384. The sky coverage of the base simulations is an octant until z ≈ 0.8 (decreasing to about 1800 deg2 at z ≈ 2.4), whereas the huge simulations offer full-sky coverage until z ≈ 2.2. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper Suprime-Cam data sets. The mock catalogues are validated against theoretical predictions for various clustering and lensing statistics, such as correlation multipoles, galaxy–shear, and shear–shear, showing excellent agreement. All products can be downloaded via a Globus endpoint (see Data Availability section).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Projected WIMP sensitivity of the LUX-ZEPLIN dark matter experiment

LUX-ZEPLIN (LZ) is a next generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with Weakly Interacting Massive Particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above $1.6 \times 10^{-48}$ cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$ discovery potential is projected reaching cross sections below the existing and projected exclusion limits of similar experiments that are currently operating. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of $2.7 \times 10^{-43}$ cm$^{2}$ ($8.1 \times 10^{-42}$ cm$^{2}$) for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Energy Distribution of the Galactic Center Excess’s Sources

The Galactic Center Excess (GCE) may yet herald the discovery of annihilating dark matter. Weighing against that conclusion are analyses showing evidence for dim point sources within the spatial structure of the emission. Because of technical limitations these analyses are purely spatial with all spectral information that could disentangle the excess from astrophysical backgrounds discarded. Here, we demonstrate that a neural network simulation-based inference approach can jointly analyze the spatial and spectra data. The addition is profound: energy information drives the putative point sources to be significantly dimmer, indicating either the GCE is truly diffuse in nature or made of an exceptionally large number of sources. Quantitatively, for our best fit background model, the excess is essentially consistent with Poisson emission as predicted by dark matter. If due to point sources, our median prediction is O(10^{5}) sources, or more than 35 000 at 90% confidence-both orders of magnitude larger than the hundreds preferred by earlier point-source analyses of the GCE, although variations allowed by background systematics could reduce the required number of sources by roughly an order of magnitude.

List, Florian↗
DOE OSTI2022

Probability of Forming Gaps in the GD-1 Stream by Close Encounters of Globular Clusters

One of the most intriguing properties of the GD-1 stellar stream is the existence of three gaps. If these gaps were formed by close encounters with dark matter subhalos, the GD-1 stream opens an exciting window through which we can see the size, mass, and velocity distributions of the dark matter subhalos in the Milky Way. However, in order to use the GD-1 stream as a probe of the dark matter substructure, we need to disprove that these gaps are not due to the perturbations from baryonic components of the Milky Way. Here we ran a large number of test-particle simulations to investigate the probability that each of the known globular clusters (GCs) can form a GD-1-like gap, by using the kinematical data of the GD-1 stream and GCs from Gaia early data release 3 and by fully taking the observational uncertainty into account. We found that the probability that all of the three gaps were formed by GCs is as low as 1.7 × 10 -5 , and the expected number of gaps formed by GCs is only 0.057 in our fiducial model. Our result highly disfavors a scenario in which GCs form the gaps. Given that other baryonic perturbers (e.g., giant molecular clouds) are even less likely to form a gap in the retrograde-moving GD-1 stream, we conclude that at least one of the gaps in the GD-1 stream was formed by dark matter subhalos if the gaps were formed by flyby perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Beyond mass: detecting secondary halo properties with galaxy-galaxy lensing

ABSTRACT Secondary halo properties beyond mass, such as the mass accretion rate (MAR), concentration, and the half mass scale, are essential in understanding the formation of large-scale structure and dark matter haloes. In this paper, we study the impact of secondary halo properties on the galaxy-galaxy lensing observable, ΔΣ. We build an emulator trained on N-body simulations to model ΔΣ and quantify the impact of different secondary parameters on the ΔΣ profile. We focus on the impact of MAR on ΔΣ. We show that a 3σ detection of variations in MAR at fixed halo mass could be achieved with the Hyper Suprime Cam survey assuming no baryonic effects and a proxy for MAR with scatter <1.5. We show that the full radial profile of ΔΣ depends on secondary properties at fixed halo mass. Consequently, an emulator that can perform full shape fitting yields better than two times improvement upon the constraints on MAR than only using the outer part of the halo. Finally, we highlight that miscentring and MAR impact the radial profile of ΔΣ in a similar fashion, implying that miscentring and MAR need to be modelled jointly for unbiased estimates of both effects. We show that present-day lensing data sets have the statistical capability to place constraints on halo MAR within our assumptions. Our analysis opens up new possibilities for observationally measuring the assembly history of the dark matter haloes that host galaxies and clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

S&TR September 2025: Computing Grand Challenge Turns 20

Livermore’s Computing Grand Challenge Program enters its 20th year with more unclassified high-performance computing (HPC) power than ever before. This unique, peer-reviewed competition awards HPC allocations on top supercomputers to multidisciplinary teams with high-impact projects. The Grand Challenge encourages researchers to innovate, pushes scientific discovery to new heights, improves the Laboratory’s HPC capabilities, and extends HPC accessibility to collaborators. Awardees must adapt to successive generations of HPC hardware and learn to run simulations at scale. The feature article spotlights three Grand Challenge teams whose research broke new ground in key scientific pursuits—the essence of dark matter, explosion-generated seismic waves, and protein interactions linked to cancer—while underscoring the importance of academic partnerships and considering the program’s future.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Antihelium-3 sensitivity for the GRAMS experiment

The Gamma-Ray and AntiMatter Survey (GRAMS) is a next-generation balloon/satellite mission utilizing a Liquid Argon Time Projection Chamber (LArTPC) detector to measure both MeV gamma rays and antinuclei produced by dark matter annihilation or decay. The GRAMS can identify antihelium-3 events based on the measurements of X-rays and charged pions from the decay of the exotic atoms, Time of Flight (TOF), energy deposition, and stopping range. This paper shows the antihelium-3 sensitivity estimation using a GEANT4 Monte Carlo simulation. For the proposed long-duration balloon (LDB) flight program (35 days x 3 flights) and future satellite mission (2-year observation/10-year observation), the sensitivities become 1.47 × 10 −7 [m 2 s sr GeV/n] −1 and 1.55 × 10 −9 [m 2 s sr GeV/n] −1 /3.10 x 10 -10 [m 2 s sr GeV/n] −1 , respectively. The results indicate that GRAMS can extensively investigate various dark matter models through the antihelium-3 measurements.

Antihelium-3↗
DOE OSTI2021

The Last Journey. II. SMACC—Subhalo Mass-loss Analysis Using Core Catalogs

Here we introduce Subhalo Mass-loss Analysis using Core Catalogs (SMACC). SMACC adds a mass model to substructure merger trees based on halo "core tracking." Our approach avoids the need for running expensive subhalo finding algorithms and instead uses subhalo mass-loss modeling to assign masses to halo cores. We present details of the SMACC methodology and demonstrate its excellent performance in describing halo substructure and its evolution. Validation of the approach is carried out using cosmological simulations at significantly different resolutions. We apply SMACC to the 1.24-trillion-particle Last Journey simulation and construct core catalogs with the additional mass information. These catalogs can be readily used as input to semianalytic models or subhalo abundance matching approaches to determine approximate galaxy distributions, as well as for in-depth studies of small-scale structure evolution.

79 ASTRONOMY AND ASTROPHYSICS↗