Unfortunately, the quantification of interaction energies within these approaches requires up to sixdimensional numerical integrations over the irregularly-shaped atomic basins, limiting the efficiency and accuracy of the approach. Addressing these limitations remains a focus for advancing real space energy decomposition schemes. An alternative strategy to study chemical interactions in complex systems relies on fragmentation methods in which the total energy of a system is expressed as a sum of contributions from all possible n-tuples of fragments, e.g., monomers, dimers, trimers and so on. As this expansion needs to be truncated in practical applications, the accuracy of such computations is strongly affected by the choice of the computational strategy used for incorporating environmental effects. Hence, a generally applicable quantum chemical method that is capable of providing a quantitative measure of the strength of individual chemical interactions to the energy of a complex chemical system is still lacking.