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

Position-dependent Voronoi probability distribution functions for matter and halos

Here, we measure the Voronoi density probability distribution function (PDF) for both dark matter and halos in N-body simulations. For the dark matter, Voronoi densities represent the matter density field smoothed on a uniform mass scale, which approximates the Lagrangian density field. For halos, the Voronoi densities contain information about the local environment of each halo. We measure the halo virial masses, the total amount of dark matter within each halo Voronoi cell, and the halo Voronoi cell volumes, and we show how halo abundances depend on these three quantities. We then study the position-dependent Voronoi density PDF, measured within finite subregions of the Universe, using separate universe simulations. We demonstrate that the spatial variation of the position-dependent PDF is due to large-scale density fluctuations, indicating that the position-dependent PDF is a biased tracer of large-scale structure. We measure this bias for the dark matter, and interpret it as the bias of regions of the Lagrangian density field that are selected based on density. For the halos, this bias can be interpreted as a form of assembly bias. We present the mapping from late-time to early-time Voronoi density for each simulation dark matter particle, which is highly stochastic. We compare the median of this stochastic map with spherical collapse calculations and discuss challenges involved in modeling the evolution of the density field on these scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

ForestFlow: predicting the Lyman-α forest clustering from linear to nonlinear scales

On large scales, the Lyman-α forest provides insights into the expansion history of the Universe, while on small scales, it imposes strict constraints on the growth history, the nature of dark matter, and the sum of neutrino masses. This work introduces ForestFlow, a novel framework that bridges the gap between large- and small-scale analyses, which have traditionally relied on distinct modeling approaches. Using conditional normalizing flows, ForestFlow predicts the two Lyman-α linear biases (b δ and b η ) and six parameters describing small-scale deviations of the three-dimensional flux power spectrum (P 3D ) from linear theory as a function of cosmology and intergalactic medium physics. These are then combined with a Boltzmann solver to make consistent predictions, from arbitrarily large scales down to the nonlinear regime, for P 3D and any other statistics derived from it. Trained on a suite of 30 fixed-and-paired cosmological hydrodynamical simulations spanning redshifts from z = 2 to 4.5, ForestFlow achieves 3 and 1.5% precision in describing P 3D and the one-dimensional flux power spectrum (P 1D ) from linear scales to k = 5 Mpc −1 and k ∥ = 4 Mpc −1 , respectively. Thanks to its conditional parameterization, ForestFlow shows similar performance for ionization histories and two ΛCDM model extensions – massive neutrinos and curvature – even though none of these are included in the training set. This framework will enable full-scale cosmological analyses of Lyman-α forest measurements from the DESI survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Using action space clustering to constrain the recent accretion history of Milky Way-like galaxies

ABSTRACT In the currently favoured cosmological paradigm galaxies form hierarchically through the accretion of satellites. Since a satellite is less massive than the host, its stars occupy a smaller volume in action space. Actions are conserved when the potential of the host halo changes adiabatically, so stars from an accreted satellite would remain clustered in action space as the host evolves. In this paper, we identify recently disrupted accreted satellites in three Milky Way-like disc galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellites through simulation snapshots. We try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the z = 0 snapshot. Even with completely error-free mock data we find that only 35 per cent (14/39) satellites are well recovered while the rest (25/39) are poorly recovered (i.e. either contaminated or split up). Most (10/14 ∼70 per cent) of the well-recovered satellites have infall times <7.1 Gyr ago and total mass >4 × 108M⊙ (stellar mass more than 1.2 × 106 M⊙, although our upper mass limit is likely to be resolution dependent). Since cosmological simulations predict that stellar haloes include a population of in situ stars, we test our ability to recover satellites when the data include 10–50 per cent in situ contamination. We find that most previously well-recovered satellites stay well recovered even with 50 per cent contamination. With the wealth of 6D phase space data becoming available we expect that cluster analysis in action space will be useful in identifying the majority of recently accreted and moderately massive satellites in the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Quantum maximum entropy closure for small flavor coherence

