Research Article

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

See allHide authors and affiliations

Science  11 Mar 2016:
Vol. 351, Issue 6278, pp. 1171-1176
DOI: 10.1126/science.aae0509

Describing dispersion forces

Dispersion or van der Waals interactions are attractive forces that arise from induced dipoles. They are not seen just in atoms and molecules but also in larger nanostructures and even macroscopic objects. Ambrosetti et al. created a qualitatively correct description of van der Waals interactions between polarizable nanostructures over a wide range of finite distances. This required delocalized electrons that have wavelike electron density fluctuations, unlike the more common approaches with dipoles fixed on atoms. Furthermore, the authors observed an enhancement in the nonlocality of the charge density response on the scale of 10 to 20 nm.

Science, this issue p. 1171


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.

The assembly of complex nanostructures and biological systems from simpler building blocks is often driven by noncovalent van der Waals (vdW) or dispersion interactions that arise from electrodynamic correlations between instantaneous charge fluctuations in matter (1, 2). The influence of vdW forces extends well beyond binding energies and encompasses the structural (3, 4), mechanical (5, 6), spectroscopic (7), and even electronic (8) signatures of condensed matter. A common way to characterize vdW interactions is by power laws in the distance D between two or more objects (e.g., atoms, molecules, nanostructures, surfaces, or solids); the most familiar is arguably the Lennard-Jones potential, which is characterized by a short-range repulsive wall with a D–12 dependence and a long-range attractive tail with a D–6 dependence. Even a slight variation in these power laws can have a profound impact on observed properties and therefore demands an accurate, physically sound theoretical description.

Thus far, both our conceptual understanding of vdW interactions and the quantitative models widely used for describing these quantum mechanical phenomena are primarily rooted in low-order intermolecular perturbation theory (IPT), wherein vdW binding originates from the interactions between transient local multipoles (9), and macroscopic Lifshitz theory (10). Although IPT-based approaches have had enormous success in describing vdW binding in (small) gas-phase molecular systems (11, 12), recent advanced experimental techniques have produced several findings that are challenging the basic assumptions of IPT and macroscopic approaches for nanostructured materials, and are strongly indicative that even our qualitative understanding of these interactions is incomplete and needs to be substantially revised (13). Examples of such experimental observations include (i) ultra–long-range vdW interactions extending up to tens of nanometers into heterogeneous dielectric interfaces (14, 15), (ii) complete screening of the vdW interaction between an atomic force microscope (AFM) tip and a SiO2 surface by the presence of one or more layers of graphene adsorbed on the surface (16), (iii) superlinear sticking power laws for the self-assembly of metallic clusters on carbon nanotubes with increasing surface area (17), and (iv) nonlinear increases in the vdW attraction between homologous molecules and an Au(111) surface as a function of molecular size (18). Satisfactory theoretical explanations for these experimental findings either require ad hoc modifications to IPT [(iii) and (iv)] or are inherently outside the domain of applicability of IPT [(i) and (ii)].

To address these issues, we note that the spatial extent of the instantaneous charge density fluctuations responsible for vdW interactions depends rather sensitively on the nature and character of the occupied-to-virtual transitions of the valence electrons in a given molecule or material. In this regard, the inherent delocalization characteristic of low-dimensional metallic systems leads to a remarkably slow decay in the interaction energies between metallic chains and layers as a function of their separation (1921), which represents a substantial modification of the “conventional” asymptotic behavior of the vdW interaction energy that dominates the binding at very large distances (e.g., beyond 10 to 20 nm in bilayer graphene) (22). Moreover, Misquitta et al. (23, 24) demonstrated that upon closure of the band gap, semiconducting nanowires may also exhibit unconventional power laws as a function of the interwire distance D, namely ~D–2 at intermediate separations, followed by asymptotic convergence to the pairwise-additive D–5 limit for large interwire separations.

Here, we revise and extend these seminal asymptotic results to the finite distance regime between nanostructures, because interactions at such distances (as opposed to asymptotic separations) determine the structural, mechanical, spectroscopic, and (opto)electronic properties at the nanoscale. By considering a series of relevant zero-, one-, and two-dimensional systems that include proteins, carbyne-like wires, nanotubes, graphenic layers, and MoS2, we show that a qualitatively correct description of vdW interactions between nonmetallic nanostructures can only be attained by accounting for the wavelike nature of charge density fluctuations in such systems. In doing so, our analysis extends well beyond the low-order particle- or fragment-based IPT paradigm and provides further insight into the roles played by dimensionality, topology, and polarizability in determining the magnitude of these fundamental forces at the nanoscale.

