Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependable computing”

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

Quantum Computing for Energy-Related Applications

Growing interest in quantum computing and simulations have created opportunities for its deployment to improve processes pertaining to energy production, distribution, and consumption. While quantum computing is considered as a paradigm shift in our basic understanding of physical computation, effective implementation of quantum computing in energy applications also depends on progress and development in the dimensions of both quantum computing hardware and quantum computing algorithms. To fully address the status and future challenges of quantum information science (QIS) applied within the energy sector, in this presentation, we firstly summarize recent advancements on the applications of quantum computing to energy infrastructure and materials, complex energy system processes, advanced manufacturing, and energy system security. Then, we will demonstrate the results of quantum computing both on a simulator and a quantum device accessing from OLCF. Our first example is to use the variational quantum eigensolver (VQE) with a unitary coupled cluster with singles and doubles (UCCSD) ansatz to simulate a series of LixHyq molecules (q=-1, 0, +1). The obtained results showed that the quantum computing VQE-UCCSD is comparable to classical CCSD for small systems like LiH with respect to full configuration interaction (FCI). Targeting on CO2 capture application, our second example is to use VQE to quantify molecular vibrational energies and reaction pathways between CO2 and a simplified amine-based solvent model—NH3 to form H2NCOOH. This research showcases quantum computing applications in the study of CO2 capture reactions.

Duan, Yuhua↗
DOE OSTI2025

Parallel Programming in MCNP6

Monte Carlo N-Particle (MCNP)1 is a general-purpose Monte Carlo particle transport code developed by Los Alamos National Laboratory (LANL). To efficiently handle long simulations, MCNP version 6 (MCNP6) supports parallel execution using two primary programming models: • Shared-memory task-based threading using OpenMP (Open Multi-Processing), and • Distributed-memory calculations using MPI (Message Passing Interface). The OpenMP and MPI programming models enable MCNP6 to scale from desktop systems to high-performance computing (HPC) clusters, allowing users to run MCNP in one of three parallel modes: • OpenMP-only, • MPI-only, and • Hybrid (MPI + OpenMP). The choice of parallelization mode depends on the underlying computer architecture and the characteristics of the simulation problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The Extended Embedded Self-Shielding Method in SCALE 6.3/Polaris

The SCALE transport lattice code, Polaris, has been previously developed to generate few-group homogenized cross sections for whole-core nodal diffusion simulators in which the embedded self-shielding method (ESSM) is used for resonance self-shielding calculations to process cross sections. Although the ESSM capability has been very successful in light-water reactor analysis, it may require enhancements in computational efficiency; treatment of spatially dependent resonance self-shielding effects; and handling of interrelated resonance effects among fuel, cladding, and control rod materials. Therefore, this study focuses on improving computational efficiency by using a Dancoff-based Wigner–Seitz approximation combined with a material-based resonance categorization, through which a spatially dependent ESSM capability is developed to accurately estimate self-shielded cross sections inside the fuel. Benchmark results show that the new capability significantly enhances computational efficiency and accuracy for spatially dependent local zones within the fuel and through depletion.

ESSM↗
DOE OSTI2026

Solidification cracking of refractory alloys: a computational and machine learning study to investigate composition-dependence for improved weldability and additive manufacturability

Large-batch numerical, CALculation of PHAse Diagrams (CALPHAD)-based solidification cracking calculations are performed and then analyzed with machine learning methods to generate models that relate chemistry of refractory alloys to cracking susceptibility. Kou’s solidification cracking index is used to study the refractory alloys including O, N, C binary mixtures with Mo, Ta, Nb, and W, the molybdenum-based TZM, Niobium-based C103, and Tantalum-based T111 and Ta-10 W, as well as hypothetical refractory ternary alloys. Findings strongly validate Kou’s Crack Susceptibility Index (CSI) against Varestraint test data for Nb- and Ta-based alloys, establishing CSI thresholds where refractory alloys with CSI < 15,000 K are likely weldable, CSI > 15,000 K are prone to cracking, and CSI > 25,000 K are likely unweldable (or unprintable). Furthermore, interstitial elements C, N, and O significantly increase crack susceptibility, with some existing material specifications coinciding with peak cracking susceptibility concentrations. Finally, machine learning-derived elemental potency factors enable rapid prediction of CSI from alloy chemistry for C103, TZM, Ta-10 W, and T-111 alloys. These results provide practical guidance for feedstock selection, powder reuse limits, and alloy specification amendments for welding and additive manufacturing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Size-Dependent Optical Band Gaps in Metal–Organic Framework Nanoparticles

