AVS2018 Session HC+SS-WeA: Theory and Dynamics of Heterogeneously Catalyzed Reactions

Wednesday, October 24, 2018 2:20 PM in Room 201A
Wednesday Afternoon

Session Abstract Book
(294KB, May 6, 2020)
Time Period WeA Sessions | Abstract Timeline | Topic HC Sessions | Time Periods | Topics | AVS2018 Schedule

Start Invited? Item
2:20 PM HC+SS-WeA-1 First-Principles Kinetic Monte Carlo Simulation of CO Oxidation on PdO(101): Role of Oxygen Vacancies
Minkyu Kim, A. Asthagiri (Ohio State University)

CO oxidation on transition metal (TM) oxide surfaces has been widely studied both experimentally and theoretically; however, a healthy debate continues on the coupling between changes in oxide phase and surface reactivity. In this study, we investigated CO oxidation over PdO(101) surface, which has been proposed to be important in Pd oxidation catalysis. In contrast to earlier studies of CO oxidation on TM oxide surfaces, we incorporate neighbor effects of oxygen vacancies on all the elementary surface processes. We employ density functional theory (DFT) to map out the kinetics of 50 elementary surface processes. We find that barriers for elementary steps such as O vacancy, O2dissociation, and CO2formation can be decreased by 20-40% in the presence of O vacancies, while O2desorption is increased by 70%. Using the DFT-derived rate constants, we have developed a lattice-based kinetic Monte Carlo (kMC) framework that can simulate CO oxidation under both ultra-high vacuum (UHV) and reaction conditions.

Initially, the kMC simulations were performed under UHV conditions (low CO partial pressure, no O2pressure) as a function of increasing surface temperature. As the CO oxidation proceeds and the surface lattice O atom coverage is depleted, the CO oxidation rate decreases rapidly at 400 K; however, the rate sharply increases at temperatures above 450 K. At 450 K, we find a new complex pathway to CO2formation that is activated in the presence of O vacancies and is the source of the increased CO oxidation rate despite the depletion of surface oxygen atoms. These kMC results match UHV isothermal experiments under similar conditions. KMC simulations at steady state conditions of low CO and O2pressures (PCO: 510-9Torr / PO2: 1.510-8Torr), show that the rate-limiting step is O2dissociation and this elementary step requires the presence of adjacent oxygen vacancies to be activated at temperatures below 500 K. Without the addition of O vacancy neighbor effects to the kMC model, the PdO(101) surface would be inactive to CO oxidation because surface oxygens cannot be healed by gaseous O2. In addition, we will discuss kMC simulations under reaction conditions at varying CO/O2partial pressures.

2:40 PM HC+SS-WeA-2 Surface Reactivity of Activated CO2
Richard van Lent (Leiden University, Netherlands); Anton Walsh, Michael Gleeson (DIFFER, Netherlands); Ludo Juurlink (Leiden University, Netherlands)

Catalytically converting CO2 into renewable fuels is a promising avenue that addresses the current fuel and energy storage challenges. Depending on the process, conversion of CO2 may involve initial breaking of the OC=O bond. If so, this highly endothermic step is likely rate limiting to the overall process. Ultra‑high vacuum (UHV) research on well-defined surfaces can provide fundamental insight into such processes, e.g. how dissociation may be aided by internal energy in CO2, the metal’s identity and the surface structure.

Internal energy has been shown to promote the highly activated dissociative adsorption of CO2 on Ni(100) [1]. Rovibrationally state-resolved measurements for CH4 on Ni(100) have shown that excitation of the ν3 antisymmetric stretch vibration has a promoting effect that approximately equals kinetic energy for Ni(100) surfaces [2,3]. For other metal surfaces and CH4 vibrations, the relative efficacy of vibrational and kinetic energies varies between 0.4 and 1.4 [2].

To start unravelling how CO2 dissociates on a surface, we combine standard supersonic molecular beam techniques with high resolution continuous wave (cw) IR laser excitation. We study the effect of ν3 antisymmetric stretch excitation of CO2 on dissociative adsorption on a Ni(711) surface.

For excitation, we use a 3.9-4.6 μm single mode optical parametric oscillator (OPO). The OPO is frequency-stabilized by locking onto the derivative of the lamb dip of the specific rovibrational transition of interest. Frequency-stabilized IR radiation is crossed with a supersonic molecular beam that impinges onto the cleaned Ni(711) surface under UHV conditions. We use methods to determine absolute reactivities for CO2 molecules with and without laser excitation. We extract rovibrational state-dependent absolute sticking probabilities as a function of kinetic energy. From this data, we determine the relative efficacy of vibrational and kinetic energy in overcoming the large activation barrier to dissociation.

[1] M. P. D’Evelyn, A. V. Hamza, G. E. Gdowski, and R. J. Madix, Surf. Sci.167, 451 (1986).