The many-body dispersion approach

In systems where electrons are well described by a localized representation of the occupied space (e.g., systems with finite band gaps), collective charge density fluctuations stem from the dynamically correlated motions of local dipolar excitations. Accordingly, we project the valence electronic response in a given nucleo-electronic system onto a set of N interacting atomic response functions, as outlined by the many-body dispersion (MBD) approach (25). This scheme has been applied with great success to the computation of polarizabilities (26) and dispersion interactions in weakly bound systems (27, 28), and has recently been extended to arbitrary nonmetallic molecules, solids, and nanostructures through the use of spatially distributed polarizabilities (25, 29). The advantage of the MBD method resides in an efficient and accurate quantum mechanical parameterization of the valence electronic response in terms of coupled atomic dipolar fluctuations (30, 31), allowing for a chemically accurate treatment of molecules and extended systems when used in conjunction with state-of-the-art exchange-correlation (XC) functionals in density functional theory (3, 4, 7, 29, 3235).

Within the MBD framework, the long-range correlation energy Embedded Image is computed via the adiabatic connection fluctuation-dissipation theorem (30, 36) asEmbedded Image (1)in which χ0 is the bare response function for a system of noninteracting atoms and χλ is the interacting response function “dressed” by the rescaled and range-separated (25, 29) Coulomb interaction, λv. [We adopt the atomic units convention of me = ħ = e = 1/(4πε0) = 1, where me is the mass of the electron, ħ is the Planck constant divided by 2π, e is the elementary charge, and ε0 is the permittivity of free space.] Within the dipole approximation, the energy expression in Eq. 1 can be computed exactly and is equivalent to diagonalization of the coupled dipolar Hamiltonian (30),Embedded Image (2)in which each atom p is characterized by a static dipole polarizability Embedded Image and a characteristic excitation frequency ωp, and dp represents the mass-weighted displacement of atom p from its equilibrium position Rp. The first two terms in this Hamiltonian correspond to the single-particle kinetic and potential energies, respectively; the last term describes the coupling between atoms p and q via the dipole-dipole interaction tensor Embedded Image, where v(Rpq) is the Coulomb interaction at the interatomic distance Rpq = |RpRq| (30). Embedded Image is then computed as the energetic difference between the eigenfrequencies of the collective modes of the fully interacting system, Embedded Image, and the characteristic excitation frequencies of the isolated atoms, ωp.

Dipole waves in one-dimensional carbyne wires

We begin our investigation into the nature of vdW interactions at the nanoscale with a detailed analysis of the collective charge density fluctuations sustained in a linear one-dimensional (1D) carbyne-like atomic wire comprising 2000 atoms (and subject to periodic boundary conditions). From Fig. 1A, we first note that such finite atomic wires are characterized by a set of very-low-energy collective eigenmodes (albeit with nonzero energy eigenvalues). In this regard, the minimum energy (gap) observed in these modes approaches the independent atom value of ωC = 0.43 Ha (Hartree units) at large C-C distances dC-C, but decreases to approximately 1/20th of this value at dC-C = 1.2 Å, the C-C bond length in acetylene. Because of charge conservation, the quantity Embedded Image (where Embedded Image and Embedded Image denote the static dipole polarizability and resonant frequency of the ith eigenmode, respectively) must be equivalent for every MBD eigenmode, as initially set by the independent atomic value of Embedded Image. Hence, the lowest-energy eigenmode can be characterized by a substantially increased polarizability (i.e., up to ~400 times that of a single C atom), which is strongly indicative of a marked delocalization of the charge density fluctuations over the entire wire. In this case, the low dimensionality and particular topology of this nanostructure is crucial for sustaining coherent delocalized fluctuations along the length of the wire. As depicted in Fig. 1, B and C, these fluctuations result from strongly enhanced (head-to-tail) dipole-dipole coupling along the longitudinal axis of the atomic wire, which leads to substantial anisotropy in the associated polarizability tensor (37).

