Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral function”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Beyond Contrast Transfer: Spectral SNR as a Finite-Dose Metric for STEM Phase Retrieval

The contrast transfer function (CTF) is widely used to evaluate phase retrieval methods in scanning transmission electron microscopy (STEM), including center-of-mass imaging, parallax imaging, direct ptychography, and iterative ptychography. However, the CTF reflects only the maximum usable signal, neglecting the effects of finite electron fluence and the Poisson-limited nature of detection. As a result, it can significantly overestimate practical performance, especially in low-dose regimes. Here, we employ the spectral signal-to-noise ratio (SSNR), as a finite-dose statistical framework to evaluate the recoverable signal as a function of spatial frequency. Using numerical reconstructions of white-noise objects, we show that center-of-mass, parallax, and direct ptychography exhibit dose-independent SSNRs, with close-form analytic expressions. In contrast, iterative ptychography exhibits a surprising dose dependence: at low fluence, its SSNR converges to that of direct ptychography; at high fluence, it saturates at a value consistent with the maximum detective quantum efficiency predicted by recent quantum Fisher information bounds. The results highlight the limitations of CTF-based evaluation and motivate SSNR as a more accurate, finite-dose metric for assessing STEM phase retrieval methods.

STEM phase retrieval↗
DOE OSTI2025

Remote quantification of Cm(III) and HNO 3 by fluorescence spectroscopy and chemometrics

A unique approach to remotely quantify Cm(III) (0–100 µg mL −1 ) in HNO 3 (1–12 M) using steady-state laser fluorescence spectroscopy and multivariate regression models was developed. Photoluminescence is amenable to remote measurements using fiber-optic cables and is sensitive to numerous lanthanide and actinide species. In-line measurements can provide feedback to support complex processing in harsh environments (e.g., hot cells) to help guide and optimize radiochemical separations. In this work, Cm(III) spectra were acquired remotely in a glove box as a function of HNO 3 concentration to better understand spectral characteristics and evaluate the utility of multivariate regression models in this system. Furthermore, the Cm(III) fluorescence peak shape, width, position, and intensity changed significantly as a function of HNO 3 concentration, likely because of the displacement of emission quenching inner-sphere water molecules and complexation with nitrate ions. Despite significant covariance and nonlinearity in the data, a D-optimal design strategy successfully minimized training set sample size and was used to build effective partial least squares regression models for Cm(III) and HNO 3 concentrations without a priori knowledge of solution conditions. Chemometrics for modeling complex fluorescence spectra are promising and may find widespread applicability for online analysis in numerous chemical systems found in the nuclear field.

Actinide↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2025

Understanding the Hydronium Cation in the Solid State: A Study in Synthetic Hydronium Uranyl Phosphate and Arsenate Mineral Systems and Their Irradiation Stability

Here, we present neutron diffraction, electronic structure calculations, and optical vibrational spectroscopic characterization of hydronium uranyl phosphate (HUP), the synthetic analog of chernikovite ((H 3 O)(UO 2 )(PO 4 )·3H 2 O) to gain insight into structural features of the hydronium cation in the solid state. HUP crystallizes in P2 1 /c instead of previously reported P4/ncc 1 at room temperature and Pccn below 302 K. Lower symmetry was required due to hydronium and interstitial water positions within the structure. Crystallographic positions of hydronium cations determined from powder neutron diffraction data are bolstered by Raman and attenuated total reflectance infrared spectroscopic measurements and are further informed by density functional theory with phonon eigenvector analysis for spectral assignments. Finally, HUP and its arsenate analog (trogerite, (H 3 O)(UO 2 )(AsO 4 )·3H 2 O), (HUAs) were studied using He 2+ irradiation as an analog for α radiolysis to investigate the irradiation stability of these phases and the stability of hydronium cations in the solid state. Dose studies were employed wherein each sample was irradiated to 5, 10, 15, 25, and 50 MGy. Structural insight regarding irradiated materials is gained using optical vibrational spectroscopy and powder X-ray diffraction. Surprising irradiation stability of HUAs was found up to 50 MGy of dose, which may have important implications for understanding and modeling the geologic stability of legacy U.

