Density Functional Theory Investigation of the Role of Cocatalytic

May 25, 2018 - Shakeel, Li, Yasin, Arif, Rehman, and Khan. 2018 57 (24), pp 8152–8159. Abstract: Efficient visible-light-driven materials are key fo...
2 downloads 0 Views 2MB Size
Subscriber access provided by Kaohsiung Medical University

Kinetics, Catalysis, and Reaction Engineering

Density Functional Theory Investigation of the Role of Cocatalytic Water in Methane Steam Reforming over Anatase TiO (101) 2

Alec Hook, Timothy P. Nuber, and Fuat E. Celik Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.8b00944 • Publication Date (Web): 25 May 2018 Downloaded from http://pubs.acs.org on May 25, 2018

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 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 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.

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 39 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

Industrial & Engineering Chemistry Research

Density Functional Theory Investigation of the Role of Cocatalytic Water in Methane Steam Reforming over Anatase TiO2 (101) Alec Hook,1 Timothy P. Nuber,2 Fuat E. Celik1,* 1

Department of Chemical and Biochemical Engineering

2

Department of Mechanical and Aerospace Engineering Rutgers, The State University of New Jersey 98 Brett Road, Piscataway, NJ 08854, USA

* Corresponding author: Tel. +1 848 445 5558 E-mail address: [email protected] (Fuat E. Celik)

1 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 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 39

Abstract The methane steam reforming (MSR) mechanism on the anatase TiO2 (101) surface has been investigated using periodic density functional theory. Reaction energies for several pathways were calculated and pathways involving high energy intermediates were eliminated from consideration. The remaining pathways all involved formaldehyde and formyl as intermediates. Activation energy barriers were calculated for these, revealing that C-H activation in methane and elementary dehydrogenation reactions in general involved very high energy transition states when calculated on clean surfaces. When water was coadsorbed on reaction surfaces as chemisorbed OH and H, hydrogen transfer reactions were facilitated by the proximity between hydrogen donor and acceptor and the strong affinity of hydroxyl for hydrogen. As the product of such hydrogen transfer reactions was physisorbed water, water activation completes the catalytic cycle and the role of water can be termed cocatalytic. Cocatalytic water lowered the most difficult dehydrogenation barriers by 2-3 eV, including reducing the methane activation barrier from 4.96 eV to 2.95 eV, making activity on the TiO2 surface relevant at high MSR reaction temperatures. Due to the high barrier for CHO dehydrogenation, the water-gas shift reaction is expected to make CO2 the preferred product, rather than CO, during MSR over metalfree TiO2.

2 ACS Paragon Plus Environment

Page 3 of 39 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

Industrial & Engineering Chemistry Research

1. Introduction Hydrogen gas has potential as a clean-burning energy carrier that produces only water as a byproduct of combustion. Hydrogen fuel, when used in a fuel cell for stationary or transportation power, can provide higher thermodynamic efficiency than conventional internal combustion. However, elemental hydrogen is not found on earth, and must be chemically separated from one of its compounds such as water or methane. Commercially, methane steam reforming (MSR) is used to convert fossil methane to hydrogen gas and carbon monoxide (as syngas) at high temperature, providing over half of global hydrogen production.1-3 The reaction is endergonic and endothermic, requiring high energy input at high temperature.4-6 Together with MSR, the water-gas shift (WGS) reaction takes places, where carbon monoxide reacts with water to form carbon dioxide and hydrogen. CH4 + H2O → CO + 3 H2

CO + H2O → CO2 + H2

Many transition metal catalysts have been investigated for the methane steam reforming reaction, including Rh, Ru, Ni, Ir, Pt and Pd. The steam reforming mechanism on metals is suggested to proceed via successive dehydrogenation of CH4 all the way to atomic carbon before undergoing oxidation. The mechanism of MSR has been investigated both experimentally and computationally using density functional theory (DFT). The initial C-H activation in methane is often found to be the rate-determining step at higher temperatures,7-12 while at lower temperature the oxidation of carbon to CO was thought to be limiting.13 More recent computational modeling has considered a wider range of possible pathways, involving several oxygenated intermediates formed by reaction of CHx fragments with surface O and OH groups, revealing that atomic carbon may not necessarily be involved in the reaction pathway. On Pt(111), hydroxymethylene and formyl where found to be energetically favorable intermediates between carbon and carbon

