o-MO2(110), M = Ru and Ir) - ACS Publications - American

Apr 17, 2017 - ABSTRACT: Direct methanol fuel cell (DMFC) is an efficient power source. However, the DMFC anodes are easily toxified by CO or other ...
0 downloads 0 Views 1MB Size
Subscriber access provided by Fudan University

Article

A First-Principles Study on CO Removing Mechanism on Pt Decorated Oxygen-Rich Anode Surfaces (Pt/o-MO(110), M = Ru and Ir) in DMFC 2

2

Chi-You Liu, Chun-Chih Chang, Jia-Jen Ho, and Elise Yu-Tzu Li J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b13051 • Publication Date (Web): 17 Apr 2017 Downloaded from http://pubs.acs.org on April 19, 2017

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

The Journal of Physical Chemistry C 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 29

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

The Journal of Physical Chemistry

A First-Principles Study on CO Removing Mechanism on Pt Decorated Oxygen-rich Anode Surfaces (Pt2/o-MO2(110), M = Ru and Ir) in DMFC

Chi-You Liu, Chun-Chih Chang, Jia-Jen Ho, Elise Y. Li* Department of Chemistry, National Taiwan Normal University No. 88, Section 4, Tingchow Road, Taipei 116, Taiwan

* Corresponding author E-mail: [email protected] Phone: (886)-2-77346219 Fax: (886)-2-2932424 1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Abstract Direct methanol fuel cell (DMFC) is an efficient power source. However, the DMFC anodes are easily toxified by CO or other hydrocarbons, which terminates the methanol oxidation reaction (MOR). The most commonly used high performance catalyst on DMFC anodes is Pt or bimetallic PtRu. In this work, we apply densityfunctional theory (DFT) to investigate the adsorption of CO and H2O on pristine Pt2/MO2(110) and the oxygen-rich Pt2/o-MO2(110) surfaces (M = Ru and Ir). We find that the application of the oxygen-rich surfaces significantly reduces the adsorption energies of CO and H2O molecules as well as the major reaction barrier (CO + OH → CO2) in the water-gas-shift-like (WGS-like) reactions forming CO2. Our detailed analyses on the electronic interaction between the catalysts and adsorbates indicate that Pt2/o-MO2(110) may be a promising DMFC anode material which reduces the poison problem, and that it may be the actual experimental system that is responsible for the observed efficient CO removal.

2

ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29

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

The Journal of Physical Chemistry

I. Introduction Direct methanol fuel cell (DMFC) represents a promising renewable power source for portable and mobile applications.1 During the operation, methanol oxidation reaction (MOR, CH3OH + H2O → CO2 + 6H+ + 6e-) and oxygen reduction reaction (ORR, 1.5O2 + 6H+ + 6e- → 3H2O) occur at the anode and the cathode, respectively. The most common catalysts for DMFC anodes are Pt-based materials.1-8 DFT computations have confirmed that the reaction barriers of MOR are low on Pt(111) surface even without the inclusion of the solvation effect.8 Despite of their high performance in MOR, many issues persist with Pt-based catalysts, such as high cost, slow oxidation rate, and poor catalyst stability.9-11 The most severe of all is CO poisoning on the Pt-based anode that occurs during the stepwise dehydrogenation of MOR, leading to a non-reversible loss of reaction activity within several cycles.1-7, 12-14 Previous studies have shown that the adsorbed CO may be detached via

coupling

with OH species in aqueous environment through the water-gas-shift-like (WGS-like) reactions,1-3, 12, 15-21 as shown in Scheme 1. It has been shown that the application of Pt/Ru bimetallic meterials1-2 or the combination of Pt/RuO23 and Pt/IrO2,4-7 while maintaining the MOR activity, leads to an effective production of OH species on the catalyst surface. Among these, RuO2 and IrO2 surfaces have been shown to exert the highest catalytic activity for water dissociation1,

22

and remain stable in an acidic

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

