催化学报  2020, Vol. 41 Issue (9): 1337-1347      DOI: 10.1016/S1872-2067(20)63580-X   PDF    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
Shunlin Li
Lili Wang
Mengmeng Wu
Yafei Sun
Xiaojuan Zhu
Ying Wan
Measurable surface d charge of Pd as a descriptor for the selective hydrogenation activity of quinoline
Shunlin Li, Lili Wang, Mengmeng Wu, Yafei Sun, Xiaojuan Zhu, Ying Wan     
The Education Ministry Key Laboratory of Resource Chemistry, Shanghai Key Laboratory of Rare Earth Functional Materials, and Department of Chemistry, Shanghai Normal University, Shanghai 200234, China
* Corresponding author. Ying Wan, Tel/Fax: +86-21-64322516; E-mail: ywan@shnu.edu.cn
This work was supported by the National Natural Science Foundation of China (21773156), the Ministry of Education of China (PCSIRT_IRT_16R49), the International Joint Laboratory on Resource Chemistry of China (IJLRC), the Shanghai Sci. & Tech. and Edu. Committee (17JC1404200, 19070502700), and the Shanghai Gaofeng & Gaoyuan Project for University Academic Program Development
Abstract: AuPd nanoalloys with tunable Pd concentrations have been synthesized and used as model catalysts. They have been directly imaged by high-angle annular dark-field scanning transmission electron microscopy and investigated by thorough analyses of their extended X-ray absorption fine structure, X-ray absorption near-edge structure, X-ray diffraction and X-ray photoelectron spectroscopy measurements. The bimetallic nanoparticles are embedded in a carbonaceous matrix and have almost an identical structure at the atomic level and the same electronic properties as AuPd bulk alloys with the same compositions. The d-electron increase at surface Pd sites is determined by the Pd concentration of the alloy. Similarly, their activation entropy and catalytic activity for the hydrogenation of quinoline is related to the Pd concentration, with Au50Pd50 the most active of the alloys investigated. An almost 11 times higher activity was achieved compared to a pure Pd catalyst. The experimentally measurable surface d charge at the Pd sites in the AuPd was found to linearly correlate with the activation entropy and catalytic activity for the hydrogenation of quinoline. The alloy structure is stable, showing negligible metal segregation, dissolution-redeposition and aggregation during the hydrogenation process which involves strong adsorption.
© 2020, Dalian Institute of Chemical Physics, Chinese Academy of Sciences.
Published by Elsevier B.V. All rights reserved.
Key words: Surface d charge    AuPd alloy    Activation entropy    Quinoline    Hydrogenation    
利用实验测量的Pd表面d带电荷密度描述喹啉催化加氢活性
李顺林, 王丽丽, 吴萌萌, 孙雅飞, 朱小娟, 万颖     
上海师范大学化学系, 资源化学教育部重点实验室, 上海市稀土功能材料重点实验室, 上海 200234
摘要:杂环化合物如喹啉等在Pd表面发生较强的化学吸附,占据活性位,导致钯炭催化剂活性低.因此,需要解决的关键科学问题之一是催化剂电子结构与表面吸附、反应活性的关联,这也是精准设计和调控催化材料结构的科学基础.本文针对喹啉选择加氢制1,2,3,4-四氢喹啉反应,利用有序介孔炭载配位结构精准调控的金钯合金作为模型催化剂,描述金属d带电子特性对吸附络合物稳定化能和催化活性的决定性影响.在原子水平确认了AuPd纳米催化剂具有与块体合金相似的结构,排除了共晶体或核壳结构.通过调变Au和Pd浓度,实现了Pd-Au配位数由0-8的精确调控.利用XPS和XANES测量了Pd位点上d带和非dsp,主要是s)带得失电子数,将其定义为d带电荷密度(d charge),发现其与Au-Pd配位数密切相关.在喹啉选择加氢反应中,Au50Pd50和Au67Pd33催化剂的活性最高,与纯Pd催化剂相比,活性提高了11倍.进一步将合金电子结构与反应活性相关联发现,合金催化剂中Pd位点d带电荷密度与反应物生成活性络合物的活化熵(ΔS0*)和转化频率(TOFPd)成线性相关.有序介孔炭载AuPd合金催化剂具有很高的稳定性,在反应中未见金属表面偏析、聚集和流失.金属位点d带电荷密度描述符为含强吸附质反应中高性能催化剂的设计提供了新原理.
关键词d带电荷密度    金钯合金    活化熵    喹啉    加氢    

1 Introduction

Metal catalysts play an important role in petrochemical, pharmaceutical, and fine chemical synthesis, etc [1]. The challenge is to identify descriptors of catalytic activity and discover how they determine the activity and selectivity of a catalytic material [2]. The Sabatier principle is primarily used to optimize the catalyst with an intermediate binding strength between the catalyst and atoms or molecules at the surface. Similar volcano-shaped curves can be obtained empirically relating the bond energy for a reaction with the position of the transition metal in the periodic table, and with the catalytic activity of the metal [3]. Therefore, determining a descriptor to represent the "bond strength" is important. The heat of formation of the metal oxide [4, 5], the "d-hole" [6, 7], or the number of d-orbital electrons [8, 9] in the early literature has been often used to understand these volcano curves.

In recent years, theoretical chemists have proposed descriptors based on theoretical models in an effort to correlate the electronic structure of the catalyst with the catalytic activity [10, 11]. Nörskov et al. [12, 13] qualitatively described the adsorption energy of adsorbate molecules on the metal surface by using the d-band center of the transition metal surface, because the properties of the chemical bond for adsorption are intrinsically related to the electronic structure of the transition metal surface. The adsorption energy of a simple adsorbed species or the enthalpy of a rate-determining step reflects the binding energy between the adsorbate and the surface of the catalyst, and should be located near the peak position of the volcano curve. The effectiveness of this theoretical descriptor has been verified in several simple reactions, illustrating the possibility of the rational design and optimization of catalysts for the related reactions. For example, on the basis of the Brönsted-Evans-Polanyi relations and adsorption energy correlations, the descriptor for the methanation reaction is the dissociation chemisorption energy of CO [14]. The resulting mode is in good quantitative agreement with experiment, showing a well-behaved volcano curve relating the experimentally measured activity to the dissociation energy, and Ru and Co are the most active transition metals [15]. Some Fe-Ni alloys which show a higher methanation activity, with a lower price of the constituents, is predicted [16].

The limitations for these descriptors are their application in a real and complex catalytic reaction system, and the difficulty in quantitative measuring them with experiment [17]. Therefore, for a fine chemical synthesis system, exploring the nanocatalyst activity descriptors that are convenient for experimental measurement has more important practical significance for the rational design of the catalyst. Very recently, we have reported the experimentally measured surface d charge at Pd sites as a descriptor that quantitatively correlates with the adsorption strength and catalytic activity of a material for the oxidation of benzyl alcohol [18]. The validity of this descriptor needs to be extended.

