Chemical reaction thresholds according to classical-limit quantum dynamics.

Classical-limit quantum dynamics is used to explain the origin of the quantum thresholds of chemical reactions from their classical dynamics when these are vibrationally nonadiabatic across the interaction region. This study is performed within the framework of an elementary model of chemical reaction that mimics the passage from the free rotation of the reagents to the bending vibration at the transition state to the free rotation of the products.

When Classical Trajectories Get to Quantum Accuracy: The Scattering of H2 on Pd(111).

When elementary reactive processes occur at such low energies that only a few states of reactants and/or products are available, quantum effects strongly manifest and the standard description of the dynamics within the classical framework fails. We show here, for H2 scattering on Pd(111), that by pseudoquantizing in the spirit of Bohr the relevant final actions of the system, along with adequately treating the diffraction-mediated trapping of the incoming wave, classical simulations achieve an unprecedented agreement with state-of-the-art quantum dynamics calculations.

Dynamics of H2 Eley-Rideal abstraction from W(110): sensitivity to the representation of the molecule-surface potential.

Dynamics of the Eley-Rideal (ER) abstraction of H2 from W(110) is analyzed by means of quasi-classical trajectory calculations. Simulations are based on two different molecule-surface potential energy surfaces (PES) constructed from Density Functional Theory results. One PES is obtained by fitting, using a Flexible Periodic London-Eyring-Polanyi-Sato (FPLEPS) functional form, and the other by interpolation through the corrugation reducing procedure (CRP). Then, the present study allows us to elucidate the ER dynamics sensitivity on the PES representation. Despite some sizable discrepancies between both H+H/W(110) PESs, the obtained projectile-energy dependence of the total ER cross sections are qualitatively very similar ensuring that the main physical ingredients are captured in both PES models. The obtained distributions of the final energy among the different molecular degrees of freedom barely depend on the PES model, being most likely determined by the reaction exothermicity. Therefore, a reasonably good agreement with the measured final vibrational state distribution is observed in spite of the pressure and material gaps between theoretical and experimental conditions.

Surface temperature effects on the dynamics of N2 Eley-Rideal recombination on W(100).

Quasiclassical trajectories simulations are performed to study the influence of surface temperature on the dynamics of a N atom colliding a N-preadsorbed W(100) surface under normal incidence. A generalized Langevin surface oscillator scheme is used to allow energy transfer between the nitrogen atoms and the surface. The influence of the surface temperature on the N(2) formed molecules via Eley-Rideal recombination is analyzed at T = 300, 800, and 1500 K. Ro-vibrational distributions of the N(2) molecules are only slightly affected by the presence of the thermal bath whereas kinetic energy is rather strongly decreased when going from a static surface model to a moving surface one. In terms of reactivity, the moving surface model leads to an increase of atomic trapping cross section yielding to an increase of the so-called hot atoms population and a decrease of the direct Eley-Rideal cross section. The energy exchange between the surface and the nitrogen atoms is semi-quantitatively interpreted by a simple binary collision model.

Dynamical reaction pathways in Eley-Rideal recombination of nitrogen from W(100).

The scattering of atomic nitrogen over a N-pre-adsorbed W(100) surface is theoretically described in the case of normal incidence off a single adsorbate. Dynamical reaction mechanisms, in particular Eley-Rideal (ER) abstraction, are scrutinized in the 0.1-3.0 eV collision energy range and the influence of temperature on reactivity is considered between 300 and 1500 K. Dynamics simulations suggest that, though non-activated reaction pathways exist, the abstraction process exhibits a significant collision energy threshold (0.5 eV). Such a feature, which has not been reported so far in the literature, is the consequence of a repulsive interaction between the impinging and the pre-adsorbed nitrogens along with a strong attraction towards the tungsten atoms. Above threshold, the cross section for ER reaction is found one order of magnitude lower than the one for hot-atoms formation. The abstraction process involves the collision of the impinging atom with the surface prior to reaction but temperature effects, when modeled via a generalized Langevin oscillator model, do not affect significantly reactivity.

Dynamics simulation of N(2) scattering onto W(100,110) surfaces: A stringent test for the recently developed flexible periodic London-Eyring-Polanyi-Sato potential energy surface.

