Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density perturbation”

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 OSTI2025

Beyond the two-point correlation: Constraining primordial non-Gaussianity with density-perturbation moments

Constraining primordial non-Gaussianity (PNG) on the large-scale cosmic structure (LSS) is an important step in understanding properties of the early Universe, specifically in distinguishing between different inflationary models. Measuring PNG relies on evaluating the scale-dependent correlations in the density field. New summary statistics beyond the two- and three-point correlation functions in configuration space and their Fourier-space counterparts, the power- and bispectrum may provide increased sensitivity. We introduce a new method for extracting the PNG signal imprinted on the LSS by using the first three Gaussian moments of the normalized correlation in density perturbations, evaluated on varying distance scales. We aim to assess this method’s sensitivity to local PNG, parameterized by f NL . We performed spherical convolutions on a range of scales on dark-matter-halo simulations to measure the scale-dependent correlations in the density field. From these, we computed the first three moments and compared them to a model expectation vector, parameterized to the second power in f NL . Our method provides about 21% improvement in sensitivity to f NL with respect to using the two-point correlation function alone. Notably, we find that the second moment alone carries nearly as much constraining power as the mean, highlighting the potential of higher order statistics. Given its simplicity and efficiency, this framework is well suited for application to current and upcoming large-scale surveys such as the Dark Energy Spectroscopic Instrument (DESI).

early universe↗
DOE OSTI2025

Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities

Perturbative quantum chromodynamics (pQCD), while valid only at densities exceeding those found in the cores of neutron stars, could provide constraints on the dense-matter equation of state (EOS). Here, in this work, we examine the impact of pQCD information on the inference of the EOS using a nonparametric framework based on Gaussian processes (GPs). We examine the application of pQCD constraints through a ``pQCD likelihood,'' and verify the findings of previous works; namely, a softening of the EOS at the central densities of the most massive neutron stars and a reduction in the maximum neutron-star mass. Although the pQCD likelihood can be easily integrated into existing EOS inference frameworks, this approach requires an arbitrary selection of the density at which the constraints are applied. The EOS behavior is also treated differently on either side of the chosen density. To mitigate these issues, we extend the EOS model to higher densities, thereby constructing a ``unified'' description of the EOS from the neutron-star crust to densities relevant for pQCD. In this approach the pQCD constraints effectively become part of the prior. Since the EOS is unconstrained by any calculation or data between the densities applicable to neutron stars and pQCD, we argue for maximum modeling flexibility in that regime. We compare the unified EOS with the traditional pQCD likelihood, and although we confirm the EOS softening, we do not see a reduction in the maximum neutron-star mass or any impact on macroscopic observables. Though residual model dependence cannot be ruled out, we find that pQCD suggests the speed of sound in the densest neutron-star cores has already started decreasing toward the asymptotic limit; we find that the speed of sound squared at the center of the most massive neutron star has an upper bound of $\sim 0.5$ at the 90% level.

equations of state of nuclear matter↗
DOE OSTI2021

New methodology for measuring electron density perturbations caused by plasma coherent modes using profile reflectometry: Magnitudes and radial profiles in DIII-D

New capabilities of fast-sweep frequency-modulated profile reflectometry are explored to measure electron density ne perturbation magnitudes and radial profiles due to plasma coherent modes in DIII-D. The first approach is based on the frequency analysis of phase perturbations associated with high frequency (~MHz) Alfvén eigenmodes (AEs). The measurement of ~5.5 MHz fast-ion-driven global Alfvén eigenmodes (GAEs) is demonstrated in a neutral beam-heated DIII-D plasma. The GAE induced a broad radial distribution of phase perturbations in the profile reflectometer data. Analysis of these data determined the effective cutoff location displacement and the estimated ne fluctuation profile. In the second approach, high resolution ne profiles are used directly to determine the radial structure of ne perturbations due to a neo-classical tearing mode. These new measurements broaden the application of profile reflectometry and advance the development of AE spectroscopy as a tool for non-invasive diagnosis of fast-ion-driven modes in DIII-D and burning plasmas such as ITER.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

The growth of density perturbations in the last ~10 billion years from tomographic large-scale structure data

In order to investigate the origin of the ongoing tension between the amplitude of matter fluctuations measured by weak lensing experiments at low redshifts and the value inferred from the cosmic microwave background anisotropies, we reconstruct the evolution of this amplitude from z ~ 2 using existing large-scale structure data. To do so, we decouple the linear growth of density inhomogeneities from the background expansion, and constrain its redshift dependence making use of a combination of 6 different data sets, including cosmic shear, galaxy clustering and CMB lensing. We analyze these data under a consistent harmonic-space angular power spectrum-based pipeline. We show that current data constrain the amplitude of fluctuations mostly in the range 0.2 < z < 0.7, where it is lower than predicted by . This difference is mostly driven by current cosmic shear data, although the growth histories reconstructed from different data combinations are consistent with each other, and we find no evidence of systematic deviations in any particular experiment. In spite of the tension with , the data are well-described by the model, albeit with a lower value of S 8 ≡ σ 8 (Ω m /0.3) 0.5 . As part of our analysis, we find constraints on this parameter of S 8 = 0.7781 ± 0.0094 (68% confidence level), reaching almost percent-level errors comparable with CMB measurements, and 3.4σ away from the value found by .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Thermodynamic properties of metaschoepite predicted from density functional perturbation theory

