Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theory and computation”

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

Machine learning guided prediction of solute segregation at coherent and semi-coherent metal/oxide interfaces

Investigation of semi-coherent metal/oxide interfaces with misfit dislocations using density functional theory (DFT) is computationally intensive to the point of being prohibitive, as it involves several hundreds to many thousands of atoms. In this study, we examined the solute segregation behavior at the Fe/Y 2 O 3 interface—a model interface for cladding applications in nuclear fission reactors—using a combination of DFT calculations and machine learning (ML) approaches. Both coherent and semi-coherent interfaces were considered. ML models were trained on DFT-calculated segregation energies to identify the key chemical, geometric and strain energy related features that govern solute segregation behavior at coherent Fe/Y 2 O 3 interfaces. Furthermore, it was found that ML models when trained on DFT calculated segregation energy of elements at a coherent interface, comprising of about a hundred-atom supercell, can predict the segregation energy of elements at a semi-coherent Fe/Y 2 O 3 interface (with multiple hundreds of atoms) at a fraction of computational cost (1/35th), with an accuracy comparable to DFT calculations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A terminology for scientific workflow systems

The term “scientific workflow” has evolved over the last two decades to encompass a broad range of compositions of interdependent compute tasks and data movements. It has also become an umbrella term for processing in modern scientific applications. Today, many scientific applications can be considered as workflows made of multiple dependent steps, and hundreds of workflow systems have been developed to manage and run these scientific workflows. However, no turnkey solution has emerged from the field to address the diversity of scientific processes and the infrastructure on which they are supposed to be implemented. Instead, new research problems requiring the execution of scientific workflows with some novel feature often lead to the development of an entirely new workflow system. A direct consequence of this situation is that many existing workflow management systems (WMSs) share some salient features, offer similar functionalities, and can manage the same categories of workflows but at the same time also have some distinct capabilities that can be important for specific applications. This situation makes researchers who develop workflows face the complex question of selecting a WMS. This selection can be driven by technical considerations, to find the system that is the most appropriate for their application and for the computing and storage resources available to them, or other factors such as reputation, adoption, strong community support, or long-term sustainability. To address this problem, a group of WMS developers and practitioners joined their efforts to produce a community-based terminology of WMSs. This paper summarizes their findings and introduces this new terminology to characterize WMSs. Furthermore, this terminology is composed of fives axes: workflow structure and characteristics, composition, orchestration, data management, and metadata capture. Each axis comprises several concepts that capture the prominent features of WMSs. Based on this terminology, this paper also presents a classification of 23 existing WMSs according to the proposed axes and terms.

Community-based terminology↗
DOE OSTI2025

Application of machine learning interatomic potentials in heterogeneous catalysis

Heterogeneous catalysts are crucial in modern societies as they promote sustainability by enabling lower-energy pathways for various chemical reactions. While Density Functional Theory (DFT) computations can provide critical insights into how heterogeneous catalysts operate at the atomic level, they are limited by computational costs and unfavorable scaling with system size. Recently, machine learning interatomic potentials (MLIPs) have emerged as a promising alternative to DFT, offering near-DFT accuracy at significantly reduced cost. Here, in this perspective, we discuss the application of MLIPs in heterogeneous catalyst modeling as a surrogate for DFT. We detail how MLIPs have been applied in thermal catalysis to probe active sites, enable studying complex metallic and nanoporous catalysts, and investigate the reconstruction of catalytic surfaces. We review the use of MLIPs in electrocatalysis and photocatalysis, emphasizing their capabilities in studying transition metal oxide surfaces and solid–liquid interfaces. We also discuss the current limitations of MLIPs, particularly their challenges with transferability and description of non-local interactions. Finally, we conclude by identifying promising and underexplored domains in which MLIPs can further advance our understanding of heterogeneous catalysts.

Catalytic surfaces↗
DOE OSTI2025

Charge Transfer and Recombination Pathways through Fullerene Guests in Porphyrin-Based MOFs

Porphyrin-based metal–organic frameworks (MOFs) offer a unique platform for building porous donor–acceptor networks that exhibit long-lived charge separation and transport upon incorporation of electron acceptor guest species. Here, in this work, porphyrin-based MOFs, PCN-222(H 2 ) and PCN-222(Zn), synthesized as nanoparticle suspensions, are successfully infiltrated with fullerene acceptor molecules, C 60 or PC 61 BM, in both polar and nonpolar solvent environments. The location and relative binding strength of these guest species are evaluated through a combination of N 2 physisorption measurements, photoluminescence quenching, and UV–vis absorption titration experiments. Semiempirical tight binding calculations are used to screen potential locations of the fullerene guest within the MOF pores, and hybrid density functional theory (DFT)-computed interaction energies confirm the energetically favorable positions. The fundamental photophysics of these donor–acceptor host–guest combinations are probed using ultrafast transient absorption spectroscopy. Sub-picosecond electron transfer involving initial exciplex population is observed, with slow charge recombination lifetimes on the order of τ ~1 ns for all systems in both dimethylformamide and 1,4-dioxane. Charge recombination occurs through population of fullerene and/or framework porphyrin triplet states depending on the porphyrin metalation status. The photophysics of the fullerene-loaded MOFs are discussed in the context of relevant porphyrin–fullerene donor–acceptor molecules to highlight the unique role of the framework environment in dictating photoinduced electron transfer and decay pathways.

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

