Design of Single Atom Catalysts

ABSTRACT Over the past decade, computational modeling based on density functional theory (DFT) calculations provides a deep insight into the catalytic mechanism of single-atom catalysts (SACs) and paves way for high-throughput screening of promising SACs. This review summarizes computational methods for the analysis of the electronic structures and catalytic performance of SACs, as well as introduces the utilization of descriptors for the computational design of SACs. We expect that future advances in computational methods will surely help to identify highly effective SACs for a wide variety of reactions. Graphical Abstract


Introduction
Widely applied in more than 90% of chemical processes, heterogeneous catalysts have attracted the most extensive attention in the industrial production of the majority of chemicals [1,2]. Traditional heterogeneous catalysts of supported nanoparticles have played an important role in the efficient production of chemicals. However, they may suffer from low atom utilization efficiency and low selectivity for desired products in a lot of reactions. To handle these drawbacks, scientists have turned their eyes to atomically dispersed metal atoms anchored on supports, which are defined as single-atom catalysts (SACs) [3]. The distinctive structures of SACs make them very promising as heterogeneous catalysts. Firstly, the identical active sites of SACs enable them to selectively catalyze only a certain kind of reaction. Moreover, as all of the atoms of SACs are exposed to reactants, their atom utilization efficiency is close to 100%, which can increase the catalytic activity per atom and reduce the consumption of metal resources [4]. In addition, the relatively simple structure of SACs provides scientists with an ideal model to explore the relationship between the active sites and their catalytic performance, paving the way for the rational design of SACs [5]. Upon these advantages, SACs have been experimentally fabricated for various reactions, showing promising prospects in catalysis [6][7][8][9][10].
Although SACs with high catalytic performance have been identified experimentally, the intrinsic factors that determine their catalytic performance remain unclear. Compared with experimental studies, computational modeling by density functional theory (DFT) calculations could offer more explicit insights into the reaction processes. Based on DFT calculations, the whole reaction cycle and the energy barriers of each elementary step can be calculated. Moreover, DFT calculations could be used to analyze the electronic structures of the active sites and reveal their structure-performance relationship. Favored with these advantages, DFT calculations have been widely applied to the study of SACs [11][12][13][14].
In this review, we summarize computational methods that open avenues for the rational design of SACs. We first introduce DFT calculations of the electronic structures of SACs. Then, we summarize calculation methods for evaluating the catalytic performance of SACs, including stability, activity, and selectivity. Finally, we describe the development of high-throughput methods for designing SACs based on descriptors, as shown in Figure 1. It is believed that these advanced computational methods and descriptors could be used to guide the experimental studies and design efficient SACs for a wide variety of reactions.

DFT calculations of SACs
Through DFT calculations, the electronic structure of SACs can be theoretically revealed. The electronic analysis is useful for the construction of descriptors based on the identification of the intrinsic properties determining the catalytic performance of SACs. Moreover, DFT calculations can predict the stability, activity, and selectivity of SACs. These predictive results are the criteria for the evaluation of the accuracy of identified descriptors. Thus, DFT calculations of the intrinsic properties and catalytic performance of SACs are the prerequisites for the computational design of SACs.

