催化学报  2017, Vol. 38 Issue (9): 1597-1602   PDF    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
Hangjia Shen
Xianyuan Wu
Dahao Jiang
Li Xiaonian
Jun Ni
Identification of active sites for hydrogenation over Ru/SBA-15 using in situ Fourier-transform infrared spectroscopy
Hangjia Shen, Xianyuan Wu, Dahao Jiang, Li Xiaonian, Jun Ni     
Institute of Industrial Catalysis, Zhejiang University of Technology, Hangzhou 310014, Zhejiang, China
* Corresponding author. Jun Ni, Tel/Fax: +86-571-88320092; E-mail: Junni@zjut.edu.cn
Foundation item: This work was supported by the National Natural Science Foundation of China (21303163), the Natural Science Foundation of Zhejiang Province (LY13B030006, LY17B060006), the Qianjiang Talent Project in Zhejiang Province (QJD1302011), and the Scientific Research Fund of Zhejiang Provincial Education Department (Y201328681)
Abstract: The active sites for hydrogenation over Ru/SBA-15 catalysts were identified using in situ Fouri-er-transform infrared spectroscopy.The amount of active sites was proportional to the interfacial circumference of the Ru particles.In contrast, the rate of hydrogen spillover from Ru to the support was inversely proportional to the size of the Ru metal particles.Consequently, a catalyst with small Ru metal particles has a high rate of hydrogen spillover but a low density of active sites, whereas one with large Ru particles has a low rate of hydrogen spillover but a high density of active sites. The formation of these active sites is probably an intermediate step in hydrogen spillover.
© 2017, Dalian Institute of Chemical Physics, Chinese Academy of Sciences.
Published by Elsevier B.V. All rights reserved.
Key words: Ru/SBA-15 catalyst     In-situ Fourier-transform infrared spectroscopy     Hydrogenation     Active site     Hydrogen spillover    
利用原位红外光谱确证Ru/SBA-15的加氢活性位点
沈行加, 吴先元, 江大好, 李小年, 倪珺     
浙江工业大学工业催化研究所, 浙江杭州 310014
摘要:Ru/SBA-15催化剂具有高的氢气活化能力,因此被广泛应用在加氢和氢解反应中.一般认为Ru/SBA-15催化剂的高活性与金属Ru的高分散有关,然而有研究发现在氧化硅载体上还存在溢流的氢,这部分溢流氢也很可能参与加氢和氢解反应.这就产生了两个关键性的问题:(1)Ru/SBA-15的催化加氢活性中心是什么,是金属Ru还是载体SBA-15;(2)在金属Ru上解离的H是如何迁移到载体上的.因此,加氢活性位点及其形成机理的确认对理解Ru/SBA-15催化剂的高活性至关重要. 原位红外光谱可从分子层面研究在工作状态的催化剂表面活性位点的状态,进而推测可能的反应机理.通过与催化剂Pd/SBA-15,Ru/Al2O3和SBA-15比较发现,在氢气氛围中Ru/SBA-15催化剂的原位红外谱图中存在一个独特的位于1996cm-1的峰.由于在Pd/SBA-15,Ru/Al2O3和SBA-15上都不存在这个峰,因此该峰的形成是金属Ru和SBA-15相互作用的结果.此外,Si-O键在位于1866 cm-1的合频峰不随氢气氛围变化而变化,因此可排除这个峰属于Si-O键振动的倍频峰.为了排除该峰的产生是由于CO的吸附,我们采用脉冲引入CO的方法,发现在低的CO覆盖率下,红外谱图中位于2068 cm-1处出现了一个CO在Ruδ+上的线性吸附峰.随着CO覆盖率增加,该峰逐渐蓝移至2075 cm-1,同时位于2132 cm-1处的峰强度增强了,这两个峰都归属于Run+(CO)x物种的振动峰.这些CO的化学吸附强度都很高,即使在He气中吹扫1 h后仍然存在,而1996 cm-1峰的形成是可逆的.此外,低CO覆盖率下生成的吸附峰(2068 cm-1)的强度低于1996 cm-1峰的强度,因此可以排除1996cm-1峰属于CO吸附峰的可能.既然1996 cm-1峰的形成是可逆的,将这个峰归属于载体上氢的可能性也可排除,因为形成载体上氢的过程是不可逆的.另外,形成1996 cm-1峰的速率还证明了这个峰不属于金属Ru上吸附的氢,因为金属Ru上氢的吸附是很快的.通过以上分析,我们推断1996 cm-1峰应该指认为在Ru和SBA-15界面处位点的红外峰. 为了证明这一点,我们制备了不同Ru负载量的Ru/SBA-15催化剂,发现这个界面处位点峰的峰面积与金属Ru颗粒在载体上形成交界面的周长成正比,而峰达到稳态所需时间与Ru颗粒大小成反比.这说明H2在金属Ru上发生解离吸附后迁移到Ru和SBA-15界面处,形成了Ru-H-Si物种.当金属Ru的颗粒比较小时,与载体形成交界面的周长小,Ru-H-Si物种的数量少,体现在红外谱图上峰的峰面积小,但解离的氢迁移到该界面所需时间变短了.当金属Ru的颗粒比较大时,与载体形成交界面的周长大,Ru-H-Si物种的数量多,1996 cm-1峰的峰面积大,但解离氢的迁移慢了.此外,H-D交换实验还证明这个界面处的位点具有加氢活性.与文献报道的孤立Si-H物种的红外峰位置比较发现,Ru-H-Si物种具有明显的峰红移现象,说明该物种中的Si-H键活性很高,这可能是由于金属Ru将电子转移至Si-H键的结果.总之,以上结果清晰地表明这个1996 cm-1峰归属为结构是Ru-H-Si的活性位点.
关键词Ru/SBA-15催化剂    原位红外光谱    加氢反应    活性位点    氢溢流    

