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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hierarchical Analysis of Halo Center in Cosmology

Ever-increasing data size raises many challenges for scientific data analysis. Particularly in cosmological N-body simulation, finding the center of a dark matter halo suffers heavily from the large computational cost associated with the large number of particles (up to 20 million). In this work, we exploit the latent structure embed in a halo, and we propose a hierarchical approach to approximate the exact gravitational potential calculation for each particle in order to more efficiently find the halo center. Tests of our method on data from N-body simulations show that in many cases the hierarchical algorithm performs significantly faster than existing methods with a desirable accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Effects of Dark Matter and Baryonic Physics on the Milky Way Subhalo Population in the Presence of the Large Magellanic Cloud

Given recent developments in our understanding of the Large Magellanic Cloud's (LMC) impact on the Milky Way's (MW) dark matter subhalo population, we compare the signatures of dark matter and baryonic physics on subhalos in MW systems with realistic LMC analogs. In particular, we study the effects of self-interacting dark matter~(SIDM), warm dark matter (WDM), and the Galactic disk on the peak maximum circular velocity ($V_{\mathrm{peak}}$) function, radial distribution, and spatial distribution of MW and LMC-associated subhalos using cosmological dark matter-only zoom-in simulations of MW+LMC systems. For a fixed abundance of subhalos expected to host dwarf galaxies ($V_{\mathrm{peak}}\gtrsim 20\ \mathrm{km\ s}^{-1}$), SIDM and WDM can produce a similar mass-dependent suppression of the subhalo $V_{\mathrm{peak}}$ function, while disk disruption is mass independent. Subhalos in the inner regions of the MW are preferentially disrupted by both self-interactions and the disk, while suppression in WDM is radially independent. The relative abundance of LMC-associated subhalos is not strongly affected by disk disruption or WDM, but is significantly suppressed in SIDM due to self-interactions with the LMC at early times and with the MW during LMC infall at late times, erasing spatial anisotropy in the MW subhalo population. Furthermore, these results provide avenues to distinguish dark matter and baryonic physics by combining properties of the MW and LMC subhalo populations probed by upcoming observations of satellite galaxies and stellar streams.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quantum simulations of dark sector showers

We consider dark sector scenarios where dark matter is accompanied by a dark photon and multiple-flavor dark fermions charged under the dark gauge group. We study quantum interference effects in dark sector jets, where multiple dark photons are emitted from high-energy dark fermions. We perform fully quantum simulations of dark sector showers and compare the results against those of the classical Monte-Carlo simulations. We find important differences in probability distributions of dark photon countings between quantum and classical computations. When the number of dark-fermion flavors is large, we find significant enhancements in large numbers of dark photon emissions. Such enhancements can provide distinguishing signals for our scenarios at particle colliders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2021

Neutrinos in N -body simulations

In the next decade, cosmological surveys will have the statistical power to detect the absolute neutrino mass scale. N-body simulations of large-scale structure formation play a central role in interpreting data from such surveys. Yet these simulations are Newtonian in nature. We provide a quantitative study of the limitations to treating neutrinos, implemented as N-body particles, in N-body codes, focusing on the error introduced by neglecting special relativistic effects. Special relativistic effects are potentially important due to the large thermal velocities of neutrino particles in the simulation box. We derive a self-consistent theory of linear perturbations in Newtonian and nonrelativistic neutrinos and use this to demonstrate that N-body simulations overestimate the neutrino free-streaming scale, and cause errors in the matter power spectrum that depend on the initial redshift of the simulations. For z i ≲ 100, and neutrino masses within the currently allowed range, this error is ≲ 0.5%, though represents an up to ~10 % correction to the shape of the neutrino-induced suppression to the cold dark matter power spectrum. Here, we argue that the simulations accurately model nonlinear clustering of neutrinos so that the error is confined to linear scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Faint galaxies and small halos : probes of galaxy formation and dark matter

