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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Axion Mass Prediction from Adaptive Mesh Refinement Cosmological Lattice Simulations

The quantum chromodynamics (QCD) axion arises as the pseudo-Goldstone mode of a spontaneously broken Abelian Peccei-Quinn (PQ) symmetry. If the scale of PQ symmetry breaking occurs below the inflationary reheat temperature and the domain wall number is unity, then there is a unique axion mass that gives the observed dark matter (DM) abundance. Computing this mass has been the subject of intensive numerical simulations for decades since the mass prediction informs laboratory experiments. Axion strings develop below the PQ symmetry-breaking temperature, and as the string network evolves, it emits axions that go on to become the DM. A key ingredient in the axion mass prediction is the spectral index of axion radiation emitted by the axion strings. We compute this index in this Letter using the most precise and accurate large-scale simulations to date of the axion-string network leveraging adaptive mesh refinement to achieve the precision that would, otherwise, require a static lattice with 262,144 3 lattice sites. We find a scale-invariant axion radiation spectrum to within 1% precision and find no evidence that the spectral index of radiation evolves with time. Accounting for axion production from strings prior to the QCD phase transition leads us to predict that the axion mass should be approximately 𝑚 𝑎 ∈ (45, 65) μ ⁢eV. However, we provide preliminary evidence that axions are produced in greater quantities from the string-domain-wall network collapse during the QCD phase transition, potentially increasing the mass prediction to as much as 300 μ ⁢eV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The DESI One-Percent Survey: Modelling the clustering and halo occupation of all four DESI tracers with U CHUU

We present results from a set of mock lightcones for the DESI One-Percent Survey, created from the UCHUU simulation. This 8 h −3 Gpc 3 N-body simulation comprises 2.1 trillion particles and provides high-resolution dark matter (sub)haloes in the framework of the Planck-based ΛCDM cosmology. Employing the subhalo abundance matching (SHAM) technique, we populated the UCHUU (sub)haloes with all four DESI tracers – Bright Galaxy Survey (BGS), luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars (QSOs) – to z = 2.1. Our method accounts for redshift evolution as well as the clustering dependence on luminosity and stellar mass. The two-point clustering statistics of the DESI One-Percent Survey generally agree with predictions from UCHUU across scales ranging from 0.3 h −1 Mpc to 100 h −1 Mpc for the BGS and across scales ranging from 5 h −1 Mpc to 100 h −1 Mpc for the other tracers. We observed some differences in clustering statistics that can be attributed to incompleteness of the massive end of the stellar mass function of LRGs, our use of a simplified galaxy-halo connection model for ELGs and QSOs, and cosmic variance. We find that at the high precision of UCHUU, the shape of the halo occupation distribution (HOD) of the BGS and LRG samples is smaller bias values, likely due to cosmic variance. The bias dependence on absolute magnitude, stellar mass, and redshift aligns with that of previous surveys. These results provide DESI with tools to generate high-fidelity lightcones for the remainder of the survey and enhance our understanding of the galaxy-halo connection.

cosmology↗
DOE OSTI2021

Performance study of charcoal-based radon reduction systems for ultraclean rare event detectors

The continuous emanation of radon due to trace amounts of uranium and thorium in detector materials introduces radon to the active detection volume of low-background rare event search detectors. 222 Rn produces a particularly problematic background in the physics region of interest by the "naked" beta decay of its 214 Pb daughter nucleus. While charcoal-based adsorption traps are expected to be effective for radon reduction in auxiliary circulation loops that service the warm components of current ton-scale detectors at slow flow rates (0.5-2 SLPM), radon reduction in the entire circulation loop at high flow rates $\mathcal{O}$(100s SLPM) is necessary to reach high sensitivity in future generation experiments. In this article we explore radon dynamics with a charcoal-based radon reduction system in the main circulation loop of time projection chamber detectors. Here, we find that even for perfect radon traps, circulation speeds of 2,000 SLPM are needed to reduce radon concentration in a 10 ton detector by 90%. This is faster by a factor of four than the highest circulation speeds currently achieved in dark matter detectors. We further find that the effectiveness of vacuum swing adsorption systems, which have been employed very successfully at reducing atmospheric radon levels in clean-rooms, is limited by the intrinsic radon activity of the charcoal adsorbent in ultra-low radon environments. Adsorbents with significantly lower intrinsic radon activity than in currently available activated charcoals would be necessary to build effective vacuum swing adsorption systems operated at room temperature for rare event search experiments. If such VSA systems are cooled to about 190 K, this requirement relaxes drastically.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

The electric field dependence of single electron emission in the PIXeY two-phase xenon detector

