(n) db-corr correct the electrostatic interaction energies such that they correspond to the full n-body electron density. Note that for the nuclear-nuclear attraction energy, a correction only appears at first order. The nonadditive kinetic and exchange-correlation energies account for non-classical interactions. In the spirit of subsystem-DFT , these contributions are evaluated using approximate, semi-local density functionals. For further details on the evaluation of the different terms in Eq. ( ), we refer to Refs. We note that the calculation of the density-based energy correction is only done once after the individual subsystem calculations required in the MBE. Therefore, the scaling of the computational effort with cluster size is the same for the eb-MBE(n) and the db-MBE(n), even though in the latter case the need for calculating the subsystem electron densities increases the prefactor .