Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spectral simulations”

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

Pore-Scale Simulation of Spectral Induced Polarization

A new model was developed to simulate spectral induced polarization (SIP) in porous media. The model is based on a pore-network approximation of the soil or rock, in which bonds between pores are represented by equivalent circuits. In the circuit model, electrolytic and surface conduction occur in parallel, with the electrolyte behaving as a single resistor, and the mineral surface behaving as a resistor and capacitor in series. Based on the pore geometry, fluid conductivity, and mineral surface conductivity, an impedance network is generated. SIP spectra are calculated by solving the sinusoidal steady-state electrical problem for a range of frequencies. The new model provides a tool to bridge the gap between (1) our fundamental understanding of polarization processes at the scale of a single pore or grains and (2) existing simulation capabilities for macroscale systems that do not account for pore-scale properties. The model is implemented in Matlab and can solve two-dimensional and three-dimensional (2D and 3D) networks. Examples for regular cubic lattices are presented, but the method and code are flexible with respect to pore geometries and allow for future leveraging of existing open-source capabilities for generation of pore networks and extraction of network geometries from micro imaging of core samples.

58 GEOSCIENCES↗

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. Furthermore, the strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multielectronic and multiatomic effects in the U O 4,5 x-ray absorption spectroscopy of uranium dioxide

In order to explain all of the spectral features observed in the U O 4,5 x-ray absorption spectrum of uranium dioxide (UO 2 ), it is necessary to include both multielectron effects and multiatomic effects. The 5d (core hole)-5f (electron) angular momentum coupling that gives rise to the giant resonance has been treated within ligand field density functional theory, and the electron scattering that generates the extended x-ray absorption fine structure has been included via the spectral simulation program FEFF: both within a UO 8 fluorite cluster picture. An atomic model is insufficient to explain all of the observed spectral features.

36 MATERIALS SCIENCE↗

Parametric description of intermittent probability distribution functions in solar wind and magnetohydrodynamic turbulence

ABSTRACT In this work, we find empirical evidence that the scale-dependent statistical properties of solar wind and magnetohydrodynamic (MHD) turbulence can be described in terms of a family of parametric probability distribution functions (PDFs) known as Normal Inverse Gaussian (NIG). Understanding these PDFs is one of the most important goals in turbulence theory, as they are inherently connected to the intermittent properties of solar wind turbulence. We investigate the properties of PDFs of Elsasser increments based on a large statistical sample from solar wind observations and high-resolution numerical simulations of MHD turbulence. In order to measure the PDFs and their corresponding properties, three experiments are presented: fast and slow solar wind for experimental data and a simulation of reduced MHD (RMHD) turbulence. Conditional statistics on a 23-yr-long sample of WIND data near 1 au and high-resolution pseudo-spectral simulation of steadily driven RMHD turbulence on a $2048^3$ mesh are used to construct scale-dependent PDFs. The empirical PDFs are fitted to NIG distributions, which depend on four free parameters. Our analysis shows that NIG distributions accurately capture the evolution of the PDFs, with scale-dependent parameters, from large scales characterized by a Gaussian distribution, turning to exponential tails within the inertial range and stretched exponentials at dissipative scales. We also show that empirically-measured NIG parameters exhibit well-defined scaling properties that are similar across the three empirical data sets, which may be indicative of universal behaviour.

Astronomy & Astrophysics↗

Rovibronic molecular line list for the $N_2(C^3Π_u–B^3Π_g)$ second positive system

Here, a line list for the N second positive system, $B^3Π_g—C^3Π_u$, has been compiled using the PGOPHER spectral simulation software. The line list extends the number of vibrational states of the $B^3Π_g$ up to v=29 and a maximum rotational state of J=150 for simulation temperatures up to 7000 K. New electronic–vibrational transition moments were calculated using refined potential energy curves and a transition dipole moment with the DUO software. Comparisons to experimental data and the SPECAIR software have been used to validate the new line list. The results are available in ASCII ExoMol .state and .trans files and as a PGOPHER input file for use in spectral analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time-Resolved X-ray Emission Spectroscopy and Resonant Inelastic X-ray Scattering Spectroscopy of Laser Irradiated Carbon

