Search NASA⌕ Search

SEARCH · Search NASA

Results for “regularization”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Interpolation of compound semiconductor alloy parameters from those of their constituents

Several methods have been proposed for interpolation of the value of physical parameters of quaternary alloys from those of their constituent ternary and binary sub-alloys. These expressions agree when non-linear bowing terms are not required; they differ in how the bowing terms of the bounding ternaries should be utilized. Common interpolation expressions for quaternaries can be generalized into two groups: (1) those that use a linear interpolation of the nearest ternary parameter values and (2) those that interpolate over binary values with a bowing term derived from the bounding ternaries. The second group of methods is equivalent to a polynomial expansion over the alloy’s interpolation space. For compound semiconductor alloys, the geometry of the composition space is the direct sum of the group-III and group-V mixture sub-spaces. The mixture sub-spaces are best described using barycentric coordinates on a regular simplex. A general polynomial expansion of the value of an alloy parameter using barycentric coordinates for the group-III and group-V simplex spaces is described along with an algorithm to generate interpolation expressions for alloys with arbitrary numbers of elements, including quinary and senary alloys. It is shown that a polynomial expansion produces values in closer agreement with the direct gap of quaternaries lattice-matched to common substrates than do approaches using an interpolation of the ternary values, despite a prominent recommendation to the contrary. Finally, a quaternary correction term is described that improves the predicted direct bandgap energies of GaInAsSb for compositions near those lattice matched to InP, InAs, and GaSb.

Olesberg, Jonathon T. [Sandia National Laboratorie↗
DOE OSTI2025

Analysis methodology of coherent oscillations in time- and angle-resolved photoemission spectroscopy

Oscillatory signals from coherently excited phonons are regularly observed in ultrafast pump–probe experiments on condensed matter samples. Electron–phonon coupling implies that coherent phonons also modulate the electronic band structure. These oscillations can be probed with energy and momentum resolution using time- and angle-resolved photoemission spectroscopy (trARPES), which reveals the orbital dependence of the electron–phonon coupling for a specific phonon mode. However, a comprehensive analysis remains challenging when multiple coherent phonon modes couple to multiple electronic bands. Complex spectral line shapes due to strong correlations in quantum materials add to this challenge. In this work, we examine how the frequency domain representation of trARPES data facilitates a quantitative analysis of coherent oscillations of the electronic bands. Here, we investigate the frequency domain representation of the photoemission intensity and the first moment of the energy distribution curves. Both quantities provide complementary information and are able to distinguish oscillations of binding energy, linewidth, and intensity. We analyze a representative trARPES dataset of the transition metal dichalcogenide WTe 2 and construct composite spectra, which intuitively illustrate how much each electronic band is affected by a specific phonon mode. We also show how a linearly chirped probe pulse can generate extrinsic artifacts that are distinct from the intrinsic coherent phonon signal.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Molecular aspect ratio effect on axial thermal transport in solution-spun carbon nanotube fibers

Neat, densely packed, and highly aligned carbon nanotube fibers (CNTFs) have appealing room-temperature axial thermal conductivity (⁠k⁠) and thermal diffusivity (⁠α⁠) for applications in lightweight heat spreading, flexible thermal connections, and thermoelectric active cooling. Although CNTFs are regularly produced from different input carbon nanotubes (CNTs), prior work has not quantified how the CNT molecular aspect ratio r (i.e., molecular length-to-diameter ratio) influences k and in well-aligned, packed CNTFs. Here, we perform self-heated steady-state and three-omega thermal measurements at room temperature on CNTF suspended in vacuum. Our results show that k increases from 150 to 380 W / mK for viscosity-averaged molecular aspect ratios increasing from r = 960 to 5600 and nanotube diameters of ~2 nm, which we attribute to the effects of thermal resistances between CNT bundles. CNTFs made with varying volume fraction Φ of constituent high-r and low-r CNT have properties that fall within or below the typical macroscopic rule-of-mixtures bounds. The thermal diffusivity α scales with k, leading to a sample-averaged volumetric heat capacity of 1.5 ± 0.3 MJ / m 3 K⁠. This work's findings that fibers made from longer CNT have larger k and α at room temperature motivate further investigation into thermal transport in solution-spun CNTF.

42 ENGINEERING↗
DOE OSTI2025

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗
DOE OSTI2025

Atomic ionization: sd energy imbalance and Perdew–Zunger self-interaction correction energy penalty in 3d atoms

To accurately describe the energetics of transition metal systems, density functional approximations (DFAs) must provide a balanced description of s- and d- electrons. One measure of this is the sd transfer error, which has previously been defined as E ( 3 d n − 1 4 s 1 ) − E ( 3 d n − 2 4 s 2 ) . Theoretical concerns have been raised about this definition due to its evaluation of excited-state energies using ground-state DFAs. A more serious concern appears to be strong correlation in the 4s 2 configuration. Here, we define a ground-state measure of the sd energy imbalance, based on the errors of s- and d-electron second ionization energies of the 3d atoms, that effectively circumvents the aforementioned problems. We find an improved performance as we move from the local spin density approximation (LSDA) to the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) to the regularized and restored Strongly Constrained and Appropriately Normed (r 2 SCAN) meta-GGA for first-row transition metal atoms. However, we find large (∼2 eV) ground-state sd energy imbalances when applying a Perdew–Zunger 1981 self-interaction correction. This is attributed to an “energy penalty” associated with the noded 3d orbitals. A local scaling of the self-interaction correction to LSDA results in a balance of s- and d-errors.

