Search NASA⌕ Search

SEARCH · Search NASA

Results for “cosmological 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dark Energy Survey Year 3 results: Cosmology from combined galaxy clustering and lensing validation on cosmological simulations

In this work, we present a validation of the Dark Energy Survey Year 3 (DES Y3) 3×2-point analysis choices by testing them on Buzzard2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak-lensing analyses. We show that the buzzard2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak-lensing statistics. We then show that our model for the 3×2-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the nonlinear matter power spectrum, nonlinear galaxy bias, and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy–galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fast Approximate Methods for Modified Gravity Cosmological Simulations

The accelerated expansion of the Universe poses a major theoretical puzzle. Although the assumption of a non-zero cosmological constant provides a minimal extension of general relativity that is consistent with observational data, many theories of modified gravity have been suggested as possible alternatives. Predictions of structure formation for these models in the fully nonlinear regime are very expensive and it is difficult, if not impossible, to explore such a huge space of models and parameters using high-resolution N-body simulations. Even in the mildly nonlinear regime, perturbative methods can become extremely complex. We explore whether simplified dynamical approximations, applicable for a certain set of cosmological probes, can be used to investigate models of modified gravity with acceptable accuracy in the latter instance. For the case of chameleon gravity, we find that these methods can indeed be used to explore the region around the baryon acoustic oscillation scale, k ~ 0:1 hMpc-1 but not much further.

cosmology↗
DOE OSTI2025

The HalfDome multi-survey cosmological simulations: N-body simulations

Upcoming cosmological surveys have the potential to reach groundbreaking discoveries on multiple fronts, including the neutrino mass, dark energy, and inflation. Most of the key science goals require the joint analysis of datasets from multiple surveys to break parameter degeneracies and calibrate systematics. To realize such analyses, a large set of mock simulations that realistically model correlated observables is required. In this paper we present the N-body component of the HalfDome cosmological simulations, designed for the joint analysis of Stage-IV cosmological surveys, such as Rubin LSST, Euclid, SPHEREx, Roman, DESI, PFS, Simons Observatory, CMB-S4, and LiteBIRD. Our 300TB initial data release includes full-sky lightcones and halo catalogs between z = 0–4 for 11 fixed cosmology realizations, as well as an additional run with local primordial non-Gaussianity (f NL = 20). The simulations evolve 61443 particles in a 3.75 h -1 Gpc box, reaching a minimum halo mass of ∼6 × 1012 h -1 M ⊙ and maximum scale of k ∼ h Mpc-1. Our data is publicly available: instructions to access the data and plans for future data releases can be found at https://halfdomesims.github.io.

Bayer, Adrian E↗
DOE OSTI2024

guppy i : a code for reducing the storage requirements of cosmological simulations

ABSTRACT As cosmological simulations have grown in size, the permanent storage requirements of their particle data have also grown. Even modest simulations present a major logistical challenge for the groups which run these boxes and researchers without access to high performance computing facilities often need to restrict their analysis to lower quality data. In this paper, we present guppy, a compression algorithm and code base tailored to reduce the sizes of dark matter-only cosmological simulations by approximately an order of magnitude. guppy is a ‘lossy’ algorithm, meaning that it injects a small amount of controlled and uncorrelated noise into particle properties. We perform extensive tests on the impact that this noise has on the internal structure of dark matter haloes, and identify conservative accuracy limits which ensure that compression has no practical impact on single-snapshot halo properties, profiles, and abundances. We also release functional prototype libraries in C, Python, and Go for reading and creating guppy data.

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 OSTI2025

Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: new models and constraints from the eBOSS measurement

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-α forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with 40963 hydrodynamical cells in a 120 Mpc (∼ 81 Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-α forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters A Lyα and n Lyα . While the best-fit Planck ΛCDM model has A Lyα = 8.79 and n Lyα = -2.363, from DR14 eBOSS data we find that A Lyα < 7.6 (95% CI) and n Lyα = -2.369 ± 0.008. The low value of A Lyα , in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-α forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find A Lyα = 9.8 ± 1.1 in accordance with Planck.

Walther, Michael↗
DOE OSTI2025

Mechanical cosmology: Simulating scalar fluctuations in expanding universes using synthetic mechanical lattices

Inspired by recent advances in observational astrophysics and continued explorations in the field of analog gravity, we discuss the prospect of simulating models of cosmology within the context of synthetic mechanical lattice experiments. We focus on the physics of expanding universe scenarios described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric. Specifically, quantizing scalar fluctuations in a background FLRW spacetime leads to a quadratic bosonic Hamiltonian with temporally varying pair production terms. Here we present a mapping that provides a one-to-one correspondence between these classes of cosmology models and feedback-coupled mechanical oscillators. As proof of principle, we then perform experiments on a synthetic mechanical lattice composed of such oscillators. We simulate two different FLRW expansion scenarios with universes dominated by vacuum energy and matter and discuss our experimental results. Published by the American Physical Society 2025

Rhyno, Brendan (ORCID:0000000235266998)↗
DOE OSTI2024

Comparative Analysis of Cosmological Simulations via Angular Power Spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use angular power spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE) [1]. After ensuring that the simulations can produce APS that match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Comparative analysis of cosmological simulations via angular power spectra

