Ab initio Quantum and ab initio Molecular Dynamics of the Dissociative Adsorption of Hydrogen on Pd(100)

Abstract

The dissociative adsorption of hydrogen on Pd(100) has been studied by ab initio quantum dynamics and ab initio molecular dynamics calculations. Treating all hydrogen degrees of freedom as dynamical coordinates implies a high dimensionality and requires statistical averages over thousands of trajectories. An efficient and accurate treatment of such extensive statistics is achieved in two steps: In a first step we evaluate the ab initio potential energy surface (PES) and determine an analytical representation. Then, in an independent second step dynamical calculations are performed on the analytical representation of the PES. Thus the dissociation dynamics is investigated without any crucial assumption except for the Born-Oppenheimer approximation which is anyhow employed when density-functional theory calculations are performed. The ab initio molecular dynamics is compared to detailed quantum dynamical calculations on exactly the same ab initio PES. The occurence of quantum oscillations in the sticking probability as a function of kinetic energy is addressed. They turn out to be very sensitive to the symmetry of the initial conditions. At low kinetic energies sticking is dominated by the steering effect which is illustrated using classical trajectories. The steering effects depends on the kinetic energy, but not on the mass of the molecules. Zero-point effects lead to strong differences between quantum and classical calculations of the sticking probability. The dependence of the sticking probability on the angle of incidence is analysed; it is found to be in good agreement with experimental data. The results show that the determination of the potential energy surface combined with high-dimensional dynamical calculations, in which all relevant degrees of freedon are taken into account, leads to a detailed understanding of the dissociation dynamics of hydrogen at a transition metal surface.