The existence of liquid carbon as an intermediate phase preceding the formation of novel carbon materials has been a point of contention for several decades. Experimental observation of such a liquid state requires nonthermal melting of solid carbon materials at various laser fluences and pulse properties. Reflectivity experiments performed in the mid-1980s reached opposing conclusions regarding the metallic or insulating properties of the purported liquid state. Time-resolved X-ray absorption studies showed shortening of C–C bonds and increasing diffraction densities, thought to evidence a liquid or glassy carbon state, respectively. Nevertheless, none of these experiments provided information on the electronic structure of the proposed liquid state. Herein, we report the results of time-resolved resonant inelastic X-ray scattering (RIXS) and time-resolved X-ray emission spectroscopy (XES) studies on amorphous carbon (a-C) and ultrananocrystalline diamond (UNCD) as a function of delay time between the irradiating pulse and X-ray probe. For both a-C and UNCD, we attribute decreases in RIXS or XES signals to transition blocking, relaxation, and finally, ablation. Increased signal at 20 ps following the irradiation of the UNCD is attributed to the probable formation of nanoscale structures in the ablation plume. Further, differences in the amount of signal observed between a-C and UNCD are explained by the difference in sample thickness and, specifically, incomplete melting of the UNCD film. Comparisons to spectral simulations based on MD trajectories at extreme conditions indicate that the carbon state in our experiments is crystalline. Normal mode analysis confirmed that symmetrical bending or stretching of the C–C bonds in the diamond lattice results in XES spectra with small intensity differences. Overall, we observed no evidence of melting to a liquid state, as determined by the lack of changes in the spectral properties for up to 100 ps delays following the melting pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating resonant spectroscopy simulations using multishifted biconjugate gradient

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide important insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Prospects for detecting UF6 hydrolysis intermediates by mass spectrometry, resonance Raman, and NQR spectroscopy

Simulations are performed to consider several spectrometric and spectroscopic candidates for elucidating the mechanism of the hydrolysis of uranium hexafluoride (UF6). This study is among the first to benchmark the def-mTZVP basis sets for actinide-containing molecules, and it is shown to be a suitable basis set for surveying geometrical structures and vibrational spectra when used in conjunction with density functional theory. An experiment is proposed coupling mass spectrometric ion selection with vibrational spectroscopy, and supporting infrared and Raman spectral simulations demonstrate that ionization blue-shifts bands and can change their qualitative features. Ultraviolet resonance Raman is shown to have good prospects for discriminating U–O–U bridged intermediates, evidence for which was observed recently [L. E. McNamara et al., J. Phys. Chem. A 130, 775–786 (2026)]. As a more speculative approach, we also consider the prospect of addressing 233U or 235U nuclei by quadrupole resonance (NQR) spectroscopy for assigning early stage intermediate complexes. In doing so, we provide a first order-of-magnitude estimate for the collision-induced NQR signal for the UF6 dimer, which is of interest for the interpretation of a fast decoherence time observed in liquid-phase nuclear magnetic resonance. Having provided critical insights into the formation and stability of UF6 hydrolysis intermediates in previous studies, this computational spectroscopy survey is expected to help guide and expedite future laboratory campaigns.

