Insight on the Interaction of Methanol-Selective ... - ACS Publications

Jul 3, 2016 - ... Polish Academy of Sciences, 44/52 Kasprzaka Street, 01-224 Warsaw, ... in the selective partial oxidation of methanol to methyl form...
0 downloads 0 Views 3MB Size
Subscriber access provided by University of Sussex Library

Article

Inside the Interaction of Methanol Selective Oxidation Intermediates with Au- or/and Pd-containing Monometallic and Bimetallic Core@Shell Catalysts Kamil Czelej, Karol Cwieka, Juan Carlos Colmenares, and Krzysztof J. Kurzydlowski Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.6b01906 • Publication Date (Web): 03 Jul 2016 Downloaded from http://pubs.acs.org on July 5, 2016

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.

Langmuir 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 36

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

Langmuir

Inside the Interaction of Methanol Selective Oxidation Intermediates with Au- or/and Pdcontaining Monometallic and Bimetallic Core@Shell Catalysts Kamil Czelej*,†, Karol Cwieka†, Juan Carlos Colmenares‡, Krzysztof J. Kurzydlowski† †

Faculty of Materials Science and Engineering, Warsaw University of Technology, 141 Woloska Str., 02-507 Warsaw, Poland



Institute of Physical Chemistry, Polish Academy of Science, 44/52 Kasprzaka Str., 01-224 Warsaw, Poland

1 ACS Paragon Plus Environment

Langmuir

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 36

ABSTRACT: Using density functional theory (DFT), the interaction of crucial molecules involved in the selective partial oxidation of methanol to methyl formate (MF), with monometallic Au, Pd and bimetallic Au/Pd, Pd/Au core@shell catalysts is systematically investigated. The core@shell structures modeled in this study consist of Au(111), Pd(111) cores covered by monolayer of Pd, Au, respectively. Our results indicate that the adsorption strength of the molecules examined as a function of catalytic surface decreases in the order of Au/Pd(111) > Pd(111) > Au(111) > Pd/Au(111) and correlate well with the d-band center model. The preadsorbtion of oxygen is found to have a positive impact on the selective partial oxidation reaction due to the stabilization of CH3OH and HCHO on the catalyst surface and the simultaneous intensification of MF desorption. Based on a dynamical matrix approach combined with statistical thermodynamic we propose a simple route for evaluating the Gibbs free energy of adsorption as a function of temperature. This method allows us to anticipate the relative temperature stability of molecules involved in the selective partial oxidation of methanol to MF in terms of catalytic surface.

2 ACS Paragon Plus Environment

Page 3 of 36

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

Langmuir

INTRODUCTION Selective partial oxidation of methanol to methyl formate (MF) over heterogeneous catalysts is an attractive and environmentally benign route for methanol conversion in industrial chemistry. The efficient production of methanol downstream products remains a challenging problem of technological importance; therefore, the design of highly efficient heterogeneous catalysts with enhanced selectivity and activity for partial oxidation of methanol is of paramount significance in academia. Currently, noble metals and their alloys13

, transition metal oxides in the form of nanoparticles4, 5, or metal-supported TiO26 are

commonly used to conduct methanol transformation to MF. However, severe limitations such as low yield and the production of a large amount of undesirable byproducts, mainly CO2 have been encountered7. To overcome these limitations worldwide attention has been focused on developing novel, highly efficient core@shell catalysts with tunable chemical composition8-18. Bimetallic nanoparticles exhibit distinctively different physico-chemical properties in comparison to their monometallic counterparts. Specifically, gold-palladium (Au-Pd) catalysts have displayed significantly increased activity in a number of chemical reactions, such as the oxidation of CO19, 20, the acetoxylation of ethylene to vinyl acetate21, 22, the direct synthesis of hydrogen peroxide from H2 and O223,

24

, the epoxidation of alkene25, and the

selective oxidation of alcohol to aldehyde26-28. The morphology of bimetallic nanoparticles depends on the synthesis technique, post-synthesis heat treatment, and the solubility of both metals. Due to unlimited solubility in Au-Pd system the atomic distribution within nanoparticles may range from homogeneous alloys to so-called core@shell structures. Commonly used synthesis techniques such as under potential deposition (UPD) enables the fabrication of a uniform shell on the core (from monolayer up to multilayer in thickness) via careful potential control29-33. Amongst a variety of Au/Pd core@shell structures Silva et al.10 3 ACS Paragon Plus Environment

Langmuir

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 36

have recently shown that Au cores with a monolayer of Pd exhibit the highest catalytic activity in selective oxidation of benzyl alcohol. By means of combined experimental and theoretical approaches they observed the volcano-like behavior in terms of a variable Au:Pd molar ratio. An increased catalytic activity was attributed to the balance between the number of active sites and the ease of product desorption. The superiority of core@monolayer-shell catalysts over uniform catalysts was also demonstrated by Zhang and Henkelman in their recent study11. Based solely on adsorbates binding energies computed within density functional theory and the linear correlation between the binding of adsorbates to the shell vs alloy-core composition, they proposed a simple approach to the design of a catalyst for a particular reaction. In the current study, we analyzed the crucial molecules involved in the selective partial oxidation of methanol to MF, interacting with Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) surfaces, by means of ab initio calculations. Interesting trends in adsorption energies for various molecules as a function of catalytic surface have been found, and subsequently discussed in term of the electronic properties of the surfaces. Further, our study demonstrates the impact of preadsorbed oxygen on the binding of adsorbates to the surfaces. Finally, using the dynamical matrix approach combined with statistical thermodynamic, we developed a simple route to evaluate the Gibbs free energy of adsorption, and hence the temperature stability of adsorbate species.

THEORETICAL METHODS In the slab-supercell approximation, the adsorption energy of molecule on a solid surface can be simply expressed as follows:  =   +   +  

(1)

The first term: 4 ACS Paragon Plus Environment

Page 5 of 36

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

Langmuir

  =  / −   +  

(2)

