What is the optimal mGGA exchange functional for solids?

The space of generalized gradient approximation (GGA) and meta-GGA (mGGA) exchange approximations is systematically explored by training 25 new functionals to produce accurate lattice parameter, cohesive energy, and bandgap predictions. The trained functionals are used to reproduce exact constraints in a data-driven way and to understand the accuracy trade-off between the mentioned properties. The functionals are compared to notable mGGA functionals to analyze how changes in the enhancement factor maps influence the accuracy of predictions. Some of the trained functionals are found to perform on par with specialized functionals for bandgaps, while outperforming them on the other two properties. The error surface of our trained functionals can serve as a soft-limit of what mGGA functionals can achieve.

Efficient and improved prediction of the band offsets at semiconductor heterojunctions from meta-GGA density functionals: A benchmark study.

Accurate theoretical prediction of the band offsets at interfaces of semiconductor heterostructures can often be quite challenging. Although density functional theory has been reasonably successful to carry out such calculations, efficient, accurate semilocal functionals are desirable to reduce the computational cost. In general, the semilocal functionals based on the generalized gradient approximation (GGA) significantly underestimate the bulk bandgaps. This, in turn, results in inaccurate estimates of the band offsets at the heterointerfaces. In this paper, we investigate the performance of several advanced meta-GGA functionals in the computational prediction of band offsets at semiconductor heterojunctions. In particular, we investigate the performance of r2SCAN (two times revised strongly constrained and appropriately normed functional), rMGGAC (revised semilocal functional based on cuspless hydrogen model and Pauli kinetic energy density functional), mTASK (modified Aschebrock and Kümmel meta-GGA functional), and local modified Becke-Johnson exchange-correlation functionals. Our results strongly suggest that these meta-GGA functionals for supercell calculations perform quite well, especially, when compared to computationally more demanding GW calculations. We also present band offsets calculated using ionization potentials and electron affinities, as well as band alignment via the branch point energies. Overall, our study shows that the aforementioned meta-GGA functionals can be used within the density functional theory framework to estimate the band offsets in semiconductor heterostructures with predictive accuracy.

Bandgap of two-dimensional materials: Thorough assessment of modern exchange-correlation functionals.

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.

WIEN2k: An APW+lo program for calculating the properties of solids.

The WIEN2k program is based on the augmented plane wave plus local orbitals (APW+lo) method to solve the Kohn-Sham equations of density functional theory. The APW+lo method, which considers all electrons (core and valence) self-consistently in a full-potential treatment, is implemented very efficiently in WIEN2k, since various types of parallelization are available and many optimized numerical libraries can be used. Many properties can be calculated, ranging from the basic ones, such as the electronic band structure or the optimized atomic structure, to more specialized ones such as the nuclear magnetic resonance shielding tensor or the electric polarization. After a brief presentation of the APW+lo method, we review the usage, capabilities, and features of WIEN2k (version 19) in detail. The various options, properties, and available approximations for the exchange-correlation functional, as well as the external libraries or programs that can be used with WIEN2k, are mentioned. References to relevant applications and some examples are also given.

Comparative study of the PBE and SCAN functionals: The particular case of alkali metals.

The SCAN meta-generalized gradient approximation (GGA) functional is known to describe multiple properties of various materials with different types of bonds with greater accuracy, compared to the widely used PBE GGA functional. Yet, for alkali metals, SCAN shows worse agreement with experimental results than PBE despite using more information about the system. In the current study, this behavior for alkali metals is explained by identifying an inner semicore region which, within SCAN, contributes to an underbinding. The inner semicore push toward larger lattice constants is a general feature but is particularly important for very soft materials, such as the alkali metals, while for harder materials the valence region dominates.

On the calculation of the bandgap of periodic solids with MGGA functionals using the total energy.

During the last few years, it has become more and more clear that functionals of the meta generalized gradient approximation (MGGA) are more accurate than GGA functionals for the geometry and energetics of electronic systems. However, MGGA functionals are also potentially more interesting for the electronic structure, in particular, when the potential is nonmultiplicative (i.e., when MGGAs are implemented in the generalized Kohn-Sham framework), which may help to get more accurate bandgaps. Here, we show that the calculation of bandgap of solids with MGGA functionals can also be done very accurately in a non-self-consistent manner. This scheme uses only the total energy and can, therefore, be very useful when the self-consistent implementation of a particular MGGA functional is not available. Since self-consistent MGGA calculations may be difficult to converge, the non-self-consistent scheme may also help to speed up the calculations. Furthermore, it can be applied to any other types of functionals, for which the implementation of the corresponding potential is not trivial.

Orbital-free approximations to the kinetic-energy density in exchange-correlation MGGA functionals: Tests on solids.

