Interpretation of van der Waals density functionals
Abstract
The nonlocal correlation energy in the van der Waals density functional (vdW-DF) method [Phys. Rev. Lett. 92, 246401 (2004); Phys. Rev. B 76, 125112 (2007); Phys. Rev. B 89, 035412 (2014)] can be interpreted in terms of a coupling of zero-point energies of characteristic modes of semilocal exchange-correlation (xc) holes. These xc holes reflect the internal functional in the framework of the vdW-DF method [Phys. Rev. B 82, 081101(2010)]. We explore the internal xc hole components, showing that they share properties with those of the generalized-gradient approximation. We use these results to illustrate the nonlocality in the vdW-DF description and analyze the vdW-DF formulation of nonlocal correlation.
pacs:
31.15.E-, 71.15.MbI Introduction
In a seminal paperra Rapcewicz and Ashcroft (RA) highlighted the connections between nonlocal correlations, the exchange-correlation (xc) hole concepthajo74 ; lape75 ; gulu76 ; lape77 ; gujolu79 ; wang91p13298 of density functional theory (DFT), and van der Waals (vdW) forces in the inhomogeneous electron gas. RA introduced a simple physical picture of vdW binding: electrons and their associated xc holes form neutral pairs in a system resembling condensed, vdW-bounded, atomic matter and experience mutual attraction of a dispersive nature.eilo30 ; lo30 ; jerry65 In the RA view, it is the local plasmon that characterizes the interaction components, i.e., electron-xc-hole pairs. The RA picture is supported by a previous study of the nonlinear response in the electron gas,ma predicting strong vdW binding from quantum-fluctuation contributions in the interaction diagram that also underpins an analysis of gradient-corrected correlation.rasolt ; lape80 ; RasoltDisc1 ; hula86 ; RasoltDisc2 In fact, the long-range interaction component is interpretedlavo87 as reflecting the small-momentum fluctuation components that are extracted to reach a generalized gradient approximation (GGA) in the early formulations.lape75 ; lape77 ; lape80 ; lameprl1981 ; lape82 ; lameprb1983 ; lavo87 ; lavo90 Together Refs. [ra, ; ma, ; lavo87, ] suggest that one can recover vdW forces in nonlocal functional theories that also incorporate the tremendous progress that GGA represents.KieronPerspective ; BeckePerspective
The Rutgers-Chalmers vdW-DF methodanlalu96 ; hurylula99 ; huanlula96 ; ryluladi00 ; rydbergthesis ; rydberg03p126402 ; dionthesis ; dion ; dionerratum ; langreth05p599 ; thonhauser ; cooper10p161104 ; lee10p081101 ; behy13 ; berland14p035412 ; bearcoleluscthhy14 ; ROPPreview allows efficient computations of the xc energy in an approximation that seamlessly incorporates nonlocal correlation effects, including vdW forces. The vdW-DF method is gaining recognition for helping to extend the success of nonempirical DFT to sparse matter.langrethjpcm2009 The vdW-DF method is free from external parameters and rests only on formal theory input,ROPPreview for the local density approximationgulu76 ; perdew92p13244 (LDA) and for gradient-corrected exchangelavo90 ; schwinger ; elliott09p1485 ; lee10p081101 in its specification of the plasmon behavior. It also includes a GGA exchange component.pewa86 ; pebuer96 ; zhya98 ; mulela09 ; cooper10p161104 ; vdwsolids ; berland14p035412 ; hamada14 The choice of vdW-DF exchange can be guided by conservation of the full (nonlocal) xc holeberland14p035412 and with such consistent-exchange vdW-DF it is possible to investigate bulk-structure and adsorption problems where interactions are in subtle competition.berland14p035412 ; bearcoleluscthhy14 The transferability of vdW-DF has also been probed via comparison with quantum Monte Carlo (QMC) results for hydrogen phasesqmctest and for water.qmctestwater The vdW-DF method was first tested in non-selfconsistent formsryluladi00 ; rydberg03p126402 ; dion ; langreth05p599 ; beloschy13 using GGA calculations of the electron density as input for a post-processing evaluation of the nonlocal correlations. With a formal derivation of forces arising from the nonlocal correlation termthonhauser and with the introduction of efficient algorithms for computing the vdW-DF energy and forces,linearscaling ; roso09 the vdW-DF method today benefits from experience in widespread sparse- and general-matter applications.langrethjpcm2009 ; bearcoleluscthhy14 ; ROPPreview ; cooper08p1304 ; berland10p134705 ; londero11p054116 ; vdwsolids ; berland11p1800 ; ni12p424203 ; ihm12p424205 ; lazic12p424215 ; vdWReady ; Mperspectives ; poloni12p4957 ; poloni14p861 ; bjorkman14
The nonempirical vdW-DF-method is built around a semilocal internal (or inner) functional,lee10p081101 ; berland14p035412 ; bearcoleluscthhy14 ; ROPPreview with xc hole , and an evaluation of a nonlocal correlation energy . The internal functional keeps local exchange and correlation together but limits gradient corrections to exchange. The internal functional was introduced as a concept in Ref. [lee10p081101, ] but it underpins all formulations of the general vdW-DF versionsrydbergthesis ; dionthesis ; dion ; ROPPreview since it serves to model the local variation in the plasmon-pole response from which is formulated. There have to date only been brief discussions of this internal functional.rydbergthesis ; lee10p081101 ; bearcoleluscthhy14
The goal of this paper is to present the vdW-DF construction formally both in terms of the internal xc hole and physical pictures. Standard vdW-DF presentationsdion ; langreth05p599 ; thonhauser ; lee10p081101 ; bearcoleluscthhy14 start more directly with a plasmon-pole representation of the response. However, there are benefits of tracing the plasmon view back to a discussion of the associated internal functional xc holes. This makes it possible to discuss the close connection that exists between vdW-DF and the GGA descriptions. Moreover, the emphasis on the internal-functional response allows us to interpret the vdW-DF nonlocal correlation energy in terms of the RA physics picture of vdW forces.ra In turn these results allow us to illustrate the mechanisms by which vdW-DF retains a collectivity and nonlocality in its description of the screened response and materials interactions.
The paper is organized as follows. In section II, we present an xc-hole based formulation of the vdW-DF framework. It is meant to give the reader an overview of the vdW-DF method in a self-contained and alternative derivation cast in the concepts that we explore in this paper. In section III we plot and discuss the internal functional xc-hole components of vdW-DF. Section IV contains a demonstration of the link between the vdW-DF nonlocal correlation energy and the RA physics picture. Finally, section V summarizes the paper, while an appendix details that the vdW-DF rests on a correct longitudinal projection in its description of the electrodynamics coupling.
Ii The vdW-DF framework
To discuss the nature of vdW-DF, we start out by noting that constraint-based GGA has been enormously successful at describing chemically connected systems, both those that have zero dimensions (such as atoms, molecules, and nanoparticles) and those with one or more macroscopic dimensions (such as wires, surfaces, sheets, and solids).KieronPerspective ; BeckePerspective We denote such systems ‘molecular-type’ even if they can be infinitely extended.
Our discussion benefits from considering two or more such molecular-type regions in close vicinity. Fig. 1 shows an example with benzene adsorbed on graphene at a separation of 3 Å.behy13 The vdW-DF method seeks to extend the GGA success (for an individual fragment) by adding an account of the nonlocal correlations that arise among several such fragments as well as inside the fragments.langrethjpcm2009 ; bearcoleluscthhy14 ; vdwsolids Only the coupling mediated by the longitudinal component of the electrodynamical interaction described in the Coulomb gauge, is considered. The Coulomb Green function is .
The starting point is the exact adiabatic connection formulalape75 ; gulu76 ; lape77 (ACF):
(1) | |||||
where is the electron-electron coupling constant and ‘Tr’ denotes a full trace over the variationnotationfootnote in the reducible density-density correlation function at imaginary frequency . The infinite self energy is given by where is the density operator. The reducible density-density correlation function relates an external potential change , at a characteristic frequency , to resulting density changes, . The electron dynamics causes any such external potential to be screened; the system can also be described by a corresponding screened potential . The irreducible density-density correlation function relates this screened potential to the same density change, . A Dyson equation relates and , Ref. [lape77, ].
The vdW-DF framework expresses the xc energy
(2) |
in terms of an approximation for an effective longitudinal dielectric function . We make three observations. First, the expression (2) is formally equivalent to the ACF,lape75 ; gulu76 ; lape77 as given in Eq. (1), since the coupling-constant integration is captured in the definition of . Second, the explicit relation is given in terms of an effective external-potential (density) response function via
(3) |
where denotes the fluctuation or plasmon propagatorpinesnozieres ; lu67 at coupling constant . The corresponding longitudinal dielectric function is
(4) |
Third, writing
(5) |
makes a rigorously defined longitudinal projection of an effective dielectric function , as further detailed in appendix A. In vdW-DF it is assumed that a scalar, but nonlocal, dielectric function, can be used in Eq. (5); there exists a demonstration that such a scalar can be constructed for any given xc energy functional (for example, the exact functional).rydbergthesis ; ROPPreview
To specify the plasmon-pole description and the xc hole components we also introduce an effective screened (density) response function given by
(7) |
This definition complies with the Lindhard-type formulationlindhard
(8) |
The screened response Eq. (II) is specified via a longitudinal projection
(9) |
of the local-field dielectric response . From Eq. (8) it is clear that a scalar approximation for also specifies the vdW-DF dielectric functional, , that enters in Eq. (5) and determines the xc energy Eq. (2) in the vdW-DF framework.
The many-body response nature of any xc functional is naturally expressed in the ACF evaluation of the xc hole
(10) |
the electron-deficiency (at ) produced around an electron at point .
For a description of the density functional it is sufficient to work with the spherically averaged xc hole
(11) |
The local xc energy per particle is directly related to this xc hole
(12) | |||||
(13) | |||||
The exact relation
(14) | |||||
links to , and hence to an integral over the xc hole.
ii.1 The vdW-DF logic
A central idea in the vdW-DF framework is to exploit that the formally exact formulation Eq. (5) has already made one instance of the electrodynamics coupling explicit. One obtains truly nonlocal effects in the xc functional even when using a semilocal GGA-type functional to specify the details of the nonlocal form of . Accordingly, in the vdW-DF method we split the total xc energy functional and associated xc holes into semilocal and nonlocal contributions,rydbergthesis
(15) | |||||
(16) |
The first term of Eq. (15) is also called the outer semilocal functional and it is given by LDA correlation and a GGA description of gradient-corrected exchange. It would in principle provide an approximate description of a typical GGA problem (i.e., an individual of the molecular-type fragment shown in Fig. 1) because gradient-corrections to exchange are typically more important than gradient-corrected correlation.pecosabu06 ; PBEsolhole We shall, for ease of discussion, also refer to such a description as being of a GGA type. At the same time we note that should not be evaluated in isolation.
The second term is viewed as a perturbation,rydbergthesis ; ROPPreview capturing nonlocal correlation energy from the coupling of plasmon poles that characterize . The formulation of is, however, in practice based on the use of an internal semilocal functional that is similar to , but with an energy per particle that decreases more rapidly at large values of the scaled density . This choice is made to avoid spurious contributions emerging from low-density regions.ra ; anlalu96 ; dionthesis ; berland14p035412 ; bearcoleluscthhy14 ; ROPPreview The construction via a GGA-type (that also just contains LDA correlation plus GGA gradient-corrected exchange) allows vdW-DF to rest exclusively on formal diagrammatic inputperdew92p13244 ; rasolt ; lavo87 ; ra while avoidingthonhauser to explicitly formulate a gradient-corrected correlation term which is a necessary but also complex step in the GGA formulations.lape80 ; lameprl1981 ; lavo87 ; pebuwa96 ; pebuer96
The vdW-DF internal functional is given by a GGA-type internal-functional xc hole and it is used to introduce an approximate scalar dielectric function via
(17) |
This form is motivated by the observation that the longitudinal projection in Eq. (5) becomes redundant in the homogeneous electron gas (HEG) limit. As indicated by the subscript ‘grad’, the definition via Eq. (17) rests on an expectationrydbergthesis ; ryluladi00 ; langreth05p599 that this simplification holds approximately true for a weakly perturbed electron gas. Effectively, we write
(18) | |||||
(19) |
The vdW-DF dielectric function (18) is used, via Eq. (9), to also determine an approximation for the full dielectric function and hence extend the account to also include nonlocal correlations. Eq. (19) is a GGA-guided ansatz for (and hence ) that describes the effective full coupling-constant integration and screening effects within vdW-DF. We note that coincides to linear order with the related approximation that was used in the vdW-DF method presentation by Dion et al.dion Refs. [rydbergthesis, ; dion, ; ROPPreview, ] suggest explicit forms for (and ), given in terms of a model plasmon dispersion at two coordinate points.
Using , we interpret the poles as an approximative specification of the collective modes of the system described by the internal functional Eq. (17), i.e., the zeros of . In the HEG limit the plasmon-pole dispersion (entering ) is the same everywhere and renders a direct specification of ; in the presence of gradients, the spatial and momentum variation in the plasmon poles, , represent instead only an approximative specification of the set of internal-functional collective modes . In any case, these plasmon modes are a direct reflection of the shape of the semilocal internal-functional xc hole, as explained in Sec. II.2.
The general-geometry vdW-DF versionsdion ; lee10p081101 ; berland14p035412 approximate the xc energy
(20) | |||||
(21) |
where the vdW-DF nonlocal correlation term
(22) | |||||
is evaluated by expanding both terms in the same plasmon-pole description :
(23) |
This quadratic expansion for has the same appearance whether cast in (as done in Refs. [rydbergthesis, ; ROPPreview, ]) or in terms of (as done in Ref. [dion, ; langrethjpcm2009, ]) because these agree to lowest order.
For given choices of the internal-functional form (and hence of plasmon poles in , Ref. rydbergthesis, ; ROPPreview, ) and of the vdW-DF form generally discards a cross-over term
(24) |
An improved alignment between and minimizes the difference between and an evaluation based on the formal ACF recast Eq. (2). Such an alignment reflects consistencyberland14p035412 between the plasmon response of the internal functional and that which characterizes , and it is beneficial because it allows the longitudinal projection [in Eqs. (2) & (5)] to leverage an automatic conservation of the full xc hole.berland14p035412 ; bearcoleluscthhy14
ii.2 The vdW-DF internal functional specification
The internal functional is semilocal and of a GGA type. It is specified by LDA exchange energy per particle, , where denotes the local Fermi wavevector, and an enhancement factor,
(25) |
The internal functional is thus fully given by the local value of the density and of the scaled density gradients, . In the vdW-DF design,dion ; thonhauser the internal functional is exclusively given by the LDA-correlation term (independent of ) and an exchange gradient enhancement (independent of ):
(26) |
In the vdW-DF1dion ; dionerratum and vdW-DF-cxbehy13 versions we stick with the Langreth-Vosko analysis for screened exchange,lavo90 giving
(27) |
specified by . In vdW-DF2,lee10p081101 formal scaling analysisschwinger ; elliott09p1485 for pure exchange yields an enhancement of curvature with . The form of is taken from Ref. [perdew92p13244, ].
The resulting energy-per-particle expression (25) provides a full specification of a vdW-DF internal xc hole inside a model that assumes a Gaussian spherical average formrydbergthesis ; ROPPreview
(28) |
The simple form enables analytical evaluation for many of the spatial integrations in the resulting description of . Also, the model form Eq. (28) ensures that is itself conserved,
(29) |
for all exchange-enhancement choices in and for any value of . Choosing in the Gaussian model Eq. (28) provides a simple relation between the inverse length scaledion
(30) |
of the model hole and the internal functional energy-per-particle variation Eq. (25). This variation in is in turn used to formulate the vdW-DF evaluation of in terms of a universal kernel , as detailed in Refs. [rydbergthesis, ; dion, ; dionerratum, ; thonhauser, ; ROPPreview, ].
An important point for our discussion and interpretation, Sec. IV, is that the shape of the internal semilocal xc hole, given by Eqs. (28) & (30), is used in vdW-DF to determine the local variation in the plasmon poles , Refs. [rydbergthesis, ; dion, ; ROPPreview, ]. The connection is made by noting that the spherical averaged xc hole also defines a natural wavevector decompositionlape75 ; lape77 ; lape80 ; lavo87 ; lavo90 for the internal functional energy per particle,
(31) | |||||
(32) |
Evaluating the imaginary frequency integral in the formal relation
(33) |
with the plasmon-pole specificationrydbergthesis ; ROPPreview for yields a wavevector decompositionlape77
(34) |
that links to the chosen description of the internal functional xc hole .
Iii Internal-functional exchange-correlation holes
In this section, we visualize the vdW-DF internal xc hole and compare it to that of the numerical-GGA constructions.pewa86 ; pebuwa96 This casts light on the nature of the vdW-DF since these internal xc holes define the vdW-DF dielectric function from which vdW-DF builds an account of truly nonlocal correlations.
Figure 2 shows the scaled-density gradient contours of a benzene dimer at binding separation — a typical molecular binding system.berland14p035412 The plot documents that a region with low-to-moderate values exists between the molecular fragments, where the over-lap of two decaying densities causes a saddle-point or trough-like behavior. Non-local correlation contributions arise from low density regionsbehy13 and from low-to-moderate values of the scaled density gradient .berland14p035412 There is no contradiction even if typically enhances exponentially outside a molecular-type region, because in vdW-DF, the binding from nonlocal correlations arises predominantly in the regions between molecular fragments.behy13 ; lazic12p424215 ; caciuc12p424214
Fig. 3 shows the internal functional enhancement factor, Eq. (26) of vdW-DF1 and vdW-DF2. For vdW-DF1, this figure corresponds to density gradient that are often most relevant for binding contributions.berland14p035412
The spherically averaged real-space xc hole of the internal functional of vdW-DF is extracted from an inverse Fourier transform of Eq. (28). This real-space internal xc hole is given by the following simple Gaussian form:
(35) | |||||
(36) |
We choose to discuss the role of the internal xc hole in a form scaled with a distance-weighted measure:gulu76
(37) |
The weighted expression (37) reflects how the shape of the xc hole and the electrodynamics coupling determine the energy-per-particle variation in the corresponding internal semilocal functional, as given by Eq. (13).
Fig. 4 shows a contour plot of the internal xc hole as defined by Eq. (26) and weighted and scaled according to Eq. (37). The plot represents the behavior of both vdW-DF1/vdW-DF-cx and vdW-DF2, since the vertical axis is simply the value of the internal functional enhancement factor, . The horizontal axis represents the scaled distance from the electron position (hole center).
Fig. 5 compares the internal functional xc hole against the Perdew-Wang (PW) xc hole model for the homogeneous electron gas.PWhole92 This model captures the salient non-oscillatory features of the correlation hole and therefore compares well with QMC calculations. In the bottom panel the non-oscillatory approximation for the exchange hole is used.PWhole92 The holes are plotted as functions of the scaled distances and weighted by the radial measure as earlier. The comparison is shown for two values of the Wigner-Seitz radius . The value of corresponds to the density between two neighboring C atoms of a benzene molecule, while the value of corresponds to the density 1.5 Å out of the benzene-plane above these two C atoms.
The PW model agrees well with the internal functional xc hole. In particular, the agreement is strikingly similar to the xc hole relying a non-oscillatory approximation for the exchange hole, with the exception of a slightly different trend with changing . In summary, vdW-DF not only keeps a good balance between local exchange and correlation contributions, in line with the DFT tradition,lape77 but also an internal xc hole form in fair agreement with QMC.
Equally interesting is the question if the vdW-DF internal-hole characterization also remains useful when applied to typical systems. These have density gradients and we need to selectively add the effects of gradient-corrected exchange in (by the vdW-DF design logic). The shape of the vdW-DF model internal xc hole should remain reasonable at values of the scaled gradient that are deemed relevant for the evaluation of the nonlocal correlation.
Numerical GGApewa86 ; pebuwa96 ; PBEsolhole is a well defined procedure to impose charge conservation and a negativity condition on the xc hole of a gradient expansionrasolt ; lape80 ; lameprl1981 ; lape82 ; RasoltDisc1 ; RasoltDisc2 ; lavo87 around a homogeneous electron gas. Since this xc hole construction can be used to derive popular xc functionals such as PW86, PBE, and PBEsol,pewa86 ; pebuer96 ; perucs08 we compare relevant components (including gradient corrected exchange but not gradient-corrected correlations) of also these to those of the vdW-DF internal functional at relevant nonzero values of the density gradient.
Figure 6 compares the vdW-DF2 internal functional xc-hole representations at and with the numerical-GGA specification of Ref. [pebuwa96, ]. For the correlation, only the LDA part of the numerical xc hole is included, since only this component is used in the construction of the internal xc-hole of vdW-DF. The vdW-DF2 version is chosen for comparison, because its internal functional has the same small- behavior as the PBE functionalpebuer96 and the numerical hole construction of Ref. [pebuwa96, ] leads to a GGA xc that resembles PBE. Conservation is built into the vdW-DF model of the internal xc hole (for any choice of ,) Eq. (29) and there is no need to enforce hard cut-offs as in the numerical GGA construction.ROPPreview
Comparing the two panels of Fig. 6 with those of Fig. 5, we see that in all cases the holes become deeper and shorter-ranged as increases. However, whereas the two holes agree fairly well for small , their shapes grow dissimilar as increases. The agreement with the numerical-GGA construction is best when (as is more relevant for larger and flatter fragments close to binding separationberland14p035412 ) we can limit the value of the scaled gradient to .
Iv Interpretation of the vdW-DF nonlocal correlation energy
In this section, we show that the vdW-DF method relies on a nonlocal-correlation formulation that can be interpreted as an implementation of the RA picture of vdW forces.ra As mentioned in the introduction, this picture sees nonlocal correlations as arising from an electrodynamical coupling of, originally independent, xc holes of a traditional semilocal functional description; in vdW-DF represented by the GGA-type internal xc functional holes . The vdW-DF nonlocal correlation term represents a counting of coupling-induced zero-point energy shifts of the characteristic plasmon modes of these internal-functional xc holes.ra ; anlalu96 ; ROPPreview We provide the interpretation by adapting the analysis that Mahan used to discuss the nature of vdW forces and detail their relation to Casimir forces.capo48 ; jerry65
To begin, we formulate the exact xc energy (2)
(38) | |||||
(39) |
where time-reversal symmetry has been used to extend the integration over to . The right-hand side of Eq. (38) can be evaluated from a contour around the complex-frequency plane that runs down the imaginary axis and closes around the half-plane of positive frequencies. This contour picks up characteristic poles of with simple residues, as is evident in the rewritejerry65
(40) |
The form (40) counts the sum of collective modes, given by . With a general specification of (beyond the approximation used in the vdW-DF versions) there will also be corrections from poles in . This second set of poles corresponds to singularities in the local-field responsejerry65 . Such singularities are normally associated with particle-hole excitations.SawadaPlas ; RasoltDisc1 ; jerry65
The nonlocal correlation term Eq. (22) is formulated as a difference between the xc energy when defined in terms of and . It therefore expresses an xc energy shift produced by an electrodynamical coupling. The general result Eq. (38) allows us to discuss the nature of this coupling-induced energy shift.
We first consider a single molecular-type fragment, i.e., one of the molecular-type regions in Fig. 1. We use and to denote the collective modes of the internal and of the full vdW-DF xc functional, given respectively by and . The contour evaluation Eq. (40) provides the formal evaluation,
(41) |
since, as motivated below, we can ignore contributions from particle-hole excitations.lape77
The vdW-DF idea of expressing all functional components through an analysis of the response in an internal (GGA-type) functional is what makes Eq. (41) relevant for analyzing . A key observation is that the specification of the vdW-DF internal functional keeps local exchange and local correlation together to allow a cancellation of terms arising from particle-hole excitations.lape77 Accordingly, Ref. [dion, ] uses simply a plasmon-pole representation for but leaves no room for singularities directly in . The same observation underpins our assumption, in Sec. II, that all singularities in the internal-functional specification,rydbergthesis ; ROPPreview , should be seen as exclusively reflecting collective (plasmon) poles of . Moreover, while Eq. (4) ensures that the reducible response function has singularities at the collective modes of , there can be no single-particle singularities in the irreducible response . This follows in the vdW-DF framework (and only there) because is set by the internal functional behavior, through and Eqs. (18) & (9).
The formal evaluation in Eq. (41) can be used to interpret the vdW-DF nonlocal correlation term, Eq. (22), as an implementation of the RA picture of vdW forces.ra The term tracks changes in characteristic plasmon modes of the system as described in vdW-DF and in a semilocal functional defined by xc functional holes (with a partial GGA character). Inclusion of the electrodynamical coupling changes the dielectric functions and hence the characteristic plasmon modes that characterize . In effect, the energy shift in Eq. (41) tracks the effects of coupling the GGA-type internal functional holes.
It is also interesting to compare the formal framework of the vdW-DF method and the random phase approximation (RPA).SawadaPlas ; Pines ; lape77 ; Furche08 ; Eshuis12 ; Ren12p7447 The RPA correlation energy islape77
(42) |
Our comparison will be based on the full representation (22) that underpins Eq. (41) and was used in an early seamless functional for layered structures.rydberg03p126402 We shall in Sec. V return to a discussion of RPA and vdW-DF, keeping in mind also that the recent vdW-DF versionsdion ; lee10p081101 ; berland14p035412 use a second-order expansion of Eq. (22) in . There are formal similarities between RPA and the vdW-DF method. The vdW-DF framework, Eqs. (2) & (5), builds on the exact ACF as does RPA; in fact, the RPA xc energy is obtained by inserting the approximation in Eq. (2). Also, the RPA correlation energy can be exactly reformulated,SawadaPlas ; Furche08
(43) |
where () denotes a RPA excitation energy as described to lowest (full) order in . The RPA interpretation as a counting of zero-point energy shifts, Eq. (43), resembles the interpretation (41) that we present for the nonlocal-correlation energy in the vdW-DF method. There are also fundamental differences. The RPA crafts from particle-hole excitations, typically given by Kohn-Sham orbitals and energies,Furche08 ; Eshuis12 ; Ren12p7447 whereas vdW-DF proceeds by asserting its response description through a plasmon model.rydberg03p126402 ; dion ; thonhauser The summation in the vdW-DF interpretation (41) is restricted to zero-point energy contributions defined by collective modes.
Additional details about the nature of the vdW-DF nonlocal correlation term can be obtained by considering the case of two molecular fragments, and , separated by a delineation surface as illustrated in Fig. 1. In the following we focus exclusively on the binding that arises in the nonlocal-correlation component of , noting that there will also be other interaction components (arising through the interplay between kinetic-energy repulsion, Coulomb terms, and in the outer semilocal functional ). The analysis applies also when, as in weak chemisorption, there is some density overlap, as we can proceed within a superposition-of-density scheme.harris85p1770 ; beloschy13 The analysis is not relevant for cases where there are also chemical bonds across the delineation surface but we defer a discussion of limitations until we can formulate this in terms of criteria on the electron-response description.
To make the discussion more specific, we cast the interpretation in terms of explicit approximations. We let (and ) denote the density of fragment when treated in isolation. These densities should be seen as DFT solutions as obtained in a vdW-DF version; note that extends into the area that the delineation surfaces assign as region and vice versa. We assume that the multicomponent density can be sufficiently approximated as a sum of fragment densities, . This is a general approximation schemeharris85p1770 which is often accurate for systems held together by dispersive forces in competition with other interactions.beloschy13 From separate densities and we can define per-fragment screened response functions and . We note that the corresponding reducible response function must have singularities at the vdW-DF collective modes for fragment . The same goes for the description of fragment . Next we introduce as the region-projected part of this irreducible response i.e., the matrix formed from by restricting both coordinates to reside in delineated region ‘’ as well as corresponding projections for the reducible response function and for .
At this stage we can discuss the limitations on the extended analysis presented below. One requirement is that we approximately retain a Dyson-like link, as in Eq. (II), among the response descriptions even when working with the fragment-projected response description:
(44) |
A second, related, requirement is that the collective modes also represent the poles of . This second requirement can be formulated as the condition that (where the determinant reflects an integration over the entire space) has the same zeros as is found for (where the determinant range is limited to the delineated region). The conditions can only hold approximately except when discussing well-separated fragments.
Notwithstanding the requirements for using a partitioning scheme, we proceed to deepen our analysis of the nonlocal correlation term. Such a scheme has also been used, for example, to extract an asymptotically exact evaluation of interactions among defects on a surface supporting a metallic surface state.lako78 ; ei78 ; hype00 ; focus00 ; hyei02trio ; berland09p155431 ; hy12p424219 The important part of the Coulomb coupling is in this problem the component of the Coulomb term that connects a point in the delineated region ‘’ with a point in the other region ‘’. Using a simple matrix factorization of (and of ) we thus obtain
(45) | |||||
(46) |
We note in passing that the result of Eq. (45) is consistent with the traditional result for the asymptotic vdW bindinganlalu96 ; hurylula99
(47) |
where denotes a dipole-dipole coupling tensor between points in separate regions and where denotes the external-field susceptibility of fragment (or fragment ). The connection between Eq. (45) and Eq. (47) is made by expanding the logarithm and noting that
(48) |
specifies the vdW-DF approximation for these susceptibilities.rydbergthesis
For a discussion of the electrodynamical coupling expressed in we provide a contour-integral evaluation of Eq. (45). Using the contour-integration formulation (40) for we now have contributions from the poles of the fragment response function . We assume that these give rise to contributions that resemble those specified by the exciton-susceptibility tensors in Ref. jerry65, . Equally important, the form Eq. (45) has regular plasmon-pole contributions given by the zeros of
(49) |
Adapting the argument presented in Ref. jerry65, , these inter-fragment collective mode correspond to a coupling between polarizability contributions defined in and by modes and . Overall the coupling contour integration leads to an approximative evaluation
(50) |
and establishes a further link between the nonlocal-correlation term in vdW-DF and the RA picture, viewing vdW forces as arising as a coupling of (semilocal, initially independent) xc holes.ra ; ma ; anlalu96
Finally, we note that we can extend a partition-based analysis also to the case when there are three (or more) molecular fragments, again adopting the analysis used for the study of electronic substrate-mediate interactions among defects on surfaces.ei78 ; hyei02trio For cases with three molecular-type fragments, denoted , , and , we find
(51) | |||||
The trio term does not naturally enter in the vdW-DF description when investigating a system with only two molecular-type regions near binding separation. The argument for approximating as exclusively connecting two points inside the same fragment breaks down if one were to partition an individual molecular-type region.
V Discussion
The influence of screening and nonadditivity effects on the vdW forces are explored in a significant body of literature, for example in Refs. [ra, ; Axilrod, ; BarashI, ; BarashII, ; SerneliusI, ; SerneliusII, ; Dobson06, ; Tkatchenko08, ; Lebegue10, ; Tkatchenko12, ; Ruzsinszky12, ; PT, ; DoGoDisp12, ; vdWscale, ; DobsonMB14, ; DobsonPerspective, ]. In this section, we discuss to what extent the vdW-DF methoddion ; lee10p081101 ; berland14p035412 can capture such effects. In particular, we first compare vdW-DF to the RPA for the correlation energy and then discuss vdW-DF in light of a recently suggested classification scheme of dispersion interaction effects that lies beyond pair-wise summations.DobsonPerspective
The vdW-DF method shares with RPA an electron-based foundation as they avoid partitioning into, for example, atomic components. The methods also share a zero-point-energy counting nature (Sec. IV), an emphasis on approximating the ACF through longitudinal dielectrical functions that comply with the continuity equation (App. A), and conservation of the associated xc hole.dion ; berland14p035412 ; bearcoleluscthhy14 ; ROPPreview One difference is that the vdW-DF is based on a fully screened response description via a plasmon-based starting point that reflects a GGA-type internal-functional xc hole, whereas in RPA one starts with independent-particle excitations.Furche08 ; Eshuis12 ; Ren12p7447 With the full expression (22), used in an early seamless vdW-DF functional,ryluladi00 ; rydberg03p126402 the vdW-DF method relies on the same machinery as RPA for systematically including screening effects, namely the Dyson equation for the density-density correlation function,lape77 as shown in Sec. II.
At the same time, the popular general-geometry vdW-DF versions and closely related variantsdion ; lee10p081101 ; cooper10p161104 ; vdwsolids ; berland14p035412 ; hamada14 rest on a second-order expansion (23) of , a step that is not used in RPA calculations. An interesting question is then how much of the screening, collectivity, and nonadditivity effects are retained after making this truncation. The question is complex and we limit the discussion to making some comments in the context of a recent perspective article by Dobson.DobsonPerspective
Dobson classifies nonadditivity effects as follows: Class ‘A’ contains the effects of bond formation. These effects are automatically included in vdW-DF. Class ‘B’ is the spectator effect; that is, the modification by an additional molecular-type fragment on the electrodynamics coupling between two molecular-type fragments. Class ‘C’ contains many-body effects that result with nondegenerate electron states and their ability to enhance the electronic response. A consequence is, for example, different asymptotic scaling laws for the vdW attraction between sheets of metals or among metallic nanotubes than between insulators.BarashI ; BarashII ; SerneliusI ; SerneliusII ; Dobson06 ; Lebegue10 ; DoGoDisp12 ; DobsonMB14 ; DobsonPerspective Class ‘B’ and ‘C’ are expected to be of greater importance for asymptotic interactions than at binding separations where there are contributions from many plasmons.Lebegue10 ; behy13 ; DobsonMB14
For a vdW-DF version to fully address nonadditivity effects of class ’B’, it requires that the evaluation proceeds with the full interaction form,ryluladi00 ; rydberg03p126402 not the expansion (23) used for the more recent and popular vdW-DF versions.dion ; lee10p081101 ; berland14p035412 However, vdW-DF reflects multipole enhancements in the binding among moleculesberland10p134705 ; berland11p1800 and image-plane formation in the binding of carbon nanotubeskleis08p205422 and in challenging physisorption problems.lee11p193408 ; lee12p424213 Image-plane effects are captured in those recent expanded vdW-DF versions through the stronger sensitivity to the low-density regions arising at surfaces than to the high-density regions of the bulk.behy13
To fully capture effects in class ’C’ one would also need to refine the vdW-DF inner-functional response model beyond a simple plasmon model relying on a GGA-based account. However, some of the energetic impact of these effects is also, in practice and at a cruder level, reflected in the modern vdW-DF versions: In low-density, highly homogeneous systems, typical of a metal surface,lee11p193408 ; lee12p424213 the vdW-DF plasmon model yields small excitation energies, strongly enhancing non-local correlation effects. On the other hand, except at edge regions, the GGA-based construction of vdW-DF does not distinguish between a limited molecular-type fragment, such as the center part of a polyaromatic hydrocarbon, and an extended fragment that has no gap, such as graphene or a metallic nanotube. This is a distinction that becomes important at asymptotic separations between fragments.Dobson06 ; Lebegue10 ; DoGoDisp12
The recent vdW-DF versionsdion ; lee10p081101 ; berland14p035412 are formulated with the expectation that the second-order expansion (23) is often sufficient in binding situations, with two molecular-type fragments in close proximity.discussHollow At binding it is important to treat truly nonlocal correlation effects and the more local/semi-local correlation effects on a same footing.rydbergthesis ; dion ; thonhauser ; behy13 In the sections above we have illustrated that vdW-DF should not be viewed merely as a summation of contributions from pairs of density points. It is rather an expression of coupling of semilocal xc holes with a finite extension and with a shape and dynamics that already reflect a GGA-type response behavior. The internal functional xc holes express a collectivity that generally extends beyond that of a single atom and represents a GGA-level of screening that we assume is often adequate for treating interfragment binding.
The good performance of vdW-DF, in particular for the most recent nonempirical versions and variants, indicates that the vdW-DF method is capable of accurately reflecting the complicated balance that can exist between general interaction contributions.lee11p193408 ; lee12p424213 ; behy13 ; bearcoleluscthhy14 For example, the recent consistent-exchange vdW-DF-cx versionberland14p035412 can correctly describe the competition between covalent and ionic bonds in ferroelectrics, and between exchange effects and ionic and vdW attraction in weak-organic chemisorption.bearcoleluscthhy14 One needs Axilrod-Tellerra ; Axilrod corrections (51) and beyondryluladi00 ; rydberg03p126402 ; Furche08 ; Tkatchenko12 to fully characterize the general dispersive interaction in systems that have three or more molecular-type regions. However, when two molecular-type regions are at their binding separation there is not generally room for a third molecule to get close and significantly influence that coupling.
Finally, for a quantitative discussion of what screening effects are retained in recent vdW-DF versions one needs to compare the results of these versions with those obtained when using the full formryluladi00 ; rydberg03p126402 under the same approximation for the internal-functional description. Two of us have led one early such explorationSvetlaCompare but no conclusion can be reached in that study because we refined the plasmon model between the layered-geometry formulation of Ref. [rydberg03p126402, ] and the launching of vdW-DF1.dion A comparison of the results based on the modern response form is beyond the present scope.
Vi Summary
Several formal properties of the vdW-DF theory have been highlighted. Specifically, we have documented how an effective internal xc hole can be viewed as a central building block in obtaining the nonlocal correlation of vdW-DF. We have documented how this internal xc hole resembles the xc hole construction that underpins standard GGA descriptions. Further, we have argued how the nonlocal correlations in vdW-DF can be interpreted as arising from the shift in collective modes induced by the electrodynamical coupling between such xc holes. This argument connects vdW-DF to the well-established RA picture of vdW-DF interactions. Finally we have compared the vdW-DF method to RPA.
By discussing the formal properties of vdW-DF and links to other theories, we hope to help build bridges that stimulate the dissemination of ideas, not just within the field of van der Waals interactions but also within the wider field of material modeling.
Acknowledgements.
The authors are grateful for many insightful discussions with David C. Langreth, who passed away in 2011, and with Bengt I. Lundqvist. The authors are furthermore grateful for input, insight, and encouragement from Prof. G. D. Mahan, primarily during his two extended visits to Chalmers, in 2009 and 2011. This work was supported by the Swedish research council (VR) and the Chalmers Materials Area of Advance theory initiative.Appendix A Role of the continuity equation in vdW-DF
This appendix details the correct longitudinal projection in electrodynamics and is based on notes by and discussions with D.C. Langreth. It serves to further motivate writing the xc energy
(52) |
so that it expresses the exact longitudinal projection. The form (52) reflects the continuity equation as well as the constituent equations of the electrodynamical response in materials.
Consider Ohm’s law for the current induced by a local field ,
(53) |
and corresponding to the induced charges . In Eq. (53), we use to denote the nonlocal conductivity tensor. In turn, this tensor corresponds to a nonlocal dielectric function
(54) |
Fourier transforming gives , and a continuity specification
(55) |
that relates the local field and the induced charge
(56) |
We infer an exact, general microscopic relation between the (longitudinal) local-field response and the dielectric tensor
(57) | |||||
(58) |
The result Eq. (57) reflects the continuity equation and identifies
(59) |
as the proper (consistent) definition of the longitudinal projection of the dielectric repose in an inhomogeneous system.
In Eq. (58) the difference takes the form of a local-field susceptibility tensor . Using to denote the corresponding external-field susceptibility we also have a relation for the external-field response
(60) |
Finally, the correct longitudinal projection of the dielectric tensor is given by
(61) |
This is demonstrated by expressing the microscopic relation in wavevector space
(62) | |||||
where we have used the specification Eq. (59). Solving for gives
(63) |
References
- (1) K. Rapcewicz and N. W. Ashcroft, Phys. Rev. B 44, 4032(R) (1991).
- (2) O. Gunnarsson and B. I. Lundqvist, Phys. Rev. B 13, 4274 (1976).
- (3) J. Harris and R. O. Jones, J. Phys. F.: Metal Phys. 4, 1170 (1974).
- (4) O. Gunnarsson, M. Jonson, and B. I. Lundqvist, Phys. Rev. B 20, 3136 (1979).
- (5) J. P. Perdew and D. C. Langreth, Sol. State. Commun. 17, 1425 (1975).
- (6) D. C. Langreth and J. P. Perdew, Phys. Rev. B 15, 2884 (1977).
- (7) Y. Wang and J. P. Perdew, Phys. Rev. B 44, 13298 (1991).
- (8) R. Eisenshitz and F. London, Z. Phys. 60, 491 (1930).
- (9) F. London, Z. Phys. 63, 245 (1930).
- (10) G. D. Mahan, J. Chem. Phys. 43, 1569 (1965).
- (11) A. C. Maggs and N. W. Ashcroft, Phys. Rev. Lett. 59, 113 (1987).
- (12) M. Rasolt and D. J. W. Geldart, Phys. Rev. Lett. 35, 1234 (1975); Phys. Rev. B 13, 1477 (1976).
- (13) D. C. Langreth and J. P. Perdew, Phys. Rev. B 21, 5469 (1980).
- (14) M. Rasolt and D. J. W. Geldart, Phys. Rev. B. 25, 5133 (1982).
- (15) C. D. Hu and D. C. Langreth, Phys. Rev. B 33, 943 (1986).
- (16) M. Rasolt and D. J. W. Geldart, Phys. Rev. B. 34, 1325 (1986).
- (17) D. C. Langreth and S. H. Vosko, Phys. Rev. Lett. 59, 497 (1987).
- (18) D. C. Langreth and S. H. Vosko, Adv. Quant. Chem. 21, 175 (1990).
- (19) D. C. Langreth and M. J. Mehl, Phys. Rev. Lett. 47, 446 (1981).
- (20) D. C. Langreth and J. P. Perdew, Phys. Rev. B 26, 2810 (1982).
- (21) D. C. Langreth and M. J. Mehl, Phys. Rev. B 28, 1809 (1983).
- (22) K. Burke, J. Chem. Phys. 136, 150901 (2012).
- (23) A. D. Becke, J. Chem. Phys. 140, 18A301 (2014).
- (24) Y. Andersson, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 76, 102 (1996).
- (25) E. Hult, H. Rydberg, B. I. Lundqvist, and D. C. Langreth, Phys. Rev. B 59, 4708 (1999).
- (26) E. Hult, Y. Andersson, B. I. Lundqvist, and D. C. Langreth, Phys. Rev. Lett. 77, 2029 (1996).
- (27) H. Rydberg, B. I. Lundqvist, D. C. Langreth, and M. Dion, Phys. Rev. B 62, 6997 (2000).
- (28) H. Rydberg, Nonlocal Correlations in Density Functional Theory, Ph.D. thesis, Chalmers University of Technology & Göteborg University, Sweden 2001.
- (29) H. Rydberg, M. Dion, N. Jacobson, E. Schröder, P. Hyldgaard, S. I. Simak, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 91, 126402 (2003).
- (30) M. Dion, van der Waals forces in density functional theory, Ph.D. thesis, Rutgers University, New Jersey, U.S.A. 2004.
- (31) M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004).
- (32) M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 95, 109902(E) (2005);
- (33) D.C. Langreth, M. Dion, H. Rydberg, E. Schröder, P. Hyldgaard, and B.I. Lundqvist, Int. J. Quantum Chem. 101, 599 (2005).
- (34) T. Thonhauser, V. R. Cooper, S. Li, A. Puzder, P. Hyldgaard, and D. C. Langreth, Phys. Rev. B 76, 125112 (2007).
- (35) V. R. Cooper, Phys. Rev. B 81, 161104(R) (2010).
- (36) K. Lee, É. D. Murray, L. Kong, B.I. Lundqvist, and D.C. Langreth, Phys. Rev. B 82, 081101(R) (2010).
- (37) K. Berland and P. Hyldgaard, Phys. Rev. B 87, 205421 (2013).
- (38) K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014).
- (39) K. Berland, C.A. Arter, V.R. Cooper, K. Lee, B.I. Lundqvist, E. Schröder, T. Thonhauser, and P. Hyldgaard, J. Chem. Phys. 140, 18A539 (2014).
- (40) K. Berland, V.R. Cooper, K. Lee, E. Schröder, T. Thonhauser, P. Hyldgaard, and B.I. Lundqvist, submitted to Reports on Progress in Physics, invited review, 2014.
- (41) D.C. Langreth, et al., J.Phys.:Condens. Matter, 21, 084203 (2009).
- (42) J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992).
- (43) J. Schwinger, Phys. Rev. A 22, 1827 (1980); Phys. Rev. A 24, 2353 (1981).
- (44) P. Elliott and K. Burke, Can. J. Chem. 87, 1485 (2009).
- (45) J. P. Perdew and Y. Wang, Phys. Rev. B 33, 8800 (1986).
- (46) J.P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
- (47) Y. Zhang and W. Yang, Phys. Rev. Lett. 80, 890 (1998).
- (48) É. D. Murray, K. Lee, and D. C. Langreth, J. Chem. Theory Comput. 5, 2754 (2009).
- (49) J. Klimeš, D. R. Bowler, and A. Michaelides, Phys. Rev. B 83, 195131 (2011).
- (50) I. Hamada, Phys. Rev. B 89, 121103 (2014).
- (51) R. C. Clay, J. McMinis, J. M. McMahon, C. Pierleoni, D. M. Ceperley, M. A. Morales, Phys. Rev. B 89, 184106 (2014).
- (52) M.A. Morales, J.R. Gergely, J. McMinis, J. M. McMahon, J. Kim, and D. M. Ceperley, J. Chem. Theory Comput. 10, 2355 (2014).
- (53) K. Berland, E. Londero, E. Schröder, and P. Hyldgaard, Phys. Rev. B 88, 045431 (2013).
- (54) A. Gulans, M. J. Puska, and R. M. Nieminen, Phys. Rev. B 79, 201105 (2009).
- (55) G. Roman-Perez and J. M. Soler, Phys. Rev. Lett. 103, 096102 (2009).
- (56) V. R. Cooper, T. Thonhauser, A. Puzder, E. Schröder, B. I. Lundqvist, and D. C. Langreth, J. Amer. Chem. Soc. 130, 1304 (2008).
- (57) K. Berland and P. Hyldgaard, J. Chem. Phys. 132, 134705 (2010).
- (58) E. Londero and E. Schröder, Phys. Rev. B 82, 054116 (2010).
- (59) K. Berland, Ø. Borck, and P. Hyldgaard, Comp. Phys. Commun. 182, 1800 (2011).
- (60) N. Nijem, P. Canepa, L. Kong, H. Wu, J. Li, T. Thonhauser and Y. J. Chabal, J. Phys.: Condens. Matter 24. 424203 (2012).
- (61) Y. Ihm, V. R. Cooper, L. Peng and J. R. Morris, J. Phys.: Condens. Matter 24. 424205 (2012).
- (62) P. Lazić, N. Atodiresei, R. Brako, B. Gumhalter, and S. Blügel, J. Phys.:Condens. Matt. 24, 424215 (2012).
- (63) T. Björkman, A. Gulans, A. V. Krasheninnikov, and R. M. Nieminen, J. Phys.: Cond. Matter 24, 424218 (2012).
- (64) J. Klimes and A. Michaelides, J. Chem. Phys. 137, 120901 (2012).
- (65) R. Poloni, B. Smit, and J. B. Neaton, J. Phys. Chem. A 116, 4957 (2012).
- (66) R. Poloni, K. Lee, R. F. Berger, B. Smit, and J. B. Neaton, J. Phys. Chem. Lett. 5, 861 (2014).
- (67) T. Björkman, “Testing several recent van der Waals density functionals for layered structures,” http://arXiv.org/abs/1406.2436.
- (68) All response functions and dielectric functions are nonlocal, i.e., functions of two spatial coordinates. The dependence is generally suppressed by considering equivalent operators . The Fourier transform is defined as
- (69) P. Nozières and D. Pines, Phys. Rev. 111, 442 (1958).
- (70) B. I. Lundqvist, Phys. Kond. Mat. 6, 193 (1967).
- (71) J. Lindhard, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 28, (1954).
- (72) J. P. Perdew, L. A. Constantin, E. Sagvolden, and K. Burke, Phys. Rev. Lett. 97, 223002 (2006).
- (73) L. A. Constantin, J. P. Perdew, and J. M. Pitarke, Phys. Rev. B 79, 075126 (2009).
- (74) J. P. Perdew, K. Burke, and Y. Wang, Phys. Rev. B 54, 16533 (1996).
- (75) V. Caciuc, N. Atodiresei, M. Callsen, P. Lazić, and S. Blügel, J. Phys.:Condens. Matt. 24, 424214 (2012).
- (76) J. P. Perdew and Y. Wang, Phys. Rev. B 46, 12947 (1992).
- (77) J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Phys. Rev. Lett. 100, 136406 (2008).
- (78) H. B. G. Casimir and D. Polder, Phys. Rev. 73, 360 (1948).
- (79) K. Sawada, Phys. Rev. 106, 372 (1957); K. Sawada, K. A. Brueckner, N. Fukuda, and R. Brout, Phys. Rev. 108, 507 (1957).
- (80) For an introduction to RPA calculations of the HEG ground-state energy, see David Pines, “Elementary excitations in Solids” (Addison-Wesley Publishing Company, Inc, Redwood City, California, 1963), Appendix C.
- (81) F. Furche, J. Chem. Phys. 129, 114105 (2008).
- (82) H. Eshuis, J.E. Bates, and F. Furche, Theor. Chem. Acc. 131, 1084 (2012).
- (83) X. Ren, P. Rinke, C. Joas, and M. Scheffler, J. Mater. Sci. 47, 7447 (2012).
- (84) J. Harris, Phys. Rev. B 31, 1770 (1985).
- (85) K. H. Lau and W. Kohn, Surf. Sci. 75, 69 (1978).
- (86) T. L. Einstein, Surf. Sci. 75, 161L (1978).
- (87) P. Hyldgaard and M. Persson, J. Phys.:Condens. Matt. 12, L13 (2000).
- (88) J. Repp, F. Moresco, G. Meyer, K.-H. Rieder, P. Hyldgaard and M. Persson, Phys. Rev. Lett. 85, 2981 (2000).
- (89) P. Hyldgaard and T. L. Einstein, Eur. Phys. Lett. 59, 265 (2002).
- (90) K. Berland, T. L. Einstein, and P. Hyldgaard, Phys. Rev. B 80, 155431 (2009).
- (91) P. Hyldgaard, J. Phys.: Condens. Matter 24, 424219 (2012).
- (92) B.M. Axilrod and E. Teller, J. Chem. Phys. 11, 299 (1943).
- (93) Y.U. Barash, Fiz. Tverd. Tela (Leningrad) 30, 1578 (1988) [Sov. Phys. Solid State 30, 1580 (1988)].
- (94) Y.U. Barash and O.I. Notysh, Zh. Exp. Tero. Fiz. 98, 502 (1990) [Sov. Phys. JETP 71, 301 (1990)].
- (95) B.E. Sernelius and P. Björk, Phys. Rev. B 57, 6592 (1998).
- (96) M. Boström and B.E. Sernelius, Phys. Rev. B 61, 2204 (2000).
- (97) J.F. Dobson, A. White, and A. Rubio, Phys. Rev. Lett. 96, 073201 (2006).
- (98) A. Tkatchenko and O.A. v. Lilienfeld, Phys. Rev. B 78, 045116 (2008).
- (99) S. Lebeque, J. Hart, T. Gould, J.G. Angyan, G. Kresse, and J.F. Dobson, Phys. Rev. Lett. 105, 196401 (2010).
- (100) A. Tkatchenko, R.A. DiStasio Jr., R. Car, and M. Scheffler, Phys. Rev. Lett. 108, 236402 (2012).
- (101) A. Ruzsinszky, J.P. Perdew, J. Tao, G.I. Csonka, and J.M. Pitarke, Phys. Rev. Lett. 109, 233203 (2012).
- (102) J.P. Perdew, J. Tao, P. Hao, A. Ruzsinszky, G.I. Csonka, and J.M. Pitarke, J. Phys.:Condens. Matter 24, 424207 (2012).
- (103) J. F. Dobson and T. Gould, J. Phys.:Condens. Matter 24, 073201 (2012).
- (104) V.V. Gobre and A. Tkatchenko, Nature Comm. 2013, DOI: 10.1038/ncomms3341.
- (105) J.F. Dobson, T. Gould, and G. Vignale, Phys. Rev. X, 021040 (2014).
- (106) J. F. Dobson, Int. J. Quantum Chem. 2014, DOI: 10.1002/qua.24635.
- (107) J. Kleis, E. Schröder, and P. Hyldgaard, Phys. Rev. B 77, 205422 (2008).
- (108) K. Lee, A. K. Kelkkanen, K. Berland, S. Andersson, D. C. Langreth, E. Schröder, B. I. Lundqvist, and P. Hyldgaard, Phys. Rev. B 84, 193408 (2011).
- (109) K. Lee, K. Berland, M. Yoon, S. Andersson, E. Schröder, P. Hyldgaard, and B. I. Lundqvist, J. Phys.:Condens. Matter 24, 424213 (2012).
- (110) The hollow C, formed by a shell with a near-metallic nature,Ruzsinszky12 ; PT supplements graphene and metallic nanotubes as examples where the modern expanded vdW-DF versions cannot capture the vdW interaction across all length scales. In C screening leads to the creation of a plasmon excitation of vanishing frequency and a per-fragment response that differs markedly from that obtained by treating the molecular fragments at the GGA level.Ruzsinszky12 ; PT When studying the vdW attraction at larger separations it is not possible to both truncate the screening and represent these hollow/extended structures as just a single molecular-type fragment with a GGA-type screening (as the recent vdW-DF versions do).Dobson06 ; Lebegue10 ; Ruzsinszky12 ; PT ; DoGoDisp12 ; DobsonMB14
- (111) S. D. Chakarova-Käck, A. Vojvodic, J. Kleis, P. Hyldgaard, and E. Schröder, New J. Phys. 12, 013017 (2010).