denotes the electronic contribution to the adsorption energy, computed as a difference between the total energy of slab with adsorbed molecule and the total energy of bare slab plus the total energy of isolated molecule in a gas phase. The second term is the dispersion energy correction, which is computed using a pairwise potential defined as: 

'   = − ∑' ) ∑( ∑&

where



-.)(

represents

the

 )(,&   *+ , !"  ,$ %

interatomic

(3)

(dipole-dipole)

dispersion

coefficient,

, )(,& = |, ), − , (,& |, where *+ , )(,&  is a damping function, used to cancel spurious interaction at short distances within the overlapping electron densities regime. The summation is over all atoms 0 and all translations of the unit cell 1 = 2 , 2 , 23 , under assumption that

4 ≠ 6 for 1 = 0. A detailed explanation of the dispersion energy correction and its evaluation can be found elsewhere34. The last term of the equation (1) is the zero-point energy, which

can be calculated in the harmonic approximation: 

 =  ∑; =  + ?B -@ CD

(5)

where the heat capacity in the harmonic approximation can be expanded in series: 5 ACS Paragon Plus Environment

Langmuir

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

∗

?B -@ CD = ∑; − DP

(10)

and hence the Gibbs free energy of adsorption: ∆h = h − h:

j

(11)

j denotes the Gibbs free energy of an ideal gas, which is reasonable estimation where h:

within low pressure regime.

COMPUTATIONAL DETAILS

6 ACS Paragon Plus Environment

Page 7 of 36

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

Langmuir

Electronic structure calculations were performed in the framework of density functional theory (DFT) as implemented in the Vienna Ab initio Simulation Package (VASP)36. The exchange-correlation energy was computed using Perdew-Burke-Ernzerhof (PBE) general gradient approximation (GGA) functional37. The projector augmented wave (PAW) method38 was employed to describe valence-core interaction. The valence electron configuration for each element was: (1s1) H, (2s22p2) C, (2s22p4) O, (5d106s1) Au and (4d95s1) Pd. An impact of nonlocal long-range interactions between adsorbed molecules and metal surfaces on the adsorption energy was verified a posteriori according to the vdW correction scheme proposed by Tkatchenko and Scheffler (DFT-TS)39. Dipole correction was considered for all calculations40. In order to simulate individual CH3OH, CH3O, CH2OH, CH3OOCH2, HCHO and CH3COOH molecules on Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) surfaces, p(3 x 3) slab geometries, corresponding to a 1/9 ML coverage, containing 4 metal layers were constructed from the optimized bulk unit cell of gold and palladium. In this work Au/Pd(111) and Pd/Au(111) denote an Au(111) surface with a monolayer of Pd shell and a Pd(111) surface with a monolayer of Au shell, respectively. Our model represents the most stable (111) facets of large core@shell nanoparticles41, 42. The slab approach was widely used to model physicochemical properties of nanoparticles, such as binding energy of adsorbates, surface reactivity as a function of chemical composition, concentration and geometry43-47. In comparison to cluster approach which is computationally demanding and tractable only for relatively small nanoparticles, the slab representation gives the possibility of studying variety of structures in a systematic way. In addition, there are few studies published that compare adsorption energies calculated using cluster and slab approaches. For instance, the adsorption energy of NO on Pd79 nanoclusters turned out to be 9% higher than on Pd(111) slab48. Another investigations report the adsorption energies of C atoms on Pd79 nanocluster to be 7 ACS Paragon Plus Environment

Langmuir

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 8 of 36

only 4% different than on Pd(111) slab49, 50. Since the typical size of metallic nanoparticles is tens of nanometers the difference in the adsorption energy ought to be even smaller; therefore, the use of slab approach is justified. A vacuum level of ~15 Å was large enough to screen the spurious interaction between periodic images of the surface. Convergence of the total energy for the slab-molecule system was assured by setting a 450 eV cutoff energy, and the Brillouin zone sampling in a 5x5x1 Monkhorst-Pack k-point grid51. For geometry optimization, the top 2 layers of metal atoms together with the molecule were allowed to relax towards a minimum of the total energy, while the bottom 2 layers were kept frozen in their lattice positions. After all forces acting on atoms dropped below 0.01 eV/Å and energy convergence criterion of 1e5 eV was reached, the relaxations were terminated. Vibrational modes and frequencies were provided by diagonalization of the Hessian matrix using the dynamical matrix method. To quantify charge transfer from the (111) surfaces to the molecules, a Bader analysis was carried out for all conformations using the code developed by Henkelman et al.52. Gas-phase structures of the molecules and their vibrational properties were calculated in a simulation box (20x20x20 Å) applying the gamma k-point sampling. The procedure presented in Section 2 was implemented in Fortran 95 environment and subsequently used to determine Gibbs free energies of adsorption.

RESULTS AND DISCUSSION Analysis of Surface Species The adsorption geometries and energies of crucial intermediates involved in the selective partial oxidation of methanol to MF on Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) are investigated in this section. The most stable conformations are displayed in Figure 1 and the adsorption energies together with geometrical parameters are listed in the Supporting Information (see Table S1). Apart from alkoxy hemiacetal (CH3OOCH2), the most stable 8 ACS Paragon Plus Environment

Page 9 of 36

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

Langmuir

adsorption site for a given molecule and the corresponding molecular orientation are roughly the same for each surface. The bond lengths and adsorption energies, however, depend strongly upon the chemical composition of the surface.

Figure 1. Optimized geometries of the most stable adsorption species involved in the selective partial oxidation of methanol to MF on Au(111), Pd(111) and core@shell Au/Pd(111), Pd/Au(111) catalysts. Table 1. Bader charges for possible intermediates of methanol selective partial oxidation to MF adsorbed on the Au(111), Pd(111), Au/Pd(111) and Pd/Au(111). The values in parentheses denote oxygen preadsorbed surfaces.

Molecule CH3OH CH2OH CH3O HCHO CH3OOCH2 CH3COOH

Pd/Au(111) 0.014 (0.104) 0.144 (0.288) -0.423 (-0.351) 0.003 (0.043) -0.411 (-0.379) -0.005 (0.011)

