Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userElimination of the asymptotic fractional dissociation problem in Kohn-Sham density functional theory using the ensemble-generalization approach
142
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Many approximations within density-functional theory spuriously predict that a many-electron system can dissociate into fractionally charged fragments. Here, we revisit the case of dissociated diatomic molecules, known to exhibit this problem when studied within standard approaches, including the local spin-density approximation (LSDA). By employing our recently proposed [E. Kraisler and L. Kronik, Phys. Rev. Lett. 110, 126403 (2013)] ensemble-generalization we find that asymptotic fractional dissociation is eliminated in all systems examined, even if the underlying exchange-correlation (xc) is still the LSDA. Furthermore, as a result of the ensemble generalization procedure, the Kohn-Sham potential develops a spatial step between the dissociated atoms, reflecting the emergence of the derivative discontinuity in the xc energy functional. This step, predicted in the past for the exact Kohn-Sham potential and observed in some of its more advanced approximate forms, is a desired feature that prevents any fractional charge transfer between the systems fragments. It is usually believed that simple xc approximations such as the LSDA cannot develop this step. Our findings show, however, that ensemble generalization to fractional electron densities automatically introduces the desired step even to the most simple approximate xc functionals and correctly predicts asymptotic integer dissociation.
rate research
Read More
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.
In numerical computations of response properties of electronic systems, the standard model is Kohn-Sham density functional theory (KS-DFT). Here we investigate the mathematical status of the simplest class of excitations in KS-DFT, HOMO-LUMO excitations. We show using concentration-compactness arguments that such excitations, i.e. excited states of the Kohn-Sham Hamiltonian, exist for $Z>N$, where $Z$ is the total nuclear charge and $N$ is the number of electrons. The result applies under realistic assumptions on the exchange-correlation functional, which we verify explicitly for the widely used PZ81 and PW92 functionals. By contrast, and somewhat surprisingly, we find using a method of Glaser, Martin, Grosse, and Thirring cite{glaser1976} that in case of the hydrogen and helium atoms, excited states do not exist in the neutral case $Z=N$ when the self-consistent KS ground state density is replaced by a realistic but easier to analyze approximation (in case of hydrogen, the true Schr{o}dinger ground state density). Implications for interpreting minus the HOMO eigenvalue as an approximation to the ionization potential are indicated.
We develop a density functional treatment of non-interacting abelian anyons, which is capable, in principle, of dealing with a system of a large number of anyons in an external potential. Comparison with exact results for few particles shows that the model captures the behavior qualitatively and semi-quantitatively, especially in the vicinity of the fermionic statistics. We then study anyons with statistics parameter $1+1/n$, which are thought to condense into a superconducting state. An indication of the superconducting behavior is the mean-field result that, for uniform density systems, the ground state energy increases under the application of an external magnetic field independent of its direction. Our density-functional-theory based analysis does not find that to be the case for finite systems of anyons, which can accommodate a weak external magnetic field through density transfer between the bulk and the boundary rather than through transitions across effective Landau levels, but the Meissner repulsion of the external magnetic field is recovered in the thermodynamic limit as the effect of the boundary becomes negligible. We also consider the quantum Hall effect of anyons, and show that its topological properties, such as the charge and statistics of the excitations and the quantized Hall conductance, arise in a self-consistent fashion.
A very specific ensemble of ground and excited states is shown to yield an exact formula for any excitation energy as a simple correction to the energy difference between orbitals of the Kohn-Sham ground state. This alternative scheme avoids either the need to calculate many unoccupied levels as in time-dependent density functional theory (TDDFT) or the need for many self-consistent ensemble calculations. The symmetry-eigenstate Hartree-exchange (SEHX) approximation yields results comparable to standard TDDFT for atoms. With this formalism, SEHX yields approximate double-excitations, which are missed by adiabatic TDDFT.
Recent work has established Moreau-Yosida regularization as a mathematical tool to achieve rigorous functional differentiability in density-functional theory. In this article, we extend this tool to paramagnetic current-density-functional theory, the most common density-functional framework for magnetic field effects. The extension includes a well-defined Kohn-Sham iteration scheme with a partial convergence result. To this end, we rely on a formulation of Moreau-Yosida regularization for reflexive and strictly convex function spaces. The optimal $L^p$-characterization of the paramagnetic current density $L^1cap L^{3/2}$ is derived from the $N$-representability conditions. A crucial prerequisite for the convex formulation of paramagnetic current-density-functional theory, termed compatibility between function spaces for the particle density and the current density, is pointed out and analyzed. Several results about compatible function spaces are given, including their recursive construction. The regularized, exact functionals are calculated numerically for a Kohn-Sham iteration on a quantum ring, illustrating their performance for different regularization parameters.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
