Jacob’s Ladder as Sketched by Escher: Assessing the Performance of Broadly Used Density Functionals on Transition Metal Surface Properties

The present work surveys the performance of various widely used density functional theory exchange–correlation (xc) functionals in describing observable surface properties of a total of 27 transition metals with face-centered cubic (fcc), body-centered cubic (bcc), or hexagonal close-packed (hcp) crystallographic structures. A total of 81 low Miller index surfaces were considered employing slab models. Exemplary xc functionals within the three first rungs of Jacob’s ladder were considered, including the Vosko–Wilk–Nusair xc functional within the local density approximation, the Perdew–Burke–Ernzerhof (PBE) functional within the generalized gradient approximation (GGA), and the Tao–Perdew–Staroverov–Scuseria functional as a meta-GGA functional. Hybrids were excluded in the survey because they are known to fail in properly describing metallic systems. In addition, two variants of PBE were considered, PBE adapted for solids (PBEsol) and revised PBE (RPBE), aimed at improving adsorption energies. Interlayer atomic distances, surface energies, and surface work functions were chosen as the scrutinized properties. A comparison with available experimental data, including single-crystal and polycrystalline values, shows that no xc functional is best at describing all of the surface properties. However, in statistical mean terms the PBEsol xc functional is advised, while PBE is recommended when considering both bulk and surface properties. On the basis of the present results, a discussion of adapting GGA functionals to the treatment of metallic surfaces in an alternative way to meta-GGA or hybrids is provided.