Let’s ignore the fact that the length scales of planets and the Sun would make the quantum mechanical considerations of DFT irrelevant. In that case, we can replace the Coulomb-interaction F = kqQ/r^2 with the gravitational analogue F=GmM/r^2. DFT could be used to generate a (time-averaged) mass density with the lowest energy for the external potential generated by the sun, by using a wave function consisting of a 1-body wave function for each planet and using a Kohn-Sham operator where the Coulomb terms have been replaced by gravitational analogues.