Van der Waals dispersion interactions play a crucial role in nature, determining the structure and functionality of molecular systems and materials, two examples being the deviation from the ideal gas law observed in real gases and the ability of geckos to stick to walls. These inherently quantum-mechanical interactions originate from the non-local electronic correlation among distant densities and their accurate description is crucial for a reliable computational modelling of matter. As a consequence of the underlying nonlocal electronic correlation, Kohn-Sham (KS) Density Functional Theory (DFT) is unable to model dispersion interactions and expensive wavefunction-based methodologies are usually required. In this connection, three-body Symmetry Adapted Perturbation Theory (SAPT) approaches have been shown to provide very accurate evaluations of dispersion energies. However, these methods are computationally challenging for large systems with scalings ranging from O(N 7 ) for HF-SAPT to O(N ) for SAPT(DFT) . To retain the attractive performances of DFT methods, several dispersion-corrected functionals have been proposed where dispersion is added via empirical pairwise London contributions. This simple and computationally appealing approach has proven successful in including dispersion effects within the DFT framework although the semiempirical nature of exchange-correlation functionals may have effects on their overall accuracy. One limitation of the meanfield pairwise approach in modelling dispersion interactions is the impossibility of capturing many-body dispersion (MBD) effects: many evidences have in fact demonstrated the importance of a non-additive treatment of dispersion interactions in modelling supramolecular complexes , clusters of atoms , one-dimensional wires , extended systems as well as in molecular crystals. The non-additive long-range character of dispersion interactions has been modelled, among the many different approaches , via a set of coupled fluctuating dipoles (CFD) or alternatively by quantum Drude oscillators , the former being preferred over the latter for computational reasons. In recent years, Tkatchenko et al. have formulated the MBD@rsSCS model based on the CFD and introducing the rangeseparation of the self-consistent screening (rsSCS) of polarizabilities which are related to each of the quantum harmonic oscillator in the system. In particular, the model relies on ab initio derived input quantities thus almost removing the presence of empirical parameters. Moreover, the self-consistent screening procedure improves the physical description of the polarizability anisotropy which strongly depends on the atomic environment.