skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on March 6, 2026

Title: Integrating density functional theory with machine learning for enhanced band gap prediction in metal oxides
In this study, we used a combination of density functional theory with Hubbard U correction (DFT+U) and machine learning (ML) to accurately predict the band gaps and lattice parameters of metal oxides: TiO2 (rutile and anatase), cubic ZnO, cubic ZnO2, cubic CeO2, and cubic ZrO2. Our results show that including Up values for oxygen 2p orbitals alongside Ud/f for metal 3d or 4f orbitals significantly enhances the accuracy of these predictions. Through extensive DFT+U calculations, we identify optimal (Up, Ud/f) integer pairs that closely reproduce experimentally measured band gaps and lattice parameters for each oxide: (8 eV, 8 eV) for rutile TiO2; (3 eV, 6 eV) for anatase TiO2; (6 eV, 12 eV) for c-ZnO; (10 eV, 10 eV) for c-ZnO2; (9 eV, 5 eV) for c-ZrO2; and (7 eV, 12 eV) for c-CeO2. Our ML analysis showed that simple supervised ML models can closely reproduce these DFT+U results at a fraction of the computational cost and generalize well to related polymorphs. Our approach builds on existing high-throughput DFT+U frameworks by providing fast pre-DFT estimates of structural properties and band gaps. Since this work does not aim to improve the underlying DFT+U method, the ML model shares its limitations. We also note that the reported values of Up strongly depend on the choice of correlated orbitals, and caution is recommended with a different choice of correlated orbitals.  more » « less
Award ID(s):
2503612 2055012
PAR ID:
10584845
Author(s) / Creator(s):
; ;
Publisher / Repository:
Royal Society of Chemistry
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
27
Issue:
10
ISSN:
1463-9076
Page Range / eLocation ID:
5338 to 5358
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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
  2. The amorphous-to-crystalline transition in TiO2 films upon annealing in air yields different polymorphs, depending on the oxygen partial pressure during the deposition of the amorphous precursor film. We further manipulate the resulting polymorph by introducing Sn into the system. By depositing a few nanometer-thick layer of metallic Sn between two layers of amorphous TiO2 prepared to yield the anatase polymorph of TiO2, we find that it results in the rutile polymorph if the content of Sn is high enough. If Sn is introduced as an oxide, no rutile is formed; anatase is by far the predominant phase (with a very small amount of brookite). This observation is consistent with scavenging of oxygen by elemental Sn at the Sn/TiO2 interfaces, stabilizing the rutile structure that can accommodate oxygen vacancies. 
    more » « less
  3. In standard doping, adding charge carrier to a compound results in a shift of the Fermi level towards the conduction band for electron doping and towards the valence band for hole doping. We discuss the curious case of antidoping, where the direction of band movements in response to doping is reversed. Specifically, p-type antidoping moves the previously occupied bands to the principal conduction band resulting in an increase of band gap energy and reduction of electronic conductivity. We find that this is a generic behavior for a class of materials: early transition and rare earth metal (e.g., Ti, Ce) oxides where the sum of composition-weighed formal oxidation states is positive; such compounds tend to form the well-known electron-trapped intermediate bands localized on the reduced cation orbitals. What is less known is that doping by a hole annihilates a single trapped electron on a cation. The latter thus becomes electronically inequivalent with respect to the normal cation in the undoped lattice, thus representing a symmetry-breaking effect. We give specific theoretical predictions for target compounds where hole antidoping might be observed experimentally: Magnéli-like phases (i.e., CeO2-x and TiO2-x) and ternary compounds (i.e., Ba2Ti6O13 and Ba4Ti12O27), and note that this unique behavior opens the possibility of unconventional control of materials conductivity by doping. 
    more » « less
  4. Abstract Kohn–Sham density functional theory (DFT)‐based searches for hypothetical catalysts are too computationally demanding for wide searches through diverse materials space. Here, the accuracy of computational alchemy schemes on carbides, nitrides, and oxides is assessed. With a single set of reference DFT calculations, computational alchemy approximates adsorbate binding energies (BEs) on a large number of hypothetical catalysts surfaces with negligible computational cost. Analogous to previous studies on metal alloys, computational alchemy predicts adsorbate BEs on rocksalt TiC(111), TiN(100), and TiO(100) materials, which have no bandgap, in close agreement with DFT results (with mean unsigned errors up to 0.33 eV). In contrast, it is found that semiconducting systems such as rutile TiO2(110), rutile SnO2(110), and rocksalt ZnO(100) can present more significant challenges. This work identifies these challenges being linked to the density of states at the Fermi level and by adding Pt dopants in the surface layer of TiO2, it is shown that computational alchemy can become more reliable with non‐transition metal systems. This remedy provides insight that promotes computational alchemy for broad searches for catalyst active sites through materials space beyond transition metal alloys. 
    more » « less
  5. TiO 2 has been extensively studied in many fields including photocatalysis, electrochemistry, optics, etc. Understanding the mechanism of the anatase–rutile phase transition (ART) process is critical for the design of TiO 2 -based high-activity photocatalysts and tuning its properties for other applications. In this work, the ART process using individual anatase micro-particles with a large percentage of (001) facets was monitored and studied. Phase concentration evolution obtained via Raman microscopy was correlated with the morphological evolution observed in scanning electron microscope (SEM) images. The ART of anatase microcrystals is dominated by surface nucleation and growth, but the ART processes of individual anatase particles are distinctive and depend on the various rutile nucleation sites. Two types of transformation pathways are observed. In one type of ART pathway, the rutile phase nucleated at a corner of an anatase microcrystal and grew in one direction along the edge of the crystal firstly followed by propagation over the rest of the microcrystal in the orthogonal direction on the surface and to the bulk of the crystal. The kinetics of the ART follows the first-order model with two distinct rate constants. The fast reaction rate is from the surface nucleation and growth, and the slow rate is from the bulk nucleation and growth. In the other type of ART pathway, multiple rutile nucleation sites formed simultaneously on different edges and corners of the microcrystal. The rutile phase spread over the whole crystal from these nucleation sites with a small contribution of bulk nucleation. Our study on the ART of individual micro-sized crystals bridges the material gap between bulk crystals and nano-sized TiO 2 particles. The anatase/rutile co-existing particle will provide a perfect platform to study the synergistic effect between the anatase phase and the rutile phase in their catalytic performances. 
    more » « less