Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning

Density functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we, in this study, use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO 3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ~1.5 eV for this class of materials, where ~1 eV comes from downward shifting the valence band and ~0.5 eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Assessing the Performance of Exchange‐Correlation Density Functionals in Describing the Iron‐Catalyzed Ammonia Synthesis System

Density functional theory (DFT) has been widely employed for elucidating mechanistic aspects of heterogeneous catalysis. However, the accuracy of DFT calculations relies heavily on selecting an exchange-correlation (XC) functional that correctly describes the electronic structure of materials involved in the reactions. This study assesses the accuracy of several XC density functionals for modeling the iron-catalyzed ammonia synthesis system. In the assessment of functional accuracy, experimental references are compared to DFT-calculated values for the formation energy of gas-phase ammonia and nitrogen, bulk Fe/Fe-nitride (γ′-Fe 4 N) lattice constants and cohesive/formation energies, and nitrogen and ammonia binding energies on Fe(100), Fe(111), Fe(110), and γ′-Fe 4 N(111). It is observed that the experimental value for each of these descriptors is accurately modeled by at least one functional. RPBE alone provides reliable estimates for both the lattice constant and cohesive energy of Fe and γ′-Fe 4 N. Temperature-programmed desorption experiments led to estimates for N and NH 3 adsorption across several Fe-based facets that are best captured by RPBE. These results highlight the importance of choosing an appropriate XC functional that accurately describes experimental systems and offer insights into effectively modeling the interactions between nitrogen and ammonia on Fe-based surfaces.

adsorption↗
DOE OSTI2024

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗
DOE OSTI2023

Assessing the source of error in the Thomas–Fermi–von Weizsäcker density functional

We investigate the source of error in the Thomas–Fermi–von Weizsäcker (TFW) density functional relative to Kohn–Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn–Sham density, the corresponding energy deviates significantly from the Kohn–Sham value. Here, we show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn–Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn–Sham solution.

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

Some Problems in Density Functional Theory

Though calculations based on density functional theory (DFT) are used remarkably widely in chemistry, physics, materials science, and biomolecular research and though the modern form of DFT has been studied for almost 60 years, some mathematical problems remain. From a physical science perspective, it is far from clear whether those problems are of major import. For context, we provide an outline of the basic structure of DFT as it is presented and used conventionally in physical sciences, note some unresolved mathematical difficulties with those conventional demonstrations, then pose several questions regarding both the time-independent and time-dependent forms of DFT that could benefit from attention in applied mathematics. Furthermore, progress on any of these would aid in development of better approximate functionals and in interpretation of DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

CpFe(CO) 2 Radical Generated from Dinuclear [CpFe(CO) 2 ] 2 and Mononuclear (Cp)(CO) 2 Fe(H): Density Functional Theory Is Accurate for One, But Not Both

Density functional theory (DFT) methods remain the most practical approach to calculating properties and reaction mechanisms of transition metal complexes. While the accuracy of DFT methods has been evaluated for some properties of mononuclear organometallic complexes there has been a general lack of evaluation for dinuclear organometallic complexes, in particular bonding changes related to reaction mechanisms. Here, this work evaluated DFT and coupled cluster methods for the accuracy of calculating the CpFe(CO) 2 radical (Fp•) generated from dinuclear [CpFe(CO) 2 ] 2 (Fp 2 ) and mononuclear [(Cp)(CO) 2 Fe(H)] (Fp-H). This transition metal radical fragment was evaluated because dinuclear complexes built with it have recently shown a variety of unique reactions but has proven challenging to accurately calculate with DFT methods. Here we show that DFT methods provide a surprising wide range of fragmentation energies for Fp 2 and lower and mid rung DFT methods as well as DLPNO–CCSD(T) perform well for this dissociation energy. The highest rung double-hybrid methods have a large range in the Fp 2 dissociation energy, and the energy greatly depends on the amount of MP2 correlation energy included. For generating Fp• from Fp-H the lower and mid rung methods that worked well for Fp 2 showed significant error. Double-hybrid methods unfortunately are only accurate for the Fe–H bond if they are very inaccurate for the Fp 2 dissociation energy. While DLPNO–CCSD(T) is not perfect, and not close to chemically accurate for the Fe–H bond, it does provide reasonable accuracy for both Fp 2 and Fp-H dissociation energies.

density functional theory↗
DOE OSTI2022

Double and Charge-Transfer Excitations in Time-Dependent Density Functional Theory

Time-dependent density functional theory has emerged as a method of choice for calculations of spectra and response properties in physics, chemistry, and biology, with its system-size scaling enabling computations on systems much larger than otherwise possible. While increasingly complex and interesting systems have been successfully tackled with relatively simple functional approximations, there has also been increasing awareness that these functionals tend to fail for certain classes of approximations. Here I review the fundamental challenges the approximate functionals have in describing double excitations and charge-transfer excitations, which are two of the most common impediments for the theory to be applied in a black-box way. At the same time, I describe the progress made in recent decades in developing functional approximations that give useful predictions for these excitations.

Chemistry↗
NASA NTRS2012

Ensemble Averaged Probability Density Function (APDF) for Compressible Turbulent Reacting Flows