Science & Technology - Other Topics↗
DOE OSTI2025

Reply to Smith and Siegel: Most lithium hops in paddlewheel-claimed conductors occur without spatially and temporally correlated anion-group rotations

We appreciate the engagement from Smith and Siegel on the topic of the “paddlewheel effect” and are pleased to see alignment regarding the absence of a paddlewheel effect, in which large-angle anion-group rotations directly propel lithium hops. In our paper, we clearly state the ambiguity surrounding the term “paddlewheel” to describe lithium transport. We believe that science is not well served by this vague term or its regular redefinition which is why we distinguish three types of anion-group rotation events: large-angle rotations (n-fold rotation returning to a rotationally invariant configuration), librations (e.g., rotational vibrations), and static changes of orientation in response to the change of Li occupancy. Here, we directly evaluate the spatial and temporal correlation between such polyanion rotations and Li hops.

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

Experimental and theoretical evidence of universality in superfluid vortex reconnections

The minimum separation between reconnecting vortices in fluids and superfluids obeys a universal scaling law with respect to time. The prereconnection and the postreconnection prefactors of this scaling law are different, a property related to irreversibility and to energy transfer and dissipation mechanisms. In the present work, we determine the temperature dependence of these prefactors in superfluid helium from experiments and a numeric model which fully accounts for the coupled dynamics of the superfluid vortex lines and the thermal normal fluid component. At all temperatures, we observe a pre- and postreconnection asymmetry similar to that observed in other superfluids and in classical viscous fluids, indicating that vortex reconnections display a universal behavior independent of the small-scale regularizing dynamics. We also numerically show that each vortex reconnection event represents a sudden injection of energy in the normal fluid. Finally we argue that in a turbulent flow, these punctuated energy injections can sustain the normal fluid in a perturbed state, provided that the density of superfluid vortices is large enough.

Reconnections↗
DOE OSTI2025

Personalized Tucker Decomposition: Modeling Commonality and Peculiarity on Tensor Data

In this paper, we propose a personalized Tucker decomposition (perTucker) to address the limitations of traditional tensor decomposition methods in capturing heterogeneity across different datasets. perTucker decomposes tensor data into shared global components and personalized local components. We introduce an order orthogonality assumption and develop a proximal gradient regularized block coordinate descent algorithm guaranteed to converge to a stationary point. The unique and common representations learned by perTucker reveal intrinsic statistical patterns in data and provide valuable information for a wide range of downstream analytics, including anomaly detection, source classification, and clustering. We demonstrate perTucker’s effectiveness through a simulation study and two case studies on solar flare detection and tonnage signal classification.

14 SOLAR ENERGY↗
DOE OSTI2025

Microfluidic cells for the 1–10 2 MPa pressure range

Thin membrane-delimited fluid cells supporting up to 1 at (0.1 MPa) of pressure are well known and commercially available for use in vacuum chambers of electron, photon, or various particle beam microscopies or spectroscopies. Hereby, we report on the development of fluid cells capable of working at 1–10 MPa, extending the analysis domain for investigating chemical, biochemical, or physical processes at pressures of interest in chemical synthesis, underwater biochemistry studies or underground geology. We explored ways to optimize cell membranes to better resist pressure beyond simply increasing the thickness or decreasing the size of the membranes, using finite element analysis and experimental validation via membrane bulging experiments and failure statistics. Fluid cell prototypes were fabricated using ∼75 nm-thick SiN x membranes, engineered to withstand 4.7 MPa (average value), compared to regular (un-engineered) membranes withstanding only 3.4 MPa (average value). The fluid cell prototypes include eight microchannels for feeding/evacuating the fluids and applying pressure into micro-reaction chambers, two electrodes for electrochemical or conduction measurements in the sample, and a possible pressure or temperature sensor, customizable for specific experiments.

