Molecular adsorption on coinage metal subnanoclusters: A DFT+D3 investigation

Gold and silver subnanoclusters with few atoms are prominent candidates for catalysis-related applications, primarily because of the large fraction of lower-coordinated atoms exposed and ready to interact with external chemical species. However, an in-depth energetic analysis is necessary to characterize the relevant terms within the molecular adsorption process that can frame the interactions within the Sabatier principle. Herein, we investigate the interaction between Ag-n and Au-n subnanoclusters (clu, n = 2-7) and N-2, NO, CO, and O-2 molecules, using scalar-relativistic density functional theory calculations within van der Waals D3 corrections. The onefold top site is preferred for all chemisorption cases, with a predominance of linear (approximate to 180 degrees) and bent (approximate to 120 degrees) molecular geometries. A larger magnitude of adsorption energy is correlated with smaller distances between molecules and clusters and with the weakening of the adsorbates bond strength represented by the increase of the equilibrium distances and decrease of molecular stretching frequencies. From the energetic decomposition, the interaction energy term was established as an excellent descriptor to classify subnanoclusters in the adsorption/desorption process concomitant with the Sabatier principle. The limiting cases: (i) weak molecular adsorption on the subnanoclusters, which may compromise the reaction activation, where an interaction energy magnitude close to 0 eV is observed (e.g., physisorption in N-2/Ag-6); and (ii) strong molecular interactions with the subnanoclusters, given the interaction energy magnitude is larger than at least one of the individual fragment binding energies (e.g., strong chemisorption in CO/Au-4 and NO/Au-4), conferring a decrease in the desorption rate and an increase in the possible poisoning rate. However, the intermediate cases are promising by involving interaction energy magnitudes between zero and fragment binding energies. Following the molecular closed-shell (open-shell) electronic configuration, we find a predominant electrostatic (covalent) nature of the physical interactions for N-2 center dot center dot center dot clu and CO center dot center dot center dot clu (O-2 center dot center dot center dot clu and NO center dot center dot center dot clu), except in the physisorption case (N-2/Ag-6) where dispersive interaction is dominant. Our results clarify questions about the molecular adsorption on subnanoclusters as a relevant mechanistic step present in nanocatalytic reactions.