Coverage Effect on the Activity of the Acetylene Semihydrogenation

Coverage Effect on the Activity of the Acetylene Semihydrogenation over Pd–Sn Catalysts: A Density Functional Theory Study ...
2 downloads 0 Views 3MB Size
Article Cite This: J. Phys. Chem. C 2018, 122, 6005−6013

pubs.acs.org/JPCC

Coverage Effect on the Activity of the Acetylene Semihydrogenation over Pd−Sn Catalysts: A Density Functional Theory Study Jiubing Zhao, Shenjun Zha, Rentao Mu, Zhi-Jian Zhao,* and Jinlong Gong* Key Laboratory for Green Chemical Technology of Ministry of Education, School of Chemical Engineering and Technology, Tianjin University; Collaborative Innovation Center of Chemical Science and Engineering, Tianjin 300072, China

Downloaded via IOWA STATE UNIV on January 6, 2019 at 21:42:16 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: The existence of acetylene impurities in ethylene feedstock is harmful to downstream polymerization reactions. The removal of acetylene can be achieved via semihydrogenation reaction, which is normally catalyzed by Pd-based catalysts. This paper describes the coverage effect on the activity of acetylene hydrogenation reactions over Pd and PdSn alloy surfaces. High-coverage models are presented to construct coverage-dependent adsorption energies of C2H2, C2H4, and H2 on Pd(111) and Pd3Sn(111) surfaces. It has been validated that the downshift of d-band center caused by preadsorbed molecules makes the adsorption weaker along with the increase of coverage, and the geometric effect can be neglected. An iterative method has been applied to predict surface coverages of reaction intermediates. Previous calculations with low-coverage models indicate that alloying Pd with late or post-transition metals, in general, enhances ethylene selectivity, accompanied with lower hydrogenation activity. However, by applying a high-coverage model, we show that the predicted hydrogenation barriers are comparable over Pd(111) and Pd3Sn(111) surfaces. single-site heterogeneous catalysts.18,19 For example, the Lindlar catalyst is widely employed in industries for hydrogenation reactions, where lead or sulfur was used to tune the hydrogenation ability of Pd. Currently, the most commonly used industrial catalyst for acetylene semihydrogenation is Agmodified Pd catalyst.20,21 Other Pd alloy catalysts such as PdCu,22 PdGa,23 and PdNi24 have also shown improved catalytic properties in acetylene semihydrogenation compared with pure Pd. In most cases, these alloying components are late and post-transition metals, which shift down the d-band center of Pd in the alloy, thus enhance the desorption of the desired product ethylene, and prevent its overhydrogenation.25−27 Actually, the difference between the hydrogenation barrier and the desorption barrier of ethylene is defined as the crucial factor influencing the selectivity.28,29 During the last two decades, density functional theory (DFT) has become a powerful tool for studying heterogeneous catalytic processes.30−32 Normally, the atomic scale analysis of elementary steps offered by DFT studies is very hard to be performed experimentally.33 In addition, the mechanistic understanding offered by DFT calculations served as guidelines for designing new catalysts with higher activity and/or selectivity. Using DFT calculations, Nørskov and co-workers

1. INTRODUCTION Ethylene is one of the most important feedstocks in the chemical industry.1 In 2016, the global production of ethylene was over 150 million tones, with most of it going to the plastic industry by polymerization reactions.2 Normally, the production of polyethylene requires an acetylene impurity content in the ethylene feedstock below 5 ppm by volume to prevent catalyst poisoning and ensure the purity of the produced polyethylene.3 However, ethylene feedstock produced by cracking can contain more than 2 vol % acetylene. Thus, it must be removed before polymerization processes.4 Most of the existing and under-construction ethylene plants use catalytic hydrogenation method to refine the cracking gas to produce the desired polymer-grade ethylene product.5 Pd-based catalysts have been recognized to selectively hydrogenate acetylene to ethylene.6,7 Previous studies have proposed that hydrogenation occurs via a Horiuti−Polanyi mechanism.8−10 Namely, H2 is dissociated on Pd,11 and the hydrogen adatom then diffuses toward the adsorbed acetylene on Pd, forming new C−H bonds. In addition, a concerted hydrogenation mechanism has also been reported for sequential hydrogenation,12 which we did not take into account. It has been reported that there were two major issues in this process:6,13−17 overhydrogenation to form ethane and oligomerization to generate green oil. The selectivity of hydrogenation reactions can be tuned by alloying Pd with another metal component or by formation of © 2018 American Chemical Society

Received: November 19, 2017 Revised: February 25, 2018 Published: March 1, 2018 6005

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

Article

The Journal of Physical Chemistry C screened about 70 different bimetallic compounds for acetylene semihydrogenation reactions.28 The optimal catalyst suggested in this study is considered to be located within a window with a suitable binding strength of the adsorbates,34−36 which is equivalent to the well-known Sabatier principle.37 When the adsorption of acetylene is too weak, the activity is low because of the poor activation of the reactants. However, the product ethylene easily desorbs from catalyst surface to increase the selectivity. As we know, too strong adsorption is easy to activate the reactants accompanying with high activity but too strong adsorption normally activates a series of side reactions, for example, overhydrogenation of ethylene, leading to low selectivity. Indeed, previous calculations have indicated that Pd alloy catalysts, for example, PdAg,38 PdFe,39 PdGa,40 and PdIn,41 usually improve the selectivity to ethylene simultaneously suppress the activity of catalysts compared with pure Pd.39,41−44 However, recent studies by Esmaeili et al. reported that the addition of Sn into Pd enhances both activity and selectivity for acetylene semihydrogenation reaction.45,46 To gain a deeper understanding on this phenomenon, we performed DFT calculations over a series of PdSn alloy surfaces, including (111) of Pd and Pd3Sn surfaces as well as PdSn(010) surface. By applying the coverage-dependent adsorption energy together with Langmuir adsorption isotherm, the adsorbate coverages have been predicted and corresponding barriers at this coverage have been calculated. Our calculation indicates that a high-coverage environment is a key factor to correctly describe the relative activity over calculated PdSn alloy surfaces. Herein, our study shed light on how the coverage effect influences the reaction kinetics.

