Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Assessing the roles of Cu- and Ag-deficient layers in chalcopyrite-based solar cells through first principles calculations

Chalcopyrites are a demonstrated material platform for realizing efficient thin-film photovoltaics, with the most well known Cu(In,Ga)Se 2 (CIGS)-based solar cells exceeding 23%. Several factors, including flexibility in tuning the absorber bandgap, enhanced surface treatments, and the electrically benign nature of common defects are responsible for the existing high performance and future promise in chalcopyrite-based photovoltaic devices. The introduction of Cu-poor phases (also known as ordered-vacancy compounds or OVCs) between the absorber and buffer layers in CIGS solar cells is known to enhance device performance; however, the overall properties and role of OVCs remain poorly understood. Using first principles calculations based on the density functional theory with screened hybrid functionals, here we explore the electronic structure and stability of OVCs and their band offsets with defect-free chalcopyrite layers in Cu- and Ag-based compounds (ABX 2 where A=Cu, Ag; B=In, Ga, Al; and X=S, Se). Using AB 3 X 5 and AB 5 X 8 stoichiometries as model OVC systems, we report on the variation of the bandgap with the A/B ratio and discuss the trends in other Cu- and Ag-based chalcopyrites beyond CuInSe 2 . We find that the valence and conduction bands are lower in energy in OVCs with respect to the parent ABX 2 chalcopyrite owing to a reduced p–d interaction between X and A atoms. We additionally perform device-level simulations to assess the implications of the results, finding that the valence band offsets of OVCs are favorable, while the conduction band offsets of chalcopyrites beyond CuInSe 2 -based absorbers may be detrimental in conventional solar cell device designs.

14 SOLAR ENERGY↗

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↗

Self-trapped hole and impurity-related broad luminescence in $β$-Ga 2 O 3

Here we explore the luminescence properties of self-trapped holes and impurity-related acceptors using one-dimensional configuration coordinate diagrams derived from hybrid functional calculations. The photoluminescence spectrum of as-grown $β$-Ga 2 O 3 typically consists of a broad band in the wavelength region from ultraviolet to green and is often dominated by an impurity independent ultraviolet band that is commonly attributed to self-trapped holes. Here, we use the self-trapped hole as a benchmark to evaluate the accuracy of the theoretical defect luminescence spectra and estimate the optical properties of Mg Ga , Be Ga , Ca Ga , Cd Ga , Zn Ga , Li Ga , and N O acceptor impurities, as well as their complexes with hydrogen donors. We also explore $V$ Ga acceptors complexed with hydrogen and Si Ga donor impurities. The results show that these defects can give rise to broad luminescence bands peaking in the infrared to visible part of the spectrum, making them potential candidates for the defect origin of broad luminescence bands in $β$-Ga 2 O 3 .

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

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↗

Prospects for n -type doping of (Al x Ga 1– x ) 2 O 3 alloys

Here we systematically explore the properties of group-IV (C, Si, Ge, and Sn) and transition metal (Hf, Zr, and Ta) dopants substituting on the cation site in (Al x Ga 1– x ) 2 O 3 (AlGO) alloys using first-principles calculations with a hybrid functional. In Ga 2 O 3 , each of these dopants acts as a shallow donor. In Al 2 O 3 , they are deep defects characterized by the formation of either DX centers or positive- U (+/0) levels. Combining our calculations of dopant charge-state transition levels with information of the AlGO alloy band structure, we estimate the critical Al composition at which each dopant transitions from being a shallow to a deep donor. We identify Si to be the most efficient dopant to achieve n -type conductivity in high Al-content AlGO alloys, acting as a shallow donor over the entire predicted stability range for AlGO solid solution alloys.

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

The role of water vapor during the synthesis of hydrogen doped In 2 O 3

Hydrogen-doped In 2 O 3 synthesized using magnetron sputtering or atomic layer deposition in the presence of water vapor has high transparency in the near infra-red region and mobility values that are almost three to four times of those in commercially used In 2 O 3 :Sn at a much lower carrier concentration. However, simple questions like what happens when water molecules enter In 2 O 3 or why the dominant charge carrier changes from an oxygen vacancy in In 2 O 3 :Sn to H + in In 2 O 3 are not clear. Using hybrid functional based density functional theory calculations, in this study we show that water molecules spontaneously split into H + and OH – . The H + interstitial bonds with an oxygen in In 2 O 3 , and the OH – interstitial occupies a vacant lattice site or is trapped by an oxygen vacancy. Thus, in agreement with experimental findings, our results suggest that the concentration of oxygen vacancies decreases as the partial pressure of water increases. The charge transition level, ϵ ( – / + ) , of a hydrogen interstitial and the Fermi level at which the formation energies of OH – and H + are equal lie above the conduction band minimum, suggesting that the H + interstitial has a lower formation energy. In addition, water molecules do not exhibit any propensity to form clusters, meaning that H + ions can diffuse and remain spatially separated from each other.

36 MATERIALS SCIENCE↗

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↗

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↗