3 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 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 39

monoxide,14 and hydroxymethylidyne played a similar role on Pd(111) and Pd(100).15 Reactive oxygen species on palladium oxides have been implicated in the formation of methoxyl on PdO(111) surfaces,16,17 and calculations over PdO(100)18 and PdO(101)19 have evoked formaldehyde as an intermediate due to the large barriers for CHx dehydrogenation. Formyl was suggested in the mechanism of methane oxidation over both PdO(100) and CeO2(111) as a low energy intermediate and precursor to CO, rather than atomic carbon.20 Typically, alumina is used as the support for metal nanoparticles in MSR, but several advantages have been suggested for using TiO2 as a support instead, including lower coke formation, wider range of feasible steam to carbon ratios, and lower reaction temperature for supported Ni particles.21-23 The Ni/TiO2 system has been investigated for reforming of alcohols as well, including methanol,24 ethanol,25 and glycerol,26 where the interaction of the metal with the support has been implicated in the determining catalyst stability. The application of TiO2 and other semiconductors in methane steam reforming is potentially interesting as a photochemical alternative to the high-temperature thermocatalytic route in current industrial practice. Such a process, if it could lower the high-temperature requirement of MSR, could potentially double the atom efficiency by reducing or eliminating the need for fuel burn for heat generation, reducing carbon dioxide emissions and improving resource sufficiency. The Yoshida group has been investigating supported precious metals such as Pt and Rh on TiO2,27 titanates,28,29 tantalates,30,31, and gallium oxide,32 for photocatalytic MSR under mild conditions. In the reaction of water and methane over Pt/TiO2, CO2 was the primary carbonaceous product, with trace amounts of CO and ethane observed, and strong evidence for the formation of surface formates. Recent experimental work in our own group has found that

4 ACS Paragon Plus Environment

Page 5 of 39 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

Industrial & Engineering Chemistry Research

metal-free anatase is active for MSR at 973 K under photocatalytic (UV) and thermocatalytic (dark) reaction conditions, yielding H2 and CO2 selectively. The WGS reaction has also been investigated both experimentally and computationally over supported metal nanoparticles, 33-41 where reducible metal oxides such as TiO2 have been implicated in the dissociation of water.42-44 The structure and behavior of water molecules on TiO2 surfaces has been probed with DFT methods to give insight into binding geometries and energies.45-48 In our previous work,49 we investigated WGS mechanisms on the (101) surface of anatase TiO2 with density functional theory calculations. The associative pathway with carboxyl as the intermediate was favored on the clean TiO2 surface. At high water coverage, adsorbed water molecules would take a cocatalytic role in hydrogen transfer steps by acting as a hydrogen acceptor or donor, reducing the activation energy barriers of such steps. This meant that the associative formate pathway was competitive with the carboxyl pathway when the cocatalytic role of water was taken into consideration. Formation of OH radicals was predicted to be restricted to photocatalytic reaction conditions due to the large formation energies, consistent with experimental reports.50-52 In the present work, DFT is applied to the MSR reaction on anatase TiO2 (101). A series of potential reaction intermediates formed from dehydrogenation, O addition, and OH addition reactions with methane and its derivatives is considered. Using binding energies and structures for these intermediates, potential energy surfaces for multiple reaction pathways are generated. The most energetically unfavorable pathways are excluded from further consideration. The remaining pathways all involve formaldehyde dehydrogenation to formyl as an intermediate, via either methoxyl or methylene intermediates. Activation energy barriers for these remaining pathways are calculated and compared. The most favorable pathways involve the formation of

5 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 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 6 of 39

