Within density functional theory (DFT), adding a Hubbard
Extensive Benchmarking of DFT+U Calculations for Predicting Band Gaps
Accurate computational predictions of band gaps are of practical importance to the modeling and development of semiconductor technologies, such as (opto)electronic devices and photoelectrochemical cells. Among available electronic-structure methods, density-functional theory (DFT) with the Hubbard U correction (DFT+U) applied to band edge states is a computationally tractable approach to improve the accuracy of band gap predictions beyond that of DFT calculations based on (semi)local functionals. At variance with DFT approximations, which are not intended to describe optical band gaps and other excited-state properties, DFT+U can be interpreted as an approximate spectral-potential method when U is determined by imposing the piecewise linearity of the total energy with respect to electronic occupations in the Hubbard manifold (thus removing self-interaction errors in this subspace), thereby providing a (heuristic) justification for using DFT+U to predict band gaps. However, it is still frequent in the literature to determine the Hubbard U parameters semiempirically by tuning their values to reproduce experimental band gaps, which ultimately alters the description of other total-energy characteristics. Here, we present an extensive assessment of DFT+U band gaps computed using self-consistent ab initio U parameters obtained from density-functional perturbation theory to impose the aforementioned piecewise linearity of the total energy. The study is carried out on 20 compounds containing transition-metal or p-block (group III-IV) elements, including oxides, nitrides, sulfides, oxynitrides, and oxysulfides. By comparing DFT+U results obtained using nonorthogonalized and orthogonalized atomic orbitals as Hubbard projectors, we find that the predicted band gaps are extremely sensitive to the type of projector functions and that the orthogonalized projectors give the most accurate band gaps, in satisfactory agreement with experimental data. This work demonstrates that DFT+U may serve as a useful method for high-throughput workflows that require reliable band gap predictions at moderate computational cost.
more »
« less
- Award ID(s):
- 1729338
- NSF-PAR ID:
- 10223643
- Date Published:
- Journal Name:
- Applied Sciences
- Volume:
- 11
- Issue:
- 5
- ISSN:
- 2076-3417
- Page Range / eLocation ID:
- 2395
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract U correction can mitigate some of the deficiencies of local and semi-local exchange-correlation functionals, while maintaining computational efficiency. However, the accuracy of DFT+U largely depends on the chosen HubbardU values. We propose an approach to determining the optimalU parameters for a given material by machine learning. The Bayesian optimization (BO) algorithm is used with an objective function formulated to reproduce the band structures produced by more accurate hybrid functionals. This approach is demonstrated for transition metal oxides, europium chalcogenides, and narrow-gap semiconductors. The band structures obtained using the BOU values are in agreement with hybrid functional results. Additionally, comparison to the linear response (LR) approach to determining U demonstrates that the BO method is superior. -
Oxides of p-block metals (e.g., indium oxide) and semimetals (e.g., antimony oxide) are of broad practical interest as transparent conductors and light absorbers for solar photoconversion due to the tunability of their electronic conductivity and optical absorption. Comparatively, these oxides have found limited applications in solar-to-hydrogen photocatalysis primarily due to their high electronegativity, which impedes electron transfer for converting protons into molecular hydrogen. We have shown recently that inserting s-block metal cations into p-block oxides is effective at lowering electronegativities while affording further control of band gaps. Here, we explain the origins of this dual tunability by demonstrating the mediator role of s-block metal cations in modulating orbital hybridization while not contributing to frontier electronic states. From this result, we carry out a comprehensive computational study of 109 ternary oxides of s- and p-block metal elements as candidate photocatalysts for solar hydrogen generation. We downselect the most desirable materials using band gaps and band edges obtained from Hubbard-corrected density-functional theory with Hubbard parameters computed entirely from first principles, evaluate the stability of these oxides in aqueous conditions, and characterize experimentally four of the remaining materials, synthesized with high phase uniformity, to assess the accuracy of computational predictions. We thus propose seven oxide semiconductors, including CsIn3O5, Sr2In2O5, and KSbO2 which, to the extent of our literature review, have not been previously considered as water-splitting photocatalysts.more » « less
-
A central problem of materials science is to determine whether a hypothetical material is stable without being synthesized, which is mathematically equivalent to a global optimization problem on a highly nonlinear and multimodal potential energy surface (PES). This optimization problem poses multiple outstanding challenges, including the exceedingly high dimensionality of the PES, and that PES must be constructed from a reliable, sophisticated, parameters-free, and thus very expensive computational method, for which density functional theory (DFT) is an example. DFT is a quantum mechanics-based method that can predict, among other things, the total potential energy of a given configuration of atoms. DFT, although accurate, is computationally expensive. In this work, we propose a novel expansion-exploration-exploitation framework to find the global minimum of the PES. Starting from a few atomic configurations, this “known” space is expanded to construct a big candidate set. The expansion begins in a nonadaptive manner, where new configurations are added without their potential energy being considered. A novel feature of this step is that it tends to generate a space-filling design without the knowledge of the boundaries of the domain space. If needed, the nonadaptive expansion of the space of configurations is followed by adaptive expansion, where “promising regions” of the domain space (those with low-energy configurations) are further expanded. Once a candidate set of configurations is obtained, it is simultaneously explored and exploited using Bayesian optimization to find the global minimum. The methodology is demonstrated using a problem of finding the most stable crystal structure of aluminum. History: Kwok Tsui served as the senior editor for this article. Funding: The authors acknowledge a U.S. National Science Foundation Grant DMREF-1921873 and XSEDE through Grant DMR170031. Data Ethics & Reproducibility Note: The code capsule is available on Code Ocean at https://codeocean.com/capsule/3366149/tree and in the e-Companion to this article (available at https://doi.org/10.1287/ijds.2023.0028 ).more » « less
-
null (Ed.)Achieving a molecular-level understanding of how the structures and compositions of metal–organic frameworks (MOFs) influence their charge carrier concentration and charge transport mechanism—the two key parameters of electrical conductivity—is essential for the successful development of electrically conducting MOFs, which have recently emerged as one of the most coveted functional materials due to their diverse potential applications in advanced electronics and energy technologies. Herein, we have constructed four new alkali metal (Na, K, Rb, and Cs) frameworks based on an electron-rich tetrathiafulvalene tetracarboxylate (TTFTC) ligand, which formed continuous π-stacks, albeit with different π–π-stacking and S⋯S distances ( d π–π and d S⋯S ). These MOFs also contained different amounts of aerobically oxidized TTFTC˙ + radical cations that were quantified by electron spin resonance (ESR) spectroscopy. Density functional theory calculations and diffuse reflectance spectroscopy demonstrated that depending on the π–π-interaction and TTFTC˙ + population, these MOFs enjoyed varying degrees of TTFTC/TTFTC˙ + intervalence charge transfer (IVCT) interactions, which commensurately affected their electronic and optical band gaps and electrical conductivity. Having the shortest d π–π (3.39 Å) and the largest initial TTFTC˙ + population (∼23%), the oxidized Na-MOF 1-ox displayed the narrowest band gap (1.33 eV) and the highest room temperature electrical conductivity (3.6 × 10 −5 S cm −1 ), whereas owing to its longest d π–π (3.68 Å) and a negligible TTFTC˙ + population, neutral Cs-MOF 4 exhibited the widest band gap (2.15 eV) and the lowest electrical conductivity (1.8 × 10 −7 S cm −1 ). The freshly prepared but not optimally oxidized K-MOF 2 and Rb-MOF 3 initially displayed intermediate band gaps and conductivity, however, upon prolonged aerobic oxidation, which raised the TTFTC˙ + population to saturation levels (∼25 and 10%, respectively), the resulting 2-ox and 3-ox displayed much narrower band gaps (∼1.35 eV) and higher electrical conductivity (6.6 × 10 −5 and 4.7 × 10 −5 S cm −1 , respectively). The computational studies indicated that charge movement in these MOFs occurred predominantly through the π-stacked ligands, while the experimental results displayed the combined effects of π–π-interactions, TTFTC˙ + population, and TTFTC/TTFTC˙ + IVCT interaction on their electronic and optical properties, demonstrating that IVCT interactions between the mixed-valent ligands could be exploited as an effective design strategy to develop electrically conducting MOFs.more » « less
-
null (Ed.)The production of hydrogen fuels, via water splitting, is of practical relevance for meeting global energy needs and mitigating the environmental consequences of fossil-fuel-based transportation. Water photoelectrolysis has been proposed as a viable approach for generating hydrogen, provided that stable and inexpensive photocatalysts with conversion efficiencies over 10% can be discovered, synthesized at scale, and successfully deployed (Pinaud et al. , Energy Environ. Sci. , 2013, 6 , 1983). While a number of first-principles studies have focused on the data-driven discovery of photocatalysts, in the absence of systematic experimental validation, the success rate of these predictions may be limited. We address this problem by developing a screening procedure with co-validation between experiment and theory to expedite the synthesis, characterization, and testing of the computationally predicted, most desirable materials. Starting with 70 150 compounds in the Materials Project database, the proposed protocol yielded 71 candidate photocatalysts, 11 of which were synthesized as single-phase materials. Experiments confirmed hydrogen generation and favorable band alignment for 6 of the 11 compounds, with the most promising ones belonging to the families of alkali and alkaline-earth indates and orthoplumbates. This study shows the accuracy of a nonempirical, Hubbard-corrected density-functional theory method to predict band gaps and band offsets at a fraction of the computational cost of hybrid functionals, and outlines an effective strategy to identify photocatalysts for solar hydrogen generation.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/app11052395