The catalytic hydrogenation of N-heterocyclic compounds is of fundamental interest in organic synthesis [19]. For example, 1, 2, 3, 4-tetrahydroquinolines (py-THQ) are ubiquitous in numerous biologically active natural products and pharmacologically relevant therapeutic agents [20, 21]. But their hydrogenation remains a scientific and technological challenge due to the potential poisoning of the catalyst by either the substrate or their reduced products [22]. The electronic structure would significantly change the bond strength between metal and adsorbates and adsorption configuration [23, 24]. For example, an amine-containing support with electron-donating properties offers an electron-rich surface of Pd nanoparticles. Catalysts prefer the adsorption of quinoline over py-THQ which in turn enhances the catalytic activity and selectivity [25]. To the best of our knowledge, this relationship has not been measured.

Here an attempt is made to correlate the experimentally measurable surface d charge at Pd sites in the AuPd nanocatalysts to the activation entropy and catalytic activity for the selective hydrogenation of quinoline. A series of AuPd nanoalloys with different Pd concentrations has been synthesized which show a volcano-shape curve for the surface d-charge gain at Pd sites with Pd concentration. The maximum occurs at 33–50 at% Pd. The activation entropy and catalytic activity follow a similar volcano-shape curve. The good qualitative agreement between the surface d charge at Pd sites and the activation entropy and catalytic activity paves the way for an understanding of the trends for these complex reactions using experimentally measurable descriptors, and further optimizing and designing novel catalysts.

2 Experimental
2.1 Synthesis of AuPd nanocatalysts supported on ordered mesoporous carbon

Ordered mesoporous carbon-supported AuPd nanoalloys were synthesized by a coordination-assisted method reported earlier using a commercial triblock copolymer (F127) as a structure-directing agent, thiol-containing silane (MPTMS) as a coordination agent, chloroauric acid and palladium chloride as metal sources, phenolic resins as the carbon source, and tetraethylorthosilicate (TEOS) as the silicon source [18]. In a typical synthesis, 5.2 g TEOS (Sigma-Aldrich, > 99 wt%), 4.9 g MPTMS (Acros Chemical Inc, 85 wt%), 5.0 g HCl (0.2 M) and 25.0 g ethanol were mixed for 30 min. Then, 8.0 g F127 (Acros Chemical Inc, EO106PO70EO106, MW= 12600 g mol‒1), 20.5 mL PdCl2 (palladium concentration in ethanol: 56.4 mmol L‒1) and 9.0 mL HAuCl4 (gold concentration in ethanol: 24.3 mmol L‒1), and 25.0 g phenolic resins, which were prepared according to the established procedure as described in Supporting Information (SI), were added to the solution in that order. The mixed solution was stirred for 2 h at 40 ℃ and was then poured into several dishes. After continuous evaporation of the ethanol in an oven at 40 ℃ for 6 h followed by a further 100 ℃ thermopolymerization in an oven for 24 h, light yellow transparent films were scraped from the dishes and were refluxed in sulfuric acid with mechanical stirring (1.0 g of solid per 100 mL of 50 wt% sulfuric acid) at 90 ℃ for 24 h. After filtration, washing with distilled water, and drying at 80 ℃ under vacuum overnight, yellow solid materials were obtained. These materials were finally calcined under a nitrogen flow with heating rates of 1 ℃ min-1 from room temperature to 350 ℃ and then 5 ℃ min‒1 to 600 ℃, where it was kept for 3 h. The obtained materials were denoted Au50Pd50, the subscripts representing the respective Au/(Au + Pd) and Pd/(Au + Pd) values in molar percentages in the mother liquor.

The total amount of metals is ~3.3 wt% for the catalyst and both the amounts and the Au:Pd ratios are close to the theoretical values. The synthesis of catalysts with different Au:Pd ratios was identical except for the use of different amounts of chloroauric acid and palladium chloride. Monometallic Au and Pd catalysts were synthesized by a similar method but using only one of the precursors.

2.2 Catalytic tests

The quinoline hydrogenation reactions were carried out in a 25 mL autoclave equipped with a Teflon tube (Parr). 5 mL of Milli-Q water, 0.5 mmol of quinoline (Admas-beta, 99 wt%) and 20 mg of the Au50Pd50 catalyst were added to the autoclave. The catalysts were treated at 180 ℃ in an oven for 6 h before use. The autoclave was purged at least three times with hydrogen and the H2 pressure was set to 1.0 MPa. The stirring speed was fixed at 800 rpm to exclude diffusion limitations, and the temperature was kept at 80 ℃. After each reaction, the catalyst was separated by hot filtration. The filtrate was extracted with 30 mL ethyl acetate and analyzed using an Agilent 7890B gas chromatograph (GC) equipped with a DB-1 capillary column and a flame ionization detector. The solid catalyst was recovered for reuse by washing with water and drying at 80 ℃ overnight under vacuum. Each test was repeated at least three times with an experimental error of ± 5%. The catalytic results were given in terms of the conversion of quinoline, yield of py-THQ, 5, 6, 7, 8-tetrahydroquinoline (bz-THQ) and decahydroquinoline (DHQ), initial reaction rate (r0), turn-over frequency (TOF). The estimated TOF on the basis of exposed Pd-involving sites (TOFPd) was taken to eliminate the effect of Au dilution. The calculations for TOFPd, activation energy (Ea), and entropy of activation (ΔS0*) were described in SI. The TOFPd value was calculated at a conversion below 20% and was reproducible to within ± 5%.

A soluble Pd trapping experiment using thiol group-modified mesoporous silica (SH-SBA-15) was performed to test the stability of the catalysts. SH-SBA-15 was added to the reaction with the substrate and catalyst (S:Pd molar ratio 35:1). The synthesis details of the SH-SBA-15 are given in the SI.

The recycled hydrogenation reaction was carried out over used catalyst under the reaction conditions described above. After the sixth run the catalyst was referred to as Au50Pd50-R. In order to keep the same amount of catalyst in each run, several parallel reactions were carried out at the same time.

2.3 Characterization

N2 adsorption-desorption isotherms were measured at 77 K using a Micromeritics TriStar II 3020 analyzer. The material was vacuum dried at 385 K for 8 h. The specific surface area (ABET) of the material was obtained by the Brunauer-Emmett-Teller (BET) method; the pore volume (V) and pore size (D) of the material were calculated using the Barrett-Joyner-Halenda (BJH) model. X-ray diffraction (XRD) measurements were performed on a Rigaku D max-3C diffractometer using Cu Kα radiation (40 kV, 20 mA, λ = 0.15408 nm). According to the Scherrer formula, the metal (111) interlayer spacing was calculated as d = λ/2sinθ, based on the (111) diffraction peak in the wide-angle XRD pattern. Transmission electron microscopy (TEM) images were recorded on a JEM 2100 microscope operating at 200 kV. A JEM 2100 microscope and a FEI Tecnai G2 F30 microscope operating at 200 kV were used for high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) imaging. The energy-dispersive X-ray spectroscopy (EDX) spectrum was obtained in the STEM mode with a focused electron beam in the sub-nanometer range. A Varian VISTAMPX inductively coupled plasma-atomic emission spectrometry (ICP-AES) was used to determine the metal content. X-ray photoelectron spectroscopy (XPS) was performed on a Perkin-Elmer PHI 5000 CESCA model instrument. The measured XPS data was calibrated with a carbon C 1s binding energy of 284.6 eV. X-ray absorption spectra of the Au L3-edge (11919 eV) and Pd K-edge (24350 eV) of samples were collected on a BL14W1 of the Shanghai Synchrotron Radiation Facility (SSRF). The data processing was performed using the ATHENA program. The Au and Pd extended X-ray absorption fine structure (EXAFS) data of the sample were fitted by a fast Fourier inverse transform (IFEFFIT) in R space using a software package.

