Quantum dynamics of dissociative chemisorption of CH4 on Ni(111): Influence of the bending vibration

Abstract

Two-dimensional, three-dimensional, and four-dimensional quantum dynamic calculations are performed on the dissociative chemisorption of CH(4) on Ni(111) using the multiconfiguration time-dependent Hartree (MCTDH) method. The potential energy surface used for these calculations is 15-dimensional (15D) and was obtained with density functional theory for points which are concentrated in the region that is dynamically relevant to reaction. Many reduced dimensionality calculations were already performed on this system, but the molecule was generally treated as pseudodiatomic. The main improvement of our model is that we try to describe CH(4) as a polyatomic molecule by including a degree of freedom describing a bending vibration in our three-dimensional and four-dimensional models. Using a polyspherical coordinate system, a general expression for the 15D kinetic energy operator is derived, which discards all the singularities in the operator and includes rotational and Coriolis coupling. We use seven rigid constraints to fix the CH(3) umbrella of the molecule to its gas phase equilibrium geometry and to derive two-dimensional, three-dimensional, and four-dimensional Hamiltonians, which were used in the MCTDH method. Only four degrees of freedom evolve strongly along the 15D minimum energy path: the distance of the center of mass of the molecule to the surface, the dissociative C[Single Bond]H bond distance, the polar orientation of the molecule, and the bending angle between the dissociative C[Single Bond]H bond and the umbrella. A selection of these coordinates is included in each of our models. The polar rotation is found to be important in determining the mode selective behavior of the reaction. Furthermore, our calculations are in good agreement with the finding of Xiang et al. [J. Chem. Phys. 117, 7698 (2002)] in their reduced dimensional calculation that the helicopter motion of the umbrella symmetry axis is less efficient than its cartwheel motion for promoting the reaction. The effect of pre-exciting the bend modes is qualitatively incorrect at higher energies, suggesting the necessity of including additional rotational and vibrational degrees of freedom in the model.