Per Hyldgaard

Van der Waals density functional: Self-consistent potential and the nature of the van der Waals bond

We derive the exchange-correlation potential corresponding to the nonlocal van der Waals density functional [M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)]. We use this potential for a self-consistent calculation of the ground state properties of a number of van der Waals complexes as well as crystalline silicon. For the latter, where little or no van der Waals interaction is expected, we find that the results are mostly determined by semilocal exchange and correlation as in standard generalized gradient approximations (GGA), with the fully nonlocal term giving little effect. On the other hand, our results for the van der Waals complexes show that the self-consistency has little effect on the atomic interaction energy and structure at equilibrium distances. This finding validates previous calculations with the same functional that treated the fully nonlocal term as a post-GGA perturbation. A comparison of our results with wave-function calculations demonstrates the usefulness of our approach. The exchange-correlation potential also allows us to calculate Hellmann-Feynman forces, hence providing the means for efficient geometry relaxations as well as unleashing the potential use of other standard techniques that depend on the self-consistent charge distribution. The nature of the van der Waals bond is discussed in terms of the self-consistent bonding charge.

Van der Waals density functional for layered structures.

To understand sparse systems, we must account for both strong local atom bonds and weak nonlocal van der Waals forces between atoms separated by empty space. A fully nonlocal functional form [Phys. Rev. B 62, 6997 (2000)]] of density-functional theory (DFT) is applied here to the layered systems graphite, boron nitride, and molybdenum sulfide to compute bond lengths, binding energies, and compressibilities. These key examples show that the DFT with the generalized-gradient approximation does not apply for calculating properties of sparse matter, while use of the fully nonlocal version appears to be one way to proceed.

van der Waals forces in density functional theory: a review of the vdW-DF method.

A density functional theory (DFT) that accounts for van der Waals (vdW) interactions in condensed matter, materials physics, chemistry, and biology is reviewed. The insights that led to the construction of the Rutgers-Chalmers van der Waals density functional (vdW-DF) are presented with the aim of giving a historical perspective, while also emphasizing more recent efforts which have sought to improve its accuracy. In addition to technical details, we discuss a range of recent applications that illustrate the necessity of including dispersion interactions in DFT. This review highlights the value of the vdW-DF method as a general-purpose method, not only for dispersion bound systems, but also in densely packed systems where these types of interactions are traditionally thought to be negligible.

Exchange functional that tests the robustness of the plasmon description of the van der Waals density functional

Is the plasmon description within the nonlocal correlation of the van der Waals density functional by Dion and coworkers (vdW-DF) robust enough to describe all exchange-correlation components? To address this question, we design an exchange functional based on this plasmon description as well as recent analysis on exchange in the large-s regime. In the regime with reduced gradients s = |del n|/2nk(F)(n) smaller than approximate to 2.5, dominating the nonlocal correlation part of the binding energy, the enhancement factor F-x(s) closely resembles the Langreth-Vosko screened exchange. In the s regime beyond, dominated by exchange, F-x(s) passes smoothly over to the revised Perdew-Wang-86 form. We term the specific exchange functional LV-PW86r, wheras the full van der Waals functional version emphasizing consistent handling of exchange is termed vdW-DF-cx. Our tests indicate that vdW-DF-cx produces accurate separations and binding energies of the S22 data set of molecular dimers as well as accurate lattice constants and bulk moduli of layered materials and tightly bound solids. These results suggest that the plasmon description within vdW-DF gives a good description of both exchange and correlation effects in the low-to-moderate s regime.

Zero-frequency current noise for the double-tunnel-junction Coulomb blockade.

We compute the zero-frequency current noise numerically and in several limits analytically for the Coulomb-blockade problem consisting of two tunnel junctions connected in series. At low temperatures over a wide range of voltages, capacitances, and resistances it is shown that the noise measures the variance in the number of electrons in the region between the two tunnel junctions. The average current, on the other hand, is linearly related to the mean number of electrons for an asymmetric pair of junctions. Thus, the noise provides additional information about transport in these devices which is not available from measuring the current alone.

Physisorption of nucleobases on graphene: a comparative van der Waals study

