To assess the extent of vdW charge density polarization in real systems, we examine two sets of large organic complexes. The first one consists of the 12 host-guest dimers from the S12L dataset . The second, referred to as the L7+ dataset, includes 11 systems: seven large complexes from the L7 dataset along with four addi- tional systems -C 60 buckyball dimer, parallel-displaced stacked dimers of circumcoronene (C3C3PD) and circumcircumcoronene (C4C4PD), and a DNA-ellipticine intercalation complex. As a reference, we also include a parallel-displaced stacked benzene dimer in both datasets (as complex '1' in S12L and C1C1PD in L7+), representing the prototypical π -π interaction. In total, we consider 24 unique systems, ranging from 24 to 240 atoms. For consistency, the atomic geometries and nomenclature follow the original publications . The additional geometries of the C 60 dimer and the DNA-ellipticine complex were taken from Ref. , while the atomic geometries of C3C3PD and C4C4PD were optimized at the PBE+MBD@rsSCS level in FHI-aims with 'tight' basis sets. The vdW-induced charge displacements, ∆ρ pol (r), were computed for all systems using both the sc-TS and MBD methods, and compared against the PBE functional. The choice of PBE is justified by its widespread use and relatively low computational cost compared to more accurate meta-GGA or hybrid functionals. Test calculations with Hartree-Fock (HF) for selected systems yielded qualitatively similar results, supporting the generality of our conclusions (see Figure ).