We present a machine-learning-augmented chemisorption model that enables fast and accurate prediction of the surface reactivity of metal alloys within a broad chemical space. Specifically, we show that artificial neural networks, a family of biologically inspired learning algorithms, trained with a set of ab initio adsorption energies and electronic fingerprints of idealized bimetallic surfaces, can capture complex, nonlinear interactions of adsorbates (e.g., *CO) on multimetallics with ∼0.1 eV error, outperforming the two-level interaction model in prediction. By leveraging scaling relations between adsorption energies of similar adsorbates, we illustrate that this integrated approach greatly facilitates high-throughput catalyst screening and, as a specific case, suggests promising {100}-terminated multimetallic alloys with improved efficiency and selectivity for CO2 electrochemical reduction to C2 species. Statistical analysis of the network response to perturbations of input features underpins our fundamental understanding of chemical bonding on metal surfaces.
Keywords: alloys; artificial neural networks; carbon dioxide reduction; density functional theory; machine learning; reactivity descriptors.