environment.5, 23 The OH species generated on the RuO2 or IrO2 surfaces would thus couple efficiently with the CO molecules adsorbed at neighboring Pt active sites, recovering the MOR performance.1, 3, 12 In addition, RuO2 and IrO2 can act as benefit species for dispersing Pt into smaller clusters or atoms on surfaces.4-5, 12, 24-26 Previous studies have shown that the facile water splitting reaction driven by RuO2 or IrO2 surfaces results in the O or OH-covered surfaces, depending on the environment.1, 3, 12, 23-24, 27 Both experimental28-29 and theoretical30-34 studies indicated that in the presence of oxygen, the oxygen-covered (also called “oxygen-rich” or “orich”) RuO2 or IrO2 surfaces might form. Crihan et al.29 observed via STM the existence of fully occupied o-rich RuO2(110) surface at room temperature. First-Principles calculations reported energy barriers as low as 0.25 and 0.29 eV, respectively, for the dissociation of the oxygen molecule into the two coordinatively unsaturated O site (Ocus) on pristine RuO2(110)31 and IrO2(110)34 surfaces. The existence of such o-rich surfaces is often overlooked regarding the catalytic performance of RuO2 or IrO2 surfaces for MOR reactions. It is evident that a thorough examination for the catalytic reaction mechanism on these surfaces is necessary. In this study, we employ the DFT study to investigate the Pt atom dispersion as well as aggregation on pristine or o-rich metal oxide surfaces (MO2(110) and oMO2(110), M = Ru and Ir), and calculate the reaction barriers of WGS reactions on 4

ACS Paragon Plus Environment

Page 4 of 29

Page 5 of 29

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

The Journal of Physical Chemistry

both Pt2/MO2(110) and Pt2/o-MO2(110) surfaces. We characterize the RuO2(110) and IrO2(110) surface properties and analyze the effect of the adsorbed oxygen atoms on the WGS-like reactions. Our detailed analyses on the electronic interaction between the catalysts and adsorbates indicate that the o-MO2(110) surfaces significantly increase the efficiency of CO release via the WGS-like reactions, and successfully explains the high catalytic performance of MOR observed in previous experiments.

II. Computational Details The density-functional theory (DFT) and the plane-wave method, as implemented in the Vienna ab initio simulation package (VASP),35-38 are employed to calculate the energies and structures of reactants, intermediates, transition structures and products of the reactions on RuO2(110) and IrO2(110). The projector-augmented-wave method (PAW)39-40 is used in conjunction with the generalized gradient-approximation (GGA) and Perdew−Wang 1991 (PW91)41-42 exchange-correlation functional, with a cutoff energy of 400 eV. The energy convergence criteria for the electronic and the ionic steps are 10-4 and 10-3 eV, respectively. The Monkhorst-Pack k-point grids43 are set as 4 x 2 x 1 and 4 x 6 x 1 for the (2 x 2) and (6 x 2)-MO2(110) (M = Ru and Ir) supercells, respectively. A vacuum space of over 12 Å is introduced to ensure negligible interactions between slabs. All layers are fully relaxed during optimization. 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

The adsorption energies are calculated based on: Eads=E(slab

Page 6 of 29

+ adsorbate)-(Eslab

+

Eadsorbate), where E(slab + adsorbate), Eslab, and Eadsorbate are the computed electronic energies of the surface with adsorbed molecule, the clean surface and the free molecule, respectively. The transition state structures are located by nudged-elastic-band (NEB) method44-46, and the profiles of potential-energy surfaces (PES) are constructed accordingly. Frequency calculations are applied to verify the adsorbed intermediates and the transition states (with only one imaginary frequency). The Bader charge analysis is used in our calculation47-49. The electron density difference diagrams are plotted at the 0.003 |e|/Bohr3 isosurfaces.

III. Result and discussion Formation of Pt2/MO2(110) and Pt2/o-MO2(110) Figure 1(a)~(c) shows the structures of the pristine (2 x 2)-MO2(110) (M = Ru and Ir) surfaces, and the (2 x 2)-MO2(110) adsorbed with one or two Pt atoms. The positions of the bridging oxygen (Obr), the three-folded oxygen (O3f), the bridging metal (Mbr), and the coordinatively unsaturated metal (Mcus) sites are labelled. The optimized RuO2(110) and IrO2(110) surfaces give rise to almost identical structures as shown in Figure 1. The ELF diagram of the sliced MO2(110) surface shows that the electrons are localized at the Obr and O3f sites (red region) in both RuO2(110) and IrO2(110) systems, 6

ACS Paragon Plus Environment

Page 7 of 29

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

The Journal of Physical Chemistry

