Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local density approximations”

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 OSTI2021

Exploring and enhancing the accuracy of interior-scaled Perdew–Zunger self-interaction correction

The Perdew–Zunger self-interaction correction (PZ-SIC) improves the performance of density functional approximations for the properties that involve significant self-interaction error (SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by local scaling self-interaction correction (LSIC) of the PZ-SIC to the local spin density approximation (LSDA). Here, we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of the atomic number Z in the asymptotic expansion of the exchange–correlation (xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator z σ , which distinguishes one-electron regions from many-electron regions. LSIC+ applied to the LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA), and almost close to the strongly constrained and appropriately normed (SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled down SIC to functionals beyond the LSDA requires a gauge transformation of the functional’s energy density. Here, the resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds.

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

Symmetry breaking and self-interaction correction in the chromium atom and dimer

Density functional approximations to the exchange–correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew–Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ~4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. In conclusion, these results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions.

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

Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions

Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu + , Zn + , Zn 2+ , and Mn + in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Furthermore, analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.

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

Study of self-interaction errors in density functional predictions of dipole polarizabilities and ionization energies of water clusters using Perdew-Zunger and locally scaled self-interaction corrected methods

Here, we studied the effects of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modelled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz. the local spin density approximation (LDA), the Perdew-Burke-Ernzherof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-Lowdin orbital SIC (FLO-SIC) framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.58 Bohr 3 with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but over-corrects, resulting in underestimated polarizabilities in SIC-LDA, -PBE, and -SCAN. Finally, we applied a recently proposed local-scaling SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute error of 0.08 Bohr 3 for LSIC-LDA and 0.06 Bohr 3 for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as an absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAE of 0.4 eV for LSIC-LDA and 0.04 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies while the LSIC-PBE ionization energies are more accurate than reported GW results.

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

Vertical detachment energies of ammonia cluster anions using self-interaction-corrected methods

Systems with weakly bound extra electrons impose great challenges to semilocal density functional approximations (DFAs), which suffer from self-interaction errors. Small ammonia clusters are one such example of weakly bound anions where the extra electron is weakly bound. We applied two self-interaction correction (SIC) schemes, viz., the well-known Perdew–Zunger and the recently developed locally scaled SIC (LSIC) with the local spin density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation (GGA), and the SCAN meta-GGA functionals to calculate the vertical detachment energies (VDEs) of small ammonia cluster anions (NH3)n–. Our results show that the LSIC significantly reduces the errors in calculations of VDE with LSDA and PBE-GGA functionals leading to better agreement with the reference values calculated with coupled cluster singles and doubles with perturbative triples [CCSD(T)]. Accurate prediction of VDE as an absolute of the highest occupied molecular orbital (HOMO) is challenging for DFAs. Our results show that VDEs estimated from the negative of HOMO eigenvalues with the LSIC-LSDA and Perdew–Zunger SIC-PBE are within 11 meV of the reference CCSD(T) results. Here, the LSIC method performs consistently well for the VDE estimates, from both the total energy differences and the absolute HOMO eigenvalues.

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

Atomic ionization: sd energy imbalance and Perdew–Zunger self-interaction correction energy penalty in 3d atoms

To accurately describe the energetics of transition metal systems, density functional approximations (DFAs) must provide a balanced description of s- and d- electrons. One measure of this is the sd transfer error, which has previously been defined as E ( 3 d n − 1 4 s 1 ) − E ( 3 d n − 2 4 s 2 ) . Theoretical concerns have been raised about this definition due to its evaluation of excited-state energies using ground-state DFAs. A more serious concern appears to be strong correlation in the 4s 2 configuration. Here, we define a ground-state measure of the sd energy imbalance, based on the errors of s- and d-electron second ionization energies of the 3d atoms, that effectively circumvents the aforementioned problems. We find an improved performance as we move from the local spin density approximation (LSDA) to the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) to the regularized and restored Strongly Constrained and Appropriately Normed (r 2 SCAN) meta-GGA for first-row transition metal atoms. However, we find large (∼2 eV) ground-state sd energy imbalances when applying a Perdew–Zunger 1981 self-interaction correction. This is attributed to an “energy penalty” associated with the noded 3d orbitals. A local scaling of the self-interaction correction to LSDA results in a balance of s- and d-errors.

Science & Technology - Other Topics↗
DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2022

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C 2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B 2 , C 2 , O 2 , and F 2 in the local spin density approximation (LSDA), the Perdew–Burke–Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C 2 . Here, this work shows that spin symmetry breaking in singlet C 2 , which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

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

Density-related properties from self-interaction corrected density functional theory calculations

The Perdew–Zunger self-interaction correction (PZ-SIC) removes unphysical electron self-interaction from calculations employing standard density functional approximations. Doing so improves many computed properties, bringing them into better agreement with experimental observations or with results from high-level quantum chemistry calculations. However, while PZ-SIC generally corrects in the right direction relative to the corresponding reference values, in many cases, it over-corrects. For this reason, scaled-down versions of PZ-SIC have been proposed and investigated. These approaches have mostly employed exterior scaling in which SIC correction terms are scaled in the same way at every point in space. Recently, a new local, or interior, scaling SIC method was proposed on non-empirical grounds to restore a property of the exact, but unknown, density functional that is broken in PZ-SIC. Here, the scaling at each point depends on the character of the charge density at that point. However, the local scaling can be done in various ways while still restoring the behavior of the exact functional. In this work, we compare and contrast the performance of various interior scaling approaches for addressing over-corrections of calculated molecular dipole moments and atomic polarizabilities and properties that reflect the nature of the electronic charge density.

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

Spin response and neutrino mean free path in neutron matter

The longitudinal and transverse density- and spin-density response functions in pure neutron matter (PNM) are derived over a wide range of densities within the time dependent local spin density approximation. The underlying density functional was derived from an Auxiliary Field Diffusion Monte Carlo computation of the equation of state of unpolarized and fully spin polarized PNM. In order to assess the dependence of the results on the specific underlying nucleon–nucleon Hamiltonian, we used both the phenomenological Argonne AV8' + UIX force, and local chiral forces up to next-to-next-to-leading-order (N 2 LO). Finally, the resulting response function has then been applied to the study of the neutrino mean free path in PNM.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Approximate bounds and temperature dependence of adiabatic connection integrands for the uniform electron gas

Thermal density functional theory is commonly used in simulations of warm dense matter, a highly energetic phase characterized by substantial thermal effects and by correlated electrons demanding quantum mechanical treatment. Methods that account for temperature dependence, such as Mermin–Kohn–Sham finite-temperature density functional theory and free energy density functional theory, are now employed with more regularity and available in many standard code packages. However, approximations from zero-temperature density functional theory are still often used in temperature-dependent simulations using thermally weighted electronic densities as an input to exchange–correlation functional approximations, a practice known to miss temperature-dependent effects in the exchange–correlation free energy of these systems. Here, in this work, the temperature-dependent adiabatic connection is demonstrated and analyzed using a well-known parameterization of the uniform electron gas free energy. Useful tools based on this formalism for analyzing and constraining approximations of the exchange–correlation at zero temperature are leveraged for the finite-temperature case. Inspired by the Lieb–Oxford inequality, which provides a lower bound for the ground-state exchange–correlation energy, bounds for the exchange–correlation at finite temperatures are approximated for various degrees of electronic correlation.

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

The exact exchange–correlation potential in time-dependent density functional theory: Choreographing electrons with steps and peaks

The time-dependent exchange–correlation potential has the unusual task of directing fictitious non-interacting electrons to move with exactly the same probability density as true interacting electrons. This has intriguing implications for its structure, especially in the non-perturbative regime, leading to step and peak features that cannot be captured by bootstrapping any ground-state functional approximation. Here, we review what has been learned about these features in the exact exchange–correlation potential of time-dependent density functional theory in the past decade or so and implications for the performance of simulations when electrons are driven far from any ground state.

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

Revealing quasi-excitations in the low-density homogeneous electron gas with model exchange–correlation kernels

Time-dependent density functional theory within the linear response regime provides a solid mathematical framework to capture excitations. The accuracy of the theory, however, largely depends on the approximations for the exchange–correlation (xc) kernels. Away from the long-wavelength (or q = 0 short wave-vector) and zero-frequency (ω = 0) limit, the correlation contribution to the kernel becomes more relevant and dominant over exchange. The dielectric function, in principle, can encompass xc effects relevant to describe low-density physics. Furthermore, besides collective plasmon excitations, the dielectric function can reveal collective electron–hole excitations, often dubbed “ghost excitons.” Besides collective excitons, the physics of the low-density regime is rich, as exemplified by a static charge-density wave that was recently found for r s > 69, and was shown to be associated with softening of the plasmon mode. These excitations are seen to be present in much higher density 2D homogeneous electron gases of r s ≳ 4. Here, in this work, we perform a thorough analysis with xc model kernels for excitations of various nature. The uniform electron gas, as a useful model of real metallic systems, is used as a platform for our analysis. We highlight the relevance of exact constraints as we display and explain screening and excitations in the low-density region.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Self-interaction correction in water–ion clusters

We study the importance of self-interaction errors in density functional approximations for various water–ion clusters. We have employed the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in conjunction with LSDA, PBE, and SCAN to describe binding energies of hydrogen-bonded water–ion clusters, i.e., water–hydronium, water– hydroxide, water–halide, as well as non-hydrogen-bonded water–alkali clusters. In the hydrogen-bonded water–ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules, because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water– ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher–level calculations. Here, the non-hydrogen-bonded water–alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies (≤ 50 meV/H 2 O) for hydrogen-bonded water–ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement.

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

Electronic structure, magnetic properties, and exchange splitting of gadolinium intermetallics

Here, we report here the electronic structure calculations of the magnetic moment, magneto-crystalline anisotropy energy (MAE), and density of states (DOS) of Gd-metal, GdAl 2 , GdAl 3 , GdCo 2 , GdCo 3 , and GdCo 5 . From the DOS, the exchange splitting energy is derived to map ferrimagnetic to paramagnetic phase transition (Curie temperature, T C ). The employed theoretical method here is based on advanced density functional theory, specifically the linearized augmented plane wave method (LAPW) within the local spin density approximation (LSDA), including spin–orbit coupling (SOC) and Hubbard model parameter (U). DOS analysis of both Gd-Al and Gd-Co systems reveals a strong correlation between exchange splitting energy and T C . The Gd-Co system exhibits an enormous exchange splitting energy among Co-3d electrons, vastly overwhelming the vanishing Gd 5d exchange splitting. Finally, these findings suggest that exchange splitting of Co-3d electrons plays a key role in the ground state magnetism and magnetic transition temperatures in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Efficient Langevin dynamics for “noisy” forces

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Additionally, stochastic approaches provide a gentler system-size dependency at the cost of introducing “noisy” forces, which could limit the efficiency of the sampling. When the forces are deterministic, the first-order Langevin dynamics (FOLD) offers efficient sampling by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator [G. Mazzola and S. Sorella, Phys. Rev. Lett. 118, 015703 (2017)]. However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure that severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