Figure 1. Top and side views of (a) Pd(111), (b) Pd3Sn(111), and (c) PdSn(010) model surfaces. The blue and pink balls denote Pd and Sn atoms, respectively. This definition is used throughout this paper.

Sads = 0.7Sgas − 3.3R

(1)

where Sads is the entropy of the adsorbate; Sgas is the entropy of the corresponding gas-phase species, obtained from thermodynamic tables;54 and R is the ideal gas constant. We assume no entropy change for surface reaction because the harmonic approximation is poor at determining entropy in low-frequency vibrational modes. It has been observed that the adsorption energy can be strongly affected by the presence of other coadsorbates on the surface,55 which has been partially considered in our work. To investigate the coverage effect on adsorption, differential adsorption energies of the reactants C2H2, C2H4, and H are calculated at coverages of 1/16, 1/8, 3/16, 1/4, 5/16, 3/8, 7/ 16, and 1/2 ML at adsorption sites, which are performed on a (4 × 4) model slab. The differential adsorption free energy of acetylene, ethylene, and hydrogen is defined as

2. MODELS AND COMPUTATIONAL DETAILS Periodic DFT calculations were carried out with the Vienna Ab initio Simulation Package.47 The calculations employed the generalized-gradient approximation (GGA) in the form of the Bayesian error estimation functional with van der Waals corrections.48,49 The interactions between the atomic cores and electrons were described by the projector augmented wave method.50 The valence wave functions were expanded by a plane wave with a cutoff energy of 400 eV. A Methfessel− Paxton smearing with a 0.15 eV width was employed to speed up the convergence, and the total energies were evaluated by extrapolating to zero broadening.51 The three slab models employed in this study are sketched in Figure 1. For Pd and Pd3Sn, we chose the flat surface (111) as the sufficient evaluated surface to investigate the properties of catalysts. For Pd/Sn alloy at a 1:1 ratio, we calculated its surface energies and found that (010) is the relative stable and ordered surface among other low-miller index surfaces (see Table S1 in the Supporting Information). Meanwhile, the segregation energies of PdSn catalysts were calculated (Table S2), and no segregation case was found from PdSn catalysts. The thickness of the employed slabs is four layers (except for five layers of PdSn(010)), with top two layers relaxed on each surface. Optimized geometries were found when the force on each relaxed atom was less than 0.02 eV/Å. A k-point mesh of 3 × 3 × 1 was used in a 4 × 4 unit cell of each model slab. Entropic contributions and zero-point energies were taken into account, which converts the DFT-calculated energy into free energy.52 The entropy of the adsorbate was calculated according the Campbell’s method with the following equation:53

Gads(θx) = Gtotal, N − Gtotal, N − 1 − Ggas

(2)

where Gtotal,N and Gtotal,N−1 are the total Gibbs free energy after adsorption of N and N − 1 adsorbates, respectively, and Ggas is the Gibbs free energy of the adsorbate in the gas phase. On PdSn(010), the adsorption free energies are always positive for both kinds of adsorbates even at the lowest coverage (1/16 ML) considered in our calculations. Thus, we just consider the coverage effect on Pd and Pd3Sn surfaces. All kinetic and thermodynamic analysis are carried out at 373 K with the 1, 0.1, and 8.9 atm partial pressures for H2, C2H2, and C2H4, respectively, which are close to the experimental conditions.23

3. RESULTS AND DISCUSSION 3.1. Adsorption of Acetylene and Ethylene on Pd/Sn Surfaces at a Low Coverage. The adsorption energies for C2H2 and C2H4 and their adsorption geometries on Pd(111), Pd3Sn(111), and PdSn(010) are shown in Table 1 and Figure 2, respectively. The adsorption energies and Gibbs free Table 1. Adsorption Energies (Eads/eV) and Gibbs Free Adsorption Energies (Gads/eV) of C2H2 and C2H4 on Surfaces Pd(111) C2H2 C2H4 6006

Pd3Sn(111)

PdSn(010)

Eads

Gads

Eads

Gads

Eads

Gads

−1.83 −0.81

−1.42 −0.40

−1.08 −0.47

−0.68 −0.05

−0.18 −0.27

0.19 0.15

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

Article

The Journal of Physical Chemistry C

Table 2. Single-Point Energy of Gas-Phase Acetylene on the (111) Surface of Pd and Pd3Sn and the Ratio of the SinglePoint Energy of N Molecular Acetylene to One Molecular Acetylene Pd(111) coverage/ML

Esingle/eV

1/16 2/16 3/16 4/16 5/16 6/16

−18.43 −36.88 −54.74 −72.96 −91.93 −109.99

Pd3Sn(111) ratio

coverage/ML

Esingle/eV

ratio

2.00 2.97 3.96 4.99 5.97

1/16 2/16 3/16 4/16 5/16 6/16

−18.93 −37.76 −56.60 −75.42 −94.21 −113.83

2.00 2.99 3.98 4.98 6.01

Figure 2. Top and side views of C2H2 and C2H4 adsorption geometries on Pd(111), Pd3Sn(111), and PdSn(010) model surfaces.

Figure 4. Plots of the Gibbs adsorption energy of acetylene on Pd(111) and Pd3Sn(111) against the d-band centers. The coverage of each point is performed in the plot.

Scheme 1. Iterative Method To Identify the Coverage of C2H2, C2H4 and H

Figure 3. (a,b) Trend of the changes in differential adsorption free energies of C2H2 and C2H4 with the coverage on Pd(111) and Pd3Sn(111) surfaces. The coverage dependence of the differential adsorption free energy of C2H2 with one molecular C2H4 coadsorbing on Pd(111) surface is illustrated by the dashed red line in (a). The relationship between the coverage and the differential adsorption free energy of C2H4 with xC2H2 coadsorbing on Pd(111) surface is shown by the blue dashed line.