Cosmology is the study of the universe's origin, large-scale structure, and dynamics. One such area of study is understanding dark energy and dark matter. Dark matter has not been observed to interact with luminous matter other than through gravitational attractions, forming large scale structures observed throughout the universe. In this area, the primary model is the Lambda Cold Dark Matter (ΛCDM) model. This standard model consists of cosmological dark energy (Λ), cold dark matter (CDM), and normal luminous matter. We use Angular Power Spectra (APS) to compare simulated ΛCDM models with observed large-scale galaxy distributions from the DECam Local Volume Exploration Survey (DELVE). After ensuring that the simulations can produce APS which match observational results, we vary cosmological parameters to constrain the allowed ranges of parameter values that are consistent with the observed data.

Feithen, William↗
DOE OSTI2024

Simulating cosmological evolution by quantum quench of an atomic Bose-Einstein condensate

In cosmological evolution, it is the homogeneous scalar field (inflaton) that drives the universe to expand isotropically and to generate standard model particles. However, to simulate cosmology, atomic gas research has focused on the dynamics of Bose-Einstein condensates (BEC) with continuously applied forces. In this paper we argue that a complementary approach needs also to be pursued; we thus consider the analog BEC experiments in a nondriven and naturally closed atomic system. We implement this using a BEC in an optical lattice which, after a quench, freely transitions from an unstable to a stable state. Furthermore, this dynamical evolution displays the counterpart “preheating,” “reheating,” and “thermalization” phases of cosmology. Importantly, our studies of these analog processes yield tractable analytic models. Additionally, of great utility to the cold atom community, such understanding elucidates the dynamics of nonadiabatic condensate preparation. Indeed, the dynamical processes discussed here are generic and in future cold atom reequilibration experiments it will be important to observe both the preheating stage, corresponding to a fragmented condensate, and the reheating stage, corresponding to a particle cloud.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Bridging the Gap between Cosmological Simulations with Graph Neural Networks and Domain Adaptation

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross dataset tasks (up to 28% better relative error and up to almost an order of magnitude better χ 2 ). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Accuracy of power spectra in dissipationless cosmological simulations

Abstract We exploit a suite of large N-body simulations (up to N = 40963) performed with abacus, of scale-free models with a range of spectral indices n, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At 1 per cent level it starts at early times from a fraction of $k_\Lambda$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, ${\sim}2{-}3 k_\Lambda$ at the very latest times we evolve to. At the $5{{\ \rm per\ cent}}$ level, accuracy extends up to wavenumbers of order $5k_\Lambda$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple n-dependence, allowing a extrapolation to place conservative bounds on the accuracy of N-body simulations of non-scale-free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small-scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_\Lambda$ may need revision in light of our convergence analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Comparing XRISM Cluster Velocity Dispersions with Predictions from Cosmological Simulations: Are Feedback Models Too Ejective?

