Effective local potentials for density and density-matrix functional approximations with non-negative screening density

A way to improve the accuracy of the spectral properties in density functional theory (DFT) is to impose constraints on the effective, Kohn-Sham (KS), local potential [J. Chem. Phys. 136, 224109 (2012)]. As illustrated, a convenient variational quantity in that approach is the "screening" or "electron repulsion" density, ρ_rep, corresponding to the local, KS Hartree, exchange and correlation potential through Poisson's equation. Two constraints, applied to this minimization, largely remove self-interaction errors from the effective potential: (i) ρ_rep integrates to N - 1, where N is the number of electrons, and (ii) ρ_rep ≥ 0 everywhere. In this work, we introduce an effective "screening" amplitude, f, as the variational quantity, with the screening density being ρ_rep = f². In this way, the positivity condition for ρ_rep is automatically satisfied, and the minimization problem becomes more efficient and robust. We apply this technique to molecular calculations, employing several approximations in DFT and in reduced density matrix functional theory. We find that the proposed development is an accurate, yet robust, variant of the constrained effective potential method.