In our implementation, the contributions of the electron-nuclei attraction and of the Coulomb repulsion are evaluated recursively order-by-order. As long as the contributions of the nonadditive kinetic and exchange-correlation energies [Eqs. ( ) and ( )] are approximate using (semi-)local density functionals, all terms in Eq. ( ) can be evaluated using the coordinates of the nuclei and the electron densities of the monomers, dimers, trimers etc. only. Therefore, while in our previous work, we applied the db-MBE only in combination with DFT calculations, it can be extended to any quantum-chemical method that provides an electron density, including wavefunction-based methods such as coupled cluster theory.