Adsorption and Vibrational Spectroscopy of CO on the Surface of MgO from Periodic Local Coupled-Cluster Theory

TL;DR

This study uses advanced periodic coupled-cluster calculations to accurately predict CO adsorption properties on MgO surfaces, validating the method against experimental vibrational frequency shifts and assessing various computational approaches.

Contribution

It introduces a cost-effective periodic local coupled-cluster approach for surface chemistry, providing highly accurate adsorption energies and vibrational frequencies.

Findings

CCSD(T) accurately predicts positive vibrational frequency shift (+14.7 cm^{-1})

MP2 and CCSD give reasonable adsorption energies and shifts

Many DFT methods overestimate adsorption energy and predict negative frequency shifts

Abstract

The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of $+14.7~\textrm{cm}^{-1}$, in excellent agreement with the experimental shift of $+14.3~\textrm{cm}^{-1}$. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by $5$ -- $15$~kJ/mol and predict a negative frequency shift of about $-20~\textrm{cm}^{-1}$, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.