Dual phase xenon detectors are widely used in experimental searches for galactic dark matter particles. The origin of single electron backgrounds following prompt scintillation and proportional scintillation signals in these detectors is not fully understood, although there has been progress in recent years. Here, in this paper, we describe single electron backgrounds in 83m Kr calibration events and their correlation with drift and extraction fields, using the Particle Identification in Xenon at Yale (PIXeY) dual-phase xenon time projection chamber. The single electron background induced by the Fowler-Nordheim (FN) effect is measured, and its electric field dependence is quantified. The photoionization of grids and impurities by prompt scintillation and proportional scintillation also contributes to the single electron background.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Comparing galaxy formation in the L-GALAXIES semi-analytical model and the IllustrisTNG simulations

ABSTRACT We perform a comparison, object by object and statistically, between the Munich semi-analytical model, L-GALAXIES, and the IllustrisTNG hydrodynamical simulations. By running L-GALAXIES on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate, and gas content of galaxies, as well as the baryonic content of subhaloes and haloes in the two models. We find that both the stellar mass functions and the stellar masses of individual galaxies agree to better than ${\sim} 0.2\,$dex. On the other hand, specific star formation rates and gas contents can differ more substantially. At z = 0, the transition between low-mass star-forming galaxies and high-mass quenched galaxies occurs at a stellar mass scale ${\sim} 0.5\,$dex lower in IllustrisTNG than that in L-GALAXIES. IllustrisTNG also produces substantially more quenched galaxies at higher redshifts. Both models predict a halo baryon fraction close to the cosmic value for clusters, but IllustrisTNG predicts lower baryon fractions in group environments. These differences are primarily due to differences in modelling feedback from stars and supermassive black holes. The gas content and star formation rates of galaxies in and around clusters and groups differ substantially, with IllustrisTNG satellites less star forming and less gas rich. We show that environmental processes such as ram-pressure stripping are stronger and operate to larger distances and for a broader host mass range in IllustrisTNG. We suggest that the treatment of galaxy evolution in the semi-analytic model needs to be improved by prescriptions that capture local environmental effects more accurately.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Debiasing with Diffusion: Probabilistic Reconstruction of Dark Matter Fields from Galaxies with CAMELS

Abstract Galaxies are biased tracers of the underlying cosmic web, which is dominated by dark matter (DM) components that cannot be directly observed. Galaxy formation simulations can be used to study the relationship between DM density fields and galaxy distributions. However, this relationship can be sensitive to assumptions in cosmology and astrophysical processes embedded in galaxy formation models, which remain uncertain in many aspects. In this work, we develop a diffusion generative model to reconstruct DM fields from galaxies. The diffusion model is trained on the CAMELS simulation suite that contains thousands of state-of-the-art galaxy formation simulations with varying cosmological parameters and subgrid astrophysics. We demonstrate that the diffusion model can predict the unbiased posterior distribution of the underlying DM fields from the given stellar density fields while being able to marginalize over uncertainties in cosmological and astrophysical models. Interestingly, the model generalizes to simulation volumes ≈500 times larger than those it was trained on and across different galaxy formation models. The code for reproducing these results can be found athttps://github.com/victoriaono/variational-diffusion-cdm✎.

Astronomy & Astrophysics↗
NASA NTRS1992

Coherent structures in the universe and the adhesion model

An adhesion model is used to study the formation process of large-scale structures due to nonlinear gravitational growth of small initial fluctuations in the universe dominated by dark matter. The model is compared with 2D N-body simulations with initial power-law spectral indices n = -2, 0, +2, and various cutoffs. It is found that the adhesion model imitates the skeleton of the structure extremely well for the parameters of the initial spectra until the stage when the nonlinear scale reaches the correlation length R(phi) of the initial gravitational potential. The model explains the origin of large-scale coherent sructures, such as superpancakes and superfilaments, as a result of coherent motion of clumps due to large-scale inhomogeneities in the initial gravitational potential. It is found that clumps of mass identified in the N-body simulations correspond to several knots in the adhesion model, which influences the way of calculating the mass distribution function. The distribution functions of velocities and masses of clumps and areas of cells in the adhesion model satisfy self-similar scaling laws of the n = 2 model.

Kofman, Lev↗
DOE OSTI2020

Axion structure formation – II. The wrath of collapse

ABSTRACT The first paper in this series showed that quantum chromodynamic axion dark matter, as a highly correlated Bose fluid, contains extra-classical physics on cosmological scales. The source of the derived extra-classical physics is exchange–correlation interactions induced by the constraints of symmetric particle exchange and interaxion correlations from self-gravitation. The paper also showed that the impact of extra-classical physics on early structure formation is marginal, as the exchange–correlation interaction is inherently non-linear. This paper continues the study of axion structure formation into the non-linear regime, considering the case of full collapse and virialization. The N-body method is chosen to study the collapse, and its algorithms are derived for a condensed Bose fluid. Simulations of isolated gravitational collapse are performed for both Bose and cold dark matter fluids using a prototype N-body code. Unique Bose structures are found to survive even the most violent collapses. Bose post-collapse features include dynamical changes to global structures, creation of new broad sub-structures, violations of classical binding energy conditions, and new fine structures. Effective models of the novel structures are constructed and possibilities for their observation are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

