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

Criteria for the Best Approximation by Simple Partial Fractions on Semi-Axis and Axis

We study uniform approximation of real-valued functions f, f(∞) = 0, on ℝ{sup +} and ℝ by real-valued simple partial fractions (the logarithmic derivatives of polynomials). We obtain a criterion for the best approximation on ℝ{sup +} and ℝ in terms of the Chebyshev alternance. This criterion is similar to the known criterion on finite segments. For the problem of approximating odd functions on ℝ we construct an alternance criterion with a weakened condition on the poles of fractions. We present a criterion for the best approximation by simple partial fractions on ℝ{sup +} and ℝ in terms of Kolmogorov. We prove analogs of the de la Vallee-Poussin alternation theorem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Global and local curvature in density functional theory

Piecewise linearity of the energy with respect to fractional electron removal or addition is a requirement of an electronic structure method that necessitates the presence of a derivative discontinuity at integer electron occupation. Semi-local exchange-correlation (xc) approximations within density functional theory (DFT) fail to reproduce this behavior, giving rise to deviations from linearity with a convex global curvature that is evidence of many-electron, self-interaction error and electron delocalization. Popular functional tuning strategies focus on reproducing piecewise linearity, especially to improve predictions of optical properties. In a divergent approach, Hubbard U-augmented DFT (i.e., DFT+U) treats self-interaction errors by reducing the local curvature of the energy with respect to electron removal or addition from one localized subshell to the surrounding system. Although it has been suggested that DFT+U should simultaneously alleviate global and local curvature in the atomic limit, no detailed study on real systems has been carried out to probe the validity of this statement. In this work, we show when DFT+U should minimize deviations from linearity and demonstrate that a “+U” correction will never worsen the deviation from linearity of the underlying xc approximation. However, we explain varying degrees of efficiency of the approach over 27 octahedral transition metal complexes with respect to transition metal (Sc–Cu) and ligand strength (CO, NH{sub 3}, and H{sub 2}O) and investigate select pathological cases where the delocalization error is invisible to DFT+U within an atomic projection framework. Finally, we demonstrate that the global and local curvatures represent different quantities that show opposing behavior with increasing ligand field strength, and we identify where these two may still coincide.

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

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

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

Technical Note: Statistical dependences between channels in radiochromic film readings. Implications in multichannel dosimetry

Purpose: This note studies the statistical relationships between color channels in radiochromic film readings with flatbed scanners. The same relationships are studied for noise. Finally, their implications for multichannel film dosimetry are discussed. Methods: Radiochromic films exposed to wedged fields of 6 MV energy were read in a flatbed scanner. The joint histograms of pairs of color channels were used to obtain the joint and conditional probability density functions between channels. Then, the conditional expectations and variances of one channel given another channel were obtained. Noise was extracted from film readings by means of a multiresolution analysis. Two different dose ranges were analyzed, the first one ranging from 112 to 473 cGy and the second one from 52 to 1290 cGy. Results: For the smallest dose range, the conditional expectations of one channel given another channel can be approximated by linear functions, while the conditional variances are fairly constant. The slopes of the linear relationships between channels can be used to simplify the expression that estimates the dose by means of the multichannel method. The slopes of the linear relationships between each channel and the red one can also be interpreted as weights in the final contribution to dose estimation. However, for the largest dose range, the conditional expectations of one channel given another channel are no longer linear functions. Finally, noises in different channels were found to correlate weakly. Conclusions: Signals present in different channels of radiochromic film readings show a strong statistical dependence. By contrast, noise correlates weakly between channels. For the smallest dose range analyzed, the linear behavior between the conditional expectation of one channel given another channel can be used to simplify calculations in multichannel film dosimetry.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies

Abstract A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗
DOE OSTI2024

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Analysis of over-magnetization of elemental transition metal solids from the SCAN density functional