Decades of research into size-dependent semiconductor optical gaps have focused on quantum confinement as the dominant mechanism. Emerging reports indicate that lattice strain─intentional or incidental─can impart optical shifts similar or greater in magnitude. Here, we report evidence of optical absorption and photoluminescence spectra of M(1,2,3-triazolate)2 (M = Mg, Cr, Mn, Fe, Co, Cu, Zn, or Cd) nanoparticles that blueshift from bulk values with decreasing particle sizes in a manner that defies explanation by conventional quantum confinement. Here, the phenomenon persists for particle sizes as large as 200 nm, whereas quantum confinement generally ceases beyond 20–30 nm diameters and follows a weaker dependence on the particle radius. Computational simulations and crystallographic analysis suggest that this behavior arises from size-dependent changes to metal–linker bonding that manifest in strain values comparable to literature reports of strain-induced optical shifts in other classes of materials. This behavior appears beyond this family of materials in other notable examples of metal–organic frameworks (MOFs), including the well-studied Cu3(trimesate)2 (CuBTC), where smaller sizes correlate with blueshifted optical gaps. Taken together, these results represent one of the few examples of size-dependent strain in crystalline materials and reinforce the emerging view that MOFs become softer materials when isolated as nanoparticles.

Electrical conductivity↗
DOE OSTI2024

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.

Science & Technology - Other Topics↗
DOE OSTI2026

A Benchmark Suite for Evaluating Scientific AI Workloads on GPUs

AI applications have been steadily increasing in the allocation portfolio among leadership computing facilities. These applications depend on deep learning frameworks with hardware acceleration and underlying software systems. With the rapid development of applications, software stacks, and hardware devices, it is essential to evaluate the performance of core operations in AI workloads for direction of optimizations and procurement of next-generation high-performance computing (HPC) infrastructures. Currently, most benchmarks lack scientific AI workloads. So, we present DeepKernelBench and the experimental results of evaluating the benchmark suite for early observations and performance comparisons on datacenter GPUs using representative workloads for scientific AI, including Attentions, General matrix multiplications, Geometrics and Fourier neural operations.

Jin, Zheming [Advanced Micro Devices (AMD)]↗
DOE OSTI2025

Shape Anisotropy-Dependent Leaking in Magnetic Neurons for Bio-Mimetic Neuromorphic Computing

Spiking neural networks seek to emulate biological computation through interconnected artificial neuron and synapse devices. Spintronic neurons can leverage magnetization physics to mimic biological neuron functions, such as integration tied to magnetic domain wall (DW) propagation in a patterned nanotrack and firing tied to the resistance change of a magnetic tunnel junction (MTJ), captured in the domain wall-magnetic tunnel junction (DW-MTJ) device. Leaking, relaxation of a neuron when it is not under stimulation, is also predicted to be implemented based on DW drift as a DW relaxes to a low energy position, but it has not been well explored or demonstrated in device prototypes. Here, in this work, we study DW-MTJ artificial neurons capable of leaky integrate-and-fire (LIF) behavior and demonstrate geometry-dependent leaking dynamics that results in repeatable, tunable LIF operation. Studying the behavior of five different device designs, we show tuning the geometry, stimulating fields and currents, and location of electrical contacts results in a wide range of neuron behavior. Additionally, implementation of an asymmetric notch allows for nonlinear pinning which increased expressivity without sacrificing leaking. The measured behavior is implemented in a simulated spiking neural network that outperforms a 1D model of continuous DW motion and approaches the performance of an ideal LIF activation function. The results show that the analog LIF capability of DW-MTJ neurons combines many desirable neuron functions into a single device, which can result in varied forms of multifunctional neuromorphic computing.

42 ENGINEERING↗
DOE OSTI2025

Shaping dark photon spectral distortions

The cosmic microwave background (CMB) spectrum is an extraordinary tool for exploring physics beyond the Standard Model. The exquisite precision of its measurement makes it particularly sensitive to small effects caused by hidden sector interactions. In particular, CMB spectral distortions can unveil the existence of dark photons which are kinetically coupled to the standard photon. In this work, we use the COBE-FIRAS dataset to derive accurate and robust limits on photon-to-dark-photon oscillations for a large range of dark photon masses, from 10$^{−10}$ to 10$^{−4}$ eV. We consider in detail the redshift dependence of the bounds, computing CMB distortions due to photon injection/removal using a Green’s function method. Our treatment improves on previous results, which had set limits studying energy injection/removal into baryons rather than photon injection/removal, or ignoring the redshift evolution of distortions. The difference between our treatment and previous ones is particularly noticeable in the predicted spectral shape of the distortions, a smoking gun signature for photon-to-dark-photon oscillations. The characterization of the spectral shape is crucial for future CMB missions, which could improve the present sensitivity by orders of magnitude, exploring regions of the dark photon parameter space that are otherwise difficult to access [inline-graphic not available: see fulltext].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗
DOE OSTI2026