An efficient method to construct the six dimensional global potential energy surface (PES) for two atoms interacting with a periodic rigid surface, the flexible periodic London-Eyring-Polanyi-Sato model, has been proposed recently. The main advantages of this model, compared to state-of-the-art interpolated ab initio PESs developed in the past, reside in its global nature along with the small number of electronic structure calculations required for its construction. In this work, we investigate to which extent this global representation is able to reproduce the fine details of the scattering dynamics of N(2) onto W(100,110) surfaces reported in previous dynamics simulations based on locally interpolated PESs. The N(2)/W(100) and N(2)/W(110) systems are chosen as benchmarks as they exhibit very unusual and distinct dissociative adsorption dynamics although chemically similar. The reaction pathways as well as the role of dynamic trapping are scrutinized. Besides, elastic/inelastic scattering dynamics including internal state and angular distributions of reflected molecules are also investigated. The results are shown to be in fair agreement with previous theoretical predictions.

DFT study of dissociative adsorption of hydrogen sulfide on Cu(111) and Au(111).

Density functional theory (DFT) is used to investigate the reaction pathways for H2S adsorption on Au(111) and Cu(111) at low coverage as well as the full decomposition of H2S on Cu(111). On both surfaces, a weakly bonded molecular state is found with the S atom bond on top sites being molecular adsorption, a nonactivated process. The H-SH dissociation process is endothermic on Au(111), and all reaction pathways present high activation energy barriers which explains the extremely low dissociation probability of H2S on defect-free Au(111) estimated from experiments. This scenario slightly changes for H2S/Cu(111): (i) dissociated configurations are energetically more favorable than the molecular state and (ii) the H-SH bond cleavage process presents a relatively small activation energy barrier. This is not inconsistent with low but nonzero reactive sticking probability of thermal H2S molecules reported in experiments. The complete energy profile for the H2S adsorption and full decomposition is compatible with the accumulation of S-adatoms observed experimentally.

Modified Shepard interpolation method applied to trapping mediated adsorption dynamics.

The modified Shepard (MS) interpolation method is applied to H(2)/Pd(111) to investigate its performance for a system for which dissociative adsorption takes place through a direct as well as an indirect (i.e. dynamic trapping) mechanism. The input data were obtained from an available accurate potential energy surface (PES) interpolated by using the corrugation reducing procedure (CRP). Dissociation probabilities obtained from classical trajectory calculations with the MS-PES are in very good agreement with the results for the CRP-PES. Thus, this study confirms the MS method as a promising tool to tackle low energy adsorption dynamics of polyatomic molecules, usually dominated by trapping.

Multiconfiguration time-dependent Hartree method applied to molecular dissociation on surfaces: H2 + Pt(111).

Four-dimensional quantum dynamics calculations are performed on the dissociative chemisorption of H(2) on Pt(111) using the multiconfiguration time-dependent Hartree method. The aim of this work is to study the performance of the multiconfiguration time-dependent Hartree method for a gas-surface reaction by comparison with the standard time-dependent wave-packet propagation method. The initial-state resolved dissociation probability of H(2) is calculated within two four-dimensional models. The first four-dimensional model treats explicitly the rotational motion of the molecule and the H(2) dissociation is studied above two different sites of the Pt(111). For this model, only a potential-energy surface of general form was available. This potential was refitted to a sum of product form to allow efficient calculations with the multiconfiguration time-dependent Hartree method. The second model focuses on the description of the center-of-mass motion parallel to the surface, the rotational motion of the molecule being frozen. These four-dimensional quantum dynamics calculations yield important insights which can help with performing full six-dimensional calculations on H(2) dissociating on Pt(111). The multiconfiguration time-dependent Hartree method is shown to be particularly efficient for computing initial-state selective dissociation probabilities for the system studied, with a good accuracy and a reduced amount of memory and computational time when compared to the standard time-dependent wave-packet method.

Application of the modified Shepard interpolation method to the determination of the potential energy surface for a molecule-surface reaction: H2 + Pt(111).

We have used a modified Shepard (MS) interpolation method, initially developed for gas phase reactions, to build a potential energy surface (PES) for studying the dissociative chemisorption of H2 on Pt(111). The aim was to study the efficiency and the accuracy of this interpolation method for an activated multidimensional molecule-surface reactive problem. The strategy used is based on previous applications of the MS method to gas phase reactions, but modified to take into account special features of molecule-surface reactions, like the presence of many similar reaction pathways which vary only slightly with surface site. The efficiency of the interpolation method was tested by using an already existing PES to provide the input data required for the construction of the new PES. The construction of the new PES required half as many ab initio data points as the construction of the old PES, and the comparison of the two PESs shows that the method is able to reproduce with good accuracy the most important features of the H2 + Pt(111) interaction potential. Finally, accuracy tests were done by comparing the results of dynamics simulations using the two different PESs. The good agreement obtained for reaction probabilities and probabilities for rotationally and diffractionally inelastic scattering shows clearly that the MS interpolation method can be used efficiently to yield accurate PESs for activated molecule-surface reactions.