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

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r2SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew-Burke-Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r2SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

Kinetic Energy Density Functionals Based on a Generalized Screened Coulomb Potential: Linear Response and Future Perspectives

We consider kinetic energy functionals that depend, beside the usual semilocal quantities (density, gradient, Laplacian of the density), on a generalized Yukawa potential, that is the screened Coulomb potential of the density raised to some power. These functionals, named Yukawa generalized gradient approximations (yGGA), are potentially efficient real-space semilocal methods that include significant non-local effects and can describe different important exact properties of the kinetic energy. In this work, we focus in particular on the linear response behavior for the homogeneous electron gas (HEG). We show that such functionals are able to reproduce the exact Lindhard function behavior with a very good accuracy, outperforming all other semilocal kinetic functionals. These theoretical advances allow us to perform a detailed analysis of a special class of yGGAs, namely the linear yGGA functionals. Thus, we show how the present approach can generalize the yGGA functionals improving the HEG linear behavior and leading to an extended formula for the kinetic functional. Moreover, testing on several jellium cluster model systems allows highlighting advantages and limitations of the linear yGGA functionals and future perspectives for the development of yGGA kinetic functionals.

Density functionals with asymptotic-potential corrections are required for the simulation of spectroscopic properties of materials

Five effects of correction of the asymptotic potential error in density functionals are identified that significantly improve calculated properties of molecular excited states involving charge-transfer character. Newly developed materials-science computational methods are used to demonstrate how these effects manifest in materials spectroscopy. Connection is made considering chlorophyll-a as a paradigm for molecular spectroscopy, 22 iconic materials as paradigms for 3D materials spectroscopy, and the VN - defect in hexagonal boron nitride as an example of the spectroscopy of defects in 2D materials pertaining to nanophotonics. Defects can equally be thought of as being "molecular" and "materials" in nature and hence bridge the relms of molecular and materials spectroscopies. It is concluded that the density functional HSE06, currently considered as the standard for accurate calculations of materials spectroscopy, should be replaced, in most instances, by the computationally similar but asymptotically corrected CAM-B3LYP functional, with some specific functionals for materials-use only providing further improvements.

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.

Analysis of the kinetic energy functional in the generalized gradient approximation

A new density functional for the total kinetic energy in the generalized gradient approximation is developed through an enhancement factor that leads to the correct behavior in the limits when the reduced density gradient tends to 0 and to infinity and by making use of the conjoint conjecture for the interpolation between these two limits, through the incorporation, in the intermediate region of constraints that are associated with the exchange energy functional. The resulting functional leads to a reasonable description of the kinetic energies of atoms and molecules when it is used in combination with Hartree-Fock densities. Additionally, in order to improve the behavior of the kinetic energy density, a new enhancement factor for the Pauli kinetic energy is proposed by incorporating the correct behavior into the limits when the reduced density gradient tends to 0 and to infinity, together with the positivity condition, and imposing through the interpolation function that the sum of its integral over the whole space and the Weiszacker energy must be equal to the value obtained with the enhancement factor developed for the total kinetic energy.

Laplacian free and asymptotic corrected semilocal exchange potential applied to the band gap of solids

It is well known that the modified semilocal exchange potentials explicitly designed for the study of solid-state band gaps are very successful in describing these properties. These exchange potentials are in principle designed either from a spherically averaged exchange hole or by satisfying the exact asymptotic conditions. In this present attempt, we use the recently developed novel technique of density matrix expansion to construct the model exchange hole potential. The proposed exchange hole potential is free from the Laplacian of density and generalized through the coordinate transformation. An improvement in the exchange energies of atoms using this potential is shown. The salient feature of the proposed semilocal potential is that it can be used within generalized Kohn-Sham formalism because of its Laplacian free representation. This modified potential is used in the framework of TBMBJ [Phys. Rev. Lett., 2009, 102, 226401] to calculate the band gaps of materials. The comparison and assessment of the newly constructed Laplacian free, asymptotically corrected semilocal potential to address the band gap problem show good agreement with the experimental band gaps and diversify the studies done in the same direction.