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

Artificial intelligence “sees” split electrons

Chemical bonds between atoms are stabilized by the exchange-correlation (xc) energy, a quantum-mechanical effect in which “social distancing” by electrons lowers their electrostatic repulsion energy. Kohn-Sham density functional theory (DFT) states that the electron density determines this xc energy, but the density functional must be approximated. Furthermore, this is usually done by satisfying exact constraints of the exact functional (making the approximation predictive), by fitting to data (making it interpolative), or both. Two exact constraints—the ensemble-based piecewise linear variation of the total energy with respect to fractional electron number and fractional electron z-component of spin —require hard-to-control nonlocality. On page 1385 of this issue, Kirkpatrick et al. have taken a big step toward more accurate predictions for chemistry through the machine learning of molecular data plus the fractional charge and spin constraints, expressed as data that a machine can learn.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H 2 + Cu(111) and H 2 + Ag(111)

Specific reaction parameter density functionals (SRP-DFs), which can describe dissociative chemisorption reactions on metals to within chemical accuracy, have so far been based on exchange functionals within the generalized gradient approximation (GGA) and on GGA correlation functionals or van der Waals correlation functionals. These functionals are capable of describing the molecule–metal surface interaction accurately, but they suffer from the general GGA problem that this can be done only at the cost of a rather poor description of the metal. Here, we show that it is possible also to construct SRP-DFs for H 2 dissociation on Cu(111) based on meta-GGA functionals, introducing three new functionals based on the “made-simple” (MS) concept. The exchange parts of the three functionals (MS-PBEl, MS-B86bl, and MS-RPBEl) are based on the expressions for the PBE, B86b, and RPBE exchange functionals. Quasi-classical trajectory (QCT) calculations performed with potential energy surfaces (PESs) obtained with the three MS functionals reproduce molecular beam experiments on H 2 , D 2 + Cu(111) with chemical accuracy. Therefore, these three non-empirical functionals themselves are also capable of describing H2 dissociation on Cu(111) with chemical accuracy. Similarly, QCT calculations performed on the MS-PBEl and MS-B86bl PESs reproduced molecular beam and associative desorption experiments on D 2 , H 2 + Ag(111) more accurately than was possible with the SRP48 density functional for H 2 + Cu(111). Also, the three new MS functionals describe the Cu, Ag, Au, and Pt metals more accurately than the all-purpose Perdew–Burke–Ernzerhof (PBE) functional. The only disadvantage we noted of the new MS functionals is that, as found for the example of H 2 + Cu(111), the reaction barrier height obtained by taking weighted averages of the MS-PBEl and MS-RPBEl functionals is tunable over a smaller range (9 kJ/mol) than possible with the standard GGA PBE and RPBE functionals (33 kJ/mol). As a result of this restricted tunability, it is not possible to construct an SRP-DF for H 2 + Ag(111) on the basis of the three examined MS meta-GGA functionals.

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

Analyzing the errors of DFT approximations for compressed water systems

We report an extensive study of the errors of density functional theory (DFT) approximations for compressed water systems. The approximations studied are based on the widely used PBE and BLYP exchange-correlation functionals, and we characterize their errors before and after correction for 1- and 2-body errors, the corrections being performed using the methods of Gaussian approximation potentials. The errors of the uncorrected and corrected approximations are investigated for two related types of water system: first, the compressed liquid at temperature 420 K and density 1.245 g/cm{sup 3} where the experimental pressure is 15 kilobars; second, thermal samples of compressed water clusters from the trimer to the 27-mer. For the liquid, we report four first-principles molecular dynamics simulations, two generated with the uncorrected PBE and BLYP approximations and a further two with their 1- and 2-body corrected counterparts. The errors of the simulations are characterized by comparing with experimental data for the pressure, with neutron-diffraction data for the three radial distribution functions, and with quantum Monte Carlo (QMC) benchmarks for the energies of sets of configurations of the liquid in periodic boundary conditions. The DFT errors of the configuration samples of compressed water clusters are computed using QMC benchmarks. We find that the 2-body and beyond-2-body errors in the liquid are closely related to similar errors exhibited by the clusters. For both the liquid and the clusters, beyond-2-body errors of DFT make a substantial contribution to the overall errors, so that correction for 1- and 2-body errors does not suffice to give a satisfactory description. For BLYP, a recent representation of 3-body energies due to Medders, Babin, and Paesani [J. Chem. Theory Comput. 9, 1103 (2013)] gives a reasonably good way of correcting for beyond-2-body errors, after which the remaining errors are typically 0.5 mE{sub h} ≃ 15 meV/monomer for the liquid and the clusters.

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

