Next-generation non-local van der Waals density functional

The fundamental ideas for a non-local density functional theory -- capable of reliably capturing van der Waals interaction -- were already conceived in the 1990s. In 2004, a seminal paper introduced the first practical non-local exchange-correlation functional called vdW-DF, which has become widely successful and laid the foundation for much further research. However, since then, the functional form of vdW-DF has remained unchanged. Several successful modifications paired the original functional with different (local) exchange functionals to improve performance and the successor vdW-DF2 also updated one internal parameter. Bringing together different insights from almost two decades of development and testing, we present the next-generation non-local correlation functional called vdW-DF3, in which we change the functional form while staying true to the original design philosophy. Although many popular functionals show good performance around the binding separation of van der Waals complexes, they often result in significant errors at larger separations. With vdW-DF3, we address this problem by taking advantage of a recently uncovered and largely unconstrained degree of freedom within the vdW-DF framework that can be constrained through empirical input, making our functional semi-empirical. For two different parameterizations, we benchmark vdW-DF3 against a large set of well-studied test cases and compare our results with the most popular functionals, finding good performance in general for a wide array of systems and a significant improvement in accuracy at larger separations. Finally, we discuss the achievable performance within the current vdW-DF framework, the flexibility in functional design offered by vdW-DF3, as well as possible future directions for non-local van der Waals density functional theory.

Non-trivial retardation effects in dispersion forces: From anomalous distance dependence to novel traps

In the study of dispersion forces, nonretarded, retarded and thermal asymptotes with their distinct scaling laws are regarded as cornerstone results governing interactions at different separations. Here, we show that when particles interact in a medium, the influence of retardation is qualitatively different, making it necessary to consider the non-monotonous potential in full. We discuss different regimes for several cases and find an anomalous behaviour of the retarded asymptote. It can change sign, and lead to a trapping potential.

Thermoelectric transport trends in group 4 half-Heusler alloys

The thermoelectric properties of 54 different group 4 half-Heusler (HH) alloys have been studied from first principles. Electronic transport was studied with density functional theory using hybrid functionals facilitated by the $mathbf{k} cdot mathbf{p}$ method, while the temperature dependent effective potential method was used for the phonon contributions to the figure of merit $ZT$. The phonon thermal conductivity was calculated including anharmonic phonon-phonon, isotope, alloy and grain-boundary scattering. HH alloys have an ${it XYZ}$ composition and those studied here are in the group 4-9-15 (Ti,Zr,Hf)(Co,Rh,Ir)(As,Sb,Bi) and group 4-10-14 (Ti,Zr,Hf)(Ni,Pd,Pt)(Ge,Sn,Pb). The electronic part of the thermal conductivity was found to significantly impact $ZT$ and thus the optimal doping level. Furthermore, the choice of functional was found to significantly affect thermoelectric properties, particularly for structures exhibiting band alignment features. The intrinsic thermal conductivity was significantly reduced when alloy and grain boundary scattering were accounted for, which also reduced the spread in thermal conductivity. It was found that sub-lattice disorder on the ${it Z}$-site, i.e. the site occupied by group 14 or 15 elements, was more effective than ${it X}$-site substitution, occupied by group 4 elements. The calculations confirmed that ZrNiSn, ZrCoSb and ZrCoBi based alloys display promising thermoelectric properties. A few other n-type and p-type compounds were also predicted to be potentially excellent thermoelectric materials, given that sufficiently high charge carrier concentrations can be achieved. This study provides insight into the thermoelectric potential of HH alloys and casts light on strategies to optimize thermoelectric performance of multicomponent alloys.

Assessment of two hybrid van der Waals density functionals for covalent and non-covalent binding of molecules