Current investigations have revealed that the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation exchange-correlation functional significantly over-magnetizes elemental Fe, Co, and Ni solids. For the paradigmatic case, bcc Fe, the error relative to experiment is ≳ 20 % . A comparative analysis of magnetization results from SCAN and its deorbitalized counterpart, SCAN-L, leads to identification of the source of the discrepancy. It is not from the difference between Kohn-Sham (SCAN-L) and generalized Kohn-Sham (SCAN) procedures. The key is the iso-orbital indicator α (the ratio of the local Pauli and Thomas-Fermi kinetic energy densities). Its deorbitalized counterpart α L has more dispersion in both spin channels with respect to magnetization in an approximate region between 0.6 and 1.2 bohrs around an Fe nucleus. The overall effect is that the SCAN switching function evaluated with α L reduces the energetic disadvantage of the down channel with respect to up compared to the original α , which in turn reduces the magnetization. This identifies the origin of the SCAN magnetization error as insensitivity of the SCAN switching function to α values in the approximate range 0.5 ≲ α ≲ 0.8 and oversensitivity for α ≳ 0.8 .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shock-induced metallization of polystyrene along the principal Hugoniot investigated by advanced thermal density functionals

To date, none of the ab-initio molecular dynamics simulations of polystyrene, often used as an ablator material in inertial confinement fusion targets, with the standard ground-state exchange-correlation (XC) functional in density-functional theory can satisfactorily agree with experiments in terms of reflectivity measurements. Here, we use recently developed thermal strongly constrained and appropriately normed Laplacian dependent meta-generalized gradient approximation XC density functional (T-SCAN-L) and thermal hybrid XC density functional (KDT0) to show that the inclusion of thermal and inhomogeneity effects is crucial for accurate prediction of structural evolution and corresponding insulator-metal transition (IMT) during shock compression. Optical reflectivity calculated as an indicator of IMT is in perfect accord with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Time-dependent current-density-functional theory of charge, energy and spin transport and dynamics in nanoscale systems (Final Report)

The focus of this research program is on developing density functional methods for the calculation of non-equilibrium phenomena of the transport type (finite wavelength, low frequency) and of the dynamical type (long wavelength, high frequency). One such method is the time-dependent current density functional theory, in which the evolution of the physical system is mimicked in a fictitious non-interacting system subjected to an effective one-particle potential, which we approximate as a function of the equilibrium density and non-equilibrium current density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Optimal Power Flow of Radial Networks and Its Variations: A Sequential Convex Optimization Approach

This paper proposes a sequential convex optimization method to solve broader classes of optimal power flow (OPF) problems over radial networks. The non-convex branch power flow equation is decomposed as a second-order cone inequality and a non-convex constraint involving the difference of two convex functions. Provided with an initial solution offered by an inexact second-order cone programming relaxation model, this approach solves a sequence of convexified penalization problems, where concave terms are approximated by linear functions and updated in each iteration. It could recover a feasible power flow solution, which usually appears to be very close, if not equal, to the global optimal one. Two variations of the OPF problem, in which non-cost related objectives are optimized subject to power flow constraints and the convex relaxation is generally inexact, are elaborated in detail. One is the maximum loadability problem, which is formulated as a special OPF problem that seeks the maximal distance to the boundary of power flow insolvability. The proposed method is shown to outperform commercial nonlinear solvers in terms of robustness and efficiency. The other is the bi-objective OPF problem. A non-parametric scalarization model is suggested, and is further reformulated as an extended OPF problem by convexifying the objective function. It provides a single trade-off solution without any subjective preference. The proposed computation framework also helps retrieve the Pareto front of the bi-objective OPF via the e-constraint method or the normal boundary intersection method. This paper also discusses extensions for OPF problems over meshed networks based on the semidefinite programming relaxation method.

bi-objective optimization↗
DOE OSTI2018

Vibrational and Thermal Properties of Oxyanionic Crystals

