Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cosmological evolution”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Constraining accuracy of the pairwise velocities in N -body simulations using scale-free models

We present a continuation of an analysis that aims to quantify resolution of N-body simulations by exploiting large (up to N = 4096 3 ) simulations of scale-free cosmologies run using abacus. Here, we focus on radial pairwise velocities of the matter field, both by direct estimation and through the cumulative two-point correlation function (using the pair conservation equation). We find that convergence at the 1 per cent level of the mean relative pairwise velocity can be demonstrated over a range of scales, evolving from a few times the grid spacing at early times to slightly below this scale at late times. We show the analysis of two different box sizes as well as from averaging results from the smaller boxes, and compare the power of the two aforementioned estimators in constraining accuracy at each scale. Down to scales of the order of the smoothing parameter, convergence is obtained at ∼5 per cent precision, and shows a behaviour indicating asymptotic stable clustering. In conclusion, we also infer for LCDM simulations conservative estimates on the evolution of the lower cut-off to resolution (at 1 and 5 per cent precision) as a function of redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Relativistic Alfvén Waves Entering Charge-starvation in the Magnetospheres of Neutron Stars

Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current j A . We derive an analytic expression for the evolution of the wavevector k and the growth of j A . In the strongly sheared regime, j A may exceed the maximum current j0 that can be supported by the background e ± plasma. We investigate these charge-starved waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave is able to propagate successfully even when κ ≡ j A /j 0 $\ll$ 1. It sustains j A by compressing and advecting the plasma along the magnetic field lines with an increasing Lorentz factor, γ ≳ κ 1/2 . The simulations show how plasma instabilities lead to gradual dissipation of the wave energy. Our results suggest that an extremely high charge-starvation parameter κ ≳ 10 4 may be required in order for this mechanism to power the observed fast radio bursts (FRBs) from SGR 1935+2154. However, cosmological FRBs with much higher luminosities are unlikely to be a result of charge-starvation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

DESI DR1 Lyα 1D power spectrum: the optimal estimator measurement

The one-dimensional power spectrum P 1D of Lyα forest offers rich insights into cosmological and astrophysical parameters, including constraints on the sum of neutrino masses, warm dark matter models, and the thermal state of the intergalactic medium. We present the measurement of P 1D using the optimal quadratic maximum likelihood estimator applied to over 300,000 Lyα quasars from Data Release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI) survey. This sample represents the largest to date for P 1D measurements and is larger than the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) by a factor of 1.7. We conduct a meticulous investigation of instrumental and analysis systematics and quantify their impact on P 1D . This includes the development of a cross-exposure estimator that eliminates the need to model the pipeline noise and has strong potential for future P 1D measurements. We also present new insights into metal contamination through the 1D correlation function. Using a fitting function we measure the evolution of the Lyα forest bias with high precision: b F (z) = (-0.218 ± 0.002) × ((1 + z)/4) 2.96±0.06 . In a companion validation paper, we substantially extend our previous suite of CCD image simulations to quantify the pipeline's exquisite performance accurately. In another companion paper, we present DR1 P 1D measurements using the Fast Fourier Transform (FFT) approach to power spectrum estimation. These two measurements produce a forest bias parameter that differs by 2.2 sigma. However, our model is simplistic, so this disagreement will be investigated in future work.

Lyman alpha forest↗
DOE OSTI2026

The Fate of Nucleated Black Holes in de Sitter Quantum Gravity

The Euclidean Nariai geometry has long been proposed as the instanton describing the nucleation of maximal-mass black holes in de Sitter space. We place this interpretation on firmer footing by showing that, once an observer is included, the gravitational path integral produces the imaginary phase required for a transition rate. As a warmup, we revisit the Hawking-Moss instanton and, as a byproduct, find that scalar fields can enhance black-hole nucleation, suggesting a quantum-gravity bound on scalar potentials with de Sitter solutions. We then study the subsequent semiclassical evolution of the nucleated black hole. We show that the previously claimed "anti-evaporation" channel is unphysical, arising from a quantum state with singular horizons. In a smooth state, the black hole instead undergoes standard thermal Hawking evaporation. We verify explicit agreement with the no-boundary state and argue that this evaporation is not subject to large quantum-gravity corrections. The nucleated black hole thus evaporates completely back to the maximally-entropic empty de Sitter vacuum, making the full process a Boltzmann fluctuation.

FOS: Physical sciences↗
DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on ultralight axion dark matter through galaxy cluster number counts

