Insight into the description of van der Waals forces for benzene adsorption on transition metal (111) surfaces

TL;DR

This study compares two advanced van der Waals-inclusive density functional theory methods to accurately model benzene adsorption on various transition metal (111) surfaces, highlighting their strengths and limitations.

Contribution

It provides a detailed comparison of interatomic vdW and nonlocal functional approaches for benzene adsorption, emphasizing the importance of electrodynamic screening and optimized exchange functionals.

Findings

Both methods accurately predict adsorption energies and equilibrium distances.

Explicit electrodynamic screening improves interatomic vdW results.

Discrepancies occur at large adsorbate-substrate separations.

Abstract

Exploring the role of van der Waals (vdW) forces on the adsorption of molecules on extended metal surfaces has become possible in recent years thanks to exciting developments in density functional theory (DFT). Among these newly developed vdW-inclusive methods, interatomic vdW approaches that account for the nonlocal screening within the bulk [V. G. Ruiz, W. Liu, E. Zojer, M. Scheffler, and A. Tkatchenko, Phys. Rev. Lett. 108, 146103 (2012)] and improved nonlocal functionals [J. Klimes, D. R. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 022201(2010)] have emerged as promising candidates to account efficiently and accurately for the lack of long-range vdW forces in most popular DFT exchange-correlation functionals. Here we have used these two approaches to compute benzene adsorption on a range of close-packed (111) surfaces upon which it either physisorbs (Cu, Ag, and Au) or chemisorbs (Rh, Pd, Ir, and Pt). We have thoroughly compared the performance between the two classes of vdW-inclusive methods and when available compared the results obtained with experimental data. By examining the computed adsorption energies, equilibrium distances, and binding curves we conclude that both methods allow for an accurate treatment of adsorption at equilibrium adsorbate-substrate distances. To this end, explicit inclusion of electrodynamic screening in the interatomic vdW scheme and optimized exchange functionals in the case of nonlocal vdW density functionals is mandatory. Nevertheless, some discrepancies are found between these two classes of methods at large adsorbate-substrate separations.