Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Probing the Chemical Kinetics of Minimalist Functional Group Gasoline Surrogates

We report surrogate mixtures are routinely used for understanding gasoline fuel combustion in engine simulations. The general trend in surrogate formulation has been to increase the number of fuel components in a mixture to better emulate real fuel properties. Recently, a new surrogate design strategy based on functional group analysis of real gasolines was proposed using a minimal number of species [minimalist functional group (MFG)-approach]. MFG surrogates (having just one or two components) could experimentally capture the ignition delay time (IDT), threshold sooting index, and smoke point of different gasoline fuels with hundreds of components. However, other combustion characteristics were not explored, and kinetic modeling of MFG surrogates was not reported. These aspects are addressed in this paper, where the combustion behavior of MFG surrogates for various gasolines was assessed by simulating IDT, jet-stirred reactor oxidation, and premixed laminar flame speeds using chemical kinetic modeling. MFG simulations were compared with experimental data of the real gasolines as well as with the more complex multicomponent (five to nine species) surrogates. This study reveals that binary MFG surrogate mixtures are capable of accurately simulating the combustion behavior of more complex gasoline fuels with hundreds of components.

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

Fermi-liquid theory and divergences of the two-particle irreducible vertex in the periodic Anderson lattice

Here we analyze the divergences of the irreducible vertex function in dynamical mean field theory, which may indicate either a nonphysical breakdown of the perturbation theory or a response to some physical phenomenon. To investigate this question, we construct a quasiparticle vertex from the diverging irreducible vertex functions. This vertex describes the scattering between quasiparticles and quasiholes in a Fermi liquid. We show that the quasparticle vertex does not diverge in the charge channel, wherein the irreducible vertex does diverge; and we show that the quasiparticle vertex does diverge in the spin channel, wherein the irreducible vertex does not diverge. This divergence occurs at the Mott transition, wherein the Fermi-liquid theory breaks down. Both the half filled Hubbard and Anderson lattices are investigated. In general, our results support that the divergences of the irreducible vertex function do not indicate a nonphysical failure of the perturbation theory. Instead, the divergences are the mathematical consequence of inverting a matrix (the local charge susceptibility) which accumulates increasingly negative diagonal elements as the Hubbard interaction suppresses charge fluctuations. Indeed, we find that the first divergences of the irreducible vertex in both Hubbard and Anderson lattices occurs near the maximum magnitude of the (negative) vertex-connected part of the charge susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Convergent momentum-space OPE and bootstrap equations in conformal field theory

General principles of quantum field theory imply that there exists an operator product expansion (OPE) for Wightman functions in Minkowski momentum space that converges for arbitrary kinematics. This convergence is guaranteed to hold in the sense of a distribution, meaning that it holds for correlation functions smeared by smooth test functions. The conformal blocks for this OPE are conceptually extremely simple: they are products of 3-point functions. We construct the conformal blocks in 2-dimensional conformal field theory and show that the OPE in fact converges pointwise to an ordinary function in a specific kinematic region. Using microcausality, we also formulate a bootstrap equation directly in terms of momentum space Wightman functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Pair breaking in superconductors with strong spin-orbit coupling

Here we study the influence of symmetry-breaking perturbations on superconductivity in multiorbital materials, with a particular focus on an external magnetic field. We introduce the field-fitness function which characterizes the pair-breaking effects of the perturbation on a given superconducting state. For even-parity superconductors we find that this field-fitness function for an external magnetic field is one, implying that the paramagnetic response is controlled only by a generalized effective g factor. For odd-parity superconductors, the interplay of the effective g factor and the field-fitness function can lead to counterintuitive results. We demonstrate this for p-wave pairing in the effective j = $\frac{3}{2}$ electronic states of the Luttinger-Kohn model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Chapter 10: Advances in Protein Engineering and its Application in Synthetic Biology

Protein engineering has been used successfully in fields ranging from medicine to food science to biofuels. Applications of protein engineering include developing antiviral peptides or other protein therapeutics, antibody engineering, designing protein-based logic circuits, engineering enzymes to be more specific or to function under industrially relevant conditions such as at higher temperatures or high/low pH, modifying cell signaling or regulatory functions, and so on. Advances in recombinant DNA, "omics," and CRISPR-Cas (clustered regularly interspaced short palindromic repeats and its associated proteins) technologies, combined with high-throughput screening facilities, will lead to improved methods for protein engineering, enabling easy modification of more proteins/enzymes for new specific applications. New methods for rational design, directed evolution, and computer-aided protein design will further accelerate the speed of protein evolution and expand the scope for protein engineering. In this chapter we discuss general protein engineering strategies and advances in engineering proteins with desired functions, focusing on the "design" and "build" part of the design-build-test-learn cycle.

BIOMASS FUELS↗
DOE OSTI2023