The dynamics of the intracluster medium (ICM), the hot plasma that fills galaxy clusters, are shaped by gravity-driven cluster mergers and feedback from supermassive black holes (SMBHs) in the cluster cores. XRISM measurements of ICM velocities in several clusters offer insights into these processes. We compare XRISM measurements for nine galaxy clusters (Virgo, Perseus, Centaurus, Hydra-A, PKS 0745–19, A2029, Coma, A2319, and Ophiuchus) with predictions from three state-of-the-art cosmological simulation suites, TNG-Cluster, the Three Hundred Project GADGET-X, and GIZMO-SIMBA, that employ different models of feedback. In cool cores, XRISM reveals systematically lower velocity dispersions than the simulations predict, with all 10 measurements below the median simulated values by a factor of 1.5–1.7 on average and all falling within the bottom 10% of the predicted distributions. The observed kinetic-to-total pressure ratio is also lower, with a median value of 2.2%, compared to the predicted 5.0%–6.5% for the three simulations. Outside the cool cores and in non-cool-core (NCC) clusters, simulations show better agreement with XRISM measurements, except for the outskirts of the relaxed, cool-core cluster A2029, which exhibits an exceptionally low kinetic pressure support (<1%), with none of the simulated systems in either of the three suites reaching such low levels. The NCC Coma and A2319 exhibit dispersions at the lower end but within the simulated spread. Our comparison suggests that the three numerical models may overestimate the kinetic effects of SMBH feedback in cluster cores. Additional XRISM observations of NCC clusters will clarify if there is a systematic tension in the gravity-dominated regime as well.

galaxy clusters↗
DOE OSTI2025

Evaluating quenching in cosmological simulations of galaxy formation with spectral covariance in the optical window

ABSTRACT Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of the observables. Here, we follow an alternative approach based on the spectral covariance in a model-independent way. We build upon previous work by Sharbaf et al. that studied the covariance of high-quality SDSS (Sloan Digital Sky Survey) continuum-subtracted spectra in a relatively narrow range of velocity dispersion ($\sigma \in [100,150]$ km s$^{-1}$). Here, the same analysis is applied to synthetic data from the eagle and IllustrisTNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral covariance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age–metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with active galactic nucleus (AGN) feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid prescriptions, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasize that this methodology only relies on the spectral covariance to assess whether simulations provide a true representation of galaxy formation.

Sharbaf, Z. (ORCID:0009000450545946)↗
DOE OSTI2022

Farpoint: A High-resolution Cosmology Simulation at the Gigaparsec Scale

Abstract In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000 h −1 Mpc) 3 and evolves close to two trillion particles, corresponding to a mass resolution of m p ∼ 4.6 × 10 7 h −1 M ⊙ . These specifications enable comprehensive investigations of the galaxy–halo connection, capturing halos down to small masses. Further, the large volume resolves scales typical of modern surveys with good statistical coverage of high-mass halos. The simulation was carried out on the GPU-accelerated system Summit, one of the fastest supercomputers currently available. We provide specifics about the Farpoint run and present an initial set of results. The high mass resolution facilitates precise measurements of important global statistics, such as the halo concentration–mass relation and the correlation function down to small scales. Selected subsets of the simulation data products are publicly available via the HACC Simulation Data Portal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Snowmass2021 Cosmic Frontier White Paper: Cosmological Simulations for Dark Matter Physics

Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CMB-S4. The time is now to build a novel simulation program to interpret observations so that we can identify novel signatures of DM microphysics across a large dynamic range of length scales and cosmic time. This white paper identifies the key elements that are needed for such a simulation program. We identify areas of growth on both the particle theory side as well as the simulation algorithm and implementation side, so that we can robustly simulate the cosmic evolution of DM for well-motivated models. We recommend that simulations include a fully calibrated and well-tested treatment of baryonic physics, and that outputs should connect with observations in the space of observables. We identify the tools and methods currently available to make predictions and the path forward for building more of these tools. A strong cosmic DM simulation program is key to translating cosmological observations to robust constraints on DM fundamental physics, and provides a connection to lab-based probes of DM physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Testing the parametric model for self-interacting dark matter using matched halos in cosmological simulations

Here, we systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in Ref.Yang et al. (2023) with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and subhalos. We focus on the accuracy of the model in predicting halo density profiles at z = 0 and the evolution of maximum circular velocity. We find the model predictions agree with the simulations within 10%–50% for most of the simulated (sub)halos, 50%–100% for extreme cases. This indicates that the model effectively captures the gravothermal evolution of the halos with very strong, velocity-dependent self-interactions. For an example application, we apply the model to study the impact of various SIDM scenarios on strong lensing perturber systems, demonstrating its utility in predicting SIDM effects for small-scale structure analyses. Our findings confirm that the model is an effective tool for mapping CDM halos into their SIDM counterparts.

Yang, Daneng [Chinese Academy of Sciences, Nanjing↗
DOE OSTI2022

Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