Dynamical screening of van der Waals interactions in nanostructured solids: Sublimation of fullerenes

Sublimation energy is one of the most important properties of molecular crystals, but it is difficult to study, because the attractive long-range van der Waals (vdW) interaction plays an important role. Here, we apply efficient semilocal density functional theory (DFT), corrected with the dynamically screened vdW interaction (DFT + vdW), the Rutgers-Chalmers nonlocal vdW-DF, and the pairwise-based dispersion-corrected DFT-D2 developed by Grimme and co-workers, to study the sublimation of fullerenes. We find that the short-range part, which accounts for the interaction due to the orbital overlap between fullerenes, is negligibly small. Our calculation shows that there exists a strong screening effect on the vdW interaction arising from the valence electrons of fullerenes. On the other hand, higher-order contributions can be as important as the leading-order term. The reasons are that (i) the surface of fullerene molecules is metallic and thus highly polarizable, (ii) the band gap of fullerene solids is small (less than 2 eV), and (iii) fullerene molecules in the solid phase are so densely packed, yielding the high valence electron density and small equilibrium intermolecular distances (the first nearest neighbor distance is only about 10 Å for C60). However, these two effects make opposite contributions, leading to significant error cancellation between these two contributions. We demonstrate that, by considering higher-order contributions and the dynamical screening, the DFT + vdW method can yield sublimation energies of fullerenes in good agreement with reference values, followed by vdW-DF and DFT-D2. The insights from this study are important for a better understanding of the long-range nature of vdW interactions in nanostructured solids.