Two hybrid van der Waals density functionals (vdW-DFs) are constructed using 25%, Fock exchange with i) the consistent-exchange vdW-DF-cx functional and ii) with the vdW-DF2 functional. The ability to describe covalent and non-covalent binding properties of molecules are assessed. For properties related to covalent binding, atomization energies (G2-1 set), molecular reaction energies (G2RC set), as well as ionization energies (G21IP set) are benchmarked against experimental reference values. We find that hybrid-vdW-DF-cx yields results that are rather similar to those of the standard non-empirical hybrid PBE0 [JCP 110, 6158 (1996)]. Hybrid vdW-DF2 follows somewhat different trends, showing on average significantly larger deviations from the reference energies, with a MAD of 14.5 kcal/mol for the G2-1 set. Non-covalent binding properties of molecules are assessed using the S22 benchmark set of non-covalently bonded dimers and the X40 set of dimers of small halogenated molecules, using wavefunction-based quantum chemistry results for references. For the S22 set, hybrid-vdW-DF-cx performs better than standard vdW-DF-cx for the mostly hydrogen-bonded systems. Hybrid-vdW-DF2 offers a slight improvement over standard vdW-DF2. Similar trends are found for the X40 set, with hybrid-vdW-DF-cx performing particularly well for binding involving the strongly polar hydrogen halides, but poorly for systems with tiny binding energies. Our study of the X40 set reveals both the potential of mixing Fock exchange with vdW-DF, but also highlights shortcomings of the hybrids constructed here. The solid performance of hybrid-vdW-DF-cx for covalent-bonded systems, as well as the strengths and issues uncovered for non-covalently bonded systems, makes this study a good starting point for developing even more precise hybrid vdW-DFs.

Interpretation of van der Waals density functionals

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.

A Harris-type van der Waals density functional scheme

Large biomolecular systems, whose function may involve thousands of atoms, cannot easily be addressed with parameter-free density functional theory (DFT) calculations. Until recently a central problem was that such systems possess an inherent sparseness, that is, they are formed from components that are mutually separated by low-electron-density regions where dispersive forces contribute significantly to the cohesion and behavior. The introduction of, for example, the van der Waals density functional (vdW-DF) method [PRL 92, 246401 (2004)] has addressed part of this sparse-matter system challenge. However, while a vdW-DF study is often as computationally efficient as a study performed in the generalized gradient approximation, the scope of large-sparse-matter DFT is still limited by computer time and memory. It is costly to self-consistently determine the electron wavefunctions and hence the kinetic-energy repulsion. In this paper we propose and evaluate an adaption of the Harris scheme [PRB 31, 1770 (1985)]. This is done to speed up non-selfconsistent vdW-DF studies of molecular-system interaction energies. Also, the Harris-type analysis establishes a formal link between dispersion-interaction effects on the effective potential for electron dynamics and the impact of including selfconsistency in vdW-DF calculations [PRB 76, 125112 (2007)].

Benchmarking van der Waals Density Functionals with Experimental Data: Potential Energy Curves for H2 Molecules on Cu(111), (100), and (110) Surfaces

Detailed physisorption data from experiment for the H_2 molecule on low-index Cu surfaces challenge theory. Recently, density-functional theory (DFT) has been developed to account for nonlocal correlation effects, including van der Waals (dispersion) forces. We show that the functional vdW-DF2 gives a potential-energy curve, potential-well energy levels, and difference in lateral corrugation promisingly close to the results obtained by resonant elastic backscattering-diffraction experiments. The backscattering barrier is found selective for choice of exchange-functional approximation. Further, the DFT-D3 and TS-vdW corrections to traditional DFT formulations are also benchmarked, and deviations are analyzed.

Temperature stability of intersubband transitions in AlN/GaN quantum wells

Temperature dependence of intersubband transitions in AlN/GaN multiple quantum wells grown with molecular beam epitaxy is investigated both by absorption studies at different temperatures and modeling of conduction-band electrons. For the absorption study, the sample is heated in increments up to $400^circ$C. The self-consistent Schrodinger-Poisson modeling includes temperature effects of the band-gap and the influence of thermal expansion on the piezoelectric field. We find that the intersubband absorption energy decreases only by $sim 6$ meV at $400^circ$C relative to its room temperature value.

A van der Waals density functional study of adenine on graphite: Single molecular adsorption and overlayer binding

The adsorption of an adenine molecule on graphene is studied using a first-principles van der Waals functional (vdW-DF) [Dion et al., Phys. Rev. Lett. 92, 246401 (2004)]. The cohesive energy of an ordered adenine overlayer is also estimated. For the adsorption of a single molecule, we determine the optimal binding configuration and adsorption energy by translating and rotating the molecule. The adsorption energy for a single molecule of adenine is found to be 711 meV, which is close to the calculated adsorption energy of the similar-sized naphthalene. Based on the single molecular binding configuration, we estimate the cohesive energy of a two-dimensional ordered overlayer. We find a significantly stronger binding energy for the ordered overlayer than for single-molecule adsorption.