Quantum angular moment transport schemes are an important avenue toward describing neutrino flavor mixing phenomena in dense astrophysical environments such as supernovae and merging neutron stars. Successful implementation will require new closure relations that go beyond those used in classical transport. In this paper, we derive the first analytic expression for a quantum M1 closure, valid in the limit of small flavor coherence, based on the maximum entropy principle. We verify that the resulting closure relation has the appropriate limits and characteristic speeds in the diffusive and free-streaming regimes. We then use this new closure in a moment linear stability analysis to search for fast flavor instabilities in a binary neutron star merger simulation and find better results as compared with previously designed, ad hoc , semiclassical closures.

astrophysical & cosmological simulations↗
DOE OSTI2025

Moving lens effect: Simulations, forecasts, and foreground mitigation

The peculiar motion of massive objects across the line of sight imprints a dipolar temperature anisotropy pattern on the cosmic microwave background known as the moving lens effect. This effect provides a unique probe of the transverse components of the peculiar velocity field, but has not yet been detected due to its small size. We implement and validate a stacking estimator for the moving lens signal using a galaxy catalog as a tracer of massive haloes combined with reconstructed velocities from the galaxy number density field. Using simulations, we forecast detection prospects for the moving lens signal from current and upcoming microwave background and galaxy surveys. Here, we demonstrate a new foreground mitigation strategy likely sufficient for current datasets, and discuss various sources of systematic error and noise. Upcoming galaxy surveys will provide high-significance statistical detections of the moving lens effect.

Astrophysical & cosmological simulations↗
DOE OSTI2025

Impact of neutrino flavor conversions on neutron star merger dynamics, ejecta, nucleosynthesis, and multimessenger signals

We present numerical relativity simulations of binary neutron star mergers incorporating neutrino flavor transformations triggered by fast flavor instability, quantum many-body effects, or potential beyond standard model physics. In both long-lived and short-lived remnant scenarios, neutrino flavor conversions modify species-dependent neutrino luminosities and mean energies, and drive the matter towards more neutron rich conditions. They produce up to more neutron rich ejecta and significantly boost the r-process yields, especially in low-density, near-equatorial outflows. We identify regions unstable to fast flavor instabilities and find that these instabilities persist despite flavor conversions. We further test the sensitivity to the equilibration timescale of the flavor conversions, finding that slower flavor conversions can interact with thermodynamic equilibration, and increase the neutron richness of the ejecta. Flavor conversions may also contribute to stronger gravitational wave and neutrino emissions, pointing to a correlation between neutrino transport and merger dynamics. Furthermore, these results highlight the potential impact of flavor conversions while motivating future work to improve on theoretical understanding of flavor instabilities in global simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2023

Collective neutrino oscillations on a quantum computer with hybrid quantum-classical algorithm

We simulate the time evolution of collective neutrino oscillations in two-flavor settings on a quantum computer. We explore the generalization of Trotter-Suzuki approximation to time-dependent Hamiltonian dynamics. The trotterization steps are further optimized using the Cartan decomposition of two-qubit unitary gates U ϵ SU(4) in the minimum number of controlled-NOT (CNOT) gates making the algorithm more resilient to the hardware noise. As a result, a more efficient hybrid quantum-classical algorithm is also explored to solve the problem on noisy intermediate-scale quantum devices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Neutrino Flavor Transformation in Neutron Star Mergers

Here, we present the first numerical relativity simulations including neutrino flavor transformations that could result from flavor instabilities, quantum many-body effects, or potential beyond standard model physics in neutron star mergers. We find that neutrino flavor transformations impact the composition and structure of the remnant, potentially leaving an imprint on the postmerger gravitational-wave signal. They also have a significant impact on the composition and nucleosynthesis yields of the ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Magnetic flux emergence in binary neutron star remnants

Using high-resolution AthenaK simulations of a twisted toroidal flux tube, we study the flux emergence of magnetic structures in the shear layer of a hot massive neutron star typical of a binary neutron star remnant. High-resolution simulations demonstrate that magnetic buoyant instabilities allow for emergence only for extremely large magnetic fields significantly exceeding 10 17 G, and more typical fields around 10 16 G are instead dominated by hydrodynamic effects. Because merger remnants tend to be stable against hydrodynamic convection, our work places strong limitations on the mechanisms by which massive binary neutron star remnants can produce the magnetically-driven outflows needed to power jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The DESI One-Percent survey: exploring the Halo Occupation Distribution of Emission Line Galaxies with AbacusSummit simulations

