where E KS [ρ I ] refers to the Kohn-Sham energy expression evaluated for the density of subsystem I, E elstat int [{ρ I }] collects all electrostatic interactions between the subsystems, and the non-additive exchange-correlation (NAXC) energy is defined in analogy to Eq. ( ) (see Sect. 4.2 for details).