The vibrational and thermal properties of dolomite and alkali chlorates and perchlorates were studied in the gradient approximation of density functional theory using the method of a linear combination of atomic orbitals (LCAO). Long-wave vibration frequencies, IR and Raman spectra, and mode Gruneisen parameters were calculated. Equation-of-state parameters, thermodynamic potentials, entropy, heat capacity, and thermal expansion coefficient were also determined. The thermal expansion coefficient of dolomite was established to be much lower than for chlorates and perchlorates. The temperature dependence of the heat capacity at T > 200 K was shown to be generally governed by intramolecular vibrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Direct computation of parameters for accurate polarizable force fields

We present an improved electronic linear response model to incorporate polarization and charge-transfer effects in polarizable force fields. This model is a generalization of the Atom-Condensed Kohn-Sham Density Functional Theory (DFT), approximated to second order (ACKS2): it can now be defined with any underlying variational theory (next to KS-DFT) and it can include atomic multipoles and off-center basis functions. Parameters in this model are computed efficiently as expectation values of an electronic wavefunction, obviating the need for their calibration, regularization, and manual tuning. In the limit of a complete density and potential basis set in the ACKS2 model, the linear response properties of the underlying theory for a given molecular geometry are reproduced exactly. A numerical validation with a test set of 110 molecules shows that very accurate models can already be obtained with fluctuating charges and dipoles. These features greatly facilitate the development of polarizable force fields.

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

Statistics of dark matter halos in the excursion set peak framework

We derive approximated, yet very accurate analytical expressions for the abundance and clustering properties of dark matter halos in the excursion set peak framework; the latter relies on the standard excursion set approach, but also includes the effects of a realistic filtering of the density field, a mass-dependent threshold for collapse, and the prescription from peak theory that halos tend to form around density maxima. We find that our approximations work excellently for diverse power spectra, collapse thresholds and density filters. Moreover, when adopting a cold dark matter power spectra, a tophat filtering and a mass-dependent collapse threshold (supplemented with conceivable scatter), our approximated halo mass function and halo bias represent very well the outcomes of cosmological N-body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Conditional probability density functional theory

Here we present conditional probability (CP) density functional theory (DFT) as a formally exact theory. In essence, CP-DFT determines the ground-state energy of a system by finding the CP density from a series of independent Kohn-Sham (KS) DFT calculations. By directly calculating CP densities, we bypass the need for an approximate XC energy functional. In this paper we discuss and derive several key properties of the CP density and corresponding CP-KS potential. Illustrative examples are used throughout to help guide the reader through the various concepts and theory presented. We explore a suitable CP-DFT approximation and discuss exact conditions, limitations, and results for selected examples.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Discrete Hahn polynomials for numerical solution of two-dimensional variable-order fractional Rayleigh–Stokes problem

The main aim of this paper is to find the numerical solutions of 2D Rayleigh–Stokes problem with the variable-order fractional derivatives in the Riemann–Liouville sense. The presented method is based on collocation procedure in combination with the new operational matrix of the variable-order fractional derivatives, in the Caputo sense, for the discrete Hahn polynomials. The main advantage of the proposed method is obtaining a global approximation for spatial and temporal discretizations, and it reduced the problem to an algebraic system, which is easier to solve. Also, the profit of approximating a continuous function by Hahn polynomials is that for computing the coefficients of the expansion, we only have to compute a summation and the calculation of coefficients is exact. The error bound for the approximate solution is estimated. Finally, we evaluate results of the presented method with other numerical methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Practical algorithms for multivariate rational approximation

We present two approaches for computing rational approximations to multivariate functions, motivated by their effectiveness as surrogate models for high-energy physics (HEP) applications. Our first approach builds on the Stieltjes process to efficiently and robustly compute the coefficients of the rational approximation. Our second approach is based on an optimization formulation that allows us to include structural constraints on the rational approximation (in particular, constraints demanding the absence of singularities), resulting in a semi-infinite optimization problem that we solve using an outer approximation approach. We present results for synthetic and real-life HEP data, and we compare the approximation quality of our approaches with that of traditional polynomial approximations.

97 MATHEMATICS AND COMPUTING↗