Bader charge of the moiety on the Au(111) Pt(111) 0.026 (0.108) 0.050 (0.081) 0.135 (0.282) 0.099 (0.174) -0.429 (-0.410) -0.372 (-0.359) 0.013 (0.044) 0.021 (0.021) -0.438 (-0.409) -0.390 (-0.368) -0.005 (-0.003) -0.006 (-0.002)

Au/Pd(111) 0.040 (0.079) 0.082 (0.170) -0.398 (-0.383) 0.010 (0.031) -0.413 (-0.418) -0.001 (-0.003)

9 ACS Paragon Plus Environment

Langmuir

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 36

CH3OH. Relaxation of the gas-phase methanol leads to equilibrium bond lengths equal to 1.10 Å for C–H, 1.43 Å for O–C, and 0.97 Å for O–H, which is in good agreement with the experimental values of 1.09 Å, 1.43 Å and 0.95 Å, respectively53. According to previous investigations54-56, the CH3OH adsorption on metallic surface is realized by donation of the lone-pair electrons from oxygen to the surface atoms, forming a weakly bound moiety, unstable upon increasing temperature of system. Slightly positive Bader charges of CH3OH molecule adsorbed on the examined surfaces, presented in Table 1, support the donation nature of adsorption. Our calculations reveal that methanol is bound to the surface top site through its oxygen atom, with a decreasing M-O distance of 2.95 Å, 2.77 Å, 2.37 Å and 2.35 Å for Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), respectively. The corresponding adsorption energies of -0.30 eV, -0.36 eV, -0.48 eV, and -0.53 eV follow the same trend (see Figure 2 and Supporting Information) and agree well with available experimental data of 0.42 eV for Au(111)57 and -0.38÷-058 eV for Pd(111)58. Previous theoretical works report the adsorption energy of CH3OH on Au(111) to be -0.15 eV59 and -0.17 eV3. The significant underestimation may be explained by negligence of long-range van der Waals interactions, which have a remarkable, over 60%, contribution to the binding energy.

CH2OH. A hydroxymethyl radical is the intermediate formed by the β-H elimination of the adsorbed methanol. It has a nonplanar structure with the most stable conformation at a top site, bound to the surface through the carbon atom. Due to the open-shell electronic configuration of CH2OH, its adsorption energy is relatively high, -1.66 eV, -1.81 eV, 2.40 eV, and -2.42 eV for Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), respectively. The corresponding M-C distances are computed to be 2.19 Å, 2.17 Å, 2.06 Å and 2.04 Å. Positive Bader charges, quantitatively larger for CH2OH than CH3OH, indicate that the delocalized

10 ACS Paragon Plus Environment

Page 11 of 36

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

Langmuir

radical electron in a π*orbital can be effectively transferred from the molecule to the metal atoms.

CH3O. A methoxy radical is formed by the α-H elimination of the adsorbed methanol molecule. The CH3O prefers to adsorb at the three-fold fcc site through its oxygen atom, with the C–O axis nearly perpendicular to the surface. The chemical bond formed between molecule and surface incorporates 3 metal atoms adjacent to the hollow site. As opposed to CH2OH, methoxy acts as an electron acceptor (see Table 1), with the computed Bader charges of ~ -0.4 eV. The adsorption energies of -1.29 eV, -1.86 eV, -2.35 eV, and -2.57 eV for Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), respectively and small van der Waals energy contributions indicate rather strong chemisorption.

HCHO. There are two possibilities of formaldehyde (HCHO) formation during the selective partial oxidation of methanol to MF: first is the α-H elimination of the adsorbed CH2OH or second is the β-H elimination of the adsorbed CH3O. Our results reveal a weak interaction between HCHO and the examined surfaces, consistently with previous theoretical predictions1, 3. In case of Au-terminated surfaces the adsorption has predominantly physical nature, whereas for Pd-terminated surfaces the contribution of dispersive and electronic energy is comparable. Formaldehyde preferentially binds to the bridge site through oxygen and hydrogen atoms, and its molecular plane is perpendicular to the surface. The calculated adsorption energies are -0.23 eV, -0.27 eV, -0.36 eV, and -0.42 eV for Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), respectively. The distance between molecule and the nearest surface atom decreases significantly from ~3.00 Å (Au-termination) to ~2.27 Å (Pd-termination).

11 ACS Paragon Plus Environment

Langmuir

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 36

CH3OOCH2. An alkoxy hemiacetal intermediate is formed as a result of the coupling reaction between electron-rich oxygen in methoxy and electron-deficient carbon in formaldehyde. Two distinguishable adsorption sites, dependent on the surface termination have been found. For Au-terminated surfaces, CH3OOCH2 binds to the bridge site through two of its oxygen atoms, whereas for Pd-terminated surfaces, the most stable adsorption site is the three-fold fcc hollow with only one oxygen involved. Negative Bader charges indicate a charge transfer from surface to the molecule. The calculated adsorption energies of 1.56 eV, -1.75 eV, -2.31 eV, and -2.53 eV for Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), respectively suggest a strong chemisorption.

CH3COOH. The final product, methyl formate, is created by abstraction of a hydrogen atom from CH3OOCH2. As can be seen in Figure 2 the interaction between MF and the examined surfaces is very weak and almost entirely governed by the van der Waals forces. The most stable conformation was found at the bridge site through oxygen and hydrogen atoms, with large (3.29–3.90 Å) distance between molecule and the nearest surface atom. This conformation affords an adsorption energies of -0.24 eV, -0.27 eV, -0.29 eV, and -0.36 eV for Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), respectively.

12 ACS Paragon Plus Environment

Page 13 of 36

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

Langmuir

