No abstract

We explore a coupled three-dimensional (3D) time-dependent wave packet formalism in hyperspherical coordinates for a 4D reactive scattering problem on the lowest adiabatic singlet surface (1(1)A') of the D(+) + H2 reaction. The coupling among the wavepackets arises through quantization of the rotation matrix, which represents the orientation of the three particles in space. The required transformation from Jacobi to hyperspherical coordinates and vice versa during initialization and projection of the wave packet on the asymptotic state(s), and the coupled equations of motion, are briefly discussed. With the long-range potential known to contribute significantly on the D(+) + H2 system, we demonstrate the workability of our approach, where the convergence profiles of the reaction probability for the reactive noncharge transfer (RNCT) process [D(+) + H2(v=0, j=0,1) → HD(v',j') + H(+)] are shown for three different collisional energies (1.7, 2.1, and 2.5 eV) with respect to the helicity (K) and total angular momentum (J) quantum numbers. The calculated reactive cross-section is presented as a function of the collision energy for two different initial states of the diatom (v = 0, j = 0, 1).

We include the effect of the phonon modes originating from the three layers of Cu(1nn) surface atoms on the dynamics of incoming molecular [H(2)(v, j)/D(2)(v, j)] degrees of freedom (DOFs) through a mean-field approach, where the surface temperature is incorporated into the effective potential by considering Bose-Einstein probability (BEP) factor for the initial state distribution of the surface modes calculated within harmonic approximation. Such time and temperature dependent effective Hamiltonian is further subdivided assuming a weak coupling between the two sets of molecular DOFs, namely, (x, y, z, Z) and (X, Y), respectively, in particular, to reduce the computational cost and the corresponding coupled quantum dynamical equations of motion have been formulated in terms of Time Dependent Discrete Variable Representation (TDDVR) approach. We demonstrate the workability of TDDVR method to investigate the scattering of H(2)(v, j) on Cu(1nn) surface by calculating the reaction probabilities and scattering cross-sections. Calculated results show that the phonon modes affect (a) the state-to-state transition probabilities of the scattered H(2) molecule substantially but chemisorption and physisorption processes negligibly and (b) the reaction probability of the incoming D(2) molecule noticeably.

We include the phonon modes originating from the three layers of Cu(100)/Cu(111) surface atoms on the dynamics of molecular [H(2)(v,j)/D(2)(v,j)] degrees of freedom (DOFs) through a mean field approach, where the surface temperature is incorporated into the effective Hamiltonian (potential) either by considering Boltzmann probability (BP) or by including the Bose-Einstein probability (BEP) factor for the initial state distribution of the surface modes. The formulation of effective potential has been carried out by invoking the expression of transition probabilities for phonon modes known from the "stochastic" treatment of linearly forced harmonic oscillator (LFHO). We perform four-dimensional (4D⊗2D) as well as six-dimensional (6D) quantum dynamics on a parametrically time and temperature-dependent effective Hamiltonian to calculate elastic/inelastic scattering cross-section of the scattered molecule for the H(2)(v,j)-Cu(100) system, and dissociative chemisorption-physisorption for both H(2)(v,j)-Cu(100) and D(2)(v,j)-Cu(111) systems. Calculated sticking probabilities by either 4D⊗2D or 6D quantum dynamics on an effective potential constructed by using BP factor for the initial state distribution of the phonon modes could not show any surface temperature dependence. In the BEP case, (a) both 4D⊗2D and 6D quantum dynamics demonstrate that the phonon modes of the Cu(100) surface affect the state-to-state transition probabilities of the scattered H(2) molecule substantially, and (b) the sticking probabilities due to the collision of H(2) on Cu(100) and D(2) on Cu(111) surfaces show noticeable and substantial change, respectively, as function of surface temperature only when the quantum dynamics of all six molecular DOFs are treated in a fully correlated manner (6D).

A recently proposed coupled three-dimensional time-dependent wave-packet formalism in hyperspherical coordinates is shown to yield accurate results for the reactive non-charge transfer process in the title system at collision energies as low as 100 K, where the lowest sheet of the accurate double many body expansion form for the singlet H3 (+) is used. The results are compared with available experimental data as well as time-independent calculations, and the agreement shown to be generally good.