To explore how vdW density polarization scales with system size, we consider a homologous sequence of polyaromatic hydrocarbon (PAH) dimers from the L7+ dataset: benzene (C1C1PD), coronene (C2C2PD), circumcoronene (C3C3PD), and circum-circumcoronene (C4C4PD). As illustrated in Figure , Q vdW exhibits a perfect linear correlation with system size. The ratio Q MBD /Q PBE reaches as high as 80 %, underscoring the importance of vdW density polarization in large aromatic systems. In contrast, Q PBE shows signs of saturation in C4C4PD, with the value approximately 30 % lower than expected from linear extrapolation. This saturation is observed for other semi-local and hybrid functionals as well (see Figure in the Supporting Information), indicating the possible limitations of these functionals in reproducing ρ(r) for larger aromatic systems.