3 Results and discussion
3.1 Catalytic performance

The selective hydrogenation of quinoline to py-THQ was first conducted over pristine palladium and gold catalysts (Fig. 1). The maximum conversion of quinoline was as low as 5% over the gold catalyst, indicating a rather low activity. This is in agreement with both theoretical calculations and experiment. High hydrogen dissociation energies were observed over the Au (111) surface [26]. Hydrogenation reactions catalyzed by Au typically occur under demanding conditions, for example, supported on reducible oxides including TiO2 and Fe2O3 and pre-adsorption of the quinoline, which enhances the dissociation of H2 [20]. An increased activity was observed for the Pd catalyst, but the selectivity was quite low, especially at high conversions. It has been reported that by using Pd/AC as the catalyst, only 81.7% selectivity to py-THQ was obtained with the byproduct being bz-THQ [27]. With the Au50Pd50 nanoalloy catalyst, significant quinoline conversion was obtained. Remarkably, no byproducts are detectable.

Fig. 1. (a) Conversion plots for the selective hydrogenation of quinoline over AuPd nanocatalysts with different Pd concentrations; and (b) selectivity of py-THQ along with conversion. For comparison, the results for monometallic Au and Pd are also provided. Reaction conditions: 20 mg of catalyst, 0.5 mmol of substrate, 5 mL of water, 80 ℃, and 1 MPa H2.

The leaching of dissolved metal ions is first excluded since these ions have been long argued as catalysing hydrogenation in solution, and the stability of the catalysts has been studied. A bare thiol-containing solid has been used to quench quinoline hydrogenation by soluble molecular metallic species. For example, Jones et al. found that catalysis by soluble Pd species for Suzuki and Mizoroki-Heck coupling reactions is quenched by the addition of SH-SiO2 either at the initial stage or during the reaction [28]. We have previously demonstrated SH-SiO2 as an effective and versatile poison compared with poly(4-vinylpyridine) and thiourea-based resin QuadraPure® TU, in quenching any catalysis by leached Pd and gold species [29, 30]. No significant difference in the kinetic profile of AuPd nanoalloys was observed in the reactions in the presence/absence of SH-SBA-15 with an S:Pd molar ratio of 35, confirming that catalysis occurs on the Pd surface rather than in solution (Fig. S1). The stability of the catalyst was studied in more detail. Quinoline conversion in the initial phase of the reaction and the end of the reaction show no substantial changes (Fig. 2) in repeated uses. These results indicate the negligible leaching of metal ions and little aggregation of the Pd particles. In addition, the strong adsorption of the reactants and products and their accumulation on the Pd surface is possibly inhibited. The solution after each run was collected, and the concentration of metal ions was below the detection limit, confirming that the catalyst has not leached [28]. These results guarantee that this is purely an investigation of surface catalysis and surface metal concentrations are almost unchanged after the reaction.

Fig. 2. The initial reaction rate (r0) and the turn over number (TON) for the Au50Pd50 catalyst in successive cycles for the selective hydrogenation of quinoline. Reaction conditions: 20 mg of catalyst, 0.5 mmol of substrate, 5 mL of water, 80 ℃, and 1 MPa H2.

To ensure a kinetic regime, the limitations of diffusion and mass transfer were eliminated. The experiments were carried out at different stirring rates and the results showed that the reaction rate remained constant in the range of 600–1000 rpm. The Madon-Boudart test was also conducted [31]. Au50Pd50 catalysts with different metal loadings, ranging from 1.2 to 2.0 wt%, were synthesized (Figs. S2 and S3). The reaction rate followed a linear relationship with the metal content in the catalyst, demonstrating that there was no mass transfer effect (Fig. 3). The metal nanoparticles were easily accessible by the reactants.

Fig. 3. The Madon-Boudart test on the reaction rate of quinoline hydrogenation using the Au50Pd50 catalyst with different metal concentrations but similar metal dispersions.

Figure 1(a) shows the variation of quinoline conversion for AuPd with different Pd concentrations. As the Pd content in the catalyst decreased from 100 at% to 80 at%, 67 at% and 50 at%, the quinoline conversion increased from 12.71% to 22.92%, 26.60% and 50.64% for one hour reactions. However, decreasing the Pd concentration to 33 at% did not result in an increase of the quinoline conversion. The quinoline conversion was reduced to approximately 36.75% and 12.50% for Au67Pd33 and Au80Pd20 for one hour reactions, respectively. A similar relationship has been observed between the activity and Pd concentration in the selective oxidation of aromatic alcohols [18]. Byproducts were detected over Au20Pd80 and Pd during the reaction. In contrast, the selectivity to py-THQ was extremely high over the Au80Pd20, Au67Pd33, Au50Pd50 and Au33Pd67 catalysts with no detectable byproducts. This result is analogous to the N-modified Pd and thermally-reduced Pd supported on carbon nanotubes [32].

It is noted that the conversion plot is nearly linear during the initial period of the reaction over the catalysts studied, similar to the kinetics for the oxidation of benzyl alcohol over AuPd nanoalloys, suggesting pseudo zero-order reaction kinetics [18]. In addition, a change of quinoline concentration in the range of 0.05–0.2 mmol mL‒1 results in an almost unchanged initial reaction rate over Au50Pd50 (Fig. S4). These phenomena imply that the reaction order for quinoline is pseudo zero order. Since the gold catalyst shows relatively low activity in hydrogenation, we simplify the active sites to Pd [20]. The reaction constant (kPd) can be deduced from r0, and an Arrhenius plot can be plotted from which Ea and pre-exponential factors (A) can be estimated (Fig. S5). As shown in Fig. 4, the activation energies follow a reverse volcano-like curve as a function of percent Pd with the minimum at Au50Pd50. The apparent entropy change (ΔS0*) reflects the activation step of the chemical reaction in the thermodynamic formulation of the rate constant and would be influenced by changes in the available electronic energy states. It has been reported that the electronic energy states may promote an electron to a localized band for adsorption of a reaction intermediate [33]. ΔS0* follows a volcano-shaped curve, where the maximum ΔS0* is also observed for the 50 at% Au catalyst. A linear relationship is observed between Ea and ΔS0* (Fig. S6), which is called the compensation effect. This effect can perhaps be explained by the change in the freedom of the system, due to the loosened bonds between the surface atoms and the adsorbates as the energy of the system increases [34]. These observations strikingly demonstrate the influence of the electronic energy states of the metal on the adsorption, and in turn on the catalytic properties.

Fig. 4. Relationship between the activation energy (Ea) and entropy of activation (ΔS0*) and the Pd concentration in AuPd nanocatalysts for quinoline hydrogenation.
3.2 Catalyst characterization

