Design of High-Performance Pd-based Alloy ... - ACS Publications

the efficient catalysts cannot be effectively achieved solely based on the well-known Sabatier analysis. In this paper, we proposed a strategy to desi...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Newcastle, Australia

Article

Design of High-Performance Pd-based Alloy Nanocatalysts for Direct Synthesis of H2O2 Haoxiang Xu, Daojian Cheng, and Yi Gao ACS Catal., Just Accepted Manuscript • Publication Date (Web): 09 Feb 2017 Downloaded from http://pubs.acs.org on February 9, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Catalysis is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Design of High-Performance Pd-Based Alloy Nanocatalysts for Direct Synthesis of H2O2 Haoxiang Xua, Daojian Chenga*, Yi Gaob,* a

International Research Center for Soft Matter, Beijing Key Laboratory of Energy Environmental Catalysis, State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, Beijing100029, China E-mail: [email protected] b

Division of Interfacial Water and Key Laboratory of Interfacial Physics and Technology,

Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China E-mail: [email protected]

ACS Paragon Plus Environment

1

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 21

ABSTRACT The direct synthesis of hydrogen peroxide (H2O2) is a promising alternative to the commercialized indirect process. However, it is still a big challenge for the development of Pdbased catalysts with outstanding activity and selectivity, because the design and optimization of the efficient catalysts cannot be effectively achieved solely based on the well-known Sabatier analysis. In this paper, we proposed a strategy to design more efficient Pd-based nanocatalysts combining the density functional theory (DFT) calculations and Sabatier analysis. The average valence electron of Pd-shell atoms is identified as the intrinsic factor for the activity and selectivity of the Pd-based nanocatalysts, which can be effectively tuned by the dopants. By introducing the dopants with suitable electronegativity, the valence electrons of Pd-shell atoms could be adjusted to the optimal range to enhance the activity and selectivity of the nanocluster simultaneously. With this strategy, Pd-W, Pd-Pb, Au-Pd-W, Au-Pd-Pb, Au-Pd-Mo and Au-PdRu are predicted as the potential candidates with the catalytic performance far exceeding the state-of-the-art experimental systems by scanning the periodic table. This work not only predicts potential Pd-based alloy nanocatalysts for direct synthesis of H2O2 for future experiments, but provides a viable way for the design of highly efficient heterogeneous catalysts in extensive applications.

KEYWORDS : direct synthesis of H2O2 , Pd-based alloys , nanoparticle , selectivity , catalyst design , density functional theory , electronegativity

ACS Paragon Plus Environment

2

Page 3 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

1. Introduction Hydrogen peroxide (H2O2) is well known to the general public due to its ubiquitous use as a bleaching agent in the pulp industry, as a disinfectant in the pharmaceutical industry and as a green oxidant in the synthesis of specialty chemicals.1-10 At present, hydrogen peroxide is manufactured industrially by an indirect process in which anthraquinones are sequentially hydrogenated and oxidized in a manner.11 However, considering the high capital and operating costs of energy-intensive separation and concentration in the indirect process,12-14 industry and academia are keen to identify an alternative direct process. To date, the vast majority of researches into the direct synthesis of H2O2 are based on supported Pd nanoparticles.15-23 However, it has not been demonstrated on commercial scale yet, since the addition of vast strong acid24 and halide15,25-26 promoters to the reaction medium promotes metal leaching and requires further purification of the H2O2 products despite suppressing the process of side reaction. Hutchings and co-workers firstly reported that Au-Pd27-32 alloy catalysts can significantly enhance the activity and selectivity in the direct synthesis of H2O2, which stimulates the extensive studies of the promotional effect by the addition of Au33-35, Pt36-38, Ag39, Ni40 to Pd catalysts both in the experiment and theory. Recently, Pd-Sn alloy catalyst is synthesized in experiment and its selectivity is characterized as the current state-of-the-art.41 Edwards et al. reported that the addition of a small amount of Pt to Au-Pd alloy catalyst significantly enhanced the activity and selectivity in the direct synthesis of H2O242. Despite of the above significant progress, the available experimental Pd-based catalysts cannot achieve the desired selectivity and reactivity simultaneously with an absent of acid and halides promoters. On the other hand, a wide range of theoretical studies have been carried out to search for new