A recent study of Mejia-Rodriguez and Trickey [Phys. Rev. A 96, 052512 (2017)] showed that the deorbitalization procedure (replacing the exact Kohn-Sham kinetic-energy density by an approximate orbital-free expression) applied to exchange-correlation functionals of the meta-generalized gradient approximation (MGGA) can lead to important changes in the results for molecular properties. For the present work, the deorbitalization of MGGA functionals is further investigated by considering various properties of solids. It is shown that depending on the MGGA, common orbital-free approximations to the kinetic-energy density can be sufficiently accurate for the lattice constant, bulk modulus, and cohesive energy. For the bandgap, calculated with the modified Becke-Johnson MGGA potential, the deorbitalization has a larger impact on the results.

Importance of the Kinetic Energy Density for Band Gap Calculations in Solids with Density Functional Theory.

Recently, exchange-correlation potentials in density functional theory were developed with the goal of providing improved band gaps in solids. Among them, the semilocal potentials are particularly interesting for large systems since they lead to calculations that are much faster than with hybrid functionals or methods like GW. We present an exhaustive comparison of semilocal exchange-correlation potentials for band gap calculations on a large test set of solids, and particular attention is paid to the potential HLE16 proposed by Verma and Truhlar. It is shown that the most accurate potential is the modified Becke-Johnson potential, which, most noticeably, is much more accurate than all other semilocal potentials for strongly correlated systems. This can be attributed to its additional dependence on the kinetic energy density. It is also shown that the modified Becke-Johnson potential is at least as accurate as the hybrid functionals and more reliable for solids with large band gaps.

Rungs 1 to 4 of DFT Jacob's ladder: Extensive test on the lattice constant, bulk modulus, and cohesive energy of solids.

A large panel of old and recently proposed exchange-correlation functionals belonging to rungs 1 to 4 of Jacob's ladder of density functional theory are tested (with and without a dispersion correction term) for the calculation of the lattice constant, bulk modulus, and cohesive energy of solids. Particular attention will be paid to the functionals MGGA_MS2 [J. Sun et al., J. Chem. Phys. 138, 044113 (2013)], mBEEF [J. Wellendorff et al., J. Chem. Phys. 140, 144107 (2014)], and SCAN [J. Sun et al., Phys. Rev. Lett. 115, 036402 (2015)] which are meta-generalized gradient approximations (meta-GGA) and are developed with the goal to be universally good. Another goal is also to determine for which semilocal functionals and groups of solids it is beneficial (or not necessary) to use the Hartree-Fock exchange or a dispersion correction term. It is concluded that for strongly bound solids, functionals of the GGA, i.e., rung 2 of Jacob's ladder, are as accurate as the more sophisticated functionals of the higher rungs, while it is necessary to use dispersion corrected functionals in order to expect at least meaningful results for weakly bound solids. If results for finite systems are also considered, then the meta-GGA functionals are overall clearly superior to the GGA functionals.

Site-selective adsorption of phthalocyanine on h-BN/Rh(111) nanomesh.

Experiment and computer simulations were conducted in order to study the adsorption of the phthalocyanine molecules H2Pc and CuPc on the h-BN/Rh(111) nanomesh. We combine STM investigations with the exploration of the potential energy surface as resulting from density functional theory calculations. Both approaches indicate a pronounced adsorption selectivity in the so called pore regions of the h-BN nanomesh, whereas the adsorption energy landscape in the pore turns out to be very shallow. This is seen by the inability to image the molecule stably at 77 K by scanning tunneling microscopy. Understanding the nature of the binding by rationalizing the site-selectivity and the mobility of the molecules is quite a challenge for both experiment and theory. In particular, we observe that the choice of the functional in the DFT description is crucial to be able to discriminate among adsorption sites that are very close in energy and to resolve low energy barriers. Our study reveals how the shape of the corrugated h-BN layer is the dominant factor that determines the subtle features of the potential energy surface for the adsorption of phthalocyanine.

Nonlocal van der Waals functionals: the case of rare-gas dimers and solids.

Recently, the nonlocal van der Waals (vdW) density functionals [M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)] have attracted considerable attention due to their good performance for systems where weak interactions are important. Since the physics of dispersion is included in these functionals, they are usually more accurate and show less erratic behavior than the semilocal and hybrid methods. In this work, several variants of the vdW functionals have been tested on rare-gas dimers (from He2 to Kr2) and solids (Ne, Ar, and Kr) and their accuracy compared to standard semilocal approximations, supplemented or not by an atom-pairwise dispersion correction [S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, J. Chem. Phys. 132, 154104 (2010)]. An analysis of the results in terms of energy decomposition is also provided.

Similarity Clustering for Representative Sets of Inorganic Solids for Density Functional Testing.