It is most important that the interpretation of catalytic results for metal alloys takes into account the possibility of phase separation [35]. TEM images of bimetallic Au-Pd catalysts supported by a carbonaceous matrix show that the metal nanoparticles are homogeneously dispersed in the matrix with an average size of approximately 3.0 nm, regardless of the metal composition (Figs. 5, 6, S7 and S8). Therefore, the size effect for the bimetallic AuPd nanoparticles on the catalytic behavior can be excluded. Hexagonally, and linearly arranged pores in large domains in the absence of distinct defects confirm the ordered 2D hexagonal mesostructure with p6m symmetry for all catalysts. HAADF-STEM images accompanied by EDX patterns were taken. The Au:Pd ratio is calculated to be close to the value of the mother liquor used for their synthesis (Table S1). Using the EDX spectra obtained with a focused electron beam in the sub-nanometer range in the STEM mode, the distributions of gold and palladium can be seen. Representative line profiles for Au67Pd33 and Au50Pd50 as examples of the different compositions in a single particle are shown in Figs. 5(d) and 6(d). These lines match well from the periphery to the centre, demonstrating that these bimetallic catalysts have uniform compositions throughout the particles with the Au-to-Pd molar ratio of approximately 2 and 1, respectively [36]. The homogeneity of the Pd distribution on the surface instead of a Pd rich phase confirms the uniform structure of the alloy.

Fig. 5. TEM images viewed from the (110) (a) and (001) (b) directions for the Au67Pd33 catalyst; HAADF-STEM image (c) and line profiles in the EDX pattern (d) of a single nanoparticle. The inset in (b) is the particle size distribution.
Fig. 6. TEM images viewed from the (110) (a) and (001) (b) directions for the Au50Pd50 catalyst; HAADF-STEM image (c) and line profiles in the EDX pattern (d) of a single nanoparticle. Inset in (b) is the particle size distribution.

Small-angle XRD (SAXRD) patterns show well-resolved diffraction peaks belonging to a 2D hexagonal mesostructure with a lattice constant of approximately 9.0 nm (Fig. 7(a)), regardless of the metal composition, in agreement with the TEM results. Typical type-IV nitrogen sorption isotherms are detected for all the catalysts studied, and are characteristic of mesoporous solids with a uniform pore size (Fig. S9). The asymmetric hysteresis loops are attributed to H2 type, different from the H1 type for a pristine ordered mesoporous silica/carbon composite (OMSC), indicates the non-uniformity of the mesopores originating from the intercalated AuPd nanoparticles. The calculated pore-size distribution curves are narrow with a most probable value of approximately 3.6 nm. The pore volumes and the BET surface areas are close to those for pristine OMSC and are approximately 0.23 cm3 g‒1 and 340 m2 g‒1 (Table S1). These phenomena indicate that monometallic or bimetallic nanoparticles show no obvious effect on the mesostructure.

Fig. 7. SAXRD (a) and WAXRD (b) patterns with a high resolution fine scans for AuPd nanocatalysts with different Pd concentrations. For comparison, the results for monometallic Au and Pd are also provided. Au50Pd50-R represents the re-used catalyst after six catalysis cycles.

Wide-angle XRD (WAXRD) patterns show wide diffraction peaks at approximately 23° which are attributed to amorphous carbon and silica, and the diffuse peak at approximately 39° is attributed to the AuPd alloy with a very small particle size. Higher resolution scans were also collected in the vicinity of the diffraction peak at approximately 39° (Fig. 7(b)) and show that the diffraction peak is symmetric, which is characteristic of a nanoalloy. A clear shift of the characteristic diffraction peak for AuPd alloys compared to monometallic Au and Pd is observed, and the peak is located between them. For example, the Au50Pd50 sample has a cell constant of 0.4017 nm, giving further evidence on the alloy structure. As a result, the formation of a monometallic core-shell structure or macroscopic phase separation is excluded. The elemental compositions, estimated using Vegard's law in WAXRD patterns, are close to the theoretical values for dilute Pd catalysts, but show the distinct deviation for Pd-rich ones, due to an obviously larger lattice constant compared to the alloy having complete solid solubility. This phenomenon implies that an ideal alloy is formed for dilute Pd catalysts, while gold segregation in micro-domains occurs for dilute Au catalysts.

Typical doublets of the Au 4f, and Pd 3d core level bands are detected for monometallic or bimetallic catalysts (Fig. 8). The doublets for Au 4f7/2 and 4f5/2 can be extensively deconvoluted, with the fitted dominant peaks at 84.30 and 88.00 eV corresponding to Au0. With increased Pd concentration, the Au 4f core levels show a gradual negative shift for AuPd nanoalloys. Similar core-level shifts have been reported for either bulk AuPd alloys or nanoparticles, which is ascribed to the Au-Pd alloying interaction and increased number of Au-Pd bonds [37-40]. A relatively large shift occurs when mixing Pd with Au, and the shift then undergoes an extremely slow change in dilute Au alloys, implying that the Au-Pd coordination number increases according to the Pd concentration in dilute Pd alloys, but may show a slight increase or even be constant in dilute Au alloys. This result is in good agreement with the XRD results, thus confirming the alloy structure.

Fig. 8. XPS spectra of the 4f level of Au (a) and the 3d level of Pd (b) for AuPd nanocatalysts with different Pd concentrations. For comparison, the spectra for monometallic Au and Pd are also given. Detailed peak fitting is given in the SI.

Since a distinct overlap is observed between the Pd 3d5/2 and the Au 4d5/2 components for the bimetallic catalysts, the Au 4d5/2 intensity was calculated from the well-resolved Au 4f5/2 intensity, and this value was subtracted from the above overlapped peak to determine the Pd 3d5/2 intensity. The extensively deconvoluted doublet in the Pd 3d region for the monometallic Pd catalyst shows that the fitted major peak is due to metallic Pd, and the minor one to a Pd–O coordination [30, 41]. It is noted that the Pd 3d core levels for AuPd nanoalloys show a gradual negative shift with decreased Pd loading due to the Au-Pd alloying interaction, and for dilute Pd alloys is relatively large. Considering the atomically dispersed Pd0 species in these alloys, which are mainly responsible for the significant XPS positive chemical shift of 1.4 eV for Pd 3d lines due to their smaller electron relaxation energy [42], the negative shifts for the Pd 3d core-level bands in dilute Pd nanoalloys are even more dominant. It is interesting that both the Au 4f and Pd 3d core level bands in bimetallic species show obvious negative shifts compared to the corresponding monometallic species. These phenomena are in accordance with bulk AuPd alloys [43], but in contrast to Cu–Pd bulk alloys [44], where the Pd shifts are positive and those for Cu are negative. Medlin et al. [45] also observed a shift of the Au 4f core level to a lower binding energy in a AuPd/C nanoalloy catalyst compared to monometallic Au/C, and agreed with that charge transfer may not be the correct explanation for the core-level shifts. Therefore, the hybridization changes between Pd and Au are the possible reason that the occupied valence shells redistribute. The estimated Au:Pd ratios in alloys with different Pd concentrations are similar to the theoretical values further confirming the alloy structure, and the surface adopting similar Pd-Au coordination with the bulk (Table S1).

