Search

To accurately describe the energetics of transition metal systems, density functional approximations (DFAs) must provide a balanced description of s- and d- electrons. One measure of this is the sd transfer error, which has previously been defined as $E(\mathrm{3d}^{n-1}\mathrm{4s}^{1}) - E(\mathrm{3d}^{n-2}\mathrm{4s}^{2})$. Theoretical concerns have been raised about this definition due to its evaluation of excited-state energies using ground-state DFAs. A more serious concern appears to be strong correlation in the ${4s}^2$ configuration. Here we define a ground-state measure of the sd energy imbalance, based on the errors of s- and d-electron second ionization energies of the 3d atoms, that effectively circumvents the aforementioned problems. We find an improved performance as we move from LSDA to PBE to r$^2$SCAN for first-row transition metal atoms. However, we find large (~ 2 eV) ground-state sd energy imbalances when applying a Perdew-Zunger 1981 self-interaction correction. This is attributed to an "energy penalty" associated with the noded 3d orbitals. A local scaling of the self-interaction correction to LSDA results in a balance of s- and d-errors.

Monolayer molybdenum disulfide ($\mathrm{MoS_2}$) under strain has many interesting properties and possible applications in technology. A recent experimental study examined the effect of strain on the bandgap of monolayer $\mathrm{MoS_2}$ on a mildly curved graphite surface, reporting that under biaxial strain with a Poisson's ratio of 0.44, the bandgap decreases at a rate of 400 meV/\% strain. In this work, we performed density functional theory (DFT) calculations for a free-standing $\mathrm{MoS_2}$ monolayer, using the generalized gradient approximation (GGA) PBE, the hybrid functional HSE06, and many-body perturbation theory with the GW approximation using PBE wavefunctions (G0W0@PBE). For the unstrained monolayer, we found a standard level of agreement for the bandgap between theory and experiment. For biaxial strain at the experimental Poisson's ratio, we found that the bandgap decreases at rates of 63 meV/\% strain (PBE), 73 meV/\% strain (HSE06), and 43 meV/\% strain (G0W0@PBE), which are significantly smaller than the experimental rate. We also found that PBE predicts a similarly smaller rate (90 meV/\% strain) for a different Poisson's ratio of 0.25. Spin-orbit correction (SOC) has little effect on the gap or its strain dependence. The strong disagreement between theory and experiment may reflect an unexpectedly strong effect of the substrate on the strain dependence of the gap. Additionally, we observed a transition from a direct to an indirect bandgap under strain, and (under an equal biaxial strain of 10\%) a semiconductor-to-metal transition, consistent with previous theoretical work.

Defects in a material can significantly tune properties and enhance utility. Hybrid functionals like HSE06 are often used to describe solids with such defects. However, geometry optimization (including accounting for effects such as Jahn-Teller distortion) using hybrid functionals is challenging for the large supercells needed for defect study. The proposed r$^2$SCAN+rVV10+U+$\mathrm{U_d}$ method, which is computationally much cheaper and faster than hybrid functionals, can successfully describe defects in materials with the proper choice of U (for the d orbitals of the host atom) and $\mathrm{U_d}$ (for those of the defect atom), as shown here for small polaron defects in layered transition-metal oxides. For a range of U and $\mathrm{U_d}$ around literature values (from solid-state reaction energies) for a given transition-metal ion and its oxidation state, we find that this approach predicts localized polaronic states in band gaps, as hybrid functionals do. The layered materials birnessite ($\mathrm{K_nMnO_2}, n=0.03 $) and $\mathrm{K_nNiO_2},n=0.03$, with one K atom intercalated between layers in a supercell, are found to have one localized occupied $\mathrm{e_g}$ polaronic state on the transition metal ion reduced by the insertion of the K atom, when the geometry is calculated as above using published U values. The expected Jahn-Teller distortion is not observed when U=$\mathrm{U_d}$=0. Layered cobalt oxide with additional potassium ions intercalated ($\mathrm{K_nCoO_2},n=1.03$) is different, due to a dramatic difference in electronic configuration of the defected Co(II) ion: A single extra K atom in the supercell leads to four localized electrons in the band gap, using standard U values, and even for U=$\mathrm{U_d}$=0.

