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

ALD with Alternative Coreactants: Which Work, Which Do Not, and Why

In this paper, we present a combined experimental and theoretical study that systematically evaluates a series of alcohols (i.e., primary, secondary, and tertiary) as coreactants for trimethylaluminum (TMA)-based atomic layer deposition (ALD) of Al 2 O 3 . We employed in situ quartz crystal microbalance techniques and ex situ X-ray photoelectron spectroscopy to probe growth across a range of temperatures from T = 120 to 285 °C. Dispersion-inclusive hybrid density functional theory was employed to identify potential reaction pathways and compute corresponding Gibbs free-energy barriers and rate coefficients. Experimentally on an Al 2 O 3 surface, sustained thin-film growth, where the rate is constant over many cycles, was only observed with tertiary alcohols (tert-butanol and 2-methyl-2-butanol) at elevated temperatures (e.g., T = 285 °C); for tert-butanol specifically, sustained growth also occurred at T = 240 °C, with no sustained growth at or below T = 210 °C. Growth using primary and secondary alcohols decayed to a negligible value after only 2–3 cycles on an Al 2 O 3 surface at all tested temperatures. Intentional trace additions of H 2 O to anhydrous alcohols restored sustained growth but at reduced deposition rates relative to pure H 2 O. Theoretical analysis supports a bimolecular alkoxy β-H elimination mechanism as the primary reaction pathway: initiated by the formation of a bound alkoxy species, the subsequent β-H elimination step is rate-determining, with a barrier that systematically decreases as the degree of the alcohol coreactant is increased from primary to secondary to tertiary. Computed rate coefficients (k) also indicate a strong temperature dependence (with k increasing by ∼4–5 orders of magnitude from T = 120 to 285 °C) and are consistent with the experimental observation that only tertiary alcohols lead to steady thin-film growth at elevated temperatures. This integrated approach establishes why only certain alcohols can sustain Al 2 O 3 ALD and delivers a predictive framework for identifying effective alcohol coreactants.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beyond the two-point correlation: Constraining primordial non-Gaussianity with density-perturbation moments

Constraining primordial non-Gaussianity (PNG) on the large-scale cosmic structure (LSS) is an important step in understanding properties of the early Universe, specifically in distinguishing between different inflationary models. Measuring PNG relies on evaluating the scale-dependent correlations in the density field. New summary statistics beyond the two- and three-point correlation functions in configuration space and their Fourier-space counterparts, the power- and bispectrum may provide increased sensitivity. We introduce a new method for extracting the PNG signal imprinted on the LSS by using the first three Gaussian moments of the normalized correlation in density perturbations, evaluated on varying distance scales. We aim to assess this method’s sensitivity to local PNG, parameterized by f NL . We performed spherical convolutions on a range of scales on dark-matter-halo simulations to measure the scale-dependent correlations in the density field. From these, we computed the first three moments and compared them to a model expectation vector, parameterized to the second power in f NL . Our method provides about 21% improvement in sensitivity to f NL with respect to using the two-point correlation function alone. Notably, we find that the second moment alone carries nearly as much constraining power as the mean, highlighting the potential of higher order statistics. Given its simplicity and efficiency, this framework is well suited for application to current and upcoming large-scale surveys such as the Dark Energy Spectroscopic Instrument (DESI).

early universe↗
DOE OSTI2024

Identification of large polarons and exciton polarons in rutile and anatase polymorphs of titanium dioxide

Titanium dioxide (TiO 2 ) is a wide-gap semiconductor with numerous applications in photocatalysis, photovoltaics, and neuromorphic computing. The unique functional properties of this material critically depend on its ability to transport charge in the form of polarons, namely narrow electron wavepackets accompanied by local distortions of the crystal lattice. It is currently well established that the most important polymorphs of TiO 2 , the rutile and anatase phases, harbor small electron polarons and small hole polarons, respectively. However, whether additional polaronic species exist in TiO 2 , and under which conditions, remain open questions. Here, we provide definitive answers to these questions by exploring the rich landscape of polaron quasiparticles in TiO 2 via recently developed ab initio techniques. In addition to the already known small polarons, we identify three species, namely a large hole polaron in rutile, a large quasi-two-dimensional electron polaron in anatase, and a large exciton polaron in anatase. These findings complete the puzzle on the polaron physics of TiO 2 and pave the way for systematically probing and manipulating polarons in a broad class of complex oxides and quantum materials.

Dai, Zhenbang↗
DOE OSTI2025

Firewalls from General Covariance