[2] L. B. F. Juurlink, D. R. Killelea, and A. L. Utz, Prog. Surf. Sci. 84, 69 (2009).

[3] B. L. Yoder, R. Bisson, and R. D. Beck, Sci. 329 , 553 (2010).

3:00 PM Invited HC+SS-WeA-3 Shining Light on Complexity: State- and Energy-Resolved Studies of Gas-Surface Reaction Dynamics and Mechanism
Arthur Utz (Tufts University)

Gas-surface reactions on a metal surface can be far more complex than their net chemical equation would suggest. The heterogeneity of the surface, and its ability to restructure, impact the availability and energetics of reaction sites, energy redistribution processes triggered by the gas-surface collision dictate transition state access, and many coupled chemical processes can influence the reaction's outcome.

The talk will survey the ability to vibrational state selected measurements of gas-surface reactivity to systematically unravel important aspects of this complexity and provide molecular-level mechanistic insights that can guide our understanding and predictions about these processes. The experiments use infrared light to excite a single quantum rotational and vibrational level of the gas-surface reagent molecule and a molecular beam to define its incident kinetic energy. In this way, we are able to make energy-resolved measurements of reactivity while systematically varying the distribution of energy among the system's energetic degrees of freedom. Computational studies by our collaborators inform and guide our interpretation of the data.

Recent results will be featured and examine the impact of a vibrationally excited reagent molecule's symmetry on surface-induced energy redistribution, the quenching lifetime of a vibrationally excited molecule on a metal surface, the role of surface atom motion on reaction energetics, the ability to use these measurements to benchmark DFT electronic structure calculations, and how isolated surface structural features, including steps and adatoms, modify reaction energetics.

3:40 PM BREAK
4:20 PM HC+SS-WeA-7 Vibration-driven Reaction of CO2 on Cu Surfaces via Eley-Rideal Type Mechanism
Junji Nakamura, Jiamei Quan, Taijun Kozarashi, Tomoyasu Mogi, Takumi Imabayashi, Kotarou Takeyasu, Takahiro Kondo (University of Tsukuba, Japan)

We have long studied the formation of formate (HCOOa) intermediates as an intermediates of methanol synthesis by hydrogenation of CO2 on Cu surfaces [1-4]. It has been suggested based on the kinetic analysis for the high pressure bulb experiments that the formate formation proceeds via Eley-Rideal type mechanism, in which CO2 directly reacts with adsorbed hydrogen atom on Cu. In order to clarify the mechanism and dynamics, we carried out CO2 molecular beam studies with the translational energy of 1.12-1.97 eV and the nozzle temperature of 800-1100 K. It was found that hot CO2 in the molecular beam reacts directly with pre-adsorbed hydrogen atoms on cold Cu(111) and Cu(100) surfaces at 120–220 K to form formate adspecies (CO2 + Ha → HCOOa). That is, even at the low surface temperatures, formate species is formed rapidly when supplying energy only to CO2. This indicates that the energy to overcome the reaction barrier comes from the hot CO2 molecule itself instead of the Cu surface. The vibrational energy of CO2 was much more effective for the reaction compared to the translational energy; and the reaction rate was independent of the surface temperature. The independence of surface temperature indicates the E-R type mechanism that the CO2 molecule is not thermally equilibrated with the Cu surfaces but directly reacts with Ha. If the impinging CO2 molecule chemisorbs on the Cu surfaces before reacting with Ha, the reaction rate should be dependent on the surface temperature, rather than the kinetic and the internal energies of CO2. We evaluated the barrier distribution factors (W) of so-called Luntz equation by analysis of the reaction probability curves, which shows significant efficacy of the vibrational energy on the reaction of CO2, which has not been observed before. DFT calculations are consistent with the experimental observations [5]. The direct reaction of CO2 with adsorbed hydrogen atom can be regarded as the first example of vibration-driven bond formation reactions on surfaces.

References

  1. J. Nakamura, Y. Choi, T. Fujitani, Top. Catal. 22, 277 (2003).

  2. H. Nakano, I. Nakamura, T. Fujitani, and J. Nakamura, J. Phys. Chem. B 105, 1355 (2001).

  3. G.Wang, Y.Morikawa, T.Matsumoto, and J.Nakamura, J. Phys. Chem. B 110, 9 (2006).

  4. J. Quan, T. Kondo, G. Wang, and J. Nakamura, Angew. Chem. Int. Ed. 56, 3496 (2017).

  5. F. Muttaqien, H. Oshima, Y. Hamamoto, K. Inagaki, I. Hamada, Y. Morikawa, Chem. Commun. 53, 9222 (2017).

4:40 PM Invited HC+SS-WeA-8 First Principles Reaction Kinetics over Metals, Oxides and Nanoparticles
Henrik Grönbeck (Chalmers University of Technology, Gothenburg, Sweden)