formaldehyde and formyl as surface intermediates. Most C-H cleavage reaction steps possess high activation energy barriers. In this context, dissociatively adsorbed water can assist in dehydrogenation steps by acting as a hydrogen acceptor, similar to the effect seen in WGS. This cocatalytic effect of chemisorbed water lowered these barriers by up to 3 eV. The water-assisted dehydrogenation mechanisms are therefore essential in the description of the activity of anatase (101) for MSR, and help interpret the metal-free activity found in high-temperature MSR experiments. 2. Methods Periodic DFT calculations were performed using the VASP code53,54 within the generalized gradient approximation (GGA-PBE)55 using projector-augmented wave (PAW)56,57 potentials. Corrections to the on-site Coulomb interactions have been accounted for through the Hubbard term (DFT+U).58 Dudarev’s approach59 was applied to increase the accuracy of the band gap and reaction kinetics with a Ueff = 4 for titanium atoms. This Ueff was chosen to align both the band gap and lattice parameters with experimental values as well as prevent electron delocalization.60–62 The single-electron wave functions were expanded using plane waves with an energy cutoff of 400 eV. All oxide slabs were based on the (101) surface of anatase TiO2, and modeled by a (2 × 2) surface unit cell with four (101) layers for a total of 32 Ti atoms and 64 O atoms. The lattice constants used for anatase TiO2 were a = 3.784 Å and c = 9.515 Å.63 A vacuum layer of 24 Å was used to separate periodic images of the slab in the z direction (normal to the surface), and a dipole correction was applied and the electrostatic potential was adjusted to ensure that interaction between the surface slab and its periodic images was negligible.64 The Brillouin zone was sampled using a (2 × 2 × 1) Gamma-centered Monkhorst-Pack k-point mesh65 following a

6 ACS Paragon Plus Environment

Page 7 of 39 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

Industrial & Engineering Chemistry Research

convergence test for adsorbate binding energies with respect to sampling mesh size. All four layers were allowed to relax in all calculations. Binding energies converged with respect to number of crystallographic layers at slab thicknesses of four (101) layers. Calculated differences between four layers (32 Ti atoms) and five layers (40 Ti atoms) were 0.7 eV. Each of these reactions’ transition states were high in energy due in part to the large distance between active sites on the anatase (101) surface. Undercoordinated oxygen atoms, which were the preferred binding sites for many adsorbates including hydrogen atoms, are 3.8 Å apart. Any adsorbates that bind to these sites through the carbon atom must place carbon and hydrogen on these two distant sites in the final state, or else place the hydrogen in an energetically unfavorable binding site as in Figure 6. The transition state must therefore span a large physical distance without any stabilization from nearby surface atoms. In the case of hydrogen transfer between formyl and carbon monoxide, the hydrogen atom in this transition state was found to have a gas-like character, with no interaction with the surface.49 Another way to describe this effect is that when both products from the C-H bond scission prefer the same site and site density is low, the competition for binding sites leads to long migration distances for the smaller daughter fragment – the H atom – and consequently large barriers. Not all dehydrogenation reactions were equally difficult. C-H bond scission in formaldehyde was only 1.62 eV. The reaction energy of only 0.10 eV was also significantly less

25 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 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 39

than with the previous four examples. This was because CH2O and CHO are the only adsorbates where both the reactant and the product interact with the surface via C and O atoms, and as such the transition state is somewhat stabilized relative to other C-H bond scission reactions. The O-H bond cleavage barrier in CH2OH was facile in contrast to all C-H cleavage reactions, with a barrier of only 0.48 eV. This is readily explained by the preferred binding structure of hydroxymethyl in Figure 4. The oxygen atom is trivalent in this structure, with the dissociating hydrogen atom dissociate pointing towards the hydrogen acceptor site, and an undercoordinated surface oxygen atom only 1.82 Å away. Unlike hydrogen atom transfer reactions including all the dehydrogenation reactions above, O and OH addition reactions are not expected to benefit from cocatalytic water, and use water as a reactant instead. This means that the barriers for methyl oxidation to methoxyl, methylene oxidation to formaldehyde, and methyl hydroxylation to hydroxymethyl are unchanged in pathways where cocatalytic water is considered. All three of these reactions had small reaction energies of +0.11 eV, -0.02 eV, and -0.21 eV respectively, and used chemisorbed water as a reactant. In methyl and methylene oxidation reactions, the OH group in chemisorbed water acts as the oxidant, transferring an O atom to the CHx species and transferring H to a nearby lattice oxygen atom. The barriers for these reactions were similar, 2.06 eV and 2.00 eV respectively. These large barriers arise due to the concerted cleavage of the O-H bond while forming the C-O bond and a new O-H bond in the transition states. Attempts to find step-wise mechanisms were unsuccessful. Taking methylene oxidation as an example, the hydrogen atom from chemisorbed water transferred to the surface oxygen atom already bonded to one H atom. The product of this elementary step is formaldehyde and a water molecule formed from a lattice oxygen atom,