1 Introduction

The use of SBA-15-supported Ru catalysts in reactions such as hydrogenation [1, 2] and hydrogenolysis [3, 4] is becoming increasingly important because of the good performance of such catalysts in H2 activation. The surface area of SBA-15 is higher than those of previously reported supports for Ru catalysts [5, 6], and it has an ordered mesoporous structure, which favors high dispersion of Ru. SBA-15 has few acidic or basic sites and did not affect the studied reactions. The high reactivity of Ru/SBA-15 is therefore normally attributed to high dispersion of small Ru particles [7]. However, strong binding of spillover hydrogen to SiO2 has been detected on SiO2-supported Ru catalysts [8]. Moreover, the addition of Ru as a promoter in SiO2-supported catalysts improves the reducibilities and activities of the catalysts; this is also ascribed to hydrogen spillover [9]. These findings emphasize the importance of hydrogen spillover in reactions involving hydrogen, and raise the following issues in understanding of catalysts under working conditions. (1) Are the active sites for hydrogenation on Ru/SBA-15 catalysts Ru metal particles or the support? (2) How does hydrogen spill over from Ru to SBA-15? However, to the best of our knowledge, such a fundamental understanding of reaction mechanisms and the nature of active sites is still lacking.

In situ or in operando experimental methods need to be used to address these issues of real working catalysts [10]. In situ infrared (IR) spectroscopy is a powerful non-destructive technique, and is particularly used to probe surface sites and reaction mechanisms at the molecular level during catalytic reactions [11]. For example, IR spectroscopy has shown that benzene adsorption increases back-bonding of Pt d electrons to CO antibonding orbitals, based on the shift in the ν(CO) vibration to lower frequencies. This indicates that benzene is adsorbed on a single Pt atom (as in π arene complexes) with global transfer of electrons to Pt [12]. In situ IR spectroscopy has also shown that during the catalytic cycle in Baeyer-Villiger oxidation on Sn-zeolite beta, a ketone is coordinated to the Lewis acid center, resulting in carbonyl group activation [13].

Here, we report that the active sites for hydrogenation reactions over Ru/SAB-15 can be identified using in situ Fourier-transform IR (FTIR) spectroscopy. The formation of these active sites is probably the intermediate step in hydrogen spillover.

2 Experimental
2.1 Catalyst preparation

SBA-15 was synthesized using a previously reported method [14]. P123 (2 g) was completely dissolved in 1.6 mol/L HCl (75 mL), and tetramethyl orthosilicate (3.3 g) was added under vigorous stirring; stirring was continued for about 20 h. The temperature was kept at 40 ℃. The resultant white mixture was transferred to a Teflon-lined autoclave and aged at 100 ℃ for 24 h without stirring. The white powder was recovered, washed, and dried at 110 ℃ for 6 h. Calcination at 550 ℃ for 10 h in air gave mesoporous SBA-15.