ACS Paragon Plus Environment

3

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 21

catalysts for the direct synthesis of H2O2 using the Sabatier analysis.43,44 However, Sabatier analysis alone cannot provide guidelines for finding new nanoalloy catalysts, since Sabatier analysis merely provide the correlation between adsorption energy and activation energy. In particular, current studies are mostly focused on the dopants of noble metals. Thus, how to design more efficient Pd-based based on earth-abundant, cost-effective materials to replace noble metals is an urgent and practical task. In this work, we proposed a new strategy to design Pd-based binary and ternary alloy nanocatalysts for direct synthesis of H2O2. Combining with density functional theory (DFT) calculations, Sabatier analysis, and analysis of electronic structures, we identified the average valance electron of Pd-shell atoms as the intrinsic effect for the activity and selectivity of the nanocatalysts, which can be tuned by the doped metals. The doped elements with slightly higher electronegativity (compared to Pd) could withdraw electrons from Pd-shell atoms to the optimal range to enhance the activity and selectivity simultaneously. By the systematic search of the periodic table, Pd-W, Pd-Pb, Au-Pd-W, Au-Pd-Pb, Au-Pd-Mo and Au-Pd-Ru are predicted as good candidates for Pd-based binary and ternary alloy nanocatalysts, which exhibit better catalytic performance potentially than the available experimental systems. Our work serves as a general example to predict new catalysts, which may also find applications in designing new heterogeneous catalysts beyond direct synthesis of H2O2.

2. Computational Details DFT calculations are performed using the PWSCF (Plane-Wave Self-Consistent Field) planewave code in the Quantum ESPRESSO package.45 All the calculations were carried out using Perdew–Burke–Ernzerhof (PBE) xc-functional46 and ultrasoft pseudopotentials.47 The energy cut-offs for the selection of the plane waves for the description of the wave function and the

ACS Paragon Plus Environment

4

Page 5 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

electronic density were set to 40 and 400 Ry, respectively. All calculations on cluster models were performed in a 15×15×15 Å3 cubic cell. The first Brillouin zone was sampled at the Gamma-point and the electronic levels were broadened though a Gaussian smearing technique with a smearing parameter of 0.002 Ry. All calculations were performed fully spin-polarized, since it is well-known that magnetic considerations can be significant for certain oxygen-based adsorbates and transitionstates.48-49 Geometric optimizations was performed with a force convergence criterion of PdPt > Pd > AuPd > SnPd > Pt > Au in both the experimental results41,42 and our DFT calculations, indicating the DFT calculations can qualitatively characterize the relative ordering of H2O2 productivity rate by experiments. Meanwhile, the comparison between H2O2 decomposition/hydrogenation rate in experiments42 and side reaction activity by DFT calculations is listed in Table S6 and plotted in Figure 3b, which confirms the accuracy of activity characterization by DFT calculations. High activity of side

reaction

on

Sn25Pd30 in

DFT

calculations

is

in

conflict

with

slow

H2O2

hydrogenation/decomposition rate in experiment because its high selectivity is attributed to tin oxide surface layer that encapsulates small Pd-rich particles (active for decomposing and hydrogenating the H2O2) rather than Sn-Pd alloy.41

ACS Paragon Plus Environment

8

Page 9 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Figure 3: (a) Hydrogen peroxide productivity in experiment and calculated As for main reaction in theory as a function of different model clusters. (b) Hydrogen peroxide hydrogenation/decomposition in experiment and calculated As for side reaction in theoryas a function of different model clusters.

Sabatier analysis provides upper bounds for reaction activity and has successfully described trends in reactivity in terms of simple descriptors, e.g. adsorption energies of certain intermediates. In this paper, Sabatier analysis are performed through a multidimensional activity volcano with ∆GO* and ∆GH*. The correlations between As of main/side reaction and ∆GO*/∆GH* are derived from scaling relations of the binding energies of various adsorbates, and from Brønsted-Evans-Polanyi (BEP) relationships59-61 (see SI including Figure S1 and S2 for full details on the correlations). The activity volcano for main reaction is shown in Figure 4. Similar

