Accurate description of the quantum dynamical surface temperature effects on the dissociative chemisorption of H from Cu(111).

Accurately describing surface temperature effects for the dissociative scattering of H on a metal surface on a quantum dynamical (QD) level is currently one of the open challenges for theoretical surface scientists. We present the first QD simulations of hydrogen dissociating on a Cu(111) surface, which accurately describe all relevant surface temperature effects, using the static corrugation model. The reaction probabilities we obtain show very good agreement with those found using quasi-classical dynamics (QCD), both for individual surface slabs and for an averaged, thus Monte Carlo sampled, set of thermally distorted surface configurations. Rovibrationally elastic scattering probabilities show a much clearer difference between the QCD and QD results, which appears to be traceable back toward thermally distorted surface configurations with very low dissociation probabilities and underlines the importance of investigating more observables than just dissociation. By reducing the number of distorted surface atoms included in the dynamical model, we also show that only including one surface atom, or even three surface atoms, is generally not enough to accurately describe the effects of the surface temperature on dissociation and elastic scattering. These results are a major step forward in accurately describing hydrogen scattering from a thermally excited Cu(111) surface and open up a pathway to better describe reaction and scattering from other relevant crystal facets, such as stepped surfaces, at moderately elevated surface temperatures where quantum effects are expected to play a more important role in the dissociation of H on Cu.