Search NASA⌕ Search

SEARCH · Search NASA

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fermi surface topology and magnetotransport properties of superconducting Pd 3 Bi 2 Se 2

Pd 3 Bi 2 Se 2 is a rare realization of a superconducting metal with a non-zero topological invariant. Here, in this study, we report the growth of high-quality single crystals of layered Pd 3 Bi 2 Se 2 with a superconducting transition at T c ≈ 0.80 K and upper critical fields of ~10 mT and ~5 mT for the in plane and out-of-plane directions, respectively. Our density functional theory (DFT) calculations reveal three pairs of doubly degenerate bands crossing the Fermi level all displaying clear three dimensional dispersion consistent with the overall low electronic anisotropy (<2). The multiband electronic nature of Pd 3 Bi 2 Se 2 is evident in magneto-transport measurements, yielding a sign changing Hall resistivity at low temperatures. The magnetoresistance is non-saturating and follows Kohler’s scaling rule. We interpret the magneto-transport data in terms of open orbits that are revealed in the DFT calculated Fermi surface. de Haas-van Alphen (dHvA) oscillation measurements using torque magnetometry on single crystals yield four frequencies for out-of-plane fields: F α = (150 ± 26) T, F β = (293 ± 10) T, F γ = (375 ± 20) T and F η = (1017 ± 12) T, with the low frequency dominating the spectrum. Through the measurement of angular dependent dHvA oscillations and DFT calculations we identify the F α frequency with an approximately ellipsoidal electron pocket centered on the L 2 point of the Brillouin zone. Lifshitz-Kosevich analysis of the dHvA oscillations reveals a small cyclotron effective mass m* = (0.11 ± 0.02)m 0 and a nontrivial Berry phase for the dominant orbit. The presence of nontrivial topology in a bulk superconductor positions Pd 3 Bi 2 Se 2 as a potential candidate for exploring topological superconductivity.

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

Role of magnetic and structural symmetry breaking in forming the Mott insulating gap in Nb 3 ⁢Cl 8

The α-phase of the gapped insulator Nb 3 Cl 8 has recently emerged as the long-sought critical testing bed for examining the importance of strong interelectronic correlation vs symmetry breaking in understanding insulation of such Mott compounds. Structural symmetry breaking detected by density functional theory (DFT) energy lowering (such as dimer formation, disproportionation, or Jahn-Teller distortions) explains insulation in both d-electron Mott-like systems and in non-d-electron cases without recourse to strong correlation. Yet, in Nb 3 Cl 8 , structural symmetry breaking alone (viz. formation of Nb trimers) fails to explain insulation, leading instead to a partially occupied metallic flat band, in contrast with experimental observations. We examine the role of magnetic symmetry breaking, noting that Nb 3 Cl 8 is an observed paramagnet (not an antiferromagnet), thus potentially carrying also short-range ordered magnetic moments. Describing the latter as a polymorphous distribution of nonzero local moments with total zero net magnetization is demonstrated to lower the DFT total energy, while gapping the system without recourse to strong correlation or long-range magnetic order. This suggests that degeneracy removal by symmetry breaking in mean-field-like approaches—either structural, or magnetic, or both—can reduce or eliminate the need for strong correlation, allowing the use of DFT for such Mott systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Light-Induced Renormalization of the Dirac Quasiparticles in the Nodal-Line Semimetal ZrSiSe

In nodal-line semimetals, linearly dispersing states form Dirac loops in the reciprocal space with a high degree of electron-hole symmetry and a reduced density of states near the Fermi level. The result is reduced electronic screening and enhanced correlations between Dirac quasiparticles. In this study we investigate the electronic structure of ZrSiSe, by combining time- and angle-resolved photoelectron spectroscopy with ab initio density functional theory (DFT) complemented by an extended Hubbard model ( DFT + U + V ) and by time-dependent DFT + U + V . We show that electronic correlations are reduced on an ultrashort timescale by optical excitation of high-energy electrons-hole pairs, which transiently screen the Coulomb interaction. Our findings demonstrate an all-optical method for engineering the band structure of a quantum material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Towards predictive band gaps for halide perovskites: Lessons from one-shot and eigenvalue self-consistent GW