Electronic structures
The performance of catalysts is in principle determined completely by electronic structures [17]. The eximious performance of SACs in a series of reactions has been elucidated from the analysis at the electronic scale [18][19][20][21]. Moreover, revealing the electronic factors that determine the catalytic performance may point out directions for tuning the catalytic performance of SACs. (Reprinted from Ref [28]. Copyright 2017 American Chemical Society). (c) Schematic of CO oxidation on MoS 2 supported Co SAC. (d) Schematic of titanium nitride supported Pt SAC selectively catalyze oxygen reduction and formic acid oxidation reaction. (Reprinted from Ref [15]. Copyright 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim). (e) Volcano relation between adsorption free energy of OOH and limiting potential of SACs for oxygen electroreduction reaction. (Reprinted from Ref [66]. 2019 American Chemical Society). (f) Relation between adsorption free energy of OH and structure descriptor φ of graphene supported transition metals. (Reprinted from Ref [27]. Copyright 2018 The Author(s), under exclusive licence to Macmillan Publishers Ltd, part of Springer Nature). (g) DFT calculated adsorption energy versus predictive adsorption energy for transition metals on oxide supports (γ-Al 2 O 3 , MgO, and MgAl 2 O 4 ). (Reprinted from Ref [16]. Copyright 2019 American Chemical Society).
The analysis of the electronic structures can be used to understand the stability and activity of SACs. As illustrated in Figure 2, the stability and activity of SACs are in a competitive relationship. The stability of SACs corresponds to the fully occupied bonding orbitals [22]. The more hybridization between the orbitals of the single atoms and the support can lead to stronger binding strength. For example, partial density of states (PDOS) of Pt@nitrogen-doped graphene shows that the discrete 5d-orbitals of Pt interact with the N-2p orbitals around the Fermi level [20]. Charge transfers from Pt to N atom, and the calculated adsorption energy of Pt (5.171 eV) clearly suggests that it is stable, which is in accordance with experiment [20]. The orbital hybridization between the single atom and the support also induces the strong metal-support interactions (SMSIs) within oxide-supported single metal atoms [23]. Due to SMSI, single atoms can stably stay on the supports. For instance, DFT calculations and experiments prove that Au single atoms are ultrastable on iron oxide nanocrystallites [24]. The Au atoms are bound to the support by strong covalent metalsupport interactions (CMSI) [24]. Bader charge analysis and calculated valence electron configurations show that the anchored Au single atoms are highly positively charged, which enables them to interact strongly with the negatively charged support [24]. The SMSI and electrostatic interaction thus account for the high stability of Au single atoms against sintering. On the other hand, the activity of SACs corresponds to the free electrons near the Fermi level because of the unsaturated coordination [22]. The free electrons determine the activity of SACs by tuning the binding strength between them and the adsorbates. As it is well-known that the activation of adsorbates can be ascribed to the moderate binding energy, namely, not too weak to make the activation energy of reactants too high or too strong to let the products hard to desorb. For instance, after dilute Ag by Cu, the activation barrier of bulk Cu for H 2 formation in methanol reforming is reduced by 0.1 eV [25]. This is partially because the dilution narrows the 3d state of Cu and contributes to the covalent bond formation when Cu interacts with hydrogen. The Cu in the AgCu alloy as the minority element is found to exhibit free-atom like d-state structure, and this phenomenon is also found for several other single-atom alloys (SAAs) -the structure of metal supported single metal atoms, pointing out a new approach to SAAs that may exhibit excellent catalytic performance [25].

