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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2021

Range-separated hybrid functionals for mixed dimensional heterojunctions: Application to phthalocyanines/MoS 2

We analyze the electronic structure and level alignment of transition-metal phthalocyanine (MPc) molecules adsorbed on two-dimensional MoS 2 employing density functional theory (DFT) calculations. We develop a procedure for multi-objective optimal tuning of parameters of range-separated hybrid functionals in these mixed-dimensional systems. Using this procedure, which leads to the asymptotically-correct exchange-correlation potential between molecule and two-dimensional material, we obtain electronic structures consistent with experimental photoemission results for both energy level alignment and electronic bandgaps, representing a significant advance compared to standard DFT methods. We elucidate the MoS 2 valence resonance with the transition-metal phthalocyanine non-frontier 3d orbitals and its dependence on the transition metal atomic number. Based on our calculations, we derive parameter-free, model self-energy corrections that quantitatively accounts for the effects of the heterogeneous dielectric environment on the electronic structure of these mixed-dimensional heterojunctions.

2D materials↗
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 OSTI2020

Radiative capture rates at deep defects from electronic structure calculations

We present a methodology to calculate radiative carrier capture coefficients at deep defects in semiconductors and insulators from first principles. Electronic structure and lattice relaxations are accurately described with hybrid density functional theory. Calculations of capture coefficients provide an additional validation of the accuracy of these functionals in dealing with localized defect states. We also discuss the validity of the Condon approximation, showing that even in the event of large lattice relaxations the approximation is accurate. We test the method on GaAs:VGa-TeAs and GaN:CN, for which reliable experiments are available, and demonstrate very good agreement with measured capture coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulated field-modulated x-ray absorption in titania

Here, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin–orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO 2 ) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t 2g peaks in the Ti L-edge was observed to red shift and the e g peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.

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

Bandgap analysis and carrier localization in cation-disordered ZnGeN 2

The bandgap of ZnGeN 2 changes with the degree of cation site disorder and is sought in light emitting diodes for emission at green to amber wavelengths. By combining the perspectives of carrier localization and defect states, we analyze the impact of different degrees of disorder on electronic properties in ZnGeN 2 , addressing a gap in current studies, which largely focus on dilute or fully disordered systems. The present study demonstrates changes in the density of states and localization of carriers in ZnGeN 2 calculated using bandgap-corrected density functional theory and hybrid calculations on partially disordered supercells generated using the Monte Carlo method. We use localization and density of states to discuss the ill-defined nature of a bandgap in a disordered material and identify site disorder and its impact on the structure as a mechanism controlling electronic properties and potential device performance. Decreasing the order parameter results in a large reduction of the bandgap. The reduction in bandgap is due, in part, to isolated, localized states that form above the valence band continuum associated with nitrogen coordinated by more zinc than germanium. The prevalence of defect states in all but the perfectly ordered structure creates challenges for incorporating disordered ZnGeN 2 into optical devices, but the localization associated with these defects provides insight into the mechanisms of electron/hole recombination in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hybrid functional calculations of electronic structure and carrier densities in rare-earth monopnictides

The structural parameters and electronic structure of rare-earth pnictides are calculated using density functional theory (DFT) with the Heyd, Scuseria, and Ernzerhof (HSE06) screened hybrid functional. Specifically, we focus on RE-V compounds, with RE=La, Gd, Er, and Lu, and V=As, Sb, and Bi, and analyze the effects of spin-orbit coupling and treating the RE4$\textit{f}$ electrons as valence electrons in the projector augmented wave approach. The results of HSE06 calculations are compared with DFT within the generalized gradient approximation (DFT-GGA) and other previous calculations. We find that all these RE-V compounds are semimetals with electron pockets at the $\textit{X}$ point and hole pockets at Γ. Whereas in DFT-GGA the carrier density is significantly overestimated, the computed carrier densities using HSE06 is in good agreement with the available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electronic and optical properties of Zn-doped β-Ga 2 O 3 Czochralski single crystals