hi-pressure↗
DOE OSTI2024

Pion valence quark distribution at physical pion mass of N f = 2 + 1 + 1 lattice QCD

Abstract We present a state-of-the-art calculation of the unpolarized pion valence-quark distribution in the framework of large-momentum effective theory (LaMET) with improved handling of systematic errors as well as two-loop perturbative matching. We use lattice ensembles generated by the MILC collaboration at lattice spacinga≈ 0.09 fm, lattice volume 64 3 × 96,N f = 2 + 1 + 1 flavors of highly-improved staggered quarks and a physical pion mass. The LaMET matrix elements are calculated with pions boosted to momentumP z ≈ 1.72 GeV with high-statistics ofO(10 6 ) measurements. We study the pion PDF in both hybrid-ratio and hybrid-regularization-independent momentum subtraction (hybrid-RI/MOM) schemes and also compare the systematic errors with and without the addition of leading-renormalon resummation (LRR) and renormalization-group resummation (RGR) in both the renormalization and lightcone matching. The final lightcone PDF results are presented in the modified minimal-subtraction scheme at renormalization scaleμ= 2.0 GeV. We show that thex-dependent PDFs are compatible between the hybrid-ratio and hybrid-RI/MOM renormalization with the same improvements. We also show that systematics are greatly reduced by the simultaneous inclusion of RGR and LRR and that these methods are necessary if improved precision is to be reached with higher-order terms in renormalization and matching.

Physics↗
DOE OSTI2024

One-shot omnidirectional pressure integration through matrix inversion

In this work, we present a method to perform 2D and 3D omnidirectional pressure integration from velocity measurements with a single-iteration matrix inversion approach. This work builds upon our previous work, where the rotating parallel ray approach was extended to the limit of infinite rays by taking continuous projection integrals of the ray paths and recasting the problem as an iterative matrix inversion problem. This iterative matrix equation is now 'fast-forwarded' to the 'infinity' iteration, leading to a different matrix equation that can be solved in a single step, thereby presenting the same computational complexity as the Poisson equation. We observe computational speedups of ~10 6 when compared to brute-force omnidirectional integration methods, enabling the treatment of grids of ~10 9 points and potentially even larger in a desktop setup at the time of publication. Further examination of the boundary conditions of our one-shot method shows that omnidirectional pressure integration implements a boundary condition where the boundary points are treated as interior points to the extent that information is available. Finally, we show how the method can be extended from the regular grids typical of particle image velocimetry to the unstructured meshes characteristic of particle tracking velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Overcoming time and complexity limitations in molecular dynamics investigations of equilibrium melting

Abstract A hybrid Monte-Carlo molecular-dynamics method for determining solidus and liquidus compositions in multicomponent systems is presented that overcomes both the time limitations in conventional molecular dynamics that prevent the evolution of distinct solid and liquid compositions via diffusion and the complexity challenge that prevents use of thermodynamic assessment in systems of many components. This hybrid method is validated in the Cu–Ni system against an independent assessment of solidus and liquidus compositions based on the regular solution model. Strategies for efficient mapping of different phase diagrams, based on the thermodynamic parameter T 0 , the temperature at which two phases of the composition X 0 have equal free energies, are presented and then demonstrated for the copper-nickel fully miscible system and the gold–silicon eutectic system. A calculation of the solidus and liquidus sampled during the equilibrium melting of equiatomic CrMnFeCoNi is performed, indicating that this method has potential to be extended to the study of many component alloys.

Au-Si↗
DOE OSTI2025

Coherent frequency combs from electrons colliding with a laser pulse

Highly coherent and powerful light sources capable of generating soft x-ray frequency combs are essential for high precision measurements and rigorous tests of fundamental physics. In this work, we derive the analytical conditions required for the emission of coherent radiation from an electron beam colliding with a laser pulse, modeled as a plane wave. These conditions are applied in a series of numerical simulations, where we show that a soft x-ray frequency comb can be produced if the electrons are regularly-spaced and sufficiently monoenergetic. High quality beams of this kind may be produced in the near future from laser-plasma interactions or linear accelerators. Furthermore, we highlight the advantageous role of employing few-cycle laser pulses in relaxing the stringent monoenergeticity requirements for coherent emission. The conditions derived here can also be used to optimize coherent emission in other frequency ranges, such as the terahertz domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MHD-induced SOL filaments and divertor heat flux striations in NSTX