Evaluation of stability
Structural stability is the prerequisite for industrial application of an effective SAC. Single atoms are often easy to aggregate into nanoparticles (NPs), which are reported to be inactive in many reactions. Another problem that influences the stability of SACs is the leaching of metal atoms. In this case, the amount of the active site reduces, leading to low activity and selectivity. Thus, stable SACs mean that single atoms stay on the surface without desorbing and diffusing.
Based on DFT calculations, both calculated thermodynamic and kinetic properties are useful in ascertaining the stability of SACs, and overall consideration is necessary [26]. Thermodynamically, the stability of SACs is evaluated by energetics. The ability of SACs against aggregation can be evaluated by comparing the adsorption energy of the metal atom (E ad ) on the support and the bulk cohesive energy (E coh ) defined by equations (1) and (2), respectively [27].
where n is the number of atoms in the bulk; E support+atom , E support , E free-atom , and E bulk are the energies of the support anchored with a metal atom, the pristine support, the isolated metal atom, and the bulk, respectively. If the adsorption energy is smaller than the bulk cohesive energy, SACs are stable against aggregation and vice versa. With the calculated value of adsorption energy, the ability of SACs against metal being leached can also be facilely evaluated, a stable SAC exhibits negative adsorption energy value. Moreover, the ability of SACs against aggregation can be described by Gibbs free energy. As shown in Figure 3(a), SACs with smaller Gibbs free energy than NPs are stable. For SACs, the chemical potential can be calculated by: where μ SAC and E support+atom represent the chemical potential and energy of SAC, respectively. E support represents the energy of the pristine support and E bulk represents the energy of one atom in bulk. For NPs, the chemical potentials can be approximated by the formation energy calculated from: where μ NP , γ, and R represent the chemical potential, the surface energy, and the curvature radius of the NPs, respectively. Ω represents the molar volume of one atom in bulk. For instance, the calculated chemical potential of Au single atoms and NPs on different supports are plotted in Figure 3(b), clearly showing the stability of single atoms compared with NPs [28]. In addition, by combined chemical potential calculation and electronic structure analysis, the stability of supported Fe SACs was investigated [29]. The doped N could promote charge transfer and create a different coordinate environment that improves the thermodynamic stability of these SACs. On the other hand, dynamically, energy barriers are used to evaluate the stability of SACs, as shown in Figure 3(a) [28,30]. If the energy barrier of single atoms aggregating to NPs is high, the SAC is thought to be stable against aggregation. Besides, a high enough energy barrier of the single atoms diffusing away from another site suggests that they are stable against being leached. Based on the calculated static energy barriers, the temperature (T)-related ability of self-diffusion of adatoms can also be evaluated. For example, for in-channel diffusion, the jump rate r j and diffusion coefficient D can be calculated by: D ¼ 1 2 a 2 r j X l l 2 p l (6) wherer 0 j represents the prefactor, ΔE j represents the activation barrier, a represents the distance between two adjacent sites and p l represents the probability of an l-sites jump. Based on these parameters, the ability of selfdiffusivity of adatoms on Cu(110), Ag(110) and Au(110) was investigated. The calculated results show that at 450 K, self-diffusion on Cu(110) is activated, while on Ag(110) and Au(110) are similar [31]. At high temperatures, the percentage of long-jumps on Cu(110) and Ag(110) is increased, which can be explained by the energy dissipation rate of the adatom and the multidimensional topology of the potential energy surface [31].
Especially, for electrochemical reactions, the ability of SACs against leaching can be determined by the dissolution potential [32]. The leaching equation for a supported n-atom cluster (M n /substrate) can be written as: The corresponding reaction-free energy at potential U is: where G(M 2+ ,aq) can be calculated from the chemical potential of bulk metals G(M,s) and experimental standard redox potential U 0 : when ΔG<0, metal leaching proceeds spontaneously, thus the dissolution potential U n for leaching one metal atom can be estimated as: Based on this method, Liu et al. [33] calculated U n for NC, Fe 2 O 3 , and graphdiyne supported Pt n and Ni n (n =1-20) clusters as well as NC and Fe 2 O 3 supported Pd n (n =1-20) clusters. The highest leaching equilibrium potential before the formation of atomically dispersed structure and the lowest equilibrium potential for destroying the structure can be used to define the boundaries of the electrochemical potential window within which the SAC is stable against metal leaching. Similarly, the electrochemical potential window for graphdiyne supported Pt 3 and Ni 3 single cluster catalysts (SCCs) are also plotted. This approach of constructing electrochemical potential windows may also be extended to other supported metal clusters, providing a strategy for the selection of stable SACs and SCCs.

Activity
High catalytic activity and selectivity are the ultimate goals for designing SACs. Guided by DFT calculations, these catalytic properties for several reactions can be predicted, including widely studied reactions such as CO oxidation [30], water gas shift [34], selective hydrogenation [35], half reactions in the fuel cell [27], electrochemical nitrogen [14] and CO 2 reduction [36] and other reactions such as aromatic nitroreduction [37], alkane dehydrogenation [38], selective dehalogenation [39] and so on. We first introduce computational methods for the evaluation of the catalytic activity.
As examples, we mainly introduce the industrial widely applied gas-solid phase reactions and electrochemical reactions.

Gas-solid phase reactions.
To theoretically evaluate the activity and selectivity of SACs, unraveling the reaction mechanisms is a crucial job. For gas-solid phase reactions, this process is accomplished by selecting the initial and final stable structures for each elementary step and the transition states between them as well as getting corresponding energies for these structures. With these energies, activation barriers can be calculated, and reaction cycles during the reaction process can be easily plotted. As shown in Figure 4, CO oxidation goes through a Mars-van Krevelen (Mvk) mechanism on Pt 1 /CeO 2 . The oxygen atom of the support participates in the catalytic cycle, after the formation of two CO 2 molecules, the support is reoxidized by an adsorbed O 2 . The potential energies of the corresponding structures are plotted within one diagram, showing the energy changes and reaction barriers along the reaction process.
where k B represents the Boltzmann constant, T represents the absolute temperature, E a represents the activation barrier, q represents partition functions and v i represents the vibrational frequency. Through this method, kinetic models for several gas-phase reactions have been established, and different SACs have been selected as potential catalysts, such as Au/TiO 2 [40], Ag/MgO, and Ag/Fs-defect MgO [45] for CO oxidation, and Pt/Mo 2 C [43] and Pt/TiO 2 [42] for water-gas shift reaction.