Highly Conjugated Porphyrin Arrays Enable Optical Resolution of Ferromagnetic and Antiferromagnetic Aligned States of the Triplet Exciton and an Incorporated Stable Radical

Well-defined photogenerated molecular spin systems have potential utility in spintronics and quantum information science (QIS). Because molecular magnetic, optical, and electronic properties can be controlled by design, diverse spin systems can be prepared at modest temperatures. Photogenerated molecular spin systems often involve states prepared from the interaction of excitons and charges. Resolving the nature of electron spin alignment in photogenerated spin states described by the coupling of a triplet exciton and a stable radical commonly relies on EPR spectroscopy. Here, we describe ethyne-bridged (porphinato)metal ( PM n ) oligomers that incorporate a macrocycle-bound Cu(II) radical center. Upon photoexcitation of such PM n arrays, a singdoublet ( 2 S 1 ) state is formed; ultrafast internal conversion (IC) then produces a tripdoublet ( 2 T 1 ) state, which undergoes intersystem crossing (ISC) to produce a tripquartet ( 4 T 1 ) state, before relaxation to the ground state ( 2 S 0 ). These highly conjugated Cu(II) radical-containing PMn arrays enable direct observation of copper porphyrin 2 T 1 → 4 T 1 ISC dynamics from the biexponential decay of the near-infrared (NIR) 2,4 T 1 → 2,4 T n transient absorption manifold. Multireference n-electron valence perturbation theory (NEVPT2) computations illuminate how PM n electronic structure controls the relaxation dynamics of these long-lived (>10 ns) electronically excited multiplet states. These studies show that highly conjugated and polarizable porphyrin arrays incorporating stable spin centers provide rare π-delocalized systems where the ferromagnetic and antiferromagnetic alignment between a triplet exciton and a stable radical are both spectrally resolved and addressable using transient optical spectroscopy at wavelengths exceeding 1 μm, providing new opportunities to QIS.

14 SOLAR ENERGY↗
DOE OSTI2026

Artificial intelligence-driven approaches for materials design and discovery

Materials design is an important component of modern science and technology, yet traditional approaches rely heavily on trial and error and can be inefficient. Computational techniques, enhanced by modern artificial intelligence, have reshaped the landscape of designing new materials. Among these approaches, inverse design has shown great promise in designing materials that meet specific property requirements. Here, in this Review, we present key computational advances in materials design over the past few decades. We follow the evolution of relevant materials design techniques, from high-throughput forward machine learning methods and evolutionary algorithms, to advanced artificial intelligence strategies such as reinforcement learning and deep generative models. We highlight the paradigm shift from conventional screening approaches to inverse generation driven by deep generative models. Finally, we discuss current challenges and future perspectives of materials inverse design. This Review may serve as a brief guide to the approaches, progress and outlook of designing future functional materials with technological relevance.

computational methods↗
DOE OSTI2025

Circuit complexity and functionality: A statistical thermodynamics perspective

Circuit complexity, defined as the minimum circuit size required for implementing a particular Boolean computation, is a foundational concept in computer science. Determining circuit complexity is believed to be a hard computational problem. Recently, in the context of black holes, circuit complexity has been promoted to a physical property, wherein the growth of complexity is reflected in the time evolution of the Einstein-Rosen bridge (“wormhole”) connecting the two sides of an anti-de Sitter “eternal” black hole. Here, we are motivated by an independent set of considerations and explore links between complexity and thermodynamics for functionally equivalent circuits, making the physics-inspired approach relevant to real computational problems, for which functionality is the key element of interest. In particular, our thermodynamic framework provides an alternative perspective on the obfuscation of programs of arbitrary length—an important problem in cryptography—as thermalization through recursive mixing of neighboring sections of a circuit, which can be viewed as the mixing of two containers with “gases of gates.” This recursive process equilibrates the average complexity and leads to the saturation of the circuit entropy, while preserving functionality of the overall circuit. The thermodynamic arguments hinge on ergodicity in the space of circuits which we conjecture is limited to disconnected ergodic sectors due to fragmentation. The notion of fragmentation has important implications for the problem of circuit obfuscation as it implies that there are circuits of same size and functionality that cannot be connected via a polynomial number of local moves. Furthermore, we argue that fragmentation is unavoidable unless the complexity classes NP and coNP coincide, a statement that implies the collapse of the polynomial hierarchy of computational complexity theory to its first level.

