Zhaoming
Xia
,
Yue
Yin
,
Jun
Li
and
Hai
Xiao
*
Department of Chemistry and Key Laboratory of Organic Optoelectronics and Molecular Engineering of the Ministry of Education, Tsinghua University, Beijing 100084, China. E-mail: haixiao@tsinghua.edu.cn
First published on 1st February 2023
Owing to limited degrees of freedom, the active sites of stable single-atom catalyst (SAC) often have one structure that is energetically much lower than other local-minimum structures. Thus, the SAC adopts one lowest-energy structure (LES) with an overwhelmingly larger proportion than any other high-energy metastable structure (HEMS), and the LES is commonly assumed to be solely responsible for the catalytic performance of an SAC. Herein, we demonstrate with SACs anchored on CeO2 that the HEMS of an SAC, even though its proportion remains several orders of magnitude lower than the LES throughout the catalytic reaction, can dictate catalysis with extraordinary activity arising from its unique coordination environment and oxidation state. Thus, we unravel the key role of HEMS-enabled catalysis in single-atom catalysis, which shakes the common assumption in the studies of SACs and urges new developments in both experiment and theory to identify and exploit catalysis via HEMSs.
An intricate yet fundamental aspect is the dynamic structural evolution of active sites coupled with catalytic mechanisms,9 and this has been well exemplified by theoretical studies of catalysis on supported clusters and nanoparticles.7,8,10,11 The active sites of these catalysts are composed of multiple atoms and thus possess large degrees of freedom, which result in ensembles of local-minimum structures (LMSs) that are energetically close12 and thus accessible to catalysis, based on the Boltzmann distribution.8 In addition, the LMSs in the accessible ensembles can vary with the reaction atmosphere, which is induced by adsorptions.13,14 These impose challenging complications on deciphering the structure–performance relationships of these catalysts.
Nevertheless, as the number of LMSs exponentially decreases with the system size,12 the scenario becomes intrinsically different when the size of the active site is pushed to the limit of a single atom, i.e., the concept of a single-atom catalyst (SAC). Since it was first proposed in 2011,15 the SAC has developed into a promising design concept for highly efficient heterogeneous catalysts,16–19 and it also serves as an ideal model for the study of the structure–performance relationship of catalytic active sites, owing to its well-defined local structure with precise tunability.16,20 This exactly arises from the fact that because of the limited degrees of freedom, the active sites of stable SACs often present one lowest-energy structure (LES) dominant in quantity, and the other local minima are high-energy metastable structures (HEMSs) in negligibly low concentrations, as illustrated in Fig. 1.
Thus, it is a sensible and common practice for both experiment and theory to assume that the LES of an SAC is solely responsible for its catalytic performance, and a clear structure–performance relationship can be established. Yet, it has been reported that SACs can reversibly switch between different structures, including cluster forms and varying coordinations,9,21,22 in which the different structures are energetically close and thus are accessible with detectable concentrations, while the HEMSs remain reasonably ignored.
Herein, we found that the HEMSs of SACs can significantly contribute to catalysis, even though they remain at negligibly low concentrations throughout the catalytic cycle, which is thus often overlooked in both experiment and theory. By investigating the local minimum structures and phase diagrams of transition metal M (M = Fe, Co, Ni, Cu, Ru, Rh, Pd, Ag, Os, Ir, Pt, Au) SACs (M1) anchored on the CeO2 (110) surface (M1/CeO2) with density functional theory (DFT) calculations, we identified that under a reductive environment, both Au1/CeO2 and Pt1/CeO2 possess a HEMS of the same type that is energetically higher than the LES by over 0.7 eV and is thus with a negligibly low concentration less than 10−5. However, we demonstrated that this HEMS of Au1/CeO2 and Pt1/CeO2 SACs plays a key role in delivering the catalysis of reduction reactions by its much higher activity arising from the unique oxidation state of M1, and it is thus the HEMS that enables the catalytic performance of SACs for these reactions.
The thresholds for converging the electronic wavefunctions and atomic structures are 10−6 eV for the energy and 0.02 eV Å−1 for the maximal force, respectively. The calculated lattice constant (5.418 Å) and Ce–O bond (2.360 Å) length of the CeO2 bulk are consistent with experimental results.32 The CeO2 (110) surface was then modeled by a p(3 × 2) 5-layer supercell with the bottom 3 layers fixed, and a 15 Å vacuum gap was included to separate the periodic images along the z-axis. A 2 × 2 × 1 Monkhorst–Pack grid was used to sample the Brillouin zone. Dipole correction along the z-axis was applied.
The Gibbs free energy corrections (including enthalpic and entropic terms) of gas-phase molecules and surface-adsorbed species were calculated with the ideal gas approximation and the harmonic approximation for the vibrational degrees of freedom using the thermochemistry modules of the Atomic Simulation Environment (ASE).33
Comparisons between electronic affinities with and without relativistic effects were obtained using the ADF program,34 which uses all-electron Slater-type orbitals as basis sets. The relativistic effects are included in the Hamiltonian via the zero-order regular approximation (ZORA).35
Microkinetic modeling was conducted using CatMAP36 to explicitly investigate the catalytic kinetics under varying conditions. The partial pressures of all the gas phase reactants were set to 1 bar, while the partial pressures of the gas phase products were set to 0.1 bar for C2H2 and 10−15 bar for HCOOH to ensure overall exothermic reactions with positive net forward reaction rates. Note that the extremely low partial pressure of HCOOH is in accordance with the experimental conditions (see the note in the ESI†).
The most stable phase, i.e., the LES of the active site, on each M1/CeO2 SAC widely changes with varying atmosphere. In particular, the coordination number of M to O (CNO) varies among different LESs. In the reductive condition represented by the pressure of H2 (pH2), the phase diagrams (Fig. S11 and S12†) indicate that all the M1/CeO2 SACs retain a particular type of LES with a CNO of 2 (labeled as the MsubO-2 type) in a wide temperature range at pH2 of 1 bar.
In the LES of the MsubO-2 type, the single-atom M is anchored at the O vacancy site and bonds to the two nearest O atoms to form an O–M–O configuration (Fig. 2a). The M atom can lose one of the M–O bonds to form an M–O configuration, and this is labeled as the MsubO-1 type (Fig. 2a). Only with M = Cu, Pd, Pt, and Au, can the M1/CeO2 SAC adopt MsubO-1 as a LMS.
![]() | ||
Fig. 2 (a) The MsubO-1 and MsubO-2 types of structures for active sites in the M1/CeO2 SACs. The two are well separated in energy by more than 0.7 eV, with the MsubO-2 type being the LES and the MsubO-1 type being the HEMS under reductive conditions. The energy gap results in the extremely low concentrations of (b) PtsubO-1 and (c) AusubO-1 in the phase diagrams, while the concentrations of (d) PtsubO-2 and (e) AusubO-2 dominate the phase diagrams of active sites in the corresponding M1/CeO2 SACs under varying temperatures and pH2 (in logarithmic scale). Note that the number of atoms for MsubO-1 and MsubO-2 are the same, while the dependence of the phase diagram on pH2 arises from the competition among all 8 types of SAC structures considered in Fig. S2–S8,† for which the numbers of O atoms are different. |
Under the reducing environment with the standard pH2, MsubO-2 is the LES, and MsubO-1 is the HEMS that is over 0.7 eV higher in energy (Fig. 2a and Table S1†). Therefore, the concentrations of HEMS (MsubO-1) are less than 10−5 in the phase diagrams of the active sites of M1/CeO2 SACs (Fig. 2b and c), and the active sites dominantly adopt the LES (MsubO-2) (Fig. 2d and e). Conventionally, the catalytic performance by such a HEMS has often been overlooked for further investigation. However, this HEMS (MsubO-1) has an unusually low CNO of only 1 (Fig. 2a), and such a highly unsaturated coordination mode of the single-atom M may spark extraordinarily high catalytic activity. Therefore, we further scrutinized the role of this HEMS of the MsubO-1 type in catalysis, in comparison with that of the LES of the MsubO-2 type.
Fig. 3a shows the free energy profiles of the reaction pathways for CO2HR catalyzed by the LES (L-pathway) and the HEMS (H-pathway) of the Au1/CeO2 SAC under standard conditions, with structures of all states illustrated in Fig. 3b. In the L-pathway, starting from the LES of the MsubO-2 type (L0), the CO2 molecule prefers to bind with the surface O atom to form a carbonate group (L1, Fig. 3c) instead of weak physisorption on the positively charged Au1 site of L0 (Fig. S13†). H2 then easily dissociates into H+ and H− on the Lewis pair composed of the carbonate group and Au1 site (L2) with a small barrier of 0.58 eV (L-TS1). Subsequently, H− attacks the carbonate group to produce HCOOH (L3), which, however, poses a formidable barrier of 3.31 eV (L-TS2).
In contrast to the L-pathway, CO2 in the H-pathway can spontaneously and directly bind to the Au1 site in the HEMS of the MsubO-1 type (H1, Fig. 3d) with a bent configuration, indicating a chemisorbed CO2. It is subsequently hydrogenated by H2via an Eley–Rideal (ER) mechanism with a barrier of 1.38 eV (H-TS1) to form the co-adsorbed COOH and H on Au1 (H2), which combine to produce HCOOH (H3) by overcoming a small barrier of 0.31 eV (H-TS2). Upon releasing HCOOH, the Au1 site transforms to L0, which recovers H0 by overcoming a barrier of 1.15 eV.
Fig. 3c and d highlights the key difference between the CO2 activation modes driven by the two types of active sites, i.e., the LES of the MsubO-2 type activates CO2 on its O site to form a carbonate group, while the HEMS of the MsubO-1 type chemisorbs CO2 directly on its Au1 site to form a bent CO2.
In addition, we investigated the direct transformation of L0 to H1 upon the adsorption of CO2 on Au. This, however, requires the C atom of CO2 to approach the positively charged Au1 site of L0, which leads to a large repulsion, and it is more favored for the C atom of CO2 to bind to the surface O to form L1. Thus, the transformation from L0 to H1 follows a detour through L1 with a barrier of 1.02 eV (LH-TS0 in Fig. 3a) referenced to L0, and this forks to a new pathway labelled as the LH-pathway.
The three pathways, in fact, compose a complex reaction network via common states, and thus, they all contribute to the kinetics of catalysis. Therefore, we performed microkinetic modeling36 of this complex reaction network. The results (Fig. 3e) show that the H-pathway dominates the contribution to the TOF, and is larger than those of the L- and LH-pathways by several orders of magnitude. Nevertheless, Fig. 3f shows that the coverages of states in the H-pathway (H0–H3) are extremely small (<10−19 at 300 K and <10−4 even at 1000 K), and thus, they remain marginal and barely detectable. In contrast, the state (L2 or L0) in the L-pathway dominates the Au1/CeO2 SAC surface at all temperatures considered. Therefore, it is in fact the HEMS with a miniscule concentration, instead of the LES, that drives the CO2HR on the Au1/CeO2 SAC, and the extraordinary catalytic ability of the HEMS arises from the unique oxidation state of its Au1 site to activate CO2, as discussed later.
Fig. 4a shows the free energy profiles of the reaction pathways for SHA catalyzed by the LES (L-pathway) and the HEMS (H-pathway) of the Pt1/CeO2 SAC under standard conditions, with the structures of all states illustrated in Fig. 4b. In the L-pathway, the LES of the MsubO-2 type (L0) prefers the physisorption of H2 (L1) to that of C2H2 and delivers the dissociative chemisorption of H2 (L2) by overcoming a barrier of 0.90 eV (L-TS1). Subsequently, C2H2 consumes the adsorbed H atoms (L3, Fig. 4c) through an ER mechanism with a small barrier of 0.28 eV (L-TS2) and produces C2H4.
In contrast to L0, H0 in the H-pathway favors strong bonding to C2H2 (H1, Fig. 4d), which is exothermic by 2.01 eV. Such strong adsorption of C2H2 can couple with the transformation of active sites and leads to a low barrier of 0.59 eV at 300 K (LH-TS1 in Fig. 4a). This forks to a new pathway labeled as the LH-pathway (Fig. 4a) that involves the coupling of C2H2 adsorption and L0–H1 site transformation. However, this low barrier rapidly increases with temperature due to the entropic contribution of gas phase C2H2, and it reaches 1.32 eV at 800 K (LH-TS1 in Fig. S15†). Thus, the transformation of the active site is decoupled from the chemisorption of C2H2, favoring the direct transformation from H0 to H1 at temperatures higher than 710 K (Fig. S15† shows the scenario at 800 K), and this recovers the H-pathway with the initial state of H0.
In the H- and LH-pathways, upon the chemisorption of C2H2 (H1), H2 then hydrogenates the adsorbed C2H2via an ER mechanism to produce C2H4 by sequentially overcoming two barriers of 0.93 eV (H-TS1) and 0.27 eV (H-TS2). Upon releasing C2H4, the H-pathway reaches L0, which recovers H0 via a barrier of 1.23 eV (TS0).
Fig. 4c and d highlights the key difference between the C2H2 activation modes driven by the two types of active sites. The LES of the MsubO-2 type cannot adsorb C2H2, and instead, it activates H2 to the adsorbed H atoms on its Pt and O sites, which then activate C2H2. In contrast, the HEMS of the MsubO-1 type can directly chemisorb C2H2 on its Pt1 site.
We further performed microkinetic modeling36 of the complex reaction network composed of all three pathways in Fig. 4a. The results in Fig. 4e show that the LH-pathway dominates the contribution to the TOF at low temperatures, and the H-pathway dominates the contribution to the TOF at high temperatures. Fig. 4f shows that the Pt1/CeO2 SAC surface is dominated by H1 at most of the temperatures because of the strong adsorption of C2H2 on the HEMS. Thus, the adsorption state derived from the HEMS is detectable in this scenario, but the HEMS itself (H0) retains a negligibly small coverage (<10−30 to 10−8). Nevertheless, the coverage of LES (L0) increases with temperature and then dominates the Pt1/CeO2 SAC surface at 1000 K. Therefore, the HEMS drives the SHA on the Pt1/CeO2 SAC at high temperatures, and LES drives the SHA at low temperatures. The extraordinary catalytic ability of the HEMS arises from the unique oxidation state of its Pt1 site to activate C2H2, as discussed below.
We investigated the OSs of M1's in the M1/CeO2 SACs. Fig. 5a shows the results of the Bader charge analysis,52,53 which indicates that the M1's in the LES (MsubO-2 type) of all M1/CeO2 SACs are positively charged, but the Pt1 and Au1 in the HEMS (MsubO-1 type) are negatively charged. The Hartree potential difference (ΔVH)54,55 in Fig. S17† also implies similar conclusions. Together with the spin population analysis (Fig. S18†), our results suggest that the OSs of Pt1 and Au1 in the HEMS are both −1.
We further calculated the 1s core level shifts of Pt1 and Au1 (Table S2†), and the results also indicated that the OSs of Pt1 and Au1 in the HEMS are −1. Although there are no clear-cut methods to determine the OS, various analyses provided strong evidence indicating that Pt1 and Au1 in the HEMS are negatively charged. This may be similar to the case of negatively charged Au SAC on CuO, which was experimentally observed by in situ scanning tunneling microscopy (STM) imaging and X-ray photoelectron spectroscopy (XPS) in a reducing environment.56
To understand why Pt and Au can reach negatively charged states, we calculated the electron affinities (EAs) of gas phase metal atoms with and without the relativistic effects (Fig. 5b). If no relativistic effects were taken into account, Pt and Au would have similar EAs to those of Pd and Ag, which were much less likely to take electrons from the support. With including the relativistic effects, Pt and Au possess the largest EAs among all metals. This arises from the contraction and thus stabilization of the 6s orbitals introduced by relativistic effects,57–60 which enables Pt and Au to hold extra electrons.
The negatively charged Au1 and Pt1 in the HEMSs serve as strong electron donors that can deliver facile activation of CO2 and C2H2. The differential electron density contour in Fig. 5c illustrates significant electron transfer from Au1 to CO2 that subsequently results in a bent CO2δ−. The bent CO2δ− forms strong bonding with Au s and d orbitals, as shown by the projected crystal orbital Hamilton population (pCOHP) analysis of Au–C in Fig. 5e, with an integrated COHP (ICOHP) of −3.4 eV. Consequently, the ICOHP of the C–O bond in CO2 increases from −18.4 eV to −13.7 eV upon adsorption, indicating the presence of greatly weakened C–O bonds.
Similarly, the differential charge density contour in Fig. 5d shows significant electron transfer from Pt1 to bonding with C2H2 that results in a bent C2H2, implying a transition from the sp hybridization of C to sp2 hybridization. The pCOHP analysis of Pt–C in Fig. 5d indicates strong bonding between C and the s and d orbitals of Pt, with an ICOHP of −2.8 eV. As a result, the C–C bond in C2H2 is weakened, as revealed by the increase in the ICOHP of C–C from −15.4 eV to −14.6 eV.
Our finding shakes the common assumption in the studies of SAC that the LES, which is much lower in energy by approximately 1 eV than any other HEMS, is responsible for the catalytic performance of the SAC. Therefore, this urges developments of new methodologies and techniques in experiment and theory to scrutinize the role of HEMS in the structure–performance relationship of SAC. Moreover, the HEMS in SACs may be further promoted and thus exploited with external fields, as in the scenarios of electrochemistry and photochemistry.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2sc06962h |
This journal is © The Royal Society of Chemistry 2023 |