actinides↗
DOE OSTI2024

Flaring Stars in a Non-targeted mm-wave Survey with SPT-3G

We present a flare star catalog from four years of non-targeted millimeter-wave survey data from the South Pole Telescope (SPT). The data were taken with the SPT-3G camera and cover a 1500-square-degree region of the sky from $20^{h}40^{m}0^{s}$ to $3^{h}20^{m}0^{s}$ in right ascension and $-42^{\circ}$ to $-70^{\circ}$ in declination. This region was observed on a nearly daily cadence from 2019-2022 and chosen to avoid the plane of the galaxy. A short-duration transient search of this survey yields 111 flaring events from 66 stars, increasing the number of both flaring events and detected flare stars by an order of magnitude from the previous SPT-3G data release. We provide cross-matching to Gaia DR3, as well as matches to X-ray point sources found in the second ROSAT all-sky survey. We have detected flaring stars across the main sequence, from early-type A stars to M dwarfs, as well as a large population of evolved stars. These stars are mostly nearby, spanning 10 to 1000 parsecs in distance. Most of the flare spectral indices are constant or gently rising as a function of frequency at 95/150/220 GHz. The timescale of these events can range from minutes to hours, and the peak $\nu L_{\nu}$ luminosities range from $10^{27}$ to $10^{31}$ erg s$^{-1}$ in the SPT-3G frequency bands.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
DOE OSTI2024

57 Fe nuclear resonance vibrational spectroscopic studies of tetranuclear iron clusters bearing terminal iron( III )–oxido/hydroxido moieties

57 Fe nuclear resonance vibrational spectroscopy (NRVS) has been applied to study a series of tetranuclear iron ([Fe 4 ]) clusters based on a multidentate ligand platform (L 3− ) anchored by a 1,3,5-triarylbenzene linker and pyrazolate or (tertbutylamino)pyrazolate ligand (PzNH t Bu − ). These clusters bear a terminal Fe(III)–O/OH moiety at the apical position and three additional iron centers forming the basal positions. The three basal irons are connected with the apical iron center via a μ 4 -oxido ligand. Detailed vibrational analysis via density functional theory calculations revealed that strong NRVS spectral features below 400 cm −1 can be used as an oxidation state marker for the overall [Fe 4 ] cluster core. The terminal Fe(III)–O/OH stretching frequencies, which were observed in the range of 500–700 cm −1 , can be strongly modulated (energy shifts of 20–40 cm −1 were observed) upon redox events at the three remote basal iron centers of the [Fe 4 ] cluster without the change of the terminal Fe(III) oxidation state and its coordination environment. Therefore, the current study provides a quantitative vibrational analysis of how the remote iron centers within the same iron cluster exert exquisite control of the chemical reactivities and thermodynamic properties of the specific iron site that is responsible for small molecule activation.

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

Effects of screening and pressure ionization on the electron broadening of spectral lines in dense plasmas

Collisions between electrons and radiating atoms broaden spectral absorption and emission lines in dense plasmas. High densities also introduce screening and pressure ionization effects that distort the wave functions of both bound and free electrons. In order to study how dense plasma effects influence the electron broadening of spectral lines, this paper incorporates electron wave functions from an average-atom (AA) model to calculate the linewidth of the B III 2⁢𝑝−2⁢𝑠 transition at 𝑇 = 10 eV for mass densities ranging from 𝜌 = 10 −4 to 0.4⁢ g⁡/cc. The calculation method uses the impact approximation, allowing the linewidth to be written in terms of electron-collision cross sections and an interference term. Compared to an otherwise identical calculation that uses Coulomb free wave functions, the AA method is found to modify both the cross sections and the resulting linewidth at sufficiently high density by introducing screening and pressure-ionized bound states. Screening lowers the cross sections at low energies and near electron excitation thresholds, while pressure-ionized bound states introduce resonances into the continuum. Thus, as the density increases, the relative linewidth between the AA and Coulomb calculations follows a general decrease because of screening, with sharp increases at various intervals due to pressure ionization. Finally, the AA results are also compared with a common approach to introduce screening through the interaction potential and reduced models that use the Bethe formula for the inelastic electron-collision cross sections.