Ultralight axions are hypothetical scalar particles that influence the evolution of large-scale structures of the Universe. Depending on their mass, they can potentially be part of the dark matter component of the Universe as candidates commonly referred to as fuzzy dark matter. While strong constraints have been established for pure fuzzy dark matter models, the more general scenario where ultralight axions constitute only a fraction of the dark matter has been limited to only a few observational probes. In this work, we use the galaxy cluster number counts obtained from the first All-Sky Survey (eRASS1) of the SRG/eROSITA mission together with gravitational weak lensing data from the Dark Energy Survey, the Kilo-Degree Survey, and the Hyper Suprime-Cam to constrain the fraction of ultralight axions in the mass range 10 −32 eV to 10 −24 eV. We put upper bounds on the ultralight axion relic density Ω a in independent logarithmic axion mass bins by performing a full cosmological parameter inference. We find an exclusion region in the intermediate ultralight axion mass regime with the tightest bounds reported so far in the mass bins around m a = 10 −27 eV with Ω a < 0.0035 and m a = 10 −26 eV with Ω a < 0.0079; both are at a 95% confidence level. When combined with cosmic microwave background probes, these bounds are tightened to Ω a < 0.0030 in the m a = 10 −27 eV mass bin and Ω a < 0.0058 in the m a = 10 −26 eV mass bin, with both at a 95% confidence level. This is the first time that constraints on ultralight axions have been obtained using the growth of structure measured by galaxy cluster number counts. These results pave the way for large surveys, which can be utilized to obtain tight constraints on the mass and relic density of ultralight axions with better theoretical modeling of the abundance of halos.

cosmological parameters↗
DOE OSTI2026

Precision Kinematic Sunyaev--Zel'dovich Measurements Across Halo Mass and Redshift with DESI DR2 and ACT DR6: Part I. Luminous Red Galaxies

We present the most precise measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect around luminous red galaxies to date, detecting the signal at $18σ$ significance in both harmonic and configuration space. Our analysis cross-correlates 2.4 million spectroscopic LRGs from the Dark Energy Spectroscopic Instrument (DESI) DR2 sample with Data Release 6 (DR6) of the Atacama Cosmology Telescope (ACT). We develop a novel harmonic-space cross-correlation approach using momentum-weighted kSZ templates, yielding nearly uncorrelated bandpowers within a framework consistent with other large-scale structure analyses. By incorporating the LRG halo occupation distribution (HOD) and its uncertainty, we convert measured galaxy gas profiles into halo gas profiles and provide generalized Navarro-Frenk-White (GNFW) fitting profiles, providing empirical targets for tuning feedback efficiency in hydrodynamical simulations and for baryonic modeling in large-scale structure analyses. We find strong evidence that gas profiles do not trace dark matter, providing direct evidence for gas redistribution beyond gravitational collapse. Comparing to hydrodynamical simulations, our measurements favor feedback efficiencies exceeding those in the Battaglia profile, suggesting more efficient gas ejection in group-scale halos than previously predicted. Splitting by redshift, we detect the kSZ signal at SNR $\approx 5$--$10$ in each of four bins and find amplitude evolution consistent with the expected decline in mean halo mass at fixed comoving number density. Splitting by stellar mass, we study the scaling of kSZ amplitude with galaxy properties. Together with BGS and ELG measurements in Paper II, these results span $0.1 \lesssim z \lesssim 1.6$ across three galaxy populations, demonstrating the potential of spectroscopic kSZ to map circumgalactic gas and constrain baryonic feedback.

