Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function approximation”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Green’s function coupled cluster formulations utilizing extended inner excitations

Here, we analyze new approximations of the Green’s function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i( n, m ) method, where the excitation level of the inner auxiliary operators ( m ) used to describe the ionization potential and electron affinity effects in the N-1 and N+1 particle spaces is higher than the excitation level ( n ) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called “n + 1” rule in this category [or the GFCC-i( n, n + 1) method], which states that in order to maintain size-extensivity of the Green’s function matrix elements, the excitation level of inner auxiliary operators X p (ω) and Y q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N 2 and CO molecules in the inner and outer valence regimes. Lastly, in comparison with the Green’s function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.

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

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Weak entanglement approximation for nuclear structure

The interacting shell model, a configuration-interaction method, is a venerable approach for low-lying nuclear structure calculations, but it is hampered by the exponential growth of its basis dimension as one increases the single-particle space and/or the number of active particles. Recent, quantum-information-inspired work has demonstrated that the proton and neutron sectors of a nuclear wave function are weakly entangled. Furthermore, the entanglement is smaller for nuclides away from N = Z, such as heavy, neutron-rich nuclides. Here, in this study, we implement a weak entanglement approximation to bipartite configuration-interaction wave functions, approximating low-lying levels by coupling a relatively small number of many-proton and many-neutron states. This truncation scheme, which we present in the context of past approaches, reduces the basis dimension by many orders of magnitude while preserving essential features of nuclear spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Formulation of the Rice-Walsh equation of state based on Shock Hugoniot data for porous metals

The dimensionless material parameter R introduced by Wu and Jing into the Rice-Walsh equation of state (EOS) has been deduced from the LASL shock Hugoniot data for porous Al and Cu. It was found that the parameter R/p decays smoothly with shock pressure p and displays small experimental scatter in the high pressure region. This finding led to the conclusion that the parameter has only a weak temperature dependence and is well approximated by a function of pressure alone, and the Grüneisen parameter should be temperature dependent under compression. The thermodynamic formulation of the Rice-Walsh EOS for Al and Cu was realized using the empirically determined function R(p) for each material and their known shock Hugoniot. It was then possible to reproduce porous shock Hugoniot for these metals. For most degrees of porosity, agreement between the porous data and the calculated Hugoniots using the empirical function described was very good. However, slight discrepancies were seen for Hugoniots with very high porosity. Two new thermal variables were introduced after further analysis, which enabled the calculation of the cold compression curve for these metals. The Grüneisen parameters along full-density and porous Hugoniot curve were calculated using a thermodynamic identity connecting R and the Grüneisen parameter. It was shown that the Grüneisen parameter is strongly temperature dependent. The present analysis suggested that the Rice-Walsh type EOS is a preferable choice for the analysis with its simple form, pressure-dependent empirical Wu-Jing parameter, and its compatibility with porous shock data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Closed-Form Approximation of the Total Variation Proximal Operator

Total variation (TV) is a widely used function for regularizing imaging inverse problems that is particularly appropriate for images whose underlying structure is piecewise constant. TV regularized optimization problems are typically solved using proximal methods, but the way in which they are applied is constrained by the absence of a closed-form expression for the proximal operator of the TV function. A closed-form approximation of the TV proximal operator has previously been proposed, but its accuracy was not theoretically explored in detail. Here, we address this gap by making several new theoretical contributions, proving that the approximation leads to a proximal operator of some convex function, it is equivalent to a gradient descent step on a smoothed version of TV, and that its error can be fully characterized and controlled with its scaling parameter. We experimentally validate our theoretical results on image denoising and sparse-view computed tomography (CT) image reconstruction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Theoretical prediction of the band offsets at the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions using the self-consistent ab initio DFT/GGA-1/2 method

