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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Candidate cosmic filament in the GJ526 field, mapped with the NIKA2 camera

Distinctive large-scale structures have been identified in the spatial distribution of optical galaxies up to redshift z ~ 1. In the more distant universe, the relationship between the dust-obscured population of star-forming galaxies observed at millimetre wavelengths and the network of cosmic filaments of dark matter apparent in all cosmological hydrodynamical simulations is still under study. Using the NIKA2 dual-band millimetre camera, we mapped a field of ~90 arcmin 2 in the direction of the star GJ526 simultaneously in its 1.15-mm and 2.0-mm continuum wavebands to investigate the nature of the quasi-alignment of five sources found ten years earlier with the MAMBO camera at 1.2 mm. We find that these sources are not clumps of a circumstellar debris disc around this star as initially hypothesized. Rather, they must be dust-obscured star-forming galaxies, or sub-millimetre galaxies (SMGs), in the distant background. The new NIKA2 map at 1.15 mm reveals a total of seven SMGs distributed in projection on the sky along a filament-like structure crossing the whole observed field. Furthermore, we show that the NIKA2 and supplemental Herschel photometric data are compatible with a model of the spectral energy distributions (SEDs) of these sources when a common redshift of 2.5 and typical values of the dust parameters for SMGs are adopted. Hence, we speculate that these SMGs might be located in a filament of the distant ‘cosmic web’. The length of this candidate cosmic filament crossing the whole map is at least 4 cMpc (comoving), and the separations between sources are between 0.25 cMpc and 1.25 cMpc at this redshift, in line with expectations from cosmological simulations. Nonetheless, further observations to determine the precise spectroscopic redshifts of these sources are required to definitively support this hypothesis of SMGs embedded in a cosmic filament of dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Self-Interacting Dark Matter and the Origin of Ultradiffuse Galaxies NGC1052-DF2 and -DF4

Observations of ultradiffuse galaxies NGC 1052-DF2 and -DF4 show they may contain little dark matter, challenging our understanding of galaxy formation. Using controlled N-body simulations, we explore the possibility that their properties can be reproduced through tidal stripping from the elliptical galaxy NGC 1052, in both cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. To explain the dark matter deficiency, we find that a CDM halo must have a very low concentration so that it can lose sufficient inner mass in the tidal field. In contrast, SIDM favors a higher and more reasonable concentration as core formation enhances tidal mass loss. Final stellar distributions in our SIDM benchmarks are more diffuse than the CDM one, and hence the former are in better agreement with the data. We further show that a cored CDM halo model modified by strong baryonic feedback is unlikely to reproduce the observations. Our results indicate that SIDM is more favorable for the formation of dark-matter-deficient galaxies.

79 ASTRONOMY AND 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 OSTI2024

AuriDESI: mock catalogues for the DESI Milky Way Survey

The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterizing remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way’s accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution lambda-cold dark matter magnetohydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum.

dynamics↗
DOE OSTI2024

A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-scale Structure Using Sequential Methods

Strong gravitational lenses are a singular probe of the Universe’s small-scale structure—they are sensitive to the gravitational effects of low-mass (<10 10 M ⊙ ) halos even without a luminous counterpart. Recent strong-lensing analyses of dark matter structure rely on simulation-based inference (SBI). Modern SBI methods, which leverage neural networks as density estimators, have shown promise in extracting the halo-population signal. However, it is unclear whether the constraints from these models are limited by the methodology or the data. In this study, we introduce an accelerator-optimized simulation pipeline that can generate lens images with realistic subhalo populations in milliseconds. Leveraging this simulator, we identify the main limitation of our fiducial SBI analysis: training set size. We then adopt a sequential neural posterior estimation (SNPE) approach, allowing us to refine the training distribution to align with the observed data. Using only one-fifth as many mock Hubble Space Telescope images, SNPE matches the constraints on the low-mass halo population produced by our best nonsequential model. Our experiments suggest that an over 3 order-of-magnitude increase in training set size and GPU hours would be required to achieve an equivalent result without sequential methods. While the full potential of the existing lens sample remains to be explored, the notable improvement in constraining power enabled by our sequential approach highlights that current constraints are limited primarily by methodology and not the data itself. Moreover, our results emphasize the need to treat training set generation and model optimization as interconnected stages of any cosmological analysis using SBI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dwarf Galaxies at Cosmic Noon: New JWST Constraints on Satellite Models and Subhalo Tidal Evolution