Sub-Doppler Cooling of a Trapped Ion in a Phase-Stable Polarization Gradient

Trapped ions provide a highly controlled platform for quantum sensors, clocks, simulators, and computers, all of which depend on cooling ions close to their motional ground state. Existing methods like Doppler, resolved sideband, and dark resonance cooling balance trade-offs between the final temperature and cooling rate. A traveling polarization gradient has been shown to cool multiple modes quickly and in parallel, but utilizing a stable polarization gradient can achieve lower ion energies, while also allowing tailorable light-matter interactions in the sub-wavelength regime. In this Letter, we demonstrate cooling of a trapped ion below the Doppler limit using a phase-stable polarization gradient created using trap-integrated photonic devices. At an axial frequency of 2⁢𝜋×1.45 MHz we achieve ⟨𝑛⟩=1.56±0.07 in 150 μ⁢s and cooling rates of ∼0.3 quanta/μ⁢s. Here, we examine ion dynamics under different polarization gradient phases, detunings, and intensities, showing reasonable agreement between experimental results and a multilevel model. Cooling is fast and power efficient, with lower average motional Fock state occupation when compared to simulated operation under the corresponding running wave configuration. Our results demonstrate a well-controlled test bed for studying the dynamics of multilevel atomic systems in a phase-stable polarization gradient.

Atom & ion cooling↗
DOE OSTI2024

Examining the possibility that normal nuclear matter is quarkyonic

The possibility that nuclear matter might be quarkyonic is considered. Quarkyonic matter is high baryon density matter that is confined but can be approximately thought of as a filled Fermi sea of quarks surrounded by a shell of nucleons. Here, nuclear matter is described by the IdylliQ sigma model for quarkyonic matter, generalizing the noninteracting IdylliQ model [Y. Fujimoto et al., Phys. Rev. Lett. 132, 112701 (2024)] to include interactions with a σ meson and a pion. When such interactions are included, we find that isospin-symmetric nuclear matter binds with acceptable values of the compressibility and other parameters for nuclear matter at saturation. The energy per nucleon and sound velocity of such matter is computed, and the isospin dependence is determined. Nuclear matter is formed at a density close to but slightly above the density at which quarkyonic matter forms. Quarkyonic matter predicts a strong depletion of nucleons in normal nuclear matter at low momentum. Finally, such a depletion for nucleon momenta k ≲ 120 MeV is shown to be consistent with electron scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Open-source library for performance-portable neutrino reaction rates: Application to neutron star mergers

A realistic and detailed description of neutrinos in binary neutron star (BNS) mergers is essential to build reliable models of such systems. To this end, we present bns_nurates, a novel open-source numerical library designed for the efficient on-the-fly computation of neutrino interactions, with particular focus on regimes relevant to BNS mergers. bns_nurates targets a higher level of accuracy and realism in the implementation of commonly employed reactions by accounting for relevant microphysics effects on the interactions, such as weak magnetism and mean field effects. It also includes the contributions of inelastic neutrino scattering off electrons and positrons and (inverse) nucleon decays. Finally, it offers a way to reconstruct the neutrino distribution function in the framework of moment-based transport schemes. As a first application, we compute both energy-dependent and energy-integrated neutrino emissivities and opacities for conditions extracted from a BNS merger simulation with m1 transport scheme. We find some qualitative differences in the results when considering the impact of the additional relevant reactions and of microphysics effects. For example, neutrino-electron/positron scattering reactions are important for the energy exchange of heavy-type neutrinos as they do not undergo semileptonic charged-current processes, when μ± are not accounted for. Moreover, weak magnetism and mean field effects can significantly modify the contribution of β processes for electron-type (anti)neutrinos, increasing at the same time the importance of (inverse) neutron decays. Here, the improved treatment for the reaction rates also modifies the conditions at which neutrinos decouple from matter in the system, potentially affecting their emission spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Sign Problem in Tensor-Network Contraction