adsorption energies of C2H2 and C2H4 on all sites are listed in Table S2. Under ultrahigh vacuum (UHV) conditions, acetylene preferentially adsorbs at 3-fold hollow sites on Pd(111),56 which was also found to be energetically most favorable in this study. While the most stable site for C2H4 adsorption is bridge, with two newly formed C−Pd bonds. At 6007

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

Article

The Journal of Physical Chemistry C

Figure 5. (a,b) Trend of changes in the differential adsorption free energies of H with the coverage on Pd(111) (6C2H2 coadsorbed) and Pd3Sn(111) (3C2H2 and 1C2H4 coadsorbed). (c,d) Trend of the changes in differential adsorption free energies of C2H2 and C2H4 with the coverage on Pd(111) (2H coadsorbed) and Pd3Sn(111) (3H coadsorbed) surfaces.

isolated Pd atoms, further weakening the binding strength of adsorbates by forcing them to interact with less Pd atoms. 3.2. Coverage-Dependent Adsorption Energy for Single Adsorbate. As shown in Table 1, the calculated adsorption energies of C2H2 and C2H4 are very strong according to the low-coverage model. If we simply apply the Langmuir adsorption model with a constant adsorption energy predicted by a low-coverage model, that is, over 1 eV of acetylene over Pd surfaces, all active sites will be occupied by the strong binding C2H2, which is much higher than the reported saturated coverage of acetylene, 0.46 ML over Pd(111) measured by the low-energy electron diffraction method.59 Because the adsorption energy, in general, depends on the coverage of adsorbates, understanding the interactions between adsorbed acetylene over Pd and PdSn alloy surfaces is important to gain a correct description of surface species population under reaction conditions. The calculated coverage-dependent differential adsorption energies are shown in Figure 3a,b by solid lines. The differential adsorption free energy is observed to be positively correlated with the coverage above a threshold and a constant when the coverage is below this threshold. The point where the differential adsorption free energy is 0 eV indicates thermodynamic equilibration between adsorbed and gas-phase species, corresponding to the saturation coverage. On the Pd(111) surface, the calculated saturation coverage of acetylene is about 0.40 ML, in agreement with experimental results of 0.46 ML.59 It is also found that the saturation coverage of ethylene is about 0.27 ML in this work, which is also very close to the previous experimental result of 0.25 ML.60 Note that the pressure used in our calculation is much higher than those of the surface science studies performed in refs 59 and 60. Therefore, the

Figure 6. Coverage of reactants on the (111) surface of Pd and Pd3Sn.

1/16 and 1/4 ML coverages, the calculated binding energies are −0.81 and −0.74 eV, respectively. Previous temperatureprogrammed desorption experiment measured the adsorption energy of ethylene under UHV conditions to be −0.61 eV on Pd(111) at a 1/4 ML coverage, which is slightly weaker than our DFT result at the same coverage,57 due to the well-known tendency of GGA method to overestimate binding energies.58 With the addition of Sn atoms to the Pd catalyst, the binding strength of adsorbates becomes weaker. The introduction of Sn tunes the electronic structure of Pd by downshifting its d-band center, which weakens the adsorption of adsorbates. Indeed, linear correlation between d-band center and adsorption energy has been observed (Figure S1 in the Supporting Information) when the acetylene binds with the same π adsorption mode on these terrace surfaces. Moreover, the surface-existing Sn atoms divided the complete Pd surface into small ensembles or 6008

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

Article

The Journal of Physical Chemistry C

might be caused by the following two factors: (a) geometric effect due to the increased steric hindrance and (b) electronic effect due to the shifts of the d-band center caused by an adsorbed molecule. Because the steric hindrance mainly comes from the repulsion between adsorbates, single-point calculations were done for the adsorbates only, that is, without a metal slab, with the optimized geometry over metal surfaces, which ignores the electronic effect between adsorbate and metal slab. The calculated single-point energies are shown in Table 2. Clearly, the ratio of the single-point energy of N acetylene molecules to one acetylene molecule is nearly N, indicating nearly no interactions between adsorbates; thus, the geometric effect can be neglected. Test calculations further indicate that there are no interactions for two molecular C2H2 with each other until the distance between them is less than 2.13 Å and the nearest adsorbates on the 0.44 ML model are 3.61 Å apart from each other. Further analysis indicates that the weakened binding strength is mainly caused by the electronic effect. A perfect linear relationship has been observed between the d-band center of the Pd atoms, which will host the next adsorbate, and its corresponding differential binding energy, as shown in Figure 4. Therefore, the change of adsorption energy with the coverage of adsorbates is mainly due to the downshift of the d-band center caused by the preadsorbed molecule, making the adsorption of the next adsorbate weaker. 3.3. Competing Adsorption of H, C2H2, and C2H4. Under reaction conditions, both acetylene and ethylene and other surface intermediates coadsorbed over catalyst surfaces. To identify the cross-interactions between different adsorbates, we first set up the high-coverage model, which accounts for the competing adsorption of C2H2 and C2H4 on the Pd(111) surface, as shown in Figure 3a by dashed lines. The dashed blue line shows the trend of the differential adsorption energy of C2H4 with coadsorbed C2H2 or vice versa against the coverage of C2H4. It almost coincides with the solid blue line as well as the case of the red dashed line, which suggests that the crossinteraction between C2H2 and C2H4 is nearly the same as the C2Hx (x = 2 or 4) self-interaction. This result is not surprising because the repulsive interaction mainly induced by the downshift of the d-band center via similar C−Pd interactions because of the existence of the neighboring adsorbates, either C2H4 or C2H2. Therefore, it is safe for us to employ the interaction models built via single-type adsorbates to describe the case of multiadsorbate systems, even for the environmental coverages, which are not covered by our explicit DFT calculations. On the basis of the relationship obtained in section 3.2, we can use the Langmuir adsorption model to estimate the coverage of acetylene and ethylene on Pd(111) and Pd3Sn(111) surfaces before any reaction occurs. First, we consider the