Figure 2. Adsorption energies of molecules involved in the selective partial oxidation of methanol to MF on clean and preadsorbed oxygen Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) surfaces. Impact of vdW and ZPE Corrections on the Adsorption Energy In order to predict the adsorption energy of molecules on solid surfaces at 0 K more accurately, one has to take into account both, the nonlocal van der Waals interactions and the zero point energy contribution. As shown in our previous work60 the incorporation of van der Waals energy may even determine stability of molecular adsorption by changing its character from endothermic to exothermic. The electronic, vdW and ZPE contributions to the adsorption energy of crucial intermediates involved in the selective partial oxidation of methanol to MF on Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) are juxtaposed together in Figure 2. In the case of weakly bound species (CH3OH, HCHO, CH3COOH) the conventional 13 ACS Paragon Plus Environment

Langmuir

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 36

DFT values of the adsorption energy are significantly underestimated; therefore, the vdW+ZPE term has to be taken into account to reflect experimental measurements. By contrast, strongly bound intermediates (CH2OH, CH3O, CH3OOCH2) are only slightly affected by the vdW+ZPE contribution. Their adsorption energy can be predicted fairy well using the conventional DFT approach exclusively.

Binding Energy Trends vs Intrinsic Properties of Surfaces Trends in the binding energy of six crucial molecules involved in the selective partial oxidation of methanol to MF were disclosed as a function of surface chemical composition. As can be seen in Figure 2, for a given molecule the adsorption energy with respect to the catalyst surface decreases as follows: Au/Pd(111) > Pd(111) > Au(111) > Pd/Au(111). In order to clarify relative binding properties of the surfaces examined, the d-projected density of states (d-DOS) together with the d-band center were depicted in Figure 3. The presented trends correlate well with the increasing d-band center of the top surface layer, defined as:

k =

? l ? l

(12)

where m is the d-DOS. For the core@shell catalysts there are two major factors affecting the d-band center: strain induced within the shell, derived from the core-shell mismatch, and charge redistribution11. The higher the difference in atomic radius of core and shell atoms the larger mismatch, and therefore the impact of strain on the d-band center shift is more pronounced. In turn, charge redistribution is a function of electronegativity difference between elements in the core and shell. The variation of Bader charges computed for the surface atoms (CX/SX vs CY/SX where X and Y are elements and C and S denotes core and shell, respectively) are negligibly small (0.01 - 0.04 e). The small charge transfer results from similar electronegativity of Pd and Au. Based on this observation we anticipate that the d-

14 ACS Paragon Plus Environment

Page 15 of 36

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

Langmuir

band center shift is dominated by the strain effect. Similar conclusions were drawn for PdAu@Pt catalyst by Zhang and Henkelman11.

Figure 3. The d-projected density of states (d-DOS) together with the d-band center calculated for the top layer of Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) bare surfaces. Total charges per atom of first and second layer of each slab are shown on the right side.

Effect of Preadsorbed Oxygen on Adsorption It is commonly known that preadsorbed oxygen on noble metal surfaces significantly influences both the reaction barriers and selectivity in catalytic oxidation of alcohols1, 3, 61. For this reason, an impact of preadsorbed oxygen on the adsorption of molecules involved in the selective partial oxidation of methanol to MF has to be highlighted. According to previous 15 ACS Paragon Plus Environment

Langmuir

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 36

theoretical data, atomic oxygen prefers binding at the fcc sites with very high adsorption energy. Our calculations provide oxygen adsorption energies of -2.74 eV, -3.47 eV, -4.68 eV and -4.73 eV for Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), respectively, confirming very strong, exothermic chemisorption. Adsorption energies of CH3OH, CH3O, CH2OH, CH3OOCH2, HCHO and CH3COOH on the oxygen precovered Pd/Au(111), Au(111), Pd(111) and Au/Pd(111) surfaces and corresponding geometries are presented in Figure 2 and Figure S1 (Supporting Information). Depends on molecule, the presence of atomic oxygen on catalyst surface may either improve or deteriorate the adsorption strength. A remarkable enhancement of binding energy is observed for the weakly bound CH3OH and HCHO. In the case of methanol, the reported adsorption energy gain is ~0.2 eV for Au-terminated surfaces and ~0.1 eV for Pd-terminated surfaces, whereas in the case of formaldehyde the values are ~0.15 eV and ~0.1 eV, respectively. In both cases, the creation of a hydrogen bond is responsible for additional stabilization effect manifesting itself in energy gain. From the point of view of catalytic performance, the additional preadsorbed oxygen stabilization effect will simultaneously intensify the CH3OH adsorption and hamper the HCHO desorption from the catalyst surface. In particular, the latter phenomenon plays an essential role in selectivity to MF, since the intense desorption of HCHO suppresses its coupling reaction with CH3O, which in turn, leads to the formation of CH3OOCH2 and subsequently MF. Furthermore, the adsorbed oxygen species act as an electron reservoir, and as can be seen in Bader analysis results (see Table 1), promote the charge transfer from moiety to the surface, increasing activation of chemical bonds in the molecule. This may explain the significant reduction of elemental reaction energy barriers, as observed for the selective partial oxidation of CH3OH using oxygen precovered conditions. The strongly chemisorbed species, such as CH2OH, CH3O or CH3OOCH2 are however only slightly affected by the presence of preadsorbed oxygen atoms. Basically, their binding energy is slightly lower in comparison to counterparts 16 ACS Paragon Plus Environment

Page 17 of 36

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

Langmuir

on clean surfaces. A positive effect of preadsorbed oxygen on catalytic performance in the selective partial oxidation of methanol to MF is also visible for CH3COOH. Due to the reduction in adsorption energies the ease of product (MF) desorption should be accordingly improved.

