Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Oxygen vacancy induced site-selective Mott transition in LaNiO 3

While defects such as oxygen vacancies in correlated materials can modify their electronic properties dramatically, understanding the microscopic origin of electronic correlations in materials with defects has been elusive. Lanthanum nickelate with oxygen vacancies, LaNiO 3 - x , exhibits the metal-to-insulator transition as the oxygen vacancy level x increases from the stoichiometric LaNiO 3 . In particular, LaNiO 2.5 exhibits a paramagnetic insulating phase, also stabilizing an antiferromagnetic state below T N ≃ 152 K . Here, we study the electronic structure and energetics of LaNiO 3 - x using first principles. We find that LaNiO 2.5 exhibits a “site-selective” paramagnetic Mott insulating state at T ≃ 290 K as obtained using density functional theory plus dynamical mean field theory ( DFT + DMFT ). The Ni octahedron site develops a Mott insulating state with strong correlations as the Ni e g orbital is half-filled while the Ni square-planar site with apical oxygen vacancies becomes a band insulator. Overall, our oxygen vacancy results cannot be explained by the pure change of the Ni oxidation state alone within the rigid band-shift approximation. Our DFT + DMFT density of states explains that the peak splitting of unoccupied states in LaNiO 3 - x measured by the experimental x-ray absorption spectra originates from two nonequivalent Ni ions in the vacancy-ordered structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Topological surface states on the nonpolar (110) and (111) surfaces of SmB 6

In order to clarify the controversial issue of the topological nature in a mixed-valent Kondo system, SmB 6 , we have explored the surface states on the nonpolar (110) surface of SmB 6 , employing both angle-resolved photoemission spectroscopy (ARPES) experiment and ab initio density-functional theory (DFT) band calculations. Based on ARPES spectroscopic fingerprints and the DFT surface band structures, we ascribe the observed spectral weights at $\bar{X}$ and $\bar{Y}$ on the (110) surface Brillouin zone to topological surface states (TSSs) of “topological insulator (TI)” nature and of “topological crystalline insulator (TCI)” nature, respectively. With varying the chemical potential, the double Dirac cones of the TCI nature exhibit a Lifshitz transition of Fermi surfaces with intriguing spin textures. Further, we have also examined the TSSs on the nearly nonpolar (111) surface of SmB 6 in connection with a recently reported ARPES result and proposed a way to probe the Dirac points that are buried in the bulk-projected bands.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Elimination of the linearization error in APW/LAPW basis set: Dirac-Kohn-Sham equations