A simple generalized gradient approximation for the noninteracting kinetic energy density functional

A simple, novel, non-empirical, constraint-based orbital-free generalized gradient approximation (GGA) non-interacting kinetic energy density functional is presented along with illustrative applications. The innovation is adaptation of constraint-based construction to the essential properties of pseudo-densities from the pseudo-potentials that are essential in plane-wave-basis ab initio molecular dynamics. This contrasts with constraining to the qualitatively different Kato-cusp-condition densities. The single parameter in the new functional is calibrated by satisfying Pauli potential positivity constraints for pseudo-atom densities. Finally, in static lattice tests on simple metals and semiconductors, the new LKT functional outperforms the previous best constraint-based GGA functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), is generally superior to a recently proposed meta-GGA, is reasonably competitive with parametrized two-point functionals, and is substantially faster.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Robustness analysis of an air heating plant and control law by using polynomial chaos

This paper presents a robustness analysis of an air heating plant with a multivariable closed-loop control law by using the polynomial chaos methodology (MPC). The plant consists of a PVC tube with a fan in the air input (that forces the air through the tube) and a mass flux sensor in the output. A heating resistance warms the air as it flows inside the tube, and a thermo-couple sensor measures the air temperature. The plant has thus two inputs (the fan's rotation intensity and heat generated by the resistance, both measured in percent of the maximum value) and two outputs (air temperature and air mass flux, also in percent of the maximal value). The mathematical model is obtained by System Identification techniques. The mass flux sensor, which is nonlinear, is linearized and the delays in the transfer functions are properly approximated by non-minimum phase transfer functions. The resulting model is transformed to a state-space model, which is used for control design purposes. The multivariable robust control design techniques used is the LQG/LTR, and the controllers are validated in simulation software and in the real plant. Finally, the MPC is applied by considering some of the system's parameters as random variables (one at a time, and the system's stochastic differential equations are solved by expanding the solution (a stochastic process) in an orthogonal basis of polynomial functions of the basic random variables. This method transforms the stochastic equations in a set of deterministic differential equations, which can be solved by traditional numerical methods (That is the MPC). Statistical data for the system (like expected values and variances) are then calculated. The effects of randomness in the parameters are evaluated in the open-loop and closed-loop pole's positions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Two-component hybrid time-dependent density functional theory within the Tamm-Dancoff approximation

We report the implementation of a two-component variant of time-dependent density functional theory (TDDFT) for hybrid functionals that accounts for spin-orbit effects within the Tamm-Dancoff approximation (TDA) for closed-shell systems. The influence of the admixture of Hartree-Fock exchange on excitation energies is investigated for several atoms and diatomic molecules by comparison to numbers for pure density functionals obtained previously [M. Kühn and F. Weigend, J. Chem. Theory Comput. 9, 5341 (2013)]. It is further related to changes upon switching to the local density approximation or using the full TDDFT formalism instead of TDA. Efficiency is demonstrated for a comparably large system, Ir(ppy){sub 3} (61 atoms, 1501 basis functions, lowest 10 excited states), which is a prototype molecule for organic light-emitting diodes, due to its “spin-forbidden” triplet-singlet transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics↗
DOE OSTI2016

Density functional theory for d- and f-electron materials and compounds

Here, the fundamental requirements for a computationally tractable Density Functional Theory-based method for relativistic f- and (nonrelativistic) d-electron materials and compounds are presented. The need for basing the Kohn–Sham equations on the Dirac equation is discussed. The full Dirac scheme needs exchange-correlation functionals in terms of four-currents, but ordinary functionals, using charge density and spin-magnetization, can be used in an approximate Dirac treatment. The construction of a functional that includes the additional confinement physics needed for these materials is illustrated using the subsystem-functional scheme. If future studies show that a full Dirac, four-current based, exchange-correlation functional is needed, the subsystem functional scheme is one of the few schemes that can still be used for constructing functional approximations.

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

Multiconfiguration Pair-Density Functional Theory: A Fully Translated Gradient Approximation and Its Performance for Transition Metal Dimers and the Spectroscopy of Re 2 Cl 8 2–

We extend the on-top density functional of multiconfiguration pair-density functional theory (MC-PDFT) to include the gradient of the on-top density as well as the gradient of the density. We find that the theory is reasonably stable to this extension; furthermore, it provides improved accuracy for molecules containing transition metals. We illustrate the extended on-top density functionals by applying them to Cr 2 , Cu 2 , Ag 2 , Os 2 , and Re 2 Cl 8 2– as well as to our previous database of 56 data for bond dissociation energies, barrier heights, reaction energies, proton affinities, and the water dimer. The performance of MC-PDFT is comparable to or better than that of CASPT2.

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

Increasing dose gradient and uniformity in small fields using modulation: Theory and prototypes for cone-based stereotactic radiosurgery

Purpose: To investigate the theoretical limits to the tradeoff between dose gradient and uniformity when modulation is used in the context of cone based SRS, and to design a prototype collimation device that allows for steeper dose gradients and/or higher target uniformity as compared to a standard circular collimator. Methods: An inverse planning optimization is performed in the context of idealized phantom geometry to determine the ideal fluence pattern that best approximates a “rect function” dose distribution. Ideal fluence patterns were approximated in a prototype device and radiochromic film dosimetry was utilized to compare the prototype device to a standard circular collimator. Results: For choices of prescription isodose lines above approximately 50%, utilizing modulation allows for an improved tradeoff between dose gradient index and dose heterogeneity index. Compensators placed within the circular collimator can achieve the necessary modulation. Conclusions: Using modulation with features on a submillimeter distance scale, it is possible to increase the dose gradient and/or uniformity in small fields.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

WE-DE-201-01: BEST IN PHYSICS (THERAPY): A Fast Multi-Target Inverse Treatment Planning Strategy Optimizing Dosimetric Measures for High-Dose-Rate (HDR) Brachytherapy

Purpose: Inverse treatment planning (ITP) for interstitial HDR brachytherapy of gynecologic cancers seeks to maximize coverage of the clinical target volumes (tumor and vagina) while respecting dose-volume-histogram related dosimetric measures (DMs) for organs at risk (OARs). Commercially available ITP tools do not support DM-based planning because it is computationally too expensive to solve. In this study we present a novel approach that allows fast ITP for gynecologic cancers based on DMs for the first time. Methods: This novel strategy is an optimization model based on a smooth DM-based objective function. The smooth approximation is achieved by utilizing a logistic function for the evaluation of DMs. The resulting nonconvex and constrained optimization problem is then optimized with a BFGS algorithm. The model was evaluated using the implant geometry extracted from 20 patient treatment plans under an IRB-approved retrospective study. For each plan, the final DMs were evaluated and compared to the original clinical plans. The CTVs were the contoured tumor volume and the contoured surface of the vagina. Statistical significance was evaluated with a one-sided paired Wilcoxon signed-rank test. Results: As did the clinical plans, all generated plans fulfilled the defined DMs for OARs. The proposed strategy showed a statistically significant improvement (p<0.001) in coverage of the tumor and vagina, with absolute improvements of related DMs of (6.9 +/− 7.9)% and (28.2 +/− 12.0)%, respectively. This was achieved with a statistically significant (p<0.01) decrease of the high-dose-related DM for the tumor. The runtime of the optimization was (2.3 +/− 2.0) seconds. Conclusion: We demonstrated using clinical data that our novel approach allows rapid DM-based optimization with improved coverage of CTVs with fewer hot spots. Being up to three orders of magnitude faster than the current clinical practice, the method dramatically shortens planning time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Quantifying the hydration structure of sodium and potassium ions: taking additional steps on Jacob's Ladder

The ability to reproduce the experimental structure of water around the sodium and potassium ions is a key test of the quality of interaction potentials due to the central importance of these ions in a wide range of important phenomena. In this study, we simulate the Na + and K + ions in bulk water using three density functional theory functionals: (1) the generalized gradient approximation (GGA) based dispersion corrected revised Perdew, Burke, and Ernzerhof functional (revPBE-D3) (2) the recently developed strongly constrained and appropriately normed (SCAN) functional (3) the random phase approximation (RPA) functional for potassium. We compare with experimental X-ray diffraction (XRD) and X-ray absorption fine structure (EXAFS) measurements to demonstrate that SCAN accurately reproduces key structural details of the hydration structure around the sodium and potassium cations, whereas revPBE-D3 fails to do so. However, we show that SCAN provides a worse description of pure water in comparison with revPBE-D3. RPA also shows an improvement for K + , but slow convergence prevents rigorous comparison. Finally, we analyse cluster energetics to show SCAN and RPA have smaller fluctuations of the mean error of ion–water cluster binding energies compared with revPBE-D3.

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

Approximate two-body generating Hamiltonian for the particle-hole Pfaffian wave function

We present two two-body Hamiltonians that approximate the exact particle-hole Pfaffian wave function with their ground states for all the system sizes where this wave function has been numerically constructed to date. The approximate wave functions have high overlap with the original and reproduce well the low-lying entanglement spectrum and structure factor. The approximate generating Hamiltonians are obtained by an optimization procedure where three to four pseudopotentials are varied in the neighbourhood of second Landau level Coulomb interaction or of a noninteracting model. They belong to a finite region in the variational space of Hamiltonians where each point approximately generates the particle-hole Pfaffian. Here we diagonalize the identified Hamiltonians for up to 20 electrons and find that for them the particle-hole Pfaffian shift appears energetically more favorable. The possibility to interpret the data in terms of composite fermions is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Chaos control of the brushless direct current motor using adaptive dynamic surface control based on neural network with the minimum weights

This paper investigates chaos control for the brushless DC motor (BLDCM) system by adaptive dynamic surface approach based on neural network with the minimum weights. The BLDCM system contains parameter perturbation, chaotic behavior, and uncertainty. With the help of radial basis function (RBF) neural network to approximate the unknown nonlinear functions, the adaptive law is established to overcome uncertainty of the control gain. By introducing the RBF neural network and adaptive technology into the dynamic surface control design, a robust chaos control scheme is developed. It is proved that the proposed control approach can guarantee that all signals in the closed-loop system are globally uniformly bounded, and the tracking error converges to a small neighborhood of the origin. Simulation results are provided to show that the proposed approach works well in suppressing chaos and parameter perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Charge localization and ordering in A 2 Mn 8 O 16 hollandite group oxides: Impact of density functional theory approaches

The phases of A 2 Mn 8 O 16 hollandite group oxides emerge from the competition between ionic interactions, Jahn-Teller effects, charge ordering, and magnetic interactions. Their balanced treatment with feasible computational approaches can be challenging for commonly used approximations in density functional theory. Three examples (A = Ag, Li, and K) are studied with a sequence of different approximate exchange-correlation functionals. Starting from a generalized gradient approximation (GGA), an extension to include van der Waals interactions and a recently proposed meta-GGA are considered. Then local Coulomb interactions for the Mn 3d electrons are more explicitly considered with the DFT + U approach. Finally, selected results from a hybrid functional approach provide a reference. Results for the binding energy of the A species in the parent oxide highlight the role of van der Waals interactions. Relatively accurate results for insertion energies can be achieved with a low-U and a high-U approach. In the low-U case, the materials are described as band metals with a high-symmetry, tetragonal crystal structure. In the high-U case, the electrons donated by A result in formation of local Mn 3+ centers and corresponding Jahn-Teller distortions characterized by a local order parameter. The resulting degree of monoclinic distortion depends on charge ordering and magnetic interactions in the phase formed. The reference hybrid functional results show charge localization and ordering. Comparison to low-temperature experiments of related compounds suggests that charge localization is the physically correct result for the hollandite group oxides studied here. Lastly, while competing effects in the local magnetic coupling are subtle, the fully anisotropic implementation of DFT + U gives the best overall agreement with results from the hybrid functional.

meta-GGA↗
DOE OSTI2021

Optical response and band structure of LiCoO 2 including electron-hole interaction effects

The optical response functions and band structures of LiCoO 2 are studied at different levels of approximation, from density functional theory (DFT) in the generalized gradient approximation (GGA) to quasiparticle self-consistent QSGW (with G for Green's function and W for screened Coulomb interaction) without and with ladder diagrams ($QSG\hat{W}$) and the Bethe Salpeter Equation (BSE) approach. The QSGW method is found to strongly overestimate the band gap and electron-hole or excitonic effects are found to be important. They lower the quasiparticle gap by only about 11~\% but the lowest energy peaks in absorption are found to be excitonic in nature. The contributions from different band to band transitions and the relation of excitons to band-to-band transitions are analyzed. Finally, the excitons are found to be strongly localized. Finally, a comparison to experimental data is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A continuous analogue of the tensor-train decomposition

We develop new approximation algorithms and data structures for representing and computing with multivariate functions using the functional tensor-train (FT), a continuous extension of the tensor-train (TT) decomposition. The FT represents functions using a tensor-train ansatz by replacing the three-dimensional TT cores with univariate matrix-valued functions. The main contribution of this paper is a framework to compute the FT that employs adaptive approximations of univariate fibers, and that is not connected to any tensorized discretization. The algorithm can be coupled with any univariate linear or nonlinear approximation procedure. We demonstrate that this approach can generate multivariate function approximations that are several orders of magnitude more accurate, for the same cost, than those based on the conventional approach of compressing the coefficient tensor of a tensor-product basis. Our approach is in the spirit of other continuous computation packages such as Chebfun, and yields an algorithm which requires the computation of “continuous” matrix factorizations such as the LU and QR decompositions of vector-valued functions. To support these developments, we describe continuous versions of an approximate maximum-volume cross approximation algorithm and of a rounding algorithm that re-approximates an FT by one of lower ranks. We demonstrate that our technique improves accuracy and robustness, compared to TT and quantics-TT approaches with fixed parameterizations, of high-dimensional integration, differentiation, and approximation of functions with local features such as discontinuities and other nonlinearities.

42 ENGINEERING↗
DOE OSTI2015

Systematic approach for simultaneously correcting the band-gap and p - d separation errors of common cation III-V or II-VI binaries in density functional theory calculations within a local density approximation

We propose a systematic approach that can empirically correct three major errors typically found in a density functional theory (DFT) calculation within the local density approximation (LDA) simultaneously for a set of common cation binary semiconductors, such as III-V compounds, (Ga or In)X with X = N,P,As,Sb, and II-VI compounds, (Zn or Cd)X, with X = O,S,Se,Te. By correcting (1) the binary band gaps at high-symmetry points , L, X, (2) the separation of p-and d-orbital-derived valence bands, and (3) conduction band effective masses to experimental values and doing so simultaneously for common cation binaries, the resulting DFT-LDA-based quasi-first-principles method can be used to predict the electronic structure of complex materials involving multiple binaries with comparable accuracy but much less computational cost than a GW level theory. This approach provides an efficient way to evaluate the electronic structures and other material properties of complex systems, much needed for material discovery and design.

Wang, Jianwei↗