I would consider a set of independent Kohn–Sham–like equations for each planet.
Each planet has its own kinetic energy and is subject to the gravitational attraction from the Sun, which plays the role of the “bare nucleus” of the system.
Each planet does not directly interact with the others but instead experiences an average gravitational potential arising from their collective presence.
Then we included a exchange–correlation potential as for electronic DFT.