A six-dimensional potential energy surface for Ru(0001)(2×2):CO

Abstract

We present a new global ground state potential energy surface (PES) for carbon monoxide at a coverage of 1/4, on a rigid Ru(0001) surface [Ru(0001)(2×2):CO]. All six adsorbate degrees of freedom are considered. For constructing the PES, we make use of more than 90 000 points calculated with periodic density functional theory using the RPBE exchange-correlation functional and an empirical van der Waals correction. These points are used for interpolation, utilizing a symmetry-adapted corrugation reducing procedure (CRP). Three different interpolation schemes with increasing accuracy have been realized, giving rise to three flavours of the CRP PES. The CRP PES yields in agreement with the DFT reference and experiments, the atop position of CO to be the most stable adsorption geometry, for the most accurate interpolation with an adsorption energy of 1.69 eV. The CRP PES shows that diffusion parallel to the surface is hindered by a barrier of 430 meV, and that dissociation is facilitated but still activated. As a first “real” application and further test of the new potential, the six-dimensional vibrational Schrödinger equation is solved variationally to arrive at fully coupled, anharmonic frequencies and vibrational wavefunctions for the vibrating, adsorbed CO molecule. Good agreement with experiment is found also here. Being analytical, the new PES opens an efficient way towards multidimensional dynamics.