Combined Experimental and Computational Efforts to Establish Ion Mobility, Solubility and Stability of Functional Liquids for Electrochemical Energy Storage

This work provides a computation-driven investigation of the stability of organic electrolytes for lithium-air batteries. Electrolyte instability is currently a key challenge that limits practical use of aprotic Li-air batteries, and the chemical processes that cause this instability are often kinetically-driven. Computational screening for kinetic stability involves the determination of reaction barriers for the numerous potential reaction mechanisms, barriers that are challenging to calculate due to the difficulty of locating transition state structures. Here we screen a broad set of substituted electrolytes for susceptibility to nucleophilic attack by superoxide. We find that carbonates are not typically expected to be stable and that sulfones are generally stable, validating literature trends. We study the effects of chemical functionalization with electron-donating and withdrawing groups and their interplay with steric factors, identifying functional groups and other chemical modifications that increase stability in these groups. User-input driven transition state identification is used for these initial calculations, and an automated computational pipeline is subsequently presented and validated as a means to perform further high-throughput searches across mechanisms and chemistries. The pipeline integrates cheminformatics-based reaction encoding, relaxed potential energy scans, and nudged elastic band calculations for an end-to-end approach to barrier calculations. We review this automated search approach and its current limitations, and discuss challenges and further work.

25 ENERGY STORAGE↗
DOE OSTI2025

Anti-symmetric barron functions and their approximation with sums of determinants

A fundamental problem in quantum physics is to encode functions that are completely anti-symmetric under permutations of identical particles. The architecture of neural network models for the electron wave function typically comprises an equivariant component followed by a summation of determinants. The recently introduced Generic Antisymmetric (GA) block is designed to enhance the expressivity of such neural wave functions, and it was found that the 2-layer GA block achieved more accurate energies than the corresponding single-determinant FermiNet architecure, suggesting its promise as a way to improve the expressivity of neural wave functions. In this paper we show how the function expressed by the 2-layer GA block can be decomposed into a sum of determinants. We formalize this result by defining the antisymmetric Barron space as a generalized version of the 2-layer GA block and providing an appromation theorem for this function class. This result can be viewed as a negative result showing that the 2-layer GA block is not more expressive than using multiple determinants.

Abrahamsen, Nilin↗
DOE OSTI2021

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗
DOE OSTI2026

Neural Scaling Laws of Deep ReLU and Deep Operator Network: A Theoretical Study

Neural scaling laws play a pivotal role in the performance of deep neural networks and have been observed in a wide range of tasks. However, a complete theoretical framework for understanding these scaling laws remains underdeveloped. In this paper, we explore the neural scaling laws for deep operator networks, which involve learning mappings between function spaces, with a focus on the Chen and Chen style architecture. These approaches, which include the popular Deep Operator Network (DeepONet), approximate the output functions using a linear combination of learnable basis functions and coefficients that depend on the input functions. We establish a theoretical framework to quantify the neural scaling laws by analyzing its approximation and generalization errors. We articulate the relationship between the approximation and generalization errors of deep operator networks and key factors such as network model size and training data size. Moreover, we address cases where input functions exhibit low-dimensional structures, allowing us to derive tighter error bounds. These results also hold for deep ReLU networks and other similar structures. Our results offer a partial explanation of the neural scaling laws in operator learning and provide a theoretical foundation for their applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Revisiting W–ZrC interfaces: A first principles study

We investigate W–ZrC interfaces using first-principles calculations based on the density functional theory. There have been theoretical and experimental studies exploring W–ZrC interfaces, however, the debate regarding the most stable interface continues to persist. In this study, we systematically simulated various W–ZrC interfaces merging W and ZrC surfaces with different orientations. Subsequently, we evaluated their stabilities and explained the corresponding stabilities in terms of the nature of bonding and charge-transfer processes at the interface. We find ZrC(111)–W(110) is the most stable interface with higher adhesive energy than the other interfaces. The additional stability associated with the ZrC(111)–W(110) results from significant interface reconstruction. Three layers of W and ZrC adjacent to the interface are involved in the charge-transfer process leading to stronger ionic bonds in ZrC(111)–W(110) as compared to the other potential candidate: ZrC(100)–W(100). The C and W atoms are found to be displaced from their symmetric position during the reconstruction process at the interface to facilitate stronger bonds with shorter W–C and W–Zr bonds in ZrC(111)–W(110) as compared to ZrC(100)–W(100). This leads to stronger covalent bonds in ZrC(111)–W(110) than that in ZrC(100)–W(100). Therefore, we conclude that the stronger covalent and ionic forces in ZrC(111)–W(110) than those in ZrC(100)–W(100) are responsible for making ZrC(111)–W(110) to be the most stable interface. In conclusion, this study addresses the long-standing question of the most stable W–ZrC interface and derives a number of implications for other W-transition metal carbide interfaces which are potential candidates for improving the mechanical properties of plasma facing materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wave-packet manipulation of He Rydberg states by a seeded free-electron laser