Thermodynamic Stability of Adsorbates In the previous sections we discussed the adsorption of CH3OH, CH3O, CH2OH, CH3OOCH2, HCHO and CH3COOH on Pd/Au(111), Au(111), Pd(111) and Au/Pd(111), taking into account electronic, dispersive and ZPE energy contributions. In principle, we assumed that molecules have in the gas phase 100% of their standard entropy, whereas adsorbates possess no entropy at all. These assumptions are meaningful within a low temperature regime but may fail remarkably as the temperature rises. To provide insight into the thermodynamic stability of adsorbates involved in the selective partial oxidation of methanol to MF, we calculated the Gibbs free energy changes as a function of temperature, consistent with the procedure described in Section 2. The results are presented in Figure 4, with dotted vertical lines indicating the temperature at which desorption occurs. The temperature stability of methanol, formaldehyde, and MF adsorbed on the surfaces examined are in reasonable agreement with available experimental data57, 62-64; nevertheless, there is one problem with the TPD results which needs to be highlighted. To be meaningful, weakly bound CH3OH, HCHO, and CH3COOH have to be detectable using the TPD experiment. In the case of adsorbate high-coverages the TPD peaks are high and sharp; therefore, they can be easily attributed to specific temperatures. However, as the coverage decreases well below 1 ML, the TPD peaks become low and smeared, making it difficult to unambiguously interpret. Since our calculations concern very low (~0.11 ML) coverage, the comparison 17 ACS Paragon Plus Environment

Langmuir

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

between the computed and the TPD data is qualitative rather than quantitative. The predicted temperature stability of strongly chemisorbed species on a given surface decreases as follow: CH3OOCH2 > CH2OH > CH3O. This order was determined under the assumption that only adsorption and desorption processes exist. In reality, however, decomposition of big molecules such as CH3OOCH2 may occur at elevated temperatures, changing the relative stability.

Figure 4. The Gibbs free energies of adsorption as a function of temperature calculated for molecules involved in the selective partial oxidation of methanol to MF on Au(111), Pd(111) and core@shell Au/Pd(111), Pd/Au(111) catalysts. Thermodynamic results for oxygen preadsorbed conditions are presented in Figure 5, in a similar manner. In comparison to the adsorption on bare surfaces, desorption temperature for weakly bound CH3OH and HCHO is shifted to higher values due to the additional stabilization effect via hydrogen bond. As we anticipated in previous paragraph, the preadsorbed oxygen intensify MF desorption by lowering desorption temperature. It is worth 18 ACS Paragon Plus Environment

Page 18 of 36

Page 19 of 36

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

Langmuir

mentioning that relative stability of molecular species as a function of catalytic surface changes with respect to temperature for the oxygen preadsorbed conditions. The main factor causing the changes arises from the entropy contribution, which in turn depends on the ground state geometry of adsorbate. As can be seen in the Supporting Information (see Figure S1), the ground state adsorption geometry for a given molecule can differ remarkably when the atomic oxygen is preadsorbed on the catalytic surface.

Figure 5. The Gibbs free energies of adsorption as a function of temperature calculated for molecules involved in the selective partial oxidation of methanol to MF on oxygen preadsorbed Au(111), Pd(111) and core@shell Au/Pd(111), Pd/Au(111) catalysts.

CONCLUSIONS Systematic DFT calculations have been employed to investigate the interaction of crucial molecules involved in the selective partial oxidation of methanol to MF, with monometallic 19 ACS Paragon Plus Environment

Langmuir

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 36

Au, Pd and bimetallic Au/Pd, Pd/Au core@shell catalysts. We have found that for a given molecule the adsorption energy with respect to the catalyst surface decreases in the order of Au/Pd(111) > Pd(111) > Au(111) > Pd/Au(111) and correlate well with the d-band center of the top surface layer. The oxygen preadsorbed on the surfaces examined can either stabilize weakly bound molecules (CH3OH and HCHO) by forming the hydrogen bonds, slightly decrease the adsorption strength of strongly chemisorbed species (CH3O, CH2OH and CH3OOCH2) or destabilize the weakly adsorbed CH3COOH, promoting its desorption. These results partially explain the huge impact of oxygen species on the selectivity and reaction yield observed previously for gold catalysts. Furthermore, we proposed simple yet robust route to evaluate the Gibbs free energy of adsorption, which may be useful in the interpretation of TPD data in particular for very low surface coverage. Test calculations carried out for the surfaces examined showed reasonable agreement with available experimental TPD results. Further investigation of the reaction barriers and effect of tuning the chemical composition of core and shell will be interesting in terms of maximizing the catalytic activity towards methanol selective oxidation to MF.

AUTHOR INFORMATION Corresponding Author * E-mail address: [email protected] Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

ACKNOWLEDGEMENTS 20 ACS Paragon Plus Environment

Page 21 of 36

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

Langmuir

Computing resources were provided by High Performance Computing facilities of the Interdisciplinary Centre for Mathematical and Computational Modeling (ICM) of the University of Warsaw under Grant No. G61-4 and Poznań Supercomputing and Networking Center (PSNC) under Grant No. 275. ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge via Internet at http://pubs.acs.org. Detailed geometries and adsorption energies Ead of stable adsorption species involved in the selective partial oxidation of methanol to MF on clean and oxygen precovered Au(111), Pd(111), Au/Pd(111) and Pd/Au(111).

REFERENCES (1) Grimme, S. Semiempirical GGA-type Density Functional Constructed With a Longrange Dispersion Correction, Journal of Computational Chemistry 2006, 27, 1787-1799. (2) Zhang, L.; Henkelman, G. Tuning the Oxygen Reduction Activity of Pd Shell Nanoparticles with Random Alloy Cores, Journal of Physical Chemistry C 2012, 116, 20860-20865. (3) Zhang, L.; Iyyamperumal, R.; Yancey, D. F.; Crooks, R. M.; Henkelman, G. Design of Pt-Shell Nanoparticles With Alloy Cores For The Oxygen Reduction Reaction,

ACS Nano 2013, 7, 9168-9172. (4) Zhang, L.; Henkelman, G. Computational Design of Alloy-Core@Shell Metal Nanoparticle Catalysts, ACS Catalysis 2015, 5, 655-660. (5) Zhang, H.; Watanabe, T.; Okumura, M.; Haruta M.; Toshima, N. Catalytically Highly 21 ACS Paragon Plus Environment

Langmuir

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 22 of 36

Active Top Gold Atom On Palladium Nanocluster, Nature Materials 2012, 11, 49-52. (6) Yancey, D. F.; Carino E. V.; Crooks, R. M. Electrochemical Synthesis and Electrocatalytic Properties of Au@Pt Dendrimer-Encapsulated Nanoparticles, Journal of

