State selective vibrational (2ν3) activation of the chemisorption of methane on Pt (111)

Accurate Reaction Probabilities for Translational Energies on Both Sides of the Barrier of Dissociative Chemisorption on Metal Surfaces

Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule-metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.

A versatile and narrow linewidth infrared radiation source for ro-vibration state-selected preparation of molecules in molecular beams

We describe the design and characterization of a versatile pulsed (5 ns, 10 Hz repetition rate) optical parametric oscillator and amplifier system capable of generating single longitudinal mode, narrow linewidth (0.01 cm-1) radiations in the wavelength range of 680-870 nm and 1380-4650 nm. Using a combination of power-normalized photoacoustic signal and a Fizeau interferometer-based wavemeter, we are able to actively stabilize the output wavenumber to within 0.005 cm-1 (3σ) over a timescale longer than 1000 s. We demonstrate an application of this system by performing ro-vibration state-selected preparation of CO in the v = 2 state, via direct overtone excitation (v = 0 → 2 at 2346 nm) and subsequent state-selected detection in an internally cold molecular beam.

Mode Selective Chemistry for the Dissociation of Methane on Efficient Ni/Pt-Bimetallic Alloy Catalysts

The mode selectivity of methane dissociation is studied on three different Ni/Pt-bimetallic alloy surfaces using a fully quantum approach based on reaction path Hamiltonian. Dissociative sticking probability depends on the...

Methane activation on single-atom Ir-doped metal nanoparticles from first principles

The breaking of the C-H bond of CH4 is of great importance, and one of the most efficient strategies in heterogeneous catalysis is to alter the electronic structure of a surface by doping it with different metal elements or controlling the stoichiometry. We present an in-depth study on methane activation on pure metal and single-atom Ir-doped alloy nanoparticles, which are constructed based on (100), (110), (111) surfaces using density functional theory (DFT) calculations. DFT results show that the dissociation barriers of CH4 on the Ir-doped alloy surfaces are about 0.3-0.4 eV, much lower than those on the pure metal surfaces (i.e., 0.6-0.8 eV). DFT-based transition state theory further reveals the rates of the first C-H activation on single-atom Ir-doped alloy nanoparticles at the early stages. Importantly, a strong temperature dependence is mainly contributed by the proportion of the exposed (110) surface. The Ir-doped Pt nanoparticle is found to be an efficient catalyst for methane activation in potential industrial applications. These important results are helpful for further designing new metal catalysts for methane activation at the atomic/molecular level.

Plasma-enhanced catalysis for the upgrading of methane: a review of modelling and simulation methods

Modelling methods and simulation works on the upgrading of methane via plasma and plasma-enhanced catalysis reviewed.

Dynamical Study of the Dissociative Chemisorption of CHD3 on Pd(111)

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 CHD₃ on Pd(111) using Born-Oppenheimer molecular dynamics while also performing a rough comparison with experimental data for CH₄ + 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 CHD₃ 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.

Dynamics in reactions on metal surfaces: A theoretical perspective

Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.

Low-energy methane scattering from Pt(111)

We have measured the temperature dependence of angular distributions of CH₄ from Pt(111) at an incident energy of 109 meV. A broad angular distribution has been observed along the two main symmetry directions, whereby the peak center shifts from the supra-specular position to the sub-specular position when the surface temperature increases from 120 K to 800 K. Different widths have been measured for the scattering patterns along the [ 1¯01 ] and the [ 2¯11 ] azimuthal directions. Based on calculations performed within the binary collision model, these differences have been ascribed to different corrugations of the CH₄-Pt(111) interaction potential along the two high-symmetry directions. This corrugation has been estimated from the model calculations to amount ∼0.03 Å, a factor of three larger than the one measured with helium diffraction.

Coherent quantum scattering of CH4 from Ni(111)

