Corroboration of theory for H + D2 --> D + HD (v' = 3, j' = 0) reactive scattering dynamics

The differential cross section (DCS) for the reaction H + D2 --> D + HD (v' = 3, j' = 0) exhibits particularly rich dynamics; in addition to the expected direct recoil backscattering feature, a surprising time-delayed forward scattering feature appears that has been attributed to glory scattering arising from nearside and farside interference. This fact leads to a complex DCS that depends strongly on the collision energy. Its accurate calculation requires a fully quantum mechanical (QM) treatment. We report improved measurements of this DCS over the collision energy range 1.55 < or = E(coll) < or = 1.82 eV. Previous measurements using the core extraction method, while generally in agreement with theory, lacked sufficient resolution to capture all of the noteworthy behavior of the system; in the present work, we use ion imaging to observe many previously unresolved features of the DCS, particularly in the forward-scattered region. Agreement with QM calculations is found at all collision energies, reconciling an earlier discrepancy between experiment and theory near E(coll) = 1.54 eV.     

Observation of partial wave structures in the integral cross section of the S((1)D2) + H2(j = 0) reaction

The integral cross section of the S((1)D(2)) + H(2)(j = 0) → SH + H reaction has been measured for the first time at collision energies from 0.820 down to 0.078 kJ mol(-1) in a high-resolution crossed beam experiment. The excitation function obtained exhibits a non-monotonic variation with collision energy and compares well with the results of high-level quantum calculations. In particular, the structures observed in the lower energy part, where only a few partial waves contribute, can be described in terms of the sequential opening of individual channels, consistent with the theoretical calculations.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.     

Disagreement between theory and experiment in the simplest chemical reaction: collision energy dependent rotational distributions for H + D2 --> HD(nu' = 3,j') + D

