Elsebeth Schröder

Van der Waals density functional for general geometries

A scheme within density functional theory is proposed that provides a practical way to generalize to unrestricted geometries the method applied with some success to layered geometries [Phys. Rev. Lett. 91, 126402 (2003)]]. It includes van der Waals forces in a seamless fashion. By expansion to second order in a carefully chosen quantity contained in the long-range part of the correlation functional, the nonlocal correlations are expressed in terms of a density-density interaction formula. It contains a relatively simple parametrized kernel, with parameters determined by the local density and its gradient. The proposed functional is applied to rare gas and benzene dimers, where it is shown to give a realistic description.

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.

Application of van der Waals Density Functional to an Extended System: Adsorption of Benzene and Naphthalene on Graphite

It is shown that it is now possible to include van der Waals (vdW) interactions via a nonempirical implementation of density functional (DF) theory to describe the correlation energy in electronic structure calculations on infinite systems of no particular symmetry. The vdW-DF theory [Phys. Rev. Lett. 92, 246401 (2004)] is applied to the adsorption of benzene and naphthalene on an infinite sheet of graphite, as well as the binding between two graphite sheets. A comparison with recent thermal-desorption data [Phys. Rev. B 69, 155406 (2004)] shows great promise for the vdW-DF method.

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.

Stacking interactions and the twist of DNA.

The importance of stacking interactions for the Twist and stability of DNA is investigated using the fully ab initio van der Waals density functional (vdW-DF). Our results highlight the role that binary interactions between adjacent sets of base pairs play in defining the sequence-dependent Twists observed in high-resolution experiments. Furthermore, they demonstrate that additional stability gained by the presence of thymine is due to methyl interactions with neighboring bases, thus adding to our understanding of the mechanisms that contribute to the relative stability of DNA and RNA. Our mapping of the energy required to twist each of the 10 unique base pair steps should provide valuable information for future studies of nucleic acid stability and dynamics. The method introduced will enable the nonempirical theoretical study of significantly larger pieces of DNA or DNA/amino acid complexes than previously possible.

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.

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.

van der Waals density functionals built upon the electron-gas tradition: Facing the challenge of competing interactions

The theoretical description of sparse matter attracts much interest, in particular for those ground-state properties that can be described by density functional theory. One proposed approach, the van der Waals density functional (vdW-DF) method, rests on strong physical foundations and offers simple yet accurate and robust functionals. A very recent functional within this method called vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014)] stands out in its attempt to use an exchange energy derived from the same plasmon-based theory from which the nonlocal correlation energy was derived. Encouraged by its good performance for solids, layered materials, and aromatic molecules, we apply it to several systems that are characterized by competing interactions. These include the ferroelectric response in PbTiO3, the adsorption of small molecules within metal-organic frameworks, the graphite/diamond phase transition, and the adsorption of an aromatic-molecule on the Ag(111) surface. Our results indicate that vdW-DF-cx is overall well suited to tackle these challenging systems. In addition to being a competitive density functional for sparse matter, the vdW-DF-cx construction presents a more robust general-purpose functional that could be applied to a range of materials problems with a variety of competing interactions.

Binding of polycyclic aromatic hydrocarbons and graphene dimers in density functional theory

An early van der Waals density functional (vdW-DF) described layered systems (such as graphite and graphene dimers) using a layer-averaged electron density in the evaluation of nonlocal correlations. This early vdW-DF version was also adapted to approximate the binding of polycyclic aromatic hydrocarbons (PAHs) (Chakarova S D and Schroder E 2005 J. Chem. Phys. 122 054102). In parallel to that PAH study, a new vdW-DF version (Dion M, Rydberg H, Schroder E, Langreth D C and Lundqvist B I 2004 Phys. Rev. Lett. 92 246401) was developed that provides accounts of nonlocal correlations for systems of general geometry. We apply here the latter vdW-DF version to aromatic dimers of benzene, naphthalene, anthracene and pyrene, stacked in sandwich (AA) structure, and the slipped-parallel (AB) naphthalene dimer. We further compare the results of the two methods as well as other theoretical results obtained by quantum-chemistry methods. We also compare calculations for two interacting graphene sheets in the AA and the AB structures and provide the corresponding graphene-from-graphite exfoliation energies. Finally, we present an overview of the scaling of the molecular–dimer interaction with the number of carbon atoms and with the number of carbon rings.