Fig. 1 Charge density fluctuations in carbyne wires.

(A) Energy eigenvalue spectrum corresponding to the 3N collective eigenmodes of the many-body dispersion (MBD) model (25, 29) plotted in ascending order for 1D carbyne-like atomic wires comprising N = 2000 C atoms (subject to periodic boundary conditions) and values of the interatomic C-C distance dC-C ranging from 1.2 Å (the C-C bond length in acetylene) to 3.0 Å. For reference, the independent atom value of ωC = 0.43 Ha is denoted by the dotted line. The inset shows a more detailed view of the low-energy portion of the spectrum, with the same units. (B) Schematic illustration of several representative low-energy (l), transverse (t), and high-energy (h) collective MBD eigenmodes, as indicated in (A). The colored spheres represent dipole displacements with respect to the equilibrium atomic positions (depicted in gray). (C) Longitudinal dipole displacements of several low-energy collective MBD eigenmodes (in arbitrary units) with respect to the wire coordinate for dC-C = 1.2 Å. Because of the inversion symmetry of this system, every eigenmode is doubly degenerate except for the lowest-energy longitudinal mode (depicted in red). Consistent with a normal mode analysis of a dipolar-coupled system with this underlying topology, the number of nodes present in a given MBD eigenmode is directly correlated with its relative location in the energy eigenvalue spectrum.

To further expand on these observations, we analyzed the low-energy spectrum of the collective MBD eigenmodes in this low-dimensional nanostructure. As illustrated in Fig. 1, B and C, these modes correspond to coherent dipolar fluctuations aligned along the entire wire with negligible components orthogonal to the longitudinal axis. Consistent with a normal mode analysis of a dipolar-coupled system with this underlying topology, the higher-energy modes correspond to polarization waves with an increasing number of nodes. At the highest energies (i.e., beyond the shoulder in the spectrum of Fig. 1A), these eigenmodes have sizable transverse components and therefore make a negligible contribution to the coherent dipolar fluctuations along the wire.

Upon compression of the carbyne wire (i.e., by varying dC-C from 3.0 to 1.2 Å), the gap in the eigenvalue spectrum becomes visibly reduced and is accompanied by a corresponding increase in the slope of the dispersion curves (Fig. 1A). A gapless metallic dispersion (21) is thus approached, but never reached, in the MBD model—a fact that we attribute to the intrinsic charge confinement of the valence atomic responses. A flat dispersion (approached only in the limit dC-C → ∞) would indicate a localization of the system response; with every eigenmode degenerate in energy, single-atom dipolar fluctuations would occur along the wire and result in a decorrelation of the atomic susceptibilities. Conversely, we observe from Fig. 1A that intrawire interactions can also induce nontrivial dispersion in the eigenvalue spectrum despite the absence of explicitly delocalized electrons, implying a marked nonlocality in the collective dipolar response within the MBD model. This collective behavior stems from a subtle interplay between the kinetic and potential energy operators in the MBD Hamiltonian (see. Eq. 2). Upon compression of the carbyne wire, the kinetic energy term increases relative to the potential energy terms, leading to large-amplitude oscillatory motions induced by the presence of many-body kinetic energy contributions. In addition, an analysis of the exact two-atom MBD solution reveals an underlying dependence of Embedded Image on the quantity Embedded Image; as such, the effect of varying the interatomic distance is expected to be qualitatively equivalent to an inversely proportional modification of the polarizability.

vdW interactions between carbyne wires

The analysis presented above already provides strong evidence of the importance of wavelike charge density fluctuations in accounting for vdW interactions in low-dimensional nonmetallic nanostructures. To investigate this aspect further, we now consider the case of two mutually interacting parallel carbyne-like wires. In Fig. 2, the resulting vdW power-law exponents are plotted as a function of the interwire separation D (for several values of dC-C); they exhibit strong deviations from the asymptotic D–5 behavior predicted by widely used vdW approaches based on pairwise additivity. This macroscopic effect goes beyond a simple renormalization of Hamaker constants and results from a theoretical account of both intra- and interwire many-body vdW interactions on an equal footing. As seen in Fig. 2, these modifications to the vdW interaction energy power laws become even more pronounced upon compression of the carbyne wire, reaching ~D–3 for interwire separations between 1 and 3 nm (with dC-C = 1.2 Å), but still not converging to the pairwise-additive asymptote at interwire separations exceeding 20 nm.

