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

Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation

Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more-stable bio-oil product. HDO reactions are performed by co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Noble-metal catalysts supported on reducible metal oxides (e.g., Pt/TiO2) are active and selective toward these HDO reactions. Griffin and co-workers showed that Pt/TiO2 catalysts promote the desired deoxygenation steps for m-cresol HDO while mitigating undesired C-C bond-breaking steps that reduce the overall value/energy density of the biofuel. Such model-compound studies to inform the design of improved catalysts for HDO chemistry are necessary to improve the overall process economics/efficiencies for biofuels production. One important class of bio-derived compounds that has not been studied extensively with respect to HDO chemistry, particularly at the atomic level, is carboxylic acids. Past experimental work indicates that Pt/C and Pt/TiO2 catalysts are selective toward C-C and C-O bond-dissociation products for acetic acid HDO, respectively. To understand the role of Pt/TiO2 active sites in this observed change in selectivity and guide catalyst development, the work in this presentation focuses on modeling the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in acetic acid HDO. Density functional theory (DFT) calculations were performed using the Vienna Ab initio Simulation Package (VASP). The exchange-correlation functional was approximated by the Perdew-Burke-Ernzerhof functional. Dispersion interactions were captured using the D3 method. Projector augmented-wave potentials described electron-ion interactions, and electron wavefunctions were expanded via a planewave basis with an energy cutoff of 500 eV. Activation barriers were calculated using the climbing image nudged elastic band method. Results and Discussion: To discern the role of Pt-metal and Pt-TiO2 interface sites in promoting acetic acid HDO chemistry, Pt(111) slab and anatase TiO2(101)-supported Pt-nanowire (PtNW/OH-TiO2) surface models were constructed, respectively. Because H2 is co-fed in HDO reactions, the anatase support was terminated with OH groups. Interfacial vacancies have been shown to facilitate Ru/TiO2-catalyzed phenol HDO; thus, an interfacial model with an OH vacancy was also considered (PtNW/OHv-TiO2). Pt-TiO2-interface sites stabilize adsorption of all studied acetic acid HDO surface intermediates relative to terrace Pt-metal sites, particularly when an interfacial-OH vacancy is present. Oxygenated species prefer to bind at the OH vacancy through the O atom, suggesting a preference for C-O over C-C bond cleavage at these sites. This hypothesis is supported by net-negative and net-positive average shifts in the reaction energy and activation energy barriers for C-O and C-C bond-breaking steps, respectively, at Pt-TiO2-interface sites relative to Pt-metal sites. Using the calculated energetics, the predicted minimum-energy pathway was determined for each surface. Pt(111) is predicted to follow decarboxylation, producing undesired methane and carbon dioxide. Conversely, PtNW/OH-TiO2 and PtNW/OHv-TiO2 are both predicted to produce desired acetaldehyde and ethane. The interfacial vacancy may also play a key role in facilitating the first C-O bond-breaking step in acetic acid HDO, lowering the barrier by 0.6 eV relative to the defect-free interface model. These results demonstrate the critical role of the Pt-TiO2 interface in the shift in acetic acid HDO selectivity experimentally observed on Pt/C and Pt/TiO2 catalysts. The results herein demonstrate the important role of the Pt-TiO2 interface and interfacial oxygen vacancies in improving the carbon efficiency for HDO reactions in CFP upgrading.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2019

An analytic expression approximating the Debye heat capacity function

It is important to have analytic expressions for critical functions in the equations of state of materials. The Debye model has been very successful in approximating the thermal energy properties of a variety of solids, but is nonanalytic. Existing approximations suffer from various shortcomings, the most common being lack of applicability over some temperature range. A new analytic and integrable functional form that closely approximates the Debye model for the heat capacity is presented. This form, based on the mean of two Einstein heat capacity functions with a low temperature correction, exhibits deviations from the Debye model smaller than typical experimental scatter in heat capacity data.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Generalized Optimized Effective Potential for Orbital Functionals and Self-Consistent Calculation of Random Phase Approximations

