Inelastic gas-surface scattering. I. Formalism

An examination of phonon-inelastic molecule-metal scattering using reduced density matrix and stochastic wave packet methods

We explore the application of reduced density matrix-based approaches to molecules interacting with the lattice vibrations of metals, an interaction responsible for the temperature dependence of many of the fundamental steps of catalysis. We avoid the use of simple models for the bath and instead use density functional theory to compute all molecule-phonon interactions and the properties of the lattice phonons, for methane scattering from Ir(111). We find that while the large metal mass leads to long bath correlation times, these are not significantly longer than the time over which the reduced density matrix changes due to interactions with the bath. We show that the neglect of memory is reasonable and the use of the Redfield equation is justified. We also show how the commonly used rotating wave approximation is far too severe for this scattering problem. A less restrictive approximation that is nearly exact for our system gives an equation of motion in the Lindblad form. As a result, the Monte Carlo wave packet methods can be used to describe gas-phonon scattering, guaranteeing positivity, and with all couplings derived from first-principles.

The trapping of methane on Ir(111): A first-principles quantum study

We implement a fully quantum mechanical study of methane trapping on Ir(111), where the phonons, the molecule-surface interaction, and the molecule-phonon coupling are all computed from first-principles. We find that both the surface corrugation and the phonon coupling vary strongly with molecular orientation and that there is a "chemical" aspect to this due to the catalytic nature of the metal. For example, molecules with reactive orientations can approach close to surface sites with low barriers to dissociation. Moreover, lattice motion can lower the barrier to dissociation, leading to unusual behavior for the phonon coupling. We find good agreement with experiment and two recent classical studies if we average our potential energy surface over several orientations of the molecule. We also find reasonable agreement with a recent study of methane diffraction, although we show that diffraction does not play a major role in trapping on the smooth Ir(111) surface and that trapping obeys normal energy scaling, consistent with experiment. We show that the trapping probability can be sensitive to the temperature at both high and low incidence energies. Relaxation and desorption of trapped particles are examined.

Hot Adatom Diffusion Following Oxygen Dissociation on Pd(100) and Pd(111): A First-Principles Study of the Equilibration Dynamics of Exothermic Surface Reactions

We augment ab initio molecular dynamics simulations with a quantitative account of phononic dissipation to study the hyperthermal adsorbate dynamics resulting from a noninstantaneous energy dissipation during exothermic surface chemical reactions. Comparing the hot adatom diffusion ensuing O_{2} dissociation over Pd(100) and Pd(111) we find experimentally accessible product end distances to form a rather misleading measure for the lifetime of this hyperthermal state. The lifetime is particularly long at Pd(111) where a random-walk-type diffusion leads only to small net displacements. A detailed phonon analysis rationalizes the slow equilibration through long-lived Rayleigh mode excitations that spatially confine the released energy within a nanoscopic "hot spot" around the impingement region.

Sticking of hydrogen on supported and suspended graphene at low temperature

The physisorption of atomic hydrogen on graphene is investigated quantum mechanically using a semiempirical model for the lattice dynamics. A thermally averaged wave packet propagation describes the motion of the H atoms with respect to the membrane. Two graphene configurations, either supported on a silicone oxide substrate or suspended over a hole in the substrate, are considered. In both cases, the phonon spectrum is modified in such a way that graphene is stabilized with respect to thermal fluctuations. The sticking probabilities of hydrogen on these stabilized membranes at 10 K are high at low collision energies, and larger than on graphite.

Sticking and desorption of hydrogen on graphite: a comparative study of different models

We study the physisorption of atomic hydrogen on graphitic surfaces with four different quantum mechanical methods: perturbation and effective Hamiltonian theories, close coupling wavepacket, and reduced density matrix propagation methods. Corrugation is included in the modeling of the surface. Sticking is a fast process which is well described by all methods. Sticking probabilities are of the order of a few percent in the collision energy range 0-25 meV, but are enhanced for collision energies close to those of diffraction resonances. Sticking also increases with surface temperature. Desorption is a slow process which involves multiphonon processes. We show, however, how to correct the close coupling wavepacket method to account for such phenomena and obtain correct time constants for initial state decay. Desorption time constants are in the range of 20-50 ps for a surface temperature of 300 K.

Reduced density matrix quantum approach for particle trapping and sticking on corrugated moving surfaces

A short time propagation algorithm for the reduced density matrix is derived to model the interaction of a quantum particle with a moving corrugated surface. The algorithm includes dissipative terms, which can be derived directly from the full Hamiltonian. The scattering of He from a corrugated Cu surface is examined as a function of incident energy and angle and the temperature of the substrate, with a focus on the nature of trapping. It is found that corrugation can make a significant contribution to trapping, even on a metal surface. Energy exchange with the phonons is shown to significantly modify the nature of diffraction mediated selective adsorption.