-
The isostructural alpha-gamma phase transition in cerium from the perspective of meta-generalized gradient approximations
Authors:
Ashesh Giri,
Chandra Shahi,
Adrienn Ruzsinszky
Abstract:
Meta-generalized gradient approximations (meta-GGAs) on the third rung of the functional hierarchy are gaining increasing relevance for the electronic structure. Meta-GGAs are constructed from numerous ingredients including the orbital kinetic energy density that make them more flexible than generalized gradient approximations (GGAs) including the heavily used PBE-GGA. Still, most meta-GGAs cope w…
▽ More
Meta-generalized gradient approximations (meta-GGAs) on the third rung of the functional hierarchy are gaining increasing relevance for the electronic structure. Meta-GGAs are constructed from numerous ingredients including the orbital kinetic energy density that make them more flexible than generalized gradient approximations (GGAs) including the heavily used PBE-GGA. Still, most meta-GGAs cope with the expected limitations of a semilocal density functional when band gaps or localization of electrons are needed. On the other hand, meta-GGAs are implicit functionals of the orbitals. This feature resembles hybrid density functionals with exact exchange. Efforts in recent years demonstrate that some meta-GGAs can rise beyond the accuracy of semilocal approximation when band gaps are computed. Cerium is an ideal testbed to challenge some recent meta-GGAs. Cerium shows an isostructural alpha - gamma phase transition with delocalized and localized f electrons in each phase, respectively. Since the phonon entropy term was found negligible in the alpha - gamma phase transition of cerium by accurate experiments, all changes in the transition are driven by electronic correlation. The correlation of f electron systems is hardly captured by semilocal approximations but the recent LAK meta-GGA with ultranonlocality steps out of the framework of conventional semilocal density functionals and delivers spectacular accuracy for the phase transition of cerium. LAK and further meta-GGAs inspired by the success of LAK can open a forefront of meta-GGAs for quantum materials with localized electrons.
△ Less
Submitted 4 June, 2025;
originally announced June 2025.
-
Effect of Strain on the Band Gap of Monolayer MoS$_2$
Authors:
Raj K. Sah,
Hong Tang,
Chandra Shahi,
Adrienn Ruzsinszky,
John P. Perdew
Abstract:
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…
▽ More
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.
△ Less
Submitted 14 November, 2024; v1 submitted 10 June, 2024;
originally announced June 2024.
-
Unconventional Error Cancellation Explains the Success of Hartree-Fock Density Functional Theory for Barrier Heights
Authors:
Bikash Kanungo,
Aaron D. Kaplan,
Chandra Shahi,
Vikram Gavini,
John P. Perdew
Abstract:
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…
▽ More
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.
△ Less
Submitted 2 November, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
A Step in the Direction of Resolving the Paradox of Perdew-Zunger Self-interaction Correction. II. Gauge Consistency of the Energy Density at Three Levels of Approximation
Authors:
Puskar Bhattarai,
Kamal Wagle,
Chandra Shahi,
Yoh Yamamoto,
Selim Romero,
Biswajit Santra,
Rajendra R. Zope,
Juan E. Peralta,
Koblar A. Jackson,
John P. Perdew
Abstract:
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 righ…
▽ More
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.
△ Less
Submitted 25 April, 2020;
originally announced April 2020.
-
Stretched or noded orbital densities and self-interaction correction in density functional theory
Authors:
Chandra Shahi,
Puskar Bhattarai,
Kamal Wagle,
Biswajit Santra,
Sebastian Schwalbe,
Torsten Hahn,
Jens Kortus,
Koblar A. Jackson,
Juan E. Peralta,
Kai Trepte,
Susi Lehtola,
Niraj K. Nepal,
Hemanadhan Myneni,
Bimal Neupane,
Santosh Adhikari,
Adrienn Ruzsinszky,
Yoh Yamamoto,
Tunna Baruah,
Rajendra R. Zope,
John P. Perdew
Abstract:
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 al…
▽ More
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.
△ Less
Submitted 2 April, 2019; v1 submitted 1 March, 2019;
originally announced March 2019.