A Combined Density Functional Theory and Microkinetics Simulation Study of Electrochemical CO2 Reduction on Ceria-Supported Bismuth

High Resolution Image Download MS PowerPoint Slide Direct electrochemical CO 2 reduction (ECR) into carbon-based fuels and chemicals is a promising way to upgrade waste CO 2 with renewable energy, contributing to closing carbon cycles and mitigating climate change. Here, we investigate the ECR of Bi–CeO 2 catalysts. Using a combination of density functional theory (DFT), artificial neural networks (ANN), genetic algorithms (GA), and microkinetics simulations, we conducted a comprehensive exploration of active sites, the CO 2 -to-formic acid (HCOOH) mechanism and the electrochemical behavior of Bi x /CeO 2 catalysts. Three representative models were investigated: (i) a Bi atom adsorption on CeO 2 (Bi 1 /CeO 2 ), (ii) a single Bi atom doped in the CeO 2 surface (Bi 1 –CeO 2 ), and (iii) a small cluster of eight Bi atoms adsorbed on CeO 2 (Bi 8 /CeO 2 ). ANN-GA was employed to identify the optimal structure of the Bi 8 clusters on the CeO 2 surface. Our investigation shows various reaction pathways for converting CO 2 to HCOOH and CO. For the structural models featuring Bi on the surface, the HCOO pathway toward HCOOH is the predominant one. Bi doping in CeO 2 predominantly favors the COOH pathway, resulting in CO as the main product. The former models leading to HCOOH exhibit higher current densities than the doped model, which mainly produces CO. Electronic structure analysis shows that stronger electron donation in the Bi 1 /CeO 2 model enhances HCOOH current densities by weakening the O–H bond and stabilizing the transition state. We discuss the kinetic differences in current density and selectivity as a function of the electrochemical potential. These findings not only elucidate various CO 2 conversion pathways, which can explain the formation of desirable HCOOH and unwanted CO, but also offer theoretical guidance for the design of electrocatalysts.