Electrochemical reactions.
For electrochemical reactions, the activity of SACs is determined by the potential-determining step (PDS), which is defined as the step with the most maximum value of free energy change through the reaction process [47]. The PDS is the thermodynamic bottleneck for a reaction and the degree of difficulty of this step well describes the activity of the SACs [48,49]. The corresponding potential of the PDS is called the limiting potential (U L ), defined as equation (14): where ΔG PDS is the free energy change of the PDS. In the case of HER, the free energy of H adsorption (ΔG H� ) is used to evaluate the activity of SACs [50][51][52]. An effective SAC for HER should have ΔG H� close to zero to lower the reaction barriers of hydrogen adsorption and desorption steps [27]. For example, by investigating ΔG H� on a series of SACs, as shown in Figure 5(a), Hossain et al. [53] proposed Co@N-doped graphene as a promising SAC for HER with the ΔG H� value close to zero, which is proved by their further theoretical analysis and experimental results. For multi-step electrochemical reactions with several intermediates, free energy diagrams along reaction pathways are plotted, as shown in Figure 5(b) for oxygen evolution reaction (OER). A SAC with lower free energy change of the PDS is proved to be able to exhibit higher activity in the reaction. Based on these analyses, the PDSs on different SACs are highlighted in Figure 5(b) and Ni embedded in nitrogen-doped holey graphene frameworks stands out with the lowest free energy change (0.42 eV) for the PDS [7]. Besides, with the value of the U L , overpotential defined as the difference between it and the equilibrium potential for an ideal catalyst can also be used to evaluate the activity of SACs for electrochemical reactions [54]. A smaller overpotential indicates higher activity of the corresponding SAC.

Selectivity
By plotting reaction cycles and free energy diagrams, the selectivity of SACs can also be readily determined. Generally, for competitive steps, the reaction would go through the step with a lower activation energy or free energy barrier. As shown in Figure 6, on Mo-based MOF, CO 2 reduction reaction (CO 2 RR) chooses the way with a lower free energy barrier for each step and results in the selective production of CH 4 . As adsorption is the first step for surface reactions, the difference between the adsorption free energy of the main adsorbates can also be used as a descriptor for evaluating their selectivity [55,56]. For instance, to reduce the competition of HER in nitrogen reduction reaction (NRR), the adsorption of H at the active sites should be suppressed. Therefore, by comparing the adsorption energy of N 2 and H, the defective graphene-based SACs such as Mo@C 3 , Nb@C 3 , V@C 3 , and Mn@N 3 have been proved to show high selectivity with smaller adsorption energy of N 2 [55]. This is in accordance with the ensemble effect in SACs that H can only be adsorbed at the top site of the single atoms [55]. Besides, Bader charge analysis and charge density difference show that the Ti in Ti@N 4 is positively charged, which also contributes to the reduced H adsorption by electrostatic repulsion [55].
On considering the whole reaction processes, the step that determines the kinetic activity of SACs in gas-solid phase reactions or PDSs in electrochemical reactions is used to evaluate the selectivity of SACs. The reaction process with lower reaction barriers will be selectively chosen. For example, in the gas-solid phase selective hydrogenation of 3-nitrostyrene, only the H atoms at the top site of the active atom can be added to substrates; thus, the energy barrier of H diffusion from the fcc-hollow site to the top site (denoted as E a-react ) is used as a descriptor for the catalytic activity [57]. Subsequently, the selectivity of hydrogenation of 3-nitrostyrene on Ni(111), Pt 1 /Ni(111), and Pt(111) was investigated. The selectivity of Pt(111) is low because the nitro groups of 3-nitrostyrene molecule interact with the support too weakly to be activated. While on Ni(111), a high E a-react value (above 0.5 eV) indicates that the hydrogenation is hard to proceed. Only the Pt 1 /Ni(111) SAC can selectivity hydrogenate 3-nitrostyrene with a low E a-react value (0.17 eV). Besides, in the electrochemical reduction of oxygen, the two-electron transfer oxygen reduction reaction (ORR) toward H 2 O 2 competes with four-electron transfer ORR toward H 2 O. At 0.6 V, the free energy change of PDSs at each active site on Fe@N doped graphene for fourelectron transfer ORR are negative, much lower than the free energy change of PDSs for the two-electron transfer processes, which are all positive [58]. Therefore, by comparing the PDSs, Fe@N doped graphene is proved to be selective for H 2 O production with four-electron transfer ORR. Moreover, for electrochemical reactions in a protic solvent, HER is easy to become a side reaction. In this case, SACs with lower free energy barriers for the main reactions than HER are thought to be selective [59].