A new self-consistent procedure for calculating the total energy with an orbital-dependent density functional approximation (DFA), the generalized optimized effective potential (GOEP), is developed in the present work. The GOEP is a nonlocal Hermitian potential that delivers the sets of occupied and virtual orbitals and minimizes the total energy. The GOEP optimization leads to the same minimum as does the orbital optimization. The GOEP method is promising as an effective optimization approach for orbital-dependent functionals, as demonstrated for the self-consistent calculations of the random phase approximation (RPA) to the correlation functionals in the particle–hole (ph) and particle–particle (pp) channels. The results show that the accuracy in describing the weakly interacting van der Waals systems is significantly improved in the self-consistent calculations. In particular, the important single excitations contribution in non-self-consistent RPA calculations can be captured self-consistently through the GOEP optimization, leading to orbital renormalization, without using the single excitations in the energy functional.

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

Excitation Energies from the Single-Particle Green’s Function with the GW Approximation

Quasi-particle energies are vital to predicting molecular ionization energies and bulk band structures. The state-of-the-art method for quasi-particle energy calculations, particularly for bulk systems, is the GW approximation. For excited state calculations, one needs to go beyond the GW approximation. The Bethe–Salpeter equation (BSE) is the commonly used approach for bulk-system excited state calculations beyond the GW approximation, which is accurate but computationally cumbersome. In this Article, we create a new method to extract excitation energies directly from the quasi-particle energies based on the GW approximation. Starting from the (N – 1)-electron system, we are able to calculate molecular excitation energies with orbital energies at the GW level for HOMO excitations. Our calculations demonstrate that this method can accurately capture low-lying local excitations as well as charge transfer excitations in many molecular systems. This approach is shown to outperform the time-dependent density functional theory (TDDFT) and are comparable with higher level excited state calculations, including the equation-of-motion couple cluster (EOM-CC) theory and the BSE, but with less computational effort. This new approach provides an efficient alternative to the BSE method for accurate excited state calculations.

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

Universal approximation of symmetric and anti-symmetric functions

In this work, we consider universal approximations of symmetric and anti-symmetric functions, which are important for applications in quantum physics, as well as other scientific and engineering computations. We give constructive approximations with explicit bounds on the number of parameters with respect to the dimension and the target accuracy ϵ. While the approximation still suffers from the curse of dimensionality, to the best of our knowledge, these are the first results in the literature with explicit error bounds for functions with symmetry or anti-symmetry constraints

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Computing partition functions in the one-clean-qubit model

We present a method to approximate partition functions of quantum systems using mixed-state quantum computation. For positive-semidefinite Hamiltonians, our method has an expected running-time that is almost linear in [M/( ε rel Z )] 2 , where M is the dimension of the quantum system, Z is the partition function, and ε rel is the relative precision. It is based on approximations of the exponential operator as linear combinations of certain operators related to block-encoding of Hamiltonians or Hamiltonian evolutions. The trace of each operator is estimated using a standard algorithm in the one-clean-qubit model. For large values of Z , our method may run faster than exact classical methods, whose complexities are polynomial in M . We also prove that a version of the partition function estimation problem within additive error is complete for the so-called DQC1 complexity class, suggesting that our method provides a superpolynomial speedup for certain parameter values. Overall, to attain a desired relative precision, we develop a classical procedure based on a sequence of approximations within predetermined additive errors that may be of independent interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Approximating Density Probability Distribution Functions Across Cosmologies

