No Arabic abstract
A detailed account of the Kohn-Sham algorithm from quantum chemistry, formulated rigorously in the very general setting of convex analysis on Banach spaces, is given here. Starting from a Levy-Lieb-type functional, its convex and lower semi-continuous extension is regularized to obtain differentiability. This extra layer allows to rigorously introduce, in contrast to the common unregularized approach, a well-defined Kohn-Sham iteration scheme. Convergence in a weak sense is then proven. This generalized formulation is applicable to a wide range of different density-functional theories and possibly even to models outside of quantum mechanics.
The Hohenberg-Kohn theorem plays a fundamental role in density functional theory, which has become a basic tool for the study of electronic structure of matter. In this article, we study the Hohenberg-Kohn theorem for a class of external potentials based on a unique continuation principle.
A Kohn-Sham (KS) inversion determines a KS potential and orbitals corresponding to a given electron density, a procedure that has applications in developing and evaluating functionals used in density functional theory. Despite the utility of KS
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.
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.
It has recently been observed [Phys. Rev. Lett. 113, 113002 (2014)] that the ground-state energy may be obtained directly as a simple sum of augmented Kohn-Sham orbital energies, where it was ascertained that the corresponding one-body shifted Kohn-Sham effective potential has appealing features. With this in mind, eigenvalue and virial constraints are deduced for approximating this potential.