Supported Ru catalysts with specific Ru loadings were prepared using a two-solvents method [15, 16]. Typically, an aqueous solution of RuCl3 (3.0 mL) was added gradually to cyclohexane (25.0 mL) containing a support material (SBA-15 or Al2O3, 1.0 g) at 25 ℃. After vigorous stirring for 15 min, the dark-brown solution became colorless and the originally white support became black. The black powder was separated by centrifugation and dried at 100 ℃ for 6 h; the product is denoted by Ru/support (wt%). Pd/SBA-15 (4%) was prepared by impregnation with PdCl2.

2.2 Catalyst characterization

In situ FTIR spectroscopy was performed using a Thermo Nicolet NEXUS 6700 spectrometer equipped with a mercury cadmium telluride (MCT) detector. A custom-made IR cell enabled spectra to be recorded at ambient temperature and at catalyst activation temperatures (up to 500 ℃). Self-supporting disks (9.0 mg) were prepared for IR studies and physisorbed molecules were removed using a He flow. FTIR spectra were collected with 32 scans at a resolution of 4 cm-1 at different time intervals until a steady state was reached. The catalyst was reduced in situ with H2 at 400 ℃ for 30 min at a ramping rate of 20 ℃/min, swept with He for 30 min, and cooled to 110 ℃ under He. These spectra were used as background for obtaining differential FTIR spectra.

The catalyst was subsequently treated using the following procedures; the gas flow rate was kept at 30 mL/min, unless otherwise stated. (a) H2 adsorption. The sample was treated in a H2 flow at 110 ℃ until a stable FTIR spectrum was obtained. Blank SBA-15 was treated in the same way for comparison. (b) H-D exchange. After procedure (a), the gas was switched to D2 to achieve H-D exchange. Spectra were collected until no change in the OH vibration peak was detected. The sample was then swept with He to expel D2 until a clean spectrum was obtained. The sample was then retreated with D2 and IR spectra of the sample without any interference from H2 were obtained. (c) CO adsorption. After the catalyst had been activated and cooled to 110 ℃ under He, CO was pulsed into the He stream three times using a six-way valve. After each pulse, the sample was swept with He for more than 10 min and spectra were recorded. After three pulses, the sample was swept with He for more than 8 h.

Transmission electron microscopy (TEM) was performed using a Tecnai G2 F30 S-Twin instrument (Philips-FEI Company) operated at 300 kV. After H2 adsorption, the samples were dispersed in ethanol and then supported on a Cu grid for TEM measurements.

3 Results and discussion

The in situ FTIR spectra in the range 2250-1850 cm-1 of the reduced samples after adsorption of H2 at 110 ℃ are shown in Fig. 1. No vibration bands from Pd/SBA-15 (4%), Ru/Al2O3 (4%), and SBA-15 were detected in this range, but a peak at 1996 cm-1 was observed for Ru/SBA-15 (4%). These results imply that formation of this band is associated with the presence of both Ru and SBA-15, but is not related to the Ru-H bonds in Ru metallic clusters, which would be detected for Ru/Al2O3. It has been reported that Ru-H species cannot be observed using IR, although adsorbed hydrogen species on various supported metals, e.g., Pt/Al2O3, Ir/Al2O3, Ni/Al2O3, Rh/Al2O3, and Ni/SiO2 have been observed [17, 18]. In addition, the intensity of the band arising from the combination tone vibration [19] of Si-O on Ru/SBA-15 at 1866 cm-1 did not change after exposure to H2, thus the intensity of the Si-O overtone band at 1985 cm-1 was constant. The possibility that the band at 1996 cm-1 can be ascribed to the Si-O overtone can therefore be ruled out.

Fig. 1. Differential FTIR spectra of catalysts exposed to H2 flow for 30 min at 110 ℃.

