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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Galaxy clusters and cold dark matter - A low-density unbiased universe?

Large-scale simulations of a universe dominated by cold dark matter (CDM) are tested against two fundamental properties of clusters of galaxies: the cluster mass function and the cluster correlation function. We find that standard biased CDM models are inconsistent with these observations for any bias parameter b. A low-density, low-bias CDM-type model, with or without a cosmological constant, appears to be consistent with both the cluster mass function and the cluster correlations. The low-density model agrees well with the observed correlation function of the Abell, Automatic Plate Measuring Facility (APM), and Edinburgh-Durham cluster catalogs. The model is in excellent agreement with the observed dependence of the correlation strength on cluster mean separation, reproducing the measured universal dimensionless cluster correlation. The low-density model is also consistent with other large-scale structure observations, including the APM angular galaxy-correlations, and for lambda = 1-Omega with the COBE results of the microwave background radiation fluctuations.

Bahcall, Neta A.↗
DOE OSTI2022

Cuspy dark matter density profiles in massive dwarf galaxies

ABSTRACT Rotation curves of galaxies probe their total mass distributions, including dark matter. Dwarf galaxies are excellent systems to investigate the dark matter density distribution, as they tend to have larger fractions of dark matter compared to higher mass systems. The core-cusp problem describes the discrepancy found in the slope of the dark matter density profile in the centres of galaxies (β*) between observations of dwarf galaxies (shallower cores) and dark matter-only simulations (steeper cusps). We investigate β* in six nearby spiral dwarf galaxies for which high-resolution CO J = 1–0 data were obtained with ALMA (Atacama Large Millimeter/submillimeter Array). We derive rotation curves and decompose the mass profile of the dark matter using our CO rotation curves as a tracer of the total potential and 4.5 $\mu$m photometry to define the stellar mass distribution. We find 〈β*〉 = 0.6 with a standard deviation of ±0.1 among the galaxies in this sample, in agreement with previous measurements in this mass range. The galaxies studied are on the high stellar mass end of dwarf galaxies and have cuspier profiles than lower mass dwarfs, in agreement with other observations. When the same definition of the slope is used, we observe steeper slopes than predicted by the FIRE and NIHAO simulations. This may signal that these relatively massive dwarfs underwent stronger gas inflows towards their centres than predicted by these simulations, that these simulations overpredict the frequency of accretion or feedback events, or that a combination of these or other effects are at work.

Cooke, Lauren H. (ORCID:0000000284337553)↗
NASA NTRS1993

Redshift space clustering of galaxies and cold dark matter model

The distorting effect of peculiar velocities on the power speturm and correlation function of IRAS and optical galaxies is studied. The observed redshift space power spectra and correlation functions of IRAS and optical the galaxies over the entire range of scales are directly compared with the corresponding redshift space distributions using large-scale computer simulations of cold dark matter (CDM) models in order to study the distortion effect of peculiar velocities on the power spectrum and correlation function of the galaxies. It is found that the observed power spectrum of IRAS and optical galaxies is consistent with the spectrum of an Omega = 1 CDM model. The problems that such a model currently faces may be related more to the high value of Omega in the model than to the shape of the spectrum. A low-density CDM model is also investigated and found to be consistent with the data.

Bahcall, Neta A.↗
DOE OSTI2025

Aggressively Dissipative Dark Dwarfs: The Effects of Atomic Dark Matter on the Inner Densities of Isolated Dwarf Galaxies

We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of cold dark matter and a strongly dissipative subcomponent. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as atomic dark matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of ∼10 10 M ⊙ systems that vary over the aDM microphysics and the dwarf’s evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts, and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark “clumps,” which are clusters of subresolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction (∼5%) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.

cold dark matter↗
DOE OSTI2020

Optimizing Simulation Parameters for Weak Lensing Analyses Involving Non-Gaussian Observables

We performed a series of numerical experiments to quantify the sensitivity of the predictions for weak lensing statistics obtained in ray-tracing dark matter (DM)-only simulations, to two hyper-parameters that influence the accuracy as well as the computational cost of the predictions: the thickness of the lens planes used to build past light cones and the mass resolution of the underlying DM simulation. The statistics considered are the power spectrum (PS) and a series of non-Gaussian observables, including the one-point probability density function, lensing peaks, and Minkowski functionals. Counterintuitively, we find that using thin lens planes (< 60 h {sup −1} Mpc on a 240 h {sup −1} Mpc simulation box) suppresses the PS over a broad range of scales beyond what would be acceptable for a survey comparable to the Large Synoptic Survey Telescope (LSST). A mass resolution of 7.2 × 10{sup 11} h {sup −1} M {sub ⊙} per DM particle (or 256{sup 3} particles in a (240 h {sup −1} Mpc){sup 3} box) is sufficient to extract information using the PS and non-Gaussian statistics from weak lensing data at angular scales down to 1′ with LSST-like levels of shape noise.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