26 ACS Paragon Plus Environment

Page 27 of 39 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

Industrial & Engineering Chemistry Research

leaving an oxygen vacancy. This reaction therefore pays the same >1 eV energy penalty for forming the vacancy by placing two H atoms on the same oxygen atom. This energy is then recovered as the two hydrogen atoms diffuse away from each other on the surface, restoring the oxygen atom to its lattice position. Therefore, the overall 2.00 eV barrier for this step includes +1.18 eV to form the formaldehyde/water/VO state and the overall reaction energy of -0.02 eV includes the -1.20 eV recovered by healing the vacancy. The reaction of water with oxygen vacancies is discussed in more detail in reference (49). Table 3. Reaction energies and activation energy barriers in eV for elementary steps in MSR. Reaction ∆Ε Ea (2) CH4  CH3 + H +0.97 4.96 (2’) a CH4 + H2O(c)  CH3 + H+ H2O(p) +0.75 2.95 (4) c CH3 + H2O(c)  CH3O + 2H +0.11 2.06 (5) CH3  CH2 + H +0.94 4.55 (5’) a CH3 + H2O(c)  CH2 + H+ H2O(p) +0.72 1.68 (8) d CH3O + H  CH2O + 2H +0.81 4.33 (8’) a,d CH3O + H + H2O(c)  CH2O + 2H + H2O(p) +0.60 1.69 (9) CH2OH  CH2O + H +0.19 0.48 (9’) a CH2OH + H2O(c)  CH2O + H + H2O(p) -0.03 4.3 eV barriers for many dehydrogenation steps, especially the initial methane activation barrier to form surface methyl (4.96 eV), essentially making all MSR pathways equally unlikely. However, the addition of chemisorbed water to the anatase surface revealed the cocatalytic role of water in reducing the barriers of difficult dehydrogenation steps. The ability of the OH group to act as a hydrogen acceptor and form physisorbed water stabilized the most difficult dehydrogenation transition states by >2 eV. This led to much more feasible apparent activation energy barriers of 2.2-2.6 eV relative to gaseous methane and steam as reactants. Though differences were small, the methoxyl pathway was slightly preferred over the methylene pathways to form formaldehyde. Even elementary reactions without hydrogen transfer benefitted from the presence of chemisorbed water. Direct oxidation and hydroxyl mechanisms using chemisorbed water as a reactant were thermodynamically favored over reactions involving removing lattice oxygen to form oxygen vacancies. These results further support the proposed active role of adsorbed water in catalyzing reactions on TiO2 surfaces. The cocatalytic activity of water has the potential to dramatically alter potential energy landscapes on anatase surfaces, changing the expected surface species and

34 ACS Paragon Plus Environment

Page 35 of 39 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

Industrial & Engineering Chemistry Research

their concentrations under reaction conditions. In particular, elementary reactions of methane, methanol, and formaldehyde derivatives that previously had prohibitively large kinetic barriers on the clean anatase surface appear to become feasible in the presence of adsorbed water.

Supporting Information Table S1: Initial, final, and transition state structures for elementary steps.

