Abstract
Single-atom alloys have recently emerged as highly active and selective alloy catalysts. Unlike pure metals, single-atom alloys escape the well-established conceptual framework developed nearly three decades ago for predicting catalytic performance. Although this offers the opportunity to explore so far unattainable chemistries, this leaves us without a simple guide for the design of single-atom alloys able to catalyse targeted reactions. Here, based on thousands of density functional theory calculations, we reveal a 10-electron count rule for the binding of adsorbates on the dopant atoms, usually the active sites, of single-atom alloy surfaces. A simple molecular orbital approach rationalizes this rule and the nature of the adsorbate–dopant interaction. In addition, our intuitive model can accelerate the rational design of single-atom alloy catalysts. Indeed, we illustrate how the unique insights provided by the electron count rule help identify the most promising dopant for an industrially relevant hydrogenation reaction, thereby reducing the number of potential materials by more than one order of magnitude.
Similar content being viewed by others
Main
Single-atom alloys (SAAs) have recently emerged as a new class of catalysts able to reach high activity and selectivity for a range of chemical reactions1,2,3. In these alloys, the active metal is dispersed as single atoms at the surface of a more inert host metal. This doping strategy significantly improves the catalytic performance of the otherwise poorly reactive, and yet highly selective, coinage metals (Cu, Ag and Au). There has been considerable work aiming at screening and predicting the unique reactivity and catalytic performance of SAAs2,4,5. Despite the accuracy of the predictions based on density functional theory (DFT), a simple physical model describing general trends is still lacking due to the special electronic structure of SAAs1,2.
Traditionally, the catalytic activity of a material can be described using the binding energy of species involved in the mechanism of the catalysed reaction. According to the Sabatier principle, species should bind neither too weakly nor too strongly to ensure optimal catalytic activity6, and the electronic properties of the catalyst determine the binding energies. For transition metal catalysts, the adsorption energy of a species linearly correlates with the energy centre of the electronic band consisting of the d-states of the metal7,8. This means that the more valence electrons the metal has, the lower the energy of the d-band centre and the weaker the bond to an adsorbate at the surface (Supplementary Fig. 1). Although the d-band model has shown great success in understanding catalytic performance over a few decades, there are known limitations. For example, SAAs and the related class of near-surface alloys (NSAs) can both escape the trends expected from the d-band model9,10,11,12. Recent developments have shown that NSAs follow stability rules and require corrections to the d-band model to fully describe their properties9,13. Yet, such an approach is not applicable to SAAs. Alternatively, the rise of machine learning has provided an efficient approach for the prediction of adsorption energies on traditional alloys and SAAs5,14,15,16. Albeit faster than and reportedly as accurate as DFT calculations, these models do not provide the fundamental physical principles that govern the stability of adsorbates on SAA surfaces. In a change of perspective, a few studies have suggested to consider SAAs as analogues of molecular systems17,18,19,20. The stability of such systems is related to the filling of discrete states (atomic orbitals or molecular orbitals (MO)) resulting in various electron-counting rules (octet, 4n + 2 aromaticity rule and so on)21,22,23.
Here, we propose a change of paradigm in the way we rationalize the reactivity of SAAs. Moving away from the traditional linear scaling relationships that have limited applicability for SAAs, we demonstrate that an MO approach, albeit somewhat simplistic, provides profound insight into the metal–adsorbate binding mechanism. Screening a variety of catalytically relevant adsorbates on a large set of SAA surfaces, we show that the adsorbate–dopant interaction is strongest when the vM valence electrons of the metal dopant (equivalently its group number in the periodic table) and the k valence electrons of the adsorbates interacting with the dopant sum up to ten: vM + k = 10. This 10-electron count rule is supported by a detailed analysis of the electronic structure of adsorbates bound to SAAs surfaces and generalized to molecular adsorbates. Finally, we demonstrate that, without expensive DFT calculations or complex machine learning models, this rule provides, experimentalists and theoreticians alike, insightful guidance for the design of SAA catalysts for targeted reactions of industrial significance, as illustrated herein for the reduction of nitrogen.
Results and discussion
Adsorption of adatoms
We start by analysing the periodic trends in the adsorption energies of key atomic adsorbates. Specifically, in Fig. 1 and Supplementary Fig. 2, we show the adsorption energies computed with DFT for atoms (O, N, C and H) adsorbed on a whole range of transition metal dopants on SAA surfaces. Unlike the adsorption on pure transition metals8,24, we do not observe the monotonic weakening of the adsorbates’ binding from left to right along the periodic table (Supplementary Fig. 1a,b), despite similar trends of the centre of the d-states of the dopants (Supplementary Fig. 3). Instead, we observe shallow W-shaped trends for adsorption energies on 3d metal dopants and deep V-shaped trends for adsorption energies on 4d and 5d dopants. Interestingly, the position of the minimum on 4d and 5d dopants depends on the adsorbate, but not the dopant’s period nor the host material. When the number of valence electrons of the adsorbate (vA) decreases (O > N > C > H), the minimum shifts to metal dopants with more valence electrons vM (to the right of the period), as if the dopant had to compensate for the fewer electrons brought by the adsorbate.
To further analyse the interplay between the number of electrons brought by the adsorbate and the dopant, we can tentatively count the total number of valence electrons of the dopant–adsorbate system, ignoring the host metal. For the dopants, we consider the number of electrons of the s- and d-orbitals of the outer shell, which corresponds to the group number of the element. For adsorbates, we consider all the electrons of the nsnp outer shell: these are the electrons traditionally considered when drawing Lewis structures. Now, if we add the valence electrons of the dopant and the adsorbate, we find that strongest binding is associated with 10 electrons for H and 12 electrons for p-block adsorbates (O, N and C). For 3d dopants, we can still recognize a preference towards binding to early transition metals for electron-rich adsorbates (O and N), and to later transition metals for adsorbates with fewer electrons (H and C). However, because of significant spin effects on magnetic 3d dopants, adsorption is weakened in the middle of the row, with a fixed extremum for Mn for all adsorbates. Indeed, if the effect of spin is unrealistically suppressed, 3d dopants behave like 4d and 5d dopants (Supplementary Fig. 4). Interestingly, these binding features on SAA surfaces (shape of the trends, role of the number of valence electrons) are distinct from pure transition metal surfaces but analogous to the binding of ligands in organometallic complexes (as shown in Supplementary Fig. 1) successfully modelled using an MO approach23.
This analogy with organometallic complexes and the electronic structure of SAAs akin to gas phase atoms17,18 have motivated us to consider an MO approach to rationalize the binding of adsorbates on SAAs (Fig. 2). Based on the symmetry point group C∞ν of the metal–adsorbate pair, we can construct an MO diagram (Supplementary Fig. 5). For H, the s-orbital can interact with the metal’s dz² orbital, generating a bonding σ and an antibonding σ* MO (Fig. 2a). The remaining four d-orbitals form the nonbonding nδ and nπ MOs. We can therefore fill up to five MOs before populating the antibonding σ*-orbital and thereby weakening the adsorbate–metal bond (Fig. 2a). Similarly, we can construct the MO diagrams for C, N and O. These adsorbates have partially filled p-orbitals that can contribute to the bond. The linear combination of the metal dz² and the adsorbate’s s and pz orbitals generate a bonding σ MO, a nonbonding nσ MO and an antibonding σ* MO. Additionally, the linear combination of dxz and dyz with the px and py, form two bonding π and two antibonding π* MOs. The remaining two d-orbitals, which cannot interact with the adsorbate’s orbitals, form the nonbonding nδ MOs. For p-block adsorbates, we can fill up to six MOs, that is, twelve electrons, before populating antibonding orbitals (Fig. 2e). This MO approach seems to accurately predict the point when the dopant–adsorbate’s binding weakens. This is under the assumption that the dopant’s s-orbital can be taken out of the picture and its electron(s) populate the states that we have just built.
To ensure that our MO approach appropriately captures the binding of adsorbates on SAAs (position of the orbitals and population), we further analysed the electronic structure of H (Fig. 2c,e,g) and N (Fig. 2d,f,h) adsorbed on Ag-based SAAs. When constructing MO diagrams, we are interested in knowing (1) where the MOs are on the energy scale (y-axis), and (2) which atomic orbitals interact to form a particular MO (that is, atomic orbitals with nonzero overlap interaction). This information can be extracted, from DFT calculations, using the electronic density of states (DOS), that is, the number of states found at a certain energy level, projected onto the d-states of the dopant (pDOS), and the crystal orbital Hamilton population (COHP) that quantifies the overlap interaction between two orbitals25,26. Fig. 2c shows the COHP analysis between the dz² orbital of Zr and the s-orbital of H as extracted from the electronic structure calculation of H adsorbed on ZrAg. The analysis shows a narrow positive peak at −2 eV for the pairwise interaction between the s and dz² orbitals; this is the expected bonding σ MO. The same plot shows antibonding σ* states as a negative broad band. The four other d-orbitals, referred to as nπ and nδ in Fig. 2a, do not interact with any orbitals of H, and thus cannot be identified in the COHP analysis that only shows pairwise interactions. Instead, they appear in the plot showing the density of nπ and nδ states (Fig. 2c, right) as partially populated, as expected from the MO diagram. When considering H adsorbed on RhAg (Fig. 2e), all the states shift to lower energies, with the σ, nπ and nδ orbitals now being doubly occupied and the σ* partially occupied. This results from Rh having more electrons than Zr. The same analysis for adsorbed N is shown in Fig. 2d,f,h. The COHP analysis shows that the σ-orbital indeed results from the interaction of the dz² orbital with the adsorbate’s s and pz orbitals. Approximately 7 eV higher up in energy are found the π orbitals (blue curve). In the same energy region, the COHP analysis shows that the dz² orbital interacts again with the adsorbate’s s and pz orbitals. This time, however, the binding contribution of the pz–dz² interaction (pink curve) cancels out the antibonding contribution of the s–dz² interaction: the superposition of these two opposite contributions corresponds to the expected nonbonding nσ MO. Overall, the proposed MO diagrams agree well with the computed electronic structures. Now that we have identified the different orbitals, we can estimate their filling. PdH (Fig. 2g) and RhN (Fig. 2h) have more than 10 and 12 valence electrons, respectively (vA + vM). Although the dopant’s s states—which interact with the host’s s states to form a band—are half-filled (Supplementary Table 1), there are more than vM − 1 electrons in the orbitals considered in Fig. 2. This can be attributed to charge transfer between the host and the dopant resulting in more electrons than perhaps expected, especially around the point of strongest binding and beyond. All the valence electrons of the metal (vM) should therefore be counted to predict when antibonding states start being populated. This whole analysis also provides insights into the binding of adsorbates on 3d metals. The electronic structure of N adsorbed on 3d dopants shows that the filling of antibonding states starts as early as Cr, the dopant for which the binding shows signs of weakening (Fig. 1) and intensifies for Co and Ni (Supplementary Fig. 6). This is consistent with the observed W-shape trend of the adsorption energies for 3d dopants and our proposed analogy to organometallic complexes.
The common point between H and p-block elements adsorbed on SAAs is that adsorption is the strongest when all the bonding and nonbonding MOs with d contributions are filled (Fig. 2g,h). The adsorption energy starts weakening when antibonding states are being populated. But why does filling nonbonding states seemingly strengthen the bond? Nonbonding states do not contribute to increased stabilization. What drives the stabilization is the lowering of the dopant’s levels and the reduction of the dopant’s radius which, together, lead to a better overlap with the adsorbate’s small orbitals at shorter distances (Supplementary Fig. 7). Now, if H and p-block elements seem to obey different counting rules, it is only because too many electrons are counted in the latter case. The nσ orbital, which can be interpreted as a lone pair located on the adsorbate outside the internuclear region, is consistently populated throughout the period and does not contribute to the binding or the saturation of the d orbitals. Ignoring this lone pair results in a universal 10-electron count rule for maximal binding of adsorbates on 4d and 5d dopants (Fig. 2g,h). It is essential to distinguish between the adsorbate’s total number of valence electrons vA and the number of valence electrons k that interact with the d states, especially when considering larger adsorbates such as molecules (Table 1). For example, CO has vA = 10 valence electrons, but, as expected from organometallic chemistry, only k = 2 electrons interact with the transition metal23,27. With this distinction in mind, the dopants with maximal binding to a specific adsorbate are easily identified as those with 10 − k valence electrons, hereafter referred to as d10−k dopants.
Extension to molecular fragments
We have shown the 10-electron count rule for the stability of atomic adsorbates on SAAs. Now, we extend this rule to molecular adsorbates. For example, CO interacts with its lone pair located on the carbon atom (Fig. 3b). NO, as a radical, can either interact with one electron in a bent geometry or three electrons in a linear geometry (Fig. 3b). If we only consider the linear geometry, the 10-electron rule predicts the strongest adsorption energies at d8 for CO and d7 for NO on 4d and 5d dopants. This is confirmed by our DFT calculations (Fig. 3a). It is important to note, however, that this rule only holds when the binding trends between the adsorbate and the dopant are dominated by the sharing of electrons, that is, covalent contributions. H2O and NH3 are known examples of adsorbates for which the binding trends are dominated by electrostatic interactions20. At the electronic level, this translates to a limited reshuffling of the electronic density of H2O upon adsorption compared with NO (Fig. 3d). For these adsorbates, the DFT-computed adsorption energies do not show the expected minima and roughly follow the steady variation of the electric charge of the dopant (Fig. 3c and Supplementary Table 2) over the 4d period (from −0.22e to +1.61e). It is only when we remove the electrostatic contribution, that we recover the expected trends for the covalent contribution with minima at d8 for NH3 and H2O (Supplementary Fig. 8). For these closed-shell molecules with highly polarized bonds (large electronegativity difference between O or N and H) as well as halogens and hydroxyl (Supplementary Fig. 9), the 10-electron count rule has less practical applicability and the atomic charge of the dopant is a much more robust descriptor of the binding20.
The significance of the 10-electron rule goes beyond understanding periodic trends in binding on SAAs: it helps identify active catalysts for targeted reactions. Let us consider the reduction of nitrogen to ammonia, a reaction of industrial relevance. On Au-based SAAs, the first elementary step, namely the hydrogenation of N2 to diazenyl NNH (Fig. 4a), was identified as rate-determining because of its endothermicity28. N2 interacts with the dopant via its lone pair and is predicted by the 10-electron count rule to have strongest binding for d8 (Ru, Os) on 4d and 5d dopants (Fig. 4c,d, black line). NNH, isolobal to NO, is predicted to have the strongest adsorption for d7 (Tc, Re) on 4d and 5d dopants. DFT calculations, shown in Fig. 4c,d, again confirm these predictions. Now, let us consider the trends in reaction energy required to go from the reactant (black curve) to the intermediate (green curve). Because of the rigidity of the periodic trends on 3d doped surfaces (W-shape with fixed maximum for Mn), the reaction energy does not significantly change when screening 3d dopants. However, for 4d and 5d dopants, the situation is different. Because of the curvatures and the different position of adsorption minima for N2 and NNH, the gap between the two curves closes when the binding energy of NNH is the strongest. We therefore predict the first hydrogenation of N2 to be most facile on d7 dopants, that is, Re and Tc. For reasons ranging from synthesizability, cost, to stability (Tc is not a stable isotope, it is only considered to test and illustrate the counting rule), it is worth considering dopants around the d7 minimum. Moving to the left of the periodic table (d7 to d6) provides viable options (Mo and W) with small thermodynamic barriers. Moving to the right (d7 to d8), however, is less interesting as the stability of the reactant N2 reaches its maximum for d8, thereby detrimentally increasing the reaction energy.
Interestingly, previous computationally demanding high-throughput and machine-learning studies28,29 had identified similar SAAs for the reduction of N2. According to a previous study, Re dopants offer the best compromise between donation and back-donation to weaken the N≡N bond29. Our theoretical framework offers an alternative explanation: this is the stability of NNH that controls the choice of the dopant regarding the most favourable reaction energetics, although donation and backdonation may, admittedly, play a role for the scission of the bond between the two nitrogen atoms in the subsequent steps of the reduction mechanism. Of course, this analysis is an early step towards the demonstration of a new catalyst for ammonia synthesis. However, this example illustrates the power of our approach: the counting rule remarkably identifies, without access to supercomputers, the most promising catalytic sites, and the degree to which we can extend the search of efficient catalysts to neighbouring dopants. Once identified, computational effort (multiscale modelling including effects such as solvation, electric fields, kinetics) or experimental effort can be specifically directed to the most promising materials.
In the previous paragraph we have only focused on the linear geometry of NNH, which is the most stable in the region of catalytic interest. We have identified two other geometries: the flat-lying and the bent geometries. The former binds to the dopant bringing five electrons and therefore helps early transition metals to satisfy the 10-electron rule, especially for dopants with fewer than 5 valence electrons. The bent geometry, interacting with one electron, is the only geometry found for Pd and Pt dopants. For these late transition metals, the bent geometry, with fewer electrons interacting with the d orbitals of the dopants, prevents too many antibonding states from being filled. The preference of early, mid and late transition metals for the flat-lying, linear and bent geometries is again a manifestation of the 10-electron count rule (Supplementary Fig. 10). Similar behaviour is well known in molecular and organometallic chemistry and was theorized using orbital correlation diagrams20,30,31.
Conclusion
In summary, species covalently bound at the surface of SAAs show non-monotonic binding energy trends that mostly depend on the nature of the dopant site. Our work clearly identifies the dopants with stronger or weaker affinity for a given adsorbate. For 3d dopants, binding is the strongest on Ti and V, or Fe and Co, regardless of the nature of the adsorbate. For 4d and 5d dopants, the trends show a single extremum with strongest binding when the adsorbate’s electrons fill and saturate the dopant’s d orbitals, hence the 10-electron count rule. In this case, the dopants with strongest binding depend on the nature of the adsorbate. We show that a simple molecular orbital approach rationalizes these rules. This furthers our previous work that had successfully identified the dopant charge as a descriptor for the binding of species with large electrostatic contributions but had originally failed at providing an electronic-level descriptor for the covalent contribution20. The fundamental principles identified here advance our understanding of chemical bonds on SAAs and provide a sought-after alternative to the d-band model for binding trends on the more reactive dopant sites (the d-band model still remains relevant for host sites). Finally, our model bridges the gap between the concepts of heterogeneous and homogeneous catalysis and, significantly so, establishes a clear conceptual guide, for theoreticians and experimentalists alike, for the design of more efficient catalysts, without the expensive or complex simulations only accessible to computational scientists.
Methods
All DFT calculations were performed using Vienna Ab Initio Simulation Package (VASP) version 5.4.4 (refs. 32,33) with the projector-augmented wave34 method to model core ionic potentials. The nonlocal optB86b-vdW exchange-correlation functional35,36 was used for all electronic structure calculations except the systems with CO and NO adsorption. The plane wave cut-off was set to 400 eV, electronic convergence set to 10−6 eV and geometric convergence criterion was less than 0.02 eV Å−1. The surfaces were modelled using a five-layer p(3×3) slab, with the bottom two layers fixed at the positions of the bulk host. The pure fcc host metals were optimized in bulk and lattice constants of 3.60 Å for Cu-based SAAs, 4.09 Å for Ag-based SAAs and 4.13 Å for Au-based SAAs were used. The slabs were separated by a 15 Å vacuum layer. For the screening of adsorption energies, a 3×3×1 Monkhorst−Pack37 k-point mesh was chosen for Brillouin-zone integration, which was shown to yield the same trends as a fully converged 13×13×1 Monkhorst-Pack k-point mesh (Supplementary Fig. 11). All adsorbates were placed at the atop position and optimized using the quasi-Newton algorithm implemented in VASP, to prevent relaxation of some adsorbates to the fcc site. Overall trends are preserved for fcc-site adsorption (Supplementary Fig. 12). The adsorption energies ΔE were calculated from the dipole corrected slabs with reference to gas-phase atoms. Spin-unrestricted calculations were performed for all SAA slabs, but for better convergence the charge density of a spin-restricted calculation was calculated first and used as starting guess for the spin-unrestricted calculation. Adsorption energies of CO and NO were calculated with the RPBE functional (see Supplementary Fig. 13 for comparison with optB86b-vdW)38.
Analysis of the electronic structure, that is, COHP analysis and pDOS, was performed using Lobster26,39,40,41,42. Starting from a VASP calculation with a higher 13×13×1 Monkhorst–Pack k-point mesh and symmetry switched off (ISYM = −1) the WAVECAR and other VASP output files where postprocessed using Lobster. Bader charges were determined using the VTST tools developed previously43.
The electronic population nj of MO j can be determined summing over pDOSi of its contributing AOs i. For bonding MOs, the integration over the energy 𝜖 is performed up to E±, the point where the COHP changes signs, that is, the point where the linear combination of AOs switches from being bonding to being antibonding (equation (1)), or up to the Fermi level EF if EF < E±. For antibonding MOs, the integration is performed from E± up to EF (equation (2)). If E± > EF, the antibonding MO is empty (null population). Similar integrations are performed for the nonbonding d and nσ orbitals (equations (3) and (4)).
Figure 3d was obtained using the VMD software (available at http://www.ks.uiuc.edu/Research/vmd/)44.
Data availability
The DFT calculations dataset used in this study is publicly available in the NOMAD Repository (https://doi.org/10.17172/NOMAD/2023.12.04-2). Source data are provided with this paper.
References
Hannagan, R. T., Giannakakis, G., Flytzani-Stephanopoulos, M. & Sykes, E. C. H. Single-atom alloy catalysis. Chem. Rev. 120, 12044–12088 (2020).
Réocreux, R. & Stamatakis, M. One decade of computational studies on single-atom alloys: is in silico design within reach? Acc. Chem. Res. 55, 87–97 (2022).
Zhang, T., Walsh, A. G., Yu, J. & Zhang, P. Single-atom alloy catalysts: structural analysis, electronic properties and catalytic activities. Chem. Soc. Rev. 50, 569–588 (2021).
Darby, M. T. et al. Elucidating the stability and reactivity of surface intermediates on single-atom alloy catalysts. ACS Catal. 8, 5038–5050 (2018).
Kumar, A. et al. Machine learning enabled screening of single atom alloys: predicting reactivity trend for ethanol dehydrogenation. ChemCatChem 14, e202101481 (2022).
Sabatier, P. Hydrogénations et déshydrogénations par catalyse. Ber. Dtsch. Chem. Ges. 44, 1984–2001 (1911).
Hammer, B. & Norskov, J. K. Why gold is the noblest of all the metals. Nature 376, 238–240 (1995).
Hammer, B. & Nørskov, J. K. Theoretical surface science and catalysis—calculations and concepts. Adv. Catal. 45, 71–129 (2000).
Saini, S., Halldin Stenlid, J. & Abild-Pedersen, F. Electronic structure factors and the importance of adsorbate effects in chemisorption on surface alloys. npj Comput. Mater. 8, 163 (2022).
Fung, V., Hu, G. & Sumpter, B. Electronic band contraction induced low temperature methane activation on metal alloys. J. Mater. Chem. A Mater 8, 6057–6066 (2020).
Monasterial, A. P., Hinderks, C. A., Viriyavaree, S. & Montemore, M. M. When more is less: nonmonotonic trends in adsorption on clusters in alloy surfaces. J. Chem. Phys. 153, 111102 (2020).
Sun, Z., Song, Z. & Yin, W.-J. Going beyond the d-band center to describe CO2 activation on single‐atom alloys. Adv. Energy Sustainability Res. 3, 2100152 (2022).
Calle-Vallejo, F., Martínez, J. I., García-Lastra, J. M., Rossmeisl, J. & Koper, M. T. M. Physical and chemical nature of the scaling relations between adsorption energies of atoms on metal surfaces. Phys. Rev. Lett. 108, 116103 (2012).
Yang, Z. & Gao, W. Applications of machine learning in alloy catalysts: rational selection and future development of descriptors. Adv. Sci. 9, 2106043 (2022).
Andersen, M., Levchenko, S. V., Scheffler, M. & Reuter, K. Beyond scaling relations for the description of catalytic materials. ACS Catal. 9, 2752–2759 (2019).
Toyao, T. et al. Toward effective utilization of methane: Machine learning prediction of adsorption energies on metal alloys. J. Phys. Chem. C 122, 8315–8326 (2018).
Greiner, M. T. et al. Free-atom-like d states in single-atom alloy catalysts. Nat. Chem. 10, 1008–1015 (2018).
Thirumalai, H. & Kitchin, J. R. Investigating the reactivity of single atom alloys using density functional theory. Top. Catal. 61, 462–474 (2018).
Spivey, T. D. & Holewinski, A. Selective interactions between free-atom-like d-states in single-atom alloy catalysts and near-frontier molecular orbitals. J. Am. Chem. Soc. 143, 11897–11902 (2021).
Réocreux, R., Sykes, E. C. H., Michaelides, A. & Stamatakis, M. Stick or Spill? Scaling relationships for the binding energies of adsorbates on single-atom alloy catalysts. J. Phys. Chem. Lett. 13, 7314–7319 (2022).
Lewis, G. N. The atom and the molecule. J. Am. Chem. Soc. 38, 762–785 (1916).
Hückel, E. Quantentheoretische Beiträge zum Benzolproblem. Z. Phys. 70, 204–286 (1931).
Housecroft, C. E. & Sharpe, A. G. Inorganic Chemistry (Pearson Education, 2012).
Gao, W. et al. Determining the adsorption energies of small molecules with the intrinsic properties of adsorbates and substrates. Nat. Commun. 11, 1196 (2020).
Hughbanks, T. & Hoffmann, R. Chains of trans-edge-sharing molybdenum octahedra: metal-metal bonding in extended systems. J. Am. Chem. Soc. 105, 3528–3537 (1983).
Dronskowski, R. & Bloechl, P. E. Crystal orbital Hamilton populations (COHP): energy-resolved visualization of chemical bonding in solids based on density-functional calculations. J. Phys. Chem. 97, 8617–8624 (1993).
Green, M. L. H. & Parkin, G. Application of the covalent bond classification method for the teaching of inorganic chemistry. J. Chem. Educ. 91, 807–816 (2014).
Zheng, G. et al. High-throughput screening of a single-atom alloy for electroreduction of dinitrogen to ammonia. ACS Appl. Mater. Interfaces 13, 16336–16344 (2021).
Dai, T., Wang, Z., Lang, X. & Jiang, Q. ‘Sabatier principle’ of d electron number for describing the nitrogen reduction reaction performance of single-atom alloy catalysts. J. Mater. Chem. A Mater. 10, 16900–16907 (2022).
Mingos, D. M. P. General bonding model for linear and bent transition metal-nitrosyl complexes. Inorg. Chem. 12, 1209–1211 (1973).
Hoffmann, R., Chen, M. M. L. & Thorn, D. L. Qualitative discussion of alternative coordination modes of diatomic ligands in transition metal complexes. Inorg. Chem. 16, 503–511 (1977).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758–1775 (1999).
Klimeš, J., Bowler, D. R. & Michaelides, A. Chemical accuracy for the van der Waals density functional. J. Phys. Condens. Matter 22, 022201 (2010).
Klimeš, J., Bowler, D. R. & Michaelides, A. Van der Waals density functionals applied to solids. Phys. Rev. B Condens. Matter Mater. Phys. 83, 195131 (2011).
Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 13, 5188–5192 (1976).
Hammer, B., Hansen, L. B. & Nørskov, J. K. Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 59, 7413–7421 (1999).
Deringer, V. L., Tchougréeff, A. L. & Dronskowski, R. Crystal orbital hamilton population (COHP) analysis as projected from plane-wave basis sets. J. Phys. Chem. A 115, 5461–5466 (2011).
Nelson, R. et al. LOBSTER: local orbital projections, atomic charges, and chemical‐bonding analysis from projector‐augmented‐wave‐based density‐functional theory. J. Comput. Chem. 41, 1931–1940 (2020).
Maintz, S., Deringer, V. L., Tchougréeff, A. L. & Dronskowski, R. Analytic projection from plane-wave and PAW wavefunctions and application to chemical-bonding analysis in solids. J. Comput. Chem. 34, 2557–2567 (2013).
Maintz, S., Deringer, V. L., Tchougréeff, A. L. & Dronskowski, R. LOBSTER: a tool to extract chemical bonding from plane-wave based DFT. J. Comput. Chem. 37, 1030–1035 (2016).
Tang, W., Sanville, E. & Henkelman, G. A grid-based Bader analysis algorithm without lattice bias. J. Phys. Condens. Matter 21, 084204 (2009).
Humphrey, W., Dalke, A. & Schulten, K. VMD: visual molecular dynamics. J. Mol. Graph. 14, 33–38 (1996).
Acknowledgements
We thank E. C. H. Sykes for discussions and feedback and members of the ICE group for discussions. Via our membership of the UK High-End Computing Materials Chemistry Consortium, which is funded by the Engineering and Physical Sciences Research Council (EP/R029431), this work used the UK Materials and Molecular Modelling Hub for computational resources, which is partially funded by Engineering and Physical Sciences Research Council (EP/T022212 and EP/P020194). We thank University College London high-performance computing (Myriad@UCL) for the use of their facilities and support. J.S. is grateful for support from a Feodor Lynen Research Fellowship from the Alexander von Humboldt Foundation. R.R., A.M. and M.S. acknowledge financial support from the Leverhulme Trust (grant RPG-2018-209), and R.R. and M.S. gratefully acknowledge funding from the European Union’s Horizon 2020 research and innovation programme under grant 814416 (ReaxPro).
Author information
Authors and Affiliations
Contributions
J.S. and R.R. performed DFT calculations. R.R. developped the MO diagrams. J.S. performed the DOS and COHP analyses. All authors analysed the compuational data. All authors wrote, read and commented on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Chemistry thanks Adam Holewinski and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Supplementary Information
Supplementary Figs. 1–13, Tables 1–3 and discussion.
Source data
Source Data Fig. 1
Binding energies for each material.
Source Data Fig. 2
DOS and orbital filling.
Source Data Fig. 3
Adsorption energies for CO, NO, H2O and NH3; atomic charges for each material.
Source Data Fig. 4
Formation energies for N2 + H and NNH.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Schumann, J., Stamatakis, M., Michaelides, A. et al. Ten-electron count rule for the binding of adsorbates on single-atom alloy catalysts. Nat. Chem. (2024). https://doi.org/10.1038/s41557-023-01424-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41557-023-01424-6