the American Chemical Society 2010, 132, 10988-10989. (7) Yancey, D. F.; Zhang, L.; Crooks, R. M.; Henkelman, G. Au@Pt Dendrimer Encapsulated Nanoparticles As Model Electrocatalysts For Comparison Of Experiment And Theory, Chemical Science 2012, 3, 1033-1040. (8) Xu, J.; White, T.; Li, P.; He, C.; Yu, J.; Yuan W.; Han, Y.-F. Biphasic Pd−Au Alloy Catalyst for Low-Temperature CO Oxidation, Journal of the American Chemical Society 2010, 132, 10398-10406. (9) Xu, B.; Haubrich, J.; Baker, T. A.; Kaxiras, E.; Friend, C. M. Theoretical Study of OAssisted Selective Coupling of Methanol on Au(111), The Journal of Physical Chemistry

C 2011, 115, 3703-3708. (10) Wang, L.; He, C.; Zhang, W.; Li, Z.; Yang, J. Methanol-Selective Oxidation Pathways on Au Surfaces: A First Princiles Study, The Journal of Physical Chemistry C 2014, 118, 17511-17520. (11) Tkatchenko, A.; Scheffler, M. Accurate Molecular Van Der Waals Interactions from Ground-State Electron Density and Free-Atom Reference Data, Physical Review Letters 2009, 102, 073005. (12) Tang, W.; Henkelman, G. Charge Redistribution In Core-Shell Nanoparticles To Promote Oxygen Reduction, The Journal of Chemical Physics 2009, 130, 194504. (13) Tang, W.; Zhang, L.; Henkelman, G. Catalytic Activity of Pd/Cu Random Alloy Nanoparticles for Oxygen Reduction, The Journal of Physical Chemistry Letters 2011, 2,

22 ACS Paragon Plus Environment

Page 23 of 36

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

Langmuir

1328-1331. (14) Strasser, P.; Koh, S.; Anniyev, T.; Greeley, J.; More, K.; Yu, C.; Liu, Z.; Kaya, S.; Norldund, D.; Ogasawara, H.; Toney, M. F.; Nilsson, A. Lattice-Strain Control Of The Activity In Dealloyed Core–Shell Fuel Cell Catalysts, Nature Chemistry 2010, 2, 454460. (15) Silva, T. A. G.; Teixeira-Neto, E.; López N.; Rossi, L. M. Volcano-like Behavior of AuPd Core-shell Nanoparticles in the Selective Oxidation of Alcohols, Scientific Reports 2014, 4, 5766. (16) Shiizaki, S.; Nagashima, I.; Matsamura, Y.; Haruta, M. A Kinetic Study Of Methanol Decomposition Catalyzed Over Plate-Type Palladium Catalyst, Catalysis Letters 1998, 56, 227-230. (17) Zhang, N.; Liu, S.; Xu, Y.-J. Recent Progress on Metal Core@Semiconductor Shell Nanocomposites as a Promising Type of Photocatalyst, Nanoscale 2012, 4, 2227-2238. (18) Zhang, Y.; Zhang, N.; Tang, Z.-R.; Xu, Y.-J. Graphene Oxide as a Surfactant and Support for In-Situ Synthesis of Au–Pd Nanoalloys with Improved Visible Light Photocatalytic Activity, The Journal of Physical Chemistry C 2014, 118, 5299-5308. (19) Sexton, B. A.; Hughes, A. E. A Comparison Of Weak Molecular Adsorption Of Organic Molecules On Clean Copper And Platinum Surfaces, Surface Science 1984, 140, 227248. (20) Serpell, C. J.; Cookson, J.; Ozkaya, D.; Beer, P. D. Core@Shell Bimetallic Nanoparticle Synthesis Via Anion Coordination, Nature Chemistry 2011, 3, 478-483. (21) Sasaki, K.; Wang, J. X.; Naohara, H.; Marinkovic, N.; More, K.; Inada, H.; Adzic, R. R. Recent Advances In Platinum Monolayer Electrocatalysts For Oxygen Reduction

23 ACS Paragon Plus Environment

Langmuir

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 24 of 36

Reaction: Scale-up Synthesis, Structure And Activity Of Pt Shells On Pd Cores,

Electrochimica Acta 2010, 55, 2645-2652. (22) Peterson, A. A.; Grabow, L. C.; Brennan, T. P.; Shong, B.; Ooi, C.; Wu, D. M.; Li, C. W.; Kushwaha, A.; Medford, A. J.; Mbuga, F.; Li, L.; Nørskov, J. K. Finite-Size Effects in O and CO Adsorption for the Late Transition Metals, Topics in Catalysis 2012, 55, 1276-1282. (23) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple, Physical Review Letters 1996, 77, 3865-3868. (24) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin Zone Integrations, Physical

Review B 1976, 13, 5188-5192. (25) Mavrikakis, M.; Barteau, M. A. Oxygenate Reaction Pathways On Transition Metal Surfaces, Journal of Molecular Catalysis A: Chemical 1998, 131, 135-147. (26) Liu, S.; Jin, P.; Zhang, D.; Hao, C.; Yang, X. Reaction Mechanism for Methanol Oxidation on Au(1 1 1): A Density Functional Theory Study, Applied Surface Science 2013, 265, 443-451. (27) Lide, D. R. CRC handbook of chemistry and physics, CRC Press, 2001. (28) Li, L.; Larsen, A. H.; Romero, N. A.; Morozov, V. A.; Glinsvad, C.; Abild-Petersen, F.; Greeley, J.; Jacobsen, K. W.; Nørskov, J. K. Investigation of Catalytic Finite-SizeEffects of Platinum Metal Clusters, The Journal of Physical Chemistry Letters 2013, 4, 222-226. (29) Landon, P.; Collier, P. J.; Papworth, A. J.; Kiely, C. J.; Hutchings, G. J. Direct Formation Of Hydrogen Peroxide From H2/O2 Using A Gold Catalyst, Chemical Communications 2002, 18, 2058-2059.

