Catalytic Role of Vacancy Diffusion in Ceria Supported Gold Catalyst

Dynamics of intrinsic defects are considered fundamental in the chemistry of reducible oxides, and their effect on catalytic reactions have been rarely reported. Herein, we propose a new Ov diffusion assisted Langmuir-Hinshelwood mechanism for CO oxidation, which may largely account for the origin of high reactivity of supported atomic gold catalysts.

][9][10][11][12][13][14][15] However, few studies are devoted to discussing their effect on the realistic catalytic reactions.This is primarily because of the complexity of the systems with multiple local minima with respect to the sites on which the excess electrons-driving the Ce 4+ (4f 0 ) to Ce 3+ (4f 1 ) reduction-localize. 8,9,16n this study, we chose ceria supported Au adatoms (monomer and dimer) as typical examples and carried out thorough densityfunctional theory investigations to explore the dynamics of intrinsic defects accompanying the redistribution of the 4f localized electronic states, and illustrated their key role in catalysing CO oxidation processes.
8][19][20] The high reactivity is generally attributed to the strong metal support interaction between ceria and gold particles. 21,22In particular, O v s, the most common defects in ceria, can engender excess charges to the system, which not only enhance the binding of gold particles, but also adjust the oxidation state of gold particles. 5,23owever, controversies still exist concerning the role of O v in gold catalysis.][26][27] On the other hand, recent studies showed that the subsurface O v was more stable than the surface O v and was the dominant species on reduced ceria surface. 28,29][32][33] However, Weststrate et al. 27 recently showed that highly dispersed gold adatoms on reduced ceria seem to bind molecular CO more weakly than oxidized ceria.Camellone et al. 5 demonstrated that the oxygen vacancy provides an inhibiting effect on the activity of supported gold catalysts but this effect weakens as the cluster size of gold increases.Therefore, an in-depth understanding of the interplay between gold adatoms and near-surface O v s is of fundamental importance to explore the nature of high reactivity of supported gold catalysts.
The calculations have been performed within the framework of DFT with generalized gradient approximation using the VASP code. 34,35][38] Previous studies indicated that multiple self-consistent solutions exist, corresponding to different occupations of the m projections associated with the f-shell to which the U parameter is applied. 39,40owever, in this study, the used Davidson-block iteration scheme with random initialization of the orbitals, for which the most extensive tests have been performed, is considerably robust 41 (see calculation details in ESI †).
To validate our methodology, the types and distributions of O v s in the bare ceria slab (Fig. S1 vacancy site itself.The calculated O v formation energy of SSSV (1.91 eV) is 0.19 eV lower than that of SSV (2.10 eV), suggesting that SSSV is energetically more stable than SSV (Table S1, ESI †).This is quantitatively consistent with previous calculations (0.18 eV). 9 In addition, our further examinations indicate that the vacancies tend to separate up to a distance equal to twice of the (1 Â 1) surface lattice parameter (Fig. S1 and Table S1, ESI †), i.e., the third-nearest neighbour in the subsurface layer, to form oxygen vacancy pairs, which is consistent with the recent experiments and calculations. 28,29rior to investigating the reaction mechanism of CO oxidation, we consider the interplay between gold adatoms and O v s (Fig. S3-S5 and Tables S3, S4, ESI †), which is generally believed to be important to promote the catalytic activity of gold.For a gold monomer on CeO 2 (111) support with one O v (SSV or SSSV), the most stable configurations are shown in Fig. S4a and c, ESI.† It is found that the gold monomer occupied SSV is much more stable than the SSSV with gold monomer binding to the surface oxygen, though SSSV is more stable without the presence of the gold monomer.Specifically, we find that the SSSV only needs to experience a small energy barrier of 0.27 eV to hop up from subsurface to surface site in the presence of a single gold adatom.During the hopping process, the number of excess electrons at Ce sites is reduced 3 to 1, suggesting that the charge state of gold adatom changes from +1 to À1 by attaining two excess electrons.These results are consistent with previous calculations. 4,5,12,23,42,43e further introduce another Au atom and an extra O v into the system to explore the effect of high concentration of vacancies on the behaviour of gold cluster.The most stable configurations for Au 2 @CeO 2 (111)-SS-SS, Au 2 @CeO 2 (111)-SS-S and Au 2 @CeO 2 (111)-S-S (SS means subsurface O v and S means surface O v ) are listed in Fig. 1.We find an interesting reconstruction process of Au dimer in the presence of O v pair.Initially, Au dimer is put on the surface bridged oxygen ions and the two O v are located in the subsurface layer.The system only has 4 excess electrons corresponding to the two O v and the charge state of Au dimer is roughly estimated to be neutral.It is shown that one of the two O v could easily hop up from the subsurface to the surface and be filled immediately by one Au atom.This process only needs an extremely small energy barrier of 0.04 eV (TS1 in Fig. 1).Simultaneously, it will reduce the oxidation state of the gold atom, which fills into the surface O v , from 0 to À1 but the other Au is oxidized from 0 to +1.Further hopping of the other O v with a barrier of 0.88 eV (see TS2 in Fig. 1) is relatively difficult, partially due to the strong Coulombic repulsion of two Au-ions when both the O v locate on the surface.These results strongly imply that under realistic conditions, not all the excess charges from O v s can have an impact on the gold particles due to the thermodynamically unfavorable Coulombic repulsion.Compared to the gold monomer, the presence of Au dimer cluster seems to improve the stability of subsurface O v .The energy difference for three distributions of O v shown in Fig. 1 is less than 0.2 eV.
Based on the most stable configurations of ceria-supported Au adatoms, we now consider the reactivity of CO oxidation to explore the potential effect of oxygen vacancy-metal interplay on the catalytic reactivity of reducible oxide supported metal particles.For Au monomer in O v , the adsorption of CO is very weak (E ad = À0.16eV) due to the full-occupied 6s orbitals of negatively charged Au, consistent with previous studies. 5Therefore, next, we only consider the reaction mechanisms of CO oxidation on ceria-supported gold dimer, where one O v locates at the subsurface layer and the other locates at the surface layer.We find a novel Langmuir-Hinshelwood (LH) mechanism assisted by the hopping of oxygen vacancy between the surface layer and subsurface layer, which exhibits significantly higher reactivity than the traditional LH mechanism.Note that the Mars-van Krevelen Mechanism was also tested for existence but did not prove to be prevalent in this study (Fig. S6, ESI †).
In the traditional L-H mechanism (see Fig. 2a), reaction only involves CO reacting with the adsorbed O 2 species and the O v s are not allowed to diffuse during the entire process.Initially, CO is adsorbed on the gold dimer with adsorption energy of À1.58 eV but O 2 molecule only weakly binds to one Au atom.The adsorption energy is À0.31 eV and the bond length of adsorbed O 2 is 1.24 Å, close to 1.23 Å in the gas phase, suggesting a physisorption of O 2 .Next, CO and O 2 move towards each other, forming an O-C-O-O* complex.This step experiences a barrier of 0.68 eV and is slightly endothermic by 0.12 eV.Subsequent formation of a gas-phase CO 2 is highly exothermic (À2.17 eV).The residual O* can directly react with an additional gas-phase CO to form CO 2 molecule, completing the reaction cycle.
In contrast, Fig. 2b shows a dynamic reaction mechanism of CO oxidation assisted by the oxygen vacancy diffusion.CO initially adsorbs to the Au dimer with high adsorption energy of À1.58 eV.Distinctly different from the traditional L-H mechanism, O 2 adsorption is significantly enhanced by 0.31 eV, when the surface O v diffuses into the subsurface by overcoming a small barrier of 0.29 eV.The bond length of the adsorbed O 2 is increased to 1.33 Å, indicating the formation of superoxide species (O 2 -).After that, the adsorbed CO and O 2 only need to experience a barrier of 0.15 eV to form an O-C-O-O* complex.Subsequent reactions After the catalytic cycle is completed, one O v will hop back to the surface layer, charging the Au cluster and stabilizing it.Therefore, the proposed O v diffusion assisted L-H mechanism is expected to exhibit much higher reactivity than the traditional mechanism, which may largely contribute to the high reactivity of ceria supported gold particles at low temperatures.
To understand the nature of the enhancement effect from O v diffusion, we first monitor the charge state of Au dimer during the catalytic cycle.It is shown that the migration of O v from surface to subsurface leads to the re-oxidation of negatively charged gold adatom, which initially occupies the O v site, implying that the strong interplay between ceria support and the gold adatoms may occur upon O v diffusion.We thus further analyse the change of spin density isosurface and the projected density of states (PDOS) during the O v diffusion processes, as shown in Fig. 3.In Fig. 3a and b, it is found that one additional 4f electron is localized on a Ce atom when the surface O v migrates to the subsurface, indicating electron transfer from gold atom to ceria.This is further confirmed by the increase of 4f states in Fig. 3c and d.Before the diffusion, the energy levels of unoccupied O 2 2p* states are higher than both saturated Au 6s states and occupied 4f states, and the O 2 is thus not able to be activated.When the O v migrates into the subsurface, the occupied Au 6s states become strongly correlated with the lattice O 2p (valence band) and as a result, the energy level of 6s gets lifted, which can release excess electrons into the Ce 4f states and unoccupied O 2 2p* states.Moreover, the resultant Au 6s empty states can interact with the occupied O 2 2p states, leading to a down shift of these states and stabilizing the adsorbed O 2 .In other words, the O v diffusion essentially activates the Au 6s states so that they can correlate well with the frontier states of adsorbed intermediates such as O 2 * and O-C-O-O* complex and activate them by charge transfer.
In summary, we have presented a thorough understanding of charge distributions resulting from O v s and proposed a newly dynamic L-H reaction mechanism for CO oxidation on ceria-supported Au dimer.The diffusion of O v between the surface and subsurface layer is shown to promote the reactivity of CO oxidation during the catalytic cycle and stabilize the Au dimer after the cycle is completed.In general, this study strongly suggests that the interplay between metal particle and the intrinsic defect of the support is vital in understating the origin of high reactivity of oxide supported gold catalysis.
, S2 and Tables S1, S2, ESI †) were further examined.The most stable configurations for single surface O v (SSV) and single subsurface O v (SSSV) on CeO 2 (111) surface contain two excess electrons both located at the next-nearest neighbour cerium positions relative to the This journal is © The Royal of Chemistry 2017