X-ray morphology of cluster-mass haloes in self-interacting dark matter

ABSTRACT We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq (1\!-\!3)\times 10^{15}{\, \rm M_\odot }$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\!-\!0.2R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(\sigma /m)= 0.5\!-\!1\, {\rm cm}^2\, {\rm g}^{-1}$ than CDM and SIDM with $(\sigma /m)=0.1\, {\rm cm}^2\, {\rm g}^{-1}$. Our statistical analysis indicates that the latter two models are disfavoured at the $68{{\ \rm per\ cent}}$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.

Shen, Xuejian (ORCID:000000026196823X)↗
DOE OSTI2025

Dark matter decay signals in cosmic filaments

Cosmic filaments form the backbone of the cosmic web, yet their properties and evolution remain uncertain. Using the EAGLE hydrodynamical simulations, we investigate the dark matter density profiles in filaments and their implications for dark matter decay signals. We show that GeV-scale dark matter particles decaying into electron-positron pairs can produce detectable radio synchrotron emission. By leveraging stacked filament radio data, we place stringent constraints on the dark matter decay lifetime, improving existing limits by up to two orders of magnitude for strong filamentary magnetic fields.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1990

Large N-body simulations of a universe dominated by cold dark matter

It is shown that pure gravitational evolution from an initial random density field with the standard cold dark matter (CDM) power spectrum does occasionally lead to filaments, walls, and voids remarkably like those seen in observed samples. A problem with the flat CDM model is the degree of streaming of galaxies relative to their velocity dispersion. A biased open CDM model seems to be more compatible with the observational data though it requires a higher amplitude of microwave background fluctuations.

Park, Changbom↗
DOE OSTI2023

On the statistical theory of self-gravitating collisionless dark matter flow

Dark matter, if it exists, accounts for five times as much as the ordinary baryonic matter. Compared to hydrodynamic turbulence, the flow of dark matter might possess the widest presence in our universe. This paper presents a statistical theory for the flow of dark matter that is compared with N-body simulations. By contrast to hydrodynamics of normal fluids, dark matter flow is self-gravitating, long-range, and collisionless with a scale-dependent flow behavior. The peculiar velocity field is of constant divergence nature on small scale and irrotational on large scale. The statistical measures, i.e., correlation, structure, dispersion, and spectrum functions, are modeled on both small and large scales, respectively. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. Incompressible and constant divergence flow share the same kinematic relations for even order correlations. The limiting correlation of velocity $\mathrm{ρ_{L}=1/2}$ on the smallest scale ( r = 0) is a unique feature of collisionless flow (⁠$\mathrm{ρ_{L}=1}$ for incompressible flow). On large scale, transverse velocity correlation has an exponential form $T_{2}∝e^{–r/r_2}$ with a constant comoving scale r 2 =21.3 Mpc/h that may be related to the horizon size at matter–radiation equality. All other correlation, structure, dispersion, and spectrum functions for velocity, density, and potential fields are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of ${S}_{2}^{1}∝r^{1/4}$. All other statistical measures can be obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for scale r between 1 and 7 Mpc/h. Divergence is negatively correlated for r > 30 Mpc/h that leads to a negative density correlation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Local Group timing argument and virial theorem mass estimators from cosmological simulations

ABSTRACT We identify Local Group (LG) analogues in the IllustrisTNG cosmological simulation, and use these to study two-mass estimators for the LG: One based on the timing argument (TA) and one based on the virial theorem (VT). Including updated measurements of the Milky Way-M31 tangential velocity and the cosmological constant, we show that the TA mass estimator slightly overestimates the true median LG-mass, though the ratio of the TA to the true mass is consistent at the approximate 90 per cent c.l.. These are in broad agreement with previous results using dark matter-only simulations. We show that the VT estimator better estimates the true LG-mass, though there is a larger scatter in the virial mass to true mass ratio relative to the corresponding ratio for the TA. We attribute the broader scatter in the VT estimator to several factors, including the predominantly radial orbits for LG satellite galaxies, which differs from the VT assumption of isotropic orbits. With the systematic uncertainties we derive, the updated measurements of the LG mass at 90 per cent c.l. are $4.75_{-2.41}^{+2.22} \times 10^{12}$ M⊙ from the TA and $2.0_{-1.5}^{+2.1} \times 10^{12}$ M⊙ from the VT. We consider the LMC’s effect on the TA and VT LG mass estimates, and do not find exact LMC–MW–M31 analogues in the Illustris simulations. However, in LG simulations with satellite companions as massive as the LMC, we find that the effect on the TA and VT estimators is small, though we need further studies on a larger sample of LMC–MW–M31 systems to confirm these results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Symphony: Cosmological Zoom-in Simulation Suites over Four Decades of Host Halo Mass

