We present a microkinetic model for CO(2) reduction (CO(2)R) on Cu(211) towards C2 products, based on energetics estimated from an explicit solvent model. We show that the differences in both Tafel slopes and pH dependence for C1 vs C2 activity arise from differences in their multi-step mechanisms. We find the depletion in C2 products observed at high overpotential and high pH to arise from the 2nd order dependence of C-C coupling on CO coverage, which decreases due to competition from the C1 pathway. We further demonstrate that CO(2) reduction at a fixed pH yield similar activities, due to the facile kinetics for CO2 reduction to CO on Cu, which suggests C2 products to be favored for CO2R under alkaline conditions. The mechanistic insights of this work elucidate how reaction conditions can lead to significant enhancements in selectivity and activity towards higher value C2 products.

Bimetallic catalysts are promising for the most difficult thermal and electrochemical reactions but modeling the many diverse active sites on polycrystalline samples is an open challenge. We present a general framework for addressing this complexity in a systematic and predictive fashion. Active sites for every stable low-index facet of a 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 sites. The activity of these sites is explored in parallel using a neural-network based surrogate model to share information between the many Density Functional Theory (DFT) relaxations, resulting in activity estimates with an order of magnitude fewer explicit DFT calculations. Sites with interesting activity were found and provide targets for follow-up calculations. This process was applied to the electrochemical reduction of CO 2 on nickel gallium bimetallics and indicated that most facets had similar activity to Ni surfaces, but a few exposed Ni sites with a very favorable on-top CO configuration. This motif emerged naturally from the predictive modeling and represents a class of intermetallic CO 2 reduction catalysts. These sites rationalize recent experimental reports of nickel gallium activity and why previous materials screens missed this exciting material. Most importantly these methods suggest that bimetallic catalysts will be discovered by studying facet reactivity and diversity of active sites more systematically.