CO was adsorbed on Ru/SBA-15 (4%) to enable the band at 1996 cm-1 to be distinguished from those for adsorbed CO species. Fig. 2(a) shows that at low coverage of the Ru surface with CO, a band at 2068 cm-1 ascribed to CO linearly adsorbed on Ruδ+ appeared [20]. This band developed together with a band at 2132 cm-1, both related to Run+(CO)x multi-carbonyls [21], as coverage with CO increased until a band at 2177 cm-1 indicating physisorbed CO appeared. Few of these chemisorbed CO molecules were removed when the system was flushed with He at 110 ℃, even after 8 h (Fig. 2(b)). In addition, the intensity of the carbonyl band at 2068 cm-1 created by the first CO pulse (Fig. 2(a)) was lower than that of the band at 1996 cm-1 shown in Fig. 1. The wavenumber difference is therefore not caused by a surface coverage effect because carbonyl bands with wavenumbers higher than 2068 cm-1 are expected to appear with increasing CO coverage. The CO pulse adsorption results prove that the band at 1996 cm-1 does not arise from any carbonyl species formed on Ru.

Fig. 2. In situ FTIR spectra of CO adsorbed on reduced Ru/SBA-15 (4%) as a function of increasing CO coverage (a) and swept with a He flow at various times (b) at 110 ℃. CO was introduced using a six-way valve (about 2 μL) with He as the carrier gas.

Formation of the band at 1996 cm-1 was reversible in the case of Ru/SBA-15 (4%). The band developed gradually when the catalyst was exposed to H2 and reached its maximum intensity after 60 min (Fig. 3(a)). Subsequent flushing with He led to the disappearance of this band, as shown in Fig. 3(b). This band therefore cannot be assigned to spillover hydrogen on SBA-15 because the migration of hydrogen from Ru to SiO2 is irreversible [8]. The slow growth of this band during exposure for 1 h also indicates that it cannot be attributed to hydrogen chemisorbed on Ru because hydrogen saturates the Ru particle surfaces in a shorter time (less than a few minutes) [8].

Fig. 3. Evolution of in situ FTIR spectra of reduced Ru/SBA-15 (4%) in H2 at 110 ℃. (a) During exposure to H2 after reduction; (b) During flushing with He after exposure to H2.

Based on the above discussion, we conclude that the band at 1996 cm-1 cannot be ascribed to adsorbed H species on either Ru metal particles or the SBA-15 support. We postulate that the adsorbed H species are located at the interface between the Ru metal particles and the support. To confirm this hypothesis, we performed experiments using Ru with various interfacial circumferences on SBA-15 by tuning the Ru metal particle size. If the amount of H species is proportional to the interfacial circumference of Ru, these species reside along the Ru-SiO2 interface on the periphery of the exposed Ru particles. Fig. 4 shows TEM images of Ru/SBA-15 catalysts with various Ru contents. The mean Ru particle sizes calculated from 100 to 150 Ru particles in these TEM images are listed in Table 1. Based on a spherical model (equations (1) and (2)), the Ru particle circumferences were calculated using equation (3); they are presented in Table 1.

(1)
(2)
(3)
Fig. 4. TEM images and size distributions of Ru particles of (a) Ru/SBA-15 (1%), (b) Ru/SBA-15 (4%), and (c) Ru/SBA-15 (10%).
Table 1
Textural properties of catalysts.

where ρ is the density of Ru (12.37 g/cm3), N is the number of Ru particles, r is the Ru particle radius, ω is the Ru loading, and C is the total circumference of the Ru particles.

Fig. 5 shows the evolution of the in situ FTIR spectra of reduced Ru/SBA-15 with various Ru loadings at 110 ℃. The IR spectra of Ru/SBA-15 (4%) are shown in Fig. 3(a). The band areas at surface saturation with H2 were calculated to be 0.578, 3.01, and 3.83 for Ru/SBA-15 (1%), Ru/SBA-15 (4%), and Ru/SBA-15 (10%), respectively (Table 1). It can be seen from Table 1 that the band area at saturation increased linearly with increasing Ru particle circumference, indicating that the species corresponding to the bands in the range 2100-1800 cm-1 reside at the interfacial sites between Ru particles and the support. For the sake of convenience, we denote these sites by Sx. A precondition of the formation of special species is the presence of both Ru metal and SiO2 in the hydrogen atmosphere, therefore the process can be represented by the proposed model shown in Scheme 1. Molecular H2 dissociates on active Ru metal sites and saturates the Ru surface within a few minutes. The dissociated hydrogen atoms then migrate from the Ru metal sites to the Sx sites between Ru and SiO2; this process takes a longer time. The migration is also dependent on the Ru particle size; the smaller the Ru particle size is, the faster the migration of H atoms from Ru metal sites to Sx is, as shown in Figs. 3 and 5.

