Rotational effects in the dissociation of H2 on metal surfaces studied by ab initio quantum-dynamics calculations

Six-dimensional quantum dynamics for dissociative chemisorption of H2 and D2 on Ag(111) on a permutation invariant potential energy surface

Quantum dynamics on a permutation invariant potential energy surface for H₂ dissociation on Ag(111) yield satisfactory agreement with experiment.

Dynamics of H2 dissociation on the 1/2 ML c(2 × 2)-Ti/Al(100) surface

The dissociation of H(2) on Ti-covered Al surfaces is relevant to the rehydrogenation and dehydrogenation of the NaAlH(4) hydrogen storage material. The energetically most stable structure for a 1/2 monolayer of Ti deposited on the Al(100) surface has the Ti atoms in the second layer with a c(2 × 2) structure, as has been confirmed by both low-energy electron diffraction and low-energy ion scattering experiments and density functional theory studies. In this work, we investigate the dynamics of H(2) dissociation on a slab model of this Ti/Al(100) surface. Two six-dimensional potential energy surfaces (PESs) have been built for this H(2) + Ti/Al(100) system, based on the density functional theory PW91 and RPBE exchange-correlation functionals. In the PW91 (RPBE) PES, the lowest H(2) dissociation barrier is found to be 0.65 (0.84) eV, with the minimum energy path occurring for H(2) dissociating above the bridge to top sites. Using both PESs, H(2) dissociation probabilities are calculated using the classical trajectory (CT), the quasi-classical trajectory (QCT), and the time-dependent wave-packet methods. We find that the QCT H(2) dissociation probabilities are in good agreement with the quantum dynamics results in the collision energy range studied up to 1.0 eV. We have also performed molecular beam simulations and present predictions for molecular beam experiments. Our molecular beam simulations show that H(2) dissociation on the 1/2 ML Ti/Al(100) surface is an activated process, and the reaction probability is found to be 6.9% for the PW91 functional and 1.8% for the RPBE at a nozzle temperature of 1700 K. Finally, we have also calculated H(2) dissociation rate constants by applying transition state theory and the QCT method, which could be relevant to modeling Ti-catalyzed rehydrogenation and dehydrogenation of NaAlH(4).

SIX-DIMENSIONAL DYNAMICS OF DISSOCIATIVE CHEMISORPTION OF H2 ON METAL SURFACES

The theory of time-dependent quantum dynamics of dissociative chemisorption of hydrogen on metal surfaces is reviewed, in the framework of electronically adiabatic scattering from static surfaces. Four implementations of the time-dependent wave packet (TDWP) method are discussed. In the direct product pseudo-spectral and the spherical harmonics pseudo-spectral methods, no use is made of the symmetry associated with the surface unit cell. This symmetry is exploited by the symmetry adapted wave packet and the symmetry adapted pseudo-spectral (SAPS) method, which are efficient for scattering at normal incidence. The SAPS method can be employed for potential energy surfaces of general form.

Comparison to experiment shows that the TDWP method yields good, but not yet excellent, quantitative accuracy for dissociation of (ν = 0, j = 0) H 2 if the calculations are based on accurately fitted density functional theory calculations that are performed using the generalized gradient approximation. The influence of the molecule's vibration (rotation) is well (reasonably well) described. The theory does not yet yield quantitatively accurate results for rovibrationally inelastic scattering. The theory has helped with the interpretation of existing experimental results, for instance, by solving a parodox regarding the corrugation of Pt(111) as seen by reacting and scattering H 2. The theory has also provided some exciting new predictions, for instance, concerning where on the surface of Cu(100) H2 reacts depending on its vibrational state. Future theoretical studies of H 2 reacting on metal surfaces will likely be aimed at validating GGAs for molecule-surface interactions, and understanding trends in collisions of H 2 with complex metal surfaces.

Ab initio based tight-binding molecular dynamics simulation of the sticking and scattering of O2∕Pt(111)

The sticking and scattering of O(2)Pt(111) has been studied by tight-binding molecular dynamics simulations based on an ab initio potential energy surface. We focus, in particular, on the sticking probability as a function of the angle of incidence and the energy and angular distributions in scattering. Our simulations provide an explanation for the seemingly paradox experimental findings that adsorption experiments suggest that the O(2)Pt(111) interaction potential should be strongly corrugated while scattering experiments indicate a rather small corrugation. The potential energy surface is indeed strongly corrugated which leads to a pronounced dependence of the sticking probability on the angle of incidence. The scattered O(2) molecules, however, experience a rather flat surface due to the fact that they are predominantly scattered at the repulsive tail of the potential.

Off-normal incidence dissociative sticking of H2 on Cu(100) studied using six-dimensional quantum calculations

Six-dimensional quantum calculations of the sticking probability for H2 hitting a Cu(100) surface with off-normal incidence are presented. The multiconfiguration time-dependent Hartree approach is employed for an efficient wave-packet propagation. The sticking probability is calculated for different initial momenta parallel to the surface. In contrast with the picture described in the literature, the sticking probability was found to depend on the parallel momentum. The results are explained by the topology of the potential-energy surface, which shows significant corrugation with a moderate variation of the barrier height with the surface site.

The reaction rate for dissociative adsorption of N2 on stepped Ru(0001): Six-dimensional quantum calculations

Quantum-mechanical calculations of the reaction rate for dissociative adsorption of N2 on stepped Ru(0001) are presented. Converged six-dimensional quantum calculations for this heavy-atom reaction have been performed using the multiconfiguration time-dependent Hartree method. A potential-energy surface for the transition-state region is constructed from density-functional theory calculations using Shepard interpolation. The quantum results are in very good agreement with the results of the harmonic transition-state theory. In contrast to the findings of previous model calculations on similar systems, the tunneling effect is found to be small.

High-dimensional quantum dynamical study of the dissociation of H2 on Pd(110)

We report the first six-dimensional quantum dynamical study of the dissociative adsorption of H(2) on a (110) surface. We have performed quantum coupled-channel calculations for the system H(2)/Pd(110) based on a potential energy surface (PES) that was derived from ab initio electronic structure calculations. In particular, we have focused on the effects of the corrugation and anisotropy of the PES on the H(2) dissociation probability. Our results agree well with the available experimental data for the sticking probability as a function of the initial kinetic energy and the angle of incidence. Because of the coupling between the anisotropy and corrugation of the potential energy surface our calculations predict an unusual rotational heating and a rather small rotational alignment in desorption.