ACS Paragon Plus Environment

9

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 21

to recent experimental work,42 Pt2Au23Pd30 locates at the top of volcano and shows very high activities in main reaction. In comparison, other clusters lie on hillside of the volcano and show lower activities towards H2O2 than Pt2Au23Pd30. Au25Pd30 shows slightly higher activity than Sn25Pd30, which can be verified by recent experimental work.41 As shown in Figure 4, Sabatier activity volcano is divided into two regions (by blue line) in terms of different predicted ratecontrolling steps. Where hydrogen binds very weakly, H2 dissociation is the rate-limiting step for a range of oxygen adsorption free energies, and the Au55 cluster falls into this region. Accordingly, OOH* formation from H* + O2* becomes the relevant rate-limiting step as the hydrogen binding free energy becomes stronger, and the other five clusters lie in this region.

Figure 4: Activity volcano for direct synthesis of H2O2. ∆GO*and ∆GH* are the adsorption free energies of oxygen and hydrogen adsorption, respectively. The activity, in kJ/mol, is plotted as a function of both ∆GO*and ∆GH*. In the above regions, the predicted rate-controlling steps in each region of descriptor space are shown. The dash area is the selective hydrogen peroxide region of the Sabatier volcano toward H2O2, which satisfy all of the selectivity reaction barrier inequality constraints in the text. The values for different model systems can be taken from their indicated positions.

ACS Paragon Plus Environment

10

Page 11 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

To determine whether the catalysts would selectively form hydrogen peroxide rather than water thermodynamically, the intersecting space of inequality constraints (see SI for full details) is delineated by shadow in Figure 4. The dashed area is based on the series of reaction barrier inequalities in Section 5 of Supporting Information, which indicates all activation barriers of elementary steps in main reaction path are lower than the corresponding competing ones in side reaction path. It is clear that Au55, Au25Pd30 and Pt2Au23Pd30 clusters conform to all activation barrier inequality constraints, which are favor for the production H2O2. In addition, Pt25Pd30, Pt55 and Pd55 clusters lie essentially out of the selectivity boundaries. These results are additionally in very good agreement with previous experiments.42,62-63 It is well-known that pure Au catalyst,63 50%Au-50%Pd catalysts62 and Pt-promoted Au-Pd catalysts42 are highly selective to hydrogen peroxide, while pure Pt/Pd catalyst and 50%Pt-50%Pd catalyst can form mixed product phases of H2O/H2O2.

3.2 Intrinsic effect of doped metal In order to get the deep understanding of intrinsic effect of the doped metal, we explored the electronic effects of the nanoclusters. Since Kohn–Sham orbital energy level diagram and d-band center is not a good indicator to explain the adsorption and activity trend of model systems (Figure S3-4), we turn to focus on charge redistribution caused by dopant. The activity of main/side reaction as a function of the average valance electron of Pd atoms in catalyst shell θ(Pd in shell) is plotted in Figure 5a. It is found that Pt2Au23Pd30 is of the highest activity for main reaction and outstanding selectivity. This could be attributed to the optimal θ for Pt2Au23Pd30 compared to Pd55, Pt25Pd30, Sn25Pd30 and Au25Pd30. DOS analysis of adsorbed O2 on Pd55, Au25Pd30, Pt25Pd30 and Pt2Au23Pd30 in Figure S5 shows that moderate electron from Pdbased alloy cluster to O-O bond can facilitate hydrogenation and suppress dissociation (detailed

ACS Paragon Plus Environment

11

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 21

explanation in SI). As been reported that there may be particle size effects on Pd particles for direct synthesis of hydrogen peroxide64, additional calculation on larger models, 79-atom cluster, are performed to explore whether the correlation in Figure 5a is robust to the effect of particle size (Figure S6), which shows the same order of the activity and selectivity for studied models as that in Figure 5a. As seen in Figure 5a, θ(Pd in shell) of Pd55 (10.12 |e|) is higher than the optimal range (9.85 |e| – 9.95 |e|), which suggests the performance of Pd55 can be promoted by the dopants which withdraw some electrons from Pd atoms. As is well-known, electronegativity quantifies the electron affinity of certain element in a chemical compound. The atom with higher electronegavitiy could withdraw electrons from the atom with lower electronegativity. For Pdbased alloys, lower electronegativity of Sn (1.96) increases θ(Pd in shell) of Sn25Pd30 (10.24 |e|) compared to Pd55 (10.12 |e|) since a few electrons of Sn atoms flow to Pd atoms. Pt (2.28) is of a slightly higher electronegativity than Pd (2.20), so θ(Pd in shell) of Pt25Pd30 is slightly reduced (10.07 |e|) compared to Pd55 (10.12 |e|) since a few electrons of Pd atoms flow to Pt atoms. Similarly, much higher electronegativity of Au (2.54) decreases more θ(Pd in shell) of Au25Pd30 (9.78 |e|) compared to Pd55 (10.12 |e|). θ(Pd in shell) of Pt2Au23Pd30 is a little more than that of Au25Pd30 since Pt has low electronegativity compared to Au.