The peculiar velocities of rich clusters in the hot and cold dark matter scenarios

We present the results of a study of the peculiar velocities of rich clusters of galaxies. The peculiar motion of rich clusters in various cosmological scenarios is of interest for a number of reasons. Observationally, one can measure the peculiar motion of clusters to greater distances than galaxies because cluster peculiar motions can be determined to greater accuracy. One can also test the slope of distance indicator relations using clusters to see if galaxy properties vary with environment. We have used N-body simulations to measure the amplitude and rms cluster peculiar velocity as a function of bias parameter in the hot and cold dark matter scenarios. In addition to measuring the mean and rms peculiar velocity of clusters in the two models, we determined whether the peculiar velocity vector of a given cluster is well aligned with the gravity vector due to all the particles in the simulation and the gravity vector due to the particles present only in the clusters. We have investigated the peculiar velocities of rich clusters of galaxies in the cold dark matter and hot dark matter galaxy formation scenarios. We have derived peculiar velocities and associated errors for the scenarios using four values of the bias parameter ranging from b = 1 to b = 2.5. The growth of the mean peculiar velocity with scale factor has been determined and compared to that predicted by linear theory. In addition, we have compared the orientation of force and velocity in these simulations to see if a program such as that proposed by Bertschinger and Dekel (1989) for elliptical galaxy peculiar motions can be applied to clusters. The method they describe enables one to recover the density field from large scale redshift distance samples. The method makes it possible to do this when only radial velocities are known by assuming that the velocity field is curl free. Our analysis suggests that this program if applied to clusters is only realizable for models with a low value of the bias parameter, i.e., models in which the peculiar velocities of clusters are large enough that the errors do not render the analysis impracticable.

Rhee, George F.↗
DOE OSTI2022

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, ΛCDM). Using the latest state-of-the-art cosmological ΛCDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass M * > 10 10 M ⊙ at redshift z = 0. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these ΛCDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at ≥12.52σ (TNG50-1) and ≥14.82σ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the ΛCDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an 8 5 times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the 5.58σ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in ΛCDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

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

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

Axions↗
DOE OSTI2024

How the Galaxy–Halo Connection Depends on Large-scale Environment

We investigate the connection between galaxies, dark matter halos, and their large-scale environments at z = 0 with Illustris TNG300 hydrodynamic simulation data. We predict stellar masses from subhalo properties to test two types of machine learning (ML) models: explainable boosting machines (EBMs) with simple galaxy environment features and E(3)-invariant graph neural networks (GNNs). The best-performing EBM models leverage spherically averaged overdensity features on 3 Mpc scales. Interpretations via SHapley Additive exPlanations also suggest that in the context of the TNG300 galaxy–halo connection, simple spherical overdensity on ∼3 Mpc scales is more important than cosmic web distance features measured using the DisPerSE algorithm. Meanwhile, a GNN with connectivity defined by a fixed linking length, L, outperforms the EBM models by a significant margin. As we increase the linking length scale, GNNs learn important environmental contributions up to the largest scales we probe (L = 10 Mpc). We conclude that 3 Mpc distance scales are most critical for describing the TNG galaxy–halo connection using the spherical overdensity parameterization, but that information on larger scales, which is not captured by simple environmental parameters or cosmic web features, can further augment these models. Our study highlights the benefits of using interpretable ML algorithms to explain models of astrophysical phenomena, and the power of using GNNs to flexibly learn complex relationships directly from data while imposing constraints from physical symmetries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Concentration–Mass relation of massive, dynamically relaxed galaxy clusters: agreement between observations and ΛCDM simulations

ABSTRACT The relationship linking a galaxy cluster’s total mass with the concentration of its mass profile and its redshift is a fundamental prediction of the Cold Dark Matter (CDM) paradigm of cosmic structure formation. However, confronting those predictions with observations is complicated by the fact that simulated clusters are not representative of observed samples where detailed mass profile constraints are possible. In this work, we calculate the symmetry-peakiness-alignment (SPA) morphology metrics for maps of X-ray emissivity from the three hundred project hydrodynamical simulations of galaxy clusters at four redshifts, and thereby select a sample of morphologically relaxed, simulated clusters, using observational criteria. These clusters have on average earlier formation times than the full sample, confirming that they are both morphologically and dynamically more relaxed than typical. We constrain the concentration–mass–redshift relation of both the relaxed and complete sample of simulated clusters, assuming power-law dependences on mass (κm) and 1 + z (κζ), finding κm = −0.12 ± 0.07 and κζ = −0.27 ± 0.19 for the relaxed subsample. From an equivalently selected sample of massive, relaxed clusters observed with Chandra, we find κm = −0.12 ± 0.08 and κζ = −0.48 ± 0.19, in good agreement with the simulation predictions. The simulated and observed samples also agree well on the average concentration at a pivot mass and redshift providing further validation of the ΛCDM paradigm in the properties of the largest gravitationally collapsed structures observed. This also represents the first clear detection of decreasing concentration with redshift, a longstanding prediction of simulations in data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A parametric model for self-interacting dark matter halos

