Vibrational control of the reaction pathway in the H + CHD3 → H2 + CD3 reaction

Higher-Order Split Operator Schemes for Solving Tetratomic Reactions Using the Time-Dependent Wave Packet Method

In this work, using the time-dependent quantum wave packet method, quite a few typical higher-order split operators (HOSOs) were for the first time applied to calculate the tetratomic reactive scattering processes in the hyperspherical coordinate. It was found that the HOSOs were hardly efficient for a tetratomic reaction calculation, unlike those for a triatomic reactive scattering calculation. We proposed an efficient HOSO with a force gradient (denoted as 2G1 in the main text) for efficiently and accurately calculating a tetratomic reaction using the quantum wave packet method. Several typical tetratomic reactions, such as H2 + OH, HF + OH, and H2 + OH+, are calculated for demonstrating the effectiveness of the proposed 2G1 in terms of (product state-resolved) reaction probability and inelastic probability, by comparing with the performance of the previously reported various HOSOs. We suggest that the 2G1 propagator could be applied to efficiently calculate a general tetratomic reaction.

High vibrational excitation of the reagent transforms the late-barrier H + HOD reaction into an early-barrier reaction

Polanyi's rules predict that a late-barrier reaction yields vibrationally cold products; however, experimental studies showed that the H2 product from the late-barrier H + H2O(|04⟩-) and H + HOD(vOH = 4) reactions is vibrationally hot. Here, we report a potential-averaged five-dimensional state-to-state quantum dynamics study for the H + HOD(vOH = 0-4) → H2 + OD reactions on a highly accurate potential energy surface with the total angular momentum J = 0. It is found that with the HOD vibration excitation increasing from vOH = 1 to 4, the product H2 becomes increasingly vibrationally excited and manifests a typical characteristic of an early barrier reaction for vOH = 3 to 4. Analysis of the scattering wave functions revealed that vibrational excitation in the breaking OH bond moves the location of dynamical saddle point from product side to reactant side, transforming the reaction into an early barrier reaction. Interestingly, pronounced oscillatory structures in the total and product vibrational-state-resolved reaction probabilities were observed for the H + HOD(vOH = 3, 4) reactions, in particular at low collision energies, which originate from the Feshbach resonance states trapped in the bending/torsion excited vibrational adiabatic potential wells in the entrance region due to van der Waals interactions.

State-to-State Dynamics in Mode-Specific Reactions of Cl + CH3D(v1-I, v1-II, and v4 = 1; |10⟩): Loss of Memory or Not

Several decades of the study of reaction dynamics culminate in the concept of mode specificity and bond selectivity in polyatomic systems. Until very recently, the main concern of those studies has been total reactivity and little attention has been paid to the mode-specific effects on the more detailed product-state and angular distributions. Conventional wisdom would anticipate that the fine detail should reveal a more pronounced mode dependency. However, a few recent studies showed that the product distributions could appear to be surprisingly insensitive to the modes of internal excitation of reagents. This counterintuitive finding led to a concept of loss of memory. Here, we present detailed experimental results in the reactions of the Cl atom with three distinct stretching-excited CH₃D(v_CH3 = 1) reagents. In conjunction with the previous reports on various aspects of this reaction, such a comprehensive set of data enables us to perform an in-depth examination of the validity of this new concept.

Stretching-mode specificity in the Cl + CH3D(v1-I, v1-II, and v4 = 1; |jK〉) reactions: dependency on the initial |jK〉 selectivity

The title reactions were studied at a collisional energy of 5.4 kcal mol^-1 in a crossed-beam product-imaging experiment. The dynamics attributes of the dominant ground-state CH₂D(0₀) and the accompanied C-D bend-excited CH₂D(6₁) products were imaged in reactions with totally 16 ro-vibrationally selected states of the CH₃D(v_i, |jK〉) reagents. We found that all three vibrational excitations yielded marked |jK〉-dependent rate-enhancements in forming the (0₀, 0/1)_s product pairs. Furthermore, for a given rotational |jK〉-mode, a vibrational-mode propensity of v₄ > v₁-I > v₁-II in rate promotion and a clear manifestation of the Fermi-phase-induced interference effect of the latter two were observed. Compared to the reactivity of the rotationless state |jK〉 = |00〉, a minute rotational-excitation of all three stretch-excited CH₃D(v_i = 1) reagents could yield significantly higher reaction rates for the product pair (0₀, 0)_s, but not so for (0₀, 1)_s. The signals in forming the (6₁, 0)_s pair were clearly notable but smaller than that of the ground-state reaction product pair, (0₀, 0)_g. An opposite propensity of v₁-II ≈ v₁-I > v₄ with a milder dependency on the initial |jK〉-states was observed. The angular distributions of the (0₀, 0)_s pairs were nearly identical for all ro-vibrationally excited reagents, displaying the typical trait for a direct abstraction of the rebound mechanism. Similar distributions were found for the (6₁, 0)_s pairs; yet, both pairs deviated substantially from the peripheral feature of the ground-state reaction pair of (0₀, 0)_g. Those of the (0₀, 1)_s pairs in reactions with v₄-excitation featured a prominent forward-peaking distribution-suggestive of a time-delayed, resonance-mediated pathway, again with little dependency on the initial |jK〉-states. As for the reactions with the two Fermi-dyads, v₁-I and v₁-II, albeit showing globally similar distributions to that for v₄, a substantial variation with the initial rotational-mode excitation could be discerned in the forward-peaking features. To unravel the impact of the Fermi-phase on the |jK〉-dependent attributes, we adopted a comparative approach by contrasting the observations in reactions with the Fermi-dyad reagents (the superposition states) to those with the pure-state reagents. Remarkable distinctions are unveiled and elucidated with the unexplained results explicitly pointed out, which call for future theoretical investigations for deeper understanding.

Multilayer Multiconfiguration Time-Dependent Hartree Study on the Mode-/Bond-Specific Quantum Dynamics of Water Dissociation on Cu(111)

In this work, full-dimensional (9D) quantum dynamics calculations on mode-/bond-specific surface scattering of a water molecule on a copper (111) rigid surface are performed through the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. To easily perform the ML-MCTDH calculations on such a triatomic molecule-surface system, we first choose specific Jacobi coordinates as a set of coordinates of water. Next, to efficiently perform the 9D ML-MCTDH wavepacket propagation, the potential energy surface is transferred to a canonical polyadic decomposition form with the aid of a Monte Carlo-based method. Excitation-specific dissociation probabilities of H₂O on Cu(111) are computed, and mode-/bond-specific dynamics are demonstrated by comparison with a probability curve computed for a water molecule in the ground state. The dependence of the dissociation probability of the initial state of H₂O is studied, and it is found that the excitation-specific dissociation probabilities can be divided into three groups. We find that the vibrationally excited states enhance the dissociation reactivity of H₂O, while the rotationally excited states hardly influence it.