Gold has been regarded as the most inert metal for a very long time, until the pioneer of Haruta [1, 2] found that the gold nanoparticles supported on some metal oxides showed surprisingly high catalytic activity for CO oxidation reaction at low temperature. Since then, the nano gold catalysis has attracted growing attention in the field of heterogeneous catalysis [3-5]. Technologically, numerous potential industrialization applications of nano gold catalysts were developed, such as synthesis of fine chemicals [6, 7], selective hydrogenations [8, 9], water gas shift reaction [10], carbon-carbon bond forming reaction [11], oxidation of organic compounds with molecular O2 [12, 13], pollution and emission control [14] and fuel cell applications [15] etc. Some of them have been commercially used, for instance, air purifier in various moto vehicles [16]. Scientifically, a huge number of theoretical and experimental studies have been devoted to understand the wide gap between the chemical inertness of bulk gold and high catalytic activity of nano gold [17-21].
The ability to catalyze CO oxidation reaction at low temperature is one of the unique properties of gold catalysts, and a dramatic increase in the number of studies has been performed to analyze the origin of the high catalytic activity and selectivity [6, 22-24]. Even so, it remains many puzzles about the simple CO oxidation reaction on gold catalysts. Experimentally, the account of surface active oxygen species is only 1%, thus it is very difficult to be characterized [25]. In addition, the activity and structures of gold nanoparticles strongly depend on the preparation methods, sizes of gold nanoparticles, and supports etc., which makes it very hard to identify the real active centers. Therefore, the importance of theoretical calculations is increasing for the elucidation of the nature of active centers in gold catalysts [4].
In the past two decades, the power of computational hardware, software gradually increased and therefore theoretical calculations became more and more popular and important to investigate structures, stabilities and activities of catalysts. Density functional theory (DFT) calculation relies only on the electron density, which significantly decreases the computational complexities. It becomes an increasingly powerful and useful tool to evaluate various systems, predict catalytic activities, and provide sufficient accuracy to be compared to experimental data [26]. So far, DFT calculation was widely used to investigate the mechanism of gold catalysis and provide fundamental physical insights into gold catalysis [17, 27, 28].
In this work, we will give a brief review of theoretical investigations on CO oxidation reaction catalyzed by gold nanoparticles, since most of the reviews focused on experimental results and the reviews of theoretical work were relatively rare [29, 30]. It should be noted that this work focused on gold nanoparticles and large gold clusters instead of small gold clusters consisting of only few gold atoms. Therefore, it was not concerned about some important effects such as some charge effects [31, 32] and even-odd effects [33] on gold etc., which are predominated on small gold clusters and vanish rapidly by increasing the cluster size [34]. In particular, special attention will be paid to the relativistic effects on gold catalysts in this work, since the relativistic effects strongly affect the color of gold, Au-Au bond strength, Au-Au bond lengths [35] and should play an important role for CO oxidation reaction.
In the CO oxidation reaction, the adsorption process of CO from the gas phase onto the gold surface is generally regarded as a key initial step of CO oxidation process. The adsorption of CO on all kinds of gold surfaces was well studied by DFT calculations. Mavrikakis et al. [36] found that the binding of CO on the flat Au(111) surface was very weak, but it became relatively strong at the step sites. Liu et al. [37] investigated CO adsorption on several gold surfaces and also found that the CO adsorption on Au(111) surface was very weak and the adsorption of CO increased considerably on surface defects. Fajín et al. [38] found that CO preferential adsorbed above the outermost surface atom at the edges of the (111) terraces in C-down fashion. The CO adsorption on bridge or hollow sites nearby the step was more favorable than that on top of gold atoms in the middle of the terraces, and the bridge sites along the step were also very stable for CO adsorption.
Barmparis et al. [39] presented a computational study of CO adsorption on gold nanoparticles and found that CO adsorption energy had a weak dependence on slab thickness but strongly depended on the exchange-correlation functional. Hussain et al. [40] investigated CO and NO adsorption on Au(111), Au(100), Au(110) and Au(310) surfaces, and found that the adsorption energy increased with the increasing degree of coordinative unsaturation of the gold atoms.
Some theoretical results showed that the calculated adsorption energies of CO on gold surfaces increased approximately linearly with decreasing coordination number (CN) of Au [41], regardless of the different functionals [42]. Furthermore, it was also found that the adsorption of NO, O2 and atomic oxygen increased with the decreasing Au CN [41, 43]. Based on this, some researcher proposed that the low-coordinated sites of gold atoms were active sites and thus gold nanoparticles with large amounts of low-coordinated sites showed high catalytic activity [18, 37].
We also observed the near-linear relationship between the CO adsorption energy and Au CN based on the DFT calculation results on Au(321), Au(711), Au(100) and Au(111), as shown in Fig. 1(a) [44]. However, we found that the linear relationship did not work anymore for CO adsorption on Au(100)-hex(1×5) surface. In the Au(100)-hex(1×5) structure, the top layer of Au(100) structure is reconstructed to an Au(111) structure, and the relative positions of surface Au(n) (n=1-6) atoms after relaxation are shown in Fig. 1(b). It is obvious that the CO adsorption energies on Au(3) and Au(4) of Au(100)-hex(1×5) surface deviate dramatically from the linear relationship, as shown in Fig. 1(a).
To trace the essential reason that the near-linear relationship between CO adsorption energy and the Au CN works well on some gold surfaces but not on the Au(100)-hex(1×5) surface, we investigated thoroughly the CO adsorption on gold surfaces by Hückel theory and DFT calculations [44]. It was found that only the neighboring gold atoms locating at the down area of Au-CO could strengthen CO adsorption, whereas the neighboring gold atoms locating at the side area of Au-CO weakened CO adsorption. It indicates that the CO adsorption is not only dependent on the Au CN, but also relies on the location of the neighboring gold atoms. The extended researches indicated that the strong relativistic effects of gold strengthened the influence from the neighboring gold atoms, because of the wide spatial extension of d orbitals of gold [28].
To gain detailed information about the O2 adsorption is a relevant issue in order to arrive at a deeper understanding of the CO oxidation reaction. The DFT results from N rskov et al. [36] showed that the O2 adsorption was extremely weak on Au(111) surface and gold clusters. Mills et al. [45] found that O2 did not bind to a planar gold cluster and only bound to gold clusters deposited on an Au(111) surface, which indicated that the surface roughness was an essential factor for O2 adsorption. Molina et al. [30] found that O2 adsorption on Au34 cluster and Au/MgO was very weak. The adsorption of molecular oxygen on the Au(321) surface was calculated by Fajín et al. [46], and it was found that the interaction between O2 and gold was very weak, and the structure with the O-O axis planar to the surface terraces was more favorable.
The calculation results of Gong et al. [47] showed that O2could hardly adsorb on flat and stepped gold thin films, but wire-structured gold could adsorb both CO and O2 rather strongly. They proposed a generalized structural model based on the wired-structured film for active gold. Günay et al. [48] investigated O2 adsorption on a gold cluster, and found that the stable O2 adsorption required the presence of an unpaired electron on gold cluster. Lee et al. [49] calculated O2 and O adsorption on Au38, and found that the most stable sites for O2 and O adsorption were bridge and hexagonal sites with the adsorption energies of -0.99 and -3.93 eV, respectively. In total, most theoretical work showed that the adsorption of O2 molecule on gold surfaces was relatively weak except for the wire-structured gold or some special gold clusters.
Liu et al. [17] calculated O2 adsorption on Au/TiO2(110) and found that O2 adsorption at the interface was highly favored compared to the adsorption on the gold surfaces. Molina et al. [50] also found that O2 adsorption at the interface of Au/TiO2 was reasonably strong. Remediakis et al. [51] reported that the adsorption site for O2 on Au/TiO2 was at the interface between gold cluster and Ti5c (surface 5-fold coordinated Ti atoms) on the surface, which was confirmed by Koga et al. [52].
Chen et al. [29] demonstrated that the interface of Au/TiO2 was far more favorable for O2 adsorption than low-coordinated gold atoms. The DFT results from Cai et al. [53] proved that O2 could only be adsorbed at the interfacial sites on Au/TiO2(110) surface. These DFT calculation results clearly showed that O2 preferred to adsorb at the interface between gold and TiO2 on Au/TiO2 catalysts. Due to the strong adsorption of O2 at the perimeter edge of gold oxide contacts, it was proposed that the perimeter sites were the active sites for CO oxidation in some studies [17, 50, 54].
Gong et al. [55] investigated O2 adsorption on Au/CeO2, and found that only the surface Ce of CeO2(111) could adsorb O2 with high stability, which showed that the O2 preferred to be adsorbed on CeO2surface than on gold surface. Han et al. [56] investigated CO and O2 adsorption on gold clusters with a 2D-3D structure on CeO2-x(111), and found that the Au CN and charges of gold atoms dominated the adsorption strength of CO and O2 and the influence of the interface of Au/CeO2was relatively weak.
In the last two decades, the O adsorption on gold surface has been investigated by DFT calculations on a variety of gold surfaces including gold single crystal and gold clusters. N rskov et al. [57] calculated that the dissociative chemisorption energies for O2 on gold surface was 0.54 eV, which implied that gold did not bind oxygen at all [43]. Shi et al. [58] investigated the adsorption of oxygen on Au(111) surface by DFT calculations, and found that atomic oxygen adsorption was weak compared to molecular oxygen. However, the DFT results from Fajín et al. [46] showed that the dissociation of O2 on Au(321) surface was exothermic. In addition, O2 temperature-programmed desorption (TPD) peaks on gold surfaces are higher than 500 K [59], indicating that the atomic oxygen adsorption on gold surface is not so weak.
Liu et al. [37] calculated atomic oxygen adsorption on a series of gold surfaces, and found that the bridge site of monatomic step edge was the most stable site for O adsorption on Au(221) and Au(211) surfaces. The 3-fold hollow site was the most stable site on Au(111) surface, which was confirmed by Zeng et al. [60]. In addition, Liu et al. [37] also found that the adsorption energy on stepped surface increased compared to the flat Au(111) surface. Fajín et al. [46, 61] studied atomic oxygen adsorption in stepped Au(321) surface and found that the atomic oxygen preferred to occupy the fcc hollow cavities near the edge of the step connecting (111) terraces.
We investigated atomic oxygen adsorption on gold surfaces by Hückel theory and DFT calculations [62] and found that the atomic oxygen preferred to locate at multifold sites on Au(111), Au(100), Au(110), Au(211) and Au(321) surfaces, which agreed well with the results reported by Fajín et al. [46]. Furthermore, we found that the linear O-Au-O structure was the most stable structure for atomic oxygen adsorption, and the corresponding structure on Au(110) surface is shown in Fig. 2.
N rskov et al. [41, 43] found that the binding energies for oxygen on gold were approximately linearly relationship with Au CN, which means that lower Au CN will induce to the stronger oxygen adsorption. However, Liu et al. [37] found that the oxygen adsorption on Au kink-I and kink-II with lower Au CN did not show stronger adsorption compared to the gold atoms at the steps, which is quite different from N rskov’s observations [41, 43].
We investigated Au-O interactions in various surroundings by DFT calculations, and found that the adsorption of oxygen on gold (Au-O bond) could be strongly weakened by the neighboring gold atoms locating in the plane perpendicular to the Au-O bond, while the neighboring gold at the downside would strengthen Au-O bond [28]. It might be used to explain the inconsistency between the observations of N rskov et al. [41, 43] and Liu et al. [37]. It indicates that the O adsorption as well as CO adsorption depends not only on the Au CN but also on the location of the neighboring gold atoms.
Since the formation of gold oxide films was believed to be related to the unique oxidation catalysis of gold [63, 64], the nature of gold oxide films has been theoretically explored in the last decades. Shi et al. [58] investigated the adsorption of oxygen on Au(111) surface by DFT calculations, and found that a thin surface-oxide-like configuration was the most favorable structure compared to atomic oxygen adsorption. In this structure, the oxygen atoms are quasi 3-fold coordinated to gold atoms and the gold atoms in the surface layer are 2-fold coordinated to oxygen atoms. The surface-oxide-like structure is illustrated in Fig. 3(a). It was expected to be stable up to about 420 K at atmospheric pressure. They also calculated the adsorption of oxygen on Au(100) and Au(110) surface, and found that a low coverage (0.1ML) on added-row reconstructed surface and 1 ML oxygen covered (2×1) missing row structure were the most stable structures on Au(100) and Au(100) surfaces, respectively [65].
Fajín et al. [61] studied the consecutive deposition of oxygen atoms on the Au(321) surface by DFT calculations, and observed that the increase of the oxygen coverage resulted in the formation of highly oxidized structures. The simulation results from Baker et al. [66] showed that oxygen chemisorbed on the surface was the primary species at lower coverage and lower temperature. At higher coverage and temperature, the gold oxide-like species became dominant. The energy barrier for the formation of the surface oxide was less than 0.1 eV. Boronate et al. [67, 68] investigated O2 dissociative adsorption on Au38, and found that O-O bond break on Au(100) facet was favored. The activation barrier was only 7.6 kcal/mol and the process was highly exothermic. They proposed that the oxygen adsorption could reach coverage of one monolayer oxide, where all oxygen atoms were forming the most stable O-Au-O linear structure.
We calculated the structures and stabilities of gold oxide films on gold surfaces in O2 atmosphere by a combination of thermodynamics and DFT calculations [69], and found that the dual chain structure (Fig. 3(b)) was more stable on most gold surfaces at lower oxygen chemical potential except for on an Au(111) surfaces, on which the graphene-like 12MR structure (Fig. 3(c)) was more stable. While at the higher oxygen chemical potential, the connected dual chain structure (Fig. 3(d)) became more stable. The phase transformation point from the dual chain structure or graphene-like 12MR structure to connected dual chain structure depended on the lattice mismatch between the connected dual chain structure and the gold surface.
In general, there are at least three reaction mechanisms for CO oxidation reaction on gold surfaces, namely Eley-Rideal (ER) mechanism, Langmuir-Hinshelwood (LH) mechanism, and O2 dissociation mechanism, as shown in Fig. 4. For the ER mechanism [27], the adsorbed CO reacts directly with gaseous O2. While for the LH mechanism [17], the adsorbed CO reacts with coadsorbed O2. Normally, an OCOO (also known as CO-O2) intermediate complex is formed and then it decomposes to CO2 and adsorbed O [70]. The remaining adsorbed O reacts with another CO molecule to generate a CO2 molecule. For the O2 dissociation mechanism, O2 dissociates to two adsorbed oxygen atoms before reacting with adsorbed CO to form CO2.
It should be noted that Liu et al. [71] proposed a trimolecular LH mechanism (3LH) for CO oxidation on gold surface. They found that the coadsorbed CO molecule at a unique triangular Au3site could act as a promoter for the bond break of an O-O bond, and spontaneously forming two CO2 as product. The coadsorbed CO on the triangular Au3 site significantly promotes the O-O bond scission in OCOO intermediate, therefore, it is a CO self-promoting oxidation reaction.
Liu et al. [17, 37] performed a DFT calculation on CO oxidation reaction on gold surfaces and found that the calculated reaction barrier involving O2 dissociation was much higher than that of CO reaction with molecular O2. Therefore, they believed that the possibility of O2 direct dissociation mechanism on gold surfaces was ruled out at low temperature, and they suggested that the reaction mechanism followed the reaction between O2 and CO through CO-O2 intermediate complex. Molina et al. [27, 72] found that CO oxidation for all gold systems proceeded via CO adsorption, trapping of O2, formation of CO-O2 intermediate, generation of CO2 and reaction of remaining atomic oxygen with gaseous CO to form CO2. They found that the direct separation or dissociation of O2 was not favorable.
Fajín et al. [38] investigated CO oxidation on Au(321) surface, and found that the predissociation of molecular oxygen on the Au(321) surface for CO oxidation was energetically less favorable than the CO direct reaction with molecular oxygen. Chang et al. [73] calculated CO oxidation on unsupported Au55 clusters and found that the reaction path consisted of an intermediate involving CO-O2 complex. Chen et al. [74] investigated the CO oxidation on a gold nanoparticle by DFT calculations, and found that Au29 nanoparticle exhibited high catalytic activity for CO oxidation. The catalytic process preferred via the LH mechanism at lower temperature, whereas the ER mechanism became more competitive at higher temperature. These DFT results clearly showed that CO direct reaction with molecular oxygen was more favorable than the predissociation of molecular oxygen before reacting with adsorbed CO on gold surfaces.
Liu et al. [17] studied CO oxidation on TiO2 supported Au by DFT calculations, and found that CO oxidation reaction occurred at the interface between gold and the oxide and the reaction barrier was very low. The key step of CO oxidation on Au/TiO2 was the adsorption of molecular O2 at the interface. The TiO2 support enhanced electron transfer from the gold to the antibonding states of O2, which significantly promoted O2 adsorption and activates O2 towards CO oxidation.
Molina et al. [50] investigated O2 adsorption and CO oxidation on gold nanoparticles supported on TiO2(110) surface, and found that O2 adsorption occurred either at the trough Ti (Ti5c) or leaning against the gold particles. The presence of gold strongly stabilized the adsorbed O2 through the electronic charge transfer from the Au to O2. The adsorbed O2 could react with adsorbed CO at the edge sites of gold particles and form CO2 easily. Molina et al. [50] found that CO-O2 intermediate complex was no longer formed as a metastable reaction intermediate on Au/TiO2, which is different from the conclusion from other researchers that the adsorbed CO reaction with an oxygen molecule normally gave formation of a CO-O2 intermediate complex [18, 37, 75, 76].
Li et al. [77] systematically studied the size and shape dependence of Aun/TiO2 on CO oxidation reaction, and found that the perimeter sites of Au-cluster/TiO2 could significantly promote the CO oxidation. The reaction occurred between the adsorbed CO at a perimeter Au site and a dangling O2 on the neighboring Ti5c site. They simulated the generation of oxygen vacancies at the interface between gold cluster and TiO2(111), and found that both of LH mechanism and Mars-van Krevelen (M-vK) mechanism were possible towards CO oxidation reaction on Au/TiO2 and the selection depended on the structural fluxionality of the gold cage clusters on TiO2.
Koga et al. [52, 54] examined CO oxidation over Au/TiO2(110) by DFT calculations and found that O2 was strongly adsorbed on the Ti5csite next to the perimeter. CO was activated by the back donation from gold, which induced to the lower reaction barrier for CO with O2 on the perimeter Ti5c than that on a more remote Ti site. Therefore, the Ti5csites next to the perimeter of Au/TiO2 played an important role for CO oxidation.
Vilhelmsen et al. [78] investigated the catalytic sites at the interface perimeter of Au/TiO2 and found that the active sites were strongly dependent on the surface direction. The gold on the [110] direction of TiO2(110) surface showed much higher catalytic activity than that on the [001] direction, and they assigned it to the too week CO binding energy along the [001] direction. Duan et al. [79] presented DFT calculations on CO oxidation mechanism on a gold rod on TiO2(110) surface along the [110] direction, and found that both LH mechanism with an adsorbed Au/O/Ti5c oxygen species and an Au-assisted M-vK mechanism with Obri were acceptable reaction mechanisms for CO oxidation reaction.
Camellone et al. [80] investigated CO oxidation mechanism on Au/CeO2 catalysts by DFT calculations with a Hubbard term (DFT+U), and found that both of Au+and Au3+were active centers for CO oxidation, but Au+would turn to Auδ-and then be deactivated with the presence of oxygen vacancies. The catalytic cycle for Au3+ on CeO2 was via formation of surface oxygen vacancies by CO reacting with surface oxygen, adsorption of molecular oxygen on oxygen vacancies, formation of oxygen adspecies and reaction with CO.
Henkelman et al. [81] investigated the CO oxidation mechanism on gold nanoparticles supported by CeO2 and proposed three reaction pathways for CO oxidation: by coadsorbed O2 at gold nanoparticles, by lattice oxygen of CeO2 (M-vK mechanism), and by O2 adsorbed at the Au-Ce3+ interface. They found that higher energy was required for the M-vK mechanism in the case of Au13 supported on CeO2, therefore, it was unfavorable and only active at high temperature.
Ghosh et al. [82] investigated small gold clusters on ceria surface during CO oxidation by DFT calculations and found that the morphology of gold clusters was influenced significantly by molecular adsorption. Their results indicated that M-vK mechanism displayed low activation energy on Au/CeO2, and they assigned it to the gold cluster fluxionality and the lability of CO2 intermediate at the interface.
Teng et al. [83] studied the structures of Au clusters on a CeO2(111) surface and found that the stability of surface oxygen vacancies depended on the size of the gold cluster and relative position of gold cluster and oxygen vacancy. Liu et al. [84] investigated the CO adsorption on Au/CeO2and found that the CeO2 support acted as an electron buffer which benefited the CO adsorption on Au/CeO2.
Song et al. [85] performed a computational DFT study of CO oxidation on a gold nanorod supported on CeO2(110), and investigated three CO oxidation mechanism: M-vK mechanism, co-adsorption mechanism with adsorbed CO and O2, and stepwise mechanism with O2 dissociation followed by CO oxidation, and found that all three mechanisms exhibited nearly similar overall reaction barriers, and M-vK mechanism and the stepwise mechanism strongly depended on the support termination of CeO2 surface.
Liu et al. [37] calculated O2 dissociation on Au(111) surface and found that the reaction barrier was very high. They found that the step sites were most active sites for O2 dissociation after examining many gold structures. However, even on the most active step sites, the estimated initial sticking coefficient was too small that O2 could not dissociate on them. Su et al. [86] confirmed that O2 adsorption on Au(111) surface was very weak and direct dissociation barrier of O2was relatively high. Fajín et al. [46] calculated the dissociation of molecular oxygen on the Au(321) surface, and found that it was exothermic and the reaction barrier was 1.00 eV. Molina et al. [72] found that O2 dissociation was not favorable on gold nanoparticles supported by MgO(100) surface.
Although many theoretical results supported that the CO oxidation via CO reaction with O2 was facile compared to O2 dissociation before reacting with CO no matter on gold nanoparticles or on Au/oxide systems [17, 18, 27, 37, 75, 76], some experimental results proved that Au/metal oxides catalysts could activate molecular oxygen directly [21, 25, 87-89]. For example, Weiher et al. [88] found that gold-oxygen complexes could be formed on active Au/TiO2 catalysts without the presence of CO by X-ray absorption near-edge structure (XANES). Widmann et al. [25] found that the active oxygen species was quite stable and should be adsorbed atomic oxygen, which was generated by dissociation of molecular oxygen. Takeda et al. [21] observed the morphology changes of gold nanoparticles on CeO2 and formation of surficial gold oxide in O2 atmosphere by transmission electron microscopy (TEM).
To understand the mechanism of O2 direct activation on gold catalysts, a great deal of theoretical fundamental studies has been performed. Roldan et al. [90] found a common pathway for O2 dissociation on gold nanoparticles by DFT calculations and concluded that it was very important for the presence of low-coordinated gold atoms and the critical size of gold nanoparticles. Boronat et al. [67, 68] investigated the adsorption and dissociation of O2 on gold surfaces, gold nanoparticles and Au/TiO2 surface, and found that the arrangement of the gold surface atoms was vital for O2 dissociation and the Au(100) facet was favored for scission of O-O bond.
Hussain et al. [91] calculated CO oxidation on extended gold surfaces and found that a site consisting of four gold atoms in a square geometry had sufficient ability to dissociate O2, and a Au38 cluster exposing this site provided the most favorable energetics for CO oxidation, which were in good agreement with Boronat et al.’s results [67, 68]. Duan et al. [79] identified Au/Ti5c site at the interface boundary to activate O2 and found that O2 dissociation with a reaction barrier of 0.5 eV was facile at room temperature. Koga et al. [92] found that O2 dissociation between two Ti5c sites was probable if the relative O2 population was small enough. Furthermore, they proposed that both of undissociated and dissociated O2 on Ti5c sites were possible active oxygen species, depending on the ratio of O2 population to the contact area.
To find out the most active sites for O2 dissociation on gold catalysts, we analyzed O2activation in an Au-O-O-Au structure by Hückel theory, and found that an increase of Au-O orbital interaction could significantly promote O2 adsorption and dissociation [93]. Since the O-Au-O structure could increase Au-O interaction remarkably [62, 28], we built a realistic model with a unique double linear O-Au-O structure at the perimeter of the Au/metal oxide interfaces to activate adsorbed O2, as shown in Fig. 5.
It was found that the O2 dissociation barriers with this structure were only 0.12 and 0.17 eV for Au/TiO2 and Au/CeO2 systems, respectively [93]. Therefore, the direct O2 dissociation mechanism is also a possible reaction pathway in Au/oxide systems. In addition, the unique double linear O-Au-O structure at the perimeter of the Au/metal oxide interfaces implies that the active centers with at least two-layer gold atoms are necessary towards CO oxidation, which agrees well with the experimentally observations by Goodman et al. [22, 94] that two-layer gold particles supported on TiO2 showed the highest catalytic activity for low-temperature CO oxidation reaction.
CO oxidation on gold has been investigated from a theoretical point of view, and some relevant conclusions about the influences of coordination number, charge transfer and relativity of gold have already been obtained. Since the experimental results showed that the gold nanoparticles exhibited higher reactivity towards CO oxidation than larger ones [22, 23], and the fraction of low-coordinated gold atoms decreased with the increasing the size of gold nanoparticles [43, 95], the low-coordinated gold atoms at the step, corner or edge were well investigated by DFT calculations.
N rskov et al. [36] found that the adsorption energy of O and CO was not dependent on the number of gold layers if there were more than two layers. They found that the steps and the strained surface were crucial for CO oxidation on gold surface and they suggested that the gold atoms on the corners and edges of gold nanoparticles were the active sites [43]. Liu et al. [37] also found that gold steps were the most active sites for CO and atomic oxygen adsorption, therefore, they believed that these steps were reaction active sites.
Xu et al. [96] investigated the O2 dissociation on strained and stepped gold surface by DFT calculations, and found that the steps and tensile enhanced the adsorption of atomic oxygen, and substantially facilitated O2 activation on a gold surface. Brodersen et al. [97] concluded that the reactivity of nano gold particles was dominated by low-coordinated corner-like atoms. Pascucci et al. [98] compared the reactivity of Cu, Ag and Au for CO oxidation by DFT calculation, and found that the reactivity of gold nanoparticles should be attributed to large number of low-coordinated sites.
The fact that low-coordinated gold atoms or tensile gold surface structures benefit the CO adsorption and O2 activation can be easily understood by the relativistic effects [28]. For the Au-O (or Au-CO) bond at the step sites or in a tensile surface, the inhibiting effects from the neighboring gold locating in the plane perpendicular to the Au-O (or Au-C) bond is weakened due to smaller Au CN or longer Au-Au bond length, therefore, the Au-O (or Au-CO) bond will be stronger. Furthermore, the reaction barrier for O2 activation decreases according to the Hückel theory [93]. As a result, the low-coordinated gold or gold in a tensile surface displays higher catalytic activity for CO oxidation.
Li et al. [77] studied CO oxidation on subnanometer gold on TiO2(110) by a combination of DFT calculations and microkinetics analysis, and found that a hollow-cage Au18/TiO2 displayed the highest catalytic activity towards CO oxidation among of Au1-4/TiO2, Au7/TiO2 and Au16-20/TiO2. In addition, the pyramidal Au18 isomer exhibited lower activity comparable to hollow-cage Au18 isomer, but it displayed much higher activity than pyramidal Au19 and Au20 isomers. They assigned the higher catalytic activities to the weaker CO adsorption on TiO2 than on gold clusters, which leaded to an increase of the ratio of O2/CO adsorption energy, an increase of the probability for O2 to occupy the Ti sites, and a decrease of the requirement for meet the critical line for CO oxidation.
The effects of charge transfer in CO oxidation reaction were deeply investigated by DFT calculations. Sanchez et al. [99] found that electron transfer from the surface to the gold cluster played an essential role in the activation of nano-size gold cluster for CO oxidation. Liu et al. [17] found that the oxide promoted the charge transfer from gold to anti-bonding orbitals of O2, which stabilized the adsorbed O2and activates O2 for CO oxidation. Koga et al. [52, 54, 100] found that the charge transfer from Au to O2 did not require direct bonding between Au and O2, and Ti5csites could activate O2 towards CO oxidation even~1 nm away from gold cluster. They also found that the charge transfer occurred even on the inversed TiO2(110)/Au(112) catalyst.
Molina et al. [72] investigated CO oxidation on gold nanoparticles supported by MgO(100) and found that the local reactivity of gold atoms was determined by Au-Au coordination and the MgO supports provided excess electrons to gold to form ionic bonds with CO-O2 intermediate. The DFT results of Guo et al. [101] showed that negative charge not only benefited the O and O2 adsorption but also lowered the dissociation barrier for O2. However, Henkelman et al. [81] found that positively or negatively charged gold ions had no direct effect on the pathway and rate of CO oxidation. Based the analysis by Hückel theory, we found that the excess electrons from gold or some electron donors benefited the O adsorption and O2 activation, but the excess electrons on surface oxygen ions in metal oxide blocked the O adsorption and O2 activation and thus hindered the CO oxidation reaction [93].
Since gold displays the strongest relativistic effects among the first 100 elements, the relativistic effects play an important role on color of gold, Au-Au bond strength, Au-Au bond lengths etc. [35]. The relativistic effects on catalytic activity of gold were discussed in some recent theoretical work. Tang et al. [102] calculated CO oxidation reaction on gold nanoclusters by all-electron relativistic DFT calculations, and found that high spin energy profiles and spin crossing played an important role in CO oxidation reaction.
Liu et al. [103] calculated the adsorption of many small molecules on helical gold nanorods by a relativistic DFT method, and found that the Au-C, Au-N, Au-S and Au-O bonds were surprisingly strong, but the Au-H bond was relatively weak. Kanoun et al. [104] calculated the activation of molecular oxygen on noble metal nanoparticles, and found that the reaction barrier for O2 dissociation on Au38 was reduced from 30.1 to 11.4 kcal/mol after involving of the relativistic effects.
We have investigated the Au-O interactions in different surroundings by DFT calculations in the scalar relativistic level [28] and found that Au-O bond could be significantly enhanced by the linear O-Au-O structure, because π-bonding play a key role in Au-O interaction. The strengthened Au-O bond in linear O-Au-O structure could be used to activate molecular O2 in CO oxidation, as shown in Fig. 5. The Au-O interactions can be strongly suppressed by the presence of neighboring gold atoms in the plane perpendicular to the Au-O bond [28]. All the interesting phenomena that some Au-O bonds especially for linear O-Au-O bonds are strengthened and some Au-O bonds are weakened can be assigned to the spatial extension of d orbitals of gold wider than other transition metals, which is induced by the strong relativistic effects of gold, as shown in Fig. 6.
The widely spatial extension of d orbitals of gold could strongly strengthen significantly the Au-X (X=C, S, O, etc.) π-bonding interactions, but could not remarkable enhance the σ-bonding interactions (for example Au-H interaction) [28]. In addition, the neighboring gold atoms, especially for the neighboring gold atoms in the plane perpendicular to the Au-X bond, would dramatically weaken the Au-X bond. As a result, gold atoms in the bulk region with higher coordination number are strongly suppressed by neighboring gold atoms and thus are chemically inert towards CO oxidation reaction or some other reactions. However, low-coordinated gold atoms abounding at the steps, corners or edges in gold nanoparticles are capable to adsorb CO and atomic oxygen, and benefit the CO oxidation reaction. The oxygen atoms locating at the perimeter of Au/oxide enlarge the differences between gold nanoparticles and bulk gold by the linear unique O-Au-O structures.
Some theoretical works about CO, O2 and O adsorption and CO oxidation on gold catalysts and the influences of coordination number, charge transfer and relativity of gold on CO oxidation were briefly reviewed in this work. It was found that CO preferentially adsorbed above the gold atoms with smaller coordination number. There is a near-linear relationship between CO adsorption energies and CN of gold in most cases except for the Au(100)-hex(1×5) surface, since CO adsorption depends on Au CN as well as the relative location of the neighboring gold atoms.
O2 adsorption on gold surfaces is normally very weak unless adsorption on wire-structured gold or some special gold clusters. The O2 adsorption is relatively strong at the interface of Au/TiO2 or on the surface of CeO2 in Au/CeO2 system. The most stable structure for the atomic oxygen adsorption relies on gold surfaces, and a linear O-Au-O structure is preferable in most cases at low coverage of oxygen. At high coverage, a surface gold oxide film can be formed on the gold surfaces and the structure of gold oxide film depends on the gold surfaces and oxygen chemical potential.
Some theoretical results supported that CO oxidation mechanism followed CO reaction with O2 via CO-O2 intermediate on gold surface, Au/TiO2 and Au/MgO system. The reaction mechanism for CO oxidation on Au/CeO2 is still under debate. Although the mechanism of O2 predissociation before reacting with CO on gold catalysts was not supported by many theoretical works, some theoretical researchers still found some ways for O2 direct dissociation with low reaction barrier, for example, a unique double linear O-Au-O structure at the perimeter of the Au/metal oxide interfaces to activate adsorbed O2.
The gold atoms with lower coordination number were regarded as active sites in some theoretical work, which could be clearly explained by the strong relativistic effects of gold. The relativistic effects could also be used to understand the unique behaviors of CO adsorption, O adsorption, O2 activation on gold surfaces and the wide gap between the chemical inertness of bulk gold and high catalytic activity of gold nanoparticles especially supported by some oxides. Due to the complexity of gold catalysts, a combination of theoretical calculations and experiments is strongly desired to understand the mechanism of low temperature CO oxidation reaction catalyzed by gold nanoparticles.