Figure 7. (a) Gibbs free energy profiles of acetylene hydrogenation to ethane on Pd(111) and Pd3Sn(111) with high-coverage models and PdSn(010). (b) Reaction pathways on Pd(111) and Pd3Sn(111) with a low-coverage model. (c) Comparison of the reaction barrier of first acetylene hydrogenation between high-coverage (0.35 ML C2H2 and 0.07 ML H on Pd(111) and 0.19 ML C2H2 and 0.13 ML H on Pd3Sn(111)) and low-coverage issue.

agreement might also be due to the accidental error cancellation caused by the pressure difference, DFT errors, and the approximation on the treatment of desorption entropy change. As expected, the differential adsorption free energy becomes less negative with increasing the coverage of adsorbates on Pd(111) and Pd3Sn(111) surfaces, as shown in Figure 3, which

Table 3. Free Energy Barrier (ΔG⧧/eV) and Reaction Free Energy (ΔG/eV) of Hydrogenation Reactions on Surfaces

Pd(111) Pd3Sn(111) PdSn(010)

θhigh θlow θhigh θlow θ

TS1

TS2

TS3

TS4

C2H2 + H

C2H3 + H

C2H4 + H

C2H5 + H

ΔG⧧

ΔG

ΔG⧧

ΔG

ΔG⧧

ΔG

ΔG⧧

ΔG

0.92 1.06 0.88 1.27 0.63

0.25 0.31 −0.18 0.24 −0.96

0.36 1.01 0.46 0.77 0.53

−0.67 −0.12 −0.71 −0.40 −0.95

0.77 0.98 0.75 1.09 0.85

−0.05 0.43 −0.01 0.44 0.07

0.77 0.90 0.64 0.74 0.32

−0.91 −0.49 −1.29 −0.83 −1.74

6009

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

Article

The Journal of Physical Chemistry C

coadsorbed H, that is, 0.07 ML, to identify whether the coadsorbed H would influence the coverage effect on C2H2 adsorption (coverage of C2H4 is too low to be calculated). As illustrated in Figure 5c, a linear relationship between adsorption energy and coverage (for the range larger than 0.19 ML) of C2H2 or C2H4 is established, and the function is shown in section S3 in the Supporting Information. The coverages of H, C2H2, and C2H4 based on multi-adsorbed species model are 0.07, 0.35, and 0.00 ML, respectively, on Pd(111), as shown in Figure 6. Note that the newly predicted coverage is very close to the input values, and these converged values are used in section 3.4. A similar procedure can be done to obtain the coverage over Pd3Sn(111) (Figure 5b,d). 3.4. Acetylene Hydrogenation Mechanism with Coverage Effect on Pd and PdxSny Alloy Surfaces. It has been suggested that a stronger C binding energy is expected in the case of a higher d-band center.30 After carefully analyzing the origin of coverage effect in section 3.2, the d-band center was affected not only by metal30,61 but also by the coadsorbed species on surfaces. Furthermore, the local environment around the reaction active sites and reactant coverage will have a meaningful impact on the intrinsic hydrogenation kinetics.21 Therefore, the predicted coverage of reactants gained in section 3.3 was applied in the analysis of acetylene hydrogenation reactions. On Pd(111), five C2H2 were preadsorbed as a reaction environment, and two C2H2 and one H were preadsorbed on Pd3Sn(111). Here, we follow the Horiuti−Polanyi mechanism for the acetylene hydrogenation. The reaction begins with the adsorption of acetylene and hydrogen, which subsequently combine to form vinyl. Vinyl reacts with a second adsorbed hydrogen atom to form ethylene. Adsorbed ethylene may either desorb to the gas phase or react further to form ethane. The high-coverage model was constructed by placing a number of coadsorbed CHCH, CH2CH2, and H in the unit cell, which is close to the coverage predicted in section 3.3. Under reaction conditions, other C2Hx species, for example, CCH2, have been reported to coexist during the reaction.62 We tested the influence of coadsorbed CCH2 on the corresponding hydrogenation barrier by replacing all coadsorbed CHCH by CCH2. The obtained CHCH hydrogenation barrier with coadsorbed CCH2 is 0.87 eV, very close to the value (0.92 eV) with coadsorbed CHCH over Pd(111). This is not surprising because the adsorbate−adsorbate interaction was mainly induced by the electronic effect, as discussed in section 3.2, and the binding mode is similar for CHCH and CCH2 over a 3-fold hollow site. Therefore, we did not explicitly include CCH2 in our high-coverage model. For CHCH2, we expect that its coverage should be low because of the following fact: the first hydrogenation step is expected to be the rate-determining step because of its high barrier. Therefore, the second hydrogenation step should be in quasi-equilibrium. Because the second hydrogenation step is strongly exothermic (Figure 7a), the coverage of the reactant (CHCH2*) should be much lower than the coverage of the product (CH2CH2*). Note that CH2CH*2 coverage is already low according to our analysis in section 3.3 (Figure 6). The corresponding energy profiles based on the coverage of reactants on the (111) surfaces of Pd and Pd3Sn and PdSn(010) are shown in Figure 7a, and the energy profiles without coverage effect are shown in Figure 7b. All the reaction free energies and free energy barriers are listed in Table 3. The weak adsorption of acetylene on PdSn(010) induces low

competing adsorption between C2H2 and C2H4. The coverages of C2H2 and C2H4 were calculated according to the following equations: 1 θ = * 1+K ( θ ) p C2H 2 * C2H 2 + K C2H4(θ*)pC2H4 θC2H2 =

θC2H4 =

(3)

K C2H2(θ )pC H * 2 2 1 + K C2H2(θ )pC H + K C2H4(θ )pC H * 2 2 * 2 4

(4)

K C2H4(θ )pC H * 2 4 1 + K C2H2(θ )pC H + K C2H4(θ )pC H * 2 2 * 2 4

(5)