Using a suite of self-similar cosmological simulations, we measure the probability distribution functions (PDFs) of real-space density, redshift-space density, and their geometric mean. We find that the real-space density PDF is well-described by a function of two parameters: ns, the spectral slope, and σL, the linear rms density fluctuation. For redshift-space density and the geometric mean of real- and redshift-space densities, we introduce a third parameter, ${s}_{L}=\sqrt{\langle {\left({{dv}}_{\mathrm{pec}}^{L}/{dr}\right)}^{2}\rangle }/H$. We find that density PDFs for the LCDM cosmology is also well-parameterized by these three parameters. As a result, we are able to use a suite of self-similar cosmological simulations to approximate density PDFs for a range of cosmologies. We make the density PDFs publicly available and provide an analytical fitting formula for them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Non-adiabatic approximations in time-dependent density functional theory: progress and prospects

Time-dependent density functional theory continues to draw a large number of users in a wide range of fields exploring myriad applications involving electronic spectra and dynamics. Although in principle exact, the predictivity of the calculations is limited by the available approximations for the exchange-correlation functional. In particular, it is known that the exact exchange-correlation functional has memory-dependence, but in practise adiabatic approximations are used which ignore this. Here we review the development of non-adiabatic functional approximations, their impact on calculations, and challenges in developing practical and accurate memory-dependent functionals for general purposes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Calculation and interpretation of classical turning surfaces in solids

Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.

Chemistry↗
DOE OSTI2023

Fast approximation to supershell partition functions: Explicit forms of the coefficients

In a previous work, we derived a formula for supershell partition functions, which are the cornerstone of the Super-Transition-Array approach to radiative-opacity calculations. The new expression takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients (denoted Γ $k$ ) for which we obtained a recursion relation. In the present short paper, as a complement of the previous one, we give an explicit formula for the coefficients Γ $k$ . Another recursion relation is also provided, as well as an alternative expression involving Bell polynomials. In conclusion, the connections with cycle indexes of permutation groups and with elementary symmetric polynomials are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Applicability of the thin-film approximation in terahertz photoconductivity measurements

Thin mesoporous photoconductive layers are critically important for efficient water-spitting solar cells. A detailed understanding of photoconductivity in these materials can be achieved via terahertz transient absorption measurements. Such measurements are commonly interpreted using the thin-film approximation. We compare this approximation with a numerical solution of the transfer function without approximations using experimental results for thin-film mesoporous tin oxide (SnO 2 ) samples which range in thickness from 3.3 to 12.6 µm. These samples were sensitized with either a ruthenium polypyridyl complex or a porphyrin dye. The two sensitizers have markedly different absorption coefficients, resulting in penetration depths of 15 µm and 1 µm, respectively. The thin-film approximation results are in good agreement with the numerical work-up for the short penetration length dye. For the longer penetration length samples, the thin-film formula fails even for thicknesses of only 3 µm ≈ λ/100. The imaginary part of the conductivity calculated using the thin-film formula was significantly larger in magnitude than the value without approximations. This discrepancy between the commonly used thin-film approximation and the numerical solution demonstrates the need for a careful analysis of the thin-film formula.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Assessment of exchange-correlation functionals on oxygen vacancy formation energies of metal oxides

Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Natively Periodic Fast Multipole Method: Approximating the Optimal Green’s Function

The fast multipole method (FMM) obeys periodic boundary conditions "natively" if it uses a periodic Green's function for computing the multipole expansion in the interaction zone of each FMM oct-tree node. One can define the "optimal" Green's function for such a method that results in the numerical solution that converges to the equivalent particle-mesh (PM) solution in the limit of sufficiently high order of multipoles. A discrete functional equation for the optimal Green's function can be derived, but is not practically useful as methods for its solution are not known. Instead, this paper presents an approximation for the optimal Green's function that is accurate to better than 10 -3 in ${L}_{\mathrm{MAX}}$ norm and 10 -4 in L 2 norm for practically useful multipole counts. Such an approximately optimal Green's function offers a practical way for implementing the FMM with periodic boundary conditions natively, without the need to compute lattice sums or to rely on hybrid FMM-PM approaches.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗
DOE OSTI2023

Benchmark Data Set of Crystalline Organic Semiconductors

