Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid density functional calculations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

First-principles calculations of defects and electron–phonon interactions: Seminal contributions of Audrius Alkauskas to the understanding of recombination processes

First-principles calculations of defects and electron–phonon interactions play a critical role in the design and optimization of materials for electronic and optoelectronic devices. The late Audrius Alkauskas made seminal contributions to developing rigorous first-principles methodologies for the computation of defects and electron–phonon interactions, especially in the context of understanding the fundamental mechanisms of carrier recombination in semiconductors. Alkauskas was also a pioneer in the field of quantum defects, helping to build a first-principles understanding of the prototype nitrogen-vacancy center in diamond, as well as identifying novel defects. Here, we describe the important contributions made by Alkauskas and his collaborators and outline fruitful research directions that Alkauskas would have been keen to pursue. Audrius Alkauskas’ scientific achievements and insights highlighted in this article will inspire and guide future developments and advances in the field.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Iodine Vacancies do not Cause Nonradiative Recombination in Halide Perovskites

The iodine vacancy (V I ) has frequently been discussed as a strong nonradiative recombination center in halide perovskites. This proposition was mainly based on the presence of charge-state transition levels in the band gap, as found in early first-principles calculations. In this work, we perform accurate hybrid-density-functional calculations for V I in CsPbI 3 , CsSnI 3 , and CsGeI 3 and find that V I does not have any transition levels in the band gap in CsPbI 3 , in contrast to the results from calculations based on semilocal functionals. The iodine vacancy V I does introduce levels in the band gap in CsSnI 3 and CsGeI 3 , but our explicitly computed nonradiative capture coefficients demonstrate that V I has a negligible impact on nonradiative recombination. Our study corrects a misunderstanding of the role of V I in the iodide-based perovskites, and shifts the focus toward identifying and mitigating actual recombination centers in order to further improve the optoelectronic performance.

14 SOLAR ENERGY↗
DOE OSTI2022

Role of carbon and hydrogen in limiting n -type doping of monoclinic ( Al x Ga 1 – x ) 2 O 3

In this study, we use hybrid density functional calculations to assess n -type doping in monoclinic ( Al x Ga 1 – x ) 2 O 3 alloys. We focus on silicon, the most promising donor dopant, and study the structural properties, formation energies, and charge-state transition levels of its various configurations. We also explore the impact of carbon and hydrogen, which are common impurities in metal-organic chemical vapor deposition (MOCVD). In Ga 2 O 3 , Si Ga is an effective shallow donor, but in Al 2 O 3 Si Al acts as a DX center with a ( + / – ) transition level in the band gap. Interstitial hydrogen acts as a shallow donor in Ga 2 O 3 but behaves as a compensating acceptor in n -type Al 2 O 3 . Interpolation indicates that Si is an effective donor in ( Al x Ga 1 – x ) 2 O 3 up to 70% Al, but it can be compensated by hydrogen already at 1% Al. We also assess the diffusivity of hydrogen and study complex formation. Si cation – H complexes have relatively low binding energies. Substitutional carbon on a cation site acts as a shallow donor in Ga 2 O 3 , but can be stable in a negative charge state in ( Al x Ga 1 – x ) 2 O 3 when x > 5 % . Substitutional carbon on an oxygen site ( C O ) always acts as an acceptor in n -type ( Al x Ga 1 – x ) 2 O 3 , but will incorporate only under relatively oxygen-poor conditions. C O – H complexes can actually incorporate more easily, explaining observations of carbon-related compensation in Ga 2 O 3 grown by MOCVD. We also investigate C cation – H complexes, finding they have high binding energies and act as compensating acceptors when x > 56 % ; otherwise the hydrogen just passivates the unintentional carbon donors. C-H complex formation explains why MOCVD-grown Ga 2 O 3 can exhibit record-low free-carrier concentrations, in spite of the unavoidable incorporation of carbon. Our study highlights that, while Si is in principle a suitable shallow donor in ( Al x Ga 1 – x ) 2 O 3 alloys up to high Al compositions, control of unintentional impurities is essential to avoid compensation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hybrid Functional Calculations for Antimony Doping in CdTe

CdTe-based solar cells are leading thin-film photovoltaic technology, with efficiencies over 22%, but still much lower than the theoretical maximum of 29%. Further improvements will rely on increasing the open-circuit voltage Voc, which, in turn, depends on carrier density and lifetime. Using hybrid density functional calculations, we investigate Sb doping of CdTe, focusing on its limitation as shallow acceptor and the formation of compensating AX center. Paying special attention to supercell size and effects of spin-orbit coupling, we predict an ionization energy of 116 meV in the dilute limit, much closer to recent experimental value of 103 meV from temperature-dependent Hall measurements in bulk single crystals, and in contrast to much larger values from previous calculations. We also find that the Sb-related AX centers are not major compensation centers in Sb-doped CdTe.