Magnetohydrodynamic (MHD) activity in NSTX has been found to induce striated heat flux patterns at the divertor. These striated heat fluxes were previously interpreted as strike-point splitting caused by stochastic magnetic lobes near the X-point. Recent analysis indicates that MHD can directly generate filaments in the scrape-off layer (SOL), and the MHD-related filaments are consistent with the rotation and the frequency of the MHD. Magnetic field line tracing from the mid-plane of the MHD-related filaments to the divertor is consistent with the location of striated heat fluxes on the divertor. Additionally, gas puff imaging (GPI) diagnostic data shows the generation of regular blobs following each edge harmonic oscillation (EHO) toroidal mode number n = 1 perturbation, which indicates that the MHD perturbation directly drives the blobs. These blobs are responsible for the divertor heat flux striations. Here, this finding provides a new framework for understanding MHD-induced edge transport and has important implications for heat-load control in NSTX-U and future spherical tokamaks.

MHD↗
DOE OSTI2026

Access and limits of RMP ELM suppression with n = 1 fields in DIII-D

This work reports on DIII-D experiments aimed at extending resonant magnetic perturbation (RMP) suppression of edge localized modes (ELMs) to n = 1 fields, where n is the toroidal mode number. Modeling of the 3D ideal MHD plasma response to the RMPs using the GPEC code is used to quantify edge and core resonant fluxes, guiding experimental strategies to increase plasma resilience against core error field penetration, optimize multicoil phasing, and explore higher q 95 operation. In DIII-D, ELM mitigation is regularly observed across a wide range of n = 1 RMP scenarios. A ∼100 ms phase of complete ELM suppression was achieved at q 95 ∼ 3.9 using an odd-parity coil configuration. The suppressed phase exhibited clear signatures of RMP ELM suppression, including the elimination of Dα spikes, increased pedestal rotation, enhanced magnetic response, and elevated broadband density turbulence. An optimized coil configuration for edge-to-core resonant flux did show increased edge resonance indicated by increased density pumpout, but did not yield RMP ELM suppression. At q 95 ∼ 5.1, a bifurcation to a grassy-like ELM regime occurred, while large type-I ELMs persisted. These results demonstrate progress in experimental access to n = 1 RMP ELM suppression in DIII-D, motivating further study for robust access. This work also highlights the potential role of 3D edge stability as well as rational surface alignment in RMP ELM suppression access, which has important implications for the use of low-n RMPs in future reactor-scale devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dark photon constraints from CMB temperature anisotropies

The resonant conversion, within the inter-galactic medium, of regular photons into dark photons amplifies the anisotropy observed in the CMB, thereby imposing stringent constraints on the existence of light dark photons. In this study, we investigate the impact of light dark photons, with masses in the range 3 × 10$^{-15}$ eV < mA' < 3 × 10$^{-12}$ eV on the power spectrum of temperature anisotropies within the cosmic microwave background (CMB) radiation utilizing the state-of-the-art large-volume FLAMINGO cosmological simulations. Our results show that using full Planck data, one can expect the existing constraints on the dark photon mixing parameter in this mass range to improve by an order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Plasma rotation and diamagnetic drift effects on the resistive wall modes in the negative triangularity tokamaks

Abstract It was found previously that the negative triangularity (NT) configuration is more MHD-unstable for low n modes than the positive triangularity (PT) case, although the situation is reversed for intermediate n modes and the NT configuration becomes more stable for intermediate n modes ( n = 3 − 10 ) (Zheng et al 2021 Nucl. Fusion 61 116014). Here, n is the toroidal mode number. In this work, we extend the studies to include the rotation effects, as well as the diamagnetic drift effects, to see how the resistive wall modes (RWMs) in the NT configuration are affected as compared with the PT configuration. This is particularly motivated by noting that the wall interface with the plasma is quite different between the NT and PT configurations. It affects the plasma rotation and diamagnetic drift effects on the low n RWM. We consider the DIII-D-NT-experiment equilibrium reconstructed by the EFIT code. Based on the equilibrium g-file, the extended equilibria are constructed with the VMEC code by varying the beta values while keeping the pressure and poloidal current flux profiles basically unchanged. The bootstrap current contribution to the equilibria is taken into account with the Sauter formula. The MHD stability is then computed using the AEGIS code with the rotation and diamagnetic drift effects taken into account. We found that, although the NT configuration is less stable for n = 1 MHD modes, the rotation and diamagnetic drift stabilization effects on RWMs are more effective in the NT configuration than in the PT one. Note that even in the PT case, the stabilization of RWMs by the rotation and kinetic effects is critical. Because the low-n RWMs in the regular NT case are more unstable, the rotation and diamagnetic drift stabilization effects found in this research are important for the NT tokamak concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