The enigmatic mechanism underlying unconventional high-temperature superconductivity, especially the role of lattice dynamics, has remained a subject of debate. Theoretical insights have long been hindered due to the lack of an accurate first-principles description of the lattice dynamics of cuprates. Recently, using the r2SCAN meta-GGA functional, we were able to achieve accurate phonon spectra of an insulating cuprate YBa$_2$Cu$_3$O$_6$, and discover significant magnetoelastic coupling in experimentally interesting Cu-O bond stretching optical modes [Ning et al., Phys. Rev. B 107, 045126 (2023)]. We extend this work by comparing PBE and r2SCAN performances with corrections from the on-site Hubbard U and the D4 van der Waals (vdW) methods, aiming at further understanding on both the materials science side and the density functional side. We demonstrate the importance of vdW and self-interaction corrections for accurate first-principles YBa2 Cu3 O6 lattice dynamics. Since r2SCAN by itself partially accounts for these effects, the good performance of r2SCAN is now more fully explained. In addition, the performances of the Tao-Mo series of meta-GGAs, which are constructed in a different way from SCAN/r2SCAN, are also compared and discussed., Comment: arXiv admin note: text overlap with arXiv:2210.06569

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase characterized by V-V dimerized structures, to a metallic rutile (R) phase above 340 Kelvin. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to better treat 3d electrons, can accurately predict the V-V dimer length. The spin-restricted method tends to overestimate the strength of the V-V bonds, resulting in a small V-V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each bond-calculation method underscores one of the two contentious mechanisms, i.e., Peierls or Mott, involved in the metal-insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations., Comment: 14 pages, 6 figures

Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally-efficient semi-local approximations for the exchange-correlation energy. The accuracy of a semi-local density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree-Fock electron densities, as known for about 30 years. The conventional explanation is that Hartree-Fock theory yields the more accurate density. This article presents a benchmark Kohn-Sham inversion of accurate coupled-cluster densities for the reaction H$_2$ + F $\rightarrow$ HHF $\rightarrow$ H + HF, and finds a strong, understandable cancellation between positive (excessively over-corrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree-Fock density functional theory. This confirms earlier conclusions [Kaplan et al., J. Chem. Theory Comput. 19, 532--543 (2023)] based on 76 barrier heights and three less reliable, but less expensive, fully-nonlocal density-functional proxies for the exact density., Comment: Revisions in response to peer comments. Addition of Supplemental Material. Addition of BCCD analysis, further KS inversion analysis

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C2 singlet ground state is known to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B2, C2, O2, and F2 in the local spin density approximation (LSDA), the Perdew-Burke-Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C2. This work shows that spin symmetry breaking in singlet C2, the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that under-binding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry-breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking, and the insights into strong correlation that symmetry-breaking can bring., Comment: 10 pages, 3 figures, 1 Table

A central aim of materials discovery is an accurate and numerically reliable description of thermodynamic properties, such as the enthalpies of formation and decomposition. The r$^2$SCAN revision of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) balances numerical stability with high general accuracy. To assess the r$^2$SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental bandgaps of more than 1,000 solids using r$^2$SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r$^2$SCAN+rVV10. We show that r$^2$SCAN achieves accuracy comparable to SCAN and often improves upon SCAN's already excellent accuracy. Whereas SCAN+rVV10 is often observed to worsen the formation enthalpies of SCAN, and makes no substantial correction to SCAN's cell volume predictions, r$^2$SCAN+rVV10 predicts marginally less-accurate formation enthalpies than r$^2$SCAN, and slightly more-accurate cell volumes than r$^2$SCAN. The average absolute errors in predicted formation enthalpies are found to decrease by a factor of 1.5 to 2.5 from the GGA level to the meta-GGA level. Smaller decreases in error are observed for decomposition enthalpies. For formation enthalpies r$^2$SCAN improves over SCAN for intermetallic systems. For a few classes of systems -- transition metals, intermetallics, weakly-bound solids, and enthalpies of decomposition into compounds -- GGAs are comparable to meta-GGAs. In total, r$^2$SCAN and r$^2$SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.