We present experimental rotational distributions for the reaction H + D2 --> HD(nu' = 3,j') + D at eight different collision energies between 1.49 and 1.85 eV. We combine a previous measurement of the state-resolved excitation function for this reaction [Ayers et al., J. Chem. Phys. 119, 4662 (2003)] with the current data to produce a map of the relative reactive cross section as a function of both collision energy and rotational quantum number (an E-j' plot). To compare with the experimental data, we also present E-j' plots resulting from both time-dependent and time-independent quantum mechanical calculations carried out on the BKMP2 surface. The two calculations agree well with each other, but they produce rotational distributions significantly colder than the experiment, with the difference being more pronounced at higher collision energies. Disagreement between theory and experiment might be regarded as surprising considering the simplicity of this system; potential causes of this discrepancy are discussed.     

Time dependent quantum dynamics study of the Ne + H2+(v=0-4)-->NeH+ + H proton transfer reaction

The Ne + H2+-->NeH+ + H proton transfer reaction was studied using the time dependent real wave packet quantum dynamics method at the helicity decoupling level, considering the H2+ molecular ion in the (v=0-4, j=0) vibrorotational states and a wide collision energy interval. The calculated reaction probabilities and reaction cross sections were in a rather good agreement with reanalyzed previous exact quantum dynamics results, where a much smaller collision energy interval was considered. Also, a quite good agreement with experimental data was found. These results suggested the adequacy of the approach used here to describe this and related systems.     

Collision-energy dependence of HD(nu'=1,j') product rotational distributions for the H + D2 reaction

Product rotational distributions for the reaction H + D2 --> HD(nu'=1,j') + D have been measured for 16 collision energies in the range of 1.43 < or = E(coll) < or = 2.55 eV. Time-dependent quantum-mechanical calculations agree well in general with the experimental results, but they consistently yield slightly colder distributions. In terms of the average energy channeled into rotation, the differences between experiment and theory amount to approximately 10% for all collision energies sampled. No peculiarity is found for E(coll)=2.55 eV at which the system has sufficient energy to access the first HD2 electronically excited state.     

Collision energy dependence of the HD(nu' = 2) product rotational distribution of the H + D2 reaction in the range 1.30-1.89 eV

An experimental and theoretical investigation of the collision energy dependence of the HD(nu' = 2,j') rotational product state distribution for the H + D2 reaction in the collision energy range of Ecol = 1.30-1.89 eV has been carried out. Theoretical results based on time-dependent and time-independent quantum mechanical methods agree nearly perfectly with each other, and the agreement with the experiment is good at low collision energies and very good at high collision energies. This behavior is in marked contrast to a previous report on the HD(nu' = 3,j') product state rotational distribution [Pomerantz et al., J. Chem. Phys. 120, 3244 (2004)] where a systematic difference between experiment and theory was observed, especially at the highest collision energies. The reason for this different behavior is not yet understood. In addition, this study employs Doppler-free spectroscopy to resolve an ambiguity in the E, F-X resonantly enhanced multiphoton ionization transition originating from the HD(nu' = 2,j' = 1) state, which is found to be caused by an accidental blending with the transition coming from the HD(nu' = 1,j' = 14) state.     

A time-dependent wave-packet quantum scattering study of the reaction H2+(v = 0-2,4,6;j = 1) + He--> HeH+ + H

The quantum scattering dynamics calculation was carried out for the titled reaction in the collision energy range of 0.0-2.4 eV with reactant H(2) (+) in the rotational state j = 1 and vibrational states v = 0-2, 4, and 6. The present time-dependent wave-packet calculation takes into account the Coriolis coupling (CC) and uses the accurate ab initio potential-energy surface of Palmieri et al. [Mol. Phys. 98, 1835 (2000)]. The importance of including the CC quantum scattering calculation has been revealed by the comparison between the CC calculation and the previous coupled state (CS) calculation. The CC total cross sections for the v = 2, 4, and 6 states show collision energy-dependent behaviors different from those based on the CS calculation. Furthermore, the collision energy dependence of the total cross sections obtained in the present CC calculation only exhibits minor oscillations, indicating that the chance is slim for reactive resonances in total cross sections to survive through the partial-wave averaging. The magnitude and profile of the CC total cross sections for v = 0-2 in the collision energy range of 0.0-2.5 eV are found to be consistent with experimental cross sections obtained recently by Tang et al. [J. Chem. Phys. 122, 164301 (2005)] after taking into account the experimental uncertainties.     

A time-dependent wave packet quantum scattering study of the reaction HD+ (v = 0 - 3;j0 = 1) + He --> HeH+(HeD+) + D(H)

Time-dependent wave packet quantum scattering (TWQS) calculations are presented for HD(+) (v = 0 - 3;j(0)=1) + He collisions in the center-of-mass collision energy (E(T)) range of 0.0-2.0 eV. The present TWQS approach accounts for Coriolis coupling and uses the ab initio potential energy surface of Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. For a fixed total angular momentum J, the energy dependence of reaction probabilities exhibits quantum resonance structure. The resonances are more pronounced for low J values and for the HeH(+) + D channel than for the HeD(+) + H channel and are particularly prominent near threshold. The quantum effects are no longer discernable in the integral cross sections, which compare closely to quasiclassical trajectory calculations conducted on the same potential energy surface. The integral cross sections also compare well to recent state-selected experimental values over the same reactant and translational energy range. Classical impulsive dynamics and steric arguments can account for the significant isotope effect in favor of the deuteron transfer channel observed for HD(+)(v<3) and low translational energies. At higher reactant energies, angular momentum constraints favor the proton-transfer channel, and isotopic differences in the integral cross sections are no longer significant. The integral cross sections as well as the J dependence of partial cross sections exhibit a significant alignment effect in favor of collisions with the HD(+) rotational angular momentum vector perpendicular to the Jacobi R coordinate. This effect is most pronounced for the proton-transfer channel at low vibrational and translational energies.     

Nonadiabatic effects in the H+D2 reaction

The state-to-state dynamics of the H+D2 reaction is studied by the reactant-product decoupling method using the double many-body expansion potential energy surface. Two approaches are compared: one uses only the lowest adiabatic sheet while the other employs both coupled diabatic sheets. Rotational distributions for the reaction H+D2 (upsilon = 0, j = 0)-->HD(upsilon' = 3, j')+D are obtained at eight different collision energies between 1.49 and 1.85 eV; no significant difference are found between the two approaches. Initial state-selected total reaction probabilities and integral cross sections are also given for energies ranging from 0.25 up to 2.0 eV with extremely small differences being observed between the two sets of results, thus showing that the nonadiabatic effects in the title reaction are negligible at least for small energies below 2.0 eV.     

H+O2(n0, j0)→HO+O及C+H2(n0, j0)→CH+H体系选态反应截面的过渡态理论计算

本文用微正则过渡态理论计算了H+O_2(n_0,j_0)→HO+O和C+H_2(n_0, j_0)→CH+H在ab initio势能面上的选态反应截面σ_(n_0,j_0); E.分析了势能面性质对反应截面的影响。计算结果表明, 在指定反应物分子的振动态n_0、转动态j_0时, 两个反应体系的反应截面随相对平动能的增加先是增加后是减小(j_0=1, n_0=0除外); 在给定相对平动能和反应物分子的转动态j_0时, 随反应物分子的振动量子数n_0的增加, 两个体系的选态反应截面均有较显著的增加, 在指定相对平动能和反应物分子的振动态n_0时, H+O_2体系的选态反应截面随j_0的变化较为复杂, 而C+H_2体系则比较简单(j_0=1除外)。对于H+O_2反应体系, 本文得到的反应截面与实验结果及准经典轨迹理论的计算结果符合得很好。     

Probing Feshbach resonances in F+H2(j=1)-->HF+H: dynamical effect of single quantum H2-rotation

Full quantum state resolved scattering of the F atom reaction with H(2)(j=0) and H(2)(j=1) was investigated at the collision energies of 0.19 and 0.56 kcalmol. Dramatic difference between the dynamics for the F+H(2)(j=0,1) reactions at both collision energies have been observed. Forward scattering HF(v(')=2) products have been observed unambiguously for the F+H(2)(j=1) reaction at low collision energies, which was attributed to the Feshbach resonances. This study provides a unique case of reaction resonances involving a rotationally excited reagent.     

State-to-state quantum dynamics of the H + HBr reaction: competition between the abstraction and exchange reactions

Quantum state-to-state dynamics for the H + HBr(υ(i) = 0, j(i) =0) reaction was studied on an accurate ab intio potential energy surface for the electronic ground state of BrH(2). Both the H + HBr → H(2) + Br abstraction reaction and the H' + HBr → H'Br + H exchange reaction were investigated up to a collision energy of 2.0 eV. It was found that the abstraction channel is dominant at lower collision energies, while the exchange channel becomes dominant at higher collision energies. The total integral cross section of the abstraction reaction at a collision energy of 1.6 eV was found to be 1.37 A?(2), which is larger than a recent quantum mechanical result (1.06 A?(2)) and still significantly smaller than the experimental value (3 ± 1 A?(2)). Meanwhile, similar to the previous theoretical study, our calculations also predicted much hotter product rotational state distributions than those from the experimental study. This suggests that further experimental investigations are highly desirable to elucidate the dynamic properties of the title reactions.     

Molecular beam scattering studies of the reaction D + H2 (v = 0) and D + H2 (v = 1) → HD + H

The reaction D + H₂ → HD + H has been investigated in two molecular beam scattering experiments. Angular and time-of-flight distributions have been measured for the initial vibrational ground state (v = 0) at a most probable collision energy of E_cm = 1.5 eV and for the first vibrational excited state (v = 1) at E_cm = 0.28 eV with the same apparatus. Results for the ground-state experiment are compared with quasiclassical trajectory calculations(QCT) on the LSTH-hypersurface transformed into the laboratory system and averaged over the apparatus distributions. The agreement isquite satisfactory. At this high collision energy the HD products are no longer scattered in a backward direction but in a wide angular region concentrated about θ = 90° in the center-of-mass system. The absolute reactive cross section has been determined and the agreement with the theoretical value from QCT calculations is within the experimental error. The high sensitivity of the experiment to different properties of the doubly differential cross section has also been demonstrated. A preliminary evaluation of the experiment with initial vibrational excitation (v = 1) shows that the HD-product molecules are preferably backward scattered and the change of internal energy is small supporting the concept of a reaction which is adiabatic with respect to the internal degrees of freedom.     

Quantum and classical studies of the O(3P) + H2(v = 0-3,j = 0) --> OH + H reaction using benchmark potential surfaces

We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.     

Time dependent quantum dynamics study of the O++H2(v=0,j=0)-->OH++H ion-molecule reaction and isotopic variants (D2,HD)

The time dependent real wave packet method using the helicity decoupling approximation was used to calculate the cross section evolution with collision energy (excitation function) of the O++H2(v=0,j=0)-->OH++H reaction and its isotopic variants with D2 and HD, using the best available ab initio analytical potential energy surface. The comparison of the calculated excitation functions with exact quantum results and experimental data showed that the present quantum dynamics approach is a very useful tool for the study of the selected and related systems, in a quite wide collision energy interval (approximately 0.0-1.1 eV), involving a much lower computational cost than the quantum exact methods and without a significant loss of accuracy in the cross sections.     

Barrier Height Effect on Cl+H2(D2) Reaction

Three-dimensional time-dependent quantum wave packet calculation was performed to study the reaction dynamics of Cl+H2(D2) on two potential energy surfaces (CW PESs). The first CW PES is with spin-orbit correction; the second is without spin-orbit correction. The integral cross-section and reaction probability as a function of collision energy are calculated in the collision energy range of 0.1 eV to 1.4 eV. For reaction of Cl with D2, the reaction section with spin-orbit correction has a shift toward the high energy because the barrier height increases. As for the reaction of Cl with H2 at low collision energy, it is more reactive on the PES with spin-orbit correction than on the low barrier height PES without spin-orbit correction, due to the tunnel effect for the reaction of the Cl with H2. When the collision energy is higher than 0.7 eV, the reactivity on the low barrier height PES is larger than that on the high barrier height PES. It is believed that the barrier height plays a very important role in the reactivity of Cl with (H2, D2). For the Cl+H2 reaction the barrier width is also very important because of the tunneling effect.     

Stereodynamics of the F + HD(v = 0, j = 1) reaction: direct vs. resonant mechanisms

The stereodynamics and mechanism of the F + HD(v = 0, j = 1) → HF (DF) + D (H) reactions have been thoroughly analysed at collision energies in the 0-160 meV range. Specifically, this study is focused on (i) the comparison between the stereodynamics of the collisions leading to HF and DF formation, and (ii) the stereodynamical fingerprints of the resonance that occurs at low collision energies in the HF channel and whose manifestation in the total cross section is greatly diminished for initial j > 0. While previous studies were limited to the analysis of integral cross sections (ICS), differential cross sections (DCS) and reaction probabilities, in the present work we have included the analysis of vectorial quantities such as the direction of the initial rotational angular momentum and internuclear axis, and their effect on reactivity. In particular, polarisation parameters (PP) and polarisation dependent differential cross sections (PDDCS), quantities that describe how the intrinsic HD rotational angular momentum and molecular axis polarisations contribute to reaction, are calculated and examined. The evolution of the PPs with the collision energy differs markedly between the two reaction channels. For the DF channel, the PP values are small and change very little in the energy range in which DF formation is appreciable. In contrast, rapid fluctuations in the magnitude and sign of the PPs are observed in the HF channel at low collision energies in and around the resonance. As the collision energy increases, direct (non-resonant) scattering prevails, and the various quantities are reasonably well accounted for by the QCT calculations, as in the case of the DF channel. The intrinsic directional information has been used to access the extent of control that can be achieved through polarisation of the HD molecule prior to collision. It was found that the same extrinsic preparation leads to very different outcomes on the HF channel DCS when the collision energy is close to the resonance. It is also shown that polarisation of the HD internuclear axis along the initial relative velocity enhances the effect of the resonance and allows its clear identification. Finally, the effect of different extrinsic preparations on the angle-velocity DCS is found to be strong, thus allowing considerable control of product angular distributions.     

An experimental and quasiclassical trajectory study of the rovibrationally state-selected reactions: HD+(v=0-15,j=1)+He-->HeH+(HeD+)+D(H)

The absolute integral cross sections for the formation of HeH+ and HeD+ from the collisions of HD+(v,j=1)+He have been examined over a broad range of vibrational energy levels v=0-13 at the center-of-mass collision energies (ET) of 0.6 and 1.4 eV using the vacuum ultraviolet (VUV) pulsed field ionization photoelectron secondary ion coincidence method. The ET dependencies of the integral cross sections for products HeH+ and HeD+ from HD+(v=0-4)+He collisions in the ET range of 0-3 eV have also been measured using the VUV photoionization guided ion beam mass spectrometric technique, in which vibrationally selected HD+(v) reactant ions were prepared via excitation of selected autoionization resonances of HD. At low total energies, a pronounced isotope effect is observed in absolute integral cross sections for the HeH++D and HeD++H channels with significant favoring of the deuteron transfer channel. As v is increased in the range of v=0-9, the integral cross sections of the HeH++D channel are found to approach those of HeD++H. The observed velocity distributions of products HeD+ and HeH+ are consistent with an impulsive or spectator-stripping mechanism. Detailed quasiclassical trajectory (QCT) calculations are also presented for HD+(v,j=1)+He collisions at the same energies of the experiment. The QCT calculations were performed on the most accurate ab initio potential energy surface available. If the zero-point energy of the reaction products is taken into account, the QCT cross sections for products HeH+ and HeD+ from HD+(v)+He are found to be significantly lower than the experimental results at ET values near the reaction thresholds. The agreement between the experimental and QCT cross sections improves with translational energy. Except for prethreshold reactivity, QCT calculations ignoring the zero-point energy in the products are generally in good agreement with experimental absolute cross sections. The experimental HeH+/HeD+ branching ratios for the HD+(v=0-9)+He collisions are generally consistent with QCT predictions. The observed isotope effects can be rationalized on the basis of differences in thermochemical thresholds and angular momentum conservation constraints.     

Vibrational and translational enhancement of endothermic reactions: K + HCl(υ = 0,1) and K + HF(υ = 0,1)

The relative integral cross section for the two endothermic reactions K + HCl(υ = 0 and 1) → KCl + H and K + HF(υ = 0 and 1) → KF + H has been measured as a function of the collision energy E using the crossed molecular beam technique. The vibrationally excited state (υ = 1) has been populated thermally by heating the beam source to temperatures around 2000 K. The variation of the collision energy from thermal up to around 2.1 eV was achieved by seeding the K-beam with various carrier gases. The molecular reaction product was detected by surface ionization in connection with a time-of-flight method. The total energy threshold of the reactions has been found to be equal to or only slightly above the corresponding endothermicities. This suggests a vanishing or very low barrier crest on the potential energy hypersurfaces which is contradictory to recent theoretical results. The inclusion of tunneling in case of K + HF leads to a negligible rise of the barrier heights. The efficacy of translational and vibrational energy in promoting the reactive process has been directly compared over a wide range of collision energies. For K + HCl the vibrational enhancement of the reactivity descends with increasing E from approximately a factor of 10 at E = 0.08 eV to around unity for E ? 0.5 eV. The good agreement of this experimental result with phase space calculations suggests that the marked enhancements are predominantly caused by the long-range attraction between reagents in connection with an “early” barrier on the potential energy surface. In case of K + HF vibrational energy is by a factor of up to 380 more favourable in promoting the reaction than the same amount of translational energy. Again, with rising collision energy its efficacy decreases but promotes the reaction still by a factor of 70 at E = 1.7 eV. Since phase space theory fails here the reaction is certainly non-statistical and we conclude that the observed large efficacy of vibrational energy is due to a “late” barrier. The proposed barrier positions for the two systems are in accordance with theoretical results.