ACS Paragon Plus Environment

12

Page 13 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

Figure 5: (a) Calculated As for main reaction and side reaction as a function of average valance electron charge of Pd atoms in cluster shell. (b) Average valance electron charge of Pd atoms in shell of Pd-base binary alloy and AuPd-base ternary alloy as a function of descriptor α based on electronegativity of doped element.

Based on the above analysis, metal elements in 4th, 5th, 6th row of the periodic table are examined systematically as the doped elements to Pd clusters and Au-Pd nanoalloys. It is noted that the atom arrangement of binary and ternary clusters are the same as Au25Pd30 and Pt2Au23Pd30, respectively, which exclude the influence of chemical ordering in order to focus on the influence of electronegativity. Here, we introduce α=(EM/EPd) to represent the effect of relative electronegativity of the doped metal (EM and EPd are the electronegativity of doped metal and Pd, respectively, which can be found in Table S7). Figure 5b gives the correlation between

ACS Paragon Plus Environment

13

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 21

θ(Pd in shell) and α. It is clear that Pd atoms have fewer valence electrons when doping elements with higher electronegativity into pure Pd cluster or Au-Pd nanoalloys and vice versa. Therefore, the exact intrinsic effect of a doped metal on reactivity of Pd-based alloy catalysts is redistribution of θ(Pd in shell) through the different electronegativity from that of Pd.

3.3 Prediction of New Pd-based Nanocatalysts The dash area in Figure 5b corresponds to the optimal θ(Pd in shell). Thus, W and Pb are predicted to be the promising substituted elements for Pd-based binary alloy catalyst. In addition, Mo, W, Pb, Ir, Os, Ru and Rh are selected to substitute Pt for improving Au-Pd ternary alloy nanocatalysts. In order to verify the prediction, Sabatier analysis is performed to estimate As for main/side reaction on Pb25Pd30, W25Pd30, Mo2Au23Pd30, Ru2Au23Pd30, Os2Au23Pd30, Ir2Au23Pd30, Rh2Au23Pd30,

Pb2Au23Pd30 and W2Au23Pd30 systems. Ag25Pd30, Os25Pd30, Rh25Pd30 and

Ag2Au23Pd30, which are not in the optimal region, are also chosen for comparison. The activity of main/side reaction of candidate Pd-based binary nanocatalysts is listed in Table S8. It illuminates that Pb25Pd30 and W25Pd30 has prior activity and selectivity toward direct synthesis of H2O2 compared with other candidates. Moreover, As for main/side reaction of candidate Pdbased ternary nanocatalysts by Sabatier analysis is summarized in Table S9, indicating that Mo2Au23Pd30,

Ru2Au23Pd30,

Os2Au23Pd30,

Ir2Au23Pd30,

Rh2Au23Pd30,

Pb2Au23Pd30

and

W2Au23Pd30 are superior to Ag2Au23Pd30 in both activity and selectivity. The predicted activity of main/side reaction on predicted Pd-based alloys are consistent with the calculated results by DFT calculation data in Table S10 and Figure S7, indicting the robustness of the rule. The results in Table S5,S6,S8,S9 are visualized in Figure 6a and 6b. We established a volcano relationship between α and the activity of main/side reaction and found that α can serve as the activity descriptor for guideline to choose appropriate doped metal to improve the activity and