The physisorption of the nucleobases adenine (A), cytosine (C), guanine (G), thymine (T), and uracil (U) on graphene is studied using several variants of the density functional theory (DFT): the generalized gradient approximation with the inclusion of van der Waals interaction (vdW) based on the TS approach (Tkatchenko and Scheffer 2009 Phys. Rev. Lett. 102 073005) and our simplified version of this approach (here called sTS), the van der Waals density functional vdW-DF (Dion et al 2004 Phys. Rev. Lett. 92 246401) and vdW-DF2 (Lee et al 2010 Phys. Rev. B 82 081101), and DFT-D2 (Grimme 2006 J. Comput. Chem. 27 1787) and DFT-D3 (Grimme et al 2010 J. Chem. Phys. 132 154104) methods. The binding energies of nucleobases on graphene are found to be in the following order: G > A > T > C > U within TS, sTS, vdW-DF, and DFT-D2, and in the following order: G > A > T ~ C > U within DFT-D3 and vdW-DF2. The binding separations are found to be different within different methods and in the following order: DFT-D2 < TS < DFT-D3 ~ vdW-DF2 < vdW-DF. We also comment on the efficiency of combining the DFT-D approach and vdW-DF to study systems with van der Waals interactions.

Potassium intercalation in graphite : A van der Waals density-functional study

Potassium intercalation in graphite is investigated by first-principles theory. The bonding in the potassium-graphite compound is reasonably well accounted for by traditional semilocal density functional theory (DFT) calculations. However, to investigate the intercalate formation energy from pure potassium atoms and graphite requires use of a description of the graphite interlayer binding and thus a consistent account of the nonlocal dispersive interactions. This is included seamlessly with ordinary DFT by a van der Waals density functional (vdW-DF) approach [Phys. Rev. Lett. 92, 246401 (2004)]. The use of the vdW-DF is found to stabilize the graphite crystal, with crystal parameters in fair agreement with experiments. For graphite and potassium-intercalated graphite structural parameters such as binding separation, layer binding energy, formation energy, and bulk modulus are reported. Also the adsorption and sub-surface potassium absorption energies are reported. The vdW-DF description, compared with the traditional semilocal approach, is found to weakly soften the elastic response.

Long-ranged adsorbate-adsorbate interactions mediated by a surface-state band

The adsorbate-adsorbate interaction mediated by Shockley surface states on, for example, the (111) faces of noble metals, yields an oscillatory form modulated by a 1/d 2 envelope at asymptotic adsorbate separations d . For this interaction we obtain a non-perturbative analytical estimate, specified by experimentally accessible Shockley-state parameters and the finite Fermi-level phase shift F 0, which characterize the standing-wave patterns observed in scanning tunnelling microscopy (STM) images. We provide explicit interaction results for the phase shift value F = - /2 suggested by the STM measurements of sulphur adsorbates on Cu(111), and we attempt to relate our results to the corresponding observations of correlations in the adsorbate distribution function.

Spin Signature of Nonlocal Correlation Binding in Metal-Organic Frameworks

We develop a proper nonempirical spin-density formalism for the van der Waals density functional (vdW-DF) method. We show that this generalization, termed svdW-DF, is firmly rooted in the single-particle nature of exchange and we test it on a range of spin systems. We investigate in detail the role of spin in the nonlocal correlation driven adsorption of H_{2} and CO_{2} in the linear magnets Mn-MOF74, Fe-MOF74, Co-MOF74, and Ni-MOF74. In all cases, we find that spin plays a significant role during the adsorption process despite the general weakness of the molecular-magnetic responses. The case of CO_{2} adsorption in Ni-MOF74 is particularly interesting, as the inclusion of spin effects results in an increased attraction, opposite to what the diamagnetic nature of CO_{2} would suggest. We explain this counterintuitive result, tracking the behavior to a coincidental hybridization of the O p states with the Ni d states in the down-spin channel. More generally, by providing insight on nonlocal correlation in concert with spin effects, our nonempirical svdW-DF method opens the door for a deeper understanding of weak nonlocal magnetic interactions.

Site determination and thermally assisted tunneling in homogenous nucleation

A combined low-temperature scanning tunneling microscopy and density functional theory study on the binding and diffusion of copper monomers, dimers, and trimers adsorbed on Cu(111) is presented. Whereas atoms in trimers are found in fcc sites only, monomers as well as atoms in dimers can occupy the fcc as well as the metastable hcp site. In fact the dimer fcc-hcp configuration is only 1.3 meV less favorable with respect to the fcc-fcc configuration. This enables a confined intracell dimer motion, which at temperatures below 5 K is dominated by thermally assisted tunneling.