Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximations and Expansions”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quarkonium Polarization Kinetic Equation from Open Quantum Systems and Effective Field Theories

Recent measurements of polarization phenomena in relativistic heavy ion collisions have aroused a great interest in understanding dynamical spin evolution of the QCD matter. In particular, the spin alignment signature of J/ψ has been recently observed in Pb-Pb collisions at LHC, which may infer nontrivial spin transport of quarkonia in quark gluon plasmas. Motivated by this, we study the spin-dependent in-medium dynamics of quarkonia by using the potential nonrelativistic QCD (pNRQCD) and the open quantum system framework. By applying the Markovian approximation and Wigner transformation, we systematically derive the Boltzmann transport equation for vector quarkonia with polarization dependence in the quantum optical limit. As opposed to the previous study for the spin-independent case where the collision terms depend on chromoelectric correlators, the new kinetic equation incorporates gauge invariant correlators of chromomagnetic fields that determine the recombination and dissociation terms with polarization dependence at the order we are working in the multipole expansion. In the quantum Brownian motion limit, the Lindblad equation with new transport coefficients defined in terms of the chromomagnetic field correlators have also been derived. Our formalism is generic and valid for both weakly-coupled and strongly-coupled quark gluon plasmas. It may be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Searching for the QCD Critical Point Along the Pseudo-critical/Freeze-out Line Using Padé-resummed Taylor Expansions of Cumulants of Conserved Charge Fluctuations

Using high-statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we construct estimators for the radius of convergence from an eighth order series expansion of the pressure as well as the number density. We show that the estimator for pressure and number density will be identical in the asymptotic limit. In the vicinity of the pseudo-critical temperature, T pc ≃ 156.5 MeV, we find the estimator of the radius of convergence to be µ B /T ≳ 3 for strangeness-neutral matter. We also present results for the pole structure of the Pad´e approximants for the pressure at non-zero values of the baryon chemical potential and show that the pole structure of the [4,4] Pad´e is consistent with not having a critical point at temperatures larger than 135 MeV and a baryon chemical potential smaller than µ B /T ~ 2.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Spin–spin interactions in defects in solids from mixed all-electron and pseudopotential first-principles calculations

Abstract Understanding the quantum dynamics of spin defects and their coherence properties requires an accurate modeling of spin-spin interaction in solids and molecules, for example by using spin Hamiltonians with parameters obtained from first principles calculations. We present a real-space approach based on density functional theory for the calculation of spin-Hamiltonian parameters, where only selected atoms are treated at the all-electron level, while the rest of the system is described with the pseudopotential approximation. Our approach permits calculations for systems containing more than 1000 atoms, as demonstrated for defects in diamond and silicon carbide. We show that only a small number of atoms surrounding the defect needs to be treated at the all-electron level, in order to obtain an overall all-electron accuracy for hyperfine and zero-field splitting tensors. We also present results for coherence times, computed with the cluster correlation expansion method, highlighting the importance of accurate spin-Hamiltonian parameters for quantitative predictions of spin dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimal parameterizations for observational constraints on thawing dark energy

Time-varying dark energy is often modeled in observational analyses through generic parameterizations of its equation of state w(z), which typically use two free parameters {w 0 , w a } to span a broad range of behaviors as a function of redshift. However, this broad range of behaviors can only approximately capture the dynamics of any given microphysical theory of dark energy. A complementary approach is to use targeted parameterizations designed to model specific classes of dynamical dark energy with greater precision. Focusing on the class of thawing dark energy, we quantify and compare the precision with which nineteen generic and targeted parameterizations can capture the dynamics of physically motivated thawing quintessence theories. We find that a targeted parameterization derived from a Padé expansion of w is the most reliable of these, producing accurate reconstructions of w(z), the expansion history H(z), and cosmological parameters such as H 0 and Ω m for a broad range of microphysical theories.

dark energy experiments↗
DOE OSTI2023

P-V-T equation of state of boron carbide

We report the P-V-T equation of state measurements of B 4 C to 50 GPa and approximately 2500 K in laser-heated diamond anvil cells. We obtain an ambient temperature, third-order Birch–Murnaghan fit to the P-V data that yields a bulk modulus K 0 of 221(2) GPa and derivative, (dK/dP) 0 of 3.3(1). These were used in fits with both a Mie–Grüneisen–Debye model and a temperature-dependent, Birch–Murnaghan equation of state that includes thermal pressure estimated by thermal expansion (α) and a temperature-dependent bulk modulus (dK 0 /dT). The ambient pressure thermal expansion coefficient (α 0 + α 1 T), Grüneisen γ(V) = γ 0 (V/V 0 ) q and volume-dependent Debye temperature, were used as input parameters for these fits and found to be sufficient to describe the data in the whole P-T range of this study.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
DOE OSTI2022