The band offsets of the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions are calculated using the density functional theory/generalized gradient approximation (DFT/GGA)-1/2 method, which takes into account the self-energy corrections and can give an approximate description to the quasiparticle characteristics of the electronic structure of semiconductors. We present the results of the ionization potential (IP)-based and interfacial offset-based band alignments. In the interfacial offset-based band alignment, to get the natural band offset, we use the surface calculations to estimate the change of reference level due to the interfacial strain. Based on the interface models and GGA-1/2 calculations, we find that the valence band maximum and conduction band minimum of ZnO, respectively, lie 0.64 eV and 0.57 eV above those of anatase TiO{sub 2}, while lie 0.84 eV and 1.09 eV below those of GaN, which agree well with the experimental data. However, a large discrepancy exists between the IP-based band offset and the calculated natural band offset, the mechanism of which is discussed. Our results clarify band alignment of the ZnO/anatase TiO{sub 2} heterojunction and show good agreement with the GW calculations for the GaN/ZnO heterojunction.

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

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2023

A unified and constructive framework for the universality of neural networks

Abstract One of the reasons why many neural networks are capable of replicating complicated tasks or functions is their universal approximation property. Though the past few decades have seen tremendous advances in theories of neural networks, a single constructive and elementary framework for neural network universality remains unavailable. This paper is an effort to provide a unified and constructive framework for the universality of a large class of activation functions including most of the existing ones. At the heart of the framework is the concept of neural network approximate identity (nAI). The main result is as follows: any nAI activation function is universal in the space of continuous functions on compacta. It turns out that most of the existing activation functions are nAI, and thus universal. The framework induces several advantages over the contemporary counterparts. First, it is constructive with elementary means from functional analysis, probability theory, and numerical analysis. Second, it is one of the first unified and constructive attempts that is valid for most of the existing activation functions. Third, it provides new proofs for most activation functions. Fourth, for a given activation and error tolerance, the framework provides precisely the architecture of the corresponding one-hidden neural network with a predetermined number of neurons and the values of weights/biases. Fifth, the framework allows us to abstractly present the first universal approximation with a favorable non-asymptotic rate. Sixth, our framework also provides insights into the developments, and hence providing constructive derivations, of some of the existing approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

MCML: Combining physical constraints with experimental data for a multi-purpose meta-generalized gradient approximation

The predictive power of density functional theory for materials properties can be improved without increasing the overall computational complexity by extending the generalized gradient approximation (GGA) for electronic exchange and correlation to density functionals depending on the electronic kinetic energy density in addition to the charge density and its gradient, resulting in a meta-GGA. Here, we propose an empirical meta-GGA model that is based both on physical constraints and on experimental and quantum chemistry reference data. Overall, the resulting optimized meta-GGA MCML yields improved surface and gas phase reaction energetics without sacrificing the accuracy of bulk property predictions of existing meta-GGA approaches.

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

Search for Higgs boson pair production in the γγ$b\overline{b} $ final state with 13 TeV pp collision data collected by the ATLAS experiment

A search is performed for resonant and non-resonant Higgs boson pair production in the γγ$b\overline{b} $ final state. The data set used corresponds to an integrated luminosity of 36.1 fb -1 of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. No significant excess relative to the Standard Model expectation is observed. The observed limit on the non-resonant Higgs boson pair cross-section is 0.73 pb at 95% confidence level. This observed limit is equivalent to 22 times the predicted Standard Model cross-section. The Higgs boson self-coupling (κ λ = λ HHH /λ HHH SM ) is constrained at 95% confidence level to -8.2 < κ λ < 13.2. For resonant Higgs boson pair production through X→HH→γγ$b\overline{b} $, the limit is presented, using the narrow-width approximation, as a function of mX in the range 260 GeV < m X < 1000 GeV. The observed limits range from 1.1 pb to 0.12 pb over this mass range. Open image in new window

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Collisional Broadening of CO{sub 2} 10{sup 0}0–00{sup 0}1 Transition Lines by O{sub 2} and N{sub 2} in the Range 300–700 K

