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 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 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 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 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 OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations

Here we extend our recently developed sparse-stochastic fragmented exchange formalism for ground-state near-gap hybrid DFT to calculate absorption spectra within linear-response time-dependent generalized Kohn-Sham DFT (LR-GKS-TDDFT) for systems consisting of thousands of valence electrons within a grid-based/plane-wave representation. A mixed deterministic/fragmented-stochastic compression of the exchange kernel, here using long-range explicit exchange functionals, provides an efficient method for accurate optical spectra. Both real-time propagation as well as frequency-resolved Casida-equation-type approaches for spectra are presented, and the method is applied to large molecular dyes.

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

Electronic properties of corundum-like Ir 2 O 3 and Ir 2 O 3 -Ga 2 O 3 alloys

In the hexagonal, corundum-like structure, α-Ga 2 O 3 has a bandgap of ~5.1 eV, which, combined with its relatively small electron effective mass, high Baliga's figure of merit, and high breakdown field, makes it a promising candidate for power electronics. Ga 2 O 3 is easy to dope n-type, but impossible to dope p-type, impeding the realization of some electronic device designs. Developing a lattice-matched p-type material that forms a high-quality heterojunction with n-type Ga 2 O 3 would open new opportunities in electronics and perhaps optoelectronic devices. In this work, we studied Ir 2 O 3 as a candidate for that purpose. Using hybrid density functional theory calculations we predict the electronic band structure of α-Ir 2 O 3 and compare that to α-Ga 2 O 3 , and study the stability and electronic properties of α-(Ir x Ga 1–x ) 2 O 3 alloys. We discuss the band offset between the two materials and compare it with recently available experimental data. We find that the Ir d bands that compose the top of the valence band in α-Ir 2 O 3 are much higher in energy than O p bands in α-Ga 2 O 3 , possibly enabling effective p-type doping. Finally, our results provide an insight into using the Ir 2 O 3 or Ir 2 O 3 -Ga 2 O 3 alloys as p-type material lattice-matched to α-Ga 2 O 3 for the realization of p–n heterojunctions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hybrid Density Functional Theory Comparison of Oxygen Release and Solvent Decomposition Kinetics on Li x NiO 2 Surfaces

High-nickel-content layered oxides are among the most promising electric vehicle battery cathode materials. However, their interfacial reactivity with electrolytes and tendency toward oxygen release (possibly yielding reactive 1 O 2 ) remain degradation concerns. Elucidating the most relevant (i.e., fastest) interfacial degradation mechanism will facilitate future mitigation strategies. We apply screened hybrid density functional (HSE06) calculations to compare the reaction kinetics of Li x NiO 2 surfaces with ethylene carbonate (EC) with those of O 2 release. Further, on both the (001) and (104) facets, EC oxidative decomposition exhibits lower activation energies than O 2 release. Our calculations, coupled with previously computed liquid-phase reaction rates of 1 O 2 with EC, strongly question the role of “reactive 1 O 2 ” species in electrolyte oxidative degradation. The possible role of other oxygen species is discussed. To deal with the challenges of modeling Li x NiO 2 surface reactivity, we emphasize a “local structure” approach instead of pursuing the global energy minimum.

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

Factors Affecting the Electron Conductivity in Single Crystal Li 7 La 3 Zr 2 O 12 and Li 7 P 3 S 11

One of the serious challenges in all solid-state Li ion batteries is neutral Li intrusion into the solid-state electrolyte that can ultimately cause catastrophic failure. One possibility for this is due to n-type electron conductivity that induces the reaction Li + + e – → Li 0 at sites where the potential is less than the Li + /Li potential. This paper reports hybrid density functional theory calculations of the electronic conductivity in two prototype single crystalline solid-state electrolytes, cubic Li 7 La 3 Zr 2 O 12 (c-LLZO) and Li 7 P 3 S 11 (LPS). The formation energies of important point defects that can affect electron conductivity are determined, and we find that the mechanism of n-type electron conductivity for both solid-state electrolytes is via “small” electron polaron hopping, where the quotes signify that substantial Li ion rearrangement is associated with the polaron formation and its migration. In both electrolytes, the formation energies for the small polarons at the Fermi energy are too high to generate measurable electron conductivity at room temperature. For c-LLZO, the concentration of electron polarons necessary to ensure charge neutrality from positively charged oxygen vacancies formed in synthesis can be significantly higher. Hence, the electron conductivity could be significant when measured with ion-blocking metal electrodes, and we discuss how the synthesis conditions could affect this magnitude. However, in the solid-state battery, these polarons are replaced by negatively charged Li vacancies so that the electron conductivity should remain minimal. For LPS single crystals, the inherent minimal electron conductivity is independent of synthesis conditions. Here we also show that the cost of forming Li 0 in bulk c-LLZO is enormous due to strain effects so that it could only potentially form at voids, grain boundaries, or around vacancy defects which relax the lattice strain.

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

Insights into the Reactivity of Brookite TiO 2 Nanorods in Liquid Water from Ab Initio Molecular Dynamics Simulations

Brookite TiO 2 , a rare natural polymorph of TiO 2 , has been reported to be an excellent photocatalyst for the production of hydrogen from water and aqueous alcohol solutions, especially when it is reduced and synthesized in the form of nanorods. Here, we investigate the reactivity of stoichiometric and reduced brookite nanorods in liquid water using ab initio molecular dynamics and hybrid density functional theory calculations. Our simulations show a much higher water dissociation fraction on reduced nanorods than on stoichiometric ones, with an accumulation of the resulting bridging hydroxyls (O br H) and terminal hydroxyls (Ti –OH) on different facets of the nanorod. ObrH groups accumulate preferentially on low-energy (210) facets, where they are stabilized by adjacent reduced Ti (Ti 3+ ) sites, while Ti –OH groups prefer to form at the four-fold coordinated Ti atoms on high-energy (010) facets. This hydroxylation pattern also favors the spatial localization of excited holes on the (010) facets. This coupling between water-induced surface chemistry and charge separation underpins the enhanced photocatalytic activity of brookite nanorods, providing useful information for the design of more efficient TiO2-based nanostructures for solar-driven hydrogen evolution.

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

Transition-Metal-Related Quantum Emitters in Wurtzite AlN and GaN

Transition-metal centers exhibit a paramagnetic ground state in wide-bandgap semiconductors and are promising for nanophotonics and quantum information processing. Specifically, there is a growing interest in discovering prominent paramagnetic spin defects that can be manipulated using optical methods. Here, we investigate the electronic structure and magneto-optical properties of Cr and Mn substitutional centers in wurtzite AlN and GaN. We use state-of-the-art hybrid density functional theory calculations to determine level structure, stability, optical signatures, and magnetic properties of these centers. The excitation energies are calculated using the constrained occupation approach and rigorously verified with the complete active space configuration interaction approach. Our simulations of the photoluminescence spectra indicate that $Cr$$^{1+}_{Al}$ in AlN and $Cr$$^{1+}_{Ga}$ in GaN are responsible for the observed narrow quantum emission near 1.2 eV. We compute the zero-field splitting (ZFS) parameters and outline an optical spin polarization protocol for $Cr$$^{1+}_{Al}$ and $Cr$$^{1+}_{Ga}$. Our results demonstrate that these centers are promising candidates for spin qubits.

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

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