Abstract We present Symphony, a compilation of 262 cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from 10 11 –10 15 M ⊙ . This compilation includes three existing simulation suites at the cluster and Milky Way–mass scales, and two new suites: 39 Large Magellanic Cloud-mass (10 11 M ⊙ ) and 49 strong-lens-analog (10 13 M ⊙ ) group-mass hosts. Across the entire host halo mass range, the highest-resolution regions in these simulations are resolved with a dark matter particle mass of ≈3 × 10 −7 times the host virial mass and a Plummer-equivalent gravitational softening length of ≈9 × 10 −4 times the host virial radius, on average. We measure correlations between subhalo abundance and host concentration, formation time, and maximum subhalo mass, all of which peak at the Milky Way host halo mass scale. Subhalo abundances are ≈50% higher in clusters than in lower-mass hosts at fixed sub-to-host halo mass ratios. Subhalo radial distributions are approximately self-similar as a function of host mass and are less concentrated than hosts’ underlying dark matter distributions. We compare our results to the semianalytic model Galacticus , which predicts subhalo mass functions with a higher normalization at the low-mass end and radial distributions that are slightly more concentrated than Symphony. We use UniverseMachine to model halo and subhalo star formation histories in Symphony, and we demonstrate that these predictions resolve the formation histories of the halos that host nearly all currently observable satellite galaxies in the universe. To promote open use of Symphony, data products are publicly available at http://web.stanford.edu/group/gfc/symphony .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Galaxies and gas in a cold dark matter universe

We use a combined gravity/hydrodynamics code to simulate the formation of structure in a random 22 Mpc cube of a cold dark matter universe. Adiabatic compression and shocks heat much of the gas to temperatures of 10 exp 6 - 10 exp 7 K, but a fraction of the gas cools radiatively to about 10 exp 4 K and condenses into discrete, highly overdense lumps. We identify these lumps with galaxies. The high-mass end of their baryonic mass function fits the form of the observed galaxy luminosity function. They retain independent identities after their dark halos merge, so gravitational clustering produces groups of galaxies embedded in relatively smooth envelopes of hot gas and dark matter. The galaxy correlation function is approximately an r exp -2.1 power law from separations of 35 kpc to 7 Mpc. Galaxy fluctuations are biased relative to dark matter fluctuations by a factor b about 1.5. We find no significant 'velocity bias' between galaxies and dark matter particles. However, virial analysis of the simulation's richest group leads to an estimated Omega of about 0.3, even though the simulation adopts Omega = 1.

Katz, Neal↗
DOE OSTI2022

A Hydro-particle-mesh Code for Efficient and Rapid Simulations of the Intracluster Medium

We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm tracks only one effective particle species, the updated version separately tracks the gas and dark matter particles, as they do not exactly trace each other on small scales. For the approximate hydrodynamics solver, the pressure term in the gas equations of motion is calculated using robust physical models. In particular, we use a dark matter halo model, ICM pressure profile, and IGM temperature–density relation, all of which can be systematically varied for parameter-space studies. We show that the HYPER simulation results are in good agreement with the halo model expectations for the density, temperature, and pressure radial profiles. Simulated galaxy cluster scaling relations for Sunyaev–Zel’dovich (SZ) and X-ray observables are also in good agreement with mean predictions, with scatter comparable to that found in hydrodynamic simulations. HYPER also produces lightcone catalogs of dark matter halos and full-sky tomographic maps of the lensing convergence, SZ effect, and X-ray emission. These simulation products are useful for testing data analysis pipelines, generating training data for machine learning, understanding selection and systematic effects, and for interpreting astrophysical and cosmological constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

AbacusSummit Cleaned Halo Catalogs

This is a data release from the AbacusSummit suite of cosmological N-body simulations. It contains dark matter halo catalogs and particle samples from those halos. The halo catalogs have been augmented (cleaned) using halo merger tree information. Descriptions of AbacusSummit, the cleaned halo catalogs, and instructions for manipulating these catalogs can be found at the following URL: abacussummit.readthedocs.io/ The halo catalog cleaning is structured as a set of auxiliary files that are ingested by the reader code at runtime. But this release contains both the original catalogs and the auxiliary files, thus forming a self-contained dataset. The files in this release are aggregated by redshift. We define sets of simulations (like all c000 sims or all hugebase sims), and then store one redshift from all simulations per tarball. This is designed to support the use-case where a user wants to analyze multiple simulations at a certain redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