Fig. 1
Fig. 1 Calculated path and corresponding energetics for the migration of O v on Au 2 @CeO 2 (111)-SS(S)-SS(S)-N.In CeO 2 (111)-SS(S)-SS(S)-N, SS(S) means subsurface (surface) O v and N means they are in the nth nearest neighbour position.The spin density isosurface is shown in dark area.Cerium, surface oxygen, and subsurface oxygen atoms are depicted in white, red, and pink, respectively.Open circles represent oxygen vacancies.For clarity, only the first OÀCeÀO trilayer (TL) is shown.
This study is supported by the National Natural Science Foundation of China (11574340), and the CAS-Shanghai Science Research Center (Grant No. CAS-SSRC-YJ-2015-01). Y. G. W. thanks the financial supported by Alexander von Humboldt Foundation.The authors also thank the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase) under Grant No. U1501501.The computational resources utilized in this research were provided by Shanghai

Fig. 2
Fig. 2 Catalytic cycle and corresponding charge cycle for the oxidation of CO catalysed by Au 2 @CeO 2 (111) via the traditional (a) and novel (b) L-H reaction mechanisms.The charges corresponding to gold adatoms are labeled as 1 and 2, respectively.Colour schemes: Au, gold; Ce, white; surface O, red; subsurface O, pink; surface O v , blue; subsurface O v , green; C, grey.