The advent of JWST has revolutionized the study of faint satellite galaxies at z ≳ 1, enabling statistical constraints on galaxy evolution and the galaxy–halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at 1 < z < 3.5 from recent JWST observations with predictions from cosmological dark-matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of ∼3.5. The JWST data and our fiducial model are consistent within 1σ–2σ across the full redshift (1 < z < 3.5) and stellar mass (M⋆ > 10 7 M ⊙ ) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST’s emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Contribution of peculiar shear motions to large-scale structure

Self-gravitating shear flow instability simulations in a cold dark matter-dominated expanding Einstein-de Sitter universe have been performed. When the shear flow speed exceeds a certain threshold, self-gravitating Kelvin-Helmoholtz instability occurs, forming density voids and excesses along the shear flow layer which serve as seeds for large-scale structure formation. A possible mechanism for generating shear peculiar motions are velocity fluctuations induced by the density perturbations of the postinflation era. In this scenario, short scales grow earlier than large scales. A model of this kind may contribute to the cellular structure of the luminous mass distribution in the universe.

Mueler, Hans-Reinhard↗
NASA NTRS1995

The least-action method, cold dark matter, and omega

Peebles has suggested an interesting technique, called the least-action method, to trace positions of galaxies back in time. This method applied on the Local Group galaxies seems to indicate that we live in an omega approximately = 0.1 universe. We have studied a cold dark matter (CDM) N-body simulation with omega = 0.2 and H = 50 km/s/Mpc and compared trajectories traced back by the least-action method with the ones given by the center of mass of the CDM halos. We show that the agreement between these sets of trajectories is at best qualitative. We also show that the line-of-sight peculiar velocities of halos are underestimated. This discrepancy is due to orphans, i.e., CDM particles which do not end up in halos. We vary the value of omega in the least-action method until the line-of-sight velocities agree with the CDM ones. The best value for this omega underestimates one of the CDM simulations by a factor of 4-5.

Dunn, A. M.↗
NASA NTRS1988

The topology of large-scale structure. II - Nonlinear evolution of Gaussian models

The evolution of non-Gaussian behavior in the large-scale universe from Gaussian initial conditions is studied. Topology measures developed in previous papers are applied to the smoothed initial, final, and biased matter distributions of cold dark matter, white noise, and massive neutrino simulations. When the smoothing length is approximately twice the mass correlation length or larger, the evolved models look like the initial conditions, suggesting that random phase hypotheses in cosmology can be tested with adequate data sets. When a smaller smoothing length is used, nonlinear effects are recovered, so nonlinear effects on topology can be detected in redshift surveys after smoothing at the mean intergalaxy separation. Hot dark matter models develop manifestly non-Gaussian behavior attributable to phase correlations, with a topology reminiscent of bubble or sheet distributions. Cold dark matter models remain Gaussian, and biasing does not disguise this.

Melott, Adrian L.↗
DOE OSTI2020

Microwave Cavity Simulation Using Ansys HFSS

The design of microwave cavity detectors for axion dark matter research is often accomplished using advanced full-wave electromagnetic simulation software tools. These tools provide a cost-effective approach to evaluate a wide variety of cavity configurations, frequency tuning mechanisms, and conductive or dielectric materials and coatings. One simulation software package used for this application is Ansys High Frequency Structure Simulator (HFSS), which is based upon the well-established finite element method. HFSS includes numerous features useful for microwave cavity design such as parametric geometry modeling, adaptive meshing algorithm, curvilinear mesh elements, driven modal and eigenmode matrix solvers, and optimization algorithms. This paper describes the use of the HFSS software to simulate microwave cavities for axion haloscope detectors, with an example tutorial for a cylindrical cavity. Excellent agreement between the simulated and analytical results is shown for the resonant frequency, quality factor, and form factor.

