{'fraction_non_alphanumeric': 0.04162827993312178, 'fraction_numerical': 0.03545228102501109, 'mean_word_length': 4.501783367749202, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Disentangling the size-dependent geometric and electronic effects of palladium nanocatalysts beyond selectivity

Site-selective blocking by ALD disentangles the geometric and electronic effects of catalysts beyond the selectivity volcano.


INTRODUCTION
Metal nanoparticles (NPs) are one of the most important materials in a broad application of heterogeneous catalysis, sensors, batteries, biology and medicine, etc., because of the notable size effects from their bulk counterparts (1,2). Geometrically, in addition to high surface area, the smaller NPs expose more low-coordination sites (LCSs) instead of high-coordination sites (HCSs), and their different local geometries considerably change chemical bond breaking and making in catalysis (3)(4)(5)(6)(7)(8)(9). Electronically, quantum size effects (QSEs) also alter significantly the energy levels of the entire particle, thereby affecting the orbital hybridization and overall charge transfer between metal and reactants (10,11). The above two distinct functionalities endow precise tailoring of the particle size, in addition to the particle composition (12,13), to be a valuable methodology to optimize catalysts (14)(15)(16)(17)(18). Nevertheless, HCSs and LCSs might have very different electronic structures (19), whereas QSEs depend sensitively on geometric structures (20). This strong entanglement of the geometric and electronic properties with the particle size makes identification of their distinct size-dependent functionalities nearly impossible. Disentangling the geometric and electronic properties of metal NPs and identifying when one or the other starts to dominate the reactivity would go beyond the established structure-activity relationship and open a new dimension for rational design of the catalysts (21,22). However, the above demand has never been met and represents a grand challenge.
Selective oxidation of alcohols to the corresponding carbonyl compounds, the key intermediates for fine chemical synthesis, is one of the most fundamental and useful industrial processes (23). Supported palladium (Pd)-based nanocatalysts have attracted significant attention because of their superior catalytic performance under mild conditions with molecular oxygen as the oxidant (23)(24)(25). However, the effects of the structural and electronic properties of Pd NPs on activity and selectivity remain controversial (17,23,25,26).
In this work, we report that the strong entanglement of the electronic and geometric structures of Pd NPs in Pd/Al 2 O 3 catalysts prompts the activity and selectivity both in volcano relationships with the Pd particle size unprecedentedly in aerobic oxidation of benzyl alcohol. To manipulate solely the geometric structure without interfering with the electronic structures, we developed a new strategy of selectively blocking the LCSs and HCSs of Pd NPs using atomic layer deposition (ALD) of Al 2 O 3 and FeO x , respectively. Combining with density functional theory (DFT) calculations, we disentangle unambiguously the geometric and electronic effects of Pd NPs on the activity and selectivity of oxidation of benzyl alcohol. We find that it is the geometric effect that controls the right side of the volcano peak at~4 nm, whereas the electronic effect controls the left side. Selectively blocking the Pd LCSs by Al 2 O 3 ALD goes beyond the selectivity volcano and results in a remarkable benzaldehyde selectivity and intrinsic activity on the 4-nm Pd NP catalyst at the same time.

RESULTS AND DISCUSSION
Solvent-free aerobic benzyl alcohol oxidation using O 2 as the oxidant for benzaldehyde was conducted at 120°C over various Pd/Al 2 O 3 catalysts, with the Pd particle size carefully tuned from~2.1 to 19.1 nm (figs. S1 and S2 and table S1). Typically, there are two major reaction pathways: the dehydrogenation path to benzaldehyde (the desired product) (Eq. 1) and the hydrogenolysis path to toluene (the major by-product at this reaction temperature) (Eq. 2) (27).
The H* in Eq. 2 is the intermediate surface hydrogen species generated during the dehydrogenation process. Strikingly, an inverse volcano relationship for the selectivity to benzaldehyde on the Pd particle size was found (Fig. 1A). On the left and right sides of this volcano, selectivity to benzaldehyde was as high as 87% (2.1 nm) and 93% (19.1 nm) at 50% conversion, while a lowest benzaldehyde selectivity of 72% occurred at the valley of the inverse volcano on the 4.2-nm Pd/Al 2 O 3 sample. Meanwhile, the activity in terms of turnover frequencies (TOFs) and specific rates also showed a traditional volcano shape with a maximum of 4.3 × 10 4 hour −1 and 1.2 × 10 4 hour −1 on the 4.2-nm Pd/Al 2 O 3 sample, respectively (Fig. 1B). At the volcano peak, corresponding apparent activation energies for formation of benzaldehyde and toluene were measured with values of 0.36 and 0.83 eV, respectively ( fig. S3). Note that the volcano relationship of TOFs on the Pd particle size is in good agreement with the literature (17).
In situ x-ray photoelectron spectroscopy (XPS) measurements showed that the Pd 3d binding energy gradually shifted to a higher position in a magnitude of 1.3 eV as Pd particle size decreased from 19.1 to 2.1 nm ( Fig. 2A). A sharp increase in Pd 3d binding energy was noticed below 4.2 nm (Fig. 2B), indicating a dramatic change in the corresponding electronic structure. To see the change in geometric structure, we estimated the HCS-to-LCS ratio by assuming all Pd NPs synthesized in a cubooctahedral shape, as often adopted in the literature (25). It is found that the resulting HCS-to-LCS ratio increases significantly from 0.93 to 16.7, with an increase in the Pd particle size from~2.1 to 19.1 nm (Fig. 2B). These results show that both the geometric (the HCS-to-LCS ratio) and electronic structures change largely with size. However, the entanglement between the geometric and electronic structures prevents differentiation of their individual contributions to the observed volcano relationships.
To disentangle the geometric and electronic effects, ALD of Al 2 O 3 and FeO x was performed on the 4.2-nm Pd/Al 2 O 3 sample to selectively block the LCSs and HCSs of Pd NPs, respectively (herein nAl 2 O 3 -Pd and nFeO x -Pd for n cycles of Al 2 O 3 and FeO x ALD, respectively) (28,29). Transmission electron microscopy (TEM) and energy-dispersive x-ray spectroscopy (EDS) elemental mapping confirmed the Al 2 O 3 and FeO x coating layers on Pd NPs (figs. S4 and S5). It was found that both the average Pd particle size ( fig. S6) and the corresponding electronic properties, measured by in situ XPS ( fig. S7), were essentially intact after overcoating. Diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) of CO chemisorption showed that the uncoated sample exhibited three peaks centered at about 1922, 1966, and 2088 cm −1 (Fig. 3A). They can be assigned to the m 2 bridgebonded CO on Pd(111) facets (HCSs), m 2 bridge-bonded CO on the Pd NPs edges/corners (LCSs), and linear CO on HCSs, respectively (28,30,31). When ALD of Al 2 O 3 was applied, the peak at 1966 cm −1 decreased much more quickly compared with the peak at 1922 cm −1 with increase in ALD cycles and vanished on 10Al 2 O 3 -Pd. The above evidence suggests firmly that the LCSs of Pd NPs in 10Al 2 O 3 -Pd were selectively blocked by the Al 2 O 3 overcoat, leaving the HCSs accessible for catalysis (inset of Fig. 3A), according to the literature (32,33). In contrast, FeO x coating showed an opposite trend by attenuating the CO peak at 1922 cm −1 much more quickly with ALD cycles (Fig. 3B), suggesting a preferential decoration of the HCSs. Semiquantitative analysis of CO DRIFT chemisorption showed that the HCS-to-LCS ratio was 4.3 for uncoated Pd/Al 2 O 3 , close to estimation of 3.0 by assuming corresponding Pd NPs in a cubooctahedral shape ( Fig. 2B and fig. S8) (25). Although this semiquantification has certain error bars, it could still properly represent the evolution of the HCS-to-LCS ratio with oxide overcoating. We found that the HCS-to-LCS ratio increased rapidly with Al 2 O 3 overcoating and became about 55 on 10Al 2 O 3 -Pd, with negligible LCSs exposed. On the contrary, the FeO x overcoating decreased this ratio considerably to~0.8 for 8FeO x -Pd.
The above results show that selective blocking of Pd NPs can tune effectively the HCS-to-LCS ratio from~0.8 to~55 without changing the particle size necessarily ( fig. S9), thus providing an ideal platform for disentangling the geometric and electronic effects. With the increase of Al 2 O 3 ALD cycles, benzaldehyde selectivity increased gradually from only 72 to almost 90% (fig. S10A), whereas the toluene selectivity decreased from about 24 to only 0.1% (fig. S10B). In contrast, with the increase of the FeO x coverage, benzaldehyde selectivity decreased from 72 to 62%, along with an increase in toluene selectivity from 22 to 37%. According to the strong support effect observed by Hutchings and co-workers (34), the Pd-FeO x interface might also influence the selectivity. However, we found that a Pd/Fe 2 O 3 catalyst with a Pd particle size of 3.5 ± 0.4 nm exhibited a benzaldehyde selectivity of 90% at 50% conversion, even higher than the Pd/Al 2 O 3 catalysts with similar particle sizes (Fig. 1A). Therefore, it turns out that the Pd-FeO x interface favors the benzaldehyde formation (35) and that the decreased benzaldehyde selectivity through FeO x overcoating is mainly due to the geometric effect through "site-selective blocking" of the HCSs of Pd NPs, rather than the formation of the Pd-FeO x interface. Plotting the selectivity of the two products together showed compellingly that with increase of the HCS-to-LCS ratio, benzaldehyde selectivity increased but the toluene selectivity decreased monotonously ( Fig. 3C). At the HCS-to-LCS ratio of~50, the corresponding TOF of 4.3 × 10 4 hour −1 remained similar with the volcano peak of uncoated 4.2-nm Pd/Al 2 O 3 ( Fig. 1B and fig. S11). Notably, toluene selectivity was nearly completely suppressed. This remarkable selectivity manipulation by oxide overcoating was also observed at different reaction temperatures and other Pd catalysts with different particle sizes and supports (figs. S12 to S17). As shown in Fig. 1, site-selective blockages of the LCSs broke down the selectivity volcano on the particle size and achieved a great performance in both activity and selectivity on 4.2-nm Pd/Al 2 O 3 . It is also worth noting that the 10Al-Pd/Al 2 O 3 catalyst demonstrated excellent recyclability and metal stability against sintering and leaching (fig. S18).
Selective blocking of Pd NPs by Al 2 O 3 and FeO x provides unambiguous evidence that dehydrogenation to benzaldehyde in benzyl alcohol oxidation takes place mainly at the HCSs, whereas hydrogenolysis to toluene occurs at the LCSs. Therefore, the geometric effect, namely, the decreased HCS-to-LCS ratio, controls the observed decrease in benzaldehyde selectivity as a decrease in the Pd particle size from~19.1 to 4.2 nm (Fig. 1A). However, when the particle size was smaller than 4.2 nm, the trend of the benzaldehyde selectivity was reversed. The resulting inverse volcano implies that the electronic effect starts to control the selectivity, as supported by a significant increase in the Pd 3d energy level (Fig. 2).
DFT calculations were performed to shed more light on the geometric and electronic effects of Pd NPs in benzyl alcohol oxidation. First, terraced Pd(111) and stepped Pd(211) were used to simulate the HCSs and LCSs of Pd NPs, respectively (figs. S19 to S21). For aerobic oxidation of benzyl alcohol to benzaldehyde and water, oxygen is mainly for promoting dehydrogenation and removing surface hydrogen. Oxygen-assisted dehydrogenation of benzyl alcohol to benzyloxy and hydroxyl (Fig. 4A, figs. S22 to S24, and table S2) is facile, with  barriers of 0.38 eV on Pd(111) and 0.44 eV on Pd(211). Further dehydrogenation of benzyloxy to benzaldehyde and hydrogen without assistance by oxygen necessarily has modest barriers of 0.51 eV on Pd(111) and 0.44 eV on Pd(211). In the following, reaction between hydroxyl and hydrogen to generate water closes the catalytic cycle. Calculated barriers are 0.75 eV on Pd(111) and 1.12 eV on Pd(211) (Fig. 4B); therein, the lower barrier for the former comes from its weaker binding toward hydroxyl, compared with the latter (−2.68 eV versus −3.04 eV) ( fig. S20), which is in agreement with previous result (36). In this cycle, benzaldehyde formation is facile at both the HCSs and LCSs, and hydrogen removal by hydroxyl is demanding and more efficient at the HCSs.
For toluene formation, the scission of the C-O bond in benzyl alcohol or its derived intermediates is the prerequisite. Compared with the scission of the C-O bond in benzyl alcohol-derived intermediates, we found that direct scission of the C-O bond of benzyl alcohol to form benzyl and hydroxyl is preferential with barriers of 1.02 eV on Pd(111) and 0.95 eV on Pd(211) (figs. S25 and S26 and table S2). Hydrogen's reaction with benzyl produces toluene with barriers of 1.16 eV on Pd(111) and 0.98 eV on Pd(211) (Fig. 4B). Compared with hydrogen's reaction with hydroxyl, hydrogen's reaction with benzyl is kinetically more favorable on Pd(211), but not on Pd(111). Namely, LCSs prefers to hydrogenate benzyl for toluene, and HCSs prefers to hydrogenate hydroxyl for water. It is therefore the competitive removal of hydrogen governing the reaction paths on HCSs and LCSs.
With an increase in the HCS-to-LCS ratio, the benzaldehyde selectivity would increase while the toluene decreases, as found exactly on the right side of the volcano relationship (Fig. 1A) and the site-selective blocking experiments (Fig. 3C).
To see the electronic effect on the small particles and rationalize the left side of the volcano relationship, we investigated the adsorption of benzyl alcohol, benzyl, and hydroxyl on freestanding Pd n clusters (n = 55, 147, and 309; figs. S27 and S28) with cubooctahedral structures (11,37). Although the change of hydroxyl binding is negligible, the binding strength of benzyl alcohol and benzyl at favorable sites increases pronouncedly with a decrease in the size (Fig. 4C), 0.40 eV stronger on Pd 55 than on Pd(211), where the calculated binding energies of benzyl alcohol and benzyl are −2.05 and −2.51 eV, respectively ( fig. S19). Changes in binding energies cannot come from the geometric effect, since the corresponding adsorption sites ( fig. S28) are similar. On the other hand, the clusters constructed have different electronic structures from Pd(211); in particular, the work function decreases systematically by~0.3 eV for Pd 55 . This indicates that the corresponding Fermi level shifts upward with a decrease in particle size, in line with upshifting the Pd 3d level to higher positions observed experimentally (Fig. 2). Accordingly, there is more electron transfer toward the adsorbates (Fig. 4C), enhancing eventually the corresponding binding strengths, a fact that was also found previously in Pt and Au (11,38). As shown in Fig. 4D and figs. S29 and S30, increasing electron transfer toward the adsorbates with a decrease in the particle size comes not only from the coordinated metal atoms contributing more electrons but also from the fact that there are more metal atoms contributing electrons. This is understandable since a change in Fermi level affects the entire particle, in contrast to the HCSs and LCSs differing only in their local geometries of the particle.
Enhanced binding strengths on small particles change significantly the activity and selectivity. Taking Pd 55 as an example, the barriers for the O-H (0.56 eV) bond and subsequent C-H (0.82 eV) bond cleavages of benzyl alcohol are higher than those on Pd(111) and Pd(211), respectively ( Fig. 4A and figs. S31 and S32). Moreover, a barrier of 1.03 eV for hydrogen removal by hydroxyl is lower than that of 1.15 eV by benzyl (Fig. 4B), opposite to the order on Pd(211). This indicates that although the toluene formation is more favorable than water at the LCSs, stronger binding of benzyl due to the pronounced electronic effect of small particles leads to subsequent hydrogenation less favorable than water formation. These rationalize well the low activity and the high selectivity to benzaldehyde on small particles on the left side of volcanos (Fig. 1). In the moderate particle size, the electronic effect is weakened, the geometric effect on the HCS-to-LCS ratio increases, and the interplay between these two effects results in the volcano of activity and selectivity on the particle size.
In summary, we have demonstrated that in aerobic oxidation of benzyl alcohol, both the activity and selectivity over Pd/Al 2 O 3 catalysts showed volcano relationships with the Pd particle size unprecedentedly, due to the strong entanglement of electronic and geometric structures. Selective blocking of Pd NPs by Al 2 O 3 and FeO x , according to DRIFTS CO chemisorption, provides a powerful tool to manipulate solely the geometric structure, which unambiguously confirms that the geometric effect controls the right side of the volcanos. On the smaller Pd particles, the electronic effect is involved significantly as indicated by XPS, thus reversing the trend of the benzaldehyde selectivity and decreasing the activity. First principles calculations further unveiled that the competitive removal of surface hydrogen via reaction with hydroxyl and benzyl intermediates on HCSs and LCSs makes the reaction structure sensitive and selective. Decreasing the particle size lowers the work function and increases the binding strength of benzyl alcohol and benzyl and the reaction barriers. These effects, in turn, alter the reaction path of hydrogen removal and significantly change the corresponding activity and selectivity, consistent with the experimental observation. Given the essential roles of the geometric and electronic effects in metal NPs, disentangling the geometric and electronic effects of metal NPs opens a new dimension for metal catalyst rational design.