Science & Technology - Other Topics↗
DOE OSTI2025

Dissociation of gas-phase anisole induced by low-energy electron interactions: understanding patterns of aromatic bond cleavage

Abstract In view of elucidating the fragmentation patterns of aromatic systems induced by low-energy electron interactions, dissociative electron attachment (DEA) to gas-phase anisole was performed. Anionic fragments resulting from this DEA process were detected by a quadrupole mass spectrometer, and ion yields of those fragments as a function of incident electron energy were rendered. Our study showed the formation of CH 3 − , HCC − , and OCH 3 − fragments, suggesting that various dissociation channels proceed out of DEA to anisole. We employed density functional theory to compute thermodynamic threshold energies for each potential dissociation channel. Those theoretical calculations supported the prediction that the CH 3 − and OCH 3 − fragments form via mechanisms of single-bond cleavage; the HCC − fragments may form through two-, three-, or four-body dissociation channels that entail hydrogen transfers and the cleavage of multiple aromatic bonds. The experimental resonance energies that form the CH 3 − , HCC − , and OCH 3 − fragments were 6.0 eV, 5.8 and 9.7 eV, and 9.8 eV, respectively. Given the classification of anisole as a monosubstituted aromatic species, our results explain generalizable patterns of electron-mediated dissociation in aromatic systems.

Finley, Jacob (ORCID:0009000606963352)↗
DOE OSTI2025

Communication Lower Bounds and Optimal Algorithms for Symmetric Matrix Computations

In this article, we focus on the communication costs of three symmetric matrix computations: (i) multiplying a matrix with its transpose, known as a symmetric rank-k update (SYRK) (ii) adding the result of the multiplication of a matrix with the transpose of another matrix and the transpose of that result, known as a symmetric rank-2k update (SYR2K) (iii) performing matrix multiplication with a symmetric input matrix (SYMM). All three computations appear in the Level 3 Basic Linear Algebra Subroutines (BLAS) and have wide use in applications involving symmetric matrices. We establish communication lower bounds for these kernels using sequential and distributed-memory parallel computational models, and we show that our bounds are tight by presenting communication-optimal algorithms for each setting. Our lower bound proofs rely on applying a geometric inequality for symmetric computations and analytically solving constrained nonlinear optimization problems. As a result, the symmetric matrix and its corresponding computations are accessed and performed according to a triangular block partitioning scheme in the optimal algorithms.