ACS Paragon Plus Environment

14

Page 15 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

selectivity of Pd-based catalysts. Compared with the previous descriptors such as d-band center65-66 and adsorption energy of reactant or intermedium derived from Sabatier analysis, descriptor α is more relevant to intrinsic effect of dopant and can be conveniently used to predict catalytic activity of Pd-based catalysts. According to Figure 6a and 6b, the best Pd-based binary alloy catalyst is predicted to be the doped metal with electronegativity slightly higher than Pd (α≈1.05-1.1). And the best PdAu-based ternary alloy catalyst is predicted to be the doped metal with electronegativity close to Pd (α≈0.9-1.1). Experimentally, Au-Pd nanocatalysts with a little Pt promoter show better catalytic performance than Pd and Au-Pd,42 fully consistent with our predication. In addition, it is found that all the predicted promising Pd-based catalysts in Figure 6a and 6b can not only improve the activity for direct synthesis of H2O2 (above the red dot line) but also enhance the selectivity (below the blue dot line), compared with Au25Pd30 which has been proved experimentally as an excellent catalyst for direct synthesis of H2O2. Sabatier activity volcano for these tested Pd-based catalysts is shown in Figure S8.

ACS Paragon Plus Environment

15

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 21

Figure 6: (a) As for main reaction and side reaction (Sabatier analysis) as a function of descriptor α based on electronegativity of doped element for Pd-base binary alloy clusters. (b) As for main reaction and side reaction (Sabatier analysis) as a function of descriptor α based on electronegativity of doped element for PdAu-base ternary alloy clusters. The red and blue dot lines represent the As for main and side reaction of Au25Pd30, respectively.

Besides, recent experiments reported that the selectivity of Pd-based alloy nanocatalysts would be improved if doped metal can form oxide surface layer and encapsulates small Pd-rich particles which is active for decomposing and hydrogenating the H2O2.41 According to reducibility of metal, Mo, Ru, Pb and W can form oxide after calcinations as same as Sn. Thus, it can be predicted that Pd-W, Pd-Pb, Au-Pd-W, Au-Pd-Pb, Au-Pd-Mo and Au-Pd-Ru are potentially good candidates for Pd-based alloy nanocatalysts, showing better catalytic performance than the

ACS Paragon Plus Environment

16

Page 17 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

available experimental systems.

4. Conclusions In summary, combining the density functional theory (DFT) calculations and Sabatier analysis, the average valence electron of Pd-shell atoms is identified as the intrinsic factor for the activity and selectivity of the Pd-based nanocatalysts, which can be effectively tuned by the electronegativity of the dopants. The dopants with suitable electronegativity withdraw electrons from Pd atoms to the optimal range to enhance the activity and selectivity simultaneously. Using descriptor α (relative electronegativity between doped metal and Pd), we reproduce the order of the activity and selectivity of the available experimental Pd-based alloy nanocatalysts, and also predict that Pd-W and Pd-Pb are promising substitution for Au-Pd, and W, Pb, Mo as well as Ru are more superior to Pt as a promoter for Au-Pd. Our calculations indicate these designed new Pd-based alloy nanocatalysts show better catalytic performance than the available experimental systems. By this work, we hope to attract the experimental confirmation of the remarkable performance of these designed new Pd-based alloy nanocatalysts for direct synthesis of H2O2 in near future. Overall, our strategy paves the avenue for the rational design of highly efficient new heterogeneous catalysts, which may also find applications beyond direct synthesis of H2O2.

Supporting Information The supporting information contains additional information regarding details of selection of model nanoparticles, details of adsorption free energy calculations, scaling relationships, reaction barrier inequality constraints, DOS analysis and sabatier analysis for predicted catalysts. This material is available free of charge via the internet at http://pubs.acs.org.

ACS Paragon Plus Environment

17

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 21

ACKNOWLEDGMENT This work is supported by the National Natural Science Foundation of China (21576008, 91634116, 91334203, 21273268, 11574340), and "Hundred People Project" from Chinese Academy of Sciences. YG also thanks Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase).