Benchmarking DFT functionals is complicated since the results highly depend on which properties and materials were used in the process. Unwanted biases can be introduced if a data set contains too many examples of very similar materials. We show that a clustering based on the distribution of density gradient and kinetic energy density is able to identify groups of chemically distinct solids. We then propose a method to create smaller data sets or rebalance existing data sets in a way that no region of the meta-GGA descriptor space is overrepresented, yet the new data set reproduces average errors of the original set as closely as possible. We apply the method to an existing set of 44 inorganic solids and suggest a representative set of seven solids. The representative sets generated with this method can be used to make more general benchmarks or to train new functionals.

Efficient Band Structure Calculation of Two-Dimensional Materials from Semilocal Density Functionals.

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 [J. Chem. Theory Comput.2020, 16, 2654] and GLLB-SC [Phys. Rev. B82, 2010, 115106] potentials are considered for their dominance in bulk band gap calculation. The performance of recently developed meta generalized gradient approximations, like TASK [Phys. Rev. Res.1, 2019, 033082] and MGGAC [Phys. Rev. B. 100, 2019, 155140], is also assessed. We find that the LMBJ potential constructed for 2D materials is not as successful as its parent functional, i.e., MBJ [Phys. Rev. Lett.102, 2009, 226401] 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.

Validation of Pseudopotential Calculations for the Electronic Band Gap of Solids.

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.

Large-Scale Benchmark of Exchange-Correlation Functionals for the Determination of Electronic Band Gaps of Solids.

We compile a large data set designed for the efficient benchmarking of exchange-correlation functionals for the calculation of electronic band gaps. The data set comprises information on the experimental structure and band gap of 472 nonmagnetic materials and includes a diverse group of covalent-, ionic-, and van der Waals-bonded solids. We used it to benchmark 12 functionals, ranging from standard local and semilocal functionals, passing through meta-generalized-gradient approximations, and several hybrids. We included both general purpose functionals, like the Perdew-Burke-Ernzerhof approximation, and functionals specifically crafted for the determination of band gaps. The comparison of experimental and theoretical band gaps shows that the modified Becke-Johnson is at the moment the best available density functional, closely followed by the Heyd-Scuseria-Ernzerhof screened hybrid from 2006 and the high-local-exchange generalized-gradient approximation.

How close are the Slater and Becke-Roussel potentials in solids?

The Becke­Roussel (BR) potential [Phys. Rev. A 1989, 39, 3761] was proposed as an approximation to the Slater potential, which is the Coulomb potential generated by the exact exchange hole. In the present work, a detailed comparison between the Slater and BR potentials in solids is presented. It is shown that the two potentials usually lead to very similar results for the electronic structure; however, in a few cases, e.g., Si, Ge, or strongly correlated systems like NiO, the fundamental band gap or magnetic properties can differ markedly. Such differences should not be neglected when the computationally expensive Slater potential is replaced by the cheap semilocal BR potential in approximations to the exact-exchange Kohn­Sham potential, such as the one proposed by Becke and Johnson [J. Chem. Phys. 2006, 124, 221101].

Hybrid functionals for solids with an optimized Hartree-Fock mixing parameter.

(Screened) hybrid functionals are being used more and more for solid-state calculations. Usually the fraction α of Hartree-Fock exchange is kept fixed during the calculation; however, there is no single (universal) value for α which systematically leads to satisfying accuracy. Instead, one could use a property of the system under consideration to determine α, and in this way the functional would be more flexible and potentially more accurate. Recently, it was proposed to use the static dielectric constant ε for the calculation of α (Shimazaki and Asai 2008 Chem. Phys. Lett. 466 91 and Marques et al 2011 Phys. Rev. B 83 035119). We explore this idea further and propose a scheme where the connection between ε and α is optimized based on experimental band gaps. ε, and thus α, is recalculated at each iteration of the self-consistent procedure. We present results for the bandgap and lattice constant of various semiconductors and insulators with this procedure. In addition, we show that this approach can also be combined with a non-self-consistent hybrid approximation to speed up the calculations considerably, while retaining an excellent accuracy in most cases.

Accurate band gaps of semiconductors and insulators with a semilocal exchange-correlation potential.

A modified version of the exchange potential proposed by Becke and Johnson [J. Chem. Phys. 124, 221101 (2006)10.1063/1.2213970] is tested on solids for the calculation of band gaps. The agreement with experiment is very good for all types of solids we considered (e.g., wide band gap insulators, sp semiconductors, and strongly correlated 3d transition-metal oxides) and is of the same order as the agreement obtained with the hybrid functionals or the GW methods. This semilocal exchange potential, which recovers the local-density approximation (LDA) for a constant electron density, mimics very well the behavior of orbital-dependent potentials and leads to calculations which are barely more expensive than LDA calculations. Therefore, it can be applied to very large systems in an efficient way.