as shown in Figure S1. The enclosure of these atoms forms two possible hollow sites, H1 and H2, the center of Obr-O3f-O3f and Obr-Obr-O3f triangle, respectively, for Pt adsorption. Pt atoms are found to preferentially attach to H1 sites when sequentially added onto the MO2(110) surface, as shown by the optimized structures in Figure 1(b) and (c). No stable H2-adsorbed Pt structures can be located. The adsorption energies of the first Pt atom is on the range of 3.7 to 4.1 eV on MO2(110) surface, as shown in Table 1. When the second Pt atom is attached (Figure 1(c)), the two Pt atoms exert strong mutual attraction and shift considerably away from the center of the H1 site (0.09 Å and 0.15 Å per Pt atom for RuO2(110) and IrO2(110), respectively.) toward each other, signifying substantial Pt aggregation in the Pt2/MO2(110) structure. A larger (6 x 2) -MO2(110) super-cell, as shown in Figure S2, is constructed to further investigate the Pt aggregation effect. Table S1 lists the co-adsorption energies of the two Pt atoms adsorbed on combinations of H1 sites (A ~ H). The two Pt atoms exhibit maximum shifting from the H1 center when adsorbed on the nearest H1 sites (site A and B in Figure S2), and the co-adsorption energies are 0.8 ~ 0.9 eV higher than other configurations. In other words, the two adsorption events cannot be treated as independent adsorbates as the co-adsorption energies (> 8 eV) are much larger than twice of the first adsorption energies, as shown in Table 1 and Table S1. We then consider the Pt adsorption on o-rich surfaces. As previously mentioned,31, 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

34

O2 molecules dissociate spontaneously with a small energy barrier (< 0.3 eV) on

RuO2(110) and IrO2(110) surfaces to form two Ocus atoms on top of Rucus and Ircus sites, as shown in Figure 1(d). Pt atoms preferentially attach to the bridging sites between the two Ocus atoms when adsorbed on the fully oxygen covered o-MO2(110) surfaces (Figure 1(e) and (f)). The adsorption energies of Pt atom on o-MO2(110) (> 5 eV) are much larger than those on pristine MO2(110) (< 4.1 eV) , as shown in Table 1. The adsorption of the second Pt atom is almost independent from the first adsorption event, as evidenced by the near identical first and second adsorption energies. The distances between the two adsorbed Pt atoms in Pt2/o-MO2(110) are over 6 Å , significantly longer than those in Pt2/MO2(110) (~2.93 Å ). In other words, Pt atoms will disperse much better on o-rich surfaces to form the Pt decorated o-rich surfaces (Ptn/o-MO2(110) surfaces). This occurrence also rationalizes previous experimental observations of better Pt dispersion with the inclusion of RuO2 and IrO2 species.4-5, 12, 24-26 The partial density of states (PDOS) analysis of the d-band states for Pt atoms in gas phase and for Pt atoms adsorbed on MO2(110) or o-MO2(110) surfaces are also performed and the results are shown in Figure 2(a)~(e). Pt2/MO2(110) shows a broad Pt d-band distribution, reflecting substantial d orbital interaction due to the short Pt-Pt distance. As mentioned previously, upon co-adsorption on MO2(110), the Pt atoms shift away from the H1 center, moving toward each other. The PDOS for the fully relaxed 8

ACS Paragon Plus Environment

Page 8 of 29

Page 9 of 29

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

The Journal of Physical Chemistry

Pt2/MO2(110) and the fixed Pt2/MO2(110), as shown in Figure S3, indicates that this additional aggregation induces a slightly broader d-band distribution. In the contrary, the Pt2/o-MO2(110) shows a much more concentrated, gas-phase like Pt d-band distribution.

This is uncommon for metal adsorption on metal oxide surfaces

considering the large adsorption energy (> 5 eV per Pt) for this system. Adsorption of CO and H2O on Pt2/MO2(110) and Pt2/o-MO2(110) Figure 3 shows the most stable adsorption structures of CO and H2O adsorbed on Pt2/RuO2(110) and Pt2/o-RuO2(110) surfaces. The optimized structures for those on Pt2/IrO2(110) and Pt2/o-IrO2(110) surfaces are very similar. The corresponding adsorption energies and relevant bond length variations for CO and H2O are listed in Table 2 and Table S2, respectively. The CO molecule attaches onto the Pt atom and forms approximately a 45-degree angle with the surface, as shown in Figure 3(a) and (b). Surprisingly, while the C-Pt bond lengths on o-rich surfaces are slightly longer (by about 0.035 Å ) than those on pristine surfaces, the adsorption energies of CO on Pt2/oMO2(110) surfaces is decreased by over 0.75 eV than those on Pt2/MO2(110) surfaces. For the adsorption of H2O, the water molecule attaches onto the Pt atom via the oxygen atom (Ow) and forms a hydrogen bond with the nearby Obr atom (Figure 3(c) and (d)), as evidenced by the typical hydrogen bond distance between the hydrogen atom (Ha) and Obr (~1.60 Å ), and the slightly lengthened Ow-Ha bond (~1.03 Å ) (Table S2). For 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

