Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional 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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Learning Functions Varying along a Central Subspace

Many functions of interest are in a high-dimensional space but exhibit low-dimensional structures. This paper studies regression of an s-Hölder function in $R^D$ which varies along a central subspace of dimension $d$ while $d \ll D$. A direct approximation of $f$ in $R^D$ with an accuracy $\varepsilon$ requires the number of samples in the order of $\varepsilon^{-(2s+D)/s}$. In this paper, we analyze the generalized contour regression (GCR) algorithm for the estimation of the central subspace and use piecewise polynomials for function approximation. GCR is among the best estimators for the central subspace, but its sample complexity is an open question. In this paper, we partially answer this questions by proving that if a variance quantity is exactly known, GCR leads to a mean squared estimation error of $O(n^{-1})$ for the central subspace. The estimation error of this variance quantity is also given in this paper. The mean squared regression error of $f$ is proved to be in the order of $(n/\log n)^{-\frac{2s}{2s+d}}$, where the exponent depends on the dimension of the central subspace instead of the ambient space . This result demonstrates that GCR is effective in learning the low-dimensional central subspace. We also propose a modified GCR with improved efficiency. Here, the convergence rate is validated through several numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Polarization options in inclusive DIS off tensor polarized deuteron

In the near future, the Jefferson Lab b 1 experiment will provide the second measurement of tensor polarized asymmetries in inclusive DIS on the deuteron. In this asymmetry, 4 independent tensor polarized structure functions contribute. This necessitates systematic approximations in the extraction of the leading twist structure function b 1 from a single tensor asymmetry measurement. Contamination from higher twist structure functions and kinematic effects is discussed here. Using a deuteron convolution model, we quantify the systematic errors from these approximations for two different choices for the target polarization direction (momentum transfer, electron beam direction). For Jefferson Lab 12 GeV kinematics, the systematic error turns out to be comparable between the two polarization options, while at higher Q 2 values the momentum transfer direction is preferred.

