Search

Dynamics of the Eley-Rideal (ER) abstraction of H2 from W(110) is analyzed by means of quasiclassical 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. [ABSTRACT FROM AUTHOR]

An approach based on reactive force fields is applied to the parametrization of potential energy surface (PES) for chemical reactions on surfaces with a benchmark system, H2/Pd(111). We show that a simple reactive force field based on the second moment approximation does not allow for obtaining reliable results of reaction dynamics for the considered system. With a more elaborate reactive force field, i.e., reactive bond order (REBO) force field, we succeeded in obtaining a reliable PES for H2/Pd(111). The accuracy of the constructed REBO force field is carefully checked through various tests including the comparison not only between energies calculated with density functional theory and those with REBO force field but also between the available results of ab initio molecular dynamics simulations and those with our force field. Moreover, our REBO force field is endowed with some transferability since the force field constructed with a database containing only information on H2/Pd(111) allows for obtaining also accurate results for H2/Pd(100) and qualitatively correct results for H2/Pd(110) without any refitting. With the help of our reactive force field, the molecular dynamics simulation for the dissociation of H2 on the considered Pd surfaces is speeded up by five orders of magnitude compared to ab initio molecular dynamics method. The demonstrated reliability and the very high computational efficiency of reactive force fields open extremely attractive perspectives for studying large-scale complex reacting systems. [ABSTRACT FROM AUTHOR]

We study the dissociative adsorption of N2 on W(100) and W(110) by means of density functional theory and classical dynamics. Working with a full six-dimensional adiabatic potential energy surface (PES), we find that the theoretical results of the dynamical problem strongly depend on the choice of approximate exchange-correlation functional for the determination of the PES. We consider the Perdew-Wang-91 [Perdew et al., Phys. Rev. B 46, 6671 (1992)] and Perdew-Burke-Ernzerhof (RPBE) [Hammer et al., Phys. Rev. B 59, 7413 (1999)] functionals and carry out a systematic comparison between the dynamics determined by the respective PESs. Even though it has been shown in earlier works that the RPBE may provide better values for the chemisorption energies, our study brings evidence that it gives rise to a PES with excessive repulsion far from the surface. [ABSTRACT FROM AUTHOR]

We have performed H2(D2) diffraction experiments on a Pd(110) surface using two different high-sensitivity set-ups. We have found that, although the total reflectivity of Pd(110) is comparable to that observed in other reactive systems, the corresponding H2(D2) diffraction patterns are quite different: no diffraction peak, including the specular one, is observed on Pd(110). This unexpected result is the consequence of dynamic trapping. Such interpretation is supported by classical dynamics calculations based on accurate ab initio potential energy surfaces. [ABSTRACT FROM AUTHOR]

We present quasiclassical dynamics calculations of H2 and D2 scattering by the NiAl(110) surface using a recently proposed six-dimensional potential-energy surface (PES) obtained from density-functional theory calculations. The results for dissociative adsorption confirm several experimental predictions using (rotationally hot) D2 beams, namely, the existence of a dissociation barrier, the small isotopic effect, the importance of vibrational enhancement, and the existence of normal energy scaling. The latter conclusion shows that normal energy scaling is not necessarily associated with weak corrugated surfaces. The results for rotationally elastic and inelastic diffractions are also in reasonable agreement with experiment, but they show that many more diffractive transitions are responsible for the observed structures than previously assumed. This points to the validity of the PES recently proposed [P. Rivière, H. F. Busnengo, and F. Martín, J. Chem. Phys. 121, 751 (2004)] to describe dissociative adsorption as well as rotationally elastic and inelastic diffractions in the H2/NiAl(110) system. [ABSTRACT FROM AUTHOR]

We present results of extensive density functional theory (DFT) calculations for H and H2 interacting with NiAl(110). Continuous representations of the full dimensional potential energy surface (PES) for the H/NiAl(110) and H2/NiAl(110) systems are obtained by interpolation of the DFT results using the corrugation reducing procedure. We find a minimum activation energy barrier of ∼300 meV for dissociative adsorption of H2, which is consistent with the energy threshold obtained in molecular beam experiments for H2 (ν=0). We explain vibrational enhancement observed in experiments as the consequence of vibrational softening in the entrance channel over the most reactive surface site. The H2/NiAl(110) PES shows a high surface site selectivity: for energies up to 0.1 eV above threshold, H2 adsorption can only take place around top-Ni sites (within a circle of radius ∼0.3 Å). A strong energetic corrugation is observed: energy barriers for dissociation vary by more than 1 eV between the most and the least reactive sites. In contrast, geometric corrugation is much less pronounced and comparable to that of low index single metal surfaces like Cu or Pt. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]