water adsorption, the adsorption energies are also reduced by 0.40 ~ 0.65 eV on the orich surfaces. To further investigate the adsorbate-surface interaction, we carried out the Bader charge analysis for the adsorption events on pristine and o-rich MO2(110) systems, as shown in Table S3~S6. For Pt2/MO2(110) and Pt2/o-MO2(110) surfaces, Pt atoms are found to carry more positive charges when adsorbed on o-rich surfaces (q ~ 0.92 |e|) than on pristine surfaces (q < 0.3 |e|) (Table S4). The excessive positive charge of the Pt atom on the o-rich surface is due to strong electron-withdrawing by the two neighboring under-saturated Ocus atoms, which accumulate an additional negative charge (Δq = -0.11 ~ -0.16 |e|) upon Pt adsorption. When CO is adsorbed to Pt2/MO2(110), the charge of Pt increases by about 0.2 |e|, indicating significant electron donation from Pt to CO (Table S5). For the Pt2/o-MO2(110) surfaces, however, the two active Ocus atoms continue to extract more negative charge (Δq = -0.11 ~ -0.13 |e|) upon subsequent CO adsorption and suppress the electron back-donation from Pt to CO. This unique electron distribution in o-rich surfaces results in a weaker Pt-CO bonding and thus lower CO adsorption energies. The above-mentioned interactions between CO adsorbates and surfaces can be directly visualized by the electron density difference before and after CO adsorption as shown in Figure 4(a)~(d). Upon CO adsorption, both the Pt2/MO2(110) and Pt2/o10

ACS Paragon Plus Environment

Page 10 of 29

Page 11 of 29

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

The Journal of Physical Chemistry

MO2(110) surfaces show a net electron flow from the Pt atom (the yellow region) towards CO (the blue region), and from neighboring oxygen atoms on the surface, Obr or Ocus on MO2(110) or o-MO2(110), respectively, to the Pt atom. Nevertheless, a more careful inspection reveals that the electron accumulation near the active Pt atom is reduced (the blue region marked by the red arrow in Figure 4(c) and (d)) for CO adsorption on the Pt2/o-MO2(110) than on the Pt2/MO2(110) surface. In other words, the Ocus atoms on the o-rich surfaces are more reluctant to share electrons with the Pt atom, which is consistent with the Bader charge analysis results. The Bader charge analysis also reveals that the H2O adsorption results in a net charge accumulation on the Ocus atom opposite to the Pt-Ow bond via the trans-effect. The electron density difference plot, as shown in Figure S4, before and after water adsorption indicates clearly the formation of the hydrogen bond (Ow-Ha-Obr) between water molecule and surfaces, as a significant amount of electron density flows from the Obr to the H atom of water can be visualized. The formation of the hydrogen bond shall facilitate the Ow-Ha bond breaking to generate OH species, as previously observed experimentally1, 3, 12. WGS reactions on Pt2/MO2(110) and Pt2/o-MO2(110) surfaces We then consider the WGS reactions to convert CO on both Pt2/MO2(110) and Pt2/o-MO2(110) surfaces. The potential energy surface (PES) diagram and the detailed 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

structures for all intermediates are shown in Figure 5 and Figure 6, respectively. First, the CO and H2O molecules are co-adsorbed on the Pt-decorated MO2 surfaces (IM1 in Figure 6). The two species are found to attach to neighboring Pt atoms in Pt2/MO2(110) system, but to the same Pt atom in Pt2/o-MO2(110) system. The adsorbed H2O then dissociates to form OH and H species (IM2), through a low reaction barrier (Ea < 0.35 eV) assisted by the hydrogen bond interaction on these surfaces. The OH species would then couple with the CO adsorbate to form OCOH intermediate on the Pt atom (IM3). This step (CO + OH → OCOH) yields the largest reaction barrier in the entire WGS reactions; the activation barriers on Pt2/RuO2(110), Pt2/IrO2(110), Pt2/o-RuO2(110) and Pt2/o-IrO2(110) surfaces are 1.93, 1.94, 1.23, and 1.46 eV, respectively. At the [OH(a) + CO(a)] intermediate state (IM2 in Figure 6), the distances between OH and CO groups (H-O…C-O) are 2.91, 2.70, 3.04, and 2.66 Å on RuO2(110), o-RuO2(110), IrO2(110), and o-IrO2(110), respectively. The OH and CO groups are a lot closer when adsorbed on the Pt2/o-MO2(110), resulting in easier combination on oxygen-rich surfaces. Note that if the CO and H2O are adsorbed on separate Pt atoms in Pt2/o-MO2(110), the reaction would be almost impossible to proceed with a distance between the adsorbates over 6 Å . The OCOH intermediate then undergoes the O-H bond scission to form CO2 and then become desorbed (IM4). To complete the catalytic cycle, the H-H coupling barriers are calculated to assess 12