Halide perovskites constitute a chemically diverse class of crystals with great promise as photovoltaic absorber materials, featuring band gaps between about 1 and 3.5 eV depending on composition. Their diversity calls for a general computational approach to predicting their band gaps. However, such an approach is still lacking. Here, we use density functional theory (DFT) and ab initio many-body perturbation theory within the GW approximation to compute the quasiparticle or fundamental band gap of a set of ten representative halide perovskites: CH 3 NH 3 PbI 3 ( MAPbI 3 ), MAPbBr 3 , CsSnBr 3 , ( MA ) 2 BiTlBr 6 , Cs 2 TlAgBr 6 , Cs 2 TlAgCl 6 , Cs 2 BiAgBr 6 , Cs 2 InAgCl 6 , Cs 2 SnBr 6 , and Cs 2 Au 2 I 6 . Comparing with recent measurements, we find that a standard generalized gradient exchange-correlation functional can significantly underestimate the experimental band gaps of these perovskites, particularly in cases with strong spin-orbit coupling (SOC) and highly dispersive band edges, to a degree that varies with composition. We show that these nonsystematic errors are inherited by one-shot G 0 W 0 and eigenvalue self-consistent GW 0 calculations, demonstrating that semilocal DFT starting points are insufficient for MAPbI 3 , MAPbBr 3 , CsSnBr 3 , ( MA ) 2 BiTlBr 6 , Cs 2 TlAgBr 6 , and Cs 2 TlAgCl 6 . On the other hand, we find that DFT with hybrid functionals leads to an improved starting point and G W 0 results in better agreement with experiment for these perovskites. Our results suggest that G W 0 with hybrid functional-based starting points are promising for predicting band gaps of systems with large SOC and dispersive bands in this technologically important class of semiconducting crystals.</

36 MATERIALS SCIENCE↗
DOE OSTI2021

Prediction of low-Z collinear and noncollinear antiferromagnetic compounds having momentum-dependent spin splitting even without spin-orbit coupling

The over 60 years old Rashba-Dresselhaus effect predicts spin-orbit coupling (SOC) induced momentum-dependent spin splitting and spin polarization in materials with noncentrosymmetric structures. Strong SOC induced effects usually require high-atomic number (Z) elements such as rare-earth elements. It has recently been pointed out that antiferromagnets could hold SOC-independent spin splitting and spin polarization. In the present work we develop the spatial and magnetic symmetry conditions enabling such antiferromagnet (AFM)-induced spin splitting, dividing the 1651 magnetic space groups into seven different spin splitting prototypes (SST-1 to SST-7). This analysis places the physics of AFM spin splitting (SST-4) within the broader context of symmetry conditions that enable the more familiar forms of spin splitting, such as ferromagnetic Zeeman effect (SST-5), nonmagnetic no spin splitting (SST-6), and the nonmagnetic Rashba and Dresselhaus effects (SST-7). The AFM-induced spin splitting and spin polarization do not necessarily require breaking of inversion symmetry or the presence of SOC, hence can exist even in centrosymmetric, low-Z light element compounds, considerably broadening the material base for spin manipulations. We use the “inverse design” approach of first formulating the target property (here, spin splitting in low-Z compounds not restricted to low symmetry structures), then derive the enabling physical design principles—the magnetic symmetry conditions—to search realizable compounds that satisfy these a priori design principles. This process uncovers 422 magnetic space groups (160 centrosymmetric and 262 noncentrosymmetric) that could hold AFM-induced, SOC-independent spin splitting and spin polarization. We then search for stable compounds following such enabling symmetries. We investigate the electronic and spin structures of some selected prototype compounds by density functional theory (DFT) and find spin textures that are different than the traditional Rashba-Dresselhaus patterns and exist even in the absence of SOC effect. We provide the DFT results for all antiferromagnetic spin splitting prototypes (SST-1, SST-2, SST-3, SST-4), and concentrate on two limits of SST-4 that are particularly unusual: When spin splitting is momentum dependent (just like the Rashba effect) but is enabled in antiferromagnets even in the absence of SOC in the Hamiltonian. This includes examples of (a) centrosymmetric SST-4A compounds (e.g., orthorhombic LaMnO 3 illustrating collinear AFM, as well as cubic NiS 2 illustrating noncollinear AFM) and (b) noncentrosymmetric SST-4B compounds (e.g., rhombohedral MnTiO 3 illustrating collinear AFM and hexagonal ScMnO 3 illustrating noncollinear AFM). The symmetry design principles outlined here, along with their transformation into an inverse design material search approach and DFT verification, could open the way to their experimental examination.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-driven approach to parameterize SCAN + U for an accurate description of 3 d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the DFT + U method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. Here we develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-U correction, and (3) DFT/DFT + U compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for 3d transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve SCAN/SCAN + U compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize DFT + U performance for other target physical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal properties of the metallic delafossite PdCo⁢O 2 : A combined experimental and first-principles study

