The Kohn–Sham
exchange–correlation potential
of a
given system can be written exactly in terms of several quantities
including the interacting and noninteracting kinetic energy densities
of the electrons. Electronic kinetic energy densities are sharply
discontinuous at the atomic nuclei in molecules, which implies that
exact exchange–correlation potentials of molecules might also
be discontinuous at the nuclear positions. This surprising possibility
is made even more likely by the fact that molecular exchange–correlation
potentials derived from wave functions within Slater-type basis sets
do have jump discontinuities. We deduce analytically the exact behavior
of the relevant quantities near atomic nuclei and conclude that, in
the basis-set limit, all jump discontinuities cancel out nontrivially
so that exact exchange–correlation potentials are continuous
after all.