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

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↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Complex Fermi–Löwdin orbital self-interaction correction

Here, this paper introduces the use of complex Fermi orbital descriptors (FODs) in the Fermi–Löwdin self-interaction-corrected density functional theory (FLOSIC). With complex FODs, the Fermi–Löwdin orbitals (FLOs) that are used to evaluate the SIC correction to the total energy become complex. Complex FLO-SIC (cFLOSIC) calculations based on the local spin density approximation produce total energies that are generally lower than the corresponding energies found with FLOSIC restricted to real orbitals (rFLOSIC). The cFLOSIC results are qualitatively similar to earlier Perdew–Zunger SIC (PZ-SIC) calculations using complex orbitals. The energy lowering stems from the exchange–correlation part of the self-interaction correction. The Hartree part of the correction is more negative in rFLOSIC. The energy difference between real and complex solutions is greater for more strongly hybridized FLOs in atoms and for FLOs corresponding to double and triple bonds in molecules. The case of N2 is examined in detail to show the differences between the real and complex FLOs. We show that the complex triple-bond orbitals are simple, and physically appealing combinations of π and σ g orbitals that have not been discussed before. Consideration of complex FODs, and resulting unitary transformations, underscores the fact that FLO centroids are not necessarily good guesses for FOD positions in a FLOSIC calculation.

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

DMC-ICE13 : Ambient and high pressure polymorphs of ice from diffusion Monte Carlo and density functional theory

Ice is one of the most important and interesting molecular crystals, exhibiting a rich and evolving phase diagram. Recent discoveries mean that there are now 20 distinct polymorphs; a structural diversity that arises from a delicate interplay of hydrogen bonding and van der Waals dispersion forces. This wealth of structures provides a stern test of electronic structure theories, with Density Functional Theory (DFT) often not able to accurately characterize the relative energies of the various ice polymorphs. Thanks to recent advances that enable the accurate and efficient treatment of molecular crystals with Diffusion Monte Carlo (DMC), we present here the DMC-ICE13 dataset; a dataset of lattice energies of 13 ice polymorphs. This dataset encompasses the full structural complexity found in the ambient and high-pressure molecular ice polymorphs, and when experimental reference energies are available, our DMC results deliver sub-chemical accuracy. Using this dataset, we then perform an extensive benchmark of a broad range of DFT functionals. Of the functionals considered, revPBE-D3 and RSCAN reproduce reference absolute lattice energies with the smallest error, while optB86b-vdW and SCAN+rVV10 have the best performance on the relative lattice energies. Our results suggest that a single functional achieving reliable performance for all phases is still missing, and that care is needed in the selection of the most appropriate functional for the desired application. The insights obtained here may also be relevant to liquid water and other hydrogen-bonded and dispersion-bonded molecular crystals.

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

Determination of single-crystal elastic moduli of Li RE F 4 ( RE =Y, Gd, and Tb) by resonant ultrasound spectroscopy

The tetragonal fluoro-scheelite Li RE F 4 compounds (RE = rare earth) have been shown to exhibit a variety of useful optical and magnetic properties. While LiYF 4 has been widely studied, many of the fundamental thermodynamic properties of other members of this family remain unknown. Here, we report the complete elastic tensors (C ij ) of single-crystalline LiYF 4 , LiTbF 4 , and LiGdF 4 using resonant ultrasound spectroscopy and density functional theory (DFT). Here we compare the results for LiYF 4 with prior experimental results using time-of-flight ultrasound methods. This is the first report, however, of the experimental elastic tensors of LiTbF 4 and LiGdF 4 . The present results point to a softening of the elastic moduli of the Li RE F 4 system when Y is replaced by the larger ionic radius of Tb or Gd. Furthermore, we find that just 0.3% doping with Nd on the Y site also leads to a slight softening of the moduli. The variation of the elastic moduli as a function of temperature up to 216°C was also measured. A nearly linear softening of all seven independent elastic moduli was observed with increasing temperature. Phonon dispersions and phonon density of states obtained by DFT support the experimental finding of a significantly higher sound velocity due to lighter Y atoms in LiYF 4 , as compared to heavier LiTbF 4 and LiGdF 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Origin of the injection-dependent emission blueshift and linewidth broadening of III-nitride light-emitting diodes

III-nitride light-emitting diodes (LEDs) exhibit an injection-dependent emission blueshift and linewidth broadening that is severely detrimental to their color purity. By using first-principles multi-scale modeling that accurately captures the competition between polarization-charge screening, phase-space filling, and many-body plasma renormalization, we explain the current-dependent spectral characteristics of polar III-nitride LEDs fabricated with state-of-the-art quantum wells. Our analysis uncovers a fundamental connection between carrier dynamics and the injection-dependent spectral characteristics of light-emitting materials. For example, polar III-nitride LEDs offer poor control over their injection-dependent color purity due to their poor hole transport and slow carrier-recombination dynamics, which forces them to operate at or near degenerate carrier densities. Designs that accelerate carrier recombination and transport and reduce the carrier density required to operate LEDs at a given current density lessen their injection-dependent wavelength shift and linewidth broadening.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Epitaxial hexagonal boron nitride with high quantum efficiency

Two-dimensional (2D) hexagonal boron nitride (h-BN) is one of the few materials showing great promise for light emission in the far ultraviolet (UV)-C wavelength, which is more effective and safer in containing the transmission of microbial diseases than traditional UV light. In this report, we observed that h-BN, despite having an indirect energy bandgap, exhibits a remarkably high room-temperature quantum efficiency (~60%), which is orders of magnitude higher than that of other indirect bandgap material, and is enabled by strong excitonic effects and efficient exciton-phonon interactions. This study offers a new approach for the design and development of far UV-C optoelectronic devices as well as quantum photonic devices employing 2D semiconductor active regions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