Metallic delafossite materials (e.g., PdCo⁢O 2 , PtCo⁢O 2 ), have attracted much recent attention due to record-high oxide conductivities, the origins of which remain unclear. Relatively little attention has been paid to their related thermal properties, however. Here, we address this via wide temperature range experimental studies of the crystal structure, thermal expansion, and specific heat of single-crystal PdCo⁢O 2 , combined with density-functional theory (DFT) calculations of the electronic and phononic densities of states, and thus thermal properties. PdCo⁢O 2 is shown to retain the $R\bar{3}m$ space group from 12 to 1000 K, exhibiting a- and c-axis thermal expansion in good quantitative agreement with DFT-based lattice dynamics calculations. The Co-O bond lengths additionally elucidate the stability of the low-spin state of the nominally Co 3+ ions, which is a notable difference between the edge-shared Co-O octahedra in PdCo⁢O 2 and the corner-shared octahedra in Co-based perovskites. Measurements of specific heat from 1.9 to 400 K provide accurate values for the Debye temperature and Sommerfeld coefficient, the phononic part being describable via a combined Debye-Einstein approach (accounting for high-frequency oxygen-related optical phonons), with excess intermediate-temperature specific heat due to a prominent low-energy peak in the phonon density of states. Most significantly, all electronic and phononic contributions to the specific heat are shown to be remarkably closely reproduced by DFT-based calculations, establishing quantitative understanding of key thermal properties of the metallic delafossite PdCo⁢O 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure and electronic tunability of acene alkylamine based layered hybrid organic-inorganic perovskites from first principles

Twelve layered hybrid organic-inorganic perovskites combining oligoacene derivatives [phenylmethylammonium (PMA), naphthylmethylammonium (NMA), anthrylmethylammonium (AMA), and tetrylmethylammonium (TMA)] and lead halides (Cl - , Br - , and I - anions) are investigated by first-principles density functional theory (DFT), showing broad, rational tunability of band gap, quantum well type, and spin-dependent energy band properties. Six compounds are known from previous syntheses and are used to devise a computational search space for likely low-energy structures. Among the six known compounds, a refined structure is identified for (NMA) 2 ⁢ PbCl 4 and a new, lower-energy structure is suggested for (AMA) 2 ⁢PbCl 4 . The DFT based search methodology is next applied to predict the likely structures for the six unknown compounds. Computationally predicted energy levels for all 12 compounds from spin-orbit coupled hybrid DFT reveal tunable type I quantum well alignments, with frontier orbitals located on the inorganic component, on the organic component, or on both, as controlled by cation and anion selection. In conclusion, several structures spontaneously break local and global inversion symmetry, causing strong spin splitting of the conduction bands [up to 0.12 eV for (PMA) 2 ⁢PbCl 4 ] and showing potential for spin-dependent transport properties and future spintronic and tunable chiroptical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hybridized kagome bands and induced Ti magnetism in RTi3Bi4 (R = Nd, Sm, Gd) kagome metals studied using angle-resolved photoemission spectroscopy and x-ray magnetic circular dichroism

Kagome materials are known for hosting emergent quantum phenomena driven by the interaction between different lattice, charge and spin orders. Here, we present a detailed angle resolved photoemission (ARPES), density functional theory (DFT) and x-ray magnetic circular dichroism (XMCD) study of the electronic and magnetic structure of RTi3Bi4 (R = Nd, Sm, Gd). ARPES and DFT demonstrate that the bulk electronic band structure is dominated by the hybridization of the Ti bands, and the weak electron-like pocket at Γ is identified as a surface state. The isotropic XAS profile of the M 4,5-edge of the rare earth is consistent with the presence of R 3+ oxidation state. Using the XMCD sum rules, backed by the atomic multiplet theory calculations, we obtain the spin and orbital magnetic moments. The Ti L2,3-edge XMCD reveals the presence of a small magnetic moment in GdTi3Bi4, presumably driven by the proximity of the Ti kagome layers to the zigzag chains of Gd, while the total magnetic moment of Gd is shared by the f and d electrons. Our combined XMCD, ARPES and DFT study brings an important piece of information to understand the spin flip transitions and anomalous Hall effect observed in the RTi3Bi4 kagome metals.

Lim, C.↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2026

Insights on carbon dioxide adsorption in a flexible 1-D coordination polymer from in situ X-ray scattering and density functional theory