24 ACS Paragon Plus Environment

Page 25 of 36

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

Langmuir

(30) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-energy Calculations Using a Plane-wave Basis Set, Physical Review B 1996, 54, 11169-11186. (31) Koenigsmann, C.; Sutter, E.; Adzic, R. R.; Wong, S. S. Size- and CompositionDependent Enhancement of Electrocatalytic Oxygen Reduction Performance in Ultrathin Palladium–Gold (Pd1–xAux) Nanowires, The Journal of Physical Chemistry C 2012, 116, 15297-15306. (32) Iyyamperumal, R.; Zhang, L.; Henkelman, G.; Crooks, R. M. Efficient Electrocatalytic Oxidation of Formic Acid Using Au@Pt Dendrimer-Encapsulated Nanoparticles,

Journal of the American Chemical Society 2013, 135, 5521-5524. (33) Hughes, M. D.; Xu, Y.-J.; Jenkins, P.; McMorn, P.; Landon, P.; Enache, D. I.; Carley, A. F.; Attard, G. A.; Hutchings, G. J.; King, F.; Stitt, E. H.; Johnston, P.; Griffin, K.; Kiely, C. J. Tunable Gold Catalysts For Selective Hydrocarbon Oxidation Under Mild Conditions, Nature 2005, 437, 1132-1135. (34) Henkelman, G.; Arnaldsson, A.; Jonsson, H. A Fast and Robust Algorithm for Bader Decomposition of Charge Density, Computational Materials Science 2006, 36, 354-360. (35) Han, Y.-F.; Wang, J.-H.; Kumar, D.; Yan, Z.; Goodwin, D. A Kinetic Study Of Vinyl Acetate Synthesis Over Pd-Based Catalysts: Kinetics Of Vinyl Acetate Synthesis Over Pd–Au/SiO2 and Pd/SiO2 Catalysts, Journal of Catalysis 2005, 232, 467-475. (36) Han, C.; Yang, X.; Gao, G.; Wang, J.; Lu, H.; Liu, J.; Tong, M.; Liang, X. Selective Oxidation of Methanol to Methyl Formate on Catalysts of Au-Ag Alloy Nanoparticles Supported on Titania Under UV Irradiation, Green Chemistry 2014, 16, 3603-3615. (37) Gong, J.; Flaherty, D. W.; Ojifinni, R. A.; White, J. M.; Mullins, C. B. Surface Chemistry of Methanol on Clean and Atomic Oxygen Pre-Covered Au(111), The Journal

of Physical Chemistry C 2008, 112, 5501-5509. 25 ACS Paragon Plus Environment

Langmuir

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 26 of 36

(38) Gao, F.; Wang, Y.; Goodman, D. W. Reaction Kinetics and Polarization-Modulation Infrared Reflection Absorption Spectroscopy (PM-IRAS) Investigation of CO Oxidation over Supported Pd, Journal of Physical Chemistry C 2010, 114, 4036-4043. (39) Enache, D. I.; Edwards, J. K., Landon, P.; Solsona-Espriu, B.; Carley, A. F.; Herzing, A. A.; Watanabe, M.; Kiely, C. J.; Knight, D. W.; Hutchings, G. J. Solvent-Free Oxidation of Primary Alcohols to Aldehydes Using Au-Pd/TiO2 Catalysts," Science 2006, 311, 362-365. (40) Edwards, J. K.; Solsona, B.; N, E. N.; Carley, A. F.; Herzing, A. A.; Kiely, C. J.; Hutchings, G. J. Switching Off Hydrogen Peroxide Hydrogenation in the Direct Synthesis Process, Science 2009, 323, 1037-1041. (41) Chen, M.; Kumar, D.; Yi, C.-W.; Goodman, D. W. The Promotional Effect of Gold in Catalysis by Palladium-Gold, Science 2005, 310, 291-293. (42) Carino, E. V.; Crooks, R. M. Characterization of Pt@Cu Core@Shell DendrimerEncapsulated Nanoparticles Synthesized by Cu Underpotential Deposition, Langmuir 2011, 27, 4227-4235. (43) Zhang, J.; Vukmirovic, M. B.; Sasaki, K.; Nilekar, A.U.; Mavrikakis, M.; Adzic, R. R. Mixed-Metal Pt Monolayer Electrocatalysts for Enhanced Oxygen Reduction Kinetics,

Journal of the American Chemical Society 2005, 127, 12480-12481. (44) Xu, Y.; Ruban, A. V.; Mavrikakis, M. Adsorption and Dissociation of O2 on Pt−Co and Pt−Fe Alloys, Journal of the American Chemical Society 2004, 126, 4717-4725. (45) Kitchin, J. R.; Reuter, K.; Scheffler, M. Alloy Surface Segregation in Reactive Environments: First-principles Atomistic Thermodynamics Study of Ag3Pd(111) in Oxygen Atmospheres, Physical Review B 2008, 77, 075437.

26 ACS Paragon Plus Environment

Page 27 of 36

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

Langmuir

(46) Zhang, J.; Vukmirovic, M. B.; Xu, Y.; Mavrikakis, M.; Adzic, R. R. Controlling the Catalytic Activity of Platinum-Monolayer Electrocatalysts for Oxygen Reduction with Different Substrates, Angewandte Chemie International Edition 2005, 44, 2132–2135. (47) Ramos-Sanchez, G.; Praserthdam, S.; Godinez-Salomon, F.; Barker, C.; Moerbe, M.; Calderon, H. A.; Lartundo, L. A.; Leyva, M. A.; Solorza-Feria, O.; Balbuena, P. B. Challenges of Modelling Real Nanoparticles: Ni@Pt Electrocatalysts for the Oxygen Reduction Reaction, Physical Chemistry Chemical Physics 2015, 17, 28286-28297. (48) Viñes, F.; Desikusumastuti, A.; Staudt, T.; Görling, A.; Libuda, J.; Neyman, K. M. A Combined Density-Functional and IRAS Study on the Interaction of NO with Pd Nanoparticles: Identifying New Adsorption Sites with Novel Properties, The Journal of

