Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical physics methods”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effect of Framework Composition and NH 3 on the Diffusion of Cu + in Cu-CHA Catalysts Predicted by Machine-Learning Accelerated Molecular Dynamics

Cu-exchanged zeolites rely on mobile solvated Cu + cations for their catalytic activity, but the role of the framework composition in transport is not fully understood. Ab initio molecular dynamics simulations can provide quantitative atomistic insight but are too computationally expensive to explore large length and time scales or diverse compositions. We report a machine-learning interatomic potential that accurately reproduces ab initio results and effectively generalizes to allow multinanosecond simulations of large supercells and diverse chemical compositions. Biased and unbiased simulations of [Cu(NH 3 ) 2 ] + mobility show that aluminum pairing in eight-membered rings accelerates local hopping and demonstrate that increased NH 3 concentration enhances long-range diffusion. The probability of finding two [Cu(NH 3 ) 2 ] + complexes in the same cage, which is key for SCR-NOx reaction, increases with Cu content and Al content but does not correlate with the long-range mobility of Cu + . Supporting experimental evidence was obtained from reactivity tests of Cu-CHA catalysts with a controlled chemical composition.

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

Au as a Surrogate for F: The Case of UAu 6 vs UF 6

Here, anion photoelectron spectroscopy and first-principles quantum chemistry are used to demonstrate to what degree Au can act as a surrogate for F in UF6 and its anion. Unlike UF6, UAu6 exhibits strong ligand–ligand, i.e., Au–Au, interactions, resulting in three low-lying isomers, two of which are three-dimensional while the third isomer has a ring-like quasi two-dimensional structure. Additionally, all the UAu6 isomers have open-shell electrons, which in nearly all cases are localized on the central U atom. As a result, the adiabatic electron affinity and vertical detachment energy are measured to be 3.05 ± 0.05 and 3.28 ± 0.05 eV, respectively, and are in very good agreement with calculations.

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

Speciation and Reactivity Control of Cu-Oxo Clusters via Extraframework Al in Mordenite for Methane Oxidation

The stoichiometric conversion of methane to methanol by Cu-exchanged zeolites can be brought to highest yields by the presence of extraframework Al and high CH 4 chemical potentials. Here, combining theory and experiments, the differences in chemical reactivity of monometallic Cu-oxo and bimetallic Cu-Al-oxo nanoclusters stabilized in zeolite mordenite (MOR) are investigated. Cu-L 3 edge X-ray absorption near-edge structure (XANES), infrared (IR), and ultraviolet–visible (UV–vis) spectroscopies, in combination with CH 4 oxidation activity tests, support the presence of two types of active clusters in MOR and allow quantification of the relative proportions of each type in dependence of the Cu concentration. Ab initio molecular dynamics (MD) calculations and thermodynamic analyses indicate that the superior performance of materials enriched in Cu-Al-oxo clusters is related to the activity of two μ-oxo bridges in the cluster. Replacing H 2 O with ethanol in the product extraction step led to the formation of ethyl methyl ether, expanding this way the applicability of these materials for the activation and functionalization of CH 4 . We show that competition between different ion-exchanged metal-oxo structures during the synthesis of Cu-exchanged zeolites determines the formation of active species, and this provides guidelines for the synthesis of highly active materials for CH 4 activation and functionalization.

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

Regioselective On-Surface Synthesis of [3]Triangulene Graphene Nanoribbons

The integration of low-energy states into bottom-up engineered graphene nanoribbons (GNRs) is a robust strategy for realizing materials with tailored electronic band structure for nanoelectronics. Low-energy zero-modes (ZMs) can be introduced into nanographenes (NGs) by creating an imbalance between the two sublattices of graphene. This phenomenon is exemplified by the family of [n]triangulenes (n ϵ $\mathbb{N}$). Here, we demonstrate the synthesis of [3]triangulene-GNRs, a regioregular one-dimensional (1D) chain of [3]triangulenes linked by five-membered rings. Hybridization between ZMs on adjacent [3]triangulenes leads to the emergence of a narrow band gap, E g,exp ~ 0.7 eV, and topological end states that are experimentally verified using scanning tunneling spectroscopy. Tight-binding and first-principles density functional theory calculations within the local density approximation corroborate our experimental observations. Our synthetic design takes advantage of a selective on-surface head-to-tail coupling of monomer building blocks enabling the regioselective synthesis of [3]triangulene-GNRs. Detailed ab initio theory provides insights into the mechanism of on-surface radical polymerization, revealing the pivotal role of Au-C bond formation/breakage in driving selectivity.

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