The above results are in good agreement with the previous reports for AuPd nanoalloy catalysts. As a consequence, we used the earlier XAFS spectra for the present catalysts. The ratio of CNAu–Au to CNAu–Pd, generated from Au L3 k3-weighted and Fourier transformed EXAFS fittings (Table S2), is very close to the Au-to-Pd molar ratio for Au-rich nanoalloys (including Au:Pd = 1), indicative of the formation of homogenous alloy structures. A deviation occurs for dilute Au catalysts, for example, the CNAu–Au to CNAu–Pd ratios for both Au33Pd67 and Au20Pd80 are much larger than the Au-to-Pd molar ratios. These results coincidently indicate a slightly different structure for dilute Au alloys compared to dilute Pd-alloys. Gold clusters are more likely to be formed in dilute gold alloys. A calculation from the EXAFS spectra of the AuPd nanoparticles near the Pd K-edge gives a reduction of CNPd–Pd from 2.8 to 0. The Pd–O contribution is dominant in Pd-rich catalysts, demonstrating that extremely small Pd clusters are formed when a small amount of Au is dissolved in Pd. Considering the surface Pd–O and Pd–Pd coordination, the coordination number for Pd clusters to that for Pd–Au is close to the Pd/Au ratio, implying a random alloy. When the Au-to-Pd ratio exceeds 2, the peaks to Pd–Pd and Pd–O shells almost disappear. The XAFS features show high similarity to AuPd nanoalloy phase [46]. Contiguous Pd clusters disappear with increasing Au loading, isolated Pd ensembles form, and eventually each Pd atom is separated by Au. An Au-rich core and a Pd-rich shell can also be excluded because Pd–Pd coordination is clearly detected for this structure when the Pd concentration is greater than 20 at% [47, 48].

The microstructural discrepancies related to the metal concentrations are caused by surface energy versus surface mixing energy, determining whether the two components mix or form islands. The segregation energy of Au–Pd is 0.15 eV for Au-host and ‒0.14 eV for Pd-host alloys [49]. As a result, Pd thermodynamically prefers a random-distribution, and individual Pd atoms form a complete solid solution with the Au-host. While Au prefers to dissolve in Pd as small Au islands [49, 50]. These results are consistence with the previous report, and suggest the existence of a complete range of stable solid solutions.

3.3 d charge descriptor

Band theory indicates that in a face-centered cubic (fcc) transition metal alloy, the narrow eg band with a high density of states will be more readily filled than the t2g state by the introduction of s electrons from the non-transition metal present in the alloy [33]. The complete filling of d-bands at a specified alloy composition results a pronounced change in the available electronic energy states.

XANES spectra of the Au L3-edge for the Au foil clearly shows the edge that probes the transition from the 2p3/2 to 5d5/2 and 5d3/2 dipole-allowed transitions, resulting in a white line at 11925 eV whose intensity reflects the number of holes in the d-band of the Au atoms (Fig. S10). On one hand, the nanosize effect causes a decrease in the line intensity because of the smaller number of atoms forming a lattice, and therefore a decrease in the d-band width [51]. A comparison of the white line intensity for a 3 nm Au catalyst with that of a Au foil, indicates an increase in the d electron count of the Au atoms in the former. Since dilute Pd alloys have similar Au coordination numbers to a pure Au catalyst, the almost loss of the white line intensity when mixing Pd with Au indicates the alloy structure. The Au L3-edge resonance is insensitive to the alloy composition, and the involvement of 20% Pd is sufficient to deplete unoccupied Au 5d5/2 states near the Fermi level. The presence of Pd has a dominant effect on the filling of the Au d-band (an increase of the d-electron density). A combination of the hybridization changes between Pd and Au is indicated by the XPS results, which can be attributed to a change in d-band hybridization such as a decreased band width and spin-orbit splitting and is not due to charge transfer, in agreement with the fact that the extent of charge transfer is small for Pd-Au alloys [52]. It has been reported that the Au-Pd alloying interaction produces a more pronounced change of the Au d electron density than the size effect does [53]. A dominant increase in the intensity of the second feature after the edge is observed with the increase of Pd content. This change is caused by the change in the distance between the two types of atom which causes an intra-atomic redistribution of charge [37, 54]. As a result, the gold is in the metallic state, and the Pd–Au d-d interaction is believed to have an important effect on the electronic structure of the alloy. In particular, with increasing Pd content the Au 5d5/2 partial density of states is broadened and shifted to a higher binding energy relative to the Fermi level [52].

The Pd XANES spectrum at the K-edge of the Pd foil exhibits a pronounced white line due to the unfilled Pd d-band, and the K-edge of the PdO exhibits an even more dominant white line with a shift of the threshold energy to a higher energy because of the higher oxidation state of Pd in PdO (Fig. S11). The hybridization-mediated 1s → 4d, dp (a) absorption transition in the pre-edge region is unresolved. The near-edge spectrum shows two resonance peaks due to 1s → 5p (p) and 1s → 4f (f)transitions. The shape of the second absorption edge reflects the extent of 4d-5p hybridization [55]. This feature is insignificant for the PdO reference material. Similar features are observed in the spectrum of the monometallic Pd nanocatalyst and in that of the Pd foil, showing a similar coordination symmetry to the Pd metal of fcc structure. However, the post-edge threshold energy shifts to a slightly higher energy for the monometallic Pd catalyst. The possible reasons for this lie in the tiny Pd particles having a Pd metal fcc structure with very short range order or the coordination of Pd with a light O atom [56]. Upon dissolution of 20 at% and 33 at% Au into Pd nanoparticles, three differences can be seen. First, a further shift of the white line energy to a higher energy is observed, in nearly identical fashion to the absorption edge of PdO, in good agreement with the increase of Pd–O coordination. Second, a positive shift in the 1s → 5p (p) transition accompanied by signal broadening is detected, which is possibly due to an electronic perturbation of the Pd atoms going from a Pd-like environment in the monometallic structure to the AuPd alloy phase. The Pd d-band becomes occupied [57]. Third, the 1s → 4f (f) transition shifts to a lower energy and its intensity is significantly reduced. This shift can be interpreted as due to an increase in the atomic separation and the 4d-5p hybridization that has been observed for the local distortion of the Pd lattice by interstitial C [58] and AuPd bimetallic nanoparticles [59]. A further decrease in Pd concentration to 50% (Au50Pd50) leads to no further shift, implying that the 4d band is filled [60]. By comparison, both the white line and p and f transitions in dilute Pd nanoalloys shift towards that for a Pd foil. As a result, the orbital character of the screening charge may change [61].

In the present AuPd nanoalloy, for a greater than 50 at% Pd content the Pd d-band becomes occupied when alloyed with Au. With a further reduction of Pd concentration, it may be conjectured that the 5sp electrons of the surface Pd atoms lose some d-character possibly due to the increased Pd–Pd bond distances and isolated Pd atoms. The d-d separation would narrow the 4d band of surface Pd sites effectively pushing its top below the Fermi level. For Au-rich alloys with a low density of states the screening efficiency will be lower [62]. Due to strong Au–Pd d-d interactions, the orbital character of the unoccupied states at the Fermi energy changes from predominantly d to predominantly sp-type as the Pd concentration in the alloy decreases [61].

