Communication: Non-additivity of van der Waals interactions between nanostructures

Electron confinement meet electron delocalization: non-additivity and finite-size effects in the polarizabilities and dispersion coefficients of the fullerenes

In this work, we used finite-field derivative techniques and density functional theory (DFT) to compute the static isotropic polarizability series (αl with l = 1, 2, 3) for the C60-C84 fullerenes and quantitatively assess the intrinsic non-additivity in these fundamental response properties. By comparing against classical models of the fullerenes as conducting spherical shells (or solid spheres) of uniform electron density, a detailed critical analysis of the derived effective scaling laws (α1 ∼ N1.2, α2 ∼ N2.0, α3 ∼ N2.7) demonstrates that the electronic structure of finite-sized fullerenes-a unique dichotomy of electron confinement and delocalization effects due to their quasi-spherical cage-like structures and encapsulated void spaces-simultaneously limits and enhances their quantum mechanical response to electric field perturbations. Corresponding frequency-dependent polarizabilities were obtained by inputting the αl series into the hollow sphere model (within the modified single frequency approximation), and used to compute the molecular dispersion coefficients (Cn with n = 6, 8, 9, 10) needed to describe the non-trivial van der Waals (vdW) interactions in fullerene-based systems. Using first-order perturbation theory in conjunction with >140 000 DFT calculations, we also computed the non-negligible zero-point vibrational contributions to α1 in C60 and C70, thereby enabling a more accurate and direct comparison between theory and experiment for these quintessential nanostructures.

Origin of the size-dependence of the equilibrium van der Waals binding between nanostructures

Nanostructures can be bound together at equilibrium by the van der Waals (vdW) effect, a small but ubiquitous many-body attraction that presents challenges to density functional theory. How does the binding energy depend upon the size or number of atoms in one of a pair of identical nanostructures? To answer this question, we treat each nanostructure as a whole object, not as a collection of atoms. Our calculations start from an accurate static dipole polarizability for each considered nanostructure, and an accurate equilibrium center-to-center distance for the pair (the latter from experiment or from the vdW-DF-cx functional). We consider the competition in each term -C_2k/d^2k (k = 3, 4, 5) of the long-range vdW series for the interaction energy, between the size dependence of the vdW coefficient C_2k and that of the 2kth power of the center-to-center distance d. The damping of these vdW terms can be negligible, but in any case, it does not affect the size dependence for a given term in the absence of non-vdW binding. To our surprise, the vdW energy can be size-independent for quasi-spherical nanoclusters bound to one another by vdW interaction, even with strong nonadditivity of the vdW coefficient, as demonstrated for fullerenes. We also show that, for low-dimensional systems, the vdW interaction yields the strongest size-dependence, in stark contrast to that of fullerenes. We illustrate this with parallel planar polycyclic aromatic hydrocarbons. The size dependences of other morphologies or bonding types lie between, as shown by sodium clusters.

Wavelike charge density fluctuations and van der Waals interactions at the nanoscale

Recent experiments on noncovalent interactions at the nanoscale have challenged the basic assumptions of commonly used particle- or fragment-based models for describing van der Waals (vdW) or dispersion forces. We demonstrate that a qualitatively correct description of the vdW interactions between polarizable nanostructures over a wide range of finite distances can only be attained by accounting for the wavelike nature of charge density fluctuations. By considering a diverse set of materials and biological systems with markedly different dimensionalities, topologies, and polarizabilities, we find a visible enhancement in the nonlocality of the charge density response in the range of 10 to 20 nanometers. These collective wavelike fluctuations are responsible for the emergence of nontrivial modifications of the power laws that govern noncovalent interactions at the nanoscale.

Non-pairwise additivity of the leading-order dispersion energy

The leading-order (i.e., dipole-dipole) dispersion energy is calculated for one-dimensional (1D) and two-dimensional (2D) infinite lattices, and an infinite 1D array of infinitely long lines, of doubly occupied locally harmonic wells. The dispersion energy is decomposed into pairwise and non-pairwise additive components. By varying the force constant and separation of the wells, the non-pairwise additive contribution to the dispersion energy is shown to depend on the overlap of density between neighboring wells. As well separation is increased, the non-pairwise additivity of the dispersion energy decays. The different rates of decay for 1D and 2D lattices of wells is explained in terms of a Jacobian effect that influences the number of nearest neighbors. For an array of infinitely long lines of wells spaced 5 bohrs apart, and an inter-well spacing of 3 bohrs within a line, the non-pairwise additive component of the leading-order dispersion energy is -0.11 kJ mol(-1) well(-1), which is 7% of the total. The polarizability of the wells and the density overlap between them are small in comparison to that of the atomic densities that arise from the molecular density partitioning used in post-density-functional theory (DFT) damped dispersion corrections, or DFT-D methods. Therefore, the nonadditivity of the leading-order dispersion observed here is a conservative estimate of that in molecular clusters.