1
|
Wei F, Lin S, Guo H. Direct or Precursor-Mediated? Mechanisms for Methane Dissociation on Pt(110)-(2 × 1) at Both Low and High Incidence Energies. JACS AU 2023; 3:2835-2843. [PMID: 37885592 PMCID: PMC10598834 DOI: 10.1021/jacsau.3c00387] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Revised: 08/24/2023] [Accepted: 08/25/2023] [Indexed: 10/28/2023]
Abstract
The activation of alkanes on metal catalysts may involve a precursor-mediated mechanism, in which impinging molecules are first trapped on the catalyst surface to form an adsorbed precursor and may undergo extensive excursion on the surface in search of an active site. A characteristic feature of such a mechanism is an increasing initial sticking probability (S0) with decreasing incidence energy at low incidence energies. Indeed, such "negative activation" was observed on the reconstructed Pt(110)-(2 × 1) surface with a missing row structure. In this paper, we describe an extensive theoretical investigation of methane dissociation on Pt(110)-(2 × 1) using a machine-learned high-dimensional potential energy surface (PES) based on a first-principles training data set. Quasi-classical trajectories (QCTs) are calculated on the PES to simulate the dissociation of both CH4 and CHD3 at various incidence energies. The agreement with the measured initial sticking probabilities is shown to be substantially improved for high incidence energies when compared to previous theoretical studies, indicating a better characterization of the dissociation barrier. Additional QCT calculations have been carried out for the trapping and diffusion of CHD3 under experimental conditions at low incidence energies. The trapping probability is shown to increase with decreasing incidence energy, consistent with the experimentally observed "negative activation" below 10 kJ/mol. The reactivity of the trapped methane is attributed to the combined effect of its nonthermal diffusion across the surface Pt rows and the lowered barrier reached by surface thermal fluctuation. These simulations shed valuable light on the microscopic dynamics of the initial and often rate-limiting step in heterogeneous catalytic processes involving alkanes.
Collapse
Affiliation(s)
- Fenfei Wei
- State
Key Laboratory of Photocatalysis on Energy and Environment, College
of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Sen Lin
- State
Key Laboratory of Photocatalysis on Energy and Environment, College
of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University
of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
2
|
Meng Q, Chen J, Ma J, Zhang X, Chen J. Adiabatic models for the quantum dynamics of surface scattering with lattice effects. Phys Chem Chem Phys 2022; 24:16415-16436. [PMID: 35766107 DOI: 10.1039/d2cp01560a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Junbo Chen
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China. .,Xi'an Modern Chemistry Research Institute, China North Industries Group Corp., Ltd., East Zhangba Road 168, 710065 Xi'an, China
| | - Jianxing Ma
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Jun Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Yangqiao Road West 155, 350002 Fuzhou, China.,Fujian Science & Technology Innovation Laboratory for Optoelectronic Information of China, Optoelectronic Industry Base at High-tech Zone, 350108 Fuzhou, China
| |
Collapse
|
3
|
Biliškov N. Infrared spectroscopic monitoring of solid-state processes. Phys Chem Chem Phys 2022; 24:19073-19120. [DOI: 10.1039/d2cp01458k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We put a spotlight on IR spectroscopic investigations in materials science by providing a critical insight into the state of the art, covering both fundamental aspects, examples of its utilisation, and current challenges and perspectives focusing on the solid state.
Collapse
Affiliation(s)
- Nikola Biliškov
- Rudjer Bošković Institute, Bijenička c. 54, 10000 Zagreb, Croatia
- Department of Chemistry, McGill University, 801 Sherbrooke St. West, Montreal, QC, H3A 0B8, Canada
| |
Collapse
|
4
|
Wei F, Smeets EWF, Voss J, Kroes GJ, Lin S, Guo H. Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Fenfei Wei
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Egidius W. F. Smeets
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Johannes Voss
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park CA 94025, USA
| | - Geert-Jan Kroes
- Leiden University, Leiden Institute of Chemistry, Einsteinweg 55, 2333 CC, Leiden, The Netherlands
| | - Sen Lin
- State Key Laboratory of Photocatalysis on Energy and Environment, College of Chemistry, Fuzhou University, Fuzhou 350002, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
5
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (Eb, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate Eb can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for Eb for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH4. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of Eb for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate Eb for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O2, H2O, NH3, CO2, and CH3OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
6
|
Fallaque JG, Ramos M, Busnengo HF, Martín F, Díaz C. Normal and off-normal incidence dissociative dynamics of O 2(v,J) on ultrathin Cu films grown on Ru(0001). Phys Chem Chem Phys 2021; 23:7768-7776. [PMID: 33000830 DOI: 10.1039/d0cp03979a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dissociative adsorption of molecular oxygen on metal surfaces has long been controversial, mostly due to the spin-triplet nature of its ground state, to possible non-adiabatic effects, such as an abrupt charge transfer from the metal to the molecule, or even to the role played by the surface electronic state. Here, we have studied the dissociative adsorption of O2 on CuML/Ru(0001) at normal and off-normal incidence, from thermal to super-thermal energies, using quasi-classical dynamics, in the framework of the generalized Langevin oscillator model, and density functional theory based on a multidimensional potential energy surface. Our simulations reveal a rather intriguing behavior of dissociative adsorption probabilities, which exhibit normal energy scaling at incidence energies below the reaction barriers and total energy scaling above, irrespective of the reaction channel, either direct dissociation, trapping dissociation, or molecular adsorption. We directly compare our results with existing scanning tunneling spectroscopy and microscopy measurements. From this comparison, we infer that the observed experimental behavior at thermal energies may be due to ligand and strain effects, as already found for super-thermal incidence energies.
Collapse
Affiliation(s)
- J G Fallaque
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | | | | | | | | |
Collapse
|
7
|
Auras SV, van Lent R, Bashlakov D, Piñeiros Bastidas JM, Roorda T, Spierenburg R, Juurlink LBF. Scaling Platinum-Catalyzed Hydrogen Dissociation on Corrugated Surfaces. Angew Chem Int Ed Engl 2020; 59:20973-20979. [PMID: 32749736 PMCID: PMC7692953 DOI: 10.1002/anie.202005616] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 07/17/2020] [Indexed: 11/07/2022]
Abstract
We determine absolute reactivities for dissociation at low coordinated Pt sites. Two curved Pt(111) single-crystal surfaces allow us to probe either straight or highly kinked step edges with molecules impinging at a low impact energy. A model extracts the average reactivity of inner and outer kink atoms, which is compared to the reactivity of straight A- and B-type steps. Local surface coordination numbers do not adequately capture reactivity trends for H2 dissociation. We utilize the increase of reactivity with step density to determine the area over which a step causes increased dissociation. This step-type specific reactive area extends beyond the step edge onto the (111) terrace. It defines the reaction cross-section for H2 dissociation at the step, bypassing assumptions about contributions of individual types of surface atoms. Our results stress the non-local nature of H2 interaction with a surface and provide insight into reactivity differences for nearly identical step sites.
Collapse
Affiliation(s)
- Sabine V Auras
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Richard van Lent
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Dima Bashlakov
- ILTPE, National Academy of Sciences of Ukraine, 47 Nauky Ave., Kharkiv, 61103, Ukraine
| | | | - Tycho Roorda
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Rick Spierenburg
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Ludo B F Juurlink
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| |
Collapse
|
8
|
Auras SV, Lent R, Bashlakov D, Piñeiros Bastidas JM, Roorda T, Spierenburg R, Juurlink LBF. Scaling Platinum‐Catalyzed Hydrogen Dissociation on Corrugated Surfaces. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202005616] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Sabine V. Auras
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Richard Lent
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Dima Bashlakov
- ILTPE, National Academy of Sciences of Ukraine 47 Nauky Ave. Kharkiv 61103 Ukraine
| | | | - Tycho Roorda
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Rick Spierenburg
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - Ludo B. F. Juurlink
- Leiden Institute of Chemistry Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| |
Collapse
|
9
|
Abstract
AbstractUnderstanding heterogeneous catalysis at the atomic level requires detailed knowledge of the reactivity of different surface sites toward specific bond breaking and bond making events. We illustrate a new method in such investigations. We use a macroscopically curved Pt single crystal containing a large variation in density of highly kinked steps of two different chiralities. Scanning tunneling microscopy maps the entire range of surface structures present on the 31° section surrounding the Pt(111) apex. Whereas most of the surface shows the expected characteristic arrays of parallel steps, hexagonally-shaped, single-atom deep pits remain after cleaning procedures near the apex. Their orientation is indicative of the different chiralities present on the two sides of the crystal’s apex. These unintended defects locally raise the surface defect concentration, but are of little consequence to subsequent reactivity measurements for $$\text {D}_2$$
D
2
dissociation and H–D exchange as probed by supersonic molecular beam techniques. We quantify absolute elementary dissociation and relative isotopic exchange rates across the surface with high spatial resolution. At low incident energies, elementary dissociation of the homonuclear isotoplogues is dominated by the kinked steps. H–D exchange kinetics depend also mostly linearly on step density. The changing ratio of D2 dissociation to H–D formation, however, suggests that anisotropic diffusion of H(D) atoms is of influence to the measured HD production rate.
Collapse
|
10
|
Gutiérrez-González A, Beck RD. Quantum state and surface-site-resolved studies of methane chemisorption by vibrational spectroscopies. Phys Chem Chem Phys 2020; 22:17448-17459. [DOI: 10.1039/d0cp03134h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Infrared spectroscopic methods enable quantum-state-specific and surface-site-selective studies of methane chemisorption on stepped platinum surfaces.
Collapse
Affiliation(s)
- Ana Gutiérrez-González
- Laboratoire de Chimie Physique Moléculaire (LCPM)
- École Polytechnique Fédérale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| | - Rainer D. Beck
- Laboratoire de Chimie Physique Moléculaire (LCPM)
- École Polytechnique Fédérale de Lausanne (EPFL)
- CH-1015 Lausanne
- Switzerland
| |
Collapse
|
11
|
Gerrits N, Chadwick H, Kroes GJ. Dynamical Study of the Dissociative Chemisorption of CHD 3 on Pd(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2019; 123:24013-24023. [PMID: 31602282 PMCID: PMC6778984 DOI: 10.1021/acs.jpcc.9b05757] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Revised: 08/20/2019] [Indexed: 06/10/2023]
Abstract
The specific reaction parameter (SRP) approach to density functional theory has been shown to model reactions of polyatomic molecules with metal surfaces important for heterogeneous catalysis in the industry with chemical accuracy. However, transferability of the SRP functional among systems in which methane interacts with group 10 metals remains unclear for methane + Pd(111). Therefore, in this work, predictions have been made for the reaction of CHD3 on Pd(111) using Born-Oppenheimer molecular dynamics while also performing a rough comparison with experimental data for CH4 + Pd(111) obtained for lower incidence energies. Hopefully, future experiments can test the transferability of the SRP functional among group 10 metals also for Pd(111). We found that the reactivity of CHD3 on Pd(111) is intermediate between and similar to either Pt(111) or Ni(111), depending on the incidence energy and the initial vibrational state distribution. This is surprising because the barrier height and experiments performed at lower incidence energies than investigated here suggest that the reactivity of Pd(111) should be similar to that of Pt(111) only. The relative decrease in the reactivity of Pd(111) at high incidence energies is attributed to site specificity of the reaction and to dynamical effects such as the bobsled effect and energy transfer from methane to the surface.
Collapse
Affiliation(s)
- Nick Gerrits
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, the Netherlands
| | - Helen Chadwick
- Department
of Chemistry, Swansea University, Singleton Park, Swansea SA2 8PP, U.K.
| | - Geert-Jan Kroes
- Gorlaeus
Laboratories, Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, the Netherlands
| |
Collapse
|
12
|
Seminara GN, Peludhero IF, Dong W, Martínez AE, Busnengo HF. Molecular Dynamics Study of Molecular and Dissociative Adsorption Using System-Specific Force Fields Based on Ab Initio Calculations: CO/Cu(110) and $$\text {CH}_4/\text {Pt(110)}$$. Top Catal 2019. [DOI: 10.1007/s11244-019-01196-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
13
|
|