Unsaturated absorption coefficients are measured using a tunable CO{sub 2} laser in pure CO{sub 2} and CO{sub 2}:O{sub 2} = 1:2 and CO{sub 2}:N{sub 2} = 1:2 gas mixtures at the central frequencies of the R(8), R(22), P(8), P(22), and P(36) CO{sub 2} 10{sup 0}0–00{sup 0}1 transition lines under a pressure of 100 Torr in the temperature range 300–700 K. The collisional broadening coefficients of the CO{sub 2} spectral lines by O{sub 2} and N{sub 2} are obtained. The temperature dependences of these coefficients are approximated by power functions with two different exponents.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Assessing entropy for catalytic processes at complex reactive interfaces

When chemical reactions are accelerated by a catalyst, entropy differences between reactants and their transient intermediates can be the driving force behind the promotion or inhibition of desired and parasitic chemical pathways. Understanding and controlling catalytic processes therefore requires both a fundamental and practicable understanding of entropy in addition to enthalpy. In unstructured media such as the vapor phase equilibrated with sparsely covered surfaces, entropy can be adequately accounted for by well-established approaches based on translational, rotational, and harmonic vibrational partition functions. However, these approximations become inadequate in more complex condensed phase environments, e.g., solid liquid interfaces of confined reaction spaces. In this chapter, we provide an overview of the state-of-art in the computational quantification of entropy and its known ramifications on catalysis. The fundamental roles of thermodynamics and kinetics in catalysis are covered in enough detail to appreciate and contextualize the computational methods employed to compute chemically accurate estimates of entropy. These methods are discussed in appropriate detail and range from the ubiquitous harmonic oscillator approximation where entropy unrelated to high frequency oscillations is typically underestimated, to enhanced free energy sampling with molecular dynamics where the desired accuracy must be weighed against the associated computational cost of obtaining it. The rising importance of machine learning and artificial intelligence in accelerating methodological progress in this field is touched upon, as well. Finally, applications, successes, and pitfalls of using these methods are provided to showcase past and present accomplishments while clarifying where improvements in both understanding and methodology are still needed.

Kollias, Loukas↗
DOE OSTI2021

A two-dimensional electron gas suspended above a neutralizing background

Highlights: • Two-dimensional electron gas above a neutralizing jellium background. • Analytic results for the value of the total energy. • Energy contribution originating from separation. • Added energy behaves as an effective kinetic energy term. We consider a two-dimensional electron gas system consisting of a layer of electrons suspended above its jellium neutralizing background counterpart. We investigate the energy of the system in the Hartree–Fock approximation as a function of separation distance between the layer of electrons and the background. We find that the energy of the system is influenced by the separation distance in such a way that the new energy contribution originating from such a separation behaves as an effective kinetic energy term. The increase of the overall energy with increase of the separation distance adversely affects the stability of electron liquid states at high density. The results are derived for a fully spin-polarized (spinless) system of electrons in the thermodynamic limit, but can easily be generalized to a spin unpolarized state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Uniform-in-phase-space data selection with iterative normalizing flows

Improvements in computational and experimental capabilities are rapidly increasing the amount of scientific data that are routinely generated. In applications that are constrained by memory and computational intensity, excessively large datasets may hinder scientific discovery, making data reduction a critical component of data-driven methods. Datasets are growing in two directions: the number of data points and their dimensionality. Whereas dimension reduction typically aims at describing each data sample on lower-dimensional space, the focus here is on reducing the number of data points. A strategy is proposed to select data points such that they uniformly span the phase-space of the data. The algorithm proposed relies on estimating the probability map of the data and using it to construct an acceptance probability. An iterative method is used to accurately estimate the probability of the rare data points when only a small subset of the dataset is used to construct the probability map. Instead of binning the phase-space to estimate the probability map, its functional form is approximated with a normalizing flow. Therefore, the method naturally extends to high-dimensional datasets. The proposed framework is demonstrated as a viable pathway to enable data-efficient machine learning when abundant data are available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Structural and Electronic Effects on the Properties of Fe 2 (dobdc) upon Oxidation with N 2 O

