Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExact-exchange Kohn-Sham potential, surface energy, and work function of jellium slabs
107
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Exact-exchange self-consistent calculations of the Kohn-Sham potential, surface energy, and work function of jellium slabs are reported in the framework of the Optimized Effective Potential (OEP) scheme of Density Functional Theory. In the vacuum side of the jellium surface and at a distance $z$ that is larger than the slab thickness, the exchange-only Kohn-Sham potential is found to be image-like ($sim -e^2/z$) but with a coefficient that differs from that of the classical image potential $V_{im}(z)=-e^2/4z$. The three OEP contributions to the surface energy (kinetic, electrostatic, and exchange) are found to oscillate as a function of the slab thickness, as occurs in the case of the corresponding calculations based on the use of single-particle orbitals and energies obtained in the Local Density Approximation (LDA). The OEP work function presents large quantum size effects that are absent in the LDA and which reflect the intrinsic derivative discontinuity of the exact Kohn-Sham potential.
rate research
Read More
The position-dependent exact-exchange energy per particle $varepsilon_x(z)$ (defined as the interaction between a given electron at $z$ and its exact-exchange hole) at metal surfaces is investigated, by using either jellium slabs or the semi-infinite (SI) jellium model. For jellium slabs, we prove analytically and numerically that in the vacuum region far away from the surface $varepsilon_{x}^{text{Slab}}(z to infty) to - e^{2}/2z$, {it independent} of the bulk electron density, which is exactly half the corresponding exact-exchange potential $V_{x}(z to infty) to - e^2/z$ [Phys. Rev. Lett. {bf 97}, 026802 (2006)] of density-functional theory, as occurs in the case of finite systems. The fitting of $varepsilon_{x}^{text{Slab}}(z)$ to a physically motivated image-like expression is feasible, but the resulting location of the image plane shows strong finite-size oscillations every time a slab discrete energy level becomes occupied. For a semi-infinite jellium, the asymptotic behavior of $varepsilon_{x}^{text{SI}}(z)$ is somehow different. As in the case of jellium slabs $varepsilon_{x}^{text{SI}}(z to infty)$ has an image-like behavior of the form $propto - e^2/z$, but now with a density-dependent coefficient that in general differs from the slab universal coefficient 1/2. Our numerical estimates for this coefficient agree with two previous analytical estimates for the same. For an arbitrary finite thickness of a jellium slab, we find that the asymptotic limits of $varepsilon_{x}^{text{Slab}}(z)$ and $varepsilon_{x}^{text{SI}}(z)$ only coincide in the low-density limit ($r_s to infty$), where the density-dependent coefficient of the semi-infinite jellium approaches the slab {it universal} coefficient 1/2.
The behavior of the surface barrier that forms at the metal-vacuum interface is important for several fields of surface science. Within the Density Functional Theory framework, this surface barrier has two non-trivial components: exchange and correlation. Exact results are provided for the exchange component, for a jellium metal-vacuum interface, in a slab geometry. The Kohn-Sham exact-exchange potential $V_{x}(z)$ has been generated by using the Optimized Effective Potential method, through an accurate numerical solution, imposing the correct boundary condition. It has been proved analytically, and confirmed numerically, that $V_{x}(zto infty)to - e^{2}/z$; this conclusion is not affected by the inclusion of correlation effects. Also, the exact-exchange potential develops a shoulder-like structure close to the interface, on the vacuum side. The issue of the classical image potential is discussed.
We construct exact Kohn-Sham potentials for the ensemble density-functional theory (EDFT) from the ground and excited states of helium. The exchange-correlation (XC) potential is compared with the quasi-local-density approximation and both single determinant and symmetry eigenstate ghost-corrected exact exchange approximations. Symmetry eigenstate Hartree-exchange recovers distinctive features of the exact XC potential and is used to calculate the correlation potential. Unlike the exact case, excitation energies calculated from these approximations depend on ensemble weight, and it is shown that only the symmetry eigenstate method produces an ensemble derivative discontinuity. Differences in asymptotic and near-ground-state behavior of exact and approximate XC potentials are discussed in the context of producing accurate optical gaps.
A long-standing puzzle in density-functional theory is the issue of the long-range behavior of the Kohn-Sham exchange-correlation potential at metal surfaces. As an important step towards its solution, it is proved here, through a rigurouos asymptotic analysis and accurate numerical solution of the Optimized-Effective-Potential integral equation, that the Kohn-Sham exact exchange potential decays as $ln(z)/z$ far into the vacuum side of an {it extended} semi-infinite jellium. In contrast to the situation in {it localized} systems, like atoms, molecules, and slabs, this dominant contribution does not arise from the so-called Slater potential. This exact-exchange result provides a strong constraint on the suitability of approximate correlation-energy functionals.
The one-electron density of a many-electron system is the ground-state density of a one-electron Schrodinger equation. The potential $v$ appearing in this Schrodinger equation can be constructed in two ways: In density functional theory (DFT), $v$ is the sum of the Kohn-Sham (KS) potential and the Pauli potential, where the latter can be explicitly expressed in terms of the KS system of non-interacting electrons. As the KS system is fictitious, this construction is only indirectly related to the interacting many-electron system. In contrast, in the exact electron factorization (EEF), $v$ is a functional of the conditional wavefunction $phi$ that describes the spatial entanglement of the electrons in the interacting system. We compare the two constructions of the potential, provide a physical interpretation of the contributions to $v$ in the EEF, and relate it to DFT. With numerical studies of one-dimensional two- and three-electron systems, we illustrate how features of $phi$ translate to the one-electron potential $v$. We show that a change in $phi$ corresponds to a repulsive contribution to $v$, and we explain step structures of $v$ with a charge transfer encoded in $phi$. Furthermore, we provide analytic formulas for the components of $v$ by using a two-state model. Our work thus presents the mapping of a many-electron system to a one-electron system from another angle and provides insights into what determines the shape of the exact one-electron potential. We expect our findings to be helpful for the search of suitable approximations in DFT and in related theories.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