Computational design of SACs
Though of fine accuracy, the predictive results given by DFT calculations in Section 2 usually cannot be used directly for the rational design of SACs, as they do not reflect the intrinsic properties of SACs. However, based on these calculated catalytic performances and electronic analysis, effective descriptors can be identified. The descriptors can act as bridges between intrinsic properties and the catalytic performance of SACs. Starting from the desired catalytic performance, the corresponding value ranges of descriptors are determined, pointing out favored intrinsic properties of effective SACs. Identifying easily calculated and accurate descriptors thus opens an approach for the high-throughput screening and computational design of SACs.

Usage of energy descriptors
In some reactions, the catalytic performances of SACs have been described by a series of energy-related parameters, which are called the energy descriptors [14,61]. The identified energy descriptors are much easier to be figured out compared to the whole reaction process, facilitating the screening of effective SACs. Besides, the correspondence between energy descriptors and catalytic performance usually can be explained by electronic factors, which can help to understand the catalytic mechanisms and prove the reliability of identified energy descriptors.
An effective way of identifying energy descriptors starts from the Sabatier principle that an efficient catalyst should bind key reaction intermediates moderately. Thus, the adsorption energies or adsorption Gibbs free energies (ΔG ad ) of key intermediates can be used as predictive descriptors. In general, for reactions containing more than one key intermediates, there may be several descriptors. At this point, the scaling relationships between their adsorption energies or adsorption Gibbs free energies provide guidance in limiting the number of independent variables and thus simplifies the descriptors. In fact, as the variations in adsorption energies and transition state energies follow the same basic physics, transition state energies also correlate with adsorption energies. Thus, scaling relations called Brønsted-Evans-Polanyi (BEP) relations exist between activation (free) energies and reaction (free) energies, which is exactly what Sabatier principle describes [54,61]. Following these relations, several predictive descriptors can be identified for different reactions, and representative examples are demonstrated as follows.
In ORR and OER, scaling relations exist between adsorption Gibbs free energies of reaction intermediates and OH, as shown in Figure 7(a), which means that Sabatier principle works for these systems [62]. The change in adsorption Gibbs free energies is explained by the d-band center. It means that the higher the ε d , the stronger the adsorption energy of intermediates [54]. Through the reaction equation of each step, reaction free energies can be expressed by the adsorption free energy of intermediates, which are in correlation with OH adsorption Gibbs free energy. Correlating with the overpotential for ORR, a contour plot with OH adsorption Gibbs free energy as energy descriptor is discovered as Figure 7(b), and candidates can then be easily chosen [54]. Similarly, for NRR, scaling relations between adsorption energies of *NNH, *NNH 2 , *NH, *NH 2 , and *N (ΔE NNH � , ΔE NNH 2 � , ΔE NH � , ΔE NH 2 � , and ΔE N � ) are discovered [61]. The change in ΔE N � is explained by projected crystal orbital Hamilton population (pCOHP) and integrated COHP (ICOHP) analyses, and the value of ΔE N � is used as the rate-limiting steps [61]. As all the rate-limiting steps are related with *NNH and *NH 2 , their Gibbs free energy (ΔG NH 2 � and ΔG NNH � ) are discovered as energy descriptors for the limiting potential, and contour plots are plotted.
Apart from using Sabatier analysis and scaling relations, several energy descriptors have been discovered through the analysis of reaction mechanisms. For instance, all the reaction mechanisms for NRR contains two 'stable to instable' transitions: *N 2 + H + + e − = *N 2 H and *NH 2 + H + + e − = *NH 3 [14]. Besides, the first and final steps of the chemisorption of N 2 and desorption of produced NH 3 , respectively, are also crucial. Hence, the free energy change of the two key elementary steps (ΔG N 2 À N 2 H and ΔG NH 2 À NH 3 ), N 2 adsorption energy (ΔE N 2 ) and NH 3 desorption free energy (ΔG NH 3 À des ) are identified as energy descriptors in NRR [14]. Another example is the identification of energy descriptors for CO oxidation on rutile TiO 2 (110) supported SACs. The reaction goes through MvK mechanism, and only the energy changes of CO reaction with the surface oxygen atom and O 2 adsorption at the vacancy site vary with different dopants. Additionally, the calculated energy change of the former step is a linear function of ΔE V while that of the latter step decays linearly with ΔE V , as shown in Figure 7(c, d). Therefore, the energy of oxygen vacancy formation is identified as an energy descriptor [63].