Abstract We propose a parametric model for studying self-interacting dark matter (SIDM) halos. The model uses an analytical density profile, calibrated using a controlled N-body SIDM simulation that covers the entire gravothermal evolution, including core-forming and -collapsing phases. By normalizing the calibrated density profile, we obtain a universal description for SIDM halos at any evolution phase. The model allows us to infer properties of SIDM halos based on their cold dark matter (CDM) counterparts. As a basic application, we only require two characteristic parameters of an isolated CDM halo atz = 0. We then extend the model to incorporate effects induced by halo mass changes, such as major mergers or tidal stripping, making it applicable to both isolated halos and subhalos. The parametric model is tested and validated using cosmological zoom-in SIDM simulations available in the literature.

Astronomy & Astrophysics↗
DOE OSTI2024

Matter power spectra in modified gravity: a comparative study of approximations and N -body simulations

ABSTRACT Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade’s largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon, and K-mouflage. In particular, we compare the predictions from N-body simulations solving the full scalar field equation, two N-body codes with approximate time integration schemes, a parametrized modified N-body implementation, and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in 2 per cent agreement at $z\le 1$ and $k\le 1~h\,{\rm Mpc}^{-1}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.

Bose, B. (ORCID:0000000319658614)↗
DOE OSTI2021

Investigating the XENON1T low-energy electronic recoil excess using NEST

The search for dark matter, the missing mass of the Universe, is one of the most active fields of study within particle physics. The XENON1T experiment recently observed a 3.5⁢𝜎 excess potentially consistent with dark matter, or with solar axions. Here, we will use the Noble Element Simulation Technique (NEST) software to simulate the XENON1T detector, reproducing the excess. We utilize different detector efficiency and energy reconstruction models, but they primarily impact sub-keV energies and cannot explain the XENON1T excess. However, using NEST, we can reproduce their excess in multiple, unique ways, most easily via the addition of 31 ± 11 37 Ar decays. Furthermore, this results in new, modified background models, reducing the significance of the excess to ≤2.2⁢𝜎 at least using non-Profile Likelihood Ratio (PLR) methods. This is independent confirmation that the excess is a real effect, but potentially explicable by known physics. Many cross-checks of our 37 Ar hypothesis are presented.

particle dark matter↗
DOE OSTI2025

Interference with gravitational instability: Hot and fuzzy dark matter

Wave or fuzzy dark matter produced with high momenta behaves in many ways like hot particle dark matter while also possessing seemingly different phenomenology due to wave interference. We develop wave perturbation theory to show that white noise density fluctuations generated by the interference of high-momenta waves are gravitationally unstable in the usual way during matter domination above the free-streaming scale and stabilize below the free-streaming scale, much like the analogous effects for massive neutrinos in hot dark matter. We verify and illustrate these effects in the density power spectra of Newtonian Schrödinger-Poisson simulations. In the cosmological context, this would cause a gradual suppression of the initial white noise isocurvature perturbations below the free-streaming scale at matter radiation equality, unlike cold dark matter isocurvature fluctuations, and virial stability of dark matter halos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

Gas dynamic simulations of galaxy formation

Results are presented from a simulation modeling the formation of a group of galaxies in a 'standard' cold, dark matter universe with delta = 1, h sub 0 = 50 km/(s(Mpc)), baryon fraction omega sub b = 0.1 and spectrum normalization sigma sub 8 = 0.6 (bias parameter b = 1.7). Initial conditions are generated within a periodic box with comoving length 16 Mpc in a manner constrained to produce a small cluster of total mass approximately 10 exp 14 solar mass. Two sets of 643 particles are used to model the dark matter and baryon fluids. Each gas particle represents 1.08 x 10 exp -8 solar mass, implying an L* galaxy is resolved by approximately 1000 particles. The system is evolved self-consistently in three dimensions using the combined N-body/hydrodynamic scheme P3MSPH up to a final redshift z = 1. Evolving to the present is prohibited by the fact that the mean density in the simulated volume is above critical and the entire volume would be going nonlinear beyond this point, We are currently analyzing another run with somewhat poorer mass resolution which was evolved to the present.

Evrard, August E.↗