I define “horizon normalcy” as the approximate validity of semiclassical gravity and effective field theory for the description of observers that approach or cross a black hole horizon. If black holes return information, then horizon normalcy must fail substantially, at least in some global states. It has been proposed that horizon normalcy persists, so long as the Hawking radiation remains in a computationally simple state. Here I argue that state-dependent horizon normalcy—independent of the underlying mechanism and independent of the class of radiation states asserted to guarantee normalcy—requires a breakdown of general covariance far from the black hole, or else horizon normalcy will depend on the infinite future of the exterior. This is because the radiation can be in different states at different events, all spacelike to the horizon crossing event whose normalcy is at stake. I discuss a related effect in AdS/CFT, and I argue that its resolution by timefolds is of no help here.

General relativity↗
DOE OSTI2026

Rolling Root Mean Square Based Multimodal Anomaly Detection for Real Time Monitoring of Smart Grid

Reliable real-time monitoring is valuable for maintaining the operational integrity of modern electrical smart grids. Deployment of heterogeneous sensing technologies in substations has enabled high-resolution, multichannel waveform monitoring, but also introduces challenges for anomaly detection due to noise, baseline drift, and modality-dependent signal characteristics. In this work, we present a computationally efficient unsupervised method for multimodal event detection based on Rolling Root Mean Square based Event Detection (RRMSED). The method is developed using in-house, field deployed sensors collecting data at a utility substation. The sensing system comprises voltage and current sensors, triaxial accelerometers, and magnetometers, collectively capturing electrical, vibrational, and magnetic waveform measurements at high temporal resolution. RRMSED operates by extracting rolling RMS energy features and their first-order temporal differences from consecutive waveform segments for each channel and then applying channel-specific statistical thresholds learned from historical data. A persistence-based exceedance logic is employed to robustly identify transient events while suppressing impulsive noise, and to provide precise temporal localization with high resolution. The framework is designed for continuous server-side operation and can be deployed in real time without requiring complex models. Experiments on simulated waveform data with known ground truth demonstrate low false positive (FP) and false negative (FN) rates. Application to real substation data shows RRMSED to identify events that are not captured by conventional monitoring indicators including fast transient detection algorithm currently deployed in the system. These results indicate that rolling RMS based features provide an effective and practical basis for real-time multimodal event detection in smart-grid substations.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗
DOE OSTI2025

efield_from_tdJ

The efield_from_tdJ.py Python tool is designed to compute the electric field from a general time-dependent current density (J) input.

Blaine, William [LANL]↗
DOE OSTI2026

LBNL_DIATOMICS v1

LBNL_DIATOMICS is a software package developed by P. Zarzycki (Lawrence Berkeley National Laboratory) for computing isotopic effects in diatomic molecules using mass-dependent spectroscopic parameters obtained experimentally or from first principles. The package provides a streamlined pipeline from the potential energy surface to the equilibrium constant of an isotopic exchange reaction.

