In the context of the modified Becke-Johnson (mBJ) potential, we recently underlined that $bar{g}$, the average of $leftvert ablarhorightvert/rho$ in the unit cell, has markedly different values in transition-metal oxides and pure transition metals [
Tran et al., J. Appl. Phys. 126, 110902 (2019)]. However, since $bar{g}$ is a constant it is not able to provide local information about a particular atom in the system. Furthermore, while $overline{g}$ can be used only for periodic bulk solids, a local (i.e., position-dependent) version would allow us to consider also low-dimensional systems and interfaces. Such a local function has been proposed by Rauch et al. [J. Chem. Theory Comput. 16, 2654 (2020)] for the local mBJ potential. Actually, a local version of $overline{g}$, or of another similar quantity like the reduced density gradient $overline{s}$, could also be used in the framework of other methods. Here, we explored the idea to use such a local function $tilde{g}$ (or $tilde{s}$), defined as the average of $g$ (or $s$) over a certain region around a transition-metal atom, to estimate the degree of on-site correlation on this atom. We found a large difference in our correlation estimators between non-correlated and correlated materials, proving its usefulness and reliability. Our estimators can subsequently be used to determine whether or not a Hubbard $U$ on-site correction in the DFT+$U$ method should be applied to a particular atom. This is particularly interesting in cases where the degree of correlation of the transition-metal atoms is not clear, like interfaces between correlated and non-correlated materials or oxygen-covered metal surfaces. In such cases, our estimators could also be used for an interpolation of $U$ between correlated and non-correlated atoms.
The density functional theory (DFT) approximations that are the most accurate for the calculation of band gap of bulk materials are hybrid functionals like HSE06, the MBJ potential, and the GLLB-SC potential. More recently, generalized gradient appro
ximations (GGA), like HLE16, or meta-GGAs, like (m)TASK, have proven to be also quite accurate for the band gap. Here, the focus is on 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 layered systems. The present work is an extension of our recent studies [Rauch et al., Phys. Rev. B 101, 245163 (2020) and Patra et al., J. Phys. Chem. C 125, xxxxx (2021)]. Due to the lack of experimental results for the band gap of 2D systems, $G_{0}W_{0}$ results were taken as reference. It is shown that the GLLB-SC potential and mTASK functional provide the band gaps that are the closest to $G_{0}W_{0}$. Following closely, the local MBJ potential has a pretty good accuracy that is similar to the accuracy of the more expensive hybrid functional HSE06.
A degenerate perturbation $kcdot p$ approach for effective mass calculations is implemented in the all-electron density functional theory (DFT) package WIEN2k. The accuracy is tested on major group IVA, IIIA-VA, and IIB-VIA semiconductor materials. T
hen, the effective mass in graphene and CuI with defects is presented as illustrative applications. For states with significant Cu-d character additional local orbitals with higher principal quantum numbers (more radial nodes) have to be added to the basis set in order to converge the results of the perturbation theory. Caveats related to a difference between velocity and momentum matrix elements are discussed in the context of application of the method to non-local potentials, such as Hartree-Fock/DFT hybrid functionals and DFT+U.
Several recent studies have shown that SCAN, a functional belonging to the meta-generalized gradient approximation (MGGA) family, leads to significantly overestimated magnetic moments in itinerant ferromagnetic metals. However, this behavior is not i
nherent to the MGGA level of approximation since TPSS, for instance, does not lead to such severe overestimations. In order to provide a broader view of the accuracy of MGGA functionals for magnetism, we extend the assessment to more functionals, but also to antiferromagnetic solids. The results show that to describe magnetism there is overall no real advantage in using a MGGA functional compared to GGAs. For both types of approximation, an improvement in ferromagnetic metals is necessarily accompanied by a deterioration (underestimation) in antiferromagnetic insulators, and vice-versa. We also provide some analysis in order to understand in more detail the relation between the mathematical form of the functionals and the results.
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 functional
s are also potentially more interesting for the electronic structure, in particular when the potential is non-multiplicative (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 be done very accurately also 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.
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_{text{xc}}$ and/or potential $v_{text{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_{text{xc}}$, is also discussed.
The nonlocal van der Waals (NL-vdW) functionals [Dion et al., Phys. Rev. Lett. 92, 246401 (2004)] are being applied more and more frequently in solid-state physics, since they have shown to be much more reliable than the traditional semilocal functio
nals for systems where weak interactions play a major role. However, a certain number of NL-vdW functionals have been proposed during the last few years, such that it is not always clear which one should be used. In this work, an assessment of NL-vdW functionals is presented. Our test set consists of weakly bound solids, namely rare gases, layered systems like graphite, and molecular solids, but also strongly bound solids in order to provide a more general conclusion about the accuracy of NL-vdW functionals for extended systems. We found that among the tested functionals, rev-vdW-DF2 [Hamada, Phys. Rev. B 89, 121103(R) (2014)] is very accurate for weakly bound solids, but also quite reliable for strongly bound solids.
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, h
owever 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.
An alternative type of approximation for the exchange and correlation functional in density functional theory is proposed. This approximation depends on a variable $u$ that is able to detect inhomogeneities in the electron density $rho$ without using
derivatives of $rho$. Instead, $u$ depends on the orbital energies which can also be used to measure how a system differs from the homogeneous electron gas. Starting from the functional of Perdew, Burke, and Ernzerhof (PBE) [Phys. Rev. Lett. 77, 3865 (1996)], a functional depending on $u$ is constructed. Tests on the lattice constant, bulk modulus, and cohesive energy of solids show that this $u$-dependent PBE-like functional is on average as accurate as the original PBE or its solid-state version PBEsol. Since $u$ carries more nonlocality than the reduced density gradient $s$ used in functionals of the generalized gradient approximation (GGA) like PBE and $alpha$ used in meta-GGAs, it will be certainly useful for the future development of more accurate exchange-correlation functionals.
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 band gap, calculated with the modified Becke-Johnson MGGA potential, the deorbitalization has a larger impact on the results.