14 SOLAR ENERGY↗
DOE OSTI2024

The deep-acceptor nature of the chalcogen vacancies in 2D transition-metal dichalcogenides

Abstract Chalcogen vacancies in the semiconducting monolayer transition-metal dichalcogenides (TMDs) have frequently been invoked to explain a wide range of phenomena, including both unintentional p-type and n-type conductivity, as well as sub-band gap defect levels measured via tunneling or optical spectroscopy. These conflicting interpretations of the deep versus shallow nature of the chalcogen vacancies are due in part to shortcomings in prior first-principles calculations of defects in the semiconducting two-dimensional TMDs that have been used to explain experimental observations. Here we report results of hybrid density functional calculations for the chalcogen vacancy in a series of monolayer TMDs, correctly referencing the thermodynamic charge transition levels to the fundamental band gap (as opposed to the optical band gap). We find that the chalcogen vacancies are deep acceptors and cannot lead to n-type or p-type conductivity. Both the (0/−1) and (−1/−2) transition levels occur in the gap, leading to paramagnetic charge states S = 1 / 2 and S = 1, respectively, in a collinear-spin representation. We discuss trends in terms of the band alignments between the TMDs, which can serve as a guide to future experimental studies of vacancy behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of defects in ultra-high gain in fast planar tin gallium oxide UV-C photodetector by MBE

Here we report ultra-high responsivity of epitaxial (Sn x Ga 1–x ) 2 O 3 (TGO) Schottky UV-C photodetectors and experimentally identified the source of gain as deep-level defects, supported by first principles calculations. Epitaxial TGO films were grown by plasma-assisted molecular beam epitaxy on (–201) oriented n-type β-Ga 2 O 3 substrates. Fabricated vertical Schottky devices exhibited peak responsivities as high as 3.5 ×10 4 A/W at –5 V applied bias under 250 nm illumination with sharp cutoff shorter than 280 nm and fast rise/fall time in milliseconds order. Hyperspectral imaging cathodoluminescence (CL) spectra were examined to find the mid-bandgap defects, the source of this high gain. Irrespective of different tin mole fractions, the TGO epilayer exhibited extra CL peaks at the green band (~2.20 eV) not seen in β-Ga 2 O 3 along with enhancement of the blue emission-band (~2.64 eV) and suppression of the UV emission-band. Based on hybrid functional calculations of the optical emission expected for defects involving Sn in β-Ga 2 O 3 , V Ga –Sn complexes are proposed as potential defect origins of the observed green and blue emission-bands. Such complexes behave as acceptors that can efficiently trap photogenerated holes and are predicted to be predominantly responsible for the ultra-high photoconductive gain in the Sn-alloyed Ga 2 O 3 devices by means of thermionic emission and electron tunneling. Regenerating the V Ga –Sn defect complexes by optimizing the growth techniques, we have demonstrated a planar Schottky UV-C photodetector of the highest peak responsivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Defect Stability in CdTe Based on Formation Energies and Migration Barriers

Native point defects are thought to play a key role in CdTe, either as compensation centers in intentionally doped material, as a source of conductivity in nominally undoped material, or as electron–hole recombination centers. Here, the discussion of their concentration and impact has often centered only on formation energies and transition levels. Using hybrid density functional calculations, including the effects of spin–orbit coupling (SOC), we discuss the stability of native point defects in CdTe based on their formation energies and migration barriers. We show that although Cd interstitials are the lowest energy donor defects, they are unstable at room temperature due to a low migration barrier. They are important for maintaining charge neutrality during growth or annealing at high temperatures, but once the material is brought to room temperature, they are not frozen in as often assumed and are expected to anneal out, leaving the other more stable defects to determine the conductivity. Taking this into account in the solution of the charge neutrality equation, we are able to predict the conductivity type and carrier concentrations that are in good agreement with experimental observations.

14 SOLAR ENERGY↗
DOE OSTI2022

Role of morphology in defect formation and photo-induced carrier instabilities in amorphous indium oxide