ACS Paragon Plus Environment

Page 12 of 29

Page 13 of 29

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

The Journal of Physical Chemistry

the probability of forming the hydrogen molecule products. The H-H coupling barriers on all surfaces are over 1 eV and are endothermic (IM5). The desorption energies of the hydrogen molecule, 0.66~1.49 eV, are also quite high (IM6). In this case, the CO conversion reactions may stop at the CO2 and the protonated surface (IM4) or proceed to the final CO2 + H2 products (IM6).

VI. Conclusion Our DFT investigation shows that the o-rich surfaces (o-MO2(110)) facilitate Pt dispersion more effectively than pristine surfaces. The calculation results indicate that the experimentally observed Pt dispersion may have occurred on the Pt2/o-MO2(110) surfaces rather than the pristine MO2(110), as previously assumed. The Pt decorated oxygen-rich surfaces (Pt2/o-MO2(110)) significantly decrease the CO adsorption energies by about 0.8 eV than pristine MO2(110) surfaces, via the participation of the Ocus sites on o-MO2(110). Furthermore, the activation energy barriers of WGS reactions are significantly reduced on Pt2/o-MO2(110). In summary, we believe that the formation of the o-rich surfaces is partially responsible for the observed high performance in catalytic reactions of DMFC anodes.

Supporting Information 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Detailed electronic analysis and structures of MO2(110) and o-MO2(110).

Acknowledgments This research is supported by the Ministry of Science and Technology (MOST) in Taiwan (MOST 105-2113-M-003-008 and MOST 102-2113-M-003-008–MY4). We thank National Center for High-performance Computing (NCHC) of Taiwan for the help on computational resources.

Reference (1) Kakati, N.; Maiti, J.; Lee, S.-H.; Jee, S.-H.; Viswanathan, B.; Yoon, Y.-S. Anode Catalysts for Direct Methanol Fuel Cells in Acidic Media: Do We Have Any Alternative for Pt or Pt-Ru? Chem. Rev. 2014, 114, 12397-12429. (2) Arico, A. S.; Srinivasan, S. Antonucci, V. DMFCs: From Fundamental Aspects to technology Development. Fuel Cells 2001, 1, 133-161. (3) Liu, R.-J.; Luo, Z.-H.; Wei, Q.-M.; Zhou, X.-S. Pt-RuO2 nanoparticles supported on diaminoanthraquinone-functionalized carbon nanotubes as efficient catalysts for methanol oxidation. Mater. Des. 2016, 94, 132-138. (4) Wang, H.; Zheng, J.; Peng, F.; Yu, H. Pt/IrO2/CNT anode catalyst with high performance for direct methanol fuel cells. Catal. Commun. 2013, 33, 34-37. (5) Chen, A.; La Russa, D. J.; Miller, B. Effect of the Iridium Oxide Thin Film on the 14

ACS Paragon Plus Environment

Page 14 of 29

Page 15 of 29

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

The Journal of Physical Chemistry