Ansys HFSS, Cavity simulation, finite element mode↗
NASA NTRS2012

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗
DOE OSTI2024

The AGORA High-resolution Galaxy Simulations Comparison Project. V. Satellite Galaxy Populations in a Cosmological Zoom-in Simulation of a Milky Way–Mass Halo

We analyze and compare the satellite halo populations at z ~ 2 in the high-resolution cosmological zoom-in simulations of a 10 12 M ⊙ target halo (z = 0 mass) carried out on eight widely used astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Arepo-t, and Gizmo) for the AGORA High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near z = 2 for each code (hereafter "z ~ 2") at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the CosmoRun simulations and the DMO simulations, we discover that each CosmoRun halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at z ~ 2, the number of satellite galaxies is significantly fewer than that of dark matter halos in all participating AGORA simulations and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called "missing satellite problem" is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in AGORA and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution (≲100 proper pc at z = 2). We also compare other properties such as the stellar mass–halo mass relation and the mass–metallicity relation. Our work highlights the value of comparison studies such as AGORA, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2008

Resolving the Formation of Protogalaxies: Feedback from the First Stars - 3

The first stars form in dark matter halos of masses 106 M as suggested by an increasing number of numerical simulations. Radiation feedback from these stars expels most of the gas from the shallow potential well of their surrounding dark matter halos.We use cosmological adaptive mesh refinement simulations that include self-consistent Population III star formation and feedback to examine the properties of assembling early dwarf galaxies. Accurate radiative transport is modeled with adaptive ray tracing. We include supernova explosions and follow the metal enrichment of the intergalactic medium. The calculations focus on the formation of several dwarf galaxies and their progenitors. In these halos, baryon fractions in 10(exp 8) Stelar Mass halos decrease by a factor of 2 with stellar feedback and by a factor of 3 with supernova explosions.We find that radiation feedback and supernova explosions increase gaseous spin parameters up to a factor of 4 and vary with time. Stellar feedback, supernova explosions, and H2 cooling create a complex, multiphase interstellar medium whose densities and temperatures can span up to 6 orders of magnitude at a given radius. The pair-instability supernovae of Population III stars alone enrich the halos with virial temperatures of 10(exp 4) K to approximately 10(exp -3) of solar metallicity.We find that 40% of the heavy elements resides in the intergalactic medium (IGM) at the end of our calculations. The highest metallicity gas exists in supernova remnants and very dilute regions of the IGM.

Wise, John H.↗
DOE OSTI2022

Inferring dark matter substructure with astrometric lensing beyond the power spectrum

Abstract Astrometry—the precise measurement of positions and motions of celestial objects—has emerged as a promising avenue for characterizing the dark matter population in our Galaxy. By leveraging recent advances in simulation-based inference and neural network architectures, we introduce a novel method to search for global dark matter-induced gravitational lensing signatures in astrometric datasets. Our method based on neural likelihood-ratio estimation shows significantly enhanced sensitivity to a cold dark matter population and more favorable scaling with measurement noise compared to existing approaches based on two-point correlation statistics. We demonstrate the real-world viability of our method by showing it to be robust to non-trivial modeled as well as unmodeled noise features expected in astrometric measurements. This establishes machine learning as a powerful tool for characterizing dark matter using astrometric data.

convolutional neural networks (1938)↗
DOE OSTI2024

Simulating Superconducting Devices in Novel Materials Using the G4CMP Simulation Toolkit

Understanding the kinematics of phonons and charge propagation in superconducting devices is crucial for conducting low-threshold dark matter searches and minimizing correlated errors in superconducting qubits. For nearly a decade, the Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, has been limited to simulating charge and phonon transport in silicon and germanium materials. In this work, we have expanded the capabilities of G4CMP to include novel substrate materials such as sapphire (Al2O3), Gallium Arsenide (GaAs), Lithium Fluoride (LiF), Calcium Tungstate (CaWO4), and Calcium Fluoride (CaF2). We demonstrate the use of this toolkit in generating phonon transport properties for these materials and compare the results with experimentally determined values where available. Additionally, we explore the performance of superconducting devices utilizing these materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Deep probabilistic direction prediction in 3D with applications to directional dark matter detectors