The One-Percent survey of the Dark Energy Spectroscopic Instrument collected ~ 270k emission line galaxies (ELGs) at 0.8 < z < 1.6. The high completeness of the sample allowed the clustering to be measured down to scales never probed before, 0.04 Mpc/h in rp for the projected 2-point correlation function (2PCF) and 0.17 Mpc/h in galaxy pair separation s for the 2PCF monopole and quadrupole. The most striking feature of the measurements is a strong signal at the smallest scales, below 0.2 Mpc/h in rp and 1 Mpc/h in s. We analyse these data in the halo occupation distribution framework. We consider different distributions for central galaxies, a standard power law for satellites with no condition on the presence of a central galaxy and explore several extensions of these models. For all considered models, the mean halo mass of the sample is found to be log10 $\langle$Mh$\rangle$ ~ 11.9. We obtain a satellite mean occupation function which agrees with physically motivated ELG models only if we introduce central-satellite conformity, meaning that the satellite occupation is conditioned by the presence of central galaxies of the same type. To achieve in addition a good modelling of the clustering between 0.1 and 1 Mpc/h in rp, we allow for ELG positioning outside of the halo virial radius and find 0.5% of ELGs residing in the outskirts of halos. Furthermore, the satellite velocity dispersion inside halos is found to be ~ 30% larger than that of the halo dark matter particles. These are the main findings of our work. Here we investigate assembly bias as a function of halo concentration, local density or local density anisotropies and observe no significant change in our results. We split the data sample in two redshift bins and report no significant evolution with redshift. Lastly, changing the cosmology in the modelling impacts only slightly our results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamical attractors in contracting spacetimes dominated by kinetically coupled scalar fields

We present non-perturbative numerical relativity simulations of slowly contracting spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear σ model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble density fluctuations fully consistent with cosmic microwave background observations. We show that the non-linear evolution rapidly approaches a homogeneous, isotropic and flat Friedmann-Robertson-Walker (FRW) geometry for a wide range of inhomogeneous and anisotropic initial conditions. Ultimately, we find, the kinetic coupling causes the evolution to deflect away from flat FRW and towards a novel Kasner-like stationary point, but in general this occurs on time scales that are too long to be observationally relevant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

FLAMINGO: combining kinetic SZ effect and galaxy–galaxy lensing measurements to gauge the impact of feedback on large-scale structure

Energetic feedback processes associated with accreting supermassive black holes can expel gas from massive haloes and significantly alter various measures of clustering on $\sim$Mpc scales, potentially biasing the values of cosmological parameters inferred from analyses of large-scale structure (LSS) if not modelled accurately. Here, we use the state-of-the-art FLAMINGO suite of cosmological hydrodynamical simulations to gauge the impact of feedback on large-scale structure by comparing to Planck + ACT stacking measurements of the kinetic Sunyaev–Zel’dovich (kSZ) effect of SDSS BOSS galaxies. We make careful like-with-like comparisons to the observations, aided by high precision KiDS and DES galaxy–galaxy lensing measurements of the BOSS galaxies to inform the selection of the simulated galaxies. In qualitative agreement with several recent studies using dark matter only simulations corrected for baryonic effects, we find that the kSZ effect measurements prefer stronger feedback than predicted by simulations which have been calibrated to reproduce the gas fractions of low redshift X-ray-selected groups and clusters. We find that the increased feedback can help to reduce the so-called $S_8$ tension between the observed and CMB-predicted clustering on small scales as probed by cosmic shear (although at the expense of agreement with the X-ray group measurements). However, the increased feedback is only marginally effective at reducing the reported offsets between the predicted and observed clustering as probed by the thermal SZ (tSZ) effect power spectrum and tSZ effect – weak lensing cross-spectrum, both of which are sensitive to higher halo masses than cosmic shear.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

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

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

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Non-parametric Lagrangian biasing from the insights of neural nets

We present a Lagrangian model of galaxy clustering bias in which we train a neural net using the local properties of the smoothed initial density field to predict the late-time mass-weighted halo field.By fitting the mass-weighted halo field in the AbacusSummit simulations at z = 0.5, we find that including three coarsely spaced smoothing scales gives the best recovery of the halo power spectrum. Adding more smoothing scales may lead to 2–5% underestimation of the large-scale power and can cause the neural net to overfit.We find that the fitted halo-to-mass ratio can be well described by two directions in the original high-dimension feature space.Projecting the original features into these two principal components and re-training the neural net either reproduces the original training result, or outperforms it with a better match of the halo power spectrum. The elements of the principal components are unlikely to be assigned physical meanings, partly owing to the features being highly correlated between different smoothing scales.Our work illustrates a potential need to include multiple smoothing scales when studying galaxy bias, and this can be done easily with machine-learning methods that can take in high dimensional input feature space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