The formation of discs in the interior of AGB stars from the tidal disruption of planets and brown dwarfs

ABSTRACT A significant fraction of isolated white dwarfs host magnetic fields in excess of a MegaGauss. Observations suggest that these fields originate in interacting binary systems where the companion is destroyed thus leaving a singular, highly magnetized white dwarf. In post-main-sequence evolution, radial expansion of the parent star may cause orbiting companions to become engulfed. During the common envelope phase, as the orbital separation rapidly decreases, low-mass companions will tidally disrupt as they approach the giant’s core. We hydrodynamically simulate the tidal disruption of planets and brown dwarfs, and the subsequent accretion disc formation, in the interior of an asymptotic giant branch star. Compared to previous steady-state simulations, the resultant discs form with approximately the same mass fraction as estimated but have not yet reached steady state and are morphologically more extended in height and radius. The long-term evolution of the disc and the magnetic fields generated therein require future study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Phase transition mechanism and bandgap engineering of Sb2S3 at gigapascal pressures

Abstract Earth-abundant antimony trisulfide (Sb 2 S 3 ), or simply antimonite, is a promising material for capturing natural energies like solar power and heat flux. The layered structure, held up by weak van-der Waals forces, induces anisotropic behaviors in carrier transportation and thermal expansion. Here, we used stress as mechanical stimuli to destabilize the layered structure and observed the structural phase transition to a three-dimensional (3D) structure. We combined in situ x-ray diffraction (XRD), Raman spectroscopy, ultraviolet-visible spectroscopy, and first-principles calculations to study the evolution of structure and bandgap width up to 20.1 GPa. The optical band gap energy of Sb 2 S 3 followed a two-step hierarchical sequence at approximately 4 and 11 GPa. We also revealed that the first step of change is mainly caused by the redistribution of band states near the conduction band maximum. The second transition is controlled by an isostructural phase transition, with collapsed layers and the formation of a higher coordinated bulky structure. The band gap reduced from 1.73 eV at ambient to 0.68 eV at 15 GPa, making it a promising thermoelectric material under high pressure.

Cui, Zhongxun (ORCID:0000000275044177)↗
DOE OSTI2023

Variational deep learning of equilibrium transition path ensembles

Here, we present a time-dependent variational method to learn the mechanisms of equilibrium reactive processes and efficiently evaluate their rates within a transition path ensemble. This approach builds off of the variational path sampling methodology by approximating the time-dependent commitment probability within a neural network ansatz. The reaction mechanisms inferred through this approach are elucidated by a novel decomposition of the rate in terms of the components of a stochastic path action conditioned on a transition. This decomposition affords an ability to resolve the typical contribution of each reactive mode and their couplings to the rare event. The associated rate evaluation is variational and systematically improvable through the development of a cumulant expansion. We demonstrate this method in both over- and under-damped stochastic equations of motion, in low-dimensional model systems, and in the isomerization of a solvated alanine dipeptide. In all examples, we find that we can obtain quantitatively accurate estimates of the rates of the reactive events with minimal trajectory statistics and gain unique insights into transitions through the analysis of their commitment probability.

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

Mechanical failure of fresh nuclear grade iron–chromium–aluminum (FeCrAl) cladding under simulated hot zero power reactivity initiated accident conditions

The reactivity initiated accident (RIA) is a postulated accident in light water reactors initiated by a rapid increase in reactivity which causes an increase in fission rate and fuel temperature. One potential mode of fuel system failure during RIA is pellet cladding mechanical interaction (PCMI) due to the rapid thermal expansion of the fuel pellet. Simulated PCMI experiments were performed by rapidly pressurizing iron-chromium-aluminum (FeCrAl) cladding tube samples using a hydraulic modified burst test system, causing the specimens to burst under a biaxial loading path. Deformation and rupture of the specimens were tracked with a telecentric lens and high-speed camera system. Outer surface strains were calculated using digital image correlation (DIC) on speckle patterns painted on the specimens’ outer surfaces. Experiments were conducted at approximately 275 °C, representative of hot zero power RIA conditions. The failure hoop strain was DIC-calculated between approximately 1.8–3.4%, corresponding to quasi-static energy depositions of approximately 110–260 calories per gram UO 2 assuming initial pellet-cladding contact. These conditions are very similar to the proposed energy deposition limit of 150 calories per gram UO 2 for unirradiated zirconium-based cladding in US Nuclear Regulatory Commission (NRC) Draft Regulatory Guide 1327. A very small strain rate dependence was observed in the data, with the magnitude of the failure hoop strain decreasing slightly with increasing strain rate. This observed dependence may be relevant because the FeCrAl cladding strain rate will be approximately 20% higher than in zirconium-based cladding due to the harder neutron spectrum and resultant shorter neutron generation time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗
DOE OSTI2025

