Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2021

Kohn-Sham density functional theory of Abelian anyons

In this work, we develop a density functional treatment of noninteracting Abelian anyons, which is capable, in principle, of dealing with a system of a large number of anyons in an external potential. Comparison with exact results for a few particles shows that the model captures the behavior qualitatively and semiquantitatively, especially in the vicinity of the fermionic statistics. We then study anyons with statistics parameter 1+1/n, which are thought to condense into a superconducting state. An indication of the superconducting behavior is the mean-field result that, for uniform density systems, the ground-state energy increases under the application of an external magnetic field independent of its direction. Our density functional theory based analysis does not find that to be the case for finite systems of anyons, which can accommodate a weak external magnetic field through density transfer between the bulk and the boundary rather than through transitions across effective Landau levels, but the “Meissner repulsion” of the external magnetic field is recovered in the thermodynamic limit as the effect of the boundary becomes negligible. We also consider the quantum Hall effect of anyons, and show that its topological properties, such as the charge and statistics of the excitations and the quantized Hall conductance, arise in a self-consistent fashion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2021

Opening band gaps of low-dimensional materials at the meta-GGA level of density functional approximations

The quasiparticle band structure can be properly described by Hedin's GW approximation (GW), at a high computational cost. For band gaps, semilocal density functionals up to the generalized gradient approximation (GGA) level cannot compete with the accuracy of hybrid-based approximations or GW. Meta-GGA density functionals with a strong dependence on the kinetic energy density ingredient can potentially give wider band gaps compared with GGAs. The recent TASK meta-GGA density functional from Aschebrock and Kümmel [ Phys. Rev. Research 1 , 033082 (2019) ], is constructed with an enhanced nonlocality in the generalized Kohn-Sham scheme and therefore harbors great opportunities for band gap prediction. Although this approximation was found to yield excellent band gaps of bulk solids, this accuracy cannot be straightforwardly transferred to low-dimensional materials. Additionally, the reduced screening of these materials results in larger band gaps compared with their bulk counterparts, as an additional barrier to overcome. In this paper we demonstrate how the alteration of this functional affects the band gaps of monolayers and nanoribbons and present accurate band gaps competing with the revised Heyd-Scuseria-Ernzerhof (HSE06) approximation. In order to achieve this goal, we have modified the TASK functional (a) by changing the tight upper bound for one- or two-electron systems ( h X 0 ) from 1.174 to 1.29 and (b) by changing the limit of the interpolation function f X ( α → ∞ ) of the TASK functional that interpolates the exchange enhancement factor F X ( s , α ) from α = 0 to 1. The resulting modified TASK (mTASK) was tested for various materials from three dimensions to two dimensions to one dimension (nanoribbons) and was compared with the results of the higher-level hybrid functional HSE06 or with the G 0 W 0 approximation within many-body perturbation theory. We find that mTASK systematically improves the band gaps and band structures of two-dimensional (2D) and 1D systems, without significantly affecting the accuracy of the original TASK for the bulk 3D materials, when compared with the Perdew-Burke-Ernzerhof (PBE) GGA and the strongly constrained and appropriately normed (SCAN) meta-GGA. We further demonstrate the applicability of mTASK by assessing the band structures of transition metal dichalcogenide nanoribbons with respect to various bending curvatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. Here, we have focused our workflow on TiH 2 as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

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 OSTI2021

Investigating 3,4-bis(3-nitrofurazan-4-yl)furoxan detonation with a rapidly tuned density functional tight binding model

In this work, we describe a machine learning approach to rapidly tune density functional tight binding models for the description of detonation chemistry in organic molecular materials. Resulting models enable simulations on the several 10s of ps scales characteristic to these processes, with “quantum-accuracy.” We use this approach to investigate early shock chemistry in 3,4-bis(3-nitrofurazan-4-yl)furoxan, a hydrogen-free energetic material known to form onion-like nanocarbon particulates following detonation. We find that the ensuing chemistry is significantly characterized by the formation of large C x N y O z species, which are likely precursors to the experimentally observed carbon condensates. Beyond utility as a means of investigating detonation chemistry, the present approach can be used to generate quantum-based reference data for the development of full machine-learned interatomic potentials capable of simulation on even greater time and length scales, i.e., for applications where characteristic time scales exceed the reach of methods including Kohn–Sham density functional theory, which are commonly used for reference data generation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