A major challenge in heterogeneous catalysis research is the determination of dominating reaction paths and kinetic bottlenecks. One reason for the challenge is the dynamic character of the kinetics, where the active sites may change with reaction conditions. Nevertheless, it is atomic scale information that allow for catalyst development beyond trial-and-error approaches. Kinetic modeling based on first principles calculations have over the past decade grown into an important tool for investigating the importance of different catalyst phases and reaction paths. In this contribution, I will discuss work where we have used density functional theory in combination with kinetic modeling to investigate catalytic reactions over metals, oxides and nanoparticles. The examples cover different aspects of kinetic modeling including determination of adsorbate entropies, importance of attractive adsorbate-adsorbate interactions and the complexity of many types of active sites.

Complete methane oxidation to carbon dioxide and water is used as one example and we have investigated the reaction over metallic Pd(100) and Pd(111) [1] as well as PdO(101) [2]. The reaction paths are shown to be markedly different on the metallic and the oxidized surfaces. We find that the catalytic activity is highly sensitive to adsorbate-adsorbate interactions which for PdO(101) are attractive owing to electronic pairing effects. This effect is crucial and common for oxide surfaces [3].

Reactions over platinum nanoparticles are investigated using a recently developed scaling relation Monte Carlo technique [4]. Taking CO oxidation as a model reaction, we find that the overall activity is determined by complex kinetic couplings. Effects of particle shape as well as internal and external strain will be discussed [5].

[1] M. Jørgensen, H. Grönbeck, ACS Catalysis, 6, 6730 (2016).

[2] M. Van den Bossche, H. Grönbeck, J. Am. Chem. Soc. 137, 37 (2015).

[3] M. Van den Bossche, H. Grönbeck, J. Phys. Chem. C 120, 8390 (2017).

[4] M. Jørgensen, H. Grönbeck, ACS Catalysis 7, 5054 (2017).

[5] M. Jørgensen, H. Grönbeck, Angew. Chem. Int Ed. (2018).

5:20 PM HC+SS-WeA-10 Formation of Pd/Ag Sandwiches, a Stable PdAg Subsurface Alloy, and the Pd Segregation induced by CO and O2, Studied with STM, Ambient-pressure XPS, and DFT
Matthijs van Spronsen (Lawrence Berkeley National Laboratory); Kaining Duanmu (UCLA); Robert Madix (Harvard University); Miquel B. Salmeron (Lawrence Berkeley National Laboratory); Philippe Sautet (UCLA); Cynthia Friend (Harvard University)

Efficient chemical production requires the use of materials able to selectively catalyze complex chemical reactions. High selectivity can be obtained by using noble metals. To improve activity, while retaining selectivity, small amounts of an active metal can be added. These species can function as active sites for, e.g., O2 activation.

A promising alloy for selective hydrogenation is AgPd. The challenge using this material as a catalyst is the larger surface free energy of Pd, compared to Ag. This difference leads to Ag enrichment of the surface, hence, more Pd, both costly and scarce, is needed to ensure that Pd is present in the surface.

Here, the stability of Pd was probed in (sub)surface alloys formed by depositing sub-monolayer amounts of Pd on pure Ag single crystals, in order to evaluate the feasibility of using surface alloys as catalysts. These surface alloys would use much less Pd, but depend critically on understanding and control of the alloy composition and stability.

Mild annealing (~400 K) under ultra-high vacuum, leads to extensive migration of Ag and the formation of a Ag capping layer on the Pd. These measurements were based on scanning tunneling microscopy and X-ray photoelectron spectroscopy, supported by modeling (density functional theory). The etching, due to the loss of Ag, and the increase in island height, due to the capping, results in three-layer-high islands, covering both the terraces and the steps of the Ag(111) surface. The absence of significant amounts of Pd in the surface was confirmed by the lack of CO adsorption at 120 K.

The exposure of the Ag/Pd/Ag(111) surface alloy to gas atmospheres, either pure CO or O2 (0.5 & 1 Torr, respectively), reverses the energetics, thus making Pd resurface. This was concluded from large core-level shifts observed in AP XPS measurements. Due to the proximity of Pd to the surface, the required temperature for Pd to segregate was rather low, in fact, for CO it was detected for 300 K. Resurfacing, and subsequent oxidation, of Pd was much slower, requiring cycling the temperature between 300 and 400 K in the O2 atmosphere, several times. The difference stems from the fact that O2 needs to be dissociated before it can form strong bonds to Pd, while CO can adsorb directly to Pd. The initial O2-dissociation activity of the Ag-capped Pd is expected to be rather small.

To conclude, precise understanding of the stability of active metals with a high surface free energy is required to control the composition of surface alloys and its potential as active/selective catalyst.

Session Abstract Book
(294KB, May 6, 2020)
Time Period WeA Sessions | Abstract Timeline | Topic HC Sessions | Time Periods | Topics | AVS2018 Schedule