In this section, we present the results for charge displacements in selected systems as calculated by the HF method and several DFAs. Figure ,b shows the ratio of vdW-displaced charge Q vdW to its PBE and HF counterparts, Q PBE and Q HF , for the six selected systems. Overall, Q HF is about 10 % larger than Q PBE , which modifies the ratio accordingly. Nevertheless, this does not change the conclusions made in the main text based on the PBE densities. For C1-C4 sequence of polyaromatic hydrocarbons (PAH), the additional calculations of charge displacement were performed with LDA, PBE0 and B3LYP xc-functionals as well as with the Hartree-Fock (HF) method. The computational settings were kept the same as described for the PBE case in the main text. The results obtained (Figure ) witness that both (semi-)local and hybrid DFAs exhibit similar saturation of displaced charge moving from C1 to C4, and this is not specific to the PBE functional. The only different case is the HF method, showing almost linear behavior.