Electrochemical Activity of Platinum Nanoparticles. Langmuir 2004, 20, 9695-9702. (6) Baglio, V.; Amin, R. S.; El-Khatib, K. M.; Siracusano, S.; D’Urso, C.; Arico, A. S. IrO2 as a promoter of Pt-Ru for methanol electro-oxidation. Phys. Chem. Chem. Phys. 2014, 16, 10414-10418. (7) Shan, C.-C.; Tsai, D.-S.; Huang, Y.-S.; Jian, S.-H.; Cheng, C.-L. Pt-Ir-IrO2NT ThinWall Electrocatalysts Derived from IrO2 Nanotubes and Their Catalytic Activities in Methanol Oxidation. Chem. Mater. 2007, 19, 424-431. (8) Sakong, S.; Groß, A. The Importance of the Electrochemical Environment in the Electro-Oxidation of Methanol on Pt(111). ACS Catal. 2016, 6, 5575-5586. (9) Antolini, E. The problem of Ru dissolution from Pt–Ru catalysts during fuel cell operation: analysis and solutions. J. Solid State Electrochem. 2011, 15, 455-472. (10) Cabello-Moreno, N.; Crabb, E. M.; Fisher, J. M.; Russell, A. E.; Thompsett, D. Impact of PtRu Anode Catalyst Degradation on DMFC MEA Performance. ECS Trans. 2008, 16, 483-496. (11) Gancs, L.; Hakim, N.; Hult, B. N.; Mukerjee, S. Dissolution of Ru from PtRu Electrocatalysts and its Consequences in DMFCs. ECS Trans. 2006, 3, 607-618. (12) Villullas, H. M.; Mattos-Costa, F. I.; Nascente, P. A. P.; Bulhoes, L. O. S. Sol-Gel Prepared Pt-Modified Oxide Layers: Synthesis, Characterization, and Electrocatalytic Activity. Chem. Mater. 2006, 18, 5563-5570. 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(13) Sundmacher, K.; Schultz, T.; Zhou, S.; Scott, K.; Ginkel, M.; Gilles, E. D. Dynamics of the direct methanol fuel cell (DMFC): experiments and model-based analysis. Chem. Eng. Sci. 2001, 56, 333-341. (14) Parasons, R.; Van Der Noot, T. The oxidation of small organic molecules: A survey of recent fuel cell related research. J. Electroanal. Chem. 1988, 257, 9. (15) Huang, W.-J.; Wang, H.-T.; Zhou, J.-G.; Wang, J.; Duchesne, P. N.; Muir, D; Zhang, P.; Han, N.; Zhao, F.-P.; Zeng, M.; Zhong, J.; Jin, C.-H.; Li, Y.-G.; Lee, S.-T.; Dai, H.J. Highly active and durable methanol oxidation electrocatalyst based on the synergy of platinum-nickel hydroxide-graphene. Nat. Commun. 2015, 6, 10035. (16) Davó-Quinoñero, A.; Navlani-García, M.; Lozano-Castelló, D.; Bueno-López, A.; Anderson, J. A. Role of Hyfroxyl Groups in the Preferential Oxidation of CO over Copper Oxide-Cerium Oxide Catalysts. ACS Catal. 2016, 6, 1723-1731. (17) Maestri, M.; Livio, D.; Beretta, A.; Groppi, G. Hierarchical Refinement of Microkinetic Models: Assessment of the Role of the WGS and r-WGS Pathways in CH4 Partial Oxidation on Rh. Ind. Eng. Chem. Sci. 2014, 53, 10914-10928. (18) Song, W.; Hensen, E. J. M. Mechanistic Aspects of the Water-Gas Shift Reaction on Isolated and Clustered Au Atoms on CeO2(110): A Density Functional Theory Study. ACS Catal. 2014, 4, 1885-1892. (19) Tang, Q.-L.; Liu, Z.-P. Identification of the Active Cu Phase in the Water-Gas Shift 16

ACS Paragon Plus Environment

Page 16 of 29

Page 17 of 29

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

The Journal of Physical Chemistry

Reaction over Cu/ZrO2 from First Principles. J. Phys. Chem. C 2010, 114, 8423-8430. (20) Liu, P.; Rodriguez, J. A. Water-gas-shift reaction on metal nanoparticles and surfaces. J. Chem. Phys. 2007, 126, 164705. (21) Lee, M.-J.; Kang, J.-S.; Kang, Y.-S.; Chung, D.-Y.; Shin, H.; Ahn, C.-Y.; Park, S.; Kim, M.-J.; Kim, S.; Lee, K.-S.; Sung, Y.-E. Understanding the Bifunctional Effect for Removal of CO Poisoning: Blend of a Platinum Nanocatalyst and Hydrous Ruthenium Oxide as a Model System. ACS Catal. 2016, 6, 2398-2407. (22) Rossmeisl, J.; Qu, Z.-W.; Zhu, H.; Kroes, G.-J.; Nørskov, J. K. Electrolysis of water on oxide surfaces. J. Electroanal. Chem. 2007, 607, 83-89. (23) Cherevko, S.; Geiger, S.; Kasian, O.; Kulyk, N.; Grote, J.-P.; Savan, A.; Shrestha, B. R.; Merzlikin, S.; Breitbach, B.; Ludwig, A.; Mayrhofer, K. J.J. Oxygen and hydrogen evolution reactions on Ru, RuO2, Ir and IrO2 thin film electrodes in acidic and alkaline electrolytes: A comparative study on activity and stability. Catal. Today 2016, 262, 170-180. (24) Park, S.; Shao, Y.; Liu, J.; Wang, Y. Oxygen electrocatalysts for water electrolyzers and reversible fuel cells: status and perspective. Energy Environ. Sci. 2012, 5, 93319344. (25) Ye, F.; Li, J.; Wang, X.; Wang, T.; Li, S.; Wei, H.; Li, Q.; Christensen, E. Electrocatalytic properties of Ti/Pt-IrO2 anode for oxygen evolution in PEM water 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

