Gold-based (Au) nanostructures are efficient catalysts for CO oxidation, hydrogen evolution (HER), and oxygen evolution (OER) reactions, but stabilizing them on graphene (Gr) is challenging due to weak affinity from delocalized \(p_{z}\) carbon orbitals. This study investigates forming metal alloys to enhance stability and catalytic performance of Au-based nanocatalysts. Using ab initio density functional theory, we characterize \({\text {M}_{(n-x)}\text {Au}_{x}}\) sub-nanoclusters (M = Ni, Pd, Pt, Cu, and Ag) with atomicities \(n=1-4\) , both in gas-phase and supported on Gr. We find that M atoms act as “anchors,” enhancing binding to Gr and modulating catalytic efficiency. Notably, \({\text {Pt}_{(n-x)}\text {Au}_{x}}\) /Gr shows improved stability, with segregation tendencies mitigated upon adsorption on Gr. The d-band center ( \(\varepsilon _{\text {d}}\) ) model indicates catalytic potential, correlating an optimal \(\varepsilon _{\text {d}}\) range of \(-1 \text { to }-2\) eV for HER and OER catalysts. Incorporating Au into \({\text{M}_n}\) adjusts \(\varepsilon _{\text {d}}\) closer to the Fermi level, especially for Group-10 alloys, offering designs with improved stability and efficiency comparable to pure Au nanocatalysts. Our methodology leveraged SimStack, a workflow framework enabling modeling and analysis, enhancing reproducibility, and accelerating discovery. This work demonstrates SimStack’s pivotal role in advancing the understanding of composition-dependent stability and catalytic properties of Au-alloy clusters, providing a systematic approach to optimize metal-support interactions in catalytic applications.