We have measured high-resolution angular distributions of methane scattered from a Ni(111) surface at incident energies between 68.9 meV and 108.6 meV. A sharp and intense specular peak has been observed, in addition to sharp features corresponding to rotationally inelastic diffraction (RID) peaks along the two main symmetry directions of Ni(111). The intensity of the most intense RID peaks is ca. 50% of the specular one. The observation of sharp, coherent elastic peaks at such low incident energies suggests that single scattering dominates over trapping at these energies, and that the depth of the van der Waals well should be lower than 60 meV. In contrast, a broad angular distribution shifted from the specular position is observed from a graphene-covered Ni(111) surface under identical incident conditions. These results open up the possibility of studying the physisorption well between CH₄ and a transition metal surface using high-resolution molecular beams.

Quantum dynamics of polyatomic dissociative chemisorption on transition metal surfaces: mode specificity and bond selectivity

Recent advances in quantum dynamical characterization of polyatomic dissociative chemisorption on accurate global potential energy surfaces are critically reviewed.

Ab Initio Molecular Dynamics Calculations versus Quantum-State-Resolved Experiments on CHD3 + Pt(111): New Insights into a Prototypical Gas-Surface Reaction

The dissociative chemisorption of methane on metal surfaces is of fundamental and practical interest, being a rate-limiting step in the steam reforming process. The reaction is best modeled with quantum dynamics calculations, but these are currently not guaranteed to produce accurate results because they rely on potential energy surfaces based on untested density functionals and on untested dynamical approximations. To help overcome these limitations, here we present for the first time statistically accurate reaction probabilities obtained with ab initio molecular dynamics (AIMD) for a polyatomic gas-phase molecule reacting with a metal surface. Using a general purpose density functional, the AIMD reaction probabilities are in semiquantitative agreement with new quantum-state-resolved experiments on CHD3 + Pt(111). The comparison suggests the use of the sudden approximation for treating the rotations even though CHD3 has large rotational constants and yields an estimated reaction barrier of 0.9 eV for CH4 + Pt(111).

An effusive molecular beam technique for studies of polyatomic gas-surface reactivity and energy transfer

An effusive molecular beam technique is described to measure alkane dissociative sticking coefficients, S(T(g), T(s); ϑ), on metal surfaces for which the impinging gas temperature, T(g), and surface temperature, T(s), can be independently varied, along with the angle of incidence, ϑ, of the impinging gas. Effusive beam experiments with T(g) = T(s) = T allow for determination of angle-resolved dissociative sticking coefficients, S(T; ϑ), which when averaged over the cos (ϑ)/π angular distribution appropriate to the impinging flux from a thermal ambient gas yield the thermal dissociative sticking coefficient, S(T). Nonequilibrium S(T(g), T(s); ϑ) measurements for which T(g) ≠ T(s) provide additional opportunities to characterize the transition state and gas-surface energy transfer at reactive energies. A resistively heated effusive molecular beam doser controls the T(g) of the impinging gas striking the surface. The flux of molecules striking the surface from the effusive beam is determined from knowledge of the dosing geometry, chamber pressure, and pumping speed. Separate experiments with a calibrated leak serve to fix the chamber pumping speed. Postdosing Auger electron spectroscopy is used to measure the carbon of the alkyl radical reaction product that is deposited on the surface as a result of alkane dissociative sticking. As implemented in a typical ultrahigh vacuum chamber for surface analysis, the technique has provided access to a dynamic range of roughly 6 orders of magnitude in the initial dissociative sticking coefficient for small alkanes on Pt(111).

Methane dissociation on Ni(111): a new understanding of the lattice effect

The nature of the lattice motion during the dissociation of methane on Ni(111) is analyzed in great detail, and various models for including lattice effects are explored. It is shown that the thermal vibrations of the lattice strongly modify the reactivity, but that the lattice motion is relatively unperturbed by the incident molecule during the collision, in contrast with several earlier predictions. Based on these studies we propose a new model to describe the effects of lattice motion, which agrees well with exact quantum calculations.

Vibrational vs. translational energy in promoting a prototype metal-hydrocarbon insertion reaction