Lutz, Jesse J. [Center for Computing Research, San↗

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗

High-Resolution Photoelectron Spectroscopy of NO3 – Vibrationally Excited Along Its ν3 Mode

The nitrate (NO3) radical has long been the subject of both experimental and theoretical studies due to its complex electronic structure resulting from vibronic interactions between its X̃2A2' and B̃2E' states. In particular, the definite assignment of the fundamental of its degenerate stretching vibration (ν3) is still under debate. Here, we report high-resolution photoelectron spectra of vibrationally pre-excited NO3- using the recently developed IR-cryo-SEVI technique. The anions are excited through infrared (IR) excitation near 1350 cm-1, accessing the ν3 and 2ν3(e') vibrational levels with band centers at 1350.5 and ∼2700 cm-1, respectively. The IR-cryo-SEVI spectrum for 2ν3 pre-excitation shows clear evidence for an intense 321 transition. From the position of this feature (30031 cm-1), the electron affinity of NO3 also determined in this work (31680 cm-1), and the IR excitation energy, we obtain a fundamental frequency of 1051 cm-1 for the ν3 fundamental of the NO3 radical. This assignment and other features in the IR-cryo-SEVI spectra are supported by spectral simulations based on a vibronic Köppel-Domcke-Cederbaum Hamiltonian. The simulations also show that nearly all features in the IR-cryo-SEVI spectra arise because of pseudo-Jahn-Teller coupling between the X̃ and B̃ states of NO3. The results and analysis presented here settle a long-standing controversy regarding the ν3 frequency of NO3.

Anions↗

Numerical study of pure rotational fs/ps CARS coherence beating at high pressure and for multi-species rotation-vibration non-equilibrium thermometry

Coherent anti-Stokes Raman scattering (CARS) has long been the gold standard for non-intrusively measuring gas temperature in reacting flows such as flames and plasmas. Recently, the development of ultrafast hybrid fs/ps CARS has enabled the exploitation of coherence beating between neighboring spectral lines to simultaneously measure rotational and vibrational temperatures from a single pure rotational spectrum. However, the influence of non-Boltzmann vibrational state distributions and limitations due to collisional dephasing at pressures greater than 1 atm remains unclear. In this work, we use spectral simulations to investigate the effects of non-Boltzmann vibrational state distributions and the applicability of coherence beating at pressures up to 10 atm. We show that short probe pulses can be leveraged to quantify non-Boltzmann vibrational state distributions of N2. Furthermore, we demonstrate that fs/ps CARS coherence beating can simultaneously provide sensitive measurements of rotational and vibrational temperatures of both O2 and N2 in air. A sensitivity analysis was conducted to qualitatively explain the accuracy and precision comparisons between probe delays.

Chen, Timothy Y. (ORCID:0000000330269767)↗

Adapting FEFF to 5f Angular Momentum Coupling

Here, it is demonstrated that the spectral simulation program FEFF can be adapted to include the effects of 5f total angular momentum coupling in the fluorite actinide dioxide systems ThO 2 , UO 2 , and PuO 2 . N 4,5 x-ray absorption spectra produced with this modified FEFF approach will be compared to the previous experimental results, obtaining a strong agreement between the two.

36 MATERIALS SCIENCE↗

Dataset describing two reference models for full-spectral lighting and daylight simulations together with implementations for two software systems

A dataset of two spectral lighting simulation reference models - one office and one factory hall - is presented. It aims to demonstrate and support full-spectral daylight and electric lighting simulations and facilitate evaluation of non-visual effects of light. The dataset includes Rhino CAD geometry, comprehensive spectral material and light source data and window system BSDF data. Example implementations in the two software tools, Radiance and OWL, enable reproducible workflows and support adoption in other software. The dataset is openly available on Zenodo. The office model reproduces Room 518 at the University of Innsbruck, including a west-facing façade and interior furnishings. The factory hall model follows the proposed geometry in the European standard 15193 for building energy performance. Interior reflectances in the office were measured in-situ using a handheld spectrometer. Exterior spectra and factory hall materials matching specified reflectances were obtained from an online spectral materials database. Glazing transmittance was derived from IGDB data using LBNL Optics/WINDOW. BSDFs for venetian blinds at various tilt angles, and for a diffusing pane adapted from the Complex Glazing Database, were generated in WINDOW. Luminaires in both models are specified with photometric files (Eulumdat/IES) and lamp spectra (Fluorescent 840, 4000 K LED). The provided example implementations (Radiance, OWL) include prepared input data and scripts to run first spectral simulations; example results are also included. The dataset is prepared to support reuse by researchers, designers and software developers for method validation, software engineering and comparison, and development of spectral metrics and controls.

Geisler-Moroder, David↗

Energy Exascale Computational Fluid Dynamics Simulations With the Spectral Element Method

Development and application of the open-source GPU-based fluid-thermal simulation code, NekRS, are described. Time advancement is based on an efficient kth-order accurate timesplit formulation coupled with scalable iterative solvers. Spatial discretization is based on the high-order spectral element method (SEM), which affords the use of fast, low-memory, matrix-free operator evaluation. Further, recent developments include support for nonconforming meshes using overset grids and for GPU-based Lagrangian particle tracking. Results of large-eddy simulations of atmospheric boundary layers for wind-energy applications as well as extensive nuclear energy applications are presented.

42 ENGINEERING↗

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Convergence of Cloud Droplet Spectral Relative Dispersion During Entrainment‐Mixing Based on Particle‐Resolved Direct Numerical Simulations

Entrainment-mixing processes critically impact cloud microphysical properties, but their effects on the relative dispersion (d) of cloud droplet size distributions (CDSDs) remain elusive. A direct numerical simulation model is initialized with different CDSDs to fill the gap. These results show that d decreases for broad CDSDs and increases for narrow ones, ultimately converging to approximately 0.5 regardless of initial CDSDs during the evaporation-dominated mixing stage. The supersaturation fluctuation and the shape of CDSDs jointly influence the convergence behavior of d. Further sensitivity tests show that the initial microphysical/dynamical/thermodynamical conditions exert negligible effects on the final converged value of d but affect the convergence rate (k). The k generally increases with increasing droplet number concentration and dissipation rate, and increases with decreasing liquid water content, relative humidity of entrained air, and mixing fraction of cloudy air. A conceptual model with two timescales is proposed; k and the timescales are negatively correlated, meaning that slow mixing and/or evaporation process results in slow convergence of d. In conclusion, this finding provides an important reference for improving understanding and parameterization of d during the entrainment-mixing processes.

54 ENVIRONMENTAL SCIENCES↗

Closing the Accuracy Gap in Tandem Photovoltaic Testing: An Accessible and Efficient Spectral Tuning Method Using LED-Based Simulators for Research Laboratories and Industry

Accurate performance calibration of multijunction (MJ) solar cells is critical for advancing this technology toward large-scale terrestrial application, yet existing testing methods developed by photovoltaics calibration laboratories remain prohibitively complex and resource intensive for most research laboratories. Current approaches rely on expensive multisource simulators and/or intricate spectral tuning algorithms, limiting accessibility and hindering standardized comparisons of emerging MJ technologies such as perovskite-based tandems. This paper introduces a streamlined spectral irradiance adjustment method for light-emitting diode (LED)-based solar simulators, which are increasingly adopted in the photovoltaics community due to their versatility and cost-effectiveness. The method we present bins LED channels into groups based on the number of junctions of the MJ photovoltaic device under test (DUT) and their corresponding band gaps and incorporates an automated tuning algorithm that eliminates the need to adjust each channel’s irradiance individually. This tuning algorithm requires the relative spectral responsivities of both the DUT and a broadband reference cell, as well as a calibrated spectroradiometer. Measurement validation across various MJ devices, including III-V and perovskite/silicon tandems, demonstrates excellent agreement within 1% with well-established xenon-tungsten multisource simulators and associated tuning algorithms. By enabling precise spectral tuning with readily available equipment and a simple tuning algorithm, our approach bridges the measurement accuracy gap between research laboratories and accredited testing facilities, fostering more reliable device comparisons and accelerating the translation of MJ technologies into real-world energy systems.

14 SOLAR ENERGY↗