This thesis addresses the intersection of dark matter and galaxy formation physics by modeling the faintest galaxies in the Universe in a cosmological context. In the first part of this work, I present an empirical model for the connection between faint galaxies and small dark matter halos. Combining this model with tailored cosmological simulations of Milky Way analogs and state-of-the-art satellite galaxy observations yields new insights into dwarf galaxy formation and hierarchical structure formation. Next, I show that these same dwarf galaxy observations place stringent constraints on microphysical dark matter properties including its warmth, self-interactions, Standard Model interactions, minimum particle mass, and formation epoch. Finally, I describe a unified framework that combines dark matter constraints from dwarf galaxies and strong gravitational lensing, paving the way for next-generation surveys of small-scale cosmic structure to deliver unprecedented insights into dark matter physics.

Nadler, Ethan Oliver↗
DOE OSTI2024

Classical field approximation of ultralight dark matter: Quantum break times, corrections, and decoherence

The classical field approximation is widely used to better understand the predictions of ultralight dark matter. Here, in this work, we use the truncated Wigner approximation method to test the classical field approximation of ultralight dark matter. This method approximates a quantum state as an ensemble of independently evolving realizations drawn from its Wigner function. The method is highly parallelizable and allows the direct simulation of quantum corrections and decoherence times in systems many times larger than have been previously studied in reference to ultralight dark matter. Our study involves simulation of systems in 1, 2, and 3 spatial dimensions. We simulate three systems, the condensation of a Gaussian random field in three spatial dimensions, a stable collapsed object in three spatial dimensions, and the merging of two stable objects in two spatial dimensions. We study the quantum corrections to the classical field theory in each case. We find that quantum corrections grow exponentially during nonlinear growth with the timescale being approximately equal to the system dynamical time. In stable systems the corrections grow quadratically. We also find that the primary effect of quantum corrections is to reduce the amplitude of fluctuations on the de Broglie scale in the spatial density. Finally, we find that the timescale associated with decoherence due to gravitational coupling to baryonic matter is at least as fast as the quantum corrections due to gravitational interactions. These results are consistent with the predictions of the classical field theory being accurate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Confinement transition to gravitational waves in the one-flavor $SU(4)$ Hyper Stealth Dark Matter theory

The thermodynamics of the $SU(4)$ gauge theory with a single flavor of fundamental quarks is analyzed on the lattice with dynamical fermion simulations, which is the low-energy sector of a realistic, strongly-interacting dark matter model -- the Hyper Stealth Dark Matter. The gravitational wave spectrum from the first-order confinement transition in the early universe is further calculated, where the effect of the dark sea quarks, which decrease the interface tension in the effective potential of the Polyakov loop, is shown numerically to lower the gravitational wave amplitude.

Ayyar, V. [Boston U.]↗
DOE OSTI2024

The SRG/eROSITA All-Sky Survey: Exploring halo assembly bias with X-ray-selected superclusters

Numerical simulations indicate that the clustering of dark matter halos is not only dependent on the halo masses but has a secondary dependence on other properties, such as the assembly history of the halo. This phenomenon, known as the halo assembly bias (HAB), has been found mostly on galaxy scales; observational evidence on larger scales is scarce. In this work, we propose a novel method for exploring HAB on cluster scales using large samples of superclusters. Leveraging the largest-ever X-ray galaxy cluster and supercluster samples obtained from the first SRG/eROSITA all-sky survey, we constructed two subsamples of galaxy clusters that consist of supercluster members and isolated clusters, respectively. After correcting for the selection effects on redshift, mass, and survey depth, we computed the excess in the concentration of the intracluster gas of isolated clusters with respect to supercluster members, defined as δc gas ≡ c gas , ISO /c gas , SC –1 , to investigate the environmental effect on the concentration of clusters, a sign of HAB on cluster scales. We find that the average gas mass concentration of isolated clusters is a few percent higher than that of supercluster members, with a maximum significance of 2.8σ. The result for δc gas varies with the overdensity ratio, f, in supercluster identification, cluster mass proxies, and mass and redshift ranges but remains positive in almost all the measurements. We measure slightly larger δc gas when adopting a higher f for supercluster identification. The δc gas is also higher for low-mass and low-redshift clusters. We performed weak lensing analyses to compare the total mass concentration of the two classes and find a similar trend in total mass concentration as obtained from the gas mass concentration. Our results are consistent with the prediction of HAB on cluster scales, where halos located in denser environments are less concentrated; this trend is stronger for halos with lower masses and at lower redshifts. These phenomena can be explained by the fact that clusters in denser environments, such as superclusters, have experienced more mergers than isolated clusters in their assembling history. This work paves the way to explore HAB with X-ray superclusters and demonstrates that large samples of superclusters with X-ray and weak-lensing data can advance our understanding of the evolution of the large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The cosmology dependence of galaxy clustering and lensing from a hybrid N -body–perturbation theory model