Here, a detailed account of the implementation of equations of the relativistic density functional theory (RDFT) using basis sets of APW/LAPW type with flexible extensions provided by local orbitals is given. Earlier discoveries of the importance of the high derivative local orbital (HDLO) extension of the APW/LAPW basis set for enhancing the accuracy of DFT calculations are confirmed using a fully relativistic approach and α – U as an example. High energy local orbitals (HELO's), however indispensable for GW calculations, are considerably less efficient in enhancing the accuracy of DFT applications. It is shown that a simplified approach to the relativistic effects, namely considering them only inside the muffin-tin (MT) spheres, produces basically identical results (as compared to fully relativistic approach) for the electronic free energy of the five materials considered in this work. By comparing the effect of the simplified approach on the electronic free energy with its effect on the electronic kinetic energy we conclude that the insensitivity of the free energy to the way we describe the relativistic effects in the interstitial region is related to the variational property of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comprehensive anisotropic linear optical properties of the Weyl semimetals TaAs and NbAs

We report that TaAs and NbAs are two of the earliest identified Weyl semimetals that possess many intriguing optical properties, such as chirality-dependent optical excitations and giant second harmonic generation (SHG). Linear and nonlinear optics have been employed as tools to probe the Weyl physics in these crystals. Here, we extend these studies to address two important aspects: determining the complete anisotropic dielectric response and exploring if and how they can reveal essential Weyl physics. We determine the complete anisotropic dielectric functions of TaAs and NbAs by combining spectroscopic ellipsometry and density functional theory (DFT). Parameterized Lorentz oscillators are reported from 1.2 - 6 eV (experiment) and 0 - 6 eV (DFT), and good agreement is shown between them. Both linear and nonlinear optical properties have been previously reported to reveal Weyl physics. We suggest that strong optical resonances from trivial bands are the likely origin of the large optical SHG previously reported at these energies. Furthermore, by comparing the contribution of a small k -space centered around the Weyl cones with the total linear dielectric function, we find that these contributions are highly anisotropic and are < 25% of the total dielectric function below 0.5 eV; abov1 eV, these contributions are negligible. Thus, the study of Weyl physics using optical techniques requires very low energies, and even there, a careful assessment is required in distinguishing the much smaller contributions of the Weyl bands from the dominant contributions of the trivial bands and Drude response to the total dielectric function.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First-principles theory for Schottky barrier physics

We develop a first-principles theory for Schottky barrier physics. The Poisson equation is solved selfconsistently with the electrostatic charge density over the entire barrier using the density functional theory (DFT) electronic structure converged locally, allowing computation of a Schottky barrier entirely from DFT involving thousands of atomic layers in the semiconductor (SC). Here, the induced charge in the bulk consists of conduction and valence band charges from doping and band bending, as well as charge from the evanescent states in the gap of the SC. The Schottky barrier height (SBH) is determined when the induced charge density and the induced electrostatic potential reach self-consistency. Tests on the GaAs-graphene and Si/Al heterostructures yield SBH, width, along with depletion and inversion layers obtained self-consistently as functions of temperature and bulk doping.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantum physical reality of polar-nonpolar oxide heterostructures

Conducting interfaces between polar and nonpolar insulating oxides, e.g., LaAlO 3 /SrTiO 3 , have generated interest for both fundamental physics and oxide-electronics applications. Current understanding is based on an amalgamation of a classical electrostatic model (polar catastrophe model) that was originally derived for semi-infinite solids and quantum density-functional-theory (DFT) results on ultrathin films. Here we report comprehensive DFT calculations that unveil a very different purely quantum physical reality. We show that, for ultrathin polar films, the interfacial dipole does not control the electrostatic potential in the polar film—the surface and interface play equal roles, and the absence or presence of centrosymmetry in the physical LAO film results in different, purely quantum mechanisms for the generation of a conducting interface, neither involving physical-charged transfer. Predictions are made that can be tested and can guide technology development.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of the electronic charge gap on LO bond-stretching phonons in undoped La 2 CuO 4 calculated using LDA + U

Typical density-functional theory calculations that wrongly predict undoped cuprates to be metallic also predict Cu-O half- and full-breathing phonon energies that are significantly softer than observed, presumably because of weak on-site Coulomb repulsion on the Cu 3d orbitals. We used DFT+U calculations with antiferromagnetic supercells of La 2 CuO 4 to establish correlation between the on-site repulsion strength, tuned via adjusting the value of U, and phonon dispersions. Here, we find that breathing and half-breathing phonons reach experimental values when U is tuned to obtain the correct optical gap and magnetic moments. We demonstrate that using distorted supercells within DFT+U is a promising framework to model phonons in undoped cuprates and other perovskite oxides with complex, interrelated structural and electronic degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Candidate structure for the H 2 -PRE phase of solid hydrogen

Experimental progress finally reached the metallic solid hydrogen phase, which was predicted by Wigner and Huntington over 80 years ago. However, the different structures in the phase diagram are still being debated due to the difficulty of diffraction experiments for high-pressured hydrogen. The determination of crystal structures under extreme condition is both of the basic condensed matter physics, and in planetary science: the behavior of giant gaseous planets (e.g., Jupiter and Saturn) strongly depends on the properties of inner high-pressured hydrogen. This work describes possible structures appearing under high pressures of 400~600 GPa. We applied a structural search using particle swarm optimization with density functional theory (DFT) to propose several candidate structures. For these structures, we performed fixed-node diffusion Monte Carlo simulations combined with DFT zero-point energy corrections to confirm their relative stability. We found P2 1 /c-8 as a promising candidate structure for the H 2 -PRE phase. P2 1 /c-8 is predicted the most stable at 400 and 500 GPa. P2 1 /c-8 reproduces qualitatively the IR spectrum peaks observed in the H 2 -PRE phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Spin-polarized imaging of strongly interacting fermions in the ferrimagnetic state of the Weyl candidate CeBi

CeBi has an intricate magnetic phase diagram whose fully polarized state has recently been suggested as a Weyl semimetal, though the role of f states in promoting strong interactions has remained elusive. Here we focus on the less-studied but also time-reversal symmetry-breaking ferrimagnetic phase of CeBi, where our density functional theory (DFT) calculations predict additional Weyl nodes near the Fermi level E F . We use spin-polarized scanning tunneling microscopy and spectroscopy to image the surface ferrimagnetic order on the itinerant Bi p states, indicating their orbital hybridization with localized Ce f states. We observe suppression of this spin-polarized signature at E F , coincident with a Fano line shape in the conductance spectra, suggesting the Bi p states partially Kondo screen the f magnetic moments, and this p – f hybridization causes strong Fermi-level band renormalization. Furthermore, the p-band flattening is supported by our quasiparticle interference measurements, which also show band splitting in agreement with DFT, painting a consistent picture of a strongly interacting magnetic Weyl semimetal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetization reversal driven by electron localization-delocalization crossover in the inverse spinel Co 2 VO 4

Neutron diffraction, magnetization, and muon spin relaxation measurements, supplemented by density functional theory (DFT) calculations are employed to unravel temperature-driven magnetization reversal in inverse spinel Co 2 VO 4 . All measurements show a second-order magnetic phase transition at T C = 168 K to a collinear ferrimagnetic phase. Neutron diffraction measurements reveal two antiparallel ferromagnetic (FM) sublattices, belonging to magnetic ions on two distinct crystal lattice sites, where the relative balance between the two sublattices determine the net FM moment in the unit cell. As the evolution of the ordered moment with temperature differs between the two sublattices, the net magnetic moment reaches a maximum at T NC = 138 K and reverses its sign at T MR = 65 K. The DFT results suggest that the underlying microscopic mechanism for the reversal is a delocalization of the unfilled 3d -shell electrons on one sublattice just below T C , followed by a gradual localization as the temperature is lowered. Furthermore, this delocalized-localized crossover is supported by muon spectroscopy results, as strong T 1 relaxation observed below T C indicates fluctuating internal fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Correlation energy of the uniform electron gas determined by ground-state conditional probability density functional theory

Conditional-probability density functional theory (CP-DFT) is a formally exact method for finding correlation energies from Kohn-Sham DFT without evaluating an explicit energy functional. We present details on how to generate accurate exchange-correlation energies for the ground-state uniform gas. We also use the exchange hole in a CP antiparallel spin calculation to extract the high-density limit. We give a highly accurate analytic solution to the Thomas-Fermi model for this problem, showing its performance relative to Kohn-Sham and may be useful at high temperatures. We explore several approximations to the CP potential. Furthermore, results are compared to accurate parameterizations for both exchange-correlation energies and holes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic structure of intrinsic defects in c -gallium nitride: Density functional theory study without the jellium approximation

Here, we report the first nonjellium, systematic, density functional theory (DFT) study of intrinsic and extrinsic defects and defect levels in zinc-blende (cubic) gallium nitride. We use the local moment counter charge (LMCC) method, the standard Perdew-Becke-Ernzerhoff (PBE) exchange-correlation potential, and two pseudopotentials, where the Ga 3$\textit{d}$ orbitals are either in the core ($d^0$) or explicitly in the valence set ($d^{10}$). We studied 64, 216, 512, and 1000 atom supercells, and demonstrated convergence to the infinite limit, crucial for delineating deep from shallow states near band edges, and for demonstrating the elimination of finite cell-size errors. Contrary to common claims, we find that exact exchange is not required to obtain defect levels across the experimental band gap. As was true in silicon, silicon carbide, and gallium arsenide, the extremal LMCC defect levels of the aggregate of defects yield an effective LMCC defect band gap that is within 10% of the experimental gap (3.3 eV) for both pseudopotentials. We demonstrate that the gallium vacancy is more complicated than previously reported. There is dramatic metastability–a nearest-neighbor nitrogen atom shifts into the gallium site, forming an antisite, nitrogen vacancy pair, which is more stable than the simple vacancy for positive charge states. Our assessment of the $d^0$ and $d^{10}$ pseudopotentials yields minimal differences in defect structures and defect levels. The better agreement of the $d^0$ lattice constant with experiment suggests that the more computationally economical $d^0$ pseudopotentials are sufficient to achieve the fidelity possible within the physical accuracy of DFT, and thereby enable calculations in larger supercells necessary to demonstrate convergence with respect to finite size supercell errors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A-type antiferromagnetic order in semiconducting EuMg 2 Sb 2 single crystals

Eu-based Zintl-phase materials EuA 2 Pn 2 (A = Mg, In, Cd, Zn; Pn = Bi, Sb, As, P) have generated significant recent interest owing to the complex interplay of magnetism and band topology. Here, we investigated the crystallographic, magnetic, and electronic properties of the layered Zintlphase single crystals of EuMg 2 Sb 2 with the trigonal CaAl 2 Si 2 crystal structure (space group $P\bar{3}m1$). Electrical resistivity measurements complemented with angle-resolved photoemission spectroscopy (ARPES) studies and density functional theory (DFT) calculations find an activated behavior with intrinsic conductivity at high temperatures indicating a semiconducting electronic ground state with a narrow energy gap of 370 meV. Magnetic susceptibility and zero-field heat capacity measurements indicate that the compound undergoes antiferromagnetic (AFM) ordering at the Néel temperature T N = 8.0(2) K. Here, zero-field neutron-diffraction measurements reveal that the AFM ordering is A-type where the Eu spins (Eu 2+ , S = 7/2) arranged in ab-plane layers are aligned ferromagnetically in the ab plane and the Eu spins in adjacent layers are aligned antiferromagnetically. Eu-moment reorientation within the ab planes in the trigonal AFM domains associated with a very weak inplane magnetic anisotropy is also evident below T N at low fields < 0.05 T. Although isostructural semimetallic EuMg 2 Bi 2 is reported to host Dirac surface states, the observation of narrow-gap semiconducting behavior in EuMg 2 Sb 2 implies a strong role of spin-orbit coupling (SOC) in tuning the electronic states of these materials. Our DFT studies also suggest, besides the SOC, the more electronegative and smaller Sb than Bi shifts the low-lying conduction bands along the Γ-A direction to higher energy, resulting in an indirect bulk band gap between the Γ and M points for EuMg 2 Sb 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Low-energy interband transition in the infrared response of the correlated metal SrVO 3 in the ultraclean limit

We studied the low-energy electronic response of the prototypical correlated metal SrVO 3 in the ultraclean and disordered limit using infrared spectroscopy and density functional theory plus dynamical mean field theory calculations (DFT+DMFT). A strong optical excitation at 70 meV is observed in the optical response of the ultraclean samples but is hidden by the low-energy Drude-like response from intraband excitations in the more disordered samples. DFT+DMFT calculations reveal that this optical excitation originates from interband transitions between the bands split by orbital off-diagonal hopping, which has often been ignored in cubic systems, such as SrVO 3 . A memory function analysis of the optical data shows that this interband transition can lead to deviations of optical self-energy from the expected Fermi-liquid behavior. Our findings demonstrate that analysis schemes employed to extract many-body effects from optical spectra may be oversimplified to study the true electronic ground state and that improvements in material quality can guide efforts to refine theoretical approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Low-energy moiré phonons in twisted bilayer van der Waals heterostructures

We develop a low-energy continuum model for phonons in twisted moiré bilayers, based on a configuration-space approach. In this approach, interatomic force constants are obtained from density functional theory (DFT) calculations of untwisted bilayers with various in-plane shifts. This allows for efficient computation of phonon properties for any small twist angle, while maintaining DFT-level accuracy. Based on this framework, we show how the low-energy phonon modes, including interlayer shearing and layer-breathing modes, vary with the twist angle. As the twist angle decreases, the frequencies of the low-energy modes are reordered and the atomic displacement fields corresponding to phonon eigenmodes break translational symmetry, developing periodicity on the moiré length scale. Here, we demonstrate the capabilities of our model by calculating the phonon properties of three specific structures: bilayer graphene, bilayer molybdenum disulfide (MoS 2 ), and molybdenum diselenide-tungsten diselenide (MoSe 2 -WSe 2 ).

2-dimensional systems↗