We report a two-dimensional pump-control-probe spectroscopic study of the dynamics of singly excited He Rydberg-state wave packets with a seeded extreme ultraviolet (XUV) free-electron laser (FEL) source. A pair of coherent XUV pulses, defined by their coarse time separation and relative phase, created and manipulated the wave packets. The He atoms were postionized by infrared (IR) pulses, and the ion yield was measured as a function of XUV phase and IR arrival time. We tagged and sorted the relative phase of the XUV pulse pair on a single-shot basis by fitting each FEL spectrum with a suitable function that accounts for nonidealities of the XUV pulse pairs, associated with the seeding process; more generally, the fit returns the time-dependent electric field of the FEL spectra. The experimental two-dimensional maps of ion yields, measured as a function of IR (probe) delay and of XUV (pump-control) phase, were compared with the solution of the first-order time-dependent Schrödinger equation for this field. Despite the fact that the experimental conditions imply strong excitation, beyond the approximations of first-order perturbation theory, the simulated map satisfactorily reproduces the experimental one for temporally well-separated pulses. We show that by selecting data at appropriate values of pump-control phase, we enhance or suppress the amplitude of chosen wave-packet components consisting of two or more Rydberg states. When the temporal overlap of the pulse pair cannot be neglected, the phase reconstruction is underdetermined, and we provide a simplified comparison between data and simulations. Published by the American Physical Society 2024

Dumergue, Mathieu (ORCID:0000000236984406)↗
DOE OSTI2023

The five-point bootstrap

We study five-point correlation functions of scalar operators in d-dimensional conformal field theories. We develop a new approach to computing the five-point conformal blocks for exchanged primary operators of arbitrary spin by introducing a generalization of radial coordinates, using an appropriate ansatz, and perturbatively solving two quadratic Casimir differential equations. We then study five-point correlators 〈σσϵσσ〉 in the critical 3d Ising model. We truncate the operator product expansions (OPEs) in the correlator by including a finite number of primary operators with conformal dimension below a cutoff ∆ ⩽ ∆ cutoff . We then compute several OPE coefficients involving ϵ and two spinning operators by demanding that the truncated correlator approximately satisfies the crossing relation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

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

Uniting Theory and Experiment to Deliver Flexible MOFS for Superior Methane (NG) Storage

The objective of the project was to use previous insights developed through synthesis and quantitative modeling of rigid metal–organic frameworks (MOFs) in an established synergistic theoretical/experimental team to create, modify, and evaluate flexible MOFs (FlexMOFs) for natural gas (NG) storage and release at practically useful pressures and transform the NG storage economy. The goal was reduced pressure absorbed natural gas (ANG) FlexMOF storage at operating pressures less than 100 bar with physisorption exploiting the favorable thermodynamics and kinetics of flexible porous material opening in response to adsorption. Hydrogen behavior in this context was also considered. The specific aim of the project was to design and develop a standard computational modeling methodology for, first, detailed atomistic retrodiction of FlexMOF gating behavior and ultimately prediction of the effects of functionalization and/or substitution on structural transitioning. The project was also geared towards the establishment of the interaction of methane with the framework and binding sites using modeling and the FlexMOF will be both internally and external validated as a SMART metric. This has value to the scientific community both from the obvious standpoint of providing a better understanding of the promising test case systems (CdIF-13 and the MIL-53(Al) series of MOFs), but also in providing an avenue of obtaining insight into these FlexMOF systems in general, which is of particular interest given the tendency of structure-function correlation to lag behind synthesis methodology (making the latter a hit or miss proposition for applications). This gap remains significant for FlexMOF systems whose gate opening behaviors complicate computational examination. The computational methodologies are relatively inexpensive in terms of both money and computational resources enhancing the general viability of these methodologies. The ultimate gain to the public will be in the application of these techniques to design systems for natural gas storage and use to cut down on green-house emissions.

03 NATURAL GAS↗
DOE OSTI2024

Open-closed string duality, branes, and topological recursion

We consider matrix models exhibiting open-closed string duality in two-dimensional string theories with various amounts of supersymmetry. In particular, a relationship between matrix models in the β = 2 Wigner-Dyson class and models in the (1 + 2Γ, 2) Altland-Zirnbauer class relates the perturbative solutions of the two systems’ string equations. Point-like operator insertions in the closed string theory are mapped to the topological expansion of the free energy in the open string theory. We compute correlation functions of macroscopic loop operators and FZZT branes in a general topological gravity background. The relationship between the topological recursion of moduli space volumes and branes is discussed by analyzing the Virasoro conditions in the matrix models.

2D Gravity↗