Fig. 2 Power laws for interwire vdW interactions.

Interwire interaction energy power-law exponents (46) for two parallel 1D carbyne-like atomic wires as a function of the interwire separation D and several values of dC-C computed using the MBD model. Over a range of relevant nanoscale distances, there are strong deviations from the well-known asymptotic power law of D–5, as predicted by widely used pairwise-additive vdW models (denoted by the dotted line).

Such power-law modifications are particularly relevant at the nanoscale and can be understood as resulting from long-range electrodynamic coupling between the collective longitudinal modes (dipole waves) delocalized along the entirety of each carbyne-like wire. These wavelike charge density fluctuations are simply absent in localized particle- or fragment-based models for vdW interactions, despite being crucial for an accurate theoretical description of these fundamental interactions in low-dimensional nonmetallic nanostructures. These findings resemble the results obtained with approximate many-body models of metallic chains that explicitly account for the delocalized wavelike nature of metallic electrons (5, 21, 23, 38), which implies that a collective dipolar response with markedly nonlocal character can also be used to predict the emergence of wavelike dynamical electron correlation in nonmetallic systems.

To better understand the influence of these fluctuations on the vdW interaction between carbyne wires, we now consider an analytical model for this interaction in more detail. At sufficiently large D, the interwire Coulomb interaction in the continuum approximation is given by w12 = 2K0(qD), where K0 is a modified (or hyperbolic) Bessel function of the second kind and q is the magnitude of the wave vector parallel to the longitudinal wire axis (21). Within the random-phase approximation, one can account for the effects of the interwire coupling on χ1, the response function of a single wire, via Embedded Image, which predicts a splitting of Embedded Image into eigenfrequencies that correspond to coupled dipolar fluctuations that are either aligned (Embedded Image) or anti-aligned (Embedded Image) with respect to the two parallel wires [i.e., Embedded Image, wherein T(q) is the interwire dipole-dipole interaction derived from w12]. Neglecting the higher-energy transverse eigenmodes that provide smaller contributions to the interwire interaction, the interaction energy between carbyne-like wires of length L can then be computed in the nonretarded regime asEmbedded Image (3)which in turn can be expanded to second order in the interaction at large D, yieldingEmbedded Image (4)A q-space decomposition of the integrand in Eq. 4 is provided in Fig. 3 and illustrates that the power law governing the interwire interaction energy has a nontrivial dependence on D that originates from a summation over individual reciprocal-space contributions, each of which corresponds to a normal mode of the system with a different characteristic decay rate. This observation is in stark contrast to the determination of the interwire interaction energy from approximate particle- or fragment-based pairwise approaches, in which each contribution has the same D–5 power-law decay. As such, this analysis is also strongly indicative that the complex power-law exponents displayed in Fig. 2 can only be correctly described by models that allow for delocalized wavelike fluctuations of the charge density. In this regard, even state-of-the-art nonlocal XC functionals (39, 40), which describe long-range electron correlation effects in a pairwise-additive fashion, would completely miss the interaction energy profiles in Fig. 2.

Fig. 3 Analytical model for interwire vdW interactions.

A reciprocal-space decomposition analysis of the integrand Embedded Image, which is required for computing the second-order interwire interaction energy in Eq. 4. Assuming linear dispersion, the decay rate of the integrand with respect to D exhibits a nontrivial dependence on the momentum q (in nm–1), which is responsible for the observed power laws. At large D, only small q contributions substantially influence Embedded Image, so that D effectively determines the scale of relevant momenta. In this plot, the integrand has been rescaled by the constant Embedded Image and ithe corresponding color scale is provided in atomic units.