Delocalization errors, such as charge-transfer and some self-interaction errors, plague computationally-efficient and otherwise-accurate density functional approximations (DFAs). Evaluating a semi-local DFA non-self-consistently on the Hartree-Fock (HF) density is often recommended as a computationally cheap remedy for delocalization errors. For sophisticated meta-GGAs like SCAN, this approach can achieve remarkable accuracy. When this HF-DFT (or DFA@HF) significantly improves over the DFA, it is often presumed that the HF density is more accurate than the self-consistent DFA density. By applying the metrics of density-corrected density functional theory (DFT), we show that HF-DFT works for barrier heights by making a localizing charge transfer error or density over-correction, thereby producing a somewhat-reliable cancellation of density- and functional-driven errors for the energy. A quantitative analysis of the charge transfer errors in a few transition states confirms this trend. We do not have the exact functional and exact densities that are needed to evaluate the exact density- and functional-driven errors for the large BH76 database of barrier heights. Instead, we have identified and used three non-local proxy functionals (the SCAN 50% global hybrid, the range-separated hybrid LC-$\omega$PBE, and SCAN-FLOSIC) and their self-consistent densities. These functionals yield reasonably accurate self-consistent barrier heights, and their self-consistent total energies are nearly piecewise linear in fractional electron number - two important points of similarity to the exact functional. We argue that density-driven errors of the energy in a self-consistent density functional calculation are second-order in the density error, and that large density-driven errors arise primarily from incorrect electron transfers over length scales larger than the diameter of an atom., Comment: v1: Submitted in preparation for the annual FLOSIC all-hands meeting. v2: Significant revisions prior to submission at JCTC. v3: Significant revisions (completely new Supp Info) in response to private discussions with colleagues and peer review (NB: arXiv limit on number of characters in abstract, so it's truncated)

Ground-state Kohn-Sham density functional theory provides, in principle, the exact ground-state energy and electronic spin-densities of real interacting electrons in a static external potential. In practice, the exact density functional for the exchange-correlation (xc) energy must be approximated in a computationally efficient way. About twenty mathematical properties of the exact xc functional are known. In this work, we review and discuss these known constraints on the xc energy and hole. By analyzing a sequence of increasingly sophisticated density functional approximations (DFAs), we argue that: (1) the satisfaction of more exact constraints and appropriate norms makes a functional more predictive over the immense space of many-electron systems; (2) fitting to bonded systems yields an interpolative DFA that may not extrapolate well to systems unlike those in the fitting set. We discuss how the class of well-described systems has grown along with constraint satisfaction, and the possibilities for future functional development., Comment: When citing this paper, please use the following: Kaplan AD, Levy M, Perdew JP. 2023. Predictive power of the exact constraints and approximate norms in density functional theory. Annu. Rev. Phys. Chem. 74. Submitted. DOI: 10.1146/annurev-physchem-062422-013259

Lieb and Oxford (1981) derived rigorous lower bounds, in the form of local functionals of the electron density, on the indirect part of the Coulomb repulsion energy. The greatest lower bound for a given electron number N depends monotonically upon N, and the N-> infinity limit is a bound for all N. These bounds have been shown to apply to the exact density functionals for the exchange- and exchange-correlation energies that must be approximated for an accurate and computationally efficient description of atoms, molecules, and solids. A tight bound on the exact exchange energy has been derived therefrom for two-electron ground states, and is conjectured to apply to all spin-unpolarized electronic ground states. Some of these and other exact constraints have been used to construct two generations of non-empirical density functionals beyond the local density approximation: the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA., Comment: 16 pages

In a data-driven paradigm, machine learning (ML) is the central component for developing accurate and universal exchange-correlation (XC) functionals in density functional theory (DFT). It is well known that XC functionals must satisfy several exact conditions and physical constraints, such as density scaling, spin scaling, and derivative discontinuity. In this work, we demonstrate that contrastive learning is a computationally efficient and flexible method to incorporate a physical constraint in ML-based density functional design. We propose a schematic approach to incorporate the uniform density scaling property of electron density for exchange energies by adopting contrastive representation learning during the pretraining task. The pretrained hidden representation is transferred to the downstream task to predict the exchange energies calculated by DFT. The electron density encoder transferred from the pretraining task based on contrastive learning predicts exchange energies that satisfy the scaling property, while the model trained without using contrastive learning gives poor predictions for the scaling-transformed electron density systems. Furthermore, the model with pretrained encoder gives a satisfactory performance with only small fractions of the whole augmented dataset labeled, comparable to the model trained from scratch using the whole dataset. The results demonstrate that incorporating exact constraints through contrastive learning can enhance the understanding of density-energy mapping using neural network (NN) models with less data labeling, which will be beneficial to generalizing the application of NN-based XC functionals in a wide range of scenarios that are not always available experimentally but theoretically justified. This work represents a viable pathway toward the machine learning design of a universal density functional via representation learning.

SCAN+rVV10 has been demonstrated to be a versatile van der Waals (vdW) density functional that delivers good predictions of both energetic and structural properties for many types of bonding. Recently, the r$^{2}$SCAN functional has been devised as a revised form of SCAN with improved numerical stability. In this work, we refit the rVV10 functional to optimize the r$^{2}$SCAN+rVV10 vdW density functional, and test its performance for molecular interactions and layered materials. Our molecular tests demonstrate that r$^{2}$SCAN+rVV10 outperforms its predecessor SCAN+rVV10 in both efficiency (numerical stability) and accuracy. This good performance is also found in lattice constant predictions. In comparison with benchmark results from higher-level theories or experiments, r$^{2}$SCAN+rVV10 yields excellent interlayer binding energies and phonon dispersions for layered materials.

We derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of $sp$ and $sd$ metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r$^2$SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r$^2$SCAN. We compare OFR2 to the r$^2$SCAN-L LL-MGGA [D. Mejia-Rodriguez and S.B. Trickey, Phys. Rev. B 102, 121109 (2020)] and show that OFR2 tends to outperform r$^2$SCAN-L for the equilibrium properties of solids, but r$^2$SCAN-L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and non-metallic condensed matter, we continue to recommend SCAN and r$^2$SCAN. Numerical performance is discussed in detail, and our work provides an outlook to machine learning., Comment: Significant revisions in response to referee comments. Under review at Physical Review Materials

The SCAN meta-GGA exchange-correlation functional [Phys. Rev. Lett. 115, 036402 (2015)] is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single orbital electron densities accurately, and another describes slowly-varying densities accurately. To conserve constraints known for the exact exchange-correlation functional, the derivatives of this interpolation vanish in the slowly-varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional's efficiency. We have recently reported a modification to the SCAN functional, termed r$^2$SCAN [J. Phys. Chem. Lett. 11, 8208 (2020)] that introduces two regularizations into SCAN which improve its numerical performance at the expense of not recovering the fourth order term of the slowly-varying density gradient expansion for exchange. Here we show the derivation of a progression of functionals (rSCAN, r++SCAN, r$^2$SCAN, and r$^4$SCAN) with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r$^2$SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r$^4$SCAN does.

The atomization energies of molecules from first-principles density functional approximations improve from the local spin-density approximation (LSDA) to the Perdew-Burke-Ernzerhof (PBE)) generalized gradient approximation (GGA) to the strongly constrained and appropriately normed (SCAN) meta-GGA, and their sensitivities to non-spherical components of the density increase in the same order. Thus, these functional advances increase density sensitivity and imitate the exact constrained search over correlated wavefunctions better than that over ensembles. The diatomic molecules studied here, singlet C2 and F2 plus triplet B2 and O2, have cylindrically symmetric densities. Because the densities of the corresponding atoms are non-spherical, the approximate Kohn-Sham potentials for the atoms have a lower symmetry than that of the external (nuclear) potential, so that the non-interacting wavefunctions are not eigenstates of the square of total orbital angular momentum, breaking a symmetry that yields a feature of the exact ground-state density. That spatial symmetry can be preserved by a non-self-consistent approach in which a self-consistent equilibrium-ensemble calculation is followed by integer re-occupation of the Kohn-Sham orbitals, as the first of several steps. The symmetry-preserving approach is different from symmetry restoration based upon projection. First-step space- (and space-spin-) symmetry preservation in atoms is shown to have a small effect on the atomization energies of molecules, quantifying earlier observations by Fertig and Kohn. Thus, the standard Kohn-Sham way of calculating atomization energies, with self-consistent symmetry breaking to minimize the energy, is justified, at least for the common cases where the molecules cannot break symmetry., Comment: 22 Pages, 8 Tables, 1 Figure

