Partially deorbitalized meta-GGA

Mejia-Rodriguez and Trickey recently proposed a procedure for removing the explicit dependence of meta-GGA exchange-correlation energy functionals Exc on the kinetic energy density τ. We present a simple modification to this approach in which the exact Kohn-Sham τ is used as input for Exc but the functional derivative of τ with respect to the density ρ, required to calculate the potential term ∫d3r′δExc∕δτ(r′)∣ρ⋅δτ(r′)∕δρ(r), is evaluated using an approximate kinetic energy density functional. This ‘half-way’ strategy ensures that the Kohn-Sham potential is a local multiplicative function (as opposed to the non-local potential of a generalized Kohn-Sham approach) while preserving the inherent non-locality of the functional itself. Electronic structure codes can be easily modified to use the new method. We validate it by quantifying the accuracy of the predicted lattice parameters, bulk moduli, magnetic moments and cohesive energies of a large set of periodic solids. An unanticipated benefit of this method is to gauge the quality of approximate kinetic energy functionals by checking if the self-consistent solution is indeed at the variational minimum.