A combination of in situ small-angle X-ray scattering (SAXS) microstructure characterization over scales ranging from 1 nm to 10 µm and powder X-ray diffraction (XRD) structure characterization under various gas/pressure/temperature conditions with density functional theory (DFT) calculations provides new insights for the CO2 sorption behavior of a one-dimensional porous coordination polymer: catena-bis­(di­benzoyl­methanato)(4,4′-bi­pyridyl)­nickel(II), denoted NiDBM-Bpy. The NiDBM-Bpy chains are held together by van der Waals forces, but the structure of guest-free NiDBM-Bpy is unsolved due to a lack of suitable crystals and high-quality powder XRD patterns. Nevertheless, SAXS and powder XRD can follow microstructural and structural changes as a function of gas pressure, composition and temperature. Both mixed-gas flow and static supercritical CO2 regimes are explored experimentally. DFT calculations are used to model the structural variation associated with XRD changes and hysteresis in the sorption isotherms. XRD and DFT calculations suggest that an orthorhombic Fddd structure with two CO2 per Ni emerges following a transition from the structure with lower molar volume and symmetry that exists without CO2 present.

Allen, Andrew J.↗
DOE OSTI2024

Rotational Millimeter-Wave Shoe Scanner Using the Discrete Fourier Transform for Backprojection-Based Image Reconstruction

An active 3D microwave / millimeter-wave shoe scanner was previously developed at the Pacific Northwest National Laboratory (PNNL) using two linear arrays scanned over a rectilinear aperture. The radar system chirps a frequency sweep from 10-40 GHz. These frequencies allow imaging through optically opaque material such as leather, rubber, plastics, and other dielectrics. The system was designed to detect concealed items in the soles of shoes while allowing people to leave their shoes on through a security checkpoint. To shrink the footprint of the system, a new iteration of the design has been developed that scans the two linear arrays over a circular aperture. This new footprint opens the possibility of it being installed in the floor of a cylindrical millimeter-wave body scanner. The backprojection-based multilayer dielectric image reconstruction developed at PNNL can easily handle arbitrary spatial sampling, accommodating the new rotational shoe scanner design. Commonly, the fast Fourier transform (FFT) is used to efficiently compute the range response from the data collected by the system as a preprocessing step to the backprojection algorithm. It was found that converting to range using the discrete Fourier transform (DFT) directly has some advantages over the FFT. For example, nonlinear and non-uniform frequency sweeps can easily be compensated for during the computation of the DFT and only the range bins of interest need to be computed and their spacing can be chosen arbitrarily. Because the range conversion step of the image reconstruction is the fastest part of the process there is very little speed penalty for using the DFT over the FFT and it can even increase the speed of image reconstruction when the ranges of interest are fewer than the total span that is calculated in the FFT.

Millimeter-wave imaging, microwave imaging, shoe s↗
DOE OSTI2023

Voltage-Dependent First-Principles Simulation of Insertion of Chloride Ions into Al/Al 2 O 3 Interfaces Using the Quantum Continuum Approximation

Experiments have shown that pitting corrosion can develop in aluminum surfaces at potentials > − 0.5 V relative to the standard hydrogen electrode (SHE). Until recently, the onset of pitting corrosion in aluminum has not been rigorously explored at an atomistic scale because of the difficulty of incorporating a voltage into density functional theory (DFT) calculations. We introduce the Quantum Continuum Approximation (QCA) which self-consistently couples explicit DFT calculations of the metal-insulator and insulator-solution interfaces to continuum Poisson-Boltzmann electrostatic distributions describing the bulk of the insulating region. By decreasing the number of atoms necessary to explicitly simulate with DFT by an order of magnitude, QCA makes the first-principles prediction of the voltage of realistic electrochemical interfaces feasible. After developing this technique, we apply QCA to predict the formation energy of chloride atoms inserting into oxygen vacancies in Al(111)/α-Al 2 O 3 (0001) interfaces as a function of applied voltage. We predict that chloride insertion is only favorable in systems with a grain boundary in the Al 2 O 3 for voltages > − 0.2 V (SHE). Here our results roughly agree with the experimentally demonstrated onset of corrosion, demonstrating QCA's utility in modeling realistic electrochemical systems at reasonable computational cost.

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

Electronic Structure Engine

INQ is an engine for electronic structure calculations. INQ is implemented as a C++ library of components designed to build a parallel, portable code for electronic structure calculations, including ground state and excited states. INQ is based on density functional theory (DFT). It can calculate ground state properties in DFT, and also excited states using time-dependent DFT (TDDFT) in real-time and linear-response. It implements different approximations for the exchange and correlation part: semi-local functionals like LDAs, GGAs and metaGGAs, and hybrid functionals that are implemented through the ACE approach for fast execution.

ANDRADE-VALENCIA, XAVIERI↗