Acknowledgements This work was supported by the National Science Foundation under grant number DGE0903675. The authors would like to thank Dr. Ashley M. Pennington for meaningful discussions on metal-free reactivity of anatase surfaces, and Jacob D. Massa for his work on building model TiO2 surfaces. References (1) Agrafiotis, C.; von Storch, H.; Roeb, M.; Sattler, C. Solar Thermal Reforming of Methane Feedstocks for Hydrogen and Syngas Production—A Review. Renewable Sustainable Energy Rev. 2014, 29, 656-682. (2) Barreto, L.; Makihira, A.; Riahi, K. The Hydrogen Economy in the 21st Century: A Sustainable Development Scenario. Int. J. Hydrog. Energy 2003, 28, 267-284. (3) Chen, H.; Ding, Y.; Cong, N.; Dou, B.; Dupont, V.; Ghadiri, M.; Williams, P. Progress in Low Temperature Hydrogen Production with Simultaneous CO2 Abatement. Chem. Eng. Res. Des. 2011, 89, 1774-1782. (4) Zanfir, M.; Gavriilidis, A. Catalytic Combustion Assisted Methane Steam Reforming in a Catalytic Plate Reactor. Chem. Eng. Sci. 2003, 58, 3947-3960. (5) Xu, J.; Froment, G. Methane Steam Reforming, Methanation and Water-Gas Shift: I. Intrinsic Kinetics. AIChE J. 1989, 35, 88-96. (6) Wang, F.; Tan, J.; Ma, L.; Leng, Y. Effects of Key Factors on Solar Aided Methane Steam Reforming in Porous Medium Thermochemical Reactor. Energy Convers. Manage. 2015, 103, 419-430. (7) Wei, J.; Iglesia, E. Structural Requirements and Reaction Pathways in Methane Activation and Chemical Conversion Catalyzed by Rhodium. J. Catal., 2004, 225, 116-127. (8) Ligthart, D.; van Santen, R.; Hensen, E. Influence of Particle Size on the Activity and Stability in Steam Methane Reforming of Supported Rh Nanoparticles. J. Catal. 2011, 280, 206-220. (9) Wei, J.; Iglesia, E. Reaction Pathways and Site Requirements for the Activation and Chemical Conversion of Methane on Ru-Based Catalysts. J. Phys. Chem. B 2004, 108, 7253-7262. (10) Wei, J.; Iglesia, E. Isotopic and Kinetic Assessment of the Mechanism of Reactions of CH4 with CO2 or H2O to Form Synthesis Gas and Carbon on Nickel Catalysts. J. Catal. 2004, 224, 370-383. (11) Wei, J.; Iglesia, E. Isotopic and Kinetic Assessment of the Mechanism of Methane Reforming and Decomposition Reactions on Supported Iridium Catalysts. Phys. Chem. Chem. Phys. 2004, 6, 3754-3759. (12) Wei, J.; Iglesia, E. Mechanism and Site Requirements for Activation and Chemical Conversion of Methane on Supported Pt Clusters and Turnover Rate Comparisons Among Noble Metals. J. Phys. Chem. B 2004, 13, 40944103.

35 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 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 36 of 39

