Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFrom Kohn-Sham to many-electron energies via step structures in the exchange-correlation potential
155
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Accurately describing excited states within Kohn-Sham (KS) density functional theory (DFT), particularly those which induce ionization and charge transfer, remains a great challenge. Common exchange-correlation (xc) approximations are unreliable for excited states owing, in part, to the absence of a derivative discontinuity in the xc energy ($Delta$), which relates a many-electron energy difference to the corresponding KS energy difference. We demonstrate, analytically and numerically, how the relationship between KS and many-electron energies leads to the step structures observed in the exact xc potential, in four scenarios: electron addition, molecular dissociation, excitation of a finite system, and charge transfer. We further show that steps in the potential can be obtained also with common xc approximations, as simple as the LDA, when addressed from the ensemble perspective. The article therefore highlights how capturing the relationship between KS and many-electron energies with advanced xc approximations is crucial for accurately calculating excitations, as well as the ground-state density and energy of systems which consist of distinct subsystems.
rate research
Read More
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 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 model a Kohn-Sham potential with a discontinuity at integer particle numbers derived from the GLLB approximation of Gritsenko et al. We evaluate the Kohn-Sham gap and the discontinuity to obtain the quasiparticle gap. This allows us to compare the Kohn-Sham gaps to those obtained by accurate many-body perturbation theory based optimized potential methods. In addition, the resulting quasiparticle band gap is compared to experimental gaps. In the GLLB model potential, the exchange-correlation hole is modeled using a GGA energy density and the response of the hole to density variations is evaluated by using the common-denominator approximation and homogeneous electron gas based assumptions. In our modification, we have chosen the PBEsol potential as the GGA to model the exchange hole, and add a consistent correlation potential. The method is implemented in the GPAW code, which allows efficient parallelization to study large systems. A fair agreement for Kohn-Sham and the quasiparticle band gaps with semiconductors and other band gap materials is obtained with a potential which is as fast as GGA to calculate.
Using a simplified one-dimensional model of a diatomic molecule, the associated interacting density and corresponding Kohn-Sham potential have been obtained analytically for all fractional molecule occupancies $N$ between 0 and 2. For the homonuclear case, and in the dissociation limit, the exact Kohn-Sham potential builds a barrier at the midpoint between the two atoms, whose strength increases linearly with $N$, with $1 < N leq 2$. In the heteronuclear case, the disociating KS potential besides the barrier also exhibits a plateau around the atom with the higher ionization potential, whose size (but not its strength) depends on $N$. An anomalous zero-order scaling of the Kohn-Sham potential with regards to the strength of the electron-electron repulsion is clearly displayed by our model; without this property both the unusual barrier and plateau features will be absent.
In this work we introduce a new semi-implicit second order correction scheme to the kinetic Kohn-Sham lattice model. The new approach is validated by performing realistic exchange-correlation energy calculations of atoms and dimers of the first two rows of the periodic table finding good agreement with the expected values. Additionally we simulate the ethane molecule where we recover the bond lengths and compare the results with standard methods. Finally, we discuss the current applicability of pseudopotentials within the lattice kinetic Kohn-Sham approach.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