REFERENCES (1) Campos-Martin, J. M.; Blanco-Brieva, G.; Fierro, J. L. G. Angew. Chem., Int. Ed. 2006, 45, 6962–6984. (2) Esmelindro, M. C.; Oestreicher, E. G.; Márquez-Alvarez, H.; Dariva, C.; Egues, S. M. S.; Fernandes, C.; Bortoluzzi, A. J.; Drago, V.; Antunes, O. A. C. J. Inorg. Biochem. 2005, 99, 2054-2061. (3) Hall, A. S.; Kondo, A.; Maeda, K.; Mallouk, T. E. J. Am. Chem. Soc. 2013, 135, 1627616279. (4) Han, Y.; He, Z.; Wang, S.; Li, W.; Zhang, J. Catal. Sci. Technol. 2015, 5, 2630–2639. (5) Morimoto, Y.; Bunno, S.; Fujieda, N.; Sugimoto, H.; Itoh, S. J. Am. Chem. Soc. 2015, 137, 5867-5870. (6) Samanta, C. Appl. Catal., A: Gen. 2008, 350, 133-149. (7) Sato, K.; Noyori, R. Chem. Chem. Ind. 1999, 52, 1166-1169. (8) Shang, H.; Zhou, H.; Zhu, Z.; Zhang, W. J. Ind. Eng. Chem. 2012, 18, 1851–1857. (9) Tan, J.; Zhang, J. S.; Lu, Y. C.; Xu, J. H.; Luo, G. S. AIChE. J. 2012, 58, 1326–1335. (10) Turnwald, S. E.; Lorier, M. A.; Wright, L. J.; Mucalo, M. R. J. Mater. Sci. Lett. 1998, 17, 1305-1307. (11) Hess, W. T.; Kroschwitz, J. I.; Howe-Grant, M. Kirk-Othmer Encyclopediaof Chemical Technology,4th ed, Wiley: New York, 1995; Vol. 13, p 961. (12) Drelinkiewicz, A.; Waksmundzka-Góra, A. J. Mol. Catal. A: Chem. 2006, 246, 167-175. (13) Garcia-Serna, J.; Moreno, T.; Biasi, P.; Cocero, M. J.; Mikkola, J. P.; Salmi, T. O. Green Chem. 2014, 45, 2320-2343. (14) Kosydar, R.; Drelinkiewicz, A.; Ganhy, J. P. Catal. Lett. 2010, 139, 105-113. (15) Choudhary, V. R.; Jana, P. Appl. Catal., A: Gen. 2007, 329, 79-85.

ACS Paragon Plus Environment

18

Page 19 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