Density functional perturbation theory (DFPT) calculations of the thermodynamic properties of metaschoepite, (UO 2 ) 8 0 2 (OH) 12 •10H 2 0, are reported. Using a recently revised crystal structure of metaschoepite, the predicted molar entropy and isobaric heat capacity are overall significantly smaller than previous calculations using an earlier orthorhombic crystal structure model. The present DFPT calculations also show large differences between the thermo-dynamic functions of metaschoepite and schoepite, which might reflect the change in phonon properties upon removal of two H 2 0 molecules per formula unit and alteration of the H-bonded interlayer water network from schoepite to metaschoepite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Computational Protocol to Evaluate Electron–Phonon Interactions Within Density Matrix Perturbation Theory

We present a computational protocol, based on density matrix perturbation theory, to obtain non-adiabatic, frequency-dependent electron–phonon self-energies for molecules and solids. Our approach enables the evaluation of electron–phonon interaction using hybrid functionals, for spin-polarized systems, and the computational overhead to include dynamical and non-adiabatic terms in the evaluation of electron–phonon self-energies is negligible. We discuss results for molecules, as well as pristine and defective solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Polymorphism and phase transitions in Na 2 U 2 O 7 from density functional perturbation theory

Polymorphism and phase transitions in sodium diuranate, Na 2 U 2 O 7 , are investigated with density functional perturbation theory (DFPT). Thermal properties of crystalline α-, β- and γ-Na 2 U 2 O 7 polymorphs are predicted from DFPT phonon calculations, i.e., the first time for the high-temperature γ-Na 2 U 2 O 7 phase (R$\bar{3}$ with combining macron]m symmetry). The standard molar isochoric heat capacities predicted within the quasi-harmonic approximation are C p (298.15K)= 219.4 and 220.9 J K -1 mol -1 mfor P21/a α-Na 2 U 2 O 7 and C2/m β-Na 2 U 2 O 7 , respectively. Gibbs free energy calculations reveal that α-Na 2 U 2 O 7 (P2 1 /a) and β-Na 2 U 2 O 7 (C2/m) are almost energetically degenerate at low temperature, with β-Na 2 U 2 O 7 becoming slightly more stable than α-Na 2 U 2 O 7 as temperature increases. Here these findings are consistent with XRD data showing a mixture of α and β phases after cooling of γ-Na 2 U 2 O 7 to room temperature and the observation of a sluggish α → β phase transition above ca. 600 K. A recently observed α-Na 2 U 2 O 7 structure with P21 symmetry is also shown to be metastable at low temperature. Based on Gibbs free energy, no direct β → γ solid-solid phase transition is predicted at high temperature, although some experiments reported the existence of such phase transition around 1348 K. This, along with recent experiments, suggests the occurrence of a multi-step process consisting of initial β-phase decomposition, followed by recrystallization into γ-phase as temperature increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Excitation of long electrostatic solitary waves in ion beam neutralization process

Unusually long electrostatic solitary waves (ESWs) are discovered in a particle-in-cell simulation study of the process of ion beam neutralization by electron emission from a filament. These ESWs are long because the density perturbation responses to the potential wells created by the ESWs are very small. Notably, the density perturbation is small because the trapped (positive) and untrapped (negative) electron density perturbations nearly compensate each other because of a non-Maxwellian electron velocity distribution in the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Observation of quasi-coherent density fluctuation in scrape-off layer enhancing boundary transport in high-βN hybrid plasmas on DIII-D

In this work, we report the observation of a quasi-coherent density fluctuation (QCF) by the Doppler backscattering system in the scrape-off layer (SOL) region of the DIII-D tokamak. This QCF is observed in high power, high performance hybrid plasmas with near double-null divertor (DND) shape during the electron cyclotron heating period. This mode is correlated with a steepened SOL density profile, and leads to significantly elevated particle and heat fluxes between ELMs. The SOL QCF is a long-wavelength ion-scale fluctuation (k θ ρ s ≈0.2-0.4 and k r ρ s ≈0.03), and propagates in the ion diamagnetic direction in the plasma frame. Its radial expanse is about 1.5–2 cm, well beyond the typical width of heat flux λ q on DIII-D. Also, the SOL QCF does not show any clear dependence on the effective SOL collisionality, and thus may raise issues on the control of plasma-material interactions in low collisionality plasmas in which the blob-induced transport is reduced. A linear simulation using BOUT++ with a 5-field reduced model is performed and compared with experimental observations. In simulation results, an interchange-like density perturbation can be driven by the SOL density gradient, and its peak location and the radial width of the density perturbation are in agreement with the experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Second-harmonic generation tensors from high-throughput density-functional perturbation theory