Catalyst preparation
Pd/Al 2 O 3 catalysts were synthesized using the incipient wetness impregnation (WI) method (28). For the synthesis of the 4.2-nm Pd/Al 2 O 3 catalyst, 0.1 g of Pd(NO 3 ) 2 ·2H 2 O (>97.7%; Aladdin Chemicals) and 0.43 g of citric acid were first dissolved in 2.6 ml of deionized water. Then, 4.4 g of spherical Al 2 O 3 powder (99.5%; NANODUR, Alfa Aesar) was added into the solution and mixed thoroughly. The mixture was dried overnight at 125°C. Last, the dried material was calcined in 10% O 2 in Ar (Nanjing Special Gases) at 300°C for 5 hours and then reduced at 250°C for another 1 hour in 10% H 2 in Ar (Nanjing Special Gases) to obtain the 4.2-nm Pd/Al 2 O 3 catalyst. Tuning the ratio of the amount of the Pd precursor to citric acid, as well as the calcination tem-perature, would allow precise control of the Pd particle size. Details can be seen in table S1.
A Pd/C catalyst was also prepared for excluding the support effect according to a procedure reported previously (39). Briefly, 0.1 mmol PdCl 2 (Aladdin Chemicals) and 0.8 mmol sodium citrate (Sinopharm Chemical Reagent Co. Ltd.) were dissolved into 150 ml of water. Carbon black (400 mg; Vulcan XC72R, Carbot Corp.) was then added. After stirring the mixture for 20 min, followed by a 30-min sonication, 15 ml of 0.1 M sodium borohydride (96%; Sinopharm Chemical Reagent Co. Ltd.) solution was added dropwise into the suspension under vigorous stirring and then aged at 25°C for 8 hours. Next, the precipitate was filtered and washed with deionized water thoroughly. The obtained materials were then dried overnight in a vacuum oven at 45°C. Last, the dried sample was calcined in 10% O 2 in Ar at 250°C for 1 hour and then reduced at 250°C for 1 hour in 10% H 2 in Ar to obtain the Pd/C catalyst.
Alumina and FeO x ALD were carried out in a viscous flow reactor (GEMStar-6 Benchtop ALD, Arradiance). Ultrahigh purity N 2 (99.999%) was used as a carrier gas at a flow rate of 200 ml/min. Al 2 O 3 ALD was performed by alternatively exposing the sample to trimethylaluminum (TMA; 99.999%; Nanjing MO Yuan Scientific Instruments & Materials Co., Ltd.) and deionized water at 200°C at different cycles (28). The timing sequence was 3, 200, 10, and 200 s for TMA exposure, N 2 purge, water exposure, and N 2 purge, respectively. Besides the 4.2-nm Pd/Al 2 O 3 samples, alumina overcoating was also applied on 5.1-nm Pd/Al 2 O 3 and Pd/C samples. The obtained alumina-coated 5.1-nm Pd/Al 2 O 3 and Pd/C samples were denoted as nAl 2 O 3 -Pd-5.1 and nAl 2 O 3 -Pd/C, respectively, where n represents the number of ALD cycles (n = 1, 3, 5, and 10). All catalysts were calcined in 10% O 2 in Ar at 300°C for 2 hours and then reduced at 200°C for 1 hour in 10% H 2 in Ar before any characterization and reaction tests. FeO x ALD was performed by alternatively exposing the sample to ferrocene (98%; Sigma-Aldrich) and ultrahigh purity O 2 (99.999%) at 150°C at different cycles (29). The timing sequence was 60, 200, 60, and 200 s for ferrocene exposure, N 2 purge, O 2 exposure, and N 2 purge, respectively. Besides the 4.2-nm Pd/Al 2 O 3 samples, FeO x overcoating was also applied on 5.1-nm Pd/Al 2 O 3 and Pd/C samples. The obtained FeO x -coated 5.1-nm Pd/Al 2 O 3 and Pd/C samples were denoted as nFeO x -Pd-5.1 and nFeO x -Pd/C, respectively, where n represents the number of ALD cycles (n = 1, 3, 5, and 10). All catalysts were calcined in 10% O 2 in Ar at 150°C for 1 hour and then reduced at 150°C for 1 hour in 10% H 2 in Ar before any characterization and reaction tests.
Catalyst characterization TEM measurements were performed on a JEM-2100F instrument operated at 200 kV to characterize the morphology of the 10Al 2 O 3 -Pd catalyst. Scanning TEM (STEM) characterization of 8FeO x -Pd catalysts and Pd/Al 2 O 3 samples with various particle sizes was carried out on an aberration-corrected high-angle annular dark-field STEM instrument at 200 kV (JEM-ARM200F, University of Science and Technology of China). Elemental mapping using EDS was performed on the same equipment. The Pd contents in various catalysts were analyzed by an inductively coupled plasma atomic emission spectrometer.
The DRIFTS CO chemisorption measurements were performed on a Nicolet iS10 spectrometer equipped with a mercury-cadmiumtelluride detector and a low-temperature reaction cell (Praying Mantis Harrick). After loading a sample into the cell, it was first calcined in 10% O 2 in Ar at 150°C for 1 hour followed by a reduction in 10% H 2 in Ar at 150°C for another 1 hour. After cooling the sample to room temperature under a continuous flow of Ar, a background spectrum was collected. Subsequently, the sample was exposed to 10% CO in Ar at a flow rate of 20 ml/min for about 30 min until saturation. Next, the sample was purged with Ar at a flow rate of 20 ml/min for another 30 min to remove the gas-phase CO, and then, the DRIFT spectrum was collected with 256 scans at a resolution of 4 cm −1 .
In situ x-ray photoemission spectroscopy (XPS) measurements were conducted at the BL10B beamline photoemission end station at the National Synchrotron Radiation Laboratory (NSRL) in Hefei, China. Briefly, the beamline is connected to a bending magnet and covers photon energies from 100 to 1000 eV. The end-station consists of four chambers, i.e., analysis chamber, preparation chamber, quick sample load lock chamber, and a high pressure reactor. The analysis chamber, with a base pressure of <5 × 10 −10 torr, is connected to the beamline and equipped with a VG Scienta R3000 electron energy analyzer and a twin anode x-ray source. The high pressure reactor houses a reaction cell where the samples can be treated with different gases up to 20 bars and simultaneously heated up to 650°C. After sample treatment, the reactor can be pumped down to high vacuum (<10 −8 torr) for sample transfer. In the current work, the samples were first treated with the flowing 10% H 2 in Ar (20 ml/min) at 150°C for 0.5 hour at ambient pressure in the high pressure reactor. Next, the samples were transferred to the analysis chamber for XPS measurements in the Pd 3d region without exposing to air. The Al 2p binding energy at 74.3 eV was used as internal reference (40).
The Pd dispersions were determined by CO pulse chemisorption, which were conducted using a Micromeritics AutoChem II chemisorption analyzer at room temperature. Before CO pulse, the catalysts were first calcined in 10% O 2 in He at 200°C for 1 hour followed by a reduction in 10% H 2 in He at 200°C for 2 hours. Then, the catalysts were cooled to room temperature in He, and CO pulses were introduced to the catalyst surface using 10% CO in He (41). A stoichiometry of CO:Pd = 1 for dispersion calculations was assumed according to literature.