DESI2024: Dark Energy Spectroscopic Instrument (DESI) survey, presentation of Year 1 results

The series of papers included in this special issue represent the first cosmological measurements from the Dark Energy Spectroscopic Instrument (DESI). These papers present the validation of the DESI clustering catalogs, the methodology for the cosmology measurements, and the full cosmology implications. These measurements were derived from a catalog of approximately 6.4 million extragalactic sources in a redshift range of 0.1 < z < 4.2. Owing to advances in fiber positioner technology and to a larger telescope, this sample was acquired in only one year of observation, yet is more than double the size of the galaxy and quasar samples from the first four phases of the Sloan Digital Sky Survey (SDSS) that spanned 20 years. With these new galaxy and quasar catalogs, the DESI collaboration completed percent-level constraints on the cosmic expansion history and the growth of structure. In turn, these new measurements provided a new platform from which to test the equation of state of dark energy, the summed mass of the three neutrino mass eigenstates, global curvature, and General Relativity. By surpassing the precision of the previous generation spectroscopic survey, the results in this special issue mark the transition in the field of cosmology to Stage-IV dark energy science.

baryon acoustic oscillations↗
DOE OSTI2020

Wide-angle effects for peculiar velocities

ABSTRACT The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as ‘wide-angle’ effects in redshift-space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work, we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in autocorrelation and in cross-correlation with the density field. Using both linear theory and the Zel'dovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low-redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration space, Fourier space, and spherical analyses; providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Realization of Non‐Equilibrium Wurtzite Structure in Heterovalent Ternary MgSiN 2 Film Grown by Reactive Sputtering

The piezoelectric and ferroelectric applications of heterovalent ternary materials are not well explored. Epitaxial MgSiN 2 films are grown at 600 °C on (111)Pt//(001)Al 2 O 3 substrates by the reactive sputtering method using metallic Mg and Si under the N 2 atmosphere. Detailed X-ray diffraction measurements and transmission electron microscopy observations revealed that the epitaxially grown films on the substrates have a hexagonal wurtzite structure with c-axis out-of-plane orientation. The random occupation of this structure by Mg and Si differs from that of the previously reported structure in which these two cations periodically occupy the cationic sites. However, the lattice spacings closely approximate those that are previously reported, irrespective of the ordering, and they are almost comparable with those of (Al 0.8 Sc 0.2 )N. The wide bandgap of >5.0 eV in deposited MgSiN 2 is compatible with that of AlN and suggests durability against the application of strong external electric fields, possibly to induce polarization switching. In addition, MgSiN 2 is shown to have piezoelectric properties with an effective d 33 value of 2.3 pm V −1 for the first time. This work demonstrates the compositional expansion of hexagonal wurtzite to heterovalent ternary nitrides for novel piezoelectric materials, whose ferroelectricity is expected.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tailoring the Electronic Structures and Spectral Properties of ZnO with Irradiation Defects Generated Under Intense Electronic Excitation: A Combined Experimental and DFT Approach

Here, the relationship between the composition of the internal defect states, spectral properties, and correlated electronic structures of wurtzite zinc oxide (ZnO) crystals under 645 MeV Xe 35+ irradiation is systematically investigated, employing experimental characterizations combined with first-principle calculations. Based on the ion irradiation-induced thermal expansion and relaxation processes, the high concentration of vacancy/interstitial defects produced from the transient disordered phase in molten track states trigger photoelectric changes, as follows: i) the generation of internal defect states effectively reduces the intrinsic bandgap (3.25 eV → 2.66 eV); ii) a large number of defective active sites inhibits the recombination between electron–hole pairs, causing dark conductance and photoconductance to increase with increasing damage levels until optimal fluence is achieved. Based on the density functional theory (DFT) with the GGA + U (GGA = generalized gradient approximation) method, the defective models associated with the different electronic structures, density of states, formation energy, and the nature of the chemical bonding are established. The narrowing of the bandgap observed experimentally and the enhancement of carrier concentration originating from the internal electron defect states are qualitatively verified, therefore laying the foundation for designing future nanoscale photoelectronic devices and microelectronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