Optical materials play a key role in enabling modern optoelectronic technologies in a wide variety of domains such as the medical or the energy sector. Among them, nonlinear optical crystals are of primary importance to achieve a broader range of electromagnetic waves in the devices. However, numerous and contradicting requirements significantly limit the discovery of new potential candidates, which, in turn, hinders the technological development. In the present work, the static nonlinear susceptibility and dielectric tensor are computed via density-functional perturbation theory for a set of 579 inorganic semiconductors. The computational methodology is discussed and the provided database is described with respect to both its data distribution and its format. Several comparisons with both experimental and ab initio results from literature allow to confirm the reliability of our data. The aim of this work is to provide a relevant dataset to foster the identification of promising nonlinear optical crystals in order to motivate their subsequent experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Towards a quantum fluid theory of correlated many-fermion systems from first principles

Correlated many-fermion systems emerge in a broad range of phenomena in warm dense matter, plasmonics, and ultracold atoms. Quantum hydrodynamics (QHD) complements first-principles methods for many-fermion systems at larger scales. We illustrate the failure of the standard Bohm potential central to QHD for strong perturbations when the density perturbation is larger than about 10^{-3} 10 − 3 of the mean density. We then extend QHD to this regime via the many-fermion Bohm potential from first-principles. This may lead to more accurate QHD simulations beyond their common application domain in the presence of strong perturbations at scales unattainable with first-principles methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Rotation of electrothermal-instability-driven overheating structure due to helically oriented surface magnetic field on a high-current-density aluminum rod

Experiments on the 1-MA, 100-ns-rise-time Mykonos Facility demonstrate rotation of electrothermal instability (ETI)-driven overheating structure on 1.00-mm-diameter, 10-nm-surface roughness, 99.999%-pure aluminum rods, which are pulsed with helically polarized surface magnetic field. Rods are machined to include pairs of 10-micron-scale quasi-hemispherical voids or “engineered defects (ED)” which provide the dominant current density perturbation from which ETI grows most rapidly. Experiments include an axial magnetic field component through the addition of a helically wound return-current electrode or “helical return can (HRC).” For a given HRC design, azimuthal field (B ɵ ) and axial field (B z ) components rise at a prescribed and fixed ratio, driving an increasing magnetic field of constant polarization at the rod's surface; most experiments generated surface magnetic field at a 15-degree field polarization angle (from horizontal) defined as ɸ B = arctan(B z /B ɵ ). ETI-driven emission patterns from individual ED are observed to rotate along ɸ B , while emission patterns from dielectric-coated ED pairs are shown to elongate and preferentially merge along ɸ B , in qualitative agreement with 3D-magnetohydrodynamic simulations. These data strongly support that for a randomized distribution of current density perturbations on a high-current density conductor, nearby perturbations will favorably merge about ɸ B , with the degree of merging increasing with current. Such observations offer fundamental new understanding of the seeding mechanisms of the helical magneto-Rayleigh Taylor (MRT) instabilities observed from axially magnetized magnetically driven imploding liners.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On anomalous transport of multi-species plasma associated with the resistive ballooning and resistive drift waves driven turbulence

Anomalous transport of multi-species plasma related to the resistive ballooning and resistive drift wave turbulence is considered in a “cold” ion approximation. It is found that similar to the resistive drift wave turbulence [see A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024); and S. I. Krasheninnikov and R. D. Smirnov, Phys. Plasmas (to be published)] the addition of the ballooning drive does not change the main features of anomalous transport of the multi-species plasma: (i) The transport of all ion species is described as a transport of the passive scalars in the turbulent field of the electrostatic potential and electron density perturbation; (ii) the density of ion species with a larger ratio of the mass to charge has the tendency to the accumulation/depletion in the vortices of plasma flow; and (iii) the cross-field transport of all plasma species (including electrons and ions) is described by the same anomalous transport coefficient.

Physics↗
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 OSTI2020

(U) Energy-Dependent Perturbations of Scattering Cross Sections Using the PERT Card in MCNP6.2

The PERT card is used for first- and second-order perturbations of material densities, isotope densities, and isotope cross sections in the MCNP6.2 Monte Carlo code. Using the PERT card for firstand second-order sensitivities is straightforward. Reference 2 only looked at density perturbations with no energy dependence. This report looks at the sensitivity of a response to scattering cross sections as a function of energy and finds that there is a bug in the PERT capability. The SENSMG multigroup neutron sensitivity code was recently given the capability of computing sensitivities to user-specified reactions whose cross sections are available in a user-supplied NJOY output file. SENSMG uses the PARTISN multigroup discrete ordinates code for the transport. The NJOY-reaction capability is similar to that of SUSD3D, which also uses PARTISN for the multigroup transport. The PERT bug was identified by comparing MCNP6.2 outputs with SENSMG outputs. The next section of this report describes the test problem. Section III presents the MCNP6.2 PERT bug and proposes a solution. Section IV is a summary and conclusions. The SENSMG and MCNP input files are listed in the appendix.

07 ISOTOPE AND RADIATION SOURCES↗