electronic excitation & ionization↗
DOE OSTI2026

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

What do physics-informed DeepONets learn? Understanding and improving training for scientific computing applications

Physics-informed deep operator networks (DeepONets) have emerged as a promising approach toward numerically approximating the solution of partial differential equations (PDEs). In this work, we aim to develop further understanding of what is being learned by physics-informed DeepONets by assessing the universality of the extracted basis functions and demonstrating their potential toward model reduction with spectral methods. Results provide clarity about measuring the performance of a physics-informed DeepONet through the decays of singular values and expansion coefficients. In addition, we propose a transfer learning approach for improving training for physics-informed DeepONets between parameters of the same PDE as well as across different, but related, PDEs where these models struggle to train well. This approach results in significant error reduction and learned basis functions that are more effective in representing the solution of a PDE.

Deep operator networks↗
DOE OSTI2024

PV Degradation Tools [SWR-20-71]

This open-source tool functions for calculating degradation of photovoltaic modules. For example, functions to calculate front and rear relative Humidity, as well as Acceleration Factors. A degradation calculation function is also being developed, considering humidity and spectral irradiances models.

Holsapple, Derek↗
DOE OSTI2025

Enhanced spectral purity of WSe 2 quantum emitters via conformal organic adlayers

Quantum emitters in solid-state materials are typically embedded in the bulk of their hosts, making their electronic transitions inaccessible to surface modification. In contrast, two-dimensional materials, with their all-surface nature, offer a platform for tuning quantum emitters via chemical functionalization. Because of its semiconducting properties that enable electrical addressability, monolayer WSe 2 is a promising candidate for quantum emission, although the complex interplay between point defects and the localized strain needed to activate quantum emission leads to poor spectral purity. Here, we demonstrate that functionalizing monolayer WSe 2 with conformal adlayers of 3,4,9,10-perylenetetracarboxylic dianhydride (PTCDA) improves quantum emission spectral purity. Optical spectroscopy reveals that PTCDA functionalization lowers defect activation energies by 10 meV and induces a 30 nm redshift in quantum emission wavelength, while preserving the bright and dark exciton energies of monolayer WSe 2 . First-principles calculations corroborate these findings, thus providing molecular-level insight into the underlying mechanism of enhanced spectral purity.

