The χ nad term serves as a mostly positive potential due to the frozen density ρ 0 that acts on the polarizable density ∆ρ and thus forces ∆ρ to be pushed away from overlapping frozen densities. As a result, volume exclusion effects due to frozen-polarizable density overlaps are emulated via the first order non-additive term ∆n ⊺ χ. However, a large portion of this non-additive component originates from kinetic energies, and an analytical form in terms of orbital-free densities is only possible for simple cases such as the Thomas-Fermi or Weizsäcker 60 functionals. Instead, we have found that a simple empirical contribution that is proportional to the frozen density works sufficiently well, and uses the form