Physical Chemistry C 2008, 112, 16539-16549. (49) Viñes, F.; Neyman, K. M.; Görling, A. Carbon on Platinum Substrates: From Carbidic to Graphitic Phases on the (111) Surface and on Nanoparticles, The Journal of Physical

Chemistry A 2009, 113, 11963-11973. (50) Viñes, F.; Loschen, C.; Illas, F.; Neyman, K. M. Edge Sites as a Gate for Subsurface Carbon in Palladium Nanoparticles, Journal of Catalysis 2009, 266, 59-63. (51) Bowker, M.; Madix, R. J. XPS, UPS and Thermal Desorption Studies of Alcohol Adsorption on Cu(110), Surface Science 1980, 95, 190-206. (52) Blöchl, P. E. Projector Augumented-wave Method, Physical Review B 1994, 50, 1795317979. (53) Bengtsson, L. Dipole Correction for Surface Supercell Calculations, Physical Review B 1999, 59, 12301-12304. (54) Xu, B.; Friend, C. M. Oxidative Coupling of Alcohols on Gold: Insights From

27 ACS Paragon Plus Environment

Langmuir

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 28 of 36

Experiments and Theory, Faraday Discussions 2001, 152, 307-320. (55) Shah, S. A. A. Forecasting Selectivity of Au-based Partial Oxidation Catalysts via Temperature Programmed Desorption Studies on the Au(111) Model Catalyst - Master Thesis 2014. (56) Francis, S. M.; Corneille, J.; Goodman, D. W.; Bowker, M. The Adsorption And Oxidation of Methanol on Thin Palladium Films, Surface Science 1996, 364, 30-38. (57) Savara, A.; Schmidt, C. M.; Geiger, F. M.; Weitz, E. Adsorption Entropies and Enthalpies and Their Implications for Adsorbate Dynamics, The Journal of Physical

Chemistry C 2009, 113, 2806-2815. (58) Ketchie, W. C.; Murayama, M.; Davis, R. J. Selective Oxidation of Glycerol Over Carbon-Supported AuPd Catalysts, Journal of Catalysis 2007, 250, 264-273. (59) Lichtenberger, J.; Lee, D.; Iglesia, E. Catalytic Oxidation of Methanol on Pd Metal and Oxide Clusters at Near-ambient Temperatures, Physical Chemistry Chemical Physics 2007, 9, 4902-4906. (60) Czelej, K.; Ćwieka, K.; Kurzydłowski, K. J. CO2 Stability on the Ni Low-index Surfaces: van der Waals Corrected DFT Analysis, Catalysis Communications 2016, 80, 33-38. (61) Liu, Z.; Zhang, R.; Wang, S.; Li, N.; Sima, R.; Liu, G.; Wu, P.; Zeng, G.; Li, S.; Sun, Y. Study, Highly Efficient and Stable Vanadia–Titania–Sulfate Catalysts for Methanol Oxidation to Methyl Formate: Synthesis and Mechanistic, The Journal of Physical

Chemistry C 2016, 120, 6591-6900. (62) Philips, K. R.; Jensen, S. C.; Baron, M.; Li, S.-C.; Friend, C. M. Sequential Photooxidation of Methanol to Methyl Formate on TiO2(110), Journal of the American

Chemical Society 2013, 135, 574-577.

28 ACS Paragon Plus Environment

Page 29 of 36

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

Langmuir

(63) Colmenares, J. C.; Lisowski, P.; Łomot, D.; Chernyayeva, O.; Lisovytskiy, D. Sonophotodeposition of Bimetallic Photocatalysts Pd–Au/TiO2: Application to Selective Oxidation of Methanol to Methyl Formate, ChemSusChem 2015, 8, 1676-1685. (64) Mullen, G. M.; Zhang, L.; Evans, E. J.; Yan, T.; Henkelman, G.; Mullins, C. B. Control Of Selectivity In Allylic Alcohol Oxidation On Gold Surfaces: The Role Of Oxygen Adatoms And Hydroxyl Species, Physical Chemistry Chemical Physics 2015, 174730174738.

29 ACS Paragon Plus Environment

Langmuir

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 30 of 36

TOC GRAPHIC

30 ACS Paragon Plus Environment

Page 31 of 36

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

Langmuir

Figure 1. Optimized geometries of the most stable adsorption species involved in the selective partial oxidation of methanol to MF on Au(111), Pd(111) and core@shell Au/Pd(111), Pd/Au(111) catalysts. 190x204mm (300 x 300 DPI)

ACS Paragon Plus Environment

Langmuir

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

Figure 2. Adsorption energies of molecules involved in the selective partial oxidation of methanol to MF on clean and preadsorbed oxygen Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) surfaces. 140x193mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 36

Page 33 of 36

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

Langmuir

Figure 3. The d-projected density of states (d-DOS) together with the d-band center calculated for the top layer of Au(111), Pd(111), Au/Pd(111) and Pd/Au(111) bare surfaces. Total charges per atom of first and second layer of each slab are shown on the right side. 190x196mm (300 x 300 DPI)

ACS Paragon Plus Environment

Langmuir

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

Figure 4. The Gibbs free energies of adsorption as a function of temperature calculated for molecules involved in the selective partial oxidation of methanol to MF on Au(111), Pd(111) and core@shell Au/Pd(111), Pd/Au(111) catalysts. 189x129mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 34 of 36

Page 35 of 36

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

Langmuir

Figure 5. The Gibbs free energies of adsorption as a function of temperature calculated for molecules involved in the selective partial oxidation of methanol to MF on oxygen preadsorbed Au(111), Pd(111) and core@shell Au/Pd(111), Pd/Au(111) catalysts. 189x132mm (300 x 300 DPI)

ACS Paragon Plus Environment

Langmuir

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

TOC Graphic 47x26mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 36 of 36