Utilizing structure descriptors
Energy descriptors have been proved to be successful for the rational design of SACs. However, they do not include the intrinsic structural factors determining catalytic properties. Hence, it is usually hard to directly tune the values of energy descriptors such as adsorption energies of SACs to engineer their performance to the top of the volcano plot. In this sense, aided by the analysis of electronic structures and other intrinsic properties, structure descriptors relating the catalytic properties with the intrinsic structural properties have been identified. In our previous work, a universal structure descriptor that can describe the activity for several reaction systems was presented for ORR, OER, and HER reactions by the graphene-based SACs [27]. A total of 112 graphenebased SACs were constructed and calculated [27]. Though ΔG OH � and ΔG H � correlate well with ORR, OER, and HER activity, they do not reflect the intrinsic characteristics of the active center. Based on the analysis of the electronic structure of the active sites, the structural parameter φ, which is defined as N and, C, respectively; n N and n C are the number of nearest N and C atoms and α is correction coefficient), was identified as a structure descriptor. As plotted in Figure 8(a,b), SACs owning a φ value of 30-32, 30-35, and around 17 or 39 are prone to exhibit higher activity than precious-metal-based commercial catalysts for ORR, OER, and HER, respectively [27]. This work opens an avenue for identifying universal structure descriptors in the usage of predicting the performance of catalysts for several reactions through one readily obtained parameter. In another work, structure descriptor Φ =V M ×E M /r M (V M , E M and, r M represent the valence electron number of metals, the electronegativity of central metals, and the radius of metal ions, respectively) was discovered which well correlates with the reaction overpotential for M-N-C SACs for CO 2 RR [64]. As the volcano relationship in Figure 8(c) shows, in the regime of Φ <5.0 the reaction favors HER while in other regimes favors CO 2 RR. Besides, in different regimes, the reaction goes through different pathways labeled η 1 (C) (CO2 → *COO → *COOH → *CO → CO) or η 2 [64]. These suggest that the descriptor Φ is an effective descriptor that can well distinguish different reaction mechanisms, intermediates, reaction pathways, and final products. In addition, the Bader charge of the metal ions is in relation to descriptor Φ, as shown in Figure 8(d), indicating its power to reflect intrinsic factors. It is expected that Φ can be used as a universal structure descriptor of M-N-C SACs for wider applications.
Other attempts have also been made to the rational design of SACs through the construction of structure descriptors. Based on the calculated magnetic moment and adsorption energy as well as the analysis of electronic structure, 0.62X-0.21ε d -5.04 was identified as a structure descriptor for O adsorption energy on graphene-based SACs [65], where X represents the electronegativity and ε d represents the d-band center of the catalytic surface. This descriptor can be used for the high-throughput screening of SACs [66].

