The effect of the exchange-correlation functional on H2 dissociation on Ru(0001)

The specific reaction parameter (SRP) approach to density functional theory (DFT) has enabled a chemically accurate description of reactive scattering experiments for activated H2-metal systems (H2 + Cu(111) and Cu(100)), but its application has not yet resulted in a similarly accurate description of non-activated or weakly activated H2-metal systems. In this study, the effect of the choice of the exchange-correlation functional in DFT on the potential energy surface and dynamics of H2 dissociation on Ru(0001), a weakly activated system, is investigated. In total, full potential energy surfaces were calculated for over 20 different functionals. The functionals investigated include functionals incorporating an approximate description of the van der Waals dispersion in the correlation functional (vdW-DF and vdW-DF2 functionals), as well as the revTPSS meta-GGA. With two of the functionals investigated here, which include vdW-DF and vdW-DF2 correlation, it has been possible to accurately reproduce molecular beam experiments on sticking of H2 and D2, as these functionals yield a reaction probability curve with an appropriate energy width. Diffraction probabilities computed with these two functionals are however too high compared to experimental diffraction probabilities, which are extrapolated from surface temperatures (Ts) ⩾ 500 K to 0 K using a Debye-Waller model. Further research is needed to establish whether this constitutes a failure of the two candidate SRP functionals or a failure of the Debye-Waller model, the use of which can perhaps in future be avoided by performing calculations that include the effect of surface atom displacement or motion, and thereby of the experimental Ts.