Single-Molecule Detection of the Encounter and Productive Electron Transfer Complexes of a Photosynthetic Reaction Center

Small, diffusible redox proteins play an essential role in electron transfer (ET) in respiration and photosynthesis, sustaining life on Earth by shuttling electrons between membrane-bound complexes via finely tuned and reversible interactions. Ensemble kinetic studies show transient ET complexes form in two distinct stages: an “encounter” complex largely mediated by electrostatic interactions, which subsequently, through subtle reorganization of the binding interface, forms a “productive” ET complex stabilized by additional hydrophobic interactions around the redox-active cofactors. Here, using single-molecule force spectroscopy (SMFS) we dissected the transient ET complexes formed between the photosynthetic reaction center-light harvesting complex 1 (RC-LH1) of Rhodobacter sphaeroides and its native electron donor cytochrome c 2 (cyt c 2 ). Importantly, SMFS resolves the distribution of interaction forces into low (~150 pN) and high (~330 pN) components, with the former more susceptible to salt concentration and to alteration of key charged residues on the RC. Thus, the low force component is suggested to reflect the contribution of electrostatic interactions in forming the initial encounter complex, whereas the high force component reflects the additional stabilization provided by hydrophobic interactions to the productive ET complex. Employing molecular dynamics simulations, we resolve five intermediate states that comprise the encounter, productive ET and leaving complexes, predicting a weak interaction between cyt c 2 and the LH1 ring near the RC-L subunit that could lie along the exit path for oxidized cyt c 2 . The multimodal nature of the interactions of ET complexes captured here may have wider implications for ET in all domains of life.

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

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Electrochemical Hydrogen Compression Cell Design (CRADA Final Report)

Electrochemical compression and purification is a high potential technology that is likely to outcompete mechanical compression for hydrogen end use applications. However, testing and understanding what is occurring within the cell is difficult as it is not readily accessible through experimentation. Thus, the main target is to increase performance and reduce cost of electrochemical hydrogen compression (EHC) and purification (EHP) cells and stacks using multiscale mathematical modeling informed by experiment. The models can provide physical insight into limiting processes and develop mitigation methods, which can be proven out through rapid sub-scale testing of the devices.

08 HYDROGEN↗
DOE OSTI2020

Uses of Sensitivity/Uncertainty Techniques for Critical Experiment Design [Slides]

This presentation discusses the programming languages TSUNAMI and TSAR. It also provides a review of experimental designs using sensitivity analysis in regard to 7uPCX fuel design, fission product experiment analysis, MIRTE, and temperature dependent experiments at the Sandia Pulsed Reactor Facility (SPRF).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1980

An improved panel method for the solution of three-dimensional leading-edge vortex flows. Volume 1: Theory document

An improved panel method for the solution of three dimensional flow and wing and wing-body combinations with leading edge vortex separation is presented. The method employs a three dimensional inviscid flow model in which the configuration, the rolled-up vortex sheets, and the wake are represented by quadratic doublet distributions. The strength of the singularity distribution as well as shape and position of the vortex spirals are computed in an iterative fashion starting with an assumed initial sheet geometry. The method calculates forces and moments as well as detail surface pressure distributions. Improvements include the implementation of improved panel numerics for the purpose of elimination the highly nonlinear effects of ring vortices around double panel edges, and the development of a least squares procedure for damping vortex sheet geometry update instabilities. A complete description of the method is included. A variety of cases generated by the computer program implementing the method are presented which verify the mathematical assumptions of the method and which compare computed results with experimental data to verify the underlying physical assumptions made by the method.

Johnson, F. T.↗
NASA NTRS1987

Design and validation of fault-tolerant flight systems