Non-d- and d-charge transfer onto atom sites can be estimated based on the XPS and XANES results (Table S3) [63]. The binding energy shift, ΔEB, between the pure metal and the alloy measured relative to the Fermi level in XPS experiments is related to the change of the Hartree-Fock one-electron energy of the core level (Δϵ(i)), the change in the work function (ΔФ), and the final relaxation state of a core hole (ΔEγ) [64]. Δϵ(i) is directly related to the charge transfer in the initial ground state, and is caused by the effects of the Coulomb repulsion energy between the core electron and the valence electron which arises from the redistribution of valence electrons, and the change in the number of the valence electrons which is associated with the charge changes on other lattice sites due to the total charge transfer in or out of the parent atomic site. A theoretical calculation on alloy systems based on the pseudopotential linear response method, shows that in the presence of the core hole △Eγ can be as large as 2 eV, which is the same order of magnitude as the initial-state effect Δϵ(i) and therefore cannot be neglected [65]. The estimate is in good agreement with the above discussion. When adding Au to Pd, the alloy has a very small amount of charge transfer (< 0.1 electron). The d charge transfer at the Pd sites is small and the compensation comes from the gaining of 0.1–0.2 sp electrons at Au sites. The decrease of Pd concentration results in a pronounced change in the available electronic energy states, and an increase in the number of d electrons at Pd sites. At approximately 33 at%–50 at% Pd the alloy contains Pd atom pairs or isolated Pd atoms, and the d-electron gain reaches the maximum value of approximately 0.28 e. For a further decrease of Pd content the total charge transfer is large, but the d-electron gain at Pd sites is reduced to some extent. From a kinetic point of view, the measured apparent entropy change is the difference between the entropy of the transition state and the reactant. The magnitude of the entropy value, in general, shows that the activated complex is an immobile molecule attached to the catalyst surface [66, 67]. As a result, the density of states close to the Fermi level of the transition metal has a large effect on ΔS0* because of its effect on the bound state of the absorbed intermediate [68]. On the other hand, the d-band theory proposes that for a given metal and adsorbent, the center of themetal d states relative to the Fermi level determines the binding energy of the adsorbent to the surface. Inspired by these works [69], we therefore attempted to relate the increased d charge at the Pd sites to the activation entropy, ΔS0*, and the TOFPd value (Fig. 9), and both showed a linear relationship. These results indicate that the Pd charge of the Group VIII-IB nanoalloy is fundamental to the adsorbate binding strength in the surface reaction. In addition, the increase in ΔS0*, which is related to the adsorption strength, is highly related to the increase in the reaction rate. Accordingly, the d charge at surface Pd sites can serve as a descriptor, which can correlate well with the catalytic activity. It should be mentioned that we simplify the relationship. On one hand, the surface changes from Pd clusters with Pd–Pd coordination, to Pd islands with Pd–Pd and Pd–Au coordination, and finally to isolated Pd with Pd–Au coordination. Comparing the activity of Au67Pd33 to Au80Pd20 both of which possess isolated Pd, the TOFPd value is higher for the former than the latter. An overloading of Au which reduces the d-electron gain leads to a decreasing activity. As a result, we simplify the effect by electronic properties instead of structures. On the other hand, the present study is limited to a strong-adsorbate-involved reaction. Quinoline is reported to adapt a di-bridge adsorption configuration on Pd cluster with a high adsorption energy[70]. Both the adsorption strength and configuration are significantly affected by the electronic properties, in particular, d-electron. A linear relationship has been found to correlate the adsorption strength and configuration, here represented by activation entropy, with the d-electron gain at Pd sites. Taking account of the complexity for surface reaction on nanoparticles, other descriptors should also play important roles. The factors are still being investigated in our lab, and deserve to be reported separately.

Fig. 9. Relationship between the d-charge gain at the Pd sites of AuPd nanoalloy catalysts and the TOFPd (red line) of AuPd nanoalloy catalysts, and the entropy of activation (ΔS0*, blue line) for the selective hydrogenation of quinoline.

The absolute thermoelectric power, S [71] of AuPd alloys is found to be highly sensitive to the presence of unoccupied d-band states. A sudden drop of the S value is expected as the d-band is filled. The results on bulk AuPd alloys indicated that the d-band filling point was very close to 50 at% Au. In the present ordered mesoporous carbon-confined AuPd nanocatalysts, the d-band filling tendency is in good agreement with that for bulk AuPd alloy. We therefore attempted to relate the thermoelectric power S to the activation entropy ΔS0* (Fig. S12). The entropy change indeed follows a similar trend with alloy composition as does the thermoelectric power, which shows a minimum at 50 at% Au. This indicates that the d-band energy position in the Group VIII-IB nanoalloy is fundamental to the activation entropy associated with the adsorbate binding strength in the surface reaction, similar to that for transition metals [72, 73]. An increased reaction rate for the "true" nanoalloy would ordinarily be expected from the decrease in Ea and the compensation effect by the decrease in ΔS0* which is related to a weaker adsorption strength.

4 Conclusions

Experimentally measurable surface d charge at Pd sites has been proposed to serve as a descriptor that quantitatively correlates to the adsorption strength and catalytic activity of AuPd alloys for the selective hydrogenation of quinoline. AuPd nanoalloys with different Pd concentrations supported on ordered mesoporous carbon have been synthesized and used as model catalysts. These nanoalloys have tunable Au–Pd, Au–Au and Pd–Pd (Pd–O) coordination numbers, and therefore tunable electronic properties. Hybridization changes between Pd and Au are the possible reason for changing the electronic properties, which in turn have a great influence on the adsorption energy. Non-d- and d-charge transfer onto the Pd atom sites has been estimated. The d-electron gain on the Pd surface has a linear relationship with the activation entropy and catalytic activity. The maximum d-charge gain at the Pd sites occurs at Pd concentrations of 33 at%–50 at%, with the alloy showing an almost 11 times higher activity than a pure Pd catalyst.

