Kohn-Sham equations replace a complex many-body system with simpler non-interacting particles that feel the influence of the global density of all other particles. For planets and the Sun we should then find the average gravitational potential from all other bodies that gives the same overall distribution to transform the problem into single-body equations.