Meta-GGA performance in solids at almost GGA cost
Daniel Mejı́a-Rodrı́guez, S. B. Trickey
Abstract
A recent modification, $\mathrm{r}^{2}\mathrm{SCAN}$, of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA) exchange-correlation functional mostly eliminates numerical instabilities and attendant integration grid sensitivities exhibited by SCAN. Here, we show that the successful deorbitalization of SCAN to SCAN-L (SCAN with density Laplacian dependence) carries over directly to yield $\mathrm{r}^{2}\mathrm{SCAN}$-L. A major benefit is that the high iteration counts that hindered the use of SCAN-L are eliminated in $\mathrm{r}^{2}\mathrm{SCAN}$-L. It therefore is a computationally much faster meta-GGA than its orbital-dependent antecedent. Validation data for molecular heats of formation, bond lengths, and vibration frequencies (G3/99X, T96-R, and T82-F test sets, respectively) and on lattice constants, and cohesive energies (for 55 solids) and bulk moduli (for 40 solids) are provided. In addition, we show that the overmagnetization of bcc Fe, hcp Co, and fcc Ni persists in $\mathrm{r}^{2}\mathrm{SCAN}$ but does not appear in $\mathrm{r}^{2}\mathrm{SCAN}$-L. Distinct from SCAN, both $\mathrm{r}^{2}\mathrm{SCAN}$ and $\mathrm{r}^{2}\mathrm{SCAN}$-L give the correct nonmagnetic ground state for bcc V.