The overall computational cost involved in the evaluation of the MBD correlation energy arises from the N ν -times solution of Dyson-like self-consistent screening equations (Eq.( )) as well as from the eigendecomposition of V MBD both of which being as O(N 3 ) in scaling. In the following part of this section we first describe a robust and linear-scaling iterative self-consistent screening algorithm as well as a linear-scaling stochastic estimation of E int based on a state-of-the-art linear algebra numerical algorithm.