The reaction Y + CH(4) --> HYCH(3) --> YCH(2) + H(2) is initiated by C-H insertion involving a 20 +/- 3 kcal/mol potential energy barrier. The reaction is studied in crossed molecular beams under two different conditions with nearly the same total energy. One experiment is carried out at a collision energy of 15.1 kcal/mol with one quantum of CH(4) antisymmetric (nu(3)) stretching vibrational excitation (8.63 kcal/mol), the other at a collision energy of 23.8 kcal/mol. The reaction cross-section for C-H stretch excited methane (sigma(s)) is found to be at least a factor of 2.2 times larger than for ground-state methane (sigma(g)) at the same total energy.

Chemical dynamics of vibrationally excited molecules: Controlling reactions in gases and on surfaces

Experimental studies of the chemical reaction dynamics of vibrationally excited molecules reveal the ability of different vibrations to control the course of a reaction. This Perspective describes those studies for the prototypical reaction of vibrationally excited methane and its isotopologues in gases and on surfaces and looks to the prospects of similar studies in liquids. The influences of vibrational excitation on the C-H bond cleavage in a single collision reaction with Cl and in dissociative adsorption on a Ni surface bear some striking similarities. Both reactions are bond-selective processes in which the initial preparation of a molecular eigenstate containing a large component of C-H stretching results in preferential cleavage of that bond. It is possible to cleave either the C-H bond or C-D bond in the reaction of Cl with CH3D, CH2D2, or CHD3 and, similarly, to use initial excitation of the C-H stretch to promote dissociation of CHD3 to CD3 and H on a Ni surface. Different vibrational modes, such as the symmetric and antisymmetric stretches in CH3D or CH4, lead to very different reactivities, and molecules with the symmetric stretching vibration excited can be as much as 10 times more reactive than ones with the antisymmetric stretch excited. The origin of this behavior lies in the change in the vibrational motion induced by the interaction with the atomic reaction partner or the surface.

Methane dissociative adsorption on the Pt(111) surface over the 300-500 K temperature and 1-10 Torr pressure ranges

The dissociative adsorption of methane on the Pt(111) surface has been investigated and characterized over the 1-10 Torr pressure and 300-500 K temperature ranges using sum frequency generation (SFG) vibrational spectroscopy and Auger electron spectroscopy (AES). At a reaction temperature of 300 K and a pressure of 1 Torr, C-H bond dissociation occurs in methane on the Pt(111) surface to produce adsorbed methyl (CH(3)) groups, carbon, and hydrogen. SFG results suggest that C-C coupling occurs at higher reaction temperatures and pressures. At 400 K, methyl groups react with adsorbed C to form ethylidyne (C(2)H(3)), which dehydrogenates at 500 K to form ethynyl (C(2)H) and methylidyne (CH) species, as shown by SFG. By 600 K, all of the ethylidyne has reacted to form the dissociation products ethynyl and methylidyne. Calculated C-H bond dissociation probabilities for methane, determined by carbon deposition measured by AES, are in the 10(-8) range and increase with increasing reaction temperature. A mechanism has been developed and is compared with conclusions from other experimental and theoretical studies using single crystals.

Using Effusive Molecular Beams and Microcanonical Unimolecular Rate Theory to Characterize CH4 Dissociation on Pt(111)

The dissociative sticking coefficient for CH4 on Pt(111) has been measured as a function of both gas temperature (Tg) and surface temperature (Ts) using effusive molecular beam and angle-integrated ambient gas dosing methods. The experimental results are used to optimize the three parameters of a microcanonical unimolecular rate theory (MURT) model of the reactive system. The MURT calculations allow us to extract transition state properties from the data as well as to compare our data directly to other molecular beam and thermal equilibrium sticking measurements. We find a threshold energy for dissociation of E0 = 52.5 +/- 3.5 kJ mol(-1). Furthermore, the MURT with an optimized parameter set provides for a predictive understanding of the kinetics of this C-H bond activation reaction, that is, it allows us to predict the dissociative sticking coefficient of CH4 on Pt(111) for any combination of Ts and Tg even if the two are not equal to one another, indeed, the distribution of molecular energy need not even be thermal. Comparison of our results to those from recent thermal equilibrium catalysis studies on CH4 reforming over Pt nanoclusters ( approximately 2 nm diam) dispersed on oxide substrates indicates that the reactivity of Pt(111) exceeds that of the Pt nanocatalysts by several orders of magnitude.