General relativistic hydrodynamic simulations of binary strange star mergers

We perform fully general-relativistic simulations of binary strange star mergers considering two different approaches for thermal effects. The first uses a cold equation of state (EOS) derived from a modified version of the MIT bag model which is then supplemented by a Γ-law correction. The second approach employs a microphysical description of the finite-temperature effects. We describe results obtained with the two treatments, highlighting the influence of thermal effects. We find that the postmerger dynamics differs significantly in the two cases, leading to quantitative differences in the postmerger gravitational-wave spectrum and ejecta mass. The peak frequency of the postmerger gravitational-wave emission is consistent with the established quasi-universal relations for binary neutron star mergers and as a result, our simulations cannot distinguish between mergers of neutron stars and those of strange stars. Our models with realistic treatment of finite-temperature effects produce a significant amount of ejecta ≳0.02 M ⊙ ​. Here, the resulting flux of strangelets near the Earth, computed assuming that all neutron star mergers are in fact strange-stars mergers and that the binary considered here is representative, is in tension with experimental upper limits. As such, our results tentatively disfavor a scenario in which strange-quark matter is the lowest energy state of matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Fluctuation–dissipation relation for a quantum Brownian oscillator in a parametrically squeezed thermal field

In this paper we study the nonequilibrium evolution of a quantum Brownian oscillator, modeling the internal degree of freedom of a harmonic atom or an Unruh–DeWitt detector, coupled to a nonequilibrium and nonstationary quantum field bath and inquire whether a fluctuation–dissipation relation (FDR) can exist after/if it approaches equilibration. This is a nontrivial issue because a squeezed field bath cannot reach equilibration and yet, as this work shows, the system oscillator indeed can, which is a necessary condition for FDRs. We discuss three different settings: (A) The bath field essentially remains in a squeezed thermal state throughout, whose squeeze parameter is a mode- and time-independent constant. This situation is often encountered in quantum optics and quantum thermodynamics. (B) The bath field is initially in a thermal state, but is subjected to a parametric process leading to mode- and time-dependent squeezing. This scenario is encountered in cosmology and dynamical Casimir effects. The squeezing in the bath in both types of processes will affect the oscillator’s nonequilibrium evolution. We show that at late times it approaches equilibration and this stationarity condition warrants the existence of a FDR. The trait of squeezing is marked by the oscillator’s effective equilibrium temperature, and the proportionality factor in the FDR is only related to the stationary component of the noise kernel of the bath field. Setting (C) is more subtle: A finite system–bath coupling strength can set the oscillator in a squeezed state even though the bath field is stationary and does not engage in any parametric process. The squeezing of the system in this case is in general time-dependent but becomes constant when the internal dynamics is fully relaxed. We begin with comments on the broad range of physical processes involving squeezed thermal baths and end with some remarks on the significance of FDRs in capturing the essence of quantum backreaction in nonequilibrium and stochastic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The electroweak phase transition: a collider target

Determining the thermal history of electroweak symmetry breaking (EWSB) is an important challenge for particle physics and cosmology. Lattice simulations indicate that EWSB in the Standard Model (SM) occurs through a crossover transition, while the presence of new physics beyond the SM could alter this thermal history. The occurrence of a first order EWSB transition would be particularly interesting, providing the needed pre-conditions for generation of the cosmic matter-antimatter asymmetry and sources for potentially observable gravitational radiation. I provide simple, generic arguments that if such an alternate thermal history exists, the new particles involved cannot be too heavy with respect to the SM electroweak temperature, nor can they interact too feebly with the SM Higgs boson. These arguments do not rely on the decoupling limit. I derive corresponding quantitative expectations for masses and interaction strengths which imply that their effects could in principle be observed (or ruled out) by prospective next generation high energy colliders. The simple, generic arguments provide a quantitative, parametric understanding of results obtained in a wide range of explicit model studies; relate them explicitly to the electroweak temperature; and delineate broad contours of collider phenomenology pertaining to a non-standard history of EWSB.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