ABSTRACT We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from N-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of N-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these non-linear spectra from redshifts z = 0 to z = 2. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h \,{\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift z ≃ 0.4 to scales of $k_{\rm max} \approx 0.6\, h\,\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy–galaxy lensing analyses with data drawn from an independent N-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming data sets from galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

WALLABY pre-pilot survey: two dark clouds in the vicinity of NGC 1395

ABSTRACT We present the Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pre-pilot observations of two ‘dark’ H i sources (with H i masses of a few times 108 $\rm {M}_\odot$ and no known stellar counterpart) that reside within 363 kpc of NGC 1395, the most massive early-type galaxy in the Eridanus group of galaxies. We investigate whether these ‘dark’ H i sources have resulted from past tidal interactions or whether they are an extreme class of low surface brightness galaxies. Our results suggest that both scenarios are possible, and not mutually exclusive. The two ‘dark’ H i sources are compact, reside in relative isolation, and are more than 159 kpc away from their nearest H i-rich galaxy neighbour. Regardless of origin, the H i sizes and masses of both ‘dark’ H i sources are consistent with the H i size–mass relationship that is found in nearby low-mass galaxies, supporting the possibility that these H i sources are an extreme class of low surface brightness galaxies. We identified three analogues of candidate primordial ‘dark’ H i galaxies within the TNG100 cosmological, hydrodynamic simulation. All three model analogues are dark matter dominated, have assembled most of their mass 12–13 Gyr ago, and have not experienced much evolution until cluster infall 1–2 Gyr ago. Our WALLABY pre-pilot science results suggest that the upcoming large-area H i surveys will have a significant impact on our understanding of low surface brightness galaxies and the physical processes that shape them.

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 OSTI2020

The structure of dissipative dark matter halos

Dissipative dark matter self-interactions can affect halo evolution and change its structure. We perform a series of controlled N-body simulations to study impacts of the dissipative interactions on halo properties. The interplay between gravitational contraction and collisional dissipation can significantly speed up the onset of gravothermal collapse, resulting in a steep inner density profile. For reasonable choices of model parameters controlling the dissipation, the collapse timescale can be a factor of 10–100 shorter than that predicted in purely elastic self-interacting dark matter. The effect is maximized when energy loss per collision is comparable to characteristic kinetic energy of dark matter particles in the halo. Our simulations provide guidance for testing the dissipative nature of dark matter with astrophysical observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A non-linear solution to the S 8 tension – II. Analysis of DES Year 3 cosmic shear

ABSTRACT Weak galaxy lensing surveys have consistently reported low values of the S8 parameter compared to the Planck lambda cold dark matter (ΛCDM) cosmology. Amon & Efstathiou used KiDS-1000 cosmic shear measurements to propose that this tension can be reconciled if the matter fluctuation spectrum is suppressed more strongly on non-linear scales than assumed in state-of-the-art hydrodynamical simulations. In this paper, we investigate cosmic shear data from the Dark Energy Survey (DES) Year 3. The non-linear suppression of the matter power spectrum required to resolve the S8 tension between DES and the Planck ΛCDM model is not as strong as inferred using KiDS data, but is still more extreme than predictions from recent numerical simulations. An alternative possibility is that non-standard dark matter contributes to the required suppression. We investigate the redshift and scale dependence of the suppression of the matter power spectrum. If our proposed explanation of the S8 tension is correct, the required suppression must extend into the mildly non-linear regime to wavenumbers $k\sim 0.2 \, h\, {\rm Mpc}^{-1}$. In addition, all measures of S8 using linear scales should agree with the Planck ΛCDM cosmology, an expectation that will be testable to high precision in the near future.

Preston, Calvin (ORCID:0000000281136916)↗
DOE OSTI2025

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey. II. Anisotropic Large-scale Coherence in Hot Gas, Galaxies, and Dark Matter

Statistics that capture the directional dependence of the baryon distribution in the cosmic web enable unique tests of cosmology and astrophysical feedback. We use constrained oriented stacking of thermal Sunyaev–Zel’dovich (tSZ) maps to measure the anisotropic distribution of hot gas 2.5–40 Mpc away from galaxy clusters embedded in massive filaments and superclusters. The cluster selection and orientation (at a scale of ∼15 Mpc) use Dark Energy Survey (DES) Year 3 data, while expanded tSZ maps from the Atacama Cosmology Telescope Data Release 6 enable a ∼3× more significant measurement of the extended gas compared to the technique’s proof-of-concept. Decomposing stacks into cosine multipoles of order m , we detect a dipole ( m = 1) and quadrupole ( m = 2) at 8σ–10σ, as well as evidence for m = 4 signal at up to 6σ, indicating sensitivity to late-time non-Gaussianity. We compare to Cardinal simulations with spherical gas models pasted onto dark matter halos. The fiducial tSZ data can discriminate between two models that deplete pressure differently in low-mass halos (mimicking astrophysical feedback), preferring higher average pressure in extended structures. However, uncertainty in the amount of cosmic infrared background contamination reduces the constraining power. Additionally, we apply the technique to DES galaxy density and weak lensing to study for the first time their oriented relationships with tSZ. In the tSZ-to-lensing relation, averaged on 7.5 Mpc (transverse) scales, we observe dependence on redshift but not shape or radial distance. Thus, on large scales, the superclustering of gas pressure, galaxies, and total matter is coherent in shape and extent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Double X/Peanut structures in barred galaxies – insights from an N -body simulation

ABSTRACT Boxy, peanut-, or X-shaped ‘bulges’ are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper, we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation, we find a stable double configuration consisting of two X/peanut structures associated with the same galactic bar – rotating with the same pattern speed – but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilizes. This feature is consistent in terms of size, strength, and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or ‘bow-tie’, morphology. It forms just after the inner peanut, and gradually extends in time (within 1–1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.

Ciambur, Bogdan C.↗
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 OSTI2022

Measurements of the ionization efficiency of protons in methane

The amount of energy released by a nuclear recoil ionizing the atoms of the active volume of detection appears “quenched” compared to an electron of the same kinetic energy. This different behavior in ionization between electrons and nuclei is described by the Ionization Quenching Factor (IQF) and it plays a crucial role in direct dark matter searches. For low kinetic energies (below 50 keV ), IQF measurements deviate significantly from common models used for theoretical predictions and simulations. We report measurements of the IQF for proton, an appropriate target for searches of Dark Matter candidates with a mass of approximately 1 GeV , with kinetic energies in between 2 keV and 13 keV in 100 mbar of methane. We used the Comimac facility in order to produce the motion of nuclei and electrons of controlled kinetic energy in the active volume, and a NEWS-G SPC to measure the deposited energy. The Comimac electrons are used as a reference to calibrate the detector with 7 energy points. A detailed study of systematic effects led to the final results well fitted by IQF (E K )=E$^α_K$ / (β+E$^α_K$) with α=0.70±0.08 and β=1.32±0.17 . In agreement with some previous works in other gas mixtures, we measured less ionization energy than predicted from SRIM simulations, the difference reaching 33% at 2 keV .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