(16) Choudhary, V. R.; Samanta, C. Catal. Lett. 2005, 99, 79-81. (17) Ghedini, E.; Menegazzo, F.; Signoretto, M.; Manzoli, M.; Pinna, F.; Strukul, G. J. Catal. 2010, 273, 266–273. (18) Hu, B.; Deng, W.; Li, R.; Zhang, Q.; Wang, Y.; Delplanque-Janssens, F.; Paul, D.; Desmedt, F.; Miquel, P. J. Catal. 2014, 319, 15-26. (19) Kim, J.; Chung, Y. M.; Kang, S. M.; Choi, C. H.; Kim, B. Y.; Kwon, Y. T.; Kim, T. J.; Oh, S. H.; Lee, C. S. ACS Catal. 2012, 2, 1042-1048. (20) Landon, P.; Collier, P. J.; Carley, A. F.; Chadwick, D.; Papworth, A. J.; Burrows, A.; Kiely, C. J.; Hutchings, G. J. Phys. Chem. Chem. Phys. 2003, 34, 1917-1923. (21) Melada, S.; Rioda, R.; Menegazzo, F.; Pinna, F.; Strukul, G. J. Catal. 2006, 239, 422-430. (22) Ouyang, L.; Tian, P. F.; Da, G. J.; Xu, X. C.; Ao, C.; Chen, T. Y.; Si, R.; Xu, J.; Han, Y. F. J. Catal. 2015, 321, 70-80. (23) Pashkova, A.; Greiner, L.; Krtschil, U.; Hofmann, C.; Zapf, R. Appl. Catal., A: Gen. 2013, 464-465, 281-287. (24) Choudhary, V. R.; Samanta, C. J. Catal. 2006, 238, 28-38. (25) Choudhary, V. R.; Jana, P. Appl. Catal., A: Gen. 2009, 352, 35-42. (26) Choudhary, V. R.; Samanta, C.; Choudhary, T. V. Catal. Commun. 2007, 8, 1310-1316. (27) Edwards, J. K.; Freakley, S. J.; Carley, A. F.; Kiely, C. J.; Hutchings, G. J. Acc. Chem. Res. 2014, 47, 845-854. (28) Edwards, J. K.; Pritchard, J. C.; Piccinini, M.; Shaw, G.; He, Q.; Carley, A. F.; Kiely, C. J.; Hutchings, G. J. J. Catal. 2012, 292, 227-238. (29) Edwards, J. K.; Solsona, B.; Ntainjua N, E.; Carley, A. F.; Herzing, A. A.; Kiely, C. J.; Hutchings, G. J. Science 2009, 323, 1037-1041. (30) Edwards, J. K.; Thomas, A.; Carley, A. F.; Herzing, A. A.; Kiely, C. J.; Hutchings, G. J. Green Chem. 2008, 10, 388-394. (31) Landon, P.; Collier, P. J.; Papworth, A. J.; Kiely, C. J.; Hutchings, G. J. Chem. Commun. 2002, 18, 2058-2059. (32) Sankar, M.; He, Q.; Morad, M.; Pritchard, J.; Freakley, S. J.; Edwards, J. K.; Taylor, S. H.; Morgan, D. J.; Carley, A. F.; Knight, D. W.; Kiely, C. J.; Hutchings, G. J. ACS Nano 2012, 6, 6600-6613. (33) Chen, M.; Kumar, D.; Yi, C. W.; Goodman, D. W. Science 2005, 310, 291-293. (34) Lopez-Sanchez, J. A.; Dimitratos, N.; Glanville, N.; Kesavan, L.; Hammond, C.; Edwards, J. K.; Carley, A. F.; Kiely, C. J.; Hutchings, G. J. Appl. Catal., A: Gen. 2011, 391, 400-406.

ACS Paragon Plus Environment

19

ACS Catalysis

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 21

(35) Ouyang, L.; Da, G. J.; Tian, P. F.; Chen, T. Y.; Liang, G. D.; Xu, J.; Han, Y. F. J. Catal. 2014, 311, 129-136. (36) Bernardotto, G.; Menegazzo, F.; Pinna, F.; Signoretto, M.; Cruciani, G.; Strukul, G. Appl. Catal., A: Gen. 2009, 358, 129-135. (37) Liu, Q.; Bauer, J. C.; Schaak, R. E.; Lunsford, J. H. Appl. Catal., A: Gen. 2008, 339, 130– 136. (38) Xu, J.; Ouyang, L.; Da, G. J.; Song, Q. Q.; Yang, X. J.; Han, Y. F. J. Catal. 2012, 285, 7482. (39) Gu, J.; Wang, S.; He, Z.; Han, Y.; Zhang, J. Catal. Sci. Technol. 2015, 6, 809-817. (40) Maity, S.; Eswaramoorthy, M. J. Mater. Chem. A 2016, 4, 3233-3237. (41) Freakley, S. J.; He, Q.; Harrhy, J. H.; Lu, L.; Crole, D. A.; Morgan, D. J.; Ntainjua, E. N.; Edwards, J. K.; Carley, A. F.; Borisevich, A. Y.; Kiely, C. J.; Hutchings, G. J. Science 2016, 351, 965-968. (42) Edwards, J. K.; Pritchard, J.; Lu, L.; Piccinini, M.; Shaw, G.; Carley, A. F.; Morgan, D. J.; Kiely, C. J.; Hutchings, G. J. Angew. Chem., Int.Ed. 2014, 53, 2381-2384. (43) Verdaguer-Casadevall, A.; Deiana, D.; Karamad, M.; Siahrostami, S.; Malacrida, P.; Hansen, T. W.; Rossmeisl, J.; Chorkendorff, I.; Stephens, I. E. L. Nano Lett. 2014, 14, 16031608. (44) Rankin, R. B.; Greeley, J. ACS Catal. 2012, 2, 2664-2672. (45) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I.; Corso, A. D.; Gironcoli, S. d., Fabris, S; Fratesi, G.; Gebauer, R.; Gerstmann, U.; Gougoussis, C.; Kokalj, A.; Lazzeri, M.; Martin-Samos, L.; Marzari, N.; Mauri, F.; Mazzarello, R.; Paolini, S.; Pasquarello, A.; Paulatto, L.; Sbraccia, C.; Scandolo, S.; Sclauzero, G.; Seitsonen, A. P.; Smogunov, A.; Umari, P.; Wentzcovitch, R. M. J. Phys. Conden. Matter. 2009, 21, 395502. (46) Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1996, 77, 3865-3868. (47) Vanderbilt, D. Phys. Rev. B 1990, 41, 7892-7895. (48) Nilekar, A. U.; Mavrikakis, M. Surf. Sci. 2008, 602, L89–L94. (49) Ojifinni, R. A.; Froemming, N. S.; Gong, J.; Pan, M.; Kim, T. S.; White, J. M.; Henkelman, G.; Mullins, C. B. J. Am. Chem. Soc. 2008, 130, 6801-6812. (50) Atkins, P. W.; De Paula, J. Physical Chemistry, 9th ed., Oxford University Press: Oxford, U.K., 2010; p 472, 922, 924, 933. (51) Henkelman, G.; Uberuaga, B. P.; Jónsson, H. J. Chem. Phys. 2000, 113, 9901-9904. (52) Ford, D. C.; Nilekar, A. U.; Xu, Y.; Mavrikakis, M. Surf. Sci. 2010, 604, 1565–1575.