Cosyn, Wim [Florida International University, Miam↗
DOE OSTI2018

Electronic optics in graphene in the semiclassical approximation

Highlights: • A rigorous semiclassical theory for the two-dimensional Dirac equation is developed. • We discuss how the use of eikonal coordinates greatly simplifies the theory. • Electronic optics in graphene is developed beyond the standard WKB approximation. • The semiclassical phase essentially affects electron focusing in graphene. - Abstract: We study above-barrier scattering of Dirac electrons by a smooth electrostatic potential combined with a coordinate-dependent mass in graphene. We assume that the potential and mass are sufficiently smooth, so that we can define a small dimensionless semiclassical parameter h≪1. This electronic optics setup naturally leads to focusing and the formation of caustics, which are singularities in the density of trajectories. We construct a semiclassical approximation for the wavefunction in all points, placing particular emphasis on the region near the caustic, where the maximum of the intensity lies. Because of the matrix character of the Dirac equation, this wavefunction contains a nontrivial semiclassical phase, which is absent for a scalar wave equation and which influences the focusing. We carefully discuss the three steps in our semiclassical approach: the adiabatic reduction of the matrix equation to aneffective scalar equation, the construction of the wavefunction using the Maslov canonical operator and the application of the uniform approximation to the integral expression for the wavefunction in the vicinity of a caustic. We consider several numerical examples and show that our semiclassical results are in very good agreement with the results of tight-binding calculations. In particular, we show that the semiclassical phase can have a pronounced effect on the position of the focus and its intensity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Fast and accurate CMB computations in non-flat FLRW universes

We present a new method for calculating CMB anisotropies in a non-flat Friedmann universe, relying on a very stable algorithm for the calculation of hyperspherical Bessel functions, that can be pushed to arbitrary precision levels. We also introduce a new approximation scheme which gradually takes over in the flat space limit and leads to significant reductions of the computation time. Our method is implemented in the Boltzmann code class. It can be used to benchmark the accuracy of the camb code in curved space, which is found to match expectations. For default precision settings, corresponding to 0.1% for scalar temperature spectra and 0.2% for scalar polarisation spectra, our code is two to three times faster, depending on curvature. We also simplify the temperature and polarisation source terms significantly, so the different contributions to the C{sub ℓ} 's are easy to identify inside the code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
DOE OSTI2020

Weighted greedy-optimal design of computer experiments for kernel-based and Gaussian process model emulation and calibration

This article is concerned with the approximation of high-dimensional functions by kernel-based methods. Motivated by uncertainty quantification, which often necessitates the construction of approximations that are accurate with respect to a probability density function of random variables, we aim at minimizing the approximation error with respect to a weighted $L^p$-norm. We present a greedy procedure for designing computer experiments based upon a weighted modification of the pivoted Cholesky factorization. The method successively generates nested samples with the goal of minimizing error in regions of high probability. Numerical experiments validate that this new importance sampling strategy is superior to other sampling approaches, especially when used with non-product probability density functions. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

An Efficient, Globally Convergent Method for Optimization Under Uncertainty Using Adaptive Model Reduction and Sparse Grids

We present a new approach to efficiently solve optimization problems constrained by partial differential equations (PDEs) with uncertain coefficients. The method leverages two sources of inexactness that trade accuracy for speed: (1) stochastic collocation based on dimension-adaptive sparse grids (SGs), which approximates the stochastic objective function with a limited number of quadrature nodes, and (2) projection-based reduced-order models (ROMs), which generate efficient approximations to PDE solutions. These two sources of inexactness lead to inexact objective function and gradient evaluations, which are managed by a trust-region method that guarantees global convergence by adaptively refining the SG and ROM until a proposed error indicator drops below a tolerance specified by trust-region convergence theory. A key feature of the proposed method is that the error indicator---which accounts for errors incurred by both the SG and ROM---must be only an asymptotic error bound, i.e., a bound that holds up to an arbitrary constant that need not be computed. This allows the approach to be applicable to a wide range of problems, including those where sharp, computable error bounds are not available; this distinguishes the proposed method from previous works. Numerical experiments performed on a model problem from optimal flow control under uncertainty verify global convergence of the method and demonstrate the method's ability to outperform previously proposed alternatives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Full self-consistency versus quasiparticle self-consistency in diagrammatic approaches: Exactly solvable two-site Hubbard model

Self-consistent solutions of Hedin's equations (HE) for the two-site Hubbard model (HM) have been studied. They have been found for three-point vertices of increasing complexity (Γ = 1 (GW approximation), Γ₁ from the first-order perturbation theory, and the exact vertex Γ E ). Comparison is made between the cases when an additional quasiparticle (QP) approximation for Green's functions is applied during the self-consistent iterative solving of HE and when QP approximation is not applied. Results obtained with the exact vertex are directly related to the present open question—which approximation is more advantageous for future implementations, GW + DMFT or QPGW + DMFT. It is shown that in a regime of strong correlations only the originally proposed GW + DMFT scheme is able to provide reliable results. Vertex corrections based on Perturbation Theory systematically improve the GW results when full self-consistency is applied. The application of QP self-consistency combined with PT vertex corrections shows similar problems to the case when the exact vertex is applied combined with QP sc. An analysis of Ward Identity violation is performed for all studied in this work's approximations and its relation to the general accuracy of the schemes used is provided.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

How well do self-interaction corrections repair the overestimation of static polarizabilities in density functional calculations?

Herein, we examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys. , 2018, 20 , 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew–Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi–Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys. , 2019, 151 , 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.

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

Min-Max Spaces and Complexity Reduction in Min-Max Expansions

Idempotent methods have been found to be extremely helpful in the numerical solution of certain classes of nonlinear control problems. In those methods, one uses the fact that the value function lies in the space of semiconvex functions (in the case of maximizing controllers), and approximates this value using a truncated max-plus basis expansion. In some classes, the value function is actually convex, and then one specifically approximates with suprema (i.e., max-plus sums) of affine functions. Note that the space of convex functions is a max-plus linear space, or moduloid. In extending those concepts to game problems, one finds a different function space, and different algebra, to be appropriate. Here we consider functions which may be represented using infima (i.e., min-max sums) of max-plus affine functions. It is natural to refer to the class of functions so represented as the min-max linear space (or moduloid) of max-plus hypo-convex functions. We examine this space, the associated notion of duality and min-max basis expansions. In using these methods for solution of control problems, and now games, a critical step is complexity-reduction. In particular, one needs to find reduced-complexity expansions which approximate the function as well as possible. We obtain a solution to this complexity-reduction problem in the case of min-max expansions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A well-scaling natural orbital theory

Here, we introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree–Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the oneparticle density matrix.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2019

Revised M11 Exchange-Correlation Functional for Electronic Excitation Energies and Ground-State Properties

The ability of Kohn-Sham density functional theory (KS-DFT) to accurately predict various types of electronic excitation energies with (necessarily approximate) exchange-correlation functionals faces several challenges. Chief among these is that valence excitations are usually inherently multiconfigurational and therefore best treated by functionals with local exchange, whereas Rydberg and charge transfer excitations are often better treated with nonlocal exchange. The question arises of whether one can optimize a functional such that all three kinds of excitations (valence, Rydberg, and charge transfer – including long-range charge transfer) are treated in a balanced and accurate way. The goal of the present work is to try to answer that question and then to optimize a functional with the best possible balanced behavior. Of the variety of functional types available, we select range-separated hybrid meta functionals because (i) range separation allows the percentage of Hartree–Fock (HF) exchange to change with interelectronic separation, and therefore, one can have 100% HF exchange at large interelectronic separations, which gives good performance for long-range charge-transfer excitations, while the range separation allows one to simultaneously have smaller values of HF exchange at small and intermediate inter-electronic separations, which give good performance for valence and Rydberg excitations and (ii) meta functionals allow one to obtain better accuracy with high HF exchange than is possible with functionals whose local part depends only on spin densities and their gradients. Furthermore, this work starts with the range-separated hybrid meta functional, M11, and re-optimizes it (with stronger smoothness restraints) against electronic excitation energies and ground-state properties to obtain a new functional called revM11 that gives good performance for all three types of electronic excitations and at the same time gives very good predictions across-the-board for ground-state properties.

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

Data driven analytic continuation for one particle spectral functions

In this proceeding, an idea is outlined suggesting a generic treatment on any type of input data for a numerical analytic continuation problem, which is needed when dynamical information is to be extracted from a calculationally convenient one particle imaginary time Green function. The quality of the resulting spectral function will rely only on the data to be treated, viz, data-driven. This is different from the Maximum Entropy or the Stochastic method which relies on an entropy term to guide convergence of the resulting spectral function.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Two-pion femtoscopic correlations in Be+Be collisions at $\sqrt{s_{\text {NN}}} = 16.84$ GeV measured by the NA61/SHINE at CERN

This paper reports measurements of two-pion femtoscopic correlations in Be+Be collisions at a beam momentum of 150$A\,\hbox {GeV}\!/\!c$ (energy available in the center-of-mass system for nucleon pair $\sqrt{s_{\text {NN}}} = 16.84$ GeV) by the NA61/SHINE experiment at the CERN SPS accelerator. The obtained momentum space correlation functions can be well described by a Lévy distributed source model. The transverse mass dependence of the Lévy source parameters is presented, and their possible theoretical interpretations are discussed. The results show that the Lévy exponent $\alpha $ is approximately constant as a function of $m_{\text {T}}$ , and far from both the Gaussian case of $\alpha = 2$ or the conjectured value at the critical endpoint, $\alpha = 0.5$. The radius scale parameter R shows a slight decrease in $m_{\text {T}}$ , which can be explained as a signature of transverse flow. Finally, an approximately constant trend of the intercept parameter $\lambda $ as a function of $m_{\text {T}}$ was observed, similar to previous NA44 S + Pb results (obtained with a Gaussian approximation, but unlike RHIC results).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Exchange interaction between the triplet exciton and the localized spin in copper-phthalocyanine

Triplet excitonic state in the organic molecule may arise from a singlet excitation and the following inter-system crossing. Especially for a spin-bearing molecule, an exchange interaction between the triplet exciton and the original spin on the molecule can be expected. In this paper, such exchange interaction in copper-phthalocyanine (CuPc, spin-1/2 ) was investigated from first-principles by using density-functional theory within a variety of approximations to the exchange correlation, ranging from local-density approximation to long-range corrected hybrid-exchange functional. The magnitude of the computed exchange interaction is in the order of meV with the minimum value (1.5 meV, ferromagnetic) given by the long-range corrected hybrid-exchange functional CAM-B3LYP. This exchange interaction can therefore give rise to a spin coherence with an oscillation period in the order of picoseconds, which is much shorter than the triplet lifetime in CuPc (typically tens of nanoseconds). This implies that it might be possible to manipulate the localized spin on Cu experimentally using optical excitation and inter-system crossing well before the triplet state disappears.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

CrI 3 revisited with a many-body ab initio theoretical approach

CrI 3 has recently been shown to exhibit low-dimensional, long-range magnetic ordering from few layers to single layers of CrI 3 . The properties of CrI 3 bulk and few layered systems are uniquely defined by a combination of short-range intralayer and long-range interlayer interactions; including strong correlations, exchange and spin orbit coupling. Unfortunately, both the long-range van der Waals interactions, which are driven by dynamic, many-body electronic correlations and the competing strong intralayer correlations, present a formidable challenging for the local or semi local mean-field approximations employed in workhorse electronic structure approaches like density functional theory. In this work, we employ a sophisticated many-body approach that can simultaneously describe long- and short-range correlations. Herein, we establish that the fixed-node diffusion Monte Carlo (FNDMC) method reproduces the experimental interlayer separation distance of bulk CrI 3 for the high-temperature, monoclinic phase, with a reliable prediction of the interlayer binding energy. We subsequently employed, the FNDMC results to benchmark the accuracy of several density functional theory exchange-correlation approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Numerical multi-scaling method to solve the linear stochastic partial differential equations

In this paper, we propose a multi-scaling method to solve a class of the stochastic partial differential equations (SPDEs) of Itô-type. By this method, we approximate the solution of the SPDEs with strongly elliptic operator on a finite domain. The method is based on B-spline wavelet approximations; some of these functions are reshaped to satisfy on boundary conditions exactly. We prove the consistency and stability of the method. Illustrative examples are included to demonstrate the validity and applicability of our proposed technique.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Supersymmetric Dynamics and Zeta-Functions

Boson, fermion, and super oscillators and (statistical) mechanism of cosmological constant; finite approximation of the zeta-function and fermion factorization of the bosonic statistical sum concidered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