We propose a spatially and temporally nonlocal exchange-correlation (xc) kernel for the spin-unpolarized fluid phase of ground-state jellium, for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact xc kernel, to accurately describe the correlation energies of bulk jellium, and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel [A. Ruzsinszky, et al., Phys. Rev. B 101, 245135 (2020)] are maintained in the new revised MCP07 (rMCP07) kernel, while others are added. The revision $f_\mathrm{xc}^\mathrm{rMCP07}(q,\omega)$ differs from MCP07 only for non-zero frequencies $\omega$. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wavevector $q$ tends to zero, the kernel has a $-4\pi \alpha(\omega)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $\alpha(\omega)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.}, Comment: (1st revision) Changes to the TC21 model, new App. F discussing numerical methods. (2nd revision) Revisions in response to peer comments. Refocused discussion around building a kernel for jellium, rather than for metals/materials. (3rd revision) Name change of kernel, TC21 --> rMCP07; added details of QMC calculations; added CP07 and Richardson-Ashcroft correlation energies

Within the framework of Kohn-Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob's ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric-isothermal ensemble. In order to generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water towards the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional.

The Perdew-Zunger self-interaction correction(PZ-SIC) improves the performance of density functional approximations(DFAs) for the properties that involve significant self-interaction error(SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by LSIC, local scaling of the PZ-SIC to the local spin density approximation(LSDA). Here we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of atomic number Z in the asymptotic expansion of the exchange-correlation(xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator z{\sigma}, which distinguishes one-electron regions from many-electron regions. LSIC+ applied to LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof(PBE) generalized gradient approximation(GGA), and almost as well as the strongly constrained and appropriately normed(SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, both fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled-down SIC to functionals beyond LSDA requires a gauge transformation of the functional's energy density. The resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds., Comment: 12 pages, 8 figures

We study the importance of self-interaction errors in density functional approximations for various water-ion clusters. We have employed the Fermi-L\"owdin orbital self-interaction correction (FLOSIC) method in conjunction with LSDA, PBE, and SCAN to describe binding energies of hydrogen-bonded water-ion clusters, \textit{i.e.}, water-hydronium, water-hydroxide, water-halide, as well as non-hydrogen-bonded water-alkali clusters. In the hydrogen-bonded water-ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules, because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water-ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher-level calculations. The non-hydrogen-bonded water-alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies ($\leq$ 50 meV/\ce{H2O}) for hydrogen-bonded water-ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement., Comment: 13 pages, 9 figures

We combine a regularized variant of the strongly constrained and appropriately normed semilocal density functional [J. Sun, A. Ruzsinszky, and J. P. Perdew, Phys. Rev. Lett. 115, 036402 (2015)] with the latest generation semi-classical London dispersion correction. The resulting density functional approximation r2SCAN-D4 has the speed of generalized gradient approximations while approaching the accuracy of hybrid functionals for general chemical applications. We demonstrate its numerical robustness in real-life settings and benchmark molecular geometries, general main group and organo-metallic thermochemistry, as well as non-covalent interactions in supramolecular complexes and molecular crystals. Main group and transition metal bond lengths have errors of just 0.8%, which is competitive with hybrid functionals for main group molecules and outperforms them for transition metal complexes. The weighted mean absolute deviation (WTMAD2) on the large GMTKN55 database of chemical properties is exceptionally small at 7.5 kcal/mol. This also holds for metal organic reactions with an MAD of 3.3 kcal/mol. The versatile applicability to organic and metal-organic systems transfers to condensed systems, where lattice energies of molecular crystals are within chemical accuracy (errors <1 kcal/mol).

The recently proposed rSCAN functional [J. Chem. Phys. 150, 161101 (2019)] is a regularized form of the SCAN functional [Phys. Rev. Lett. 115, 036402 (2015)] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN., Comment: To appear in the Journal of Physical Chemistry Letters. Updates to the main text for clarity regarding the interpolation functions. TPSS data added throughout for comparison. VASP real-space convergence test, determination of regularization parameter eta, additional figures and data added to Supplementary Material

Strong correlations within a symmetry-unbroken ground-state wavefunction can show up in approximate density functional theory as symmetry-broken spin-densities or total densities, which are sometimes observable. They can arise from soft modes of fluctuations (sometimes collective excitations) such as spin-density or charge-density waves at non-zero wavevector. In this sense, an approximate density functional for exchange and correlation that breaks symmetry can be more revealing (albeit less accurate) than an exact functional that does not. The examples discussed here include the stretched H$_2$ molecule, antiferromagnetic solids, and the static charge-density wave/Wigner crystal phase of a low-density jellium. It is shown that (and in what sense) the static charge density wave is a soft plasmon.

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-$Z$ expansions for the corresponding energies of neutral atoms with atomic number $Z$ and electron number $N=Z$, which are correct to leading order ($-0.221 Z^{5/3}$ and $-0.021 Z \ln Z$ respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to $E_x(N,Z) \approx -0.354 N^{2/3} Z$ (as known before only for $Z \gg N \gg 1$) and $E_c \approx -0.02 N \ln N$. These asymptotic estimates are most correct for atomic ions with large $N$ and $Z \gg N$, but we find that they are qualitatively and semi-quantitatively correct even for small $N$ and for $N \approx Z$. The large-$N$ asymptotic behavior of the energy is pre-figured in small-$N$ atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to $Z$ in the $Z\to \infty$ limit for any fixed $N$., Comment: This work has been accepted for publication at the Journal of Chemical Physics. Revisions: new Appendix A (former Appendix A is now Appendix B) discussing exact Kohn-Sham perturbation series for Ec. Added material discussing the Becke 1988 functional. More discussion of non-empirical functionals' recovery of the asymptotic series, and their accuracy in predicting atomic/molecular energies

Classical turning surfaces of Kohn-Sham potentials, separating classically-allowed regions (CARs) from classically-forbidden regions (CFRs), provide a useful and rigorous approach to understanding many chemical properties of molecules. Here we calculate such surfaces for several paradigmatic solids. Our study of perfect crystals at equilibrium geometries suggests that CFRs are absent in metals, rare in covalent semiconductors, but common in ionic and molecular crystals. A CFR can appear at a monovacancy in a metal. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators. This classical picture is confirmed in the limits of extreme uniform compression of the internuclear distances, where all materials become metals without CFRs, and extreme expansion, where all materials become insulators with disconnected and widely-separated CARs around the atoms., Comment: Added supplemental information (63 pages), was missing in original submission. Minor typo corrections in Tables I and III for Eps_HO - vs(r) column (and CFR % volume for Pt monovacancy) for PBE data only

The Perdew-Zunger(PZ) self-interaction correction (SIC) was designed to correct the one-electron limit of any approximate density functional for the exchange-correlation (xc) energy, while yielding no correction to the exact functional. Unfortunately, it spoils the slowly-varying-in-space limits of the uncorrected approximate functionals, where those functionals are right by construction. The right limits can be restored by locally scaling down the energy density of the PZ SIC in many-electron regions, but then a spurious correction to the exact functional would be found unless the self-Hartree and exact self-xc terms of the PZ SIC energy density were expressed in the same gauge. Only the local density approximation satisfies the same-gauge condition for the energy density, which explains why the recent local-scaling SIC (LSIC) is found here to work excellently for atoms and molecules only with this basic approximation, and not with the more advanced generalized gradient approximations (GGAs) and meta-GGAs, which lose the Hartree gauge via simplifying integrations by parts. The transformation of energy density that achieves the Hartree gauge for the exact xc functional can also be applied to approximate functionals. Doing so leads to a simple scaled-down self-interaction (sdSIC) correction that is typically much more accurate than PZ SIC in tests for many molecular properties (including equilibrium bond lengths). The present work shows unambiguously that the largest errors of PZ SIC applied to standard functionals at three levels of approximation can be removed by restoring their correct slowly-varying-density limits. It also confirms the relevance of these limits to atoms and molecules., Comment: 21 pages, 2 figures

According to time-dependent density functional theory, the exact exchange-correlation kernel f$_{xc}$(n, q, $\omega$) determines not only the ground-state energy but also the excited-state energies/lifetimes and time-dependent linear density response of an electron gas of uniform density n $=$ 3/(4$\pi$r$^3_s$). Here we propose a parametrization of this function based upon the satisfaction of exact constraints. For the static ($\omega$ = 0) limit, we modify the model of Constantin and Pitarke at small wavevector q to recover the known second-order gradient expansion, plus other changes. For all frequencies $\omega$ at q $=$ 0, we use the model of Gross, Kohn, and Iwamoto. A Cauchy integral extends this model to complex $\omega$ and implies the standard Kramers-Kronig relations. A scaling relation permits closed forms for not only the imaginary but also the real part of f$_{xc}$ for real $\omega$. We then combine these ingredients by damping out the $\omega$ dependence at large q in the same way that the q dependence is damped. Away from q $=$ 0 and $\omega$ $=$ 0, the correlation contribution to the kernel becomes dominant over exchange, even at r$_s$ $=$ 4, the valence electron density of metallic sodium. The resulting correlation energy from integration over imaginary $\omega$ is essentially exact. The plasmon pole of the density response function is found by analytic continuation of f$_{xc}$ to $\omega$ just below the real axis, and the resulting plasmon lifetime first decreases from infinity and then increases as q grows from 0 toward the electron-hole continuum. A static charge-density wave is found for r$_s$ $>$ 69, and shown to be associated with softening of the plasmon mode. The exchange-only version of our static kernel confirms Overhauser's 1968 prediction that correlation enhances the charge-density wave., Comment: 20 pages including 11 figures

Self-interaction (SI) error, which results when exchange-correlation contributions to the total energy are approximated, limits the reliability of many density functional approximations. The Perdew-Zunger SI correction (PZSIC), when applied in conjunction with the local spin density approximation (LSDA), improves the description of many properties, but overall, this improvement is limited. Here we propose a modification to PZSIC that uses an iso-orbital indicator to identify regions where local SI corrections should be applied. Using this local-scaling SIC (LSIC) approach with LSDA, we analyze predictions for a wide range of properties including, for atoms, total energies, ionization potentials, and electron affinities, and for molecules, atomization energies, dissociation energy curves, reaction energies, and reaction barrier heights. LSIC preserves the results of PZSIC-LSDA for properties where it is successful and provides dramatic improvements for many of the other properties studied. Atomization energies calculated using LSIC are better than those of the Perdew, Burke, and Ernzerhof (PBE) generalized gradient approximation (GGA) and close to those obtained with the Strongly Constrained and Appropriately Normed (SCAN) meta-GGA. LSIC also restores the uniform gas limit for the exchange energy that is lost in PZSIC-LSDA. Further performance improvements may be obtained by an appropriate combination or modification of the local scaling factor and the particular density functional approximation., Comment: 12 pages, 6 figures, to be published in Journal of Chemical Physics

The random phase approximation (RPA) is exact for the exchange energy of a many-electron ground state, but RPA makes the correlation energy too negative by about 0.5 eV/electron. That large short-range error, which tends to cancel out of iso-electronic energy differences, is largely corrected by an exchange-correlation kernel, or (as in RPA+) by an additive local or semilocal correction. RPA+ is by construction exact for the homogeneous electron gas, and it is also accurate for the jellium surface. RPA+ often gives realistic total energies for atoms or solids in which spin-polarization corrections are absent or small. RPA and RPA+ also yield realistic singlet binding energy curves for H2 and N2, and thus RPA+ yields correct total energies even for spin-unpolarized atoms with fractional spins and strong correlation, as in stretched H2 or N2. However, RPA and RPA+ can be very wrong for spin-polarized one-electron systems (especially for stretched H2+), and also for the spin-polarization energies of atoms. The spin-polarization energy is often a small part of the total energy of an atom, but important for ionization energies, electron affinities, and the atomization energies of molecules. Here we propose a computationally efficient generalized RPA+ (gRPA+) that changes RPA+ only for spin-polarized systems by making gRPA+ exact for all one-electron densities, in the same simple semilocal way that the correlation energy densities of many meta-generalized gradent approximations are made self-correlation free. By construction, gRPA+ does not degrade the exact RPA+ description of jellium. gRPA+ is found to greatly improve upon RPA and RPA+ for the ionization energies and electron affinities of light atoms. Many versions of RPA with an approximate exchange-correlation kernel fail to be exact for all one-electron densities, and they can also be self-interaction corrected in this way.

Semi-local approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely-related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semi-local approximation makes that approximation exact for all one-electron ground- or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semi-local approximations satisfy, and suggesting the need for a generalized SIC. These conclusions are supported by calculations for one-electron densities, and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGA's) and meta-GGA's, it reduces and often worsens the atomization energies of molecules. Thus PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here in particular to the SCAN meta-GGA, for which the correlation part is already self-interaction-free. That property makes SCAN a natural first candidate for a generalized SIC.

Semi-local density functionals for the exchange-correlation energy of a many-electron system cannot be exact for all one-electron densities. In 1981, Perdew and Zunger (PZ) subtracted the fully-nonlocal self-interaction error orbital-by-orbital, making the corrected functional exact for all collections of separated one-electron densities, and making no correction to the exact functional. Although the PZ self-interaction correction (SIC) eliminates many errors of semi-local functionals, it is often worse for equilibrium properties of sp-bonded molecules and solids. Non-empirical semi-local functionals are usually designed to be exact for electron gases of uniform density, and thus also make 0% error for neutral atoms in the limit of large atomic number Z, but PZ SIC is not so designed. For localized SIC orbitals, we show analytically that the LSDA-SIC correlation energy per electron of the uniform gas in the high-density limit makes an error of -50% in the spin-unpolarized case, and -100% in the fully-spin-polarized case. Then we extrapolate from the Ne, Ar, Kr, and Xe atoms to estimate the relative errors of the PZ SIC exchange-correlation energies (with localized SIC orbitals) in the limit of large atomic number: about +5.5% for the local spin density approximation (LSDA-SIC), and about -3.5% for nonempirical generalized gradient (PBE-SIC) and meta-generalized gradient (SCAN-SIC) approximations. The SIC errors are considerably larger than those that have been estimated for LSDA-SIC by approximating the localized SIC orbitals for the uniform gas, and may explain the errors of PZ SIC for equilibrium properties, opening the door to a generalized SIC that is more widely accurate.

Adsorption of the molecule CO on metallic surfaces is an important unsolved problem in Kohn-Sham density functional theory (KS-DFT). We present a detailed study of carbon monoxide adsorption on fcc (111) surfaces of 3d, 4d and 5d metals using nonempirical semilocal density functionals for the exchange-correlation energy: the local-density approximation (LDA), two generalized gradient approximations or GGAs (PBE and PBEsol), and a meta-GGA (SCAN). The typical error pattern (as found earlier for free molecules and for free transition metal surfaces), in which results improve from LDA to PBE or PBEsol to SCAN, due to the satisfaction of more exact constraints, is not found here. Instead, for CO adsorption on transition metal surfaces, we find that, while SCAN overbinds much less than LDA, it overbinds slightly more than PBE. Moreover, the tested functionals often predict the wrong adsorption site, as first pointed out for LDA and GGA in the ``CO/Pt (111) puzzle". This abnormal pattern leads us to suspect that the errors of PBE and SCAN for this problem are density-driven self-interaction errors associated with incorrect charge transfer between molecule and metal surface. We point out that, by the variational principle, overbinding by an approximate functional would be reduced if that functional were applied not to its selfconsistent density for the adsorbed system but to an exact or more correct density for that system. Finally, we show for CO on Pt(111) that the site preference is corrected and the adsorption energy is improved for the PBE functional by using not the selfconsistent PBE density but a PBE+U density. The resulting correction to the PBE total energy is much larger for the adsorbed system than for its desorbed components, showing that the error is in the density of the adsorbed system. This seems to solve the ``CO/Pt (111) puzzle", in principle if not fully in practice.

Nanostructures can be bound together at equilibrium by the van der Waals (vdW) effect, a small but ubiquitous many-body attraction that presents challenges to density functional theory. How does the binding energy depend upon the size or number of atoms in one of a pair of identical nanostructures? To answer this question, we treat each nanostructure properly as a whole object, not as a collection of atoms. Our calculations start from an accurate static dipole polarizability for each considered nanostructure, and an accurate equilibrium center-to-center distance for the pair (the latter from experiment, or from the vdW-DF-cx functional). We consider the competition in each term $-C_{2k}/d^{2k}$ ($k=3, 4, 5$) of the long-range vdW series for the interaction energy, between the size dependence of the vdW coefficient $C_{2k}$ and that of the $2k$-th power of the center-to-center distance $d$. The damping of these vdW terms can be negligible, but in any case it does not affect the size dependence for a given term in the absence of non-vdW binding. To our surprise, the vdW energy can be size-independent for quasi-spherical nanoclusters bound to one another by vdW interaction, even with strong nonadditivity of the vdW coefficient, as demonstrated for fullerenes. We also show that, for low-dimensional systems, the vdW interaction yields the strongest size-dependence, in stark contrast to that of fullerenes. We illustrate this with parallel planar polycyclic aromatic hydrocarbons. Other cases are between, as shown by sodium clusters., Comment: 7 pages, 1 figure

The equilibrium binding energy is an important factor in the design of materials and devices. However, it presents great computational challenges for materials built up from nanostructures. Here we investigate the binding-energy scaling law from first-principles calculations. We show that the equilibrium binding energy per atom between identical nanostructures can scale up or down with nanostructure size. From the energy scaling law, we predict finite large-size limits of binding energy per atom. We find that there are two competing factors in the determination of the binding energy: Nonadditivities of van der Waals coefficients and center-to-center distance between nanostructures. To uncode the detail, the nonadditivity of the static multipole polarizability is investigated. We find that the higher-order multipole polarizability displays ultra-strong intrinsic nonadditivity, no matter if the dipole polarizability is additive or not., Comment: 13 pages, 4 figures, 7 tables

Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice I{\it h} at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner.

The Perdew-Zunger self-interaction correction cures many common problems associated with semilocal density functionals, but suffers from a size-extensivity problem when Kohn-Sham orbitals are used in the correction. Fermi-L\"{o}wdin-orbital self-interaction correction (FLOSIC) solves the size-extensivity problem, allowing its use in periodic systems and resulting in better accuracy in finite systems. Although the previously published FLOSIC algorithm [J. Chem. Phys. 140, 121103 (2014)] appears to work well in many cases, it is not fully self-consistent. This would be particularly problematic for systems where the occupied manifold is strongly changed by the correction. In this paper we demonstrate a new algorithm for FLOSIC to achieve full self-consistency with only marginal increase of computational cost. The resulting total energies are found to be lower than previously reported non-self-consistent results., Comment: 8 pages, 5 figures

We have computed the surface energies, work functions, and interlayer surface relaxations of clean (111), (110), and (100) surfaces of Al, Cu, Ru, Rh, Pd, Ag, Pt, and Au. Many of these metallic surfaces have technological or catalytic applications. We compare experimental reference values to those of the local density approximation (LDA), the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), the PBEsol (PBE for solids) GGA, the SCAN meta-GGA, and SCAN+rVV10 (SCAN with a long-range van der Waals or vdW correction). The closest agreement with uncertain experimental values is achieved by the simplest density functional (LDA) and by the most sophisticated general-purpose one (SCAN+rVV10). The long-range vdW interaction increases the surface energies by about 10%, and the work functions by about 1%. LDA works for metal surfaces through a stronger-than-usual error cancellation. PBE yields the most-underestimated and presumably least accurate surface energies and work functions. Surface energies within the random phase approximation (RPA) are also reported. Interlayer relaxations from different functionals are in reasonable agreement with one another, and usually with experiment.

Originating from a broken spatial inversion symmetry, ferroelectricity is a functionality of materials with an electric dipole that can be switched by external electric fields. Spontaneous polarization is a crucial ferroelectric property, and its amplitude is determined by the strength of polar structural distortions. Density functional theory (DFT) is one of the most widely used theoretical methods to study ferroelectric properties, yet it is limited by the levels of approximations in electron exchange-correlation. On the one hand, the local density approximation (LDA) is considered to be more accurate for the conventional perovskite ferroelectrics such as BaTiO3 and PbTiO3 than the generalized gradient approximation (GGA), which suffers from the so-called super-tetragonality error. On the other hand, GGA is more suitable for hydrogen-bonded ferroelectrics than LDA, which largely overestimates the strength of hydrogen bonding in general. We show here that the recently developed general-purpose strongly constrained and appropriately normed (SCAN) meta-GGA functional significantly improves over the traditional LDA/GGA for structural, electric, and energetic properties of diversely-bonded ferroelectric materials with a comparable computational effort, and thus enhances largely the predictive power of DFT in studies of ferroelectric materials. We also address the observed system-dependent performances of LDA and GGA for ferroelectrics from a chemical bonding point of view., Comment: 23 pages, 7 figures, 5 tables, 137 references

The structural and energetic properties of layered materials propose a challenge to density functional theory with common semilocal approximations to the exchange-correlation. By combining the most-widely used semilocal generalized gradient approximation (GGA), Perdew--Burke--Ernzerhof (PBE), with the revised Vydrov--Van Voorhis non-local correlation functional (rVV10), both excellent structural and energetic properties of 28 layered materials were recovered with a judicious parameter selection. We term the resulting functional as PBE+rVV10L with "L" denoting for layered materials. Such combination is not new, and involves only refitting a single global parameter, however the resulting excellency suggests such corrected PBE for many aspects of theoretical studies on layered materials. For comparison, we also present the results for PBE+rVV10 where the parameter is determined by the 22 interaction energies between molecules.

We present a new paradigm for the design of exchange-correlation functionals in density-functional theory. Electron pairs are correlated explicitly by means of the recently developed second order Bethe-Goldstone equation (BGE2) approach. Here we propose a screened BGE2 (sBGE2) variant that efficiently regulates the coupling of a given electron pair. sBGE2 correctly dissociates H$_2$ and H$_2^+$, a problem that has been regarded as a great challenge in density-functional theory for a long time. The sBGE2 functional is then taken as a building block for an orbital-dependent functional, termed ZRPS, which is a natural extension of the PBE0 hybrid functional. While worsening the good performance of sBGE2 in H$_2$ and H$_2^{+}$, ZRPS yields a remarkable and consistent improvement over other density functionals across various chemical environments from weak to strong correlation., Comment: 5 pages, 2 figures

The fundamental energy gap of a periodic solid distinguishes insulators from metals and characterizes low-energy single-electron excitations. But the gap in the band-structure of the exact multiplicative Kohn-Sham (KS) potential substantially underestimates the fundamental gap, a major limitation of KS density functional theory. Here we give a simple proof of a new theorem: In generalized KS theory (GKS), the band gap of an extended system equals the fundamental gap for the approximate functional if the GKS potential operator is continuous and the density change is delocalized when an electron or hole is added. Our theorem explains how GKS band gaps from meta-generalized gradient approximations (meta-GGAs) and hybrid functionals can be more realistic than those from GGAs or even from the exact KS potential. The theorem also follows from earlier work. The band edges in the GKS one-electron spectrum are also related to measurable energies. A linear chain of hydrogen molecules, solid aluminum arsenide, and solid argon provide numerical illustrations.

The uniform electron gas and the hydrogen atom play fundamental roles in condensed matter physics and quantum chemistry. The former has an infinite number of electrons uniformly distributed over the neutralizing positively-charged background, and the latter only one electron bound to the proton. The uniform electron gas was used to derive the local spin density approximation (LSDA) to the exchange-correlation functional that undergirds the development of the Kohn-Sham density functional theory. We show here that the ground-state exchange-correlation energies of the hydrogen atom and many other 1- and 2-electron systems are modeled surprisingly well by a different local spin density approximation (LSDA0). Our LSDA0 is constructed to satisfy exact constraints, but agrees surprisingly well with the exact results for a uniform two-electron density in a finite, curved three-dimensional space. We also apply LSDA0 to excited or noded 1-electron densities. Furthermore, we show that the locality of an orbital can be measured by the ratio between the exact exchange energy and its optimal lower bound.

Unlike the local density approximation (LDA) and the generalized gradient approximation (GGA), calculations with meta-generalized gradient approximations (meta-GGA) are usually done according to the generalized Kohn-Sham (gKS) formalism. The exchange-correlation potential of the gKS equation is non-multiplicative, which prevents systematic comparison of meta-GGA bandstructures to those of the LDA and the GGA. We implement the optimized effective potential (OEP) of the meta-GGA for periodic systems, which allows us to carry out meta-GGA calculations in the same KS manner as for the LDA and the GGA. We apply the OEP to several meta-GGAs, including the new SCAN functional [Phys. Rev. Lett. 115, 036402 (2015)]. We find that the KS gaps and KS band structures of meta-GGAs are close to those of GGAs. They are smaller than the more realistic gKS gaps of meta-GGAs, but probably close to the less-realistic gaps in the band structure of the exact KS potential, as can be seen by comparing with the gaps of the EXX+RPA OEP potential. The well-known grid sensitivity of meta-GGAs is much more severe in OEP calculations., Comment: 10 pages, 10 figures