Solving the Self-Interaction Problem in Kohn–Sham Density Functional Theory: Application to Atoms

In previous work, we proposed a computational methodology that addresses the elimination of the self-interaction error from the Kohn–Sham formulation of the density functional theory. We demonstrated how the exchange potential can be obtained, and presented results of calculations for atomic systems up to Kr carried out within a Cartesian coordinate system. In this paper, we provide complete details of this self-interaction free method formulated in spherical coordinates based on the explicit equidensity basis ansatz. We prove analytically that derivatives obtained using this method satisfy the Virial theorem for spherical orbitals, where the problem can be reduced to one dimension. We present the results of calculations of ground-state energies of atomic systems throughout the periodic table carried out within the exchange-only mode.