SEARCH · Search NASA

Results for “Spectral moments”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Nonequilibrium spectral moment sum rules of the Holstein–Hubbard model

Abstract We derive a general procedure for evaluating the n th derivative of a time-dependent operator in the Heisenberg representation and employ this approach to calculate the zeroth to third spectral moment sum rules of the retarded electronic Green’s function and self-energy for a system described by the Holstein–Hubbard model allowing for arbitrary spatial and time variation of all parameters (including spatially homogeneous electric fields and parameter quenches). For a translationally invariant (but time-dependent) Hamiltonian, we also provide sum rules in momentum space. The sum rules can be applied to various different phenomena like time-resolved angle-resolved photoemission spectroscopy and benchmarking the accuracy of numerical many-body calculations. This work also corrects some errors found in earlier work on simpler models.

Physics↗

DOE OSTI2022

Constraints on ion velocity distributions from fusion product spectroscopy

Abstract Recent inertial confinement fusion experiments have shown primary fusion spectral moments which are incompatible with a Maxwellian velocity distribution description. These results show that an ion kinetic description of the reacting ions is necessary. We develop a theoretical classification of non-Maxwellian ion velocity distributions using the spectral moments. At the mesoscopic level, a monoenergetic decomposition of the velocity distribution reveals there are constraints on the space of spectral moments accessible by isotropic distributions. General expressions for the directionally dependent spectral moments of anisotropic distributions are derived. At the macroscopic level, a distribution of fluid element velocities modifies the spectral moments in a constrained manner. Experimental observations can be compared to these constraints to identify the character and isotropy of the underlying reactant ion velocity distribution and determine if the plasma is hydrodynamic or kinetic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Directional Variations in Tidal Flow Multifractality and Intermittency

We explored distinct directional variations of multifractal and intermittent characteristics of ebb and flood flow velocities at nodule point, WA, tidal energy site, and complementary inspection on the East River, NY using scaling exponents of the structure function, distribution flatness, detrending moving average (DMA) analysis, multifractal detrended fluctuation analysis (MF-DFA), and high-order spectral moments. Our findings reveal that tidal flow presents higher long-range dependence (LRD) and intermittent levels for the ebb flow sections, whereas the flood flow sections exhibit a higher degree of multifractality and greater sensitivity to larger magnitude of turbulent fluctuations. We demonstrate that long-range dependence predominantly contributes to multifractal behavior in both ebb and flood flows, as evidenced by the significantly reduced multifractal spectrum width for temporally randomly permuted time series. Moreover, spectral kurtosis analysis uncovers a higher intermittent level across all frequency scales for ebb flow sections and reveals a distinct pattern of tidal flow intermittency differing from the monotonically increasing intermittent level observed in wall-bounded and grid turbulence. Finally, we demonstrate that DMA, MF-DFA, and high-order spectral moments provide more comprehensive insights than structure function scaling exponents and PDF flatness methods.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗

DOE OSTI2026

How Proton Incorporation Reshapes Lattice Dynamics In BaSnO 3 ‐Type Proton Conductors

Proton conduction in acceptor-doped perovskites is fundamentally a vibronic process: mobile H + and D + do not move independently, but dynamically co-vibrate with the surrounding oxygen–metal framework. Direct experimental evidence for this behavior is presented using in situ 119⁢ 𝑆⁢𝑛 nuclear resonance vibrational spectroscopy (NRVS) on hydrated, deuterated, and dry 𝐵⁢𝑎⁢𝑆⁢𝑛 1−𝑥⁢ 𝑌 𝑥 ⁢𝑂 3−𝛿 . Hydration induces systematic redistributions in the Sn-projected phonon density of states (PDOS), including an upshift of the first spectral moment by about 0.4 meV, indicating a stiffening of the extended Sn–O network. H/D isotopic substitution leaves the Sn-projected PDOS largely unchanged, with only subtle isotope-dependent spectral reweighting, demonstrating that protonic degrees of freedom are not localized oscillators but are embedded in collective lattice modes. These results are rationalized using a classical coupled proton–phonon oscillator model that links the observed PDOS variations to changes in effective force constants and vibrational mass terms. The model captures how H + and D + participate in cooperative lattice dynamics rather than forming isolated OH/OD entities. Overall, NRVS probes proton–lattice coupling in ceramic proton conductors and quantitatively describes how protonic defects modulate host lattice dynamics to enable phonon-assisted long-range proton transport.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Faster spectral density calculation using energy moments

Abstract Accurate predictions of inclusive scattering cross sections in the linear response regime require efficient and controllable methods to calculate the spectral density in a strongly-correlated many-body system. In this work we reformulate the recently proposed Gaussian Integral Transform technique in terms of Fourier moments of the system Hamiltonian which can be computed efficiently on a quantum computer. One of the main advantages of this framework is that it allows for an important reduction of the computational cost by exploiting previous knowledge about the energy moments of the spectral density. For a simple model of medium mass nucleus like $$^{40}$$ 40 Ca and target energy resolution of 1 MeV we find an expected speed-up of $$\approx 125$$ ≈ 125 times for the calculation of the giant dipole response and of $$\approx 50$$ ≈ 50 times for the simulation of quasi-elastic electron scattering at typical momentum transfers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Spectral evidence for local-moment ferromagnetism in the van der Waals metals $\mathrm{Fe_3GaTe_2}$ and $\mathrm{Fe_3GeTe_2}$

Magnetism in two-dimensional (2D) materials has attracted considerable attention recently for both fundamental understanding of magnetism and its tunability towards device applications. The isostructural Fe 3 GeTe 2 and Fe 3 GaTe 2 are two members of the Fe-based van der Waals (vdW) ferromagnet family, but exhibit very different Curie temperatures (T C ) of 210 and 360 K, respectively. In this report by using angle-resolved photoemission spectroscopy and density functional theory, we systematically compare the electronic structures of the two compounds. Qualitative similarities in the Fermi surface can be found between the two compounds, with expanded hole pockets in Fe 3 GaTe 2 suggesting additional hole carriers compared to Fe 3 GeTe 2 . Interestingly, we observe almost no band shift in Fe 3 GaTe 2 across its T C of 360 K, compared to a small shift in Fe 3 GeTe 2 across its T C of 210 K. The weak temperature-dependent evolution strongly deviates from the expectations of an itinerant Stoner mechanism. Our results suggest that itinerant electrons have minimal contributions to the enhancement of T C in Fe 3 GaTe 2 compared to Fe 3 GeTe 2 , and that the nature of ferromagnetism in these Fe-based vdW ferromagnets must be understood with considerations of the electron correlations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗

DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗

DOE OSTI2023

Hadronization dynamics from the spectral representation of the gauge invariant quark propagator

We study the Dirac decomposition of the gauge invariant quark propagator, whose imaginary part describes the hadronization of a quark as this interacts with the vacuum, and relate each of its coefficients to a specific sum rule for the chiral-odd and chiral-even quark spectral functions. Working at first in light-like axial gauge, we obtain a new sum rule for the spectral function associated to the gauge fixing vector, and show that its second moment is in fact equal to zero. Then, we demonstrate that the first moment of the chiral-odd quark spectral function is equal in any gauge to the so-called inclusive jet mass, which is related to the mass of the particles produced in the hadronization of a quark. Finally, we present a gauge-dependent formula that connects the second moment of the chiral-even quark spectral function to invariant mass generation and final state rescattering in the hadronization of a quark.

Costa, Caroline S. R.↗

DOE OSTI2025

Micro Rain Radar Pro Data at the Argonne Testbed for Multiscale Observational Science obtained during the CROCUS Urban Integrated Field Laboratory

The Micro Rain Radar Pro (MRR-PRO) is a vertically pointing Ka-band Doppler radar designed to capture the fine-scale structure and evolution of precipitation. By recording the full Doppler spectrum at high temporal and spatial resolution, the MRR-PRO provides insight into both hydrometeor fall velocities and precipitation microphysics. From these spectra, key moments—reflectivity, mean Doppler velocity, spectral width, and rainfall rate—are derived and stored alongside the raw spectral data in CF/Radial 1.4-compliant files. Deployed at the Argonne Testbed for Multiscale Observational Studies (ATMOS) since November 2024, the MRR-PRO delivers vertical profiles at 70 m range resolution extending up to 4.5 km above ground level. These observations enable detailed analyses of precipitation type, intensity, and vertical structure, supporting process-level studies of cloud and precipitation dynamics in diverse weather regimes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗

DOE OSTI2024

Data Science Enabled Enabled Discovery of Superconductors (Final Progress Report)

This Final Technical Report describes efforts by 4 PIs at the University of Florida (Peter Hirschfeld, Richard Hennig, Greg Stewart and James Hamlin), over the period September 2019-August 2023, to use data science and machine learning techniques to discover new conventional superconductors. The PIs constructed a discovery loop with two theorists and two experimentalists to: develop algorithms to machine learn descriptors correlating strongly with the critical temperature Tc (PI's Peter Hirschfeld, UF Physics and Richard Hennig, UF Materials Science and En), synthesize and measure properties of promising materials, and feed back the knowledge gained into the prediction algorithm. This work was motivated by the theoretical prediction and experimental discovery of high-pressure, high-pressure hydride superconductors, and to find ways to recreate the high critical temperatures in these systems at ambient pressure. Highlights from the grant include: 1) a new equation for Tc in terms of moments of the electron-phonon spectral function, improving on the so-called Allen-Dynes equation (1975); 2) study of the metastable A15 superconductor Nb3Si, formed under explosive compression at ~1000GPa to determine the kinetic barrier to the ground state structure; 3) the development of ultra-fast machine-learned atomic potentials for molecular dynamics, and 4) the discovery of superconductivity at 19K in WB2 arising from metastable defect structures in the crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

DOE OSTI2024

A Regional Phase Amplitude Model of 2-D Attenuation for North America

We analyzed seismic attenuation patterns across the North America using over 70,000,000 Lg wave amplitudes recorded by the various networks at frequencies of 0.1-32 Hz. Our inversion solved for laterally varying attenuation, site terms, moments, and apparent stress following Phillips et al., (2016). The inversion was anchored by independently constrained: corner frequencies (via coda spectral ratios) to control the attenuation-stress tradeoff and moment measurements of teleseismic (GCMT, USGS) and regional (St. Louis University and UC Berkeley) earthquakes to provided absolute scaling. The quality factor (Q) shows clear regional patterns: low values in coastal, volcanic, and tectonically active regions, and high values in stable areas like the Great Plains and major plateaus throughout North America. These 2-D Q models enable improved regional source characterization, magnitude estimation, and yield determination.

58 GEOSCIENCES↗

DOE OSTI2025

Exact spectral gaps of random one-dimensional quantum circuits

The spectral gap of local random quantum circuits is a fundamental property that determines how close the moments of the circuit's unitaries match those of a Haar random distribution. When studying spectral gaps, it is common to bound these quantities using tools from statistical mechanics or via quantum information-based inequalities. Here, by focusing on the second moment of one-dimensional unitary circuits where nearest-neighboring gates act on sets of qudits (with open and closed boundary conditions), we show that one can exactly compute the associated spectral gaps. Indeed, having access to their functional form allows us to prove several important results, such as the fact that the spectral gap for closed boundary condition is exactly the square of the gap for open boundaries, as well as improve on previously known bounds for approximate design convergence. Finally, we verify our theoretical results by numerically computing the spectral gap for systems of up to 70 qubits, as well as comparing them to gaps of random orthogonal and symplectic circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Anapole moment of neutrinos and radioactive sources near liquid xenon detectors

We show that placing a radioactive source such as 51 Cr near a liquid xenon detector may allow us to detect the contribution induced by the anapole moment to neutrino-electron scattering in the Standard Model (SM) at the 1 − 2⁢𝜎 level. Although the anapole moment of neutrinos induces a scattering rate with the same spectral shape as the neutral and charged current contributions, exposures of ∼ 60 ton × source run at XENONnT or XLZD may be enough to accumulate sufficient statistics for a detection. We also discuss a simple model where the anapole moment of neutrinos is enhanced or decreased with respect to the SM expectation, further demonstrating how a potential measurement of the anapole moment of neutrinos would allow us to constrain new physics.

dark matter direct detection↗