Flight systems must be validated to show that they are consistent with the requirements of their intended applications. While high reliability is difficult to validate, the additional complexity of fault tolerance further compounds the validation problem. The objective of NASA’s research is to develop a methodology for designing validatable fault-tolerant systems. Under the design-for-validation philosophy, emphasis is placed on developing validation methods that can be incorporated into the design process right from the start and design methods and guidance which, while incorporating fault tolerance, can assure validatability. This paper examines the statistical issues of validating highly reliable, fault tolerant system. There are many problems associated with traditional methods of designing and validating these potentially complex hardware and software systems. Useful design-for-validation methods, which include structured specification and design methodologies, mathematical proof techniques, analytical modeling, simulation and emulation, and physical testing, are discussed. Important design issues associated with fault tolerance are presented along with the related validation concerns which must be addressed. Experience has shown that synchronization and Byzantine resilience must accompany fault tolerance. Other design attributes associated with fault tolerance may be used by a designer on the basis of cost, weight, performance, and validation considerations.

Computer systems↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Coarse Graining Discrete Element Method Information in Particle-in-Cell Length Scales Using a Machine Learning Approach

This report details the development of a machine learning (ML)-driven framework to coarse-grain inter-particle collision dynamics from high-fidelity Discrete Element Method (DEM) simulations to Particle-in-Cell (PIC) scales for gas-solid systems. Traditional PIC models, while computationally efficient, rely on empirical granular stress formulations that fail to capture the full complexity of collision physics, particularly the heterogeneity in particle dynamics. This study adopts a bottom-up approach, integrating insights from DEM simulations to improve the physical fidelity and interpretability of PIC-scale models.

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

Discontinuous Galerkin sparse grids methods for time domain Maxwell’s equations

In this paper, we present a time-explicit sparse grid discontinuous Galerkin method for solving the three-dimensional time-domain Maxwell equations. The conservation properties and convergence rates are established for different choices of numerical fluxes. The convergence rates are proved theoretically and then verified by several numerical examples. Even though our scheme does not preserve the divergence, but by implying the higher order polynomial, one can observe the same convergence rate as numerical solution. Several numerical tests are presented to validate these conclusions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Extracting Inelastic Scattering Cross Sections for Finite and Aperiodic Materials from Electronic Dynamics Simulations

Explicit time-dependent electronic structure theory methods are increasingly prevalent in the areas of condensed matter physics and quantum chemistry, with the broad-band optical absorptivity of molecular and small condensed-phase systems nowadays routinely studied with such approaches. Here, in this paper, it is demonstrated that electronic dynamics simulations can similarly be employed to study cross sections for the scattering-induced electronic excitations probed in nonresonant inelastic X-ray scattering and momentum-resolved electron energy loss spectroscopies. A method is put forth for evaluating the electronic dynamic structure factor, which involves the application of a momentum boost-type perturbation and transformation of the resulting reciprocal space density fluctuations into the frequency domain. Good agreement is first demonstrated between the dynamic structure factor extracted from these electronic dynamics simulations and the corresponding transition matrix elements from linear response theory. The method is then applied to some extended (quasi)one-dimensional systems, for which the wave vector becomes a good quantum number in the thermodynamic limit. Finally, the dispersion of many-body excitations in a series of hydrogen-terminated graphene flakes (and twisted bilayers thereof) is investigated to highlight the utility of the presented approach for capturing morphology-dependent effects in the inelastic scattering cross sections of nanostructured and/or noncrystalline materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Nonlinear ARMA models for the D(st) index and their physical interpretation

Time series models successfully reproduce or predict geomagnetic activity indices from solar wind parameters. A method is presented that converts a type of nonlinear filter, the nonlinear Autoregressive Moving Average (ARMA) model to the nonlinear damped oscillator physical model. The oscillator parameters, the growth and decay, the oscillation frequencies and the coupling strength to the input are derived from the filter coefficients. Mathematical methods are derived to obtain unique and consistent filter coefficients while keeping the prediction error low. These methods are applied to an oscillator model for the Dst geomagnetic index driven by the solar wind input. A data set is examined in two ways: the model parameters are calculated as averages over short time intervals, and a nonlinear ARMA model is calculated and the model parameters are derived as a function of the phase space.

Vassiliadis, D.↗