A van der Waals density functional for solids

Torbjörn Björkman
DOI: https://doi.org/10.1103/PhysRevB.86.165109
2012-09-14
Abstract:The recent non-local correlation functional of Vydrov and van Voorhis[J. Chem. Phys. 133, 244103 (2010)] is investigated and two new versions of the functional are suggested as being appropriate for describing van der Waals interactions in solids. A refitting of the original functional is demonstrated to result in very accurate interlayer binding energies for weakly bonded layered solids. A VV10 functional based on the generalized gradient approximation by Armiento and Mattsson[Phys. Rev. B 72, 085108 (2005)], while performing slightly worse for interlayer binding is highly successful in describing the equilibrium geometries of both weakly bonded and close packed solids.
Materials Science,Chemical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to improve the existing non - local exchange - correlation (XC) functionals so that they can describe the van der Waals (vdW) interactions in solids more accurately. Specifically, the authors studied the non - local correlation functional proposed by Vydrov and van Voorhis (VV10) and proposed two new functional versions to better describe the vdW interactions in solids. ### Specific Problems and Background 1. **Limitations of Existing Methods** - The VV10 functional was originally designed to describe the vdW interactions between molecules, but in weakly - bonded layered solids, it significantly overestimates the inter - layer binding energy. - Previous vdW correction methods cannot accurately predict the geometric structure and binding energy simultaneously in weakly - bonded layered solids, indicating that there is room for improvement in these methods. 2. **Improvement Goals** - Refit the parameters of the VV10 functional to make it more suitable for describing the vdW interactions in solids. - Ensure that the new functionals can provide more accurate inter - layer binding energy and equilibrium geometric structure in weakly - bonded layered solids while maintaining a good description of strongly - bonded solids. ### Method Overview The authors selected three generalized gradient approximation (GGA) functionals as parent functionals: PW86R, AM05, and PBEsol. By evaluating the performance of these functionals, PW86R and AM05 were finally selected as the basis for refitting the VV10 functional. The refitted functionals are called PW86R - VV10sol and AM05 - VV10sol respectively. ### Key Formulas The energy expression of the functional is divided into two parts: \[ E_c = E_c^0+E_c^{nl} \] where \(E_c^0\) is the local correlation energy and \(E_c^{nl}\) is the non - local correlation energy, in the form of: \[ E_c^{nl}=\frac{\hbar^2}{2}\iint dr \, dr' \, n(r)\Phi(r, r')n(r') \] Here, \(n(r)\) is the electron density and \(\Phi(r, r')\) describes the density - density interaction. In the VV10 framework, \(\Phi(r, r')\) depends on the local response parameter \(\omega_0(r)\) and the short - range damping parameter \(\kappa(r)\), as well as the local bandgap \(\omega_g(r)\). ### Conclusion By refitting the VV10 functional, the authors obtained two new functional versions: PW86R - VV10sol and AM05 - VV10sol. In particular, AM05 - VV10sol performs excellently in weakly - bonded layered solids and can simultaneously provide accurate inter - layer binding energy and equilibrium geometric structure, which is suitable for most solid - state calculations involving vdW interactions.