Machine-learning aided catalyst design
Utilizing descriptors has indeed promoted the rational design of SACs. However, up to now, universal descriptors that are suitable for various types of SACs and reaction systems are still insufficient. As a result, a vast amount of DFT calculations still need to be carried out to screen different descriptors for each reaction system. At this point, machine-learning based methods can help to break the bottleneck [67]. Based on DFT calculations providing basic data, statistical learning could directly point out the accurate correlations between the intrinsic factors and catalytic performance, giving effective function forms to predict the catalytic properties. Recently, successful attempts by combined DFT and statistical learning or machinelearning methods have been reported to study SACs [68].
An example is to identify different descriptors that can predict interaction strengths between single atoms and the oxide support by combined DFT calculations and statistical learning approach [68]. Through DFT calculations, the oxide formation enthalpy (ΔH f, ox ) of metal adatom and the oxygen vacancy formation energy (ΔE vac ) of support correlate well with the binding energy of single metal atoms, which can be explained by the electronic interactions as introduced in Section 2.1 [68]. In the following, the feature space was expanded to containing 333,932 descriptors that include parameters of the atomic properties of the adatom and support. Through these descriptors, binding energy was predicted, and their comparison to DFT calculated data is shown in Figure 9(a,b). By LASSO+l 0 analysis, the robustnesses of the descriptors are examined, and results show that the most accurate descriptors giving the lowest RMSE are made up of |ΔH f, ox /ΔE vac | multiplied by a parameter describing the property of surface metal or together with the adatom. This work explained the trends by intrinsic factors and provided a successful model for screening metal/support combinations that produce stable SACs against sintering. Besides, machine-learning has been used to predict the activity and selectivity of SACs. Guo et al. [66] employed a feature set containing 8 intrinsic descriptors by machine-learning and found that the oxide formation enthalpy (H f, ox ) and the adjusted electron numbers of d/p orbital (e dp' ) show the highest feature importance of 0.59 and 0.20, respectively, in determining ΔG (O*) thus selectivity of certain two-dimensional SACs, as shown in Figure 9(c). Their predicted ΔG (O*) data by machine-learning versus DFT calculated data are plotted in Figure 9(d) [66]. In addition, Lin et al. [69] developed three models for ORR, OER, and HER, respectively, using machine-learning, which successfully predicted their limiting potentials on graphene-supported SACs. Similarly, Fung et al. [70] identified equation (15) to accurately predict ΔG H on M-N-C-based SACs, which is a general descriptor for HER activity.

Summary and outlook
In conclusion, SACs emerge as a new candidate that can cut down noble metal consumption in catalysis. They have shown exceptional performance for many important reactions, and their relatively simple structure makes them ideal models to discover the underlying mechanisms for catalysis. To reduce the time and economical-consuming work of complicated experiments, computational efforts have been paid to investigate the SACs. Density functional theory (DFT) calculations can reveal the properties of the active sites, and help to understand the exceptional performance of SACs. In addition, the catalytic performance including stability, activity, and selectivity of SACs can also be predicted by DFT calculations. Here, based on these DFT calculated information, we highlight the identification and usage of energy and structure descriptors to describe the performance of SACs which can guide the high-throughput screening of SACs. Energy descriptors are identified as a series of energy-related parameters that can well reflect the DFT calculated catalytic performance of SACs, which are more readily to be figured out than the whole reaction processes. Compared with energy descriptors, structure descriptors can directly relate the performance of SACs with intrinsic structural factors, making the properties of active sites more explicit and thus directly guide the rational design of SACs. Moreover, the advanced machine-learning methods have also been used for the identification of effective descriptors. Through machine-learning, effective functional forms can be readily constructed consisting of parameters selected from a large feature set. With the development of computational techniques, we are expecting the combined DFT and machine-learning methods will provide us with general strategies for identifying universal and effective structure descriptors to guide the rational design of SACs. Apart from the effective usage in the design of SACs, computational methods are also expected to be of great help to the research progress of other atomic-scale dispersed catalysts, such as double atom catalysts (DACs), single cluster catalysts (SCCs), and so on [71,72]. For example, Guo et al. [73] performed large-scale DFT calculations and obtained an activity map for NRR on bimetal catalysts supported by two-dimensional expanded phthalocyanine. The adsorption energy of N 2 H* was identified as a descriptor in the map, which has been used for the screening of promising bimetal catalysts. Besides, first-principles calculations have been applied for the demonstration of in situ dynamic generation of multinuclear sites from single atoms implicated as active sites in O 2 activation [74]. The dynamic formation of multinuclear sites is discovered as distinct from traditional heterogeneous catalysts, providing a more thorough understanding of the mechanisms of catalysis. With the rapid development of computational skills, it is expected that universal descriptors and computational methods will open new avenues for the rational design of catalysts in wider applications.