Abstract We present the first method to probabilistically predict 3D direction in a deep neural network model. The probabilistic predictions are modeled as a heteroscedastic von Mises-Fisher distribution on the sphere S 2 , giving a simple way to quantify aleatoric uncertainty. This approach generalizes the cosine distance loss which is a special case of our loss function when the uncertainty is assumed to be uniform across samples. We develop approximations required to make the likelihood function and gradient calculations stable. The method is applied to the task of predicting the 3D directions of electrons, the most complex signal in a class of experimental particle physics detectors designed to demonstrate the particle nature of dark matter and study solar neutrinos. Using simulated Monte Carlo data, the initial direction of recoiling electrons is inferred from their tortuous trajectories, as captured by the 3D detectors. For 40 keV electrons in a 70% He 30% CO 2 gas mixture at STP, the new approach achieves a mean cosine distance of 0.104 (26 ∘ ) compared to 0.556 (64 ∘ ) achieved by a non-machine learning algorithm. We show that the model is well-calibrated and accuracy can be increased further by removing samples with high predicted uncertainty. This advancement in probabilistic 3D directional learning could increase the sensitivity of directional dark matter detectors.

Computer Science↗
DOE OSTI2025

Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. II. cosmological results

Here, we present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years of the Dark Energy Survey. We implement the following: (1) second and third moments; (2) wavelet phase harmonics; (3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven 𝑤 Cold Dark Matter (𝑤⁢ CDM) cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the figure of merit (𝑆 8 , Ω m ) is improved by 70% (Λ ⁢CDM) and 90% (𝑤 ⁢CDM). When all the summary statistics are combined, we achieve a 2% constraint on the amplitude of fluctuations parameter 𝑆 8 ≡ 𝜎 8 ⁢(Ω m /0.3) 0.5 , obtaining 𝑆 8 = 0.794 ±0.017 (Λ⁢ CDM) and 𝑆 8 = 0.817 ±0.021 (𝑤 ⁢CDM), and a ∼10% constraint on Ω m , obtaining Ω m =0.259 ±0.025 (Λ ⁢CDM) and Ω m = 0.273 ±0.029 (𝑤⁢ CDM). In the context of the 𝑤⁢ CDM scenario, these statistics also strengthen the constraints on the parameter 𝑤, obtaining 𝑤 <−0.72. The constraints from different statistics are shown to be internally consistent (with a 𝑝-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the first three years of the Dark Energy Survey data, finding good consistency; we also compare with results from external datasets, such as planck constraints from the cosmic microwave background, finding statistical agreement, with discrepancies no greater than <2.2⁢𝜎.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The persistence of large scale structures. Part I. Primordial non-Gaussianity

Abstract We develop an analysis pipeline for characterizing the topology of large scale structure and extracting cosmological constraints based on persistent homology . Persistent homology is a technique from topological data analysis that quantifies the multiscale topology of a data set, in our context unifying the contributions of clusters, filament loops, and cosmic voids to cosmological constraints. We describe how this method captures the imprint of primordial local non-Gaussianity on the late-time distribution of dark matter halos, using a set of N-body simulations as a proxy for real data analysis. For our best single statistic, running the pipeline on several cubic volumes of size 40 (Gpc/h) 3 , we detect f NL loc =10 at 97.5% confidence on ~ 85% of the volumes. Additionally wetest our ability to resolve degeneracies betweenthe topological signature of f NL loc and variation of σ 8 and argue that correctly identifying nonzero f NL loc in this case is possible via an optimal template method. Our method relies on information living at $\mathcal{O}$(10) Mpc/h, a complementary scale with respect to commonly used methods such as the scale-dependent bias in the halo/galaxy power spectrum. Therefore, while still requiring a large volume, our method does not require sampling long-wavelength modes to constrain primordial non-Gaussianity. Moreover, our statistics are interpretable: we are able to reproduce previous results in certain limits and we make new predictions for unexplored observables, such as filament loops formed by dark matter halos in a simulation box.

Astronomy & Astrophysics↗