Minimization of Eq. ( ) with respect to the electron density ρ K of subsystem K (subject to the condition that ρ K can be represented by a non-interacting system with Kohn-Sham orbitals {ψ i K }) while keeping all other subsystem electron densities fixed leads to Kohn-Sham-like equations for subsystem K,