Ab initio molecular dynamics liquid-quench simulations and hybrid density functional calculations are performed to model the effects of room-temperature atomic fluctuations and photo-illumination on the structural and electronic properties of amorphous sub-stoichiometric In 2 O 2.96 . A large configurational ensemble is employed to reliably predict the distribution of localized defects as well as their response to the thermal and light activation. The results reveal that the illumination effects on the carrier concentration are greater in amorphous configurations with shorter In–O bond length and reduced polyhedral sharing as compared to the structures with a more uniform morphology. Further, the obtained correlation between the photo-induced carrier density and the reduction in the number of fully coordinated In-atoms implies that metal oxides with a significant fraction of crystalline/amorphous interfaces would show a more pronounced response to illumination. Photo-excitation also produces In–O 2 –In defects that have not been previously found in sub-stoichiometric amorphous oxides; these defects are responsible for carrier instabilities due to overdoping.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

The role of native defects and their diffusivity in CdTe

CdTe is a leading thin-film technology for solar-cell applications. Basic features such as the optimum band gap of 1.5 eV and simple manufacturing process have been driving the research on improving the current record efficiency of 22% for CdTe solar cells toward the theoretical limit of 32%. Low p-type doping and short carrier lifetimes in the CdTe absorber have been suggested as the main limiting factors and attributed to defect compensation and carrier recombination. Understanding the role of point defects and impurities in CdTe is crucial to solving these problems. Despite many years of research, the sources of compensation and the microscopic recombination mechanisms are still unclear. By using hybrid density functional calculations with spin-orbit coupling and large supercells, we investigate the electronic properties, the formation energies, and the diffusivities of the native defects in CdTe. As possible sources of compensation for p-type conductivity, we discuss the structural and electronic properties of Cd interstitial, Te vacancies, and TeCd antisites. The stability of each defect is discussed in terms of both the defect formation energy and the migration barrier. We also address possible mechanisms that lead to non-radiative recombination.

14 SOLAR ENERGY↗
DOE OSTI2023

Oxides related to cadmium telluride solar cells

Polycrystalline cadmium telluride (CdTe) is a leading material in photovoltaic technology due to its high absorption coefficient and near-optimum bandgap of 1.44 eV. It is known that CdTe film processing can promote surface oxidation depending on the growth environment and upon exposure to different oxidation conditions. For example, CdTeO3 forms when the CdTe film is treated in heated dry air, while in humid air, CdTe2O5 is detected. These oxides feature tellurium in the oxidation state +4 compared to the +2 state in CdTe. Other possible relevant oxides are CdO, TeO2, and TeO3. Using hybrid density functional calculations, we studied the electronic structure of these oxide materials and their band alignment to CdTe, which are essential parameters in the characterization of the interfaces at grain boundaries. The goal is to understand their stability and possible effects on passivating grain boundaries. The results are compared to the available experimental data.

14 SOLAR ENERGY↗
DOE OSTI2024

Self-compensation of group-V acceptors in CdTe

Cadmium Telluride is at the core of an important thin-film technology for photovoltaics that is already commercially available, yet the CdTe-based solar cell efficiency remains limited at 22%, well below the theoretical limit of ~30%. Increasing the hole concentration is crucial for higher efficiency, and group-V elements such as As, P, and Sb are potential dopants as they are shallow acceptors. Nevertheless, group-V doped p-type CdTe often exhibits low doping activation, and the compensation source remains debated. Here, we performed hybrid density functional calculations with spin-orbit coupling to investigate possible sources of hole compensation in group-V doped CdTe. First, regarding possible self-compensation of the group-V dopants, we find that the formation of AX centers is unlikely since they are found to be unstable relative to the shallow acceptor forms. However, if the group-V dopants come in during growth (such as dimer molecules As2, P2, and Sb2), we find that the impurity atoms, which would occupy nearest neighbor sites, maintain the V-V bonds, limiting the hole density. For the native defects, our study reveals that Cd interstitial is the lowest energy donor defect in p-type CdTe. Still, it has a small migration barrier of 0.5 eV, making it highly mobile and unstable at room temperature. The Te vacancy is the next lowest formation energy donor. The migration barrier of 1.4 eV indicates that the Te vacancy is stable at room temperature. The antisite CdTe is also a donor, with low formation energy and stable at room temperature, potentially limiting the hole concentration. Our results, therefore, shed light on possible compensation centers and some guidance on how to avoid them.

14 SOLAR ENERGY↗
DOE OSTI2022

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this research, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

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

Photoelectron spectra of early 3 d -transition metal dioxide molecular anions from GW calculations

