Improving the Accuracy of Modelling CO 2 Electroreduction on Copper Using Many-Body Perturbation Theory.

Copper (Cu) remains the most important metal catalyst for the carbon dioxide reduction reaction (CO 2 RR) into C 2 products. Due to limited evidence from in situ experiments, mechanistic studies are often performed in the framework of density functional theory (DFT), using functionals at the generalized gradient approximation (GGA) level, which have fundamental difficulties to correctly describe CO adsorption and surface stability. We employ the adiabatic connection fluctuation dissipation theorem within the random phase approximation (RPA), in combination with the linearized Poisson-Boltzmann equation to describe solvation effects, to investigate the mechanism of CO 2 RR on the Cu(100) facet. Qualitatively different from the DFT-GGA results, RPA results propose the formation of *OCCHO as the potential determining step towards C 2 products. The results suggest that it is important to use more accurate methods like RPA when modeling reactions involving multiple CO-related species like CO 2 RR.