Effusive Molecular Beam Study of C2H6 Dissociation on Pt(111)

The dissociative sticking coefficient for C2H6 on Pt(111) has been measured as a function of both gas temperature (Tg) and surface temperature (Ts) using effusive molecular beam and angle-integrated ambient gas dosing methods. A microcanonical unimolecular rate theory (MURT) model of the reactive system is used to extract transition state properties from the data as well as to compare our data directly with supersonic molecular beam and thermal equilibrium sticking measurements. We report for the first time the threshold energy for dissociation, E0 = 26.5 +/- 3 kJ mol(-1). This value is only weakly dependent on the other two parameters of the model. A strong surface temperature dependence in the initial sticking coefficient is observed; however, the relatively weak dependence on gas temperature indicates some combination of the following (i) not all molecular excitations are contributing equally to the enhancement of sticking, (ii) that strong entropic effects in the dissociative transition state are leading to unusually high vibrational frequencies in the transition state, and (iii) energy transfer from gas-phase rovibrational modes to the surface is surprisingly efficient. In other words, it appears that vibrational mode-specific behavior and/or molecular rotations may play stronger roles in the dissociative adsorption of C2H6 than they do for CH4. The MURT with an optimized parameter set provides for a predictive understanding of the kinetics of this C-H bond activation reaction, that is, it allows us to predict the dissociative sticking coefficient of C2H6 on Pt(111) for any combination of Ts and Tg even if the two are not equal to one another.

Microcanonical unimolecular rate theory at surfaces. III. Thermal dissociative chemisorption of methane on Pt(111) and detailed balance

A local hot spot model of gas-surface reactivity is used to investigate the state-resolved dynamics of methane dissociative chemisorption on Pt(111) under thermal equilibrium conditions. Three Pt surface oscillators, and the molecular vibrations, rotations, and the translational energy directed along the surface normal are treated as active degrees of freedom in the 16-dimensional microcanonical kinetics. Several energy transfer models for coupling a local hot spot to the surrounding substrate are developed and evaluated within the context of a master equation kinetics approach. Bounds on the thermal dissociative sticking coefficient based on limiting energy transfer models are derived. The three-parameter physisorbed complex microcanonical unimolecular rate theory (PC-MURT) is shown to closely approximate the thermal sticking under any realistic energy transfer model. Assuming an apparent threshold energy for CH(4) dissociative chemisorption of E(0)=0.61 eV on clean Pt(111), the PC-MURT is used to predict angle-resolved yield, translational, vibrational, and rotational distributions for the reactive methane flux at thermal equilibrium at 500 K. By detailed balance, these same distributions should be observed for the methane product from methyl radical hydrogenation at 500 K in the zero coverage limit if the methyl radicals are not subject to side reactions. Given that methyl radical hydrogenation can only be experimentally observed when the CH(3) radicals are kinetically stabilized against decomposition by coadsorbed H, the PC-MURT was used to evaluate E(0) in the high coverage limit. A high coverage value of E(0)=2.3 eV adequately reproduced the experimentally observed methane angular and translational energy distributions from thermal hydrogenation of methyl radicals. Although rigorous application of detailed balance arguments to this reactive system cannot be made because thermal decomposition of the methyl radicals competes with hydrogenation, approximate applicability of detailed balance would argue for a strong coverage dependence of E(0) with H coverage--a dependence not seen for methyl radical hydrogenation on Ru(0001), but not yet experimentally explored on Pt(111).