Fig. 5. Evolution of in situ FTIR spectra of reduced catalysts in H2 at 110 ℃. (a) Ru/SBA-15 (1%); (b) Ru/SBA-15 (10%).
Scheme 1. Hydrogen spillover to active sites on periphery of Ru on SBA-15 in H2.

The reactivities of these H species were evaluated using H-D exchange experiments. When D2 was introduced, the intensities of the hydroxyl bands at 3740 and 3660 cm-1 decreased, whereas the intensities of the deuteroxyl bands at 2760 and 2630 cm-1 increased, as shown in Fig. 6(a). Some of the hydroxyl groups corresponding to the band at 3660 cm-1 remained even after a few hours in D2. No changes in the bands at 2100-1800 cm-1 were observed. After flushing with a He flow for more than 1 h, the adsorbed deuterium had been completely swept away, as in the case of adsorbed hydrogen in He. We then introduced D2 again to investigate the nature of the bands in the range 2100-1800 cm-1. As shown in Fig. 6(b), a band centered at 1987 cm-1 was formed in D2 but it was red shifted relative to the band at 1994 cm-1 in Fig. 5(b), and of much lower intensity than that in H2; the former change can be ascribed to the deuterium isotope effect [22]. No replacement of H species by D species occurred during H-D exchange (Fig. 6(a)), therefore the much weaker band for D species cannot be attributed to the deuterium isotope effect alone. It mainly arises because D species are less stable than H species.

Fig. 6. In situ FTIR spectra of Ru/SBA-15 (10%) during H-D exchange (a) and on exposure to D2 at 110 ℃ after complete removal of any adsorbed H species (b).

Based on the above analysis, we propose a possible structure in which an H atom is coordinated to both Ru and Si, as shown in Scheme 2. The formation of this structure is probably the intermediate step in the hydrogen spillover process. The red shift in the Si-H stretching band [23] from 2100 to 1996 cm-1 is probably caused by the presence of Ru-H bonds from which electrons can be transferred to Si-H species. This electron transfer is responsible for the higher reactivity of Ru-H-Si species compared with those of isolated Si-H species reported in the literature [24-26]. It should be noted that hydrogenation may also take place on the Ru metal particle surfaces, but these active sites are more important when high selectivity toward targeted products is desired, especially when multi-bond hydrogenation is required [27].

Scheme 2. Proposed structure of active sites on Ru/SBA-15.
4 Conclusions

In situ FTIR spectroscopy showed the presence of special adsorbed H species on reduced Ru/SBA-15 catalysts but not on Pd/SBA-15, Ru/Al2O3, or SBA-15. These H species were differentiated from adsorbed CO species and were reproducible in H2. The formation rate of these H species depended on the Ru metal particle size, with small Ru particles giving faster formation of these H species. The amount of these H species was proportional to the Ru interfacial circumference, indicating that these species reside along the Ru-SiO2 interface on the periphery of exposed Ru particles. H-D exchange experiments suggested that these H species were active in hydrogenation. A possible structure for these active sites was proposed.