In this paper, we present a concept of the averaged probability density function (APDF) for studying compressible turbulent reacting flows. The APDF is defined as an ensemble average of the fine grained probability density function (FG-PDF) with a mass density weighting. It can be used to exactly deduce the mass density weighted, ensemble averaged turbulent mean variables. The transport equation for APDF can be derived in two ways. One is the traditional way that starts from the transport equation of FG-PDF, in which the compressible Navier- Stokes equations are embedded. The resulting transport equation of APDF is then in a traditional form that contains conditional means of all terms from the right hand side of the Navier-Stokes equations except for the chemical reaction term. These conditional means are new unknown quantities that need to be modeled. Another way of deriving the transport equation of APDF is to start directly from the ensemble averaged Navier-Stokes equations. The resulting transport equation of APDF derived from this approach appears in a closed form without any need for additional modeling. The methodology of ensemble averaging presented in this paper can be extended to other averaging procedures: for example, the Reynolds time averaging for statistically steady flow and the Reynolds spatial averaging for statistically homogeneous flow. It can also be extended to a time or spatial filtering procedure to construct the filtered density function (FDF) for the large eddy simulation (LES) of compressible turbulent reacting flows.

Shih, Tsan-Hsing↗
DOE OSTI2024

Density functional theory of material design: fundamentals and applications—II

Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.

36 MATERIALS SCIENCE↗
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 OSTI2022

Workhorse minimally empirical dispersion-corrected density functional with tests for weakly bound systems: r 2 SCAN + rVV 10

SCAN+rVV10 has been demonstrated to be a versatile van der Waals (vdW) density functional that delivers good predictions of both energetic and structural properties for many types of bonding. Recently, the r 2 SCAN functional was devised as a revised form of SCAN with improved numerical stability. In this work, we refit the rVV10 functional to optimize the r 2 SCAN+rVV10 vdW density functional and test its performance for molecular interactions and layered materials. Our molecular tests demonstrate that r 2 SCAN+rVV10 outperforms its predecessor SCAN+rVV10 in both efficiency (numerical stability) and accuracy. This good performance is also found in lattice-constant predictions. In comparison with benchmark results from higher-level theories or experiments, r 2 SCAN+rVV10 yields excellent interlayer binding energies and phonon dispersions for layered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Meta-GGA exchange-correlation free energy density functional to increase the accuracy of warm dense matter simulations

We discuss strategies for thermalization of the ground-state meta-generalized gradient approximation (meta-GGA) exchange-correlation (XC) functionals. A simple but accurate scheme is implemented via universal additive thermal correction to XC using a perturbative-like self-consistent approach. The additive correction with explicit temperature dependence is applied to the ground-state deorbitalized, strongly constrained and appropriately normed (SCAN-L) meta-GGA XC leading to thermal XC functional denoted here as T-SCAN-L. Thermal T-SCAN-L meta-GGA functional shows significant improvement in density functional theory calculation accuracy for warm dense matter by a factor of 3 to 10, achieving unprecedented accuracy of total pressure between a few tenths and 1% when compared to traditional XC functionals, as demonstrated by the comparison to pathintegral Monte Carlo simulations for helium equation of state. Furthermore, the T-SCAN-L calculations of dc conductivity of warm dense aluminum also give better agreement with experiments over other XC functionals such as PBE and SCAN-L.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Binding of radionuclides and surrogate to 18-crown-6 ether by density functional theory

For this work, we use density functional theory to investigate the interactions of cerium, americium, and curium cations with crown ethers. Our calculations reveal that the modeled structure of cerium integrated within the crown ether is in good agreement with experimental data, with the negative binding energy indicating that capturing the cerium nitrates is thermodynamically favorable. Our results demonstrate that crown ethers can also bind americium and curium, providing insights into the potential applications of crown ether in radionuclide sequestration. Finally, we explore the impact of the skeleton modification of different crown ethers through by substitution of nitrogen atoms in the core of the crown ether for oxygen atoms and find that this structural modification significantly increases the radionuclide binding energies. These findings provide insights on the potential for the use of organic linkers such as crown ethers to address the urgent needs in radionuclide sequestration, separation and sensing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Field-induced magnetization of δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of magnetization, induced by an external magnetic field, for δ-phase plutonium. The fully relativistic electronic structure accounts for Zeeman splitting effects through a Hamiltonian that couples the magnetic field to both spin and orbital magnetic moments. The electronic-structure model is further improved by an extension to DFT in terms of an orbital-orbital coupling via the conventional orbital-polarization method, as has routinely been done for plutonium. The response to the applied magnetic field is shown to be weak in the DFT model, with induced magnetic moments of the order of in magnetic fields up to 30 T. These results are in accord with recent assessments from x-ray magnetic circular dichroism measurements in magnetic field on δ-plutonium. Somewhat surprisingly, a model assuming δ-plutonium to be absent magnetic moments shows stronger response to an external static magnetic field than models allowing for antiferromagnetism or magnetic disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. In this paper, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

Modal density function and number of propagating modes in ducts

Often raised questions in duct sound propagation studies involve the total number of propagating modes, the number of propagating radial modes for a particular spinning lobe number, and the number of modes possible between two given values of cutoff ratio or eigenvalue. These questions can be answered approximately by using the modal distribution function which is the integral of the modal density function for ducts in a manner similar to that previously published for architectural acoustics. The modal density functions are derived for rectangular and circular ducts with a uniform steady flow. Results from this continuous theory are compared to the actual (discrete) modal distributions.

Rice, E. J.↗