electrolysis. Int. J. Hydrogen. Energy. 2010, 35, 8049-8055. (26) Yao, W.; Yang, J.; Wang, J.; Nuli, Y. Chemical deposition of platinum nanoparticles on iridium oxide for oxygen electrode of unitized regenerative fuel cell. Electrochem. Commun. 2007, 9, 1029-1034. (27) Hansen, H. A.; Man, I. C.; Studt, F.; Abild-Pedersen, F.; Bligaard, T.; Rossmeisl, J. Electrochemical chlorine evolution at rutile oxide (110) surfaces. Phys. Chem. Chem. Phys. 2010, 12, 283-290. (28) Over, H.; Knapp, M.; Lundgren, E.; Seitsonen, A. P.; Schmid, M.; Varga, P. Visualization of Atomic Processes on Ruthenium Dioxide using Scanning Tunneling Microscopy. Chem. Phys. Chem. 2004, 5, 167-174. (29) Crihan, D.; Knapp, M.; Seitsonen, A. P.; Over, H. Comment on “Interaction of Hydrogen with RuO2(110) Surfaces: Activity Differences between Various Oxygen Species” J. Phys. Chem. B 2006, 110, 22947. (30) Wang, H.; Schneider, W. F. Effects of coverage on the structures, energetics, and electronics of oxygen adsorption on RuO2(110). J. Chem. Phys. 2007, 127, 064706. (31) Wang, H.; Schneider, W. F.; Schmidt, D. Intermediates and Spectators in O2 Dissociation at the RuO2(110) Surface. J. Phys. Chem. C 2009, 113, 15266-15273. (32) Atmaca, D. O.; Düzenli, D.; Ozbek, M. O.; Onal, I. A density functional theory study of propylene epoxidation on RuO2(110) surface. Appl. Surf. Sci. 2016, 385, 9918

ACS Paragon Plus Environment

Page 18 of 29

Page 19 of 29

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

The Journal of Physical Chemistry

105. (33) Kim, Y. D.; Seitsonen, A. P.; Wendt, S.; Wang, J.; Fan, C.; Jacobi, K.; Over, H.; Ertl, G. Characterization of Various Oxygen Species on an Oxide Surface: RuO2(110). J. Phys. Chem. B. 2001, 105, 3752-3758. (34) Chung, W.-H.; Wang, C.-C.; Tsai, D.-S.; Jiang, J.-C.; Cheng, Y.-C.; Fan, L.-J.; Yang, Y.-W.; Huang, Y.-S. Deoxygenation of IrO2(110) surface: Core-level spectroscopy and density functional theory calculation. Surf. Sci. 2010, 604, 118-124. (35) Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558-561. (36) Kresse, G.; Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metelamorphous-semiconductor transition in germanium. Phys. Rev. B 1994, 49, 1425114269. (37) Kresse, G.; Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comp. Mater. Sci. 1996, 6, 15-50. (38) Kresse, G.; Hafner, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169-11186. (39) Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 1795317979. 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(40) Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758-1775. (41) Perdew, J. P.; Wang, Y. Accurate and simple analytic representation of the electrongas correlation energy. Phys. Rev. B 1992, 45, 13244. (42) Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Penderson, M. R.; Singh, D. J.; Fiolhais, C. Atoms, molecules, soilds, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 1992, 46, 6671. (43) Monkhorst, H. J.; Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188. (44) Ulitsky, A.; Elber, R. A new technique to calculate steepest descent paths in flexible polyatomic systems. J. Chem. Phys. 1990, 92, 1510. (45) Mills, G.; Jónsson, H.; Schenter, G. K. Reversible work transition state theory: application to dissociative adsorption of hydrogen. Surf. Sci. 1995, 324, 305. (46) Henkelman, G.; Uberuaga, B. P.; Jónsson, H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. J. Chem. Phys. 2000, 113, 9901. (47) Bader, R. F. W. A quantum theory of molecular structure and its applications. Chem. Rev. 1991, 91, 893-928. 20

ACS Paragon Plus Environment

Page 20 of 29

Page 21 of 29

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

The Journal of Physical Chemistry