To gain a more detailed understanding of the power-law exponents displayed in Fig. 2, we note that the asymptotic behavior of the interwire interaction energy is largely influenced by Embedded Image, the single-wire dispersion appearing in the integrand of Eq. 4. In fact, a first estimate of the decay rate of Embedded Image can be obtained by fitting the dispersion law with a single exponent—that is, as Embedded Image—which directly leads to the finding that Embedded Image upon variable substitution of q′ = q/D. Physically speaking, the quantity K0(qD) governs the interwire interaction by introducing a momentum cutoff at q ≈ 1/D, and the resultant integration over the infinitesimal contributions Embedded Image up to this momentum cutoff controls the scaling of Embedded Image with respect to D. In particular, a linear plasmon dispersion leads (up to logarithmic corrections) to a D–2 decay (21) in metallic wires, whereas the slightly sublinear dispersion observed here for dC-C = 1.2 Å (see the inset in Fig. 1A after the initial energetic gap) is consistent with a ~D–2.8 power law and is in qualitative agreement with the numerical evidence provided by the MBD model in Fig. 2. Despite this initial deviation from D–5, the power law at large distances again tends toward the pairwise-additive asymptote, although deviations still persist at interwire separations greater than 20 nm. A direct account of this phenomenon follows from Fig. 3: Because of the decay of K0(qD) with respect to the quantity qD, only very small values of q will contribute to the interwire interaction energy at large D (see Eq. 4). In fact, only those wave vectors near the gap [in which ω(q) ≈ constant] will be of relevance in recovering the asymptotic D–5 limit. Hence, the pairwise-additive limit can only be approached in the presence of a flat dispersion, wherein the localization of the system response to single-atom dipolar fluctuations occurring along the individual wires would actually validate the fundamental particle- or fragment-based assumption that the vdW interaction energy can be described as a summation over induced atomic dipole contributions.

Note that the results of the above analysis are quite general and are not simply intrinsic properties of the MBD coupled dipolar Hamiltonian (Eq. 2). The same conclusions are observed when the qualitatively different tight binding (TB) model (23) is used for the single wire response. Within this approach, the individual atomic wires are described by a two-site nearest-neighbor interaction that can be tuned (by manipulating the β1 and β2 free-parameter space) to reproduce both the metallic (symmetric interactions via β21 = 1) and insulating (asymmetric/single-sided interactions via β21 = 0) limits. A second-order perturbative treatment of the interwire interaction energy leads to power-law exponents that describe these two respective limits—that is, ~D–2 for metallic nanowires and ~D–5 for insulating nanowires with flat energy bands—as depicted in Fig. 4. By analogy to the MBD model, the TB response becomes markedly nonlocal in the metallic limit, the regime in which the largest deviations from the pairwise-additive ~D–5 power law are expected. Conversely, a substantial localization of the response is recovered in the insulating limit, in which the pairwise-additive approximation becomes asymptotically valid. Comparisons between the power-law exponents provided in Figs. 2 and 4 demonstrate semiquantitative agreement between the qualitatively different MBD and TB models, which is strongly indicative of the universality of the influence of wavelike charge density fluctuations on vdW interactions between nanostructures.

Fig. 4 The tight-binding model for interwire vdW interactions.

Interwire interaction energy power-law exponents (46) are shown for two parallel 1D carbyne-like atomic wires as a function of the interwire separation D computed using the tight-binding (TB) model (23, 46). With dC-C = 2.0 Å, the two limiting cases of insulating (β21 = 0) and metallic (β21 = 1) nanowires were considered, along with an intermediate case (β21 = 0.96). As previously observed with the qualitatively different MBD model (see Fig. 2), there are strong deviations from the asymptotic D–5 power law (denoted by the dotted line) over a wide range of relevant nanoscale distances.

vdW interactions between complex nanostructures

Although the results above were obtained for model 1D carbyne-like wires, our calculations suggest that collective charge density fluctuations will also strongly influence the behavior of strained hydrocarbon polymers and inorganic wires made from polarizable elements such as phosphorus, silicon, or germanium. To further elucidate the roles of dimensionality, topology, and polarizability on determining vdW interaction power laws across distances relevant at the nanoscale, we carried out a study of layered two-dimensional (2D) materials, including graphenic sheets of varying polarizability as well as bilayer MoS2. Although a smaller reduction in the power-law exponent is expected in bilayer graphene because of its higher dimensionality, a D–3.5 decay was found for interlayer distances of 0.6 to 3.0 nm, with convergence to the conventional D–4 power law only achieved for interlayer distances well beyond 10 nm (Fig. 5). In general, the interlayer interaction power-law exponents for bilayer graphene were intermediate between D–3 and the pairwise-additive D–4 asymptote, which are induced by gapless π states and are expected to be valid in the case of finite-gap doped graphene (41, 42). Increasing the polarization response of the graphenic sheets to effectively model highly polarizable 2D materials leads to substantial delocalization in the system response; for Embedded Image = 50 bohr3, the power law exhibited an extended plateau between D–2.5 and D–3 at interlayer separations from 1.5 nm to well beyond 10 nm—distances that are again quite relevant to interactions at the nanoscale.