θ* is the free site coverage; θC2H2 and θC2H4 are the coverages of C2H2 and C2H4, respectively, px is the partial pressure of the adsorbate x (x = C2H2 and C2H4), and Kx(θ*) is the coveragedependent equilibrium constant of x (x = C2H2 and C2H4) adsorption, which can be related to the differential adsorption free energies of adsorbates: ⎡ Gads,C H (θ ) ⎤ 2 2 * ⎥ K C2H2(θ ) = exp⎢ − * ⎢⎣ ⎥⎦ RT

(6)

⎡ Gads,C H (θ ) ⎤ 2 4 * ⎥ K C2H4(θ ) = exp⎢ − * ⎢⎣ ⎥⎦ RT

(7)

where R is the ideal gas constant and T is the reaction temperature. Gads,C2H2(θ*) and Gads,C2H4(θ*) are the differential adsorption free energies of C2H2 and C2H4, respectively, against free site coverage, which were fitted by the values shown in Figure 3. The corresponding fitted functions of Gads,C2H2(θ*) and Gads,C2H4(θ*) of Pd(111) and Pd3Sn(111) surfaces are shown in section S3 in the Supporting Information. By applying the coverage-dependent adsorption energy (see details in Supporting Information section S3), we obtained the coverage of C2H2 (0.43 ML) and C2H4 (0.00 ML) on Pd(111) and C2H2 (0.20 ML) and C2H4 (0.07 ML) on Pd3Sn(111). Because the reactants of acetylene hydrogenation contain H, we should also consider the coverage of H on each surface. An iterative strategy has been proposed in this study to calculate the coverages of C2H2, C2H4, and H, by constructing a highcoverage model of H or C2Hx with coexistence of fixed C2Hx or H coverage (the extrinsic coverage gained by last iteration). The detail of the iterative method is shown in Scheme 1. Our iterative method starts from the results of coverages of C2H2 and C2H4 on the (111) surfaces of Pd and Pd3Sn because of their strong binding nature over Pd and Pd3Sn. On Pd(111), we first carried out the adsorption of six molecular C2H2 (0.375 ML) on the surface as the reference of hydrogen adsorption (reference models are shown in Figure 5a). As shown in Figure 5a, we obtain the H adsorption energy as a function of its coverage, with six coadsorbed C2H2 (see section S3 in the Supporting Information). Combining the above relationship of C2H2 and C2H4, the predicted coverages of H, C2H2, and C2H4 are 0.07, 0.37, and 0.00 ML on Pd(111). Recall that the coverage-dependent adsorption energy was built under conditions of θH = 0 ML, θC2H2 = 0.43 ML, and θC2H4 = 0.00 ML, which are different from the final conditions. Because a non-negligible H coverage has been predicted, in the next step, the interaction models were developed with 6010

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

The Journal of Physical Chemistry C



surface coverage, and our (4 × 4) unit cell, corresponding to 1/ 16 ML coverage, is suitable to describe it. Comparing the reaction barrier of each step in acetylene hydrogenation at high coverage and low coverage, hydrogenation reaction becomes much easier with weakened adsorption of reactants at high coverage. Note that kinetic Monte Carlo simulations have also found the decrease of hydrogenation barrier with increasing coverage.21 Previous calculations indicated that alloying late or posttransition metals, for example, PtSn,42 PtGa,63 PdAg,38 and PdIn,41 in general, lowers the d-band center of Pd or Pt and thus suppresses its hydrogenation/dehydrogenation activity. Our calculation at a low coverage is consistent with previous findings (Figure 7b).64 However, this is not the case in our high-coverage calculations. The barrier of the rate-determining step on Pd(111) at a high coverage was 0.92 eV, only 0.14 eV lower than the value calculated at a low coverage. However, the high-coverage barrier on Pd3Sn(111) is only 0.88 eV, much lower than the low-coverage case (1.27 eV). Clearly, at a high coverage, the first hydrogenation barrier becomes comparable on both Pd(111) and Pd3Sn(111). Because higher reaction barrier reflects lower rate constants, the qualitative activity trends in our work considering the coverage effect show slightly increasing activity with Sn addition, consistent with the experimental observations.45 According to Figure 7b,c, lowcoverage calculations, which have been applied for most of the previous DFT studies,38,65 could not capture this important feature because of the relatively high hydrogenation barrier over Pd3Sn(111). Our high-coverage model also predicts similar barriers for the hydrogenation of ethylene to ethane on Pd(111) and Pd3Sn(111) (Figure S3 in the Supporting Information).