β-Ga 2 O 3 has several soluble deep acceptors that impart insulating behavior. In this study, we investigate Zn doping (0.25at.%) in bulk Czochralski and vertical gradient freeze β-Ga 2 O 3 . Representative crystals were assessed for orientation (electron backscatter diffraction and Raman spectroscopy), purity (glow discharge mass spectrometry and secondary ion mass spectrometry), optical properties (ultraviolet to near infrared absorption), and electrical properties (resistivity and current–voltage). Purity measurements indicate that Zn evaporation is insufficient to inhibit doping of Zn into β-Ga 2 O 3 . Hybrid functional calculations show Zn substitutes nearly equally on tetrahedral and octahedral sites, with less than ~0.1eV preference for the octahedral (Ga II ) site. Furthermore, calculations show that Zn Ga acts as a deep acceptor with trapping levels ~1.3 and ~0.9eV above the valence band for one and two holes, respectively. The solubility and electronic behavior of Zn dopants are consistent with measured concentrations >1×10 18 atoms/cm 3 and electrical measurements that show resistivity 10 11 –10 13 Ωcm, with no p-type conduction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ti- and Fe-related charge transition levels in $β$ – Ga 2 O 3

Here, deep-level transient spectroscopy measurements on $β$ – Ga 2 O 3 crystals reveal the presence of three defect signatures labeled $E$ 2a , $E$ 2b , and $E$ 3 with activation energies at around 0.66 eV, 0.73 eV, and 0.95 eV below the conduction band edge. Using secondary ion mass spectrometry, a correlation between the defect concentration associated with $E$ 3 and the Ti concentration present in the samples was found. Particularly, it is found that $E$ 3 is the dominant Ti-related defect in $β$ – Ga 2 O 3 and is associated with a single Ti atom. This finding is further corroborated by hybrid functional calculations that predict Ti substituting on an octahedral Ga site, denoted as Ti GaII , to be a good candidate for $E$ 3 . Moreover, the deep level transient spectroscopy results show that the level previously labeled $E$ 2 and attributed to Fe substituting on a gallium site (FeGa) consists of two overlapping signatures labeled E2a and E2b. We tentatively assign $E$ 2a and $E$ 2b to Fe substituting for Ga on a tetrahedral or an octahedral site, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of Heat Treatments in H 2 and Ar on the E 1 center in β -Gs 2 O 3

In this report the influence of heat treating n-type bulk β-Ga 2 O 3 in hydrogen (H2) and argon (Ar) gases on the presence of the defect level commonly labeled as E 1 was studied. Fourier transform-infrared spectroscopy confirms that hydrogen (H) is incorporated into β-Ga 2 O 3 during H 2 annealing at 900 °C. Deep-level transient spectroscopy measurements reveal that the concentration of the E 1 level is promoted by the introduction of H, in contrast to what is observed in samples heat-treated in an Ar flow. We further find the E 1 level to be stable against heat treatments at 650 K, both with and without an applied reverse-bias voltage. Potential candidates for the defect origin of E 1 are investigated using hybrid-functional calculations, and three types of defect complexes involving H are found to exhibit charge-state transition levels compatible with E 1 , including substitutional H at one of the threefold coordinated O sites, Ga-substitutional shallow donor impurities passivated by H, and certain configurations of singly hydrogenated Ga–O divacancies. Among these types, only the latter exhibit H binding energies that are consistent with the observed thermal stability of E 1 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso–meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4- phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi’s golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

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

Revisiting the link between magnetic properties and chemisorption at graphene nanoribbon zigzag edge

Graphene has received tremendous interest in both chemical and physical fields. Among different edges of the graphene system, the zigzag edge terminated graphene nanoribbons (ZGNRs) show unique magnetic properties in the antiferromagnetic (AFM) state. However, to date, the understanding of ZGNR chemical properties is mainly based on the partial radical concept, and in previous studies, the energy differences between the ferromagnetic (FM) and AFM states are smaller than experimental evidence. Here, we report that the strongly constrained and appropriately normed functional gives a significantly larger energy difference, which matches the experimental observation. Furthermore, utilizing the energetics in the large difference case, we propose a conceptual supplement to the previous partial radical concept: the overall stabilization of the AFM state compared to the nonmagnetic (NM) state consists of two parts that affect the adsorption energy conversely. The NM-FM energy differences will strengthen the adsorption, being in line with the previous partial radical concept. The FM-AFM energy differences will instead weaken the adsorption. Here, we perform calculations of H, OH, and LiS radical adsorption energies on ZGNRs to show that this weakening effect is numerically non-negligible: at least a ~0.2 eV difference in the adsorption energies is found. We expect that this refinement of the partial radical concept can provide a more comprehensive understanding of the chemical properties of ZGNRs. The differences in adsorption energies for the H, OH, and LiS radicals found here lead to significant changes in the predicted reactivity of the ZGNR models.

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 OSTI2022

Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions

Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu + , Zn + , Zn 2+ , and Mn + in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Furthermore, analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