Ligand removal energetics control CO 2 electroreduction selectivity on atomically precise, ligated alloy nanoclusters
Atomically precise, thiolate-protected gold nanoclusters (TPNCs) exhibit remarkable catalytic performance for the electrochemical reduction of carbon dioxide (CO 2 R) to CO. The origin of their high CO 2 R activity and selectivity has been attributed to partial ligand removal from the thiolate-covered surfaces of TPNCs to expose catalytically active sulfur atoms. Recently, heterometal doped (alloy) TPNCs have been shown to exhibit enhanced CO 2 R activity and selectivity compared to their monometallic counterparts. However, systematic studies on the effect of doping (metal type and location on TPNC) on active site exposure and CO 2 R activity are missing in literature. Herein, we apply Density Functional Theory calculations to investigate the effect of heterometal (Pt, Pd, Hg and Cd) doping of Au 25 (SR) 18 TPNC on the active site exposure and CO 2 R activity and selectivity. In this study, we reveal that doping significantly modifies relevant TPNC electronic properties, such as electron affinity, while also altering partial ligand removal and carboxyl (*COOH) intermediate formation energies. Furthermore, we demonstrate that changing the dopant (e.g. Hg) position can change the selectivity of the TPNC towards CO (g) or H 2(g) formation, highlighting the importance of dopant locations in TPNC-based CO 2 R. Most notably, we report a universal (i.e. capturing different dopant types and positions) linear trend between the ligand removal energy and i) the *COOH formation energy, as well as, ii) the hydrogen (*H) formation energy on the different alloy TPNCs. Thus, utilizing the ligand removal energy as a descriptor for CO 2 RR activity and selectivity, our work opens new avenues for accelerated computational screening of different alloy TPNCs for electrocatalytic CO 2 R applications.