Article

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.7b11394. Surface energy data; adsorption properties of acetylene and ethylene on each surface; fitted coverage-dependent adsorption energies; high-coverage data; coverage effect on Pd(100) and Pd3Sn(100); additional discussion of reaction barrier on each surface; and coverage effect on Pd(100) and Pd3Sn(100) (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (Z.-J.Z.). *E-mail: [email protected] (J.G.). ORCID

Zhi-Jian Zhao: 0000-0002-8856-5078 Jinlong Gong: 0000-0001-7263-318X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS



REFERENCES

This work was financially supported by the National Natural Science Foundation of China (nos. 21506149, 21525626, and 21676181) and the Program of Introducing Talents of Discipline to Universities (no. B06006).

(1) Chen, Y.; Vlachos, D. G. Hydrogenation of Ethylene and Dehydrogenation and Hydrogenolysis of Ethane on Pt (111) and Pt (211): A Density Functional Theory Study. J. Phys. Chem. C 2010, 114, 4973−4982. (2) Wang, K. L.-C.; Li, H.; Ecker, J. R. Ethylene Biosynthesis and Signaling Networks. Plant Cell 2002, 14, S131−S151. (3) Yang, B.; Burch, R.; Hardacre, C.; Headdock, G.; Hu, P. Origin of the Increase of Activity and Selectivity of Nickel Doped by Au, Ag, and Cu for Acetylene Hydrogenation. ACS Catal. 2012, 2, 1027−1032. (4) Sheth, P. A.; Neurock, M.; Smith, C. M. First-Principles Analysis of the Effects of Alloying Pd with Ag for the Catalytic Hydrogenation of Acetylene-Ethylene Mixtures. J. Phys. Chem. B 2005, 109, 12449− 12466. (5) He, Y.; Fan, J.; Feng, J.; Luo, C.; Yang, P.; Li, D. Pd Nanoparticles on Hydrotalcite as an Efficient Catalyst for Partial Hydrogenation of Acetylene: Effect of Support Acidic and Basic Properties. J. Catal. 2015, 331, 118−127. (6) Jin, Y.; Datye, A. K.; Rightor, E.; Gulotty, R.; Waterman, W.; Smith, M.; Holbrook, M.; Maj, J.; Blackson, J. The Influence of Catalyst Restructuring on the Selective Hydrogenation of Acetylene to Ethylene. J. Catal. 2001, 203, 292−306. (7) Mei, D.; Sheth, P.; Neurock, M.; Smith, C. First-Principles-Based Kinetic Monte Carlo Simulation of the Selective Hydrogenation of Acetylene over Pd(111). J. Catal. 2006, 242, 1−15. (8) Horiuti, J.; Polanyi, M. A Catalysed Reaction of Hydrogen with Water. Nature 1933, 132, 819. (9) Horiuti, J.; Polanyi, M. Catalysed Reaction of Hydrogen with Water and the Nature of over-Voltage. Nature 1933, 132, 931. (10) Horiuti, J.; Polanyi, M. Catalytic Interchange of Hydrogen between Water and Ethylene and between Water and Benzene. Nature 1934, 134, 377−378. (11) Mitsui, T.; Rose, M. K.; Fomin, E.; Ogletree, D. F.; Salmeron, M. Dissociative Hydrogen Adsorption on Palladium Requires Aggregates of Three or More Vacancies. Nature 2003, 422, 705−707.

4. CONCLUSIONS First-principles calculations were used to understand the effects of alloying Pd with Sn on the binding energies of various C2Hx intermediates. Three models with different PdSn mole ratios are constructed to represent PdSn alloy surfaces. The addition of Sn to Pd weakens the binding energies of all C2 intermediates studied herein, which is due to both electronic and geometric effects. Through the electronic structure calculation, it has been found that the d-band center of Pd downshifts with alloying of Pd with Sn. According to the geometric effects, Pd atoms were separated by Sn, which forces the adsorbates to interact with less Pd atoms. Therefore, the binding energies become progressively weaker with the increase of Sn composition. By investigating the coverage effect on adsorption energy, we found that the adsorption energy decreases with the increase of coverage. While analyzing the origin of coverage effect, geometric effect is found to be minimum. The electronic effect is mainly caused by the downshift of the d-band center with the increase of coverage. On the basis of the function of adsorption energy of reactants, we further obtained the surface coverages and reaction barriers of the hydrogenation reaction calculated in high coverage of reactants, which was found to be lower than that in low coverage. Moreover, the rate-limiting step becomes more feasible on Pd3Sn(111) than Pd(111) taking the coverage effect into consideration, which is opposite to the results in which the coverage effect was not considered. Our calculations suggest that a high-coverage model is essential to obtain a more accurate description of the reaction kinetics, when strong binding species exist over catalyst surfaces. 6011

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

Article

The Journal of Physical Chemistry C (12) García-Melchor, M.; Bellarosa, L.; López, N. Unique reaction path in heterogeneous catalysis: the concerted semi-hydrogenation of propyne to propene on CeO2. ACS Catal. 2014, 4, 4015−4020. (13) Yang, B.; Burch, R.; Hardacre, C.; Hu, P.; Hughes, P. Mechanistic Study of 1,3-Butadiene Formation in Acetylene Hydrogenation over the Pd-Based Catalysts Using Density Functional Calculations. J. Phys. Chem. C 2014, 118, 1560−1567. (14) Zhou, H.; Yang, X.; Li, L.; Liu, X.; Huang, Y.; Pan, X.; Wang, A.; Li, J.; Zhang, T. Pdzn Intermetallic Nanostructure with Pd−Zn−Pd Ensembles for Highly Active and Chemoselective Semi-Hydrogenation of Acetylene. ACS Catal. 2016, 6, 1054−1061. (15) Zhang, J.; Sui, Z.; Zhu, Y.-A.; Chen, D.; Zhou, X.; Yuan, W. Composition of the Green Oil in Hydrogenation of Acetylene over a Commercial Pd-Ag/Al2O3catalyst. Chem. Eng. Technol. 2016, 39, 865− 873. (16) Spanjers, C. S.; Held, J. T.; Jones, M. J.; Stanley, D. D.; Sim, R. S.; Janik, M. J.; Rioux, R. M. Zinc Inclusion to Heterogeneous Nickel Catalysts Reduces Oligomerization During the Semi-Hydrogenation of Acetylene. J. Catal. 2014, 316, 164−173. (17) Yang, B.; Burch, R.; Hardacre, C.; Hu, P.; Hughes, P. Selective Hydrogenation of Acetylene over Pd−Boron Catalysts: A Density Functional Theory Study. J. Phys. Chem. C 2014, 118, 3664−3671. (18) Vilé, G.; Albani, D.; Almora-Barrios, N.; López, N.; PérezRamírez, J. Advances in the design of nanostructured catalysts for selective hydrogenation. ChemCatChem 2016, 8, 21−33. (19) Kyriakou, G.; Boucher, M. B.; Jewell, A. D.; Lewis, E. A.; Lawton, T. J.; Baber, A. E.; Tierney, H. L.; Flytzani-Stephanopoulos, M.; Sykes, E. C. H. Isolated metal atom geometries as a strategy for selective heterogeneous hydrogenations. Science 2012, 335, 1209− 1212. (20) López, N.; Vargas-Fuentes, C. Promoters in the Hydrogenation of Alkynes in Mixtures: Insights from Density Functional Theory. Chem. Commun. 2012, 48, 1379−1391. (21) Mei, D.; Neurock, M.; Smith, C. M. Hydrogenation of Acetylene−Ethylene Mixtures over Pd and Pd−Ag Alloys: FirstPrinciples-Based Kinetic Monte Carlo Simulations. J. Catal. 2009, 268, 181−195. (22) Leviness, S.; Nair, V.; Weiss, A. H.; Schay, Z.; Guczi, L. Acetylene Hydrogenation Selectivity Control on PdCu/Al2 O 3 Catalysts. J. Mol. Catal. 1984, 25, 131−140. (23) Armbrüster, M.; Wowsnick, G.; Friedrich, M.; Heggen, M.; Cardoso-Gil, R. Synthesis and Catalytic Properties of Nanoparticulate Intermetallic Ga-Pd Compounds. J. Am. Chem. Soc. 2011, 133, 9112− 9118. (24) Trimm, D. L.; Liu, I. O. Y.; Cant, N. W. The Selective Hydrogenation of Acetylene over a Ni/SiO2 Catalyst in the Presence and Absence of Carbon Monoxide. Appl. Catal., A 2010, 374, 58−64. (25) Norskov, J. K.; Abild-Pedersen, F.; Studt, F.; Bligaard, T. Density Functional Theory in Surface Chemistry and Catalysis. Proc. Natl. Acad. Sci. U.S.A. 2011, 108, 937−943. (26) Nykänen, L.; Honkala, K. Density Functional Theory Study on Propane and Propene Adsorption on Pt(111) and PtSn Alloy Surfaces. J. Phys. Chem. C 2011, 115, 9578−9586. (27) Xin, H.; Linic, S. Communications: Exceptions to the d-Band Model of Chemisorption on Metal Surfaces: The Dominant Role of Repulsion between Adsorbate States and Metal D-States. J. Chem. Phys. 2010, 132, 221101. (28) Studt, F.; Abild-Pedersen, F.; Bligaard, T.; Sørensen, R. Z.; Christensen, C. H.; Nørskov, J. K. Identification of Non-Precious Metal Alloy Catalysts for Selective Hydrogenation of Acetylene. Science 2008, 320, 1320−1322. (29) Studt, F.; Abild-Pedersen, F.; Bligaard, T.; Sørensen, R. Z.; Christensen, C. H.; Nørskov, J. K. On the Role of Surface Modifications of Palladium Catalysts in the Selective Hydrogenation of Acetylene. Angew. Chem. 2008, 120, 9439−9442. (30) Zhao, Z.-J.; Chiu, C.-c.; Gong, J. Molecular Understandings on the Activation of Light Hydrocarbons over Heterogeneous Catalysts. Chem. Sci. 2015, 6, 4403−4425.

(31) Nørskov, J. K.; Bligaard, T.; Rossmeisl, J.; Christensen, C. H. Towards the computational design of solid catalysts. Nat. Chem. 2009, 1, 37−46. (32) Corma, A. Heterogeneous catalysis: understanding for designing, and designing for applications. Angew. Chem., Int. Ed. 2016, 55, 6112−6113. (33) Hammer, B.; Nørskov, J. K. Theoretical Surface Science and CatalysisCalculations and Concepts. Adv. Catal. 2000, 45, 71−129. (34) Bligaard, T.; Nørskov, J. K.; Dahl, S.; Matthiesen, J.; Christensen, C. H.; Sehested, J. The Brønsted−Evans−Polanyi Relation and the Volcano Curve in Heterogeneous Catalysis. J. Catal. 2004, 224, 206−217. (35) Nørskov, J. K.; Bligaard, T.; Logadottir, A.; Bahn, S.; Hansen, L. B.; Bollinger, M.; Bengaard, H.; Hammer, B.; Sljivancanin, Z.; Mavrikakis, M.; et al. Universality in Heterogeneous Catalysis. J. Catal. 2002, 209, 275−278. (36) Yang, B.; Burch, R.; Hardacre, C.; Headdock, G.; Hu, P. Understanding the Optimal Adsorption Energies for Catalyst Screening in Heterogeneous Catalysis. ACS Catal. 2013, 4, 182−186. (37) Che, M. Nobel Prize in Chemistry 1912 to Sabatier: Organic Chemistry or Catalysis? Catal. Today 2013, 218−219, 162−171. (38) Li, Q.; Song, L.; Pan, L.; Zhuang, X.; Ling, M.; Duan, L. Ensemble and Ligand Effects on the Acetylene Adsorption on Ordered PdXAg1-X/Pd(100) Surface Alloys Investigated by Periodic DFT Study. Phys. Chem. Chem. Phys. 2013, 15, 20345−20353. (39) Benedetti, A.; Fagherazzi, G.; Pinna, F.; Rampazzo, G.; Selva, M.; Strukul, G. The Influence of a Second Metal Component (Cu, Sn, Fe) on Pd/SiO2 Activity in the Hydrogenation of 2,4-Dinitrotoluene. Catal. Lett. 1991, 10, 215−223. (40) Armbrüster, M.; Kovnir, K.; Behrens, M.; Teschner, D.; Grin, Y.; Schlögl, R. Pd−Ga Intermetallic Compounds as Highly Selective Semihydrogenation Catalysts. J. Am. Chem. Soc. 2010, 132, 14745− 14747. (41) Wu, Z.; Wegener, E. C.; Tseng, H.-T.; Gallagher, J. R.; Harris, J. W.; Diaz, R. E.; Ren, Y.; Ribeiro, F. H.; Miller, J. T. Pd−In Intermetallic Alloy Nanoparticles: Highly Selective Ethane Dehydrogenation Catalysts. Catal. Sci. Technol. 2016, 6, 6965−6976. (42) Yang, M.-L.; Zhu, Y.-A.; Zhou, X.-G.; Sui, Z.-J.; Chen, D. FirstPrinciples Calculations of Propane Dehydrogenation over PtSn Catalysts. ACS Catal. 2012, 2, 1247−1258. (43) Huang, W.; McCormick, J.; Lobo, R.; Chen, J. Selective Hydrogenation of Acetylene in the Presence of Ethylene on ZeoliteSupported Bimetallic Catalysts. J. Catal. 2007, 246, 40−51. (44) Xu, Z.-F.; Wang, Y. Effects of Alloyed Metal on the Catalysis Activity of Pt for Ethanol Partial Oxidation: Adsorption and Dehydrogenation on Pt3M (M = Pt, Ru, Sn, Re, Rh, and Pd). J. Phys. Chem. C 2011, 115, 20565−20571. (45) Esmaeili, E.; Mortazavi, Y.; Khodadadi, A. A.; Rashidi, A. M.; Rashidzadeh, M. The Role of Tin-Promoted Pd/MWNTs Via the Management of Carbonaceous Species in Selective Hydrogenation of High Concentration Acetylene. Appl. Surf. Sci. 2012, 263, 513−522. (46) Esmaeili, E.; Rashidi, A. M.; Khodadadi, A. A.; Mortazavi, Y.; Rashidzadeh, M. Palladium−Tin Nanocatalysts in High Concentration Acetylene Hydrogenation: A Novel Deactivation Mechanism. Fuel Process. Technol. 2014, 120, 113−122. (47) Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15−50. (48) Mortensen, J. J.; Kaasbjerg, K.; Frederiksen, S. L.; Nørskov, J. K.; Sethna, J. P.; Jacobsen, K. W. Bayesian error estimation in densityfunctional theory. Phys. Rev. Lett. 2005, 95, 216401. (49) Wellendorff, J.; Lundgaard, K. T.; Møgelhøj, A.; Petzold, V.; Landis, D. D.; Nørskov, J. K.; Bligaard, T.; Jacobsen, K. W. Density functionals for surface science: Exchange-correlation model development with Bayesian error estimation. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, 235149. (50) Blöchl, P. E.; Jepsen, O.; Andersen, O. K. Improved tetrahedron method for Brillouin-zone integrations. Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 49, 16223−16233. 6012

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013

Article

The Journal of Physical Chemistry C (51) Methfessel, M.; Paxton, A. T. High-Precision Sampling for Brillouin-Zone Integration in Metals. Phys. Rev. B: Condens. Matter Mater. Phys. 1989, 40, 3616−3621. (52) Zhao, Z.-J.; Kulkarni, A.; Vilella, L.; Nørskov, J. K.; Studt, F. Theoretical Insights into the Selective Oxidation of Methane to Methanol in Copper-Exchanged Mordenite. ACS Catal. 2016, 6, 3760−3766. (53) Campbell, C. T.; Á rnadóttir, L.; Sellers, J. R. V. Kinetic Prefactors of Reactions on Solid Surfaces. Z. Phys. Chem. 2013, 227, 1435−1454. (54) Bruno, T. J.; Haynes, W. M.; Lide, D. R. CRC Handbook of Chemistry and Physics, 95th ed.; CRC Press: Boca Raton, FL, 2014. (55) Lausche, A. C.; Medford, A. J.; Khan, T. S.; Xu, Y.; Bligaard, T.; Abild-Pedersen, F.; Nørskov, J. K.; Studt, F. On the Effect of Coverage-Dependent Adsorbate−Adsorbate Interactions for Co Methanation on Transition Metal Surfaces. J. Catal. 2013, 307, 275−282. (56) Dunphy, J. C.; Rose, M.; Behler, S.; Ogletree, D. F.; Salmeron, M.; Sautet, P. Acetylene Structure and Dynamics on Pd (111). Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 57, R12705. (57) Pallassana, V.; Neurock, M.; Lusvardi, V. S.; Lerou, J. J.; Kragten, D. D.; van Santen, R. A. A Density Functional Theory Analysis of the Reaction Pathways and Intermediates for Ethylene Dehydrogenation over Pd(111). J. Phys. Chem. B 2002, 106, 1656−1669. (58) Zhao, Z.-J.; Moskaleva, L. V.; Aleksandrov, H. A.; Basaran, D.; Rösch, N. Ethylidyne Formation from Ethylene over Pt(111): A Mechanistic Study from First-Principle Calculations. J. Phys. Chem. C 2010, 114, 12190−12201. (59) Ormerod, R. M.; Lambert, R. M.; Bennett, D. W.; Tysoe, W. T. Temperature Programmed Desorption of Co-Adsorbed Hydrogen and Acetylene on Pd (111). Surf. Sci. 1995, 330, 1−10. (60) Stacchiola, D.; Calaza, F.; Zheng, T.; Tysoe, W. T. Hydrocarbon Conversion on Palladium Catalysts. J. Mol. Catal. A: Chem. 2005, 228, 35−45. (61) Yu, W.; Porosoff, M. D.; Chen, J. G. Review of Pt-Based Bimetallic Catalysis: From Model Surfaces to Supported Catalysts. Chem. Rev. 2012, 112, 5780−5817. (62) Podkolzin, S. G.; Alcala, R.; de Pablo, J. J.; Dumesic, J. A. Monte Carlo Simulations of Reaction Kinetics for Ethane Hydrogenolysis over Pt. J. Phys. Chem. B 2002, 106, 9604−9612. (63) Wang, T.; Jiang, F.; Liu, G.; Zeng, L.; Zhao, Z.-j.; Gong, J. Effects of Ga Doping on Pt/CeO2-Al2O3 Catalysts for Propane Dehydrogenation. AIChE J. 2016, 62, 4365−4376. (64) Yang, B.; Burch, R.; Hardacre, C.; Headdock, G.; Hu, P. Influence of Surface Structures, Subsurface Carbon and Hydrogen, and Surface Alloying on the Activity and Selectivity of Acetylene Hydrogenation on Pd Surfaces: A Density Functional Theory Study. J. Catal. 2013, 305, 264−276. (65) Ma, L.-L.; Lv, C.-Q.; Wang, G.-C. A DFT Study and MicroKinetic Analysis of Acetylene Selective Hydrogenation on Pd-Doped Cu(111) Surfaces. Appl. Surf. Sci. 2017, 410, 154−165.

6013

DOI: 10.1021/acs.jpcc.7b11394 J. Phys. Chem. C 2018, 122, 6005−6013