(13) Jones, G.; Jakobsen, J.G.; Shim, S.S.; Kleis, J.; Andersson, M.P.; Rossmeisl, J.; Abild-Pedersen, F.; Bligaard, T.; Helveg, S.; Hinnemann, B.; Rostrup-Nielsen, J.R.; Chorkendorff, I.; Sehested, J.; Nørskov, J.K. First Principles Calculations and Experimental Insight into Methane Steam Reforming over Transition Metal Catalysts. J. Catal. 2008, 259, 147-160. (14) Qi, W.; Ran, J.; Wang, R.; Du, X.; Shi, J.; Ran, M. Kinetic Mechanism of Effects of Hydrogen Addition on Methane Catalytic Combustion over Pt(111) Surface: A DFT Study with Cluster Modeling. Comput. Mater. Sci. 2016, 111, 430-442. (15) Jørgensen, M.; Grönbeck, H. First-Principles Microkinetic Modeling of Methane Oxidation over Pd(100) and Pd(111). ACS Catal. 2016, 6, 6730-6738. (16) Chin, Y.; Iglesia, E. Elementary Steps, the Role of Chemisorbed Oxygen, and the Effects of Cluster Size in Catalytic CH4–O2 Reactions on Palladium. J. Phys. Chem. C 2011, 115, 17845−17855. (17) Chin, Y.; Buda, C.; Neurock, M.; Iglesia, E. Consequences of Metal–Oxide Interconversion for C–H Bond Activation during CH4 Reactions on Pd Catalysts J. Am. Chem. Soc. 2013, 135, 15425−15442. (18) Dianat, A.; Seriani, N.; Ciacchi, L.C.; Bobeth, M.; Cuniberti, G. DFT Study of Reaction Processes of Methane Combustion on PdO(100). Chem. Phys. 2014, 443, 53-60. (19) Van den Bossche, M.; Grönbeck, H. Methane Oxidation over PdO(101) Revealed by First-Principles Kinetic Modeling. J. Am. Chem. Soc., 2015, 137, 12035-12044. (20) Mayernick, A.D.; Janik, M.J. Methane Oxidation on Pd-Ceria: A DFT Study of the Mechanism over PdXCe1xO2, Pd, and PdO. J. Catal. 2011, 278, 16-25. (21) Shinde, V.M.; Madras, G. Catalytic Performance of Highly Dispersed Ni/TiO2 for Dry and Steam Reforming of Methane. RSC Adv. 2014, 4, 4817-4826. (22) Kho, E.T.; Scott, J.; Amal, R. Ni/TiO2 for Low Temperature Steam Reforming of Methane. Chem. Eng. Sci. 2016, 140, 161-170. (23) Li, K.; Jia, L.; Wang, X.; Pu, J.; Chi, B.; Li, J. Enhanced Methane Steam Reforming Activity and Electrochemical Performance of Ni0.9Fe0.1-supported Solid Oxide Fuel Cells with Infiltrated Ni-TiO2 Particles. Sci. Rep. 2016, 6, 35981. (24) Chiarello, G.L.; Aguirre, M.H.; Selli, E. Hydrogen Production by Photocatalytic Steam Reforming of Methanol on Noble Metal-modified TiO2. J. Catal. 2010, 273, 182-190. (25) Nichele, V.; Signoretto, M.; Menegazzo, F.; Rossetti, I.; Cruciani, G.; Hydrogen Production by Ethanol Steam Reforming: Effect of the Synthesis Parameters on the Activity of Ni/TiO2 Catalysts. Int. J. Hydrog. Energy 2014, 39, 4252–4258. (26) Rossetti, I.; Gallo, A.; Dal Santo, V.; Bianchi, C.L.; Nichele, V.; Signoretto, M.; Finocchio, E.; Ramis, G.; Di Michele, A. Nickel Catalysts Supported Over TiO2, SiO2 and ZrO2 for the Steam Reforming of Glycerol. ChemCatChem 2013, 5, 294-306. (27) Yoshida, H.; Hirao, K.; Nishimoto, J.I.; Shimura, K.; Kato, S.; Itoh, H.; Hattori, T. Hydrogen Production from Methane and Water on Platinum Loaded Titanium Oxide Photocatalysts. J. Phys. Chem. C 2008, 112, 55425551. (28) Shimura, K.; Yoshida, H. Hydrogen Production from Water and Methane over Pt-loaded Calcium Titanate Photocatalyst. Energy Environ. Sci. 2010, 3, 615-617. (29) Shimura, K.; Kawai, H.; Yoshida, T.; Yoshida, H. Bifunctional Rhodium Cocatalysts for Photocatalytic Steam Reforming of Methane over Alkaline Titanate. ACS Catal. 2012, 2, 2126-2134. (30) Yoshida, H.; Kato, S.; Hirao, K.; Nishimoto, J.I.; Hattori, T. Photocatalytic Steam Reforming of Methane Over Platinum-loaded Semiconductors for Hydrogen Production. Chem. Lett. 2007, 36, 430-431. (31) Shimura, K.; Kato, S.; Yoshida, T.; Itoh, H.; Hattori, T.; Yoshida, H. Photocatalytic Steam Reforming of Methane Over Sodium Tantalate. J. Phys. Chem. C 2010, 114, 3493-3503. (32) Shimura, K.; Yoshida, T.; Yoshida, H. Photocatalytic Activation of Water and Methane over Modified Gallium Oxide for Hydrogen Production. J. Phys. Chem. C 2010, 114, 11466-11474. (33) Fu, Q.; Saltsburg, H.; Flytzani-Stephanopoulos, M. Active Nonmetallic Au and Pt Species on Ceria-based Water-gas Shift Catalysts. Science, 2003, 301, 935-938. (34) Pistonesi, C.; Juan, A.; Irigoyen, B.; Amadeo, N. Theoretical and Experimental Study of Methane Steam Reforming Reactions Over Nickel Catalyst. Appl. Surf. Sci. 2007, 253, 4427-4437. (35) Li, Y.; Fu, Q.; Flytzani-Stephanopoulos, M. Low-temperature Water-gas Shift Reaction over Cu- and Ni-loaded Cerium Oxide Catalysts. Appl. Catal., B, 2000, 27, 179-191. (36) Rodríguez, J.A.; Ma, S.; Liu, P.; Hrbek, J.; Evans, J.; Perez, M. Activity of CeOx and TiOx Nanoparticles Grown on Au (111) in the Water-gas Shift Reaction. Science, 2007, 318, 1757-1760.