This work reports a Benchmark Data set of Crystalline Organic Semiconductors to test calculations of the structural and electronic properties of these materials in the solid state. The data set contains 67 crystals consisting of mostly rigid molecules with a single dominant conformer, covering the majority of known structural types. The experimental crystal structure is available for the entire data set, whereas zero-temperature unit cell volume can be reliably estimated for a subset of 28 crystals. Using this subset, we benchmark r 2 SCAN-D3 and PBE-D3 density functionals. Then, for the entire data set, we benchmark approximate density functional theory (DFT) methods, including GFN1-xTB and DFTB3(3ob-3-1), with various dispersion corrections against r 2 SCAN-D3. Our results show that r 2 SCAN-D3 geometries are accurate within a few percent, which is comparable to the statistical uncertainty of experimental data at a fixed temperature, but the unit cell volume is systematically underestimated by 2% on average. The several times faster PBE-D3 provides an unbiased estimate of the volume for all systems except for molecules with highly polar bonds, for which the volume is substantially overestimated in correlation with the underestimation of atomic charges. Considered approximate DFT methods are orders of magnitude faster and provide qualitatively correct but overcompressed crystal structures unless the dispersion corrections are fitted by unit cell volume.

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

The role of decomposition reactions in assessing first-principles predictions of solid stability

Abstract The performance of density functional theory approximations for predicting materials thermodynamics is typically assessed by comparing calculated and experimentally determined enthalpies of formation from elemental phases, Δ H f . However, a compound competes thermodynamically with both other compounds and their constituent elemental forms, and thus, the enthalpies of the decomposition reactions to these competing phases, Δ H d , determine thermodynamic stability. We evaluated the phase diagrams for 56,791 compounds to classify decomposition reactions into three types: 1. those that produce elemental phases, 2. those that produce compounds, and 3. those that produce both. This analysis shows that the decomposition into elemental forms is rarely the competing reaction that determines compound stability and that approximately two-thirds of decomposition reactions involve no elemental phases. Using experimentally reported formation enthalpies for 1012 solid compounds, we assess the accuracy of the generalized gradient approximation (GGA) (PBE) and meta-GGA (SCAN) density functionals for predicting compound stability. For 646 decomposition reactions that are not trivially the formation reaction, PBE (mean absolute difference between theory and experiment (MAD) = 70 meV/atom) and SCAN (MAD = 59 meV/atom) perform similarly, and commonly employed correction schemes using fitted elemental reference energies make only a negligible improvement (~2 meV/atom). Furthermore, for 231 reactions involving only compounds (Type 2), the agreement between SCAN, PBE, and experiment is within ~35 meV/atom and is thus comparable to the magnitude of experimental uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A compendium of human gene functions derived from evolutionary modelling

A comprehensive, computable representation of the functional repertoire of all macromolecules encoded within the human genome is a foundational resource for biology and biomedical research. The Gene Ontology Consortium has been working towards this goal by generating a structured body of information about gene functions, which now includes experimental findings reported in more than 175,000 publications for human genes and genes in experimentally tractable model organisms 1,2 . Here, we describe the results of a large, international effort to integrate all of these findings to create a representation of human gene functions that is as complete and accurate as possible. Specifically, we apply an expert-curated, explicit evolutionary modelling approach to all human protein-coding genes. This approach integrates available experimental information across families of related genes into models that reconstruct the gain and loss of functional characteristics over evolutionary time. The models and the resulting set of 68,667 integrated gene functions cover approximately 82% of human protein-coding genes. The functional repertoire reveals a marked preponderance of molecular regulatory functions, and the models provide insights into the evolutionary origins of human gene functions. We show that our set of descriptions of functions can improve the widely used genomic technique of Gene Ontology enrichment analysis. The experimental evidence for each functional characteristic is recorded, thereby enabling the scientific community to help review and improve the resource, which we have made publicly available.

59 BASIC BIOLOGICAL SCIENCES↗