-
Exploring control of the emergent exciton insulator state in 1T-TiSe$_2$ monolayer by state-of-the-art theory models
Authors:
Hong Tang,
Li Yin,
Gábor I. Csonka,
Adrienn Ruzsinszky
Abstract:
The layered transition metal dichalcogenide 1T-TiSe$_2$ is of great research interest, having intriguing properties of charge density waves (CDW) and superconductivity under doping or pressurizing. The monolayer form of 1T-TiSe$_2$ also shows a CDW with a higher transition temperature T_c than the bulk, indicating a stronger CDW interaction. By using the meta-generalized gradient approximation (me…
▽ More
The layered transition metal dichalcogenide 1T-TiSe$_2$ is of great research interest, having intriguing properties of charge density waves (CDW) and superconductivity under doping or pressurizing. The monolayer form of 1T-TiSe$_2$ also shows a CDW with a higher transition temperature T_c than the bulk, indicating a stronger CDW interaction. By using the meta-generalized gradient approximation (metaGGA)-based model Bethe-Salpeter Equation (BSE) and many-body perturbation GW+BSE methods, we calculate the exciton binding energies and electron energy loss spectrum (EELS) for the 1T-TiSe$_2$ monolayer under different in-plane biaxial strains. We find that even without strain the 1T-TiSe$_2$ monolayer can have negative exciton energies at the Brillouin zone boundary point M, with a binding energy larger than the gap. The calculated EELS reinforces this picture, indicating EI (exciton insulator) states in 1T-TiSe$_2$ monolayer even without strain. The Wannier-Mott formula calculations of exciton binding energy corroborate results from GW+BSE. Small compressive strains enhance the EI state, and for tensile strains slightly less than 3%, the EI state in this monolayer persists. At large tensile strains, the material makes a transition to a normal semiconductor. Our results provide important information for understanding the quantum nature of this two-dimensional (2D) material. Our results from the standard G0W0@PBE+SOC+U+BSE approach are not qualitatively different from those of a more computationally efficient metaGGA-based SCAN+SOC+U+mBSE+$f_{xc}^{loc}$ approach that employs a model BSE.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Tunable band gaps and optical absorption properties of bent MoS$_2$ nanoribbons
Authors:
Hong Tang,
Bimal Neupane,
Santosh Neupane,
Shiqi Ruan,
Niraj K. Nepal,
Adrienn Ruzsinszky
Abstract:
The large tunability of band gaps and optical absorptions of armchair MoS$_2$ nanoribbons of different widths under bending is studied using density functional theory and many-body perturbation GW and Bethe-Salpeter equation approaches. We find that there are two critical bending curvatures, and the non-edge and edge band gaps generally show a non-monotonic trend with bending. The non-degenerate e…
▽ More
The large tunability of band gaps and optical absorptions of armchair MoS$_2$ nanoribbons of different widths under bending is studied using density functional theory and many-body perturbation GW and Bethe-Salpeter equation approaches. We find that there are two critical bending curvatures, and the non-edge and edge band gaps generally show a non-monotonic trend with bending. The non-degenerate edge gap splits show an oscillating feature with ribbon width n, with a period delta_n=3, due to quantum confinement effects. The complex strain patterns on the bent nanoribbons control the varying features of band structures and band gaps that result in varying exciton formations and optical properties. The binding energy and the spin singlet-triplet split of the exciton forming the lowest absorption peak generally decrease with bending curvatures. The large tunability of optical properties of bent MoS$_2$ nanoribbons is promising and will find applications in tunable optoelectronic nanodevices.
△ Less
Submitted 31 October, 2022; v1 submitted 9 September, 2021;
originally announced September 2021.
-
Self-interaction corrected Kohn-Sham effective potentials using the density-consistent effective potential method
Authors:
Carlos M. Diaz,
Luis Basurto,
Santosh Adhikari,
Yoh Yamamoto,
Adrienn Ruzsinszky,
Tunna Baruah,
Rajendra R. Zope
Abstract:
Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn-Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal…
▽ More
Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn-Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. In this work, we adapt and implement the density-consistent effective potential(DCEP) method of Kohut, Ryabinkin, and Staroverov to obtain SI corrected local effective potentials from the SI corrected Fermi-Löwdin orbitals and density in the FLOSIC scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO-LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.
△ Less
Submitted 28 July, 2021;
originally announced July 2021.
-
First-principles wavevector- and frequency-dependent exchange-correlation kernel for jellium at all densities
Authors:
Aaron D. Kaplan,
Niraj K. Nepal,
Adrienn Ruzsinszky,
Pietro Ballone,
John P. Perdew
Abstract:
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…
▽ More
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,ω)$ differs from MCP07 only for non-zero frequencies $ω$. 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πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.}
△ Less
Submitted 22 November, 2021; v1 submitted 5 July, 2021;
originally announced July 2021.
-
Progress towards understanding ultranonlocality through the wavevector and frequency dependence of approximate exchange-correlation kernels
Authors:
Niraj K. Nepal,
Aaron D. Kaplan,
J. M. Pitarke,
Adrienn Ruzsinszky
Abstract:
In the framework of time-dependent density functional theory (TDDFT), the exact exchange-correlation (xc) kernel $f_{xc}(n,q,ω)$ determines the ground-state energy, excited-state energies, lifetimes, and the time-dependent linear density response of any many-electron system. The recently developed MCP07 xc kernel $f_{xc}(n,q,ω)$ of A. Ruzsinszky et al. [Phys. Rev. B 101, 245135 (2020)] yields exce…
▽ More
In the framework of time-dependent density functional theory (TDDFT), the exact exchange-correlation (xc) kernel $f_{xc}(n,q,ω)$ determines the ground-state energy, excited-state energies, lifetimes, and the time-dependent linear density response of any many-electron system. The recently developed MCP07 xc kernel $f_{xc}(n,q,ω)$ of A. Ruzsinszky et al. [Phys. Rev. B 101, 245135 (2020)] yields excellent uniform electron gas (UEG) ground-state energies and plausible plasmon lifetimes. As MCP07 is constructed to describe $f_{xc}$ of the UEG, it cannot capture optical properties of real materials. To verify this claim, we follow Nazarov et al. [Phys. Rev. Lett. 102, 113001 (2009)] to construct the long-range, dynamic xc kernel, $\lim_{q\to 0}f_{xc}(n,q,ω) = -α(ω)e^2/q^2$, of a weakly inhomogeneous electron gas, using MCP07 and other common xc kernels. The strong wavevector and frequency dependence of the "ultranonlocality" coefficient $α(ω)$ is demonstrated for a variety of simple metals and semiconductors. We examine how imposing exact constraints on an approximate kernel shapes $α(ω)$. Comparisons to kernels derived from correlated-wavefunction calculations are drawn.
△ Less
Submitted 30 July, 2021; v1 submitted 24 May, 2021;
originally announced May 2021.
-
Exploring and enhancing the accuracy of interior-scaled Perdew-Zunger self-interaction correction
Authors:
Puskar Bhattarai,
Biswajit Santra,
Kamal Wagle,
Yoh Yamamoto,
Rajendra R. Zope,
Adrienn Ruzsinszky,
Koblar Alan Jackson,
John P. Perdew
Abstract:
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…
▽ More
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σ, 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.
△ Less
Submitted 7 January, 2021;
originally announced January 2021.
-
The Fermi-Löwdin self-interaction correction for ionization energies of organic molecules
Authors:
Santosh Adhikari,
Biswajit Santra,
Shiqi Ruan,
Puskar Bhattarai,
Niraj K. Nepal,
Koblar A. Jackson,
Adrienn Ruzsinszky
Abstract:
(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest-occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. The inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-…
▽ More
(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest-occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. The inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction corrections on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zσ, only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.
△ Less
Submitted 11 August, 2020;
originally announced August 2020.
-
Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories
Authors:
John P. Perdew,
Adrienn Ruzsinszky,
Jianwei Sun,
Niraj K. Nepal,
Aaron D. Kaplan
Abstract:
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…
▽ More
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.
△ Less
Submitted 23 July, 2020;
originally announced July 2020.
-
Constraint-based Wavevector- and Frequency-dependent Exchange-Correlation Kernel of the Uniform Electron Gas
Authors:
Adrienn Ruzsinszky,
Niraj K. Nepal,
J. M. Pitarke,
John P. Perdew
Abstract:
According to time-dependent density functional theory, the exact exchange-correlation kernel f$_{xc}$(n, q, $ω$) 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$π$r$^3_s$). Here we propose a parametrization of this function based upon the satisfaction of exact const…
▽ More
According to time-dependent density functional theory, the exact exchange-correlation kernel f$_{xc}$(n, q, $ω$) 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$π$r$^3_s$). Here we propose a parametrization of this function based upon the satisfaction of exact constraints. For the static ($ω$ = 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 $ω$ at q $=$ 0, we use the model of Gross, Kohn, and Iwamoto. A Cauchy integral extends this model to complex $ω$ 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 $ω$. We then combine these ingredients by damping out the $ω$ dependence at large q in the same way that the q dependence is damped. Away from q $=$ 0 and $ω$ $=$ 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 $ω$ is essentially exact. The plasmon pole of the density response function is found by analytic continuation of f$_{xc}$ to $ω$ 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.
△ Less
Submitted 20 April, 2020;
originally announced April 2020.
-
Examining the order-of-limits problem and lattice constant performance of the Tao--Mo Functional
Authors:
James W Furness,
Niladri Sengupta,
Jinliang Ning,
Adrienn Ruzsinszky,
Jianwei Sun
Abstract:
In their recent communication [Phys. Rev. Lett., 117, 073001 (2016)] Tao and Mo presented a semi-local density functional derived from the density matrix expansion of the exchange hole localised by a general coordinate transformation. We show that the order-of-limits problem present in the functional, dismissed as harmless in the original publication, causes severe errors in predicted phase transi…
▽ More
In their recent communication [Phys. Rev. Lett., 117, 073001 (2016)] Tao and Mo presented a semi-local density functional derived from the density matrix expansion of the exchange hole localised by a general coordinate transformation. We show that the order-of-limits problem present in the functional, dismissed as harmless in the original publication, causes severe errors in predicted phase transition pressures. We also show that the claim that lattice volume prediction accuracy exceeds that of existing similar functionals was based on comparison to reference data that misses anharmonic zero-point expansion and consequently overestimates accuracy. By highlighting these omissions, we give a more accurate assessment of the Tao-Mo functional and show a simple route to resolving the problems.
△ Less
Submitted 3 April, 2020;
originally announced April 2020.
-
Revealing the role of van der Waals interactions in thiophene adsorption on copper surfaces
Authors:
Abhirup Patra,
Adrienn Ruzsinszky
Abstract:
Accurate modeling of electronic and structural properties of organic molecule-metal interfaces are challenging problems because of the complicated electronic distribution of molecule and screening of charges at the metallic surface. This is also the reason why the organic/inorganic system can be engineered for several applications by fine-tuning the metallic work function. Here, we use density-fun…
▽ More
Accurate modeling of electronic and structural properties of organic molecule-metal interfaces are challenging problems because of the complicated electronic distribution of molecule and screening of charges at the metallic surface. This is also the reason why the organic/inorganic system can be engineered for several applications by fine-tuning the metallic work function. Here, we use density-functional theory (DFT) calculations with different level of functional approximations for a systematic study of thiophene interacting with Copper surfaces. In particular, we considered adsorbed structures with the thiophene molecule seated on the top site, with the S atom of the molecule located on the top of a Cu atom. In this work, we find that the weak chemisorption hypothesis of thiophene binding on the copper surface is well justified by the two meta-GGAs-based approximations, SCAN and SCAN+rVV10. PBE-GGA and TM meta-GGA describe it as a physisorption phenomenon by significantly underestimating the adsorption energies. Calculated adsorption energy curves reveal that non-local dispersion interaction between the molecule and metallic surface predominantly controls the bonding mechanism and thus, modifies the copper's work function. Our results imply that semi-local functionals without any kind of van der Waals (vdW) correction can often misinterpret this as physisorption, while, a fortuitous error cancellation can give a right description of this adsorption picture for a wrong reason as in the case of SCAN. The calculated density of states of the adsorbed molecule shows that the long-range vdW correction of SCAN+rVV10 causes more than enough hybridization between the \textit{p} orbitals of S atom and the copper \textit{d}-bands and therefore overestimates the adsorption energies by an average of 16\%.
△ Less
Submitted 23 July, 2019;
originally announced July 2019.
-
A simple self-interaction correction to RPA-like correlation energies
Authors:
Tim Gould,
Adrienn Ruzsinszky,
John P. Perdew
Abstract:
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 b…
▽ More
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.
△ Less
Submitted 12 May, 2019;
originally announced May 2019.
-
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.
-
The key role of antibonding electron transfer in surface chemisorption and heterogeneous catalysis
Authors:
Liping Yu,
Qimin Yan,
Adrienn Ruzsinszky
Abstract:
The description of the chemical bond between a solid surface and an atom or a molecule is the fundamental basis for understanding a broad range of scientific problems in heterogeneous catalysis, semiconductor device fabrication, and fuel cells. Widespread understandings are based on the molecular orbital theory and focused on the degree of filling of antibonding surface-adsorbate states that weake…
▽ More
The description of the chemical bond between a solid surface and an atom or a molecule is the fundamental basis for understanding a broad range of scientific problems in heterogeneous catalysis, semiconductor device fabrication, and fuel cells. Widespread understandings are based on the molecular orbital theory and focused on the degree of filling of antibonding surface-adsorbate states that weaken bonding on surfaces. The unoccupied antibonding surface-adsorbate states are often tacitly assumed to be irrelevant. Here, we show that most antibonding states become unoccupied because the electrons that would occupy these antibonding states are transferred to the lower-energy Fermi level. Such antibonding electron transfer goes beyond molecular orbital theory. It leads to an energy gain that largely controls the trends of surface adsorption strength and can serve as a primary descriptor for bonding on surfaces. This finding is illustrated from the first-principles study of hydrogen adsorption on MoS$_2$ surfaces. A clear linear relationship between the energies of antibonding electron transfer and hydrogen adsorption is identified. The hydrogen evolution reaction on MoS$_2$ is found to originate from the in-gap states induced by sulfur vacancies or edges. The effects of surface inhomogeneity (e.g., defects, step edges) on surface catalysis can be understood from the corresponding different energies gained from antibonding electron transfer. The emerging picture also offers a physically different explanation for the well-known $d$-band theory for hydrogen adsorption on transition metal surfaces.
△ Less
Submitted 5 December, 2018;
originally announced December 2018.
-
Understanding Band Gaps of Solids in Generalized Kohn-Sham Theory
Authors:
John P. Perdew,
Weitao Yang,
Kieron Burke,
Zenghui Yang,
Eberhard K. U. Gross,
Matthias Scheffler,
Gustavo E. Scuseria,
Thomas M. Henderson,
Igor Ying Zhang,
Adrienn Ruzsinszky,
Haowei Peng,
Jianwei Sun,
Egor Trushin,
Andreas Görling
Abstract:
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 theor…
▽ More
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.
△ Less
Submitted 7 March, 2017; v1 submitted 24 August, 2016;
originally announced August 2016.
-
Benchmark tests of a strongly constrained semilocal functional with a long-range dispersion correction
Authors:
J. G. Brandenburg,
J. E. Bates,
A. Ruzsinszky,
J. Sun,
J. P. Perdew
Abstract:
The strongly constrained and appropriately normed (SCAN) semilocal density functional [J. Sun, A. Ruzsinszky, J. P. Perdew \textit{Phys. Rev. Lett.} {\bf 115}, 036402 (2015)] obeys all 17 known exact constraints for meta-generalized-gradient approximations (meta-GGA) and includes some medium range correlation effects. Long-range London dispersion interactions are still missing, but can be accounte…
▽ More
The strongly constrained and appropriately normed (SCAN) semilocal density functional [J. Sun, A. Ruzsinszky, J. P. Perdew \textit{Phys. Rev. Lett.} {\bf 115}, 036402 (2015)] obeys all 17 known exact constraints for meta-generalized-gradient approximations (meta-GGA) and includes some medium range correlation effects. Long-range London dispersion interactions are still missing, but can be accounted for via an appropriate correction scheme. In this study, we combine SCAN with an efficient London dispersion correction and show that lattice energies of simple organic crystals can be improved with the applied correction by 50\%. The London-dispersion corrected SCAN meta-GGA outperforms all other tested London-dispersion corrected meta-GGAs for molecular geometries. Our new method delivers mean absolute deviations (MADs) for main group bond lengths that are consistently below 1\,pm, rotational constants with MADs of 0.2\%, and noncovalent distances with MADs below 1\%. For a large database of general main group thermochemistry and kinetics, it also delivers a weighted mean absolute deviation below 4 kcal/mol, one of the lowest for long-range corrected meta-GGA functionals. We also discuss some consequences of numerical sensitivity encountered for meta-GGAs.
△ Less
Submitted 23 May, 2016;
originally announced May 2016.
-
Strongly Constrained and Appropriately Normed Semilocal Density Functional
Authors:
Jianwei Sun,
Adrienn Ruzsinszky,
John P. Perdew
Abstract:
The ground-state energy, electron density, and related properties of ordinary matter can be computed efficiently when the exchange-correlation energy as a functional of the density is approximated semilocally. We propose the first meta-GGA (meta-generalized gradient approximation) that is fully constrained, obeying all 17 known exact constraints that a meta-GGA can. It is also exact or nearly exac…
▽ More
The ground-state energy, electron density, and related properties of ordinary matter can be computed efficiently when the exchange-correlation energy as a functional of the density is approximated semilocally. We propose the first meta-GGA (meta-generalized gradient approximation) that is fully constrained, obeying all 17 known exact constraints that a meta-GGA can. It is also exact or nearly exact for a set of appropriate norms, including rare-gas atoms and nonbonded interactions. This SCAN (strongly constrained and appropriately normed) meta-GGA achieves remarkable accuracy for systems where the exact exchange-correlation hole is localized near its electron, and especially for lattice constants and weak interactions.
△ Less
Submitted 23 June, 2015; v1 submitted 12 April, 2015;
originally announced April 2015.
-
Gedanken Densities and Exact Constraints in Density Functional Theory
Authors:
John P. Perdew,
Adrienn Ruzsinszky,
Jianwei Sun,
Kieron Burke
Abstract:
Approximations to the exact density functional for the exchange-correlation energy of a many-electron ground state can be constructed by satisfying constraints that are universal, i.e., valid for all electron densities. Gedanken densities are designed for the purpose of this construction, but need not be realistic. The uniform electron gas is an old gedanken density. Here, we propose a spherical t…
▽ More
Approximations to the exact density functional for the exchange-correlation energy of a many-electron ground state can be constructed by satisfying constraints that are universal, i.e., valid for all electron densities. Gedanken densities are designed for the purpose of this construction, but need not be realistic. The uniform electron gas is an old gedanken density. Here, we propose a spherical two-electron gedanken density in which the dimensionless density gradient can be an arbitrary positive constant wherever the density is non-zero. The Lieb-Oxford lower bound on the exchange energy can be satisfied within a generalized gradient approximation (GGA) by bounding its enhancement factor or simplest GGA exchange-energy density. This enhancement-factor bound is well known to be sufficient, but our gedanken density shows that it is also necessary. The conventional exact exchange-energy density satisfies no such local bound, but energy densities are not unique, and the simplest GGA exchange-energy density is not an approximation to it. We further derive a strongly and optimally tightened bound on the exchange enhancement factor of a two-electron density, which is satisfied by the local density approximation but is violated by all published GGA's or meta-GGA's. Finally, some consequences of the non-uniform density-scaling behavior for the asymptotics of the exchange enhancement factor of a GGA or meta-GGA are given.
△ Less
Submitted 23 April, 2014; v1 submitted 23 March, 2014;
originally announced March 2014.
-
Self Interaction Corrected Density Functional Theory with Unitary Invariance: Applications to Molecules
Authors:
Mark R Pederson,
Adrienn Ruzsinszky,
John P. Perdew
Abstract:
Standard spin-density functionals for the exchange-correlation energy of a many-electron ground state make serious self-interaction errors which can be corrected by the Perdew-Zunger self-interaction correction (SIC). We propose a size-extensive construction of SIC orbitals which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional. The SIC orbital…
▽ More
Standard spin-density functionals for the exchange-correlation energy of a many-electron ground state make serious self-interaction errors which can be corrected by the Perdew-Zunger self-interaction correction (SIC). We propose a size-extensive construction of SIC orbitals which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional. The SIC orbitals are constructed from a unitary transformation that is explicitly dependent on the non-interacting one-particle density matrix. When this SIC is applied to the local spin-density approximation, improvements are found for the atomization energies of molecules.
△ Less
Submitted 26 January, 2014;
originally announced January 2014.
-
Density Functionals that Recognize Covalent, Metallic, and Weak Bonds
Authors:
Jianwei Sun,
Bing Xiao,
Yuan Fang,
Robin Haunschild,
Pan Hao,
Adrienn Ruzsinszky,
Gabor I. Csonka,
Gustavo E. Scuseria,
John P. Perdew
Abstract:
Computationally-efficient semilocal approximations of density functional theory at the level of the local spin density approximation (LSDA) or generalized gradient approximation (GGA) poorly describe weak interactions. We show improved descriptions for weak bonds (without loss of accuracy for strong ones) from a newly-developed semilocal meta-GGA (MGGA), by applying it to molecules, surfaces, and…
▽ More
Computationally-efficient semilocal approximations of density functional theory at the level of the local spin density approximation (LSDA) or generalized gradient approximation (GGA) poorly describe weak interactions. We show improved descriptions for weak bonds (without loss of accuracy for strong ones) from a newly-developed semilocal meta-GGA (MGGA), by applying it to molecules, surfaces, and solids. We argue that this improvement comes from using the right MGGA dimensionless ingredient to recognize all types of orbital overlap.
△ Less
Submitted 3 June, 2013; v1 submitted 22 March, 2013;
originally announced March 2013.
-
Effect of orbital-overlap dependence in density functionals
Authors:
Jianwei Sun,
Bing Xiao,
Adrienn Ruzsinszky
Abstract:
The semilocal meta generalized gradient approximation (MGGA) for the exchange-correlation functional of Kohn-Sham (KS) density functional theory can yield accurate ground-state energies simultaneously for atoms, molecules, surfaces, and solids, due to the inclusion of kinetic energy density as an input. We study for the first time the effect and importance of the dependence of MGGA on the kinetic…
▽ More
The semilocal meta generalized gradient approximation (MGGA) for the exchange-correlation functional of Kohn-Sham (KS) density functional theory can yield accurate ground-state energies simultaneously for atoms, molecules, surfaces, and solids, due to the inclusion of kinetic energy density as an input. We study for the first time the effect and importance of the dependence of MGGA on the kinetic energy density through the dimensionless inhomogeneity parameter, $α$, that characterizes the extent of orbital overlap. This leads to a simple and wholly new MGGA exchange functional, which interpolates between the single-orbital regime, where $α=0$, and the slowly varying density regime, where $α\approx 1$, and then extrapolates to $α\to \infty$. When combined with a variant of the Perdew-Burke-Erzerhof (PBE) GGA correlation, the resulting MGGA performs equally well for atoms, molecules, surfaces, and solids.
△ Less
Submitted 3 May, 2012; v1 submitted 10 March, 2012;
originally announced March 2012.