Kohn-Sham (KS) density functional theory (DFT) is a very efficient method for calculating various properties of solids as, for instance, the total energy, the electron density, or the electronic band structure. The KS-DFT method leads to rather fast calculations, however the accuracy depends crucially on the chosen approximation for the exchange and correlation (xc) functional E xc and/or potential v xc . Here, an overview of xc methods to calculate the electronic band structure is given, with the focus on the so-called semilocal methods that are the fastest in KS-DFT and allow to treat systems containing up to thousands of atoms. Among them, there is the modified Becke-Johnson potential that is widely used to calculate the fundamental band gap of semiconductors and insulators. The accuracy for other properties like the magnetic moment or the electron density, that are also determined directly by v xc , is also discussed.

The density-functional theory (DFT) approximations that are the most accurate for the calculation of bandgap of bulk materials are hybrid functionals, such as HSE06, the modified Becke–Johnson (MBJ) potential, and the GLLB-SC potential. More recently, generalized gradient approximations (GGAs), such as HLE16, or meta-GGAs, such as (m)TASK, have also proven to be quite accurate for the bandgap. Here, the focus is on two-dimensional (2D) materials and the goal is to provide a broad overview of the performance of DFT functionals by considering a large test set of 298 2D systems. The present work is an extension of our recent studies [T. Rauch, M. A. L. Marques, and S. Botti, Phys. Rev. B 101, 245163 (2020); Patra et al., J. Phys. Chem. C 125, 11206 (2021)]. Due to the lack of experimental results for the bandgap of 2D systems, G0W0 results were taken as reference. It is shown that the GLLB-SC potential and mTASK functional provide the bandgaps that are the closest to G0W0. Following closely, the local MBJ potential has a pretty good accuracy that is similar to the accuracy of the more expensive hybrid functional HSE06.

Nowadays pseudopotential (PP) density functional theory calculations constitute the standard approach to tackle solid-state electronic problems. These rely on distributed PP tables that were built from all-electron atomic calculations using few popular semilocal exchange-correlation functionals, while PPs based on more modern functionals, such as meta-generalized gradient approximation and hybrid functionals, or for many-body methods, such as GW, are often not available. Because of this, employing PPs created with inconsistent exchange-correlation functionals has become a common practice. Our aim is to quantify systematically the error in the determination of the electronic band gap when cross-functional PP calculations are performed. To this end, we compare band gaps obtained with norm-conserving PPs or the projector-augmented wave method with all-electron calculations for a large data set of 473 solids. We focus, in particular, on density functionals that were designed specifically for band gap calculations. On average, the absolute error is about 0.1 eV, yielding absolute relative errors in the 5–10% range. Considering that typical errors stemming from the choice of the functional are usually larger, we conclude that the effect of choosing an inconsistent PP is rather harmless for most applications. However, we find specific cases where absolute errors can be larger than 1 eV or others where relative errors can amount to a large fraction of the band gap.

The experimental and theoretical realization of two-dimensional (2D) materials is of utmost importance in semiconducting applications. Computational modeling of these systems with satisfactory accuracy and computational efficiency is only feasible with semilocal density functional theory methods. In the search for the most useful method in predicting the band gap of 2D materials, we assess the accuracy of recently developed semilocal exchange–correlation (XC) energy functionals and potentials. Though the explicit forms of exchange–correlation (XC) potentials are very effective against XC energy functionals for the band gap of bulk solids, their performance needs to be investigated for 2D materials. In particular, the LMBJ [ 32097004 J. Chem. Theory Comput. 2020 16 2654 ] and GLLB-SC [ 115106 Phys. Rev. B 2010 82 ] potentials are considered for their dominance in bulk band gap calculation. The performance of recently developed meta generalized gradient approximations, like TASK [ 033082 Phys. Rev. Res. 2019 1 ] and MGGAC [ 155140 Phys. Rev. B 2019 100 ], is also assessed. We find that the LMBJ potential constructed for 2D materials is not as successful as its parent functional, i.e., MBJ [ 226401 19658882 Phys. Rev. Lett. 2009 102 ] in bulk solids. Due to a contribution from the derivative discontinuity, the band gaps obtained with GLLB-SC are in a certain number of cases, albeit not systematically, larger than those obtained with other methods, which leads to better agreement with the quasi-particle band gap obtained from the GW method. The band gaps obtained with TASK and MGGAC can also be quite accurate.

The DFT-1/2 method in density functional theory [L. G. Ferreira et al., Phys. Rev. B 78, 125116 (2008)] aims to provide accurate band gaps at the computational cost of semilocal calculations. The method has shown promise in a large number of cases, however some of its limitations or ambiguities on how to apply it to covalent semiconductors have been pointed out recently [K.-H. Xue et al., Comput. Mater. Science 153, 493 (2018)]. In this work, we investigate in detail some of the problems of the DFT-1/2 method with a focus on two classes of materials: covalently bonded semiconductors and transition-metal oxides. We argue for caution in the application of DFT-1/2 to these materials, and the condition to get an improved band gap is a spatial separation of the orbitals at the valence band maximum and conduction band minimum.