36 ACS Paragon Plus Environment

Page 37 of 39 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

Industrial & Engineering Chemistry Research

(37) Si, R.; Flytzani-Stephanopoulos, M. Shape and Crystal-Plane Effects of Nanoscale Ceria on the Activity of AuCeO2 Catalysts for the Water–Gas Shift Reaction. Angew. Chem. 2008, 120, 2926-2929. (38) Fu, Q.; Weber, A.; Flytzani-Stephanopoulos, M. Nanostructured Au–CeO2 Catalysts for Low-temperature Water–gas Shift. Catal. Lett. 2001, 77, 87-95. (39) Panagiotopoulou, P.; Kondarides, D.I. Effect of Morphological Characteristics of TiO2-supported Noble Metal Catalysts on their Activity for the Water–gas Shift Reaction. J. Catal. 2004, 225, 327-336. (40) Rodríguez, J. A.; Evans, J.; Graciani, J.; Park, J. B.; Liu, P.; Hrbek, J.; Sanz, J. F. High Water-Gas Shift Activity in TiO2 (110) Supported Cu and Au Nanoparticles: Role of the Oxide and Metal Particle Size. J. Phys. Chem. C 2009, 113, 7364-7370. (41) Park, J.B.; Graciani, J.; Evans, J.; Stacchiola, D.; Senanayake, S.D.; Barrio, L.; Liu, P.; Sanz, J.F.; Hrbek, J.; Rodriguez, J.A. Gold, Copper, and Platinum Nanoparticles Dispersed on CeOx/TiO2 (110) Surfaces: High Water-Gas Shift Activity and the Nature of the Mixed-Metal Oxide at the Nanometer Level. J. Am. Chem. Soc. 2009, 132, 356-363. (42) Henderson, M.A. Structural Sensitivity in the Dissociation of Water on TiO2 Single-Crystal Surfaces. Langmuir 1996, 12, 5093-5098. (43) Lindan, P.J.D.; Harrison, N.M.; Holender, J. M.; Gillan, M.J. First-principles Molecular Dynamics Simulation of Water Dissociation on TiO2 (110). Chem. Phys. Lett. 1996, 261, 246-252. (44) Bikondoa, O.; Pang, C.L.; Ithnin, R.; Muryn, C.A.; Onishi, H.; Thornton, G. Direct Visualization of DefectMediated Dissociation of Water on TiO2 (110). Nat. Mater. 2006, 5, 189-192. (45) Tilocca, A.; Selloni, A. Structure and Reactivity of Water Layers on Defect-Free and Defective Anatase TiO2(101) Surfaces. J. Phys. Chem. B. 2004, 108, 4743-4751. (46) Aschauer, U.; He, Y.; Cheng, H.; Li, S.-C.; Diebold, U.; Selloni, A. Influence of Subsurface Defects on the Surface Reactivity of TiO2: Water on Anatase (101). J. Phys. Chem. C 2010, 114, 1278-1284. (47) Sun, C.; Liu, L.-M.; Selloni, A.; Lu, G.Q.(M.); Smith, S.C. Titania-water Interactions: A Review of Theoretical Studies. J. Mater. Chem. 2010, 20, 10319-10334. (48) Aschauer, U.J.; Tilocca, A.; Selloni, A. Ab Initio Simulations of the Structure of Thin Water Layers on Defective Anatase TiO2(101) Surfaces. Int. J. Quantum Chem. 2015, 115, 1250–1257. (49) Hook, A.; Celik, F.E. Density Functional Theory Investigation of the Role of Cocatalytic Water in Water Gas Shift Reaction over Anatase TiO2 (101). Ind. Eng. Chem. Res. 2018, 57, 6830-6841. (50) Geng, Z.; Chen, X.; Yang, W.; Gui, Q.; Xu, C.; Dai, D.; Yang, X. Highly Efficient Water Dissociation on Anatase TiO2(101). J. Phys. Chem. C. 2016, 120, 26807-26813. (51) Quanjun, X.; Yu, J.; Wong, P.K. Quantitative Characterization of Hydroxyl Radicals Produced by Various Photocatalysts. J. Colloid Interface Sci. 2011, 357, 163-167. (52) Hirakawa, T.; Nosaka. Y. Properties of O2•-and OH• Formed in TiO2 Aqueous Suspensions by Photocatalytic Reaction and the Influence of H2O2 and Some Ions. Langmuir 2002, 18, 3247-3254. (53) Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using A Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15-50. (54) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using A PlaneWave Basis Set. Phys. Rev. B 1996, 54, 11169-11186. (55) Perdew, J.P; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1997, 78, 1396. (56) Blöchl, P. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979. (57) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758-1775. (58) Anisimov, V.I.; Korotin, M.A.; Mylnikova, A.S.; Kozhevnikov, A.V.; Korotin, D.M.; Lorenzana, J. Computation of Stripes in Cuprates within the LDA+U Method, Phys. Rev. B 2004, 70, 172501. (59) Dudarev, S.L.; Botton, G.A.; Savrasov, S.Y.; Humphreys, C.J;. Sutton, A.P. Electron-energy-loss Spectra and the Structural Stability of Nickel Oxide: An LSDA+U Study. Phys. Rev. B 1998, 57, 1505. (60) Allen, J.; Watson, G. Occupation Matrix Control of D- And F-Electron Localisations Using DFT + U. Phys. Chem. Chem. Phys. 2014, 16, 21016-21031. (61) Hu, Z.; Metiu, H. Choice of U for DFT+U Calculations for Titanium Oxides. J. Phys. Chem. C 2011, 115, 5841-5845. (62) Morgan, B.; Watson, G. A DFT+U Description of Oxygen Vacancies at the TiO2 Rutile (110) Surface. Surf. Sci. 2007, 601, 5034-5041.