Zarzycki, Piotr [Lawrence Berkeley National Labora↗
DOE OSTI2024

Quantum Computing and Simulations for Energy Applications

While quantum computing (QC) is considered as a paradigm shift in our basic understanding of physical computation, effective implementation of QC in energy applications also depends on progress and development in the dimensions of both QC hardware and algorithms. To fully address the status and future challenges of QC applied within the energy sector, in this presentation, we firstly summarize recent advancements on the applications of QC to energy infrastructure and materials, complex energy system processes, advanced manufacturing, and energy system security. Then, we will demonstrate the results of QC performed both on a simulator and a quantum device targeting on energy-related applications.

Paudel, Hari P.↗
DOE OSTI2025

Micro-structural features and material properties impact on adhesive metal joints via computational modeling and machine learning

The quality of structural bonding in practical applications depends on various factors arising from materials, pre-processing conditions, and manufacturing. Understanding how these factors influence bonding performance and determining their relative importance are of significant interest. Thus, this study evaluates the effects of microstructural features and material properties on the structural strength of adhesively-bonded metal joints at the submillimeter scale, utilizing a combination of Finite Element Modeling (FEM) and Machine Learning (ML) with Gradient Boosting Regression (GBR). The microstructural features include adhesive thickness, internal voids within the adhesive, adherend-adhesive interfacial voids, void size and volume fraction, and surface roughness. The material properties include the constitutive behavior of the adhesive, as well as the adherend-adhesive interfacial strength and fracture energy. The changes in structural strength and morphologies of the bonded metal structures with respect to different microstructural features and material properties were clarified by FEM. By further leveraging ML-GBR, the sequence of importance of these factors affecting bonding performance across various scenarios was summarized. This work provides valuable insights into the development of improved structural bonding for adhesive joints in industries such as automotive , aerospace, and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Water effects on NH 3 -SCR over Cu-based small-pore zeolite catalysts: A review

The selective catalytic reduction (SCR) of NO x with NH 3 as the reductant over small-pore Cu-zeolite catalysts is a leading strategy for controlling emissions from diesel engines and mobile sources. While the hydrothermal stability of Cu-zeolites such as Cu-SSZ-13, Cu-SAPO-34, and Cu-SSZ-39 has been extensively studied, the multifaceted role of water vapor under operating conditions has only recently received systematic attention. Water impacts multiple aspects of SCR catalysis, including NO x conversion, N 2 O formation, Cu ion speciation, redox kinetics, and surface reactions. This review comprehensively summarizes recent advances in understanding the effect of water on NH 3 -SCR over Cu-based catalysts. For Cu-SSZ-13, water vapor induces both inhibitory and promotive effects depending on temperature, Cu loading, and Si/Al ratio by altering NH 3 adsorption, Cu mobility, and redox half-cycle kinetics. Studies have revealed that water can suppress ammonium nitrate formation at low temperatures and enhance nitrate pathways at higher temperatures. Over Cu-SAPO-34, water promotes NO conversion by improving Cu(II) reducibility and strengthening Brønsted acidity. First-principles calculations confirm that water-coordinated Cu species lower O 2 activation barriers. In Cu-SSZ-39, water was found to facilitate Cu reduction, migration, and O 2 activation, resulting in superior water tolerance and unexpected catalytic promotion. Both experimental and computational studies highlight that water vapor exerts complex, catalyst-dependent effects on NH 3 -SCR performance, underscoring the need for tailored catalyst designs that leverage beneficial water interactions while minimizing inhibitory impacts.

Ammonia selective catalytic reduction↗
DOE OSTI2026

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗
DOE OSTI2025

Genetic programming for the nuclear many-body problem: a guide

Genetic Programming (GP) is an evolutionary algorithm that generates computer programs, or mathematical expressions, to solve complex problems. In this Guide, we demonstrate how to use GP to develop surrogate models to mitigate the computational costs of modeling atomic nuclei with ever increasing complexity. The computational burden escalates when uncertainty quantification is pursued, or when observables must be globally computed for thousands of nuclei. By studying three models in which the mean field depends on the total particle density self-consistently, we show that by constructing reduced order models supported by GP one can speed up many-body computations by several orders of magnitude with a negligible loss in accuracy.

dimensionality reduction↗
DOE OSTI2025

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗
DOE OSTI2025

A Modified Model Dielectric Function for Analyzing Optical Spectra of InGaN Nanofilms on Sapphire Substrates

Due to a lower InN bandgap energy Eg~0.7 eV, InxGa1−xN/Sapphire epifilms are considered valuable in the development of low-dimensional heterostructure-based photonic devices. Adjusting the composition x and thickness d in epitaxially grown films has offered many possibilities of light emission across a wide spectral range, from ultraviolet through visible into near-infrared regions. Optical properties have played important roles in making semiconductor materials useful in electro-optic applications. Despite the efforts to grow InxGa1−xN/Sapphire samples, no x- and d-dependent optical studies exist for ultrathin films. Many researchers have used computationally intensive methods to study the electronic band structures Ejk→, and subsequently derive optical properties. By including inter-band transitions at critical points from Ejk→, we have developed a semiempirical approach to comprehend the optical characteristics of InN, GaN and InxGa1−xN. Refractive indices of InxGa1−xN and sapphire substrate are meticulously integrated into a transfer matrix method to simulate d- and x-dependent reflectivity RE and transmission TE spectra of nanostructured InxGa1−xN/Sapphire epifilms. Analyses of RE and TE have offered accurate x-dependent shifts of energy gaps for InxGa1−xN (x = 0.5, 0.7) in excellent agreement with the experimental data.

Talwar, Devki N. (ORCID:0000000202723195)↗
DOE OSTI2024

The Transient Multi-Level method for Monte Carlo reactor statics calculations

The Transient Multi-Level (TML) method is applied to a time-dependent Monte Carlo transport solver to offload some of the computational burden of the expensive Monte Carlo solve to lower-order Coarse Mesh Finite Difference (CMFD) and Exact Point Kinetics Equations (EPKE) solvers via factorization of the neutron flux at the transport and CMFD levels using the Predictor Corrector Quasi-Static Method (PCQM). The Monte Carlo transient is solved by a modified fission source iteration scheme that introduces a single transient source bank. The method is implemented in the production-level Monte Carlo code, Shift, and verified with prescribed reactivity ramps from the two-dimensional version of the C5G7-TD reactor benchmark. The results show that, as compared to other quasi-static methods, the TML reduces the stochastic noise inherent to the transient Monte Carlo solver by factors of ~2 to 6 for various norm comparisons of the reactor power amplitude. Finally, the TML additionally reduces the number of Monte Carlo evaluations needed to simulate the transient, leading to roughly an order of magnitude improvement in CPU time relative to the standard PCQM for the problems tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

First principles investigation of dopants and defect complexes in CdSe$_x$Te$_{1-x}$

Se alloying is a common approach to improve the performance of CdTe solar cells by tuning the bandgap, defect levels, and carrier density. A fundamental understanding of these improvements, specifically the effect of Se alloying on the behavior of defects and dopants in CdTe, remains unclear. Here, in this work, we present a density functional theory (DFT) study of point defect energetics in CdTe and CdSe x Te 1-x with x = 0.25, leading to a comparison of how native defects, dopants (As and Cu), impurities (Cl and O), and related defect complexes behave in CdTe vs CdSe x Te 1-x . Our calculations, performed by combining semi-local and nonlocal hybrid functionals, show a general lowering of the formation energies of native defects as well as substitutional defects formed by As and Cl upon Se addition. For successful p-type doping with As, destabilizing Cl-based defects in the CdSeTe lattice would be essential. We find evidence for some low-energy defect complexes of As, Cl, and O in CdSe 0.25 Te 0.75 . The computed defect formation energies further enable estimates of temperature-dependent defect concentrations and self-consistent Fermi levels. A comparison of defect energetics with the energies of impurity phases reveals that As, Cu, Cl, and O overwhelmingly prefer being segregated to unwanted As 2 O 5 , AsCl 3 , Cd 2 AsCl 2 , and CuO x phases rather than remain at defect sites, but such segregation is less likely to happen in CdSe 0.25 Te 0.75 than in CdTe. Overall, our work presents a list of likely defects and complexes in CdTe and Se-incorporated CdTe, paving the way to explain and mitigate limited dopant activation in experimental observations.

CdTe↗
DOE OSTI2024

Molecular Modes Elucidate the Nuclear Magnetic Resonance Relaxation of Viscous Fluids

The Bloembergen, Purcell, and Pound (BPP) theory of nuclear magnetic resonance (NMR) relaxation in fluids dating back to 1948 continues to be the linchpin in interpreting NMR relaxation data in applications ranging from characterizing fluids in porous media to medical imaging (MRI). The BPP theory is founded on assuming molecules are hard spheres with 1 H– 1 H dipole pairs reorienting randomly; assumptions that are severe in light of modern understanding of liquids. Nevertheless, it is intriguing to this day that the BPP theory was consistent with the original experimental data for glycerol, a hydrogen-bonding molecular fluid for which the hard-sphere-rigid-dipole assumption is inapplicable. To better understand this incongruity, atomistic molecular simulations are used to compute 1 H NMR T 1 relaxation dispersion (i.e., frequency dependence) in two contrasting cases: glycerol, and a (non hydrogen-bonding) viscosity standard. At high viscosities, simulations predict distinct functional forms of T 1 for glycerol compared to the viscosity standard, in agreement with modern measurements, yet both in contrast to BPP theory. The cause of these departures from BPP theory is elucidated, without assuming any relaxation models and without any free parameters, by decomposing the simulated T 1 response into dynamic molecular modes for both intramolecular and intermolecular interactions. The decomposition into dynamic molecular modes provides an alternative framework to understand the physics of NMR relaxation for viscous fluids.

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

Convolutional Non-Homogeneous Poisson Process and its Application to Wildfire Ignition Risk Quantification for Power Delivery Networks

To quantify wildfire ignition risks on power delivery networks, the current practice predominantly relies on the empirically calculated fire danger indices, which may not well capture the effects of dynamically changing environmental factors. This article proposes a spatio-temporal point process model, known as the Convolutional Non-homogeneous Poisson Process (cNHPP), and applies the model to quantify wildfire ignition risks for power delivery networks. The proposed model captures both the current (i.e., instantaneous) and cumulative (i.e., historical) effects of key environmental processes (i.e., covariates) on wildfire risks, as well as the spatio-temporal dependency among different segments of the power delivery network. The computation and interpretation of the intensity function are thoroughly investigated. We apply the proposed approach to estimate wildfire ignition risks on major transmission lines in California, using historical fire data, meteorological and vegetation data obtained from the National Oceanic and Atmospheric Administration and National Aeronautics and Space Administration. Here, a comprehensive comparison study is performed to show the applicability and predictive capability of the proposed approach.

Non-homogeneous Poisson Process↗