Qu, F. J. [KIPAC, Menlo Park; Stanford U., HEPL; S↗
DOE OSTI2020

Biased tracers in redshift space in the EFTofLSS with exact time dependence

We study the effect of the Einstein-de Sitter (EdS) approximation on the one-loop power spectrum of galaxies in redshift space in the Effective Field Theory of Large-Scale Structure. The dark matter density perturbations and velocity divergence are treated with exact time dependence. Splitting the density perturbation into its different temporal evolutions naturally gives rise to an irreducible basis of biases. While, as in the EdS approximation, at each time this basis spans a seven-dimensional space, this space is a slightly different one, and the difference is captured by a single calculable time- and vec k-dependent function. We then compute the redshift-space galaxy one-loop power spectrum with the EdS approximation (P EdS-approx ) and without (P Exact ). For the monopole we find P 0 Exact /P 0 Ed S-approx ~ 1.003 and for the quadrupole P 2 Exact /P 2 EdS-approx ~ 1.007 at z=0.57, and sharply increasing at lower redshifts. Finally, we show that a substantial fraction of the effect remains even after allowing the bias coefficients to shift within a physically allowed range. This suggests that the EdS approximation can only fit the data to a level of precision that is roughly comparable to the precision of the next generation of cosmological surveys. Furthermore, we find that implementing the exact time dependence formalism is not demanding and is easily applicable to data. Both of these points motivate a direct study of this effect on the cosmological parameters.

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 OSTI2023

Hickson-like compact groups inhabiting different environments

ABSTRACT Although compact groups (CGs) of galaxies have been envisioned as isolated extremely dense structures in the Universe, it is accepted today that many of them could be not as isolated as thought. In this work, we study Hickson-like CGs identified in the Sloan Digital Sky Survey Data Release 16 to analyse these systems and their galaxies when embedded in different cosmological structures. To achieve this goal, we identify several cosmological structures where CGs can reside: nodes of filaments, loose groups, filaments, and cosmic voids. Our results indicate that 45 per cent of CGs do not reside in any of these structures, i.e. they can be considered non-embedded or isolated systems. Most of the embedded CGs are found inhabiting loose groups and nodes, while there are almost no CGs residing well inside cosmic voids. Some physical properties of CGs vary depending on the environment they inhabit. CGs in nodes show the largest velocity dispersions, the brightest absolute magnitude of the first-ranked galaxy, and the smallest crossing times, while the opposite occurs in non-embedded CGs. When comparing galaxies in all the environments and galaxies in CGs, CGs show the highest fractions of red/early-type galaxy members in most of the absolute magnitudes ranges. The variation between galaxies in CGs inhabiting one or another environment is not as significant as the differences caused by belonging or not to a CG. Our results suggest a plausible scenario for galaxy evolution in CGs in which both large-scale and local environments play essential roles.

Astronomy & Astrophysics↗
DOE OSTI2021

Field moment expansion method for interacting bosonic systems

In this work, we introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Atacama Cosmology Telescope: Probing the baryon content of SDSS DR15 galaxies with the thermal and kinematic Sunyaev-Zel’dovich effects

We present high signal-to-noise measurements (up to 12$\sigma$) of the average thermal Sunyaev Zel'dovich (tSZ) effect from optically selected galaxy groups and clusters and estimate their baryon content within a 2.1$^\prime$ radius aperture. Sources from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) DR15 catalog overlap with 3,700 sq. deg. of sky observed by the Atacama Cosmology Telescope (ACT) from 2008 to 2018 at 150 and 98 GHz (ACT DR5), and 2,089 sq. deg. of internal linear combination component-separated maps combining ACT and $\it{Planck}$ data (ACT DR4). The corresponding optical depths, $\bar{\tau}$, which depend on the baryon content of the halos, are estimated using results from cosmological hydrodynamic simulations assuming an AGN feedback radiative cooling model. We estimate the mean mass of the halos in multiple luminosity bins, and compare the tSZ-based $\bar{\tau}$ estimates to theoretical predictions of the baryon content for a Navarro-Frenk-White profile. We do the same for $\bar{\tau}$ estimates extracted from fits to pairwise baryon momentum measurements of the kinematic Sunyaev-Zel'dovich effect (kSZ) for the same data set obtained in a companion paper. We find that the $\bar{\tau}$ estimates from the tSZ measurements in this work and the kSZ measurements in the companion paper agree within $1\sigma$ for two out of the three disjoint luminosity bins studied, while they differ by 2-3$\sigma$ in the highest luminosity bin. The optical depth estimates account for one third to all of the theoretically predicted baryon content in the halos across luminosity bins. Potential systematic uncertainties are discussed. The tSZ and kSZ measurements provide a step towards empirical Compton-$\bar{y}$-$\bar{\tau}$ relationships to provide new tests of cluster formation and evolution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Pulsar timing anomalies: a window into baryon number violation

Here, we investigate the influence of a specific class of slow Baryon Number Violation (BNV)—one that induces quasi-equilibrium evolution — on pulsar spin characteristics. This work reveals how BNV can potentially alter observable parameters, including spin-down rates, the second derivative of spin frequency, and braking indices of pulsars. Moreover, we demonstrate that BNV could lead to anomalies in pulsar timing, along with a wide array of braking indices, both positive and negative. In addition, we examine the possibility of pulsar spin-up due to BNV, which may result in a novel mechanism for the revival of “dead” pulsars. We conclude by assessing the sensitivity required for future pulsar timing efforts to detect such BNV effects, thus highlighting the potential for pulsars to serve as laboratories for testing fundamental physics.

cosmology↗
DOE OSTI2024

Primordial black holes from bubble collisions during a first-order phase transition

We study the possibility of production of primordial black holes (PBHs) from bubble collisions during a first-order phase transition. While typical colliding bubbles are small and irrelevant for PBH production, we find that those that can produce PBHs must have a macroscopically thick fluid shell and have been born much before the typical nucleation time. Particularly large uncertainties arise from an exponential sensitivity of the nucleation rate on the required duration of bubble growth which depends on the details of the collisions and the evolution of the spacetime metric toward the end of the phase transition. We introduce a few parameters to be obtained from future numerical simulation to represent those unknowns, and estimate the PBH abundance in an Abelian Higgs benchmark model and show that it can be significant. We predict an approximately monochromatic PBH mass spectrum, and find regions in the parameter space where the PBHs can constitute entire dark matter or even overclose the Universe. Our result thus shows that models with a first-order phase transition can be constrained by overabundant PBHs or null results of other PBH searches. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Full calibration of the tomographic redshift distribution from the HSC PDR3 Shape Catalog with DESI

The calibration of tomographic redshift distributionsis essential for cosmological analysis of weak lensing data.In this work, we calibrate all four tomographic bins of the Hyper Suprime Camera (HSC) weak lensing catalog with the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 and 2 using the clustering redshifts technique. We include z > 1.2 redshift sources such as emission line galaxies (ELG) and quasars (QSO) sources in our calibration, which were not available in the previous HSC calibration (Rau et al. (2022), Mon. Not. Roy. Astron. Soc. 524 (2023) 5109), allowing a complete calibration of all the redshift bins. We find the first tomographic bin exhibits a small shift towards low redshifts. The second bin is in good agreement with the photometric calibration, while third and fourth bin exhibit a shift towards higher redshifts. However, these shifts are considerably smaller than the shifts obtained in the HSC Year 3 cosmic shear analyses. We evaluate the impact of galaxy bias and magnification effects from all the samples on the measurements, finding them to be small, and we propose corrections to reduce them further. Specifically, we relax the assumption of linear bias and only assume no redshift evolution of the cross-correlation coefficient, allowing us to leverage smaller clustering scales. We model the redshift distributions with splines and compare our results to previous analyses as well as to other parameterizations found in literature. For the two high-redshift tomographic bins, we find the shifts to higher redshifts with respect to the measurements performed in Rau+2022 to be Δz$_{3}$ =-0.039$^{+0.020}$$_{-0.021}$ and Δz$_{4}$ = -0.048$^{+0.012}$$_{-0.012}$.

Choppin de Janvry, J. [LBL, Berkeley; UC, Berkeley↗
DOE OSTI2023

Photometric Objects Around Cosmic Webs (PAC) Delineated in a Spectroscopic Survey. IV. High-precision Constraints on the Evolution of the Stellar–Halo Mass Relation at Redshift z < 0.7

Taking advantage of the Photometric objects Around Cosmic webs method developed in Paper I, we measure the excess surface density $\overline{n}$ 2 w p of the photometric objects around spectroscopic objects down to stellar masses 10 8.0 M ⊙ , 10 9.2 M ⊙ , and 10 9.8 M ⊙ in the redshift ranges of z s < 0.2, 0.2 < z s < 0.4, and 0.5 < z s < 0.7, respectively, using data from the DESI Legacy Imaging Surveys and spectroscopic samples from the Sloan Digital Sky Survey (i.e., the Main, LOWZ, and CMASS samples). We model the measured in an N-body simulation, using the abundance matching method, and we constrain the stellar–halo mass relations (SHMRs) in the three redshift ranges to percent levels. With the accurate modeling, we demonstrate that the stellar mass scatter for a given halo mass is nearly a constant, and that the empirical form of Behroozi et al. describes the SHMR better than the double-power-law form at low mass. Our SHMR accurately captures the downsizing of massive galaxies from z s = 0.7, while it also indicates that small galaxies are still growing faster than their host halos. The galaxy stellar mass functions (GSMFs) from our modeling are in perfect agreement with the model-independent measurements in Paper III, although the current work extends the GSMFs to a much smaller stellar mass. Based on the GSMFs and the SHMRs, we derive the stellar mass completeness and halo occupation distributions for the LOWZ and CMASS samples, which are useful for correctly interpreting their cosmological measurements, such as galaxy–galaxy lensing and redshift space distortion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Shaping dark photon spectral distortions

The cosmic microwave background (CMB) spectrum is an extraordinary tool for exploring physics beyond the Standard Model. The exquisite precision of its measurement makes it particularly sensitive to small effects caused by hidden sector interactions. In particular, CMB spectral distortions can unveil the existence of dark photons which are kinetically coupled to the standard photon. In this work, we use the COBE-FIRAS dataset to derive accurate and robust limits on photon-to-dark-photon oscillations for a large range of dark photon masses, from 10$^{−10}$ to 10$^{−4}$ eV. We consider in detail the redshift dependence of the bounds, computing CMB distortions due to photon injection/removal using a Green’s function method. Our treatment improves on previous results, which had set limits studying energy injection/removal into baryons rather than photon injection/removal, or ignoring the redshift evolution of distortions. The difference between our treatment and previous ones is particularly noticeable in the predicted spectral shape of the distortions, a smoking gun signature for photon-to-dark-photon oscillations. The characterization of the spectral shape is crucial for future CMB missions, which could improve the present sensitivity by orders of magnitude, exploring regions of the dark photon parameter space that are otherwise difficult to access [inline-graphic not available: see fulltext].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