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

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.

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.