Fig. 5 Power laws for interlayer vdW interactions.

Interlayer interaction energy power-law exponents (46) for two parallel 2D graphenic layers and MoS2 as a function of the interlayer separation D computed using the MBD model. To investigate the behavior of these power laws in highly polarizable 2D materials, we considered several values of Embedded Image for the graphenic layers, ranging from 12 bohr3 (the static atomic dipole polarizability of a C atom) to 50 bohr3. As found above for the case of interacting parallel 1D nanowires over a range of relevant nanoscale distances, there are substantial deviations from the asymptotic D–4 power law predicted by pairwise-additive vdW approaches (denoted by the dotted line).

As a final set of examples, we extend our observations regarding the crucial importance of many-body vdW interactions (43) to the description of complex nanostructures such as carbon nanotubes and biological systems. As Fig. 6 demonstrates, the decay rate of the MBD interaction energy between a nanowire and the 1MC5 protein (human glutathione-dependent formaldehyde dehydrogenase) is much slower than the prediction of the pairwise-additive approximation. In addition, the ratio between the MBD and pairwise interaction energies, Embedded Image, is quite large and grows as a function of the wire-protein separation, demonstrating a nontrivial coupling between delocalized charge density fluctuations sustained across the protein and along the wire.

Fig. 6 Power laws for vdW interactions between complex nanostructures.

Interaction energy power-law exponents (46) computed using the MBD model for two parallel (3,3) carbon nanotubes and a wire-protein nanostructure (depicted in insets) as a function of the respective center-of-mass separations D. As observed above for interacting 1D nanowires and 2D layered materials, we report substantial deviations in the decay rate of the MBD interaction energy in such nanostructures when compared against the pairwise-additive asymptote of D–5 (denoted by the dotted line). The ratio of the MBD and pairwise (12) interaction energies Embedded Image are provided for select values of D; the growth of this ratio with D is another indication of the nontrivial coupling between delocalized charge density fluctuations sustained across these nanostructures.


Our analysis demonstrates the possibilities for engineering the vdW interaction energy power laws in low-dimensional nanostructures by modifying their underlying dimensionality, topology, and response properties. Moreover, our findings suggest that complex power laws could also play an important role in numerous low-dimensional biologically relevant systems, including phospholipid aggregates and bilayers (44) or even the subnanometer spatula-shaped structures that determine the peculiar pedal adhesion in the gecko (45). Our study provides strong evidence that the ubiquitous vdW forces between polarizable nonmetallic nanostructures can be more completely understood in terms of collective interactions between wavelike charge density fluctuations, rather than simply a summation over pairwise interactions between instantaneous particle- or fragment-like dipolar fluctuations.

In a manner analogous to metallic systems, the marked nonlocality of the valence electronic response in low-dimensional polarizable nanostructures is reflected in the emergence of collective modes that span the entire system. It is these delocalized wavelike charge density fluctuations that are responsible for the nontrivial interaction energy power laws observed herein that substantially deviate from the predictions of standard pairwise-additive vdW approximations. As such, these wavelike fluctuations govern the magnitude of vdW interactions at large, but finite, distances of relevance at the nanoscale. Our findings reveal a smooth transition from a system consisting of independent atoms to the collective fully interacting limit, thereby providing a potential pathway for tuning the fundamental noncovalent vdW interactions responsible for the assembly of complex polarizable nanostructures.

Supplementary Materials

Theoretical Methods

Supplementary Text

Tables S1 to S3

References and Notes

  1. See supplementary materials on Science Online.
Acknowledgments: Supported by a startup grant from Cornell University (R.A.D.), European Research Council Starting Grant VDW-CMAT, and DFG/SFB-951 HIOS project A10. A.A. thanks F. Toigo and P. L. Silvestrelli for useful discussions. This research used resources of the Argonne Leadership Computing Facility at Argonne National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under contract DE-AC02-06CH11357.

Stay Connected to Science

Navigate This Article