Photoelectron spectra of early 3 d -transition metal dioxide anions, Sc O 2 - , Ti O 2 - , V O 2 - , Cr O 2 - , and Mn O 2 - , are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point ( G 0 W 0 @PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV–0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew–Burke–Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G 0 W 0 @PBE0, the G -only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3 d -transition metal dioxide anions. Adding eigenvalue self-consistency also in W , on the other hand, worsens the agreement with experiment overall. Overall, our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

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

Optimized Auxiliary Functions for Robust Mitigation of Finite-Size Errors in Periodic Hybrid Density Functional Theory

When calculating properties of periodic systems at the thermodynamic limit (TDL), the dominant source of finite size error (FSE) arises from the long-range Coulomb interaction, and can manifest as a slowly converging quadrature error when approximating an integral in the reciprocal space by a finite sum. The singularity subtraction (SS) method offers a systematic approach for reducing this quadrature error and thus the FSE. Here, in this work, we first investigate the performance of the SS method in the simplest setting, aiming at reducing the FSE in exact exchange calculations by subtracting the Coulomb contribution with a single, adjustable Gaussian auxiliary function. We demonstrate that a simple fitting method can robustly estimate the optimal Gaussian width and leads to rapid convergence toward the TDL. Furthermore, we suggest new forms of the auxiliary function, whose optimal parameters could also be determined through least-squares fitting. For a range of semiconductors and insulators, the proposed auxiliary functions achieve robust, millihartree-level accuracy in hybrid density functional theory calculations, including cases with sparse k-meshes and large basis sets.

Quiton, Stephen Jon [University of California, Ber↗
DOE OSTI2021

Analysis of defects in In 2 O 3 :H synthesized in presence of water vapor and hydrogen gas mixture

Using hybrid functional-based density functional theory calculations, we analyze the structure and kinetics of defects formed in two competing synthesis routes to prepare hydrogen-doped In 2 O 3 films, using a hydrogen and oxygen gas mixture vs synthesis in the presence of water vapor. For both of these synthesis routes, we find that H + is the dominant defect species: when the Fermi level is close to the conduction band, H + has a lower formation energy than other intrinsic or extrinsic defects. Our results also suggest that water molecules spontaneously split into H + (which occupies octahedral voids) and OH – interstitials (which occupies vacant oxygen lattice sites or oxygen vacancies). From the analysis of the binding energies between these different defects, we conclude that these defects do not cluster and are most likely to stay spatially distributed throughout the films. In addition, the sum of formation energies of an oxygen (i.e., O$^{2–}_i$) and a H + interstitial is close to the formation energy of a OH – interstitial, meaning that water molecules are completely split into 2H + and O$^{2–}_i$ at the synthesis conditions. Further, in the presence of H 2 + O 2 gas mixture, oxygen interstitials occupy oxygen vacancies while hydrogen interstitials occupy vacant oxygen lattice sites and form bonds with lattice oxygens. Our analysis of the defect equilibria suggests that the hydrogen content in films synthesized in the presence of water vapor is higher than films synthesized in the presence of a hydrogen gas mixture. At high dopant concentrations, a hydrogen bond network is formed in the system and this leads to large distortions in the lattice.

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

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

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

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure—An experimental and theoretical study

The 57Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α→ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth’s inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe2O3. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

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

Carbon complexes in highly C-doped GaN

We investigate the properties of heavily C-doped GaN grown by hydride vapor phase epitaxy using both optical experiments and hybrid density functional theory calculations. Previous work has established that carbon acceptors (C N ) give rise to a yellow luminescence band near 2.2 eV along with a blue luminescence band near 2.9 eV. Photoluminescence measurements show the yellow band shifting as a function of carbon concentration, suggesting a change in the behavior of carbon species as carbon content increases. With hybrid density functional theory we calculate the electrical and optical behavior of carbon centers containing multiple carbon impurities, which may arise in heavily doped material. We compare the behavior of these complexes to the isolated centers, and find that the dicarbon donor-acceptor (C Ga –C N ) complex is a candidate to explain the shift in the yellow luminescence peak. Tricarbon complexes have high formation energies and modest binding energies, and also give rise to optical transitions that are inconsistent with the observed spectra. We also identify the split dicarbon interstitial on the gallium site as a low-energy species with a large binding energy that may act to compensate carbon acceptors. Here, local vibrational modes are calculated for carbon impurity centers, and we compare these results to recent experiments. Dicarbon and tricarbon complexes involving C Ga and C N exhibit modes that are only slightly higher than the isolated species, while carbon interstitials and related complexes give rise to vibrational modes much higher than C Ga and C N .

36 MATERIALS SCIENCE↗