References
[1]
L. T. Mik, E. Cséfalvay, Á Cséfalvay, Chem. Rev., 2018, 118, 505-613. DOI:10.1021/acs.chemrev.7b00395
[2]
A. J. Medford, A. Vojvodic, J. S. HummelshØj, J. Voss, F. Abild-Pedersen, F. Studt, T. Bligaard, A. Nilsson, J. K. Nörskov, J. Catal., 2015, 328, 36-42. DOI:10.1016/j.jcat.2014.12.033
[3]
M. Che, Catal. Today, 2013, 218–219, 162–171.
[4]
S. J. Henley, J. D. Carey, S. R. P. Silva, Phys. Rev. B, 2005, 72, 195408.. DOI:10.1103/PhysRevB.72.195408
[5]
D. Walsh, L. Arcelli, T. Ikoma, J. Tanaka, S. Mann, Nat. Mater., 2003, 2, 386-390. DOI:10.1038/nmat903
[6]
D. R. Penn, Phys. Rev. Lett., 1979, 42, 921-925. DOI:10.1103/PhysRevLett.42.921
[7]
F. C. Zhang, T. M. Rice, Phys. Rev. B, 1988, 37, 3759-3761. DOI:10.1103/PhysRevB.37.3759
[8]
J. I. Musher, J. Am. Chem. Soc., 1972, 94, 1370-1371. DOI:10.1021/ja00759a063
[9]
J. Kanamori, Prog. Theor. Phys., 1963, 30, 275-289. DOI:10.1143/PTP.30.275
[10]
W. Kohn, A. D. Becke, R. G. Parr, J. Phys. Chem., 1996, 100, 12974-12980. DOI:10.1021/jp960669l
[11]
J. Greeley, J. K. Nörskov, M. Mavrikakis, Annu. Rev. Phys. Chem., 2002, 53, 319-348. DOI:10.1146/annurev.physchem.53.100301.131630
[12]
J. K. Nörskov, T. Bligaard, J. Rossmeisl, C. H. Christensen, Nat. Chem., 2009, 1, 37-46. DOI:10.1038/nchem.121
[13]
B. Hammer, J. K. Nörskov, Nature, 1995, 376, 238-240. DOI:10.1038/376238a0
[14]
T. Bligaard, J. K. Nörskov, S. Dahl, J. Matthiesen, C. H. Christensen, J. Sehested, J. Catal., 2004, 224, 206-217. DOI:10.1016/j.jcat.2004.02.034
[15]
J. K. Nörskov, T. Bligaard, B. Hvolbæk, F. Abild-Pedersen, I. Chorkendorff, C. H. Christensen, Chem. Soc. Rev., 2008, 37, 2163-2171. DOI:10.1039/b800260f
[16]
M. P. Andersson, T. Bligaard, A. Kustov, K. E. Larsen, J. Greeley, T. Johannessen, C. H. Christensen, J. K. Nörskov, J. Catal., 2006, 239, 501-506. DOI:10.1016/j.jcat.2006.02.016
[17]
T. Bligaard, J. K. Nörskov, In: Chemical Bonding at Surfaces and Interfaces, A. Nillson, L. G. M. Petterson, J. K. Norskov eds. Amsterdam, 2008, 255–321.
[18]
X. Zhu, Q. Guo, Y. Sun, S. Chen, J. Q. Wang, M. Wu, W. Fu, Y. Tang, X. Duan, D. Chen, Y. Wan, Nat. Commun., 2019, 10, 1428. DOI:10.1038/s41467-019-09421-5
[19]
S. Díez-González, N. Marion, S. P. Nolan, Chem. Rev., 2009, 109, 3612-3676. DOI:10.1021/cr900074m
[20]
D. Ren, L. He, L. Yu, R. S. Ding, Y. M. Liu, Y. Cao, H. Y. He, K. N. Fan, J. Am. Chem. Soc., 2012, 134, 17592-17598. DOI:10.1021/ja3066978
[21]
V. Sridharan, P. A. Suryavanshi, J. C. Menéndez, Chem. Rev., 2011, 111, 7157-7259. DOI:10.1021/cr100307m
[22]
H. Arakawa, M. Aresta, J. N. Armor, M. A. Barteau, E. J. Beckman, A. T. Bell, J. E. Bercaw, C. Creutz, E. Dinjus, D. A. Dixon, K. Domen, D. L. DuBois, J. Eckert, E. Fujita, D. H. G. Gibson, A. William, D. W. Goodman, J. Keller, G. J. Kubas, H. H. Kung, J. E. Lyons, L. E. Manzer, T. J. Marks, K. Morokuma, K. M. Nicholas, R. Periana, L. Que, J. Rostrup-Nielson, W. M. H. Sachtler, L. D. Schmidt, A. Sen, G. A. Somorjai, P. C. Stair, B. R. Stults, W. Tumas, Chem. Rev., 2001, 101, 953-996. DOI:10.1021/cr000018s
[23]
J. R. Kitchin, J. K. Nörskov, M. A. Barteau, J. G. Chen, J. Chem. Phys., 2004, 120, 10240-10246. DOI:10.1063/1.1737365
[24]
P. Strasser, S. Koh, T. Anniyev, J. Greeley, K. More, C. Yu, Z. Liu, S. Kaya, D. Nordlund, H. Ogasawara, M. F. Toney, A. Nilsson, Nat. Chem., 2010, 2, 454-460. DOI:10.1038/nchem.623
[25]
M. Guo, C. Li, Q. Yang, Catal. Sci. Technol., 2017, 7, 2221-2227. DOI:10.1039/C7CY00394C
[26]
S. Venkatachalam, T. Jacob, Phys. Chem. Chem. Phys., 2009, 11, 3263-3270. DOI:10.1039/b900250b
[27]
K. Li, B. Hao, M. Xiao, Y. Kuang, H. Shang, J. Ma, Y. Liao, H. Mao, Appl. Surf. Sci., 2019, 478, 176-182. DOI:10.1016/j.apsusc.2019.01.225
[28]
J. M. Richardson, C. W. Jones, J. Catal., 2007, 251, 80-93. DOI:10.1016/j.jcat.2007.07.005
[29]
S. Chen, L. Meng, B. Chen, W. Chen, X. Duan, X. Huang, B. Zhang, H. Fu, Y. Wan, ACS Catal., 2017, 7, 2074-2087. DOI:10.1021/acscatal.6b02720
[30]
L. Duan, R. Fu, B. Zhang, W. Shi, S. Chen, Y. Wan, ACS Catal., 2016, 6, 1062-1074. DOI:10.1021/acscatal.5b02147
[31]
R. J. Madon, M. Boudart, Ind. Eng. Chem. Fun., 1982, 21, 438-447. DOI:10.1021/i100008a022
[32]
Y. Gong, P. Zhang, X. Xu, Y. Li, H. Li, Y. Wang, J. Catal., 2013, 297, 272-280. DOI:10.1016/j.jcat.2012.10.018
[33]
H. Wise, J. Catal., 1968, 10, 69-72. DOI:10.1016/0021-9517(68)90224-8
[34]
H. F. Cordes, J. Phys. Chem., 1968, 72, 2185-2189. DOI:10.1021/j100852a052
[35]
W. M. H. Sachtler, R. Jongepier, J. Catal., 1965, 4, 665-671. DOI:10.1016/0021-9517(65)90266-6
[36]
G. J. Hutchings, C. J. Kiely, Acc. Chem. Res., 2013, 46, 1759-1772. DOI:10.1021/ar300356m
[37]
S. Marx, A. Baiker, J. Phys. Chem. C., 2009, 113, 6191-6201. DOI:10.1021/jp808362m
[38]
A. Venezia, Appl. Catal. A, 2003, 251, 359-368.
[39]
A. M. Venezia, V. La Parola, V. Nicol, G. Deganello, J. Catal., 2002, 212, 56-62. DOI:10.1006/jcat.2002.3778
[40]
J. Pritchard, L. Kesavan, M. Piccinini, Q. He, R. Tiruvalam, N. Dimitratos, J. A. Lopez–Sanchez, A. F. Carley, J. K. Edwards, C. J. Kiely, G. J. Hutchings, Langmuir, 2010, 26, 16568-16577. DOI:10.1021/la101597q
[41]
L. Duan, R. Fu, Z. Xiao, Q. Zhao, J. Q. Wang, S. Chen, Y. Wan, ACS Catal., 2015, 5, 575-586. DOI:10.1021/cs501356h
[42]
J. C. Vedrine, M. Dufaux, C. Naccache, B. Imelik, J. Chem, . Soc, Faraday Trans. 1, 1978, 74, 440-449.
[43]
P. A. P. Nascente, S. G. C. de Castro, R. Landers, G. G. Kleiman, Phys. Rev. B, 1991, 43, 4659-4666. DOI:10.1103/PhysRevB.43.4659
[44]
N. Mårtensson, R. Nyholm, H. Cal, J. Hedman, B. Johansson, Phys. Rev. B, 1981, 24, 1725-1738. DOI:10.1103/PhysRevB.24.1725
[45]
M. B. Griffin, A. A. Rodriguez, M. M. Montermore, J. R. Monnier, C. T. Williams, J. W. Medlin, J. Catal., 2013, 307, 111-120. DOI:10.1016/j.jcat.2013.07.012
[46]
A. F. Lee, C. J. Baddeley, C. Hardacre, R. M. Ormerod, R. M. Lambert, G. Schmid, H. West, J. Phys. Chem., 1995, 99, 6096-6102. DOI:10.1021/j100016a053
[47]
M. R. Knecht, M. G. Weir, A. I. Frenkel, R. M. Crooks, Chem. Mater., 2008, 20, 1019-1028. DOI:10.1021/cm0717817
[48]
P. Dash, T. Bond, C. Fowler, W. Hou, N. Coombs, R. W. J. Scott, J. Phys. Chem. C, 2009, 113, 12719-12730. DOI:10.1021/jp9037182
[49]
A. Christensen, A. V. Ruban, P. Stoltze, K. W. Jacobsen, H. L. Skriver, J. K. Nörskov, F. Besenbacher, Phys. Rev. B, 1997, 56, 5822-5834. DOI:10.1103/PhysRevB.56.5822
[50]
R. Marchal, A. Genest, S. Krüger, N. Rösch, J. Phys. Chem. C, 2013, 117, 21810-21822. DOI:10.1021/jp4061686
[51]
S. Wang, J. Wang, Q. Zhao, D. Li, J. Q. Wang, M. Cho, H. Cho, O. Terasaki, S. Chen, Y. Wan, ACS Catal., 2015, 5, 797-802. DOI:10.1021/cs501896c
[52]
S. N. Reifsnyder, H. H. Lamb, J. Phys. Chem. B, 1999, 103, 321-329. DOI:10.1021/jp982893b
[53]
F. Liu, D. Wechsler, P. Zhang, Chem. Phys. Lett., 2008, 461, 254-259. DOI:10.1016/j.cplett.2008.07.029
[54]
M. Ø, ; . Pedersen, S. Helveg, A. Ruban, I. Stensgaard, E. Lægsgaard, J. K. Nörskov, F. Besenbacher, Surf. Sci., 1999, 426, 395-409. DOI:10.1016/S0039-6028(99)00385-4
[55]
R. J. Davis, S. M. Landry, J. A. Horsley, M. Boudart, Phys. Rev. B, 1989, 39, 10580-10583. DOI:10.1103/PhysRevB.39.10580
[56]
J. Goetz, M. A. Volpe, A. M. Sica, C. E. Gigola, R. Touroude, J. Catal., 1995, 153, 86-93. DOI:10.1006/jcat.1995.1110
[57]
H. Xin, A. Vojvodic, J. Voss, J. K. Nörskov, F. Abild-Pedersen, Phys. Rev. B, 2014, 89, 115114. DOI:10.1103/PhysRevB.89.115114
[58]
J. A. McCaulley, Phys. Rev. B., 1993, 47, 4873-4879. DOI:10.1103/PhysRevB.47.4873
[59]
W. Jones, R. Su, P. P. Wells, Y. Shen, N. Dimitratos, M. Bowker, D. Morgan, B. B. Iversen, A. Chutia, F. Besenbacher, G. Hutchings, Phys. Chem. Chem. Phys., 2014, 16, 26638-44. DOI:10.1039/C4CP04693E
[60]
J. E. Muller, O. Jepsen, O. K. Andersen, J. W. Wilkins, Phys. Rev. Lett., 1978, 40, 720-722. DOI:10.1103/PhysRevLett.40.720
[61]
W. Olovsson, C. Göransson, L. V. Pourovskii, B. Johansson, I. A. Abrikosov, Phys. Rev. B, 2005, 72, 064203. DOI:10.1103/PhysRevB.72.064203
[62]
E. G. Allison, G. C. Bond, Catal. Rev., 1972, 7, 233-289. DOI:10.1080/01614947208062259
[63]
E. J. Cho, S. Lee, S. J. Oh, M. Han, Y. S. Lee, C. N. Whang, Phys. Rev. B, 1995, 52, 16443-16450. DOI:10.1103/PhysRevB.52.16443
[64]
E. Choi, S. J. Oh, M. Choi, Phys. Rev. B, 1991, 43, 6360-6368. DOI:10.1103/PhysRevB.43.6360
[65]
N. J. Castellani, D. B. Leroy, Z. Phys. B, 1988, 71, 315-319. DOI:10.1007/BF01312491
[66]
W. Tu, Y. H. Chin, Angew. Chem. Int. Ed., 2014, 53, 12148-12152. DOI:10.1002/anie.201405232
[67]
W. Tu, M. Ghoussoub, C. V. Singh, Y.-H. C. Chin, J. Am. Chem. Soc., 2017, 139, 6928-6945. DOI:10.1021/jacs.7b01632
[68]
W. Chen, D. Li, C. Peng, G. Qian, X. Duan, D. Chen, X. Zhou, J. Catal., 2017, 356, 186-196. DOI:10.1016/j.jcat.2017.10.016
[69]
B. Hammer, J. K. Nörskov, Adv. Catal., 2000, 45, 71-129.
[70]
G. Santarossa, M. Lannuzzi, A. Vargas, A. Baiker, ChemPhysChem, 2008, 9, 401-413. DOI:10.1002/cphc.200700534
[71]
T. Rowland, N. E. Cusack, R. G. Ross, J. Phys. F:Met. Phys., 1974, 4, 2189-2202. DOI:10.1088/0305-4608/4/12/015
[72]
H. Falsig, B. Hvolbaek, I. S. Kristensen, T. Jiang, T. Bligaard, C. H. Christensen, J. K. Nörskov, Angew. Chem. Int. Ed., 2008, 47, 4835-4839. DOI:10.1002/anie.200801479
[73]
N. Lopez, J. K. Nörskov, J. Am. Chem. Soc., 2002, 124, 11262-11263. DOI:10.1021/ja026998a