Statistical Mechanics Improves Density Functional Theory for Electrocatalytic Metal Surface Properties and CO Binding Trends

Herein, we combine Density Functional Theory (DFT) with ab initio molecular dynamics at ambient conditions to calculate surface properties of and CO binding trends to the electrocatalytically relevant Pt(111), Cu(111), Ag(111), and Au(111) metals. By accounting for thermal effects both the RPBE GGA and the B97M-rV meta-GGA functionals show better metal surface strain trends, CO binding energetics, and the correct low-coordination site preference for the Cu(111) and Pt(111) surfaces. The importance of dispersion is evident at finite temperatures, with only B97M-rV exhibiting binding energetics in agreement with experiment for the Ag(111) and Au(111) surfaces, and is found to consist of a mixture of chemisorbed and physisorbed CO species. Overall this work begins to address a longstanding theoretical discrepancy with experimental observations on these important metals relevant to the electrochemical CO2 reduction reaction through the full accounting of thermal fluctuations.