References
[1] C. Newman, X. B. Zhou, B. Goundie, I. T. Ghampson, R. A. Pollock, Z. Ross, M. C. Wheeler, R. W. Meulenberg, R. N. Austin, B. G. Frederick, Appl. Catal. A, 2014, 477: 64–74. DOI:10.1016/j.apcata.2014.02.030
[2] Y. F. Niu, H. Wang, X. L. Zhu, Z. Q. Song, X. N. Xie, X. Liu, J. Y. Han, Q. F. Ge, Microporous Mesoporous Mater., 2014, 198: 215–222. DOI:10.1016/j.micromeso.2014.07.030
[3] Y. M. Li, L. Ma, H. M. Liu, D. H. He, Chin. J. Catal., 2014, 35: 677–683. DOI:10.1016/S1872-2067(14)60032-2
[4] W. J. Li, L. M. Ye, J. Chen, X. P. Duan, H. Q. Lin, Y. Z. Yuan, Catal. Today, 2015, 251: 53–59. DOI:10.1016/j.cattod.2014.10.044
[5] H. Y. Liang, G. Z. Chen, S. Desinan, R. Rosei, F. Rosei, D. L. Ma, Int. J. Hydrogen Energy, 2012, 37: 17921–17927. DOI:10.1016/j.ijhydene.2012.09.026
[6] S. Akbayrak, S. Tanyıld, İ. Morkan, S. Özkar, Int. J. Hydrogen Energy, 2014, 39: 9628–9637. DOI:10.1016/j.ijhydene.2014.04.091
[7] Q. L. Yao, Z. H. Lu, K. K. Yang, X. S. Chen, M. H. Zhu, Sci. Rep., 2015, 5: 15186–15196. DOI:10.1038/srep15186
[8] D. O. Uner, M. Pruski, T. S. King, J. Catal., 1995, 156: 60–64. DOI:10.1006/jcat.1995.1231
[9] H. F. Xiong, Y. H. Zhang, K. Liew, J. L. Li, Fuel Process. Technol., 2009, 90: 237–246. DOI:10.1016/j.fuproc.2008.08.014
[10] B. R. Cuenya, Acc. Chem. Res., 2012, 46: 1682–1691.
[11] A. Vimont, F. Thibault-Starzyk, M. Daturi, Chem. Soc. Rev., 2010, 39: 4928–4950. DOI:10.1039/b919543m
[12] J. M. Basset, G. Dalmai-Imelik, M. Primet, R. Mutin, J. Catal., 1975, 37: 22–36. DOI:10.1016/0021-9517(75)90130-X
[13] A. Corma, L. T. Nemeth, M. Renz, S. Valencia, Nature, 2001, 412: 423–425. DOI:10.1038/35086546
[14] D. Y. Zhao, J. L. Feng, Q. S. Huo, N. Melosh, G. H. Frederickson, B. F. Chmelka, G. D. Stucky, Science, 1998, 279: 548–552. DOI:10.1126/science.279.5350.548
[15] M. Imperor, - Clerc, D. Bazin, M. D. Appay, P. Beaunier, A. Davidson, Chem. Mater., 2004, 16: 1813–1821. DOI:10.1021/cm035353m
[16] J. L. Liu, L. J. Zhu, Y. Pei, J. H. Zhuang, H. Li, H. X. Li, M. H. Qiao, K. N. Fan, Appl. Catal. A, 2009, 353: 282–287. DOI:10.1016/j.apcata.2008.10.056
[17] H. Ishikawa, J. N. Kondo, K. Domen, J. Phys. Chem. B, 1999, 103: 3229–3234. DOI:10.1021/jp9842852
[18] Z. Pei, T. H. Fang, S. D. Worley, J. Phys. Chem., 1995, 99: 3663–3670. DOI:10.1021/j100011a038
[19] M. Kantcheva, S. Sayan, Catal. Lett., 1999, 60: 27–38. DOI:10.1023/A:1019082218590
[20] G. H. Yokomizo, C. Louis, A. T. Bell, J. Catal., 1989, 120: 1–14. DOI:10.1016/0021-9517(89)90245-5
[21] K. Hadjiivanov, J. C. Lavalley, J. Lamotte, F. Mauge, J. Saint-Just, M. Che, J. Catal., 1998, 176: 415–425. DOI:10.1006/jcat.1998.2038
[22] C. K. Costello, J. H. Yang, H. Y. Law, Y. Wang, J. N. Lin, L. D. Marks, M. C. Kung, H. H. Kung, Appl. Catal. A, 2003, 243: 15–24. DOI:10.1016/S0926-860X(02)00533-1
[23] C. Copéret, D. P. Estes, K. Larmier, K. Searles, Chem. Rev., 2016, 116: 8463–8505. DOI:10.1021/acs.chemrev.6b00082
[24] C. Morterra, M. J. D. Low, J. Phys. Chem., 1969, 73: 321–326. DOI:10.1021/j100722a008
[25] C. Morterra, M. J. D. Low, J. Phys. Chem., 1969, 73: 327–333. DOI:10.1021/j100722a009
[26] A. Van Meerbeek, A. Jelli, J. J. Fripiat, J. Catal., 1977, 46: 320–325. DOI:10.1016/0021-9517(77)90215-9
[27] G. Kennedy, G. Melaet, H. L. Han, W. T. Ralston, G. A. Somorjai, ACS Catal., 2016, 6: 7140–7147. DOI:10.1021/acscatal.6b01640