Here, we report electronic, vibrational, and magnetic properties, together with their structural dependences, for the metal–organic framework Fe 2 (dobdc) (dobdc 4 – = 2,5-dioxido-1,4-benzenedicarboxylate) and its derivatives, Fe 2 (O) 2 (dobdc) and Fe 2 (OH) 2 (dobdc)—species arising in the previously proposed mechanism for the oxidation of ethane to ethanol using N 2 O as an oxidant. Magnetic susceptibility measurements reported for Fe 2 (dobdc) in an earlier study and reported in the current study for Fe II 0.26 [Fe III (OH)] 1.74 (dobdc)(DMF) 0.15 (THF) 0.22 , which is more simply referred to as Fe 2 (OH) 2 (dobdc), were used to confirm the computational results. Theory was also compared to experiment for infrared spectra and X-ray powder diffraction structures. Structural and magnetic properties were computed by using Kohn-Sham density functional theory both with periodic boundary conditions and with cluster models. In addition, we studied the effects of different treatments of the exchange interactions on the magnetic coupling parameters by comparing several approaches to the exchange-correlation functional: generalized gradient approximation (GGA), GGA with empirical Coulomb and exchange integrals for 3d electrons (the GGA+U method), nonseparable gradient approximation with empirical Coulomb and exchange integrals for 3d electrons (GAM+U), hybrid GGA, meta-GGA, and hybrid meta- GGA. We found the coupling between the metal centers to be ferromagnetic in the case of Fe 2 (dobdc) and antiferromagnetic in the cases of Fe 2 (O) 2 (dobdc) and Fe 2 (OH) 2 (dobdc). The shift in magnetic coupling behavior correlates with the changing electronic structure of the framework, which derives from both structural and electronic changes that occur upon metal oxidation and addition of the charge-balancing oxo and hydroxo ligands.

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

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium–Ligand Bonds

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy D e and bond length r e are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15–0.19. For the single-reference methods, Kohn–Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the D e and r e of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core–valence terms; this results in highly accurate dissociation energies and good bond lengths. In conclusion, averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

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

Extremely Slow Trap-Mediated Hole Relaxation in Room-Temperature InP/ZnSe/ZnS Quantum Dots

Time-resolved photoluminescence (PL) and transient absorption (TA) spectroscopy have been used to elucidate the hole relaxation dynamics in high-quality InP/ZnSe/ZnS quantum dots (QDs) at room temperature. In QDs having an In/P ratio close to 1:1, hole cooling is fast compared to the ~40 ps time-correlated photon-counting instrument response. However, in QDs that luminesce at wavelengths greater than around 600 nm and have excess indium, a significant fraction of the PL exhibits a slower rise time. In this case, the PL rises over a range of timescales varying from close to the instrument response to approximately 300–400 ps. The time constant of the slowest component increases with increasing InP core size and the fraction of the PL that has the slow rise time increases with excess indium in the particle. In all cases, a slow rise time and any spectral evolution on that timescale are absent in the TA results, indicating that the slow dynamics may be assigned to relaxation in the valence band. QDs emitting at wavelengths below 560 nm show very different dynamics, exhibiting a barely detectable <50 ps rise time. Based on these results and measured elemental compositions, the slow hole relaxation is assigned to the transient population of indium-based traps in the ZnSe shell that are above the valence band edge only for the larger, red-emitting QDs. Here, the relative trap to band edge relaxation rates can be understood in terms of effective mass approximation hole wave function calculations, which show that the relaxation rate is largely determined by the extent to which the hole wave function in the lowest exciton penetrates into the ZnSe shell, overlapping the trapped-hole wave functions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