Effect of Surface Modes on the Six-Dimensional Molecule–Surface Scattering Dynamics of H2–Cu(100) and D2–Cu(111) Systems

On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K

We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.

Effect of surface temperature on quantum dynamics of H2 on Cu(111) using a chemically accurate potential energy surface

The effect of surface atom vibrations on H₂ scattering from a Cu(111) surface at different temperatures is being investigated for hydrogen molecules in their rovibrational ground state (v = 0, j = 0). We assume weakly correlated interactions between molecular degrees of freedom and surface modes through a Hartree product type wavefunction. While constructing the six-dimensional effective Hamiltonian, we employ (a) a chemically accurate potential energy surface according to the static corrugation model [M. Wijzenbroek and M. F. Somers, J. Chem. Phys. 137, 054703 (2012)]; (b) normal mode frequencies and displacement vectors calculated with different surface atom interaction potentials within a cluster approximation; and (c) initial state distributions for the vibrational modes according to Bose-Einstein probability factors. We carry out 6D quantum dynamics with the so-constructed effective Hamiltonian and analyze sticking and state-to-state scattering probabilities. The surface atom vibrations affect the chemisorption dynamics. The results show physically meaningful trends for both reaction and scattering probabilities compared to experimental and other theoretical results.

Stochastic wave packet approach to nonadiabatic scattering of diatomic molecules from metals

In this contribution, we present a quantum dynamical approach to study inelastic scattering of diatomic molecules from metal surfaces at normal incidence. The dissipative dynamics obeys a stochastic Schrödinger equation describing the time-evolution of the system as a piecewise deterministic process. Energy exchange between the molecular vibrational degrees of freedom and the metal electrons is represented using operators in tensor product form, which are coupled via anharmonic transition rates calculated from first-order perturbation theory. Full dimensional observables are obtained by averaging over simulations in 4D-including the internal stretch, the distance to the surface, and the orientation angles-at different surface sites. The method is applied to the state-resolved scattering of vibrationally excited NO from Au(111), revealing important channels for quantized energy relaxation.

Quantum and classical dynamics of reactive scattering of H2 from metal surfaces

State-of-the-art theoretical models allow nowadays an accurate description of H₂/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.

The effect of phonon modes and electron–hole pair couplings on molecule–surface scattering processes

The effect of phonon modes and electron–hole pair (elhp) couplings at different surface temperature on D 2(v = 0, 1; j = 0)– Cu (111) collision has been explored by assuming weakly correlated interactions between molecular Degrees of Freedoms (DOFs) with surface modes and elhp excitations through a Hartree product type wavefunction, where the initial state distributions for the phonon modes and the elhp couplings are incorporated by using Bose–Einstein and Fermi–Dirac probability factors, respectively. We carry out four (4D⊗2D)- and six (6D)- dimensional quantum dynamics on such an effective Hamiltonian, and depict the calculated sticking/transition probabilities and energy transfer from molecule to the surface. The phonon modes slightly affect the sticking probability by broadening the profile, but the transition probability are substantially changed with respect to the rigid surface. On the contrary, the inclusion of elhp coupling along with phonon modes does not change the results much compared to the only phonon case.

The direct and precursor mediated dissociation rates of H2 on a Ni(111) surface

The dissociation and recombination rates of physisorbed H2, and the direct and steady state dissociation (i.e., the precursor mediated dissociation) rates of gas phase H2 on Ni(111), as well as the corresponding kinetic isotope effects, are calculated using the quantum instanton method, together with path integral Monte Carlo and adaptive umbrella sampling techniques. All these rates except the recombination one first decrease and then increase with the increasing temperature, and their minimum values appear at about 250, 300 and 250 K, respectively. These non-monotonic behaviors reveal that the quantum effect of H2 should be very remarkable at low temperatures. The steady state rates are smaller than the direct rates at low temperatures, however, they become larger than the direct ones at high temperatures, these two kinds of rates become equal at about 400 and 300 K on the rigid and quantum lattices, respectively. The quantum motion of the lattice can enhance the direct and steady state rates, and it increases the steady state rate much more than the direct one, for instance, the direct and steady state rates on the quantum lattice are 1.30 and 2.08 times larger than that on the rigid one at 300 K. The calculated kinetic isotope effects are much larger than 1, which reveals that H2 always has a larger rate than that of D2, and the direct process predicts much larger kinetic isotope effects than the steady state process at low temperatures. In addition, the kinetic isotope effects are not affected by the lattice motion.