We investigate how the computational difficulty of contracting tensor networks depends on the sign structure of the tensor entries. Using results from computational complexity, we observe that the approximate contraction of tensor networks with only positive entries has lower computational complexity as compared to tensor networks with general real or complex entries. This raises the question of how this transition in computational complexity manifests itself in the hardness of different tensor-network-contraction schemes. We pursue this question by studying random tensor networks with varying bias toward positive entries. First, we consider contraction via Monte Carlo sampling and find that the transition from hard to easy occurs when the tensor entries become predominantly positive; this can be understood as a tensor-network manifestation of the well-known negative-sign problem in quantum Monte Carlo. Second, we analyze the commonly used contraction based on boundary tensor networks. The performance of this scheme is governed by the number of correlations in contiguous parts of the tensor network (which by analogy can be thought of as entanglement). Remarkably, we find that the transition from hard to easy—i.e., from a volume-law to a boundary-law scaling of entanglement—already occurs for a slight bias of the tensor entries toward a positive mean, scaling inversely with the bond dimension D , and thus the problem becomes easy the earlier the larger D occurs. This is in contrast both to expectations and to the behavior found in Monte Carlo contraction, where the hardness at fixed bias increases with the bond dimension. To provide insight into this early breakdown of computational hardness and the accompanying entanglement transition, we construct an effective classical statistical-mechanical model that predicts a transition at a bias of the tensor entries of 1 / D , confirming our observations. We conclude by investigating the computational difficulty of computing expectation values of tensor-network wave functions (projected entangled-pair states, PEPSs) and find that in this setting, the complexity of entanglement-based contraction always remains low. We explain this by providing a local transformation that maps PEPS expectation values to a positive-valued tensor network. This not only provides insight into the origin of the observed boundary-law entanglement scaling but also suggests new approaches toward PEPS contraction based on positive decompositions. Published by the American Physical Society 2025

Chen, Jielun (ORCID:0000000178411545)↗
DOE OSTI2026

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗
DOE OSTI2025

PyPTkinetics

This basic research code computes the time-dependent phase transition between two solid crystalline phases under dynamic conditions.

Blaschke, Daniel N.↗
DOE OSTI2026

braggedgemodeling

braggedgemodeling (bem) is an open-source Python package for modeling neutron Bragg-edge imaging. It computes the wavelength-dependent total neutron cross-section of a material (coherent and incoherent elastic, coherent and incoherent inelastic scattering, and absorption) from its crystal structure, and implements the March-Dollase texture model and the Jorgensen peak profile, supporting quantitative analysis of energy-resolved neutron imaging data (phase, stress/strain, and texture). Published in the Journal of Open Source Software (2018).

Lin, Jiao [Oak Ridge National Laboratory (ORNL), O↗
DOE OSTI2026

Dielectric Resonator Design for Low Power and Low Temperature Microwave Plasma

Waveguide-based microwave plasmas generally operate at high temperatures (2000 - 6000K)[1], making it difficult to directly interface solid materials with the plasma without significant thermal damage. Dielectric microwave resonators (DMRs), long studied for wave-based manipulation of electromagnetic radiation for telecom and optics, can focus radiation to extremely small mode volumes, creating intense localized fields with low-power input.[2] This phenomenon can be used for applications ranging from efficient plasma electronics to near-ambient plasma-materials interactions. Such DMR-based plasmas have been demonstrated a handful of times in the literature, but the majority of research towards this utilize the lowest frequency resonance mode.[3], [4], [5] By carefully controlling the geometry of cylindrical resonators, a variety of electromagnetic modes can be excited. In this work, COMSOL Multiphysics simulations are used to study the electric field enhancement and absorption properties of CaTiO3 DMRs as a function of geometry and excitation frequency. Whereas previous studies have utilized the HEM111 resonance frequency to drive low power plasma excitation, we find that higher order resonance frequencies are more effective at field enhancement and result in less power loss within the dielectric material, hence less wasted heating. The effectiveness of these modes is also geometry dependent and can be computationally optimized for plasma generation. Complementing these computational efforts, we demonstrate a new closed-system reactor design built in a WR-650 waveguide and experimentally demonstrate the formation of atmospheric argon microwave plasma using < 30 W input power on DMR dimers. We observe a shifting resonance frequency as the DMRs heat in response to microwave excitation and develop a Python-based lock-in mechanism to effectively track the DMR resonance over time, leading to stable plasma operation. We use infrared thermal imaging to monitor the temperature of the DMR dimers and surrounding quartz chamber, demonstrating thermal temperatures < 60 degreesC. Finally, we utilize optical emission spectroscopy (OES) to probe the plasma properties as a function of the resonance mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