Al Daas, Hussam [Rutherford Appleton Laboratory, D↗
DOE OSTI2025

Tunable Few-Layer van der Waals Crystals and Heterostructures as Emerging Energy and Quantum Materials (Final Technical Report)

2D and layered (van der Waals) semiconductors offer extraordinary opportunities for manipulating optically excited charge carriers, many-body excitations, and non-charge based quantum numbers. To date, research has focused on a limited group of materials, mostly transition metal dichalcogenides in the monolayer limit. Other van der Waals semiconductors, and especially few-layer to multilayer crystals and their heterostructures, carry large potential for the discovery of phenomena of interest for future energy and information technologies. But they remain largely unexplored, often due to a lack of access to high-quality materials and approaches for measuring their properties at the relevant scales. The goal of this project was to develop an EPSCoR-State/National Laboratory Partnership that addresses the challenges of preparing high-quality van der Waals semiconductors and of probing their structure, composition, and especially their optoelectronic and photonic properties, near the atomic scale using electron microscopy techniques. A central component of the project was the development of advanced methods for electron microscopy and electron-excited spectroscopy, taking advantage of unique samples as well as leading capabilities and expertise at the partner institutions. Efforts to advance leading-edge techniques was supported by ancillary developments, such as precision sample preparation for electron microscopy/spectroscopy, coordinated chemical imaging, and analytical electron microscopy. Experiments in materials synthesis and technique development were closely linked to theory and computation. The results obtained under this project yielded multifaceted benefits to the involved partners and their institutions, DOE-BES, the wider scientific community, and society at large, particularly in the State of Nebraska through dividends from knowledge and human capital generated under the project.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Testing Classical Properties from Quantum Data

Many properties of Boolean functions can be tested far more efficiently than the function itself can be learned. However, this dramatic advantage often disappears when testers are limited to random samples of ƒ instead of adaptively chosen queries to f. In this work we investigate the quantum version of this restriction: quantum algorithms that test properties of a Boolean function f solely from copies of either the function state |ƒ⟩ ∝ ∑ x |x, ƒ(x)⟩ or the phase state |(-1) ƒ ⟩ ∝ ∑ x (-1) ƒ(x) |x⟩. For monotonicity, symmetry, and triangle-freeness, we show passive quantum testers are unboundedly or super-polynomially better than their classical passive testing counterparts. They are competitive with classic query -based testers in each case. Our new testers use techniques beyond quantum Fourier sampling, and it turns out this is necessary: we show a certain class of bent functions can be tested from 𝒪(1) function states but has a sample complexity lower bound of 2 Ω(n) for any tester relying exclusively on Fourier and classical samples. Our passive quantum testers are competitive with classical query -based testers, but this isn't universal: we exhibit a testing problem that can be solved from 𝒪(1) classical queries but requires Ω(2 n/2 ) function state copies. The Forrelation problem provides a separation of the same magnitude in the opposite direction, so we conclude that quantum data and classical queries are "maximally incomparable" resources for testing. We also begin the study of lower bounds for testing from quantum data. For quantum monotonicity testing, we prove that the ensembles of [Goldreich et al., 2000; Black, 2024], which give exponential lower bounds for classical sample-based testing, do not yield any nontrivial lower bounds for testing from quantum data. New insights specific to quantum data will be required for proving copy complexity lower bounds for testing in this model.

Boolean Functions↗
NASA NTRS1926

Navy propeller section characteristics as used in propeller design

This report contains artificial aerodynamic characteristics of a set of propeller sections to be used in designing propellers by means of the blade element theory. Characteristics computed from model propeller tests for a single section are extended to cover sections of Navy propeller sections at high Reynolds Number in the variable density tunnel of the NACA.

Weick, Fred E↗
NASA NTRS1953

A method for calculating the lift and center of pressure of wing-body-tail combinations at subsonic, transonic, and supersonic speeds

A method is presented for calculating the lift and pitching-moment characteristics of circular cylindrical bodies in combination with triangular, rectangular, or trapezoidal wings or tails through the subsonic, transonic, and supersonic speed ranges. The method covers unbanked wings, sweptback leading edges or sweptforward trailing edges, low angles of attack, and the effects of wing and tail incidence. The wing-body interference is handled by the method presented in NACA RM's A51J04 and A52B06, and the wing-tail interference is treated by assuming one completely rolled-up vortex per wing panel and evaluating the tail load by strip theory. A computing table and set of design charts are presented which reduce the calculations to routine operations. Comparison is made between the estimated and experimental characteristics for a large number of wing-body and wing-body-tail combinations. Generally speaking, the lifts were estimated to within plus-or-minus 10 percent and the centers of pressure were estimated to within plus-or-minus 0.02 of the body length. The effect of wing deflection on wing-tail interference at supersonic speeds was not correctly predicted for triangular wings with supersonic leading edges.

Nielsen, Jack N↗
DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Quantum Monte Carlo calculation of 𝛿 NS in 10 C using an effective field theory approach

Here, we compute radiative corrections to the superallowed 𝛽 decay of 10 C in an effective field theory approach using nuclear matrix elements obtained from quantum Monte Carlo calculations. These corrections are an important ingredient in the extraction of the Cabibbo-Kobayashi-Masakawa quark mixing matrix element 𝑉 𝑢⁢𝑑 , and the role of this work is to illuminate the uncertainties arising from nuclear structure. Our results provide good agreement with both the traditional extraction of 𝑉 𝑢⁢𝑑 , as well as with a more recent evaluation performed using the no-core shell model and a dispersion formalism. The dominant uncertainty in this approach is the presence of two unknown low-energy constants that enter into the relevant nuclear matrix elements. Future determinations of these low-energy constants—either from QCD or modeling them with two nucleon amplitudes—would improve the precision of the extraction in this formalism.

beta decay↗
DOE OSTI2025

Fault-tolerant operation and materials science with neutral atom logical qubits

We report on the fault-tolerant operation of logical qubits on a neutral atom quantum computer, with logical performance surpassing physical performance for multiple circuits including Bell state preparation (12x error reduction), random circuits (15x), and a prototype Anderson Impurity Model ground state solver for materials science applications (up to 6x, non-fault-tolerantly). The logical qubits are implemented via the [[4, 2, 2]] code (C 4 ). Our work constitutes the first complete realization of the benchmarking protocol proposed by Gottesman 2016 demonstrating results consistent with fault tolerance. In light of recent advances on applying concatenated C 4 /C 6 detection codes to achieve error correction with high code rates and thresholds, our work can be regarded as a building block towards a practical scheme for fault tolerant quantum computation. Our demonstration of a materials science application with logical qubits particularly demonstrates the immediate value of these techniques on current experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