(48) Henkelman, G.; Arnaldsson, A.; Jόnsson, H. A fast and robust algorithm for Bader decomposition of charge density. Comput. Mater. Sci. 2006, 36, 354-360. (49) Sanville, E.; Kenny, S. D.; Smith, R.; Henkelman, G. J. Improved grid-based algorithm for Bader charge allocation. J. Comput. Chem. 2007, 28, 899-908.

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Scheme 1. The WGS-like reaction mechanism. Symbol (a) refers to the adsorbed state.

Figure 1. Optimized structures of the (a) RuO2(110), (b) Pt/RuO2(110), (c) Pt2/RuO2(110), (d) o-RuO2(110), (e) Pt/o-RuO2(110), and (f) Pt2/o-RuO2(110) surface. The definitions of Mcus, O3f, Obr, Mbr, and Ocus sites are indicated. The colors of the elements Ru, O, and Pt are Persian green, red, and blue, respectively.

22

ACS Paragon Plus Environment

Page 22 of 29

Page 23 of 29

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

The Journal of Physical Chemistry

Figure 2. Partial density of state (PDOS) projection of the d-band structures of two Pt atoms in (a) gas phase, (b) Pt2/RuO2(110), (c) Pt2/IrO2(110), (d) Pt2/o-RuO2(110), and (e) Pt2/o-IrO2(110), respectively.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 3. Calculated structures (side view) of CO adsorbed on (a) Pt2/RuO2(110) and (b) Pt2/o-RuO2(110) surfaces, and H2O adsorbed on (c) Pt2/RuO2(110) and (d) Pt2/oRuO2(110) surfaces. The colors of elements Ru, O, Pt, C, and H are Persian green, red, blue, gray, and white, respectively.

24

ACS Paragon Plus Environment

Page 24 of 29

Page 25 of 29

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

The Journal of Physical Chemistry

Figure 4. Calculated electron density difference diagrams (Isosurface = 0.003 |e|/Bohr3) before and after CO adsorption on (a) Pt2/RuO2(110), (b) Pt2/IrO2(110), (c) Pt2/oRuO2(110), and (d) Pt2/o-IrO2(110) surfaces. Blue and yellow represent charge accumulation and depletion regions, respectively. Red arrows indicate the region with a comparatively less electron accumulation upon CO adsorption for o-rich systems.

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 5. Calculated potential energy surface (PES) diagram for the WGS reaction on Pt2/RuO2(110), Pt2/IrO2(110), Pt2/o-RuO2(110) and Pt2/o-IrO2(110) surfaces. The intermediates (IMn, n=1-6) and the reaction barriers (in eV) for the transition states (TSn, n=1-5) for each reaction step are marked. The optimized structures for each intermediate are shown in Figure 6.

26

ACS Paragon Plus Environment

Page 26 of 29

Page 27 of 29

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

The Journal of Physical Chemistry

Figure 6. The calculated reaction intermediates of WGS reactions on (a) Pt2/RuO2(110), (b) Pt2/IrO2(110), (c) Pt2/o-RuO2(110), and (d) Pt2/o-IrO2(110) surfaces. The colors of the elements Ru, Ir, O, Pt, C, and H are Persian green, Lavender, red, blue, gray, and white, respectively.

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 29

Table 1. The adsorption energies (∆E, in eV) of one or two Pt atoms on RuO2(110), IrO2(110), o-RuO2(110), and o-IrO2(110) surfaces. Adsorption Energy ∆E (eV)

Adsorbate

a

RuO2(110)

IrO2(110)

Pt (H1)

-3.68

-4.09

2Pt (H1)

-8.46a

-9.40a

o-RuO2(110)

o-IrO2(110)

First Pt

-5.44

-5.98

Second Pt

-5.31

-5.87

The co-adsorption energy of two Pt atoms

Table 2. The adsorption energies (∆E, in eV) of CO and H2O molecules on Pt2/RuO2(110), Pt2/o-RuO2(110), Pt2/IrO2(110), and Pt2/o-IrO2(110) surfaces. Adsorbate

Adsorption Energy ∆E (eV) Pt2/RuO2(110)

Pt2/IrO2(110)

CO

-2.86

-2.84

H2O

-1.75

-1.67

Pt2/o-RuO2(110)

Pt2/ o-IrO2(110)

CO

-1.99

-2.07

H2O

-1.09

-1.24

28

ACS Paragon Plus Environment

Page 29 of 29

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

The Journal of Physical Chemistry

TOC Graphic The CO removal reaction catalyzed by Pt2/o-MO2(110) (M = Ru and Ir).

29

ACS Paragon Plus Environment