MultiCAM: a multivariable framework for connecting the mass accretion history of haloes with their properties

ABSTRACT Models that connect galaxy and halo properties often summarize a halo’s mass accretion history (MAH) with a single value and use this value as the basis for predictions. However, a single-value summary fails to capture the complexity of MAHs and information can be lost in the process. We present MultiCAM, a generalization of traditional abundance matching frameworks, which can simultaneously connect the full MAH of a halo with multiple halo and/or galaxy properties. As a first case study, we apply MultiCAM to the problem of connecting dark matter halo properties to their MAHs in the context of a dark matter-only simulation. While some halo properties, such as concentration, are more strongly correlated to the early-time mass growth of a halo, others, like the virial ratio, have stronger correlations with late-time mass growth. This highlights the necessity of considering the impact of the entire MAH on halo properties. For most of the halo properties we consider, we find that MultiCAM models that use the full MAH achieve higher accuracy than conditional abundance matching models that use a single epoch. We also demonstrate an extension of MultiCAM that captures the covariance between predicted halo properties. This extension provides a baseline model for applications where the covariance between predicted properties is important.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Cold dark matter. 2: Spatial and velocity statistics

We examine high-resolution gravitational N-body simulations of the omega = 1 cold dark matter (CDM) model in order to determine whether there is any normalization of the initial density fluctuation spectrum that yields acceptable results for galaxy clustering and velocities. Dense dark matter halos in the evolved mass distribution are identified with luminous galaxies; the most massive halos are also considered as sites for galaxy groups, with a range of possibilities explored for the group mass-to-light ratios. We verify the earlier conclusions of White et al. (1987) for the low-amplitude (high-bias) CDM model-the galaxy correlation function is marginally acceptable but that there are too many galaxies. We also show that the peak biasing method does not accurately reproduce the results obtained using dense halos identified in the simulations themselves. The Cosmic Background Explorer (COBE) anisotropy implies a higher normalization, resulting in problems with excessive pairwise galaxy velocity dispersion unless a strong velocity bias is present. Although we confirm the strong velocity bias of halos reported by Couchman & Carlberg (1992), we show that the galaxy motions are still too large on small scales. We find no amplitude for which the CDM model can reconcile simultaneously and galaxy correlation function, the low pairwise velocity dispersion, and the richness distribution of groups and clusters. With the normalization implied by COBE, the CDM spectrum has too much power on small scales if omega = 1.

Gelb, James M.↗
DOE OSTI2021

`SkyPy`: A package for modelling the Universe

SkyPy is an open-source Python package for simulating the astrophysical sky. It comprises a library of physical and empirical models across a range of observables and a command-line script to run end-to-end simulations. The library provides functions that sample realisations of sources and their associated properties from probability distributions. Simulation pipelines are constructed from these models using a YAML-based configuration syntax, while task scheduling and data dependencies are handled internally and the modular design allows users to interface with external software. SkyPy is developed and maintained by a diverse community of domain experts with a focus on software sustainability and interoperability. By fostering development, it provides a framework for correlated simulations of a range of cosmological probes including galaxy populations, large scale structure, the cosmic microwave background, supernovae and gravitational waves. Version 0.4 implements functions that model various properties of galaxies including luminosity functions, redshift distributions and optical photometry from spectral energy distribution templates. Future releases will provide additional modules, for example, to simulate populations of dark matter halos and model the galaxy-halo connection, making use of existing software packages from the astrophysics community where appropriate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Condensation and evaporation of boson stars

Abstract Axion-like particles, including the QCD axion, are well-motivated dark matter candidates. Numerical simulations have revealed coherent soliton configurations, also known as boson stars, in the centers of axion halos. We study evolution of axion solitons immersed into a gas of axion waves with Maxwellian velocity distribution. Combining analytical approach with controlled numerical simulations we find that heavy solitons grow by condensation of axions from the gas, while light solitons evaporate. We deduce the parametric dependence of the soliton growth/evaporation rate and show that it is proportional to the rate of the kinetic relaxation in the gas. The proportionality coefficient is controlled by the product of the soliton radius and the typical gas momentum or, equivalently, the ratio of the gas and soliton virial temperatures. We discuss the asymptotics of the rate when this parameter is large or small.

Physics↗