ACS Paragon Plus Environment

20

Page 21 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

(53) Ham, H. C.; Hwang, G. S.; Han, J.; Nam, S. W.; Lim, T. H. J. Phys. Chem. C 2009, 113, 12943-12945. (54) Ham, H. C.; Hwang, G. S.; Han, J.; Nam, S. W.; Lim, T. H. J. Phys. Chem. C 2010, 114, 14922-14928. (55) Li, J.; Ishihara, T.; Yoshizawa, K. J. Phys. Chem. C 2011, 115, 25359-25367. (56) Li, L.; Yoshizawa, K. Catal. Today 2014, 248, 142–148. (57) Staykov, A.; Kamachi, T.; Ishihara, T.; Yoshizawa, K. J. Phys. Chem. C 2008, 112, 1950119505. (58) Todorovic, R.; Meyer, R. J. Catal. Today 2011, 160, 242-248. (59) Bligaard, T.; Nørskov, J. K.; Dahl, S.; Matthiesen, J.; Christensen, C. H.; Sehested, J. J. Catal. 2004, 224, 206-217. (60) Ferrin, P.; Simonetti, D.; Kandoi, S.; Kunkes, E.; Dumesic, J. A.; Nørskov, J. K.; Mavrikakis, M. J. Am. Chem. Soc. 2009, 131, 5809-5815. (61) Nørskov, J. K.; Bligaard, T.; Logadottir, A.; Bahn, S.; Hansen, L. B.; Bollinger, M.; Bengaard, H.; Hammer, B.; Sljivancanin, Z.; Mavrikakis, M.; Xu, Y.; Dahl, S.; Jacobsen, C. J. H. J. Catal. 2002, 209, 275-278. (62) Edwards, J. K.; Carley, A. F.; Herzing, A. A.; Kiely, C. J.; Hutchings, G. J. Faraday Discuss. 2008, 138, 225-239. (63) Ishihara, T.; Ohura, Y.; Yoshida, S.; Hata, Y.; Nishiguchi, H.; Takita, Y. Appl. Catal., A: Gen. 2005, 291, 215-221. (64) Wilson, N. M.; Flaherty, D. W. J. Am. Chem. Soc. 2016, 138, 574-586. (65) Hammer, B.; Norskov, J. K. Nature 1995, 376, 238-240. (66) Kitchin, J. R.; Nørskov, J. K.; Barteau, M. A.; Chen, J. G. J. Chem. Phys. 2004, 120, 10240-10246.

TOC Figure

ACS Paragon Plus Environment

21