37 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 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 38 of 39

(63) Howard, C.; Sabine, T.; Dickson, F. Structural and Thermal Parameters for Rutile and Anatase. Acta Crystallogr., Sect. B: Struct. Sci. 1991, 47, 462-468. (64) Neugebauer, J.; Scheffler, M. Adsorbate-Substrate and Adsorbate-Adsorbate Interactions of Na and K Adlayers on Al(111). Phys. Rev. B 1992, 46, 16067-16080. (65) Monkhorst, H.; Pack, J. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188-5192. (66) Henkelman, G.; Uberuaga, B.P.; Jonsson, H.J. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths, J. Chem. Phys. 2000, 113, 9901– 9904. (67) Ravi, M.; Ranocchiari, M. van Bokhoven, J.A. The Direct Catalytic Oxidation of Methane to Methanol-A Critical Assessment. Angew. Chem., Int. Ed. 2017, 56, 16464-16483. (68) Scanlon, D.O.; Dunnill, C.W.; Buckeridge, J.; Shevlin, S.A.; Logsdail, A.J.; Woodley, S.M.; Catlow, C.R.A.; Powell, M.J.; Palgrave, R.G.; Parkin, I.P.; Watson, G.W.; Keal, T.W.; Sherwood, P.; Walsh, A.; Sokol, A.A. Band Alignment of Rutile and Anatase TiO2, Nature Mater., 2013, 12, 798-801.

38 ACS Paragon Plus Environment

Page 39 of 39 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

Industrial & Engineering Chemistry Research

Table of Contents Graphic

39 ACS Paragon Plus Environment