Catalyst evaluation
Solvent-free aerobic oxidation of benzyl alcohol using molecular O 2 was conducted in a batch-type reactor, as described in our previous work (42). The amount of catalyst used for the reaction test was normalized to have the same Pd content of 1.6 mmol. Typically, 5 ml of benzyl alcohol (48 mmol) and the normalized Pd catalyst were coadded into a 25-ml three-necked glass flask equipped with a reflux condenser. Before reaction, the system was first charged with O 2 by bubbling ultrahigh purity O 2 (99.999%; Nanjing Special Gases) at a flow rate of 15 ml/min for 20 min to remove air. Under the continuous flow of O 2 , the reactor was quickly immersed into a silicon oil bath at the desired temperature (90, 120, or 140°C) to initiate the reaction. During each reaction, the mixture was vigorously stirred at a rate of 1500 rpm to exclude any mass transfer limitation. The effect of oxygen flow rates on the possible selectivity change was also ruled out. Last, the reaction products were analyzed using a Shimadzu GC-2014 gas chromatograph, equipped with an Rtx-1 capillary column and an autoinjector. The TOFs were calculated using Eq. 3 TOF ¼ moles of benzly alcohol consumed moles of total active sites Â reaction time ð3Þ Here, the moles of active sites were normalized to the number of exposed surface Pd sites calculated by the cubooctahedral "magic cluster" model. On alumina-coated samples, the number of exposed surface Pd sites was determined by the Pd dispersion data (table S3). The TOFs were calculated from the rate of the intitial 60-min reaction, unless otherwise noted. The specific rates were also calculated on the basis of Eq. 4 Specific rate ¼ moles of benzly alcohol consumed moles of total Pd atoms Â reaction time ð4Þ Theoretical calculations DFT calculations were performed with the Vienna ab initio Simulation Package (43,44). The exchange-correlation interaction was described by optPBE functional with van der Waals correction (45)(46)(47). The Kohn-Sham equations were solved in a plane wave basis set with a kinetic energy cutoff of 400 eV. To model Pd(111) and Pd(211) surfaces, the four-layer slab models were used with (4 by 4) and (2 by 4) unit cells, respectively. The freestanding Pd 55 , Pd 147 , and Pd 309 clusters were modeled using the cubeoctahedral structure based on the literature (11,38). A (3 by 3 by 1) k-point mesh and G-point were used to sample the surface Brillouin zone for the surface and cluster, respectively. A 12-Å vacuum was introduced with correction of dipole moment between the repeated slabs/clusters. During optimization, the bottom two layers of the slab were fixed, while the remaining atoms of the slab, all the atoms of the clusters and the adsorbates were relaxed until the residual force was less than 0.02 eV/Å. The barriers of the elementary steps were calculated by the climbing image nudged elastic band (CI-NEB) method (48), and the transition states were confirmed by single imaginary frequency.

SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/ content/full/5/1/eaat6413/DC1 Supplementary Text Fig. S1. Morphologies of Pd/Al 2 O 3 catalysts with various particle sizes.   Fig. S24. Energy profiles of the direct and oxidative dehydrogenation paths on Pd surfaces. Fig. S25. Energy profiles of the hydrogenolysis path to toluene formation on Pd surfaces. Fig. S26. Energy profiles of C-O scission of (C 6 H 5 )CH 2 OH*, (C 6 H 5 )CHOH*, and (C 6 H 5 )CH 2 O* on Pd surfaces. Fig. S27. Models of freestanding Pd 55 , Pd 147 , and Pd 309 clusters. Fig. S28. The optimized most stable geometries of (C 6 H 5 )CH 2 OH*, (C 6 H 5 )CH 2 *, and OH* on freestanding Pd 55 , Pd 147 , and Pd 309 clusters. Fig. S29. Comparison of the side view of charge density difference upon adsorption of benzyl on Pd 55 , Pd 147 , Pd 309 , and Pd(211) surfaces. Fig. S30. Comparison of the top view of charge density difference upon adsorption of benzyl on Pd 55 , Pd 147 , Pd 309 , and Pd(211) surfaces. Fig. S31. Optimized transition state structures of the key steps involved in benzaldehyde and toluene formations on Pd 55 cluster. Fig. S32. Energy profiles of the oxidative dehydrogenation and hydrogenolysis paths on Pd 55 cluster. Table S1. Conditions for synthesis of Pd/Al 2 O 3 samples with various Pd particle sizes using the WI method.