Ananth, Riddhi [Northwestern Univ., Evanston, IL (↗
DOE OSTI2024

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spectral Bounds on Hyperbolic 3-Manifolds: Associativity and the Trace Formula

We constrain the low-energy spectra of Laplace operators on closed hyperbolic manifolds and orbifolds in three dimensions, including the standard Laplace--Beltrami operator on functions and the Laplacian on powers of the cotangent bundle. Our approach employs linear programming techniques to derive rigorous bounds by leveraging two types of spectral identities. The first type, inspired by the conformal bootstrap, arises from the consistency of the spectral decomposition of the product of Laplace eigensections, and involves the Laplacian spectra as well as integrals of triple products of eigensections. We formulate these conditions in the language of representation theory of PSL 2 (C) and use them to prove upper bounds on the first and second Laplacian eigenvalues. The second type of spectral identities follows from the Selberg trace formula. We use them to find upper bounds on the spectral gap of the Laplace--Beltrami operator on hyperbolic 3-orbifolds, as well as on the systole length of hyperbolic 3-manifolds, as a function of the volume. Further, we prove that the spectral gap λ 1 of the Laplace--Beltrami operator on all closed hyperbolic 3-manifolds satisfies λ 1 < 47.32. Along the way, we use the trace formula to estimate the low-energy spectra of a large set of example orbifolds and compare them with our general bounds, finding that the bounds are nearly sharp in several cases.

Bonifacio, James [University of Mississippi, MS (U↗
DOE OSTI2026

Solid state 59 Co NMR study of a Np(VII) compound

Np(VII) compounds with [Co(NH 3 ) 6 ] 3+ cations were synthesized and structurally examined using powder X-ray diffraction. Multiple phases were observed, consisting of octahedral [Co(NH 3 ) 6 ] 3+ cations, discrete tetragonal bipyramidal [NpO 4 (OH) 2 ] 3− anions, and waters of hydration. Electric field gradient tensors at Co sites were measured by solid state 59 Co nuclear magnetic resonance (NMR) spectroscopy and compared with theoretical calculations. The relative contributions of the chemical shift and electric field tensors to the NMR lineshape were determined by recording spectra at field strengths of 7.04 and 11.74 Tesla. Further, the evolution of structure and morphology as a function of sample age and the effects on NMR spectral parameters has also been investigated. These results demonstrate the use of NMR at multiple fields to expand understanding of the stability and electronic structure of high valent neptunium compounds.

heptavalent neptunium↗
DOE OSTI2025

Excitation Laser Energy Dependence of the Gap-Mode TERS Spectra of WS 2 and MoS 2 on Silver

In this work we present a systematic study of the dependence of the gap mode tip-enhanced Raman scattering (TERS) response of the mono- and bi-layer WS 2 and MoS 2 on silver as a function of the excitation laser energy in a broad spectral range from 473nm to 830 nm. For this purpose, we collected consecutive TERS maps of the same area in the sample containing mono-and bi-layer regions with the same TERS probe with 6 different excitation lasers. To decrease the number of collected TERS maps, we used for the first time to the best of our knowledge, concurrent excitation and collection with two lasers simultaneously. We found that the E 2g /A 1g peak intensity ratio for the bilayer WS 2 @Ag and the ratio of the A’/A 1g peak intensity of the out-of-plane mode for the mono- and the bilayer, change in a significantly non-monotonous way as the excitation laser energy is swept from 1.58 eV to 2.62 eV. The former ratio increases at energies corresponding to A and B excitons (~2.0 eV and 2.4 eV correspondingly) in bilayer WS 2 . The absolute intensity of the A’ peak in the monolayer, and correspondingly the A’/A 1g ratio, is surprisingly high at lower excitation energies, but dips dramatically at the energy corresponding to the A exciton, being restored partially in between A and B excitons, but still showing the descending trend as the excitation laser energy increases. A somewhat similar picture was observed in mono- and bi-layers of MoS 2 @Ag, though the existing set of excitation lasers did not match the excitonic profile of this material as nicely as for the case of WS 2 . We attribute the observed behavior to the presence of intermediate (Fano resonance) or strong (Rabi splitting) coupling between the excitons in transition metal dichalcogenides (TMDs) and the plasmons in the tip-substrate nanocavity. This is akin to the so-called Fano (Rabi) transparency experimentally observed in far field scattering from transition metal dichalcogenides between two plasmonic metals. The possibility of the formation of intermediate/strong coupling between the excitonic resonances in TMDs and the nanocavity re-evaluates the role of various resonances in the gap-mode TERS and should become an important factor to be considered by TERS practitioners during planning the experiments. Lastly, based on observed phenomena and its explanation, we propose the “ideal” substrate for efficient TERS and tip enhanced photoluminescence (TEPL) measurements.

Raman spectroscopy↗
DOE OSTI2024

Fragment-based initialization for quantum subspace methods

Here, we present a novel quantum-classical algorithm called LAS-QKSD for multireference systems, by combining a classical localized active space (LAS) fragment-based multireference algorithm with the quantum Krylov subspace diagonalization (QKSD) method for quantum computers. The algorithm uses wave function information from a LAS self-consistent field (LASSCF) calculation to prepare an initial state with better overlap with the target ground state than the Hartree-Fock state. This is coupled with the use of QKSD to ultimately converge to the exact energy, providing faster convergence than starting from the Hartree-Fock state. Fragmentation has the two-fold benefit of fewer configurations on the classical side of the algorithm as well as fewer state preparation gates on the quantum side. First, we compare the LAS-QKSD method to the classical LASSCF method and to QKSD with a Hartree-Fock initial state. We then examine ways to load the LASSCF wave function using direct initialization and a QKSD-motivated spectral filtering approach. Finally, using a bimetallic complex, we show that the LAS-QKSD method is an efficient alternative to highly expensive complete active space SCF (CASSCF) calculations on strongly correlated systems.

D'Cunha, Ruhee↗