Influence of collision energy and reagent vibrational excitation on the dynamics of the reaction H + LiH

We present a detailed quasiclassical trajectory (QCT) study of the dynamics corresponding to the reaction H + LiH proceeding via depletion and H‐exchange paths on a new potential energy surface of the electronic ground state. The effects of collision energy and reagent initial vibrational excitation on the reaction probability and cross sections are studied over a wide range of collision energies. The QCT‐calculated reaction probability and cross sections are in good agreement with previous time‐dependent wave packet results. More importantly, we found that the vibrational excitation of LiH molecule inhibits the LiH depletion reaction, whereas it promotes the H‐exchange reaction. In addition, the differential cross sections calculated for the depletion reaction at different collision energies and excitation states indicate a strong forward scattering of the product molecule H₂. © 2013 Wiley Periodicals, Inc.     

Quantum study of Eley-Rideal reaction and collision induced desorption of hydrogen atoms on a graphite surface. I. H-chemisorbed case

Collision induced (CI) processes involving hydrogen atoms on a graphite surface are studied quantum mechanically within the rigid, flat surface approximation, using a time-dependent wave packet method. The Eley-Rideal (ER) reaction and collision induced desorption (CID) cross sections are obtained with the help of two propagations which use different sets of coordinates, a "product" and a "reagent" set. Several adsorbate-substrate initial states of the target H atom in the chemisorption well are considered, and CI processes are studied over a wide range of projectile energy. Results show that (i) the Eley-Rideal reaction is the major reactive outcome and (ii) CID cross sections do not exceed 4 A2 and present dynamic thresholds for low values of the target vibrational quantum number. ER cross sections show oscillations at high energies which cannot be reproduced by classical and quasiclassical trajectory calculations. They are related to the vibrational excitation of the reaction products, which is a rather steep decreasing function of the collision energy. This behavior causes a selective population of the low-lying vibrational states and allows the quantization of the product molecular states to manifest itself in a collisional observable. A peak structure in the CID cross section is also observed and is assigned to the selective population of metastable states of the transient molecular hydrogen.     

Fully converged integral cross sections of collision induced dissociation, four-center, and single exchange reactions, and accuracy of the centrifugal sudden approximation in H2 + D2 reaction

The initial state selected time-dependent wave packet method was employed to calculate the integral cross sections for the H(2) + D(2) reaction with and without the centrifugal sudden (CS) approximation by including all important K (the projection of the total angular momentum on the body-fixed axis) blocks. With a full-dimensional model, the first fully converged coupled-channel (CC) cross sections for different competitive processes from the ground rotational state were obtained: collision induced dissociation (CID), four-center (4C) reaction and single exchange (SE) reaction. The effect of the total angular momentum J on the reaction dynamics of H(2) + D(2) and the accuracy of the CS approximation have also been studied. It was found that the CID and SE processes occur in a wide range of J values while the 4C process can only take place in a narrow window of J values. For this reason, the CC cross section for the 4C channel is merely comparable to the SE channel. A comparison of the integral cross sections from CC and CS calculations showed that the CS approximation works well for the CID process but not for the 4C and SE processes, and the discrepancy between the CC and CS cross sections grows larger as the translational energy and/or the vibrational energy increase(s).     

Quantum dynamics of H + LiH reaction and its isotopic variants

Time-dependent quantum wave packet dynamics study is carried out to investigate the initial state selected channel specific reactivity of H + LiH collisional system on a new and more accurate ab initio potential energy surface developed by Wernli et al. [J. Phys. Chem. A 113, 1121 (2009)]. The H + LiH reaction proceeds through LiH depletion and H-exchange paths. While the former path is highly exoergic (by ～2.258 eV), the latter path is thermoneutral. State selected and energy resolved integral reaction cross sections and thermal rate constants are reported and compared with the literature data. The reactivity of the LiH depletion channel is found to be greater than the H-exchange channel. Rotational excitation of the reagent LiH molecule causes a decrease of reactivity of both the channels. On the other hand, the vibrational excitation of the reagent LiH decreases the reactivity of the LiH depletion channel and increases the reactivity of the H-exchange channel. The effect of isotopic substitution (H by D) on the reaction dynamics is also examined.     

Cross sections and rate constants for OH + H2 reaction on three different potential energy surfaces for ro-vibrationally excited reagents

A systematic study of the reagent ro-vibrational excitations in H(2) + OH reaction is presented on three different potential energy surfaces using the multiconfiguration time-dependent Hartree method. An exact form of the kinetic energy operator including Coriolis coupling has been used. Coupled channel results on WDSE surface for vibrational excitation of H(2) produce very large cross sections in accordance with the previous approximate results. The rate constant obtained for H(2)(v = 1) at 300 K on the YZCL2 surface shows an excellent agreement with the most recent experimental result. Quantum dynamical results for ro-vibrational excitation of reagents obtained on the WSLFH surface show similar behavior to previous quasiclassical trajectory studies. The integral cross sections obtained for excited reagent rotations exhibit contrasting trends on the three surfaces. The effects are explained considering the different orientations of the transition state structure and the individual surface characteristics.     

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.     

Bimolecular reactions of vibrationally excited molecules. Roaming atom mechanism at low kinetic energies

Quasiclassical trajectory calculations have been performed for the H + H'X(v) → X + HH' abstraction and H + H'X(v) → XH + H' (X = Cl, F) exchange reactions of the vibrationally excited diatomic reactant at a wide collision energy range extending to ultracold temperatures. Vibrational excitation of the reactant increases the abstraction cross sections significantly. If the vibrational excitation is larger than the height of the potential barrier for reaction, the reactive cross sections diverge at very low collision energies, similarly to capture reactions. The divergence is quenched by rotational excitation but returns if the reactant rotates fast. The thermal rate coefficients for vibrationally excited reactants are very large, approach or exceed the gas kinetic limit because of the capture-type divergence at low collision energies. The Arrhenius activation energies assume small negative values at and below room temperature, if the vibrational quantum number is larger than 1 for HCl and larger than 3 for HF. The exchange reaction also exhibits capture-type divergence, but the rate coefficients are larger. Comparisons are presented between classical and quantum mechanical results at low collision energies. At low collision energies the importance of the exchange reaction is enhanced by a roaming atom mechanism, namely, collisions leading to H atom exchange but bypassing the exchange barrier. Such collisions probably have a large role under ultracold conditions. The calculations indicate that for roaming to occur, long-range attractive interaction and small relative kinetic energy in the chemical reaction at the first encounter are necessary, which ensures that the partners can not leave the attractive well. Large orbital angular momentum of the primary products (equivalent to large rotational excitation in a unimolecular reaction) is favorable for roaming.     

Dynamical properties of S(3P) + HD reaction on 13A″state and their quantum wavepacket calculation

The time‐dependent wavepacket method is used to study the reaction dynamics of S(³P) + HD (v = 0, 1, 2) on the adiabatic 1³A″ potential energy surface constructed by Han and coworkers [J. Chem. Phys. 2012, 136, 094308]. The reaction probabilities and integral cross sections as a function of collision energy are obtained and discussed. The results calculated by using the CC and the CS approximation have been compared, which suggests that for this direct abstraction reaction, the cheaper CS approximation calculation is valid enough in the quantum calculation. The investigation also shows that the reaction probabilities and integral cross sections tend to increase with collision energy. By analyzing the v‐dependent behavior of the integral cross sections, the significant effect of the vibrational excitation of HD is found. Also found in the calculation is a significant resonance feature in the reaction probabilities versus collision energy. © 2014 Wiley Periodicals, Inc.     

H+ versus D+) transfer from HOD+ to N2: mode- and bond-selective effects

Reactions of HOD(+) with N(2) have been studied for HOD(+) in its ground state and with one quantum of excitation in each of its vibrational modes: (001)--predominately OH stretch, 0.396 eV, (010)--bend, 0.153 eV, and (100)--predominately OD stretch, 0.293 eV. Integral cross sections and product recoil velocities were recorded for collision energies from threshold to 4 eV. The cross sections for both H(+) and D(+) transfer rise slowly from threshold with increasing collision energy; however, all three vibrational modes enhance reaction much more strongly than equivalent amounts of collision energy and the enhancements remain large even at high collision energy, where the vibration contributes less than 10% of the total energy. Excitation of the OH stretch enhances H(+) transfer by a factor of ～5, but the effect on D(+) transfer is only slightly larger than that from an equivalent increase in collision energy, and smaller than the effect from the much lower energy bend excitation. Similarly, OD stretch excitation strongly enhances D(+) transfer, but has essentially no effect beyond that of the additional energy on H(+) transfer. The effects of the two stretch vibrations are consistent with the expectation that stretching the bond that is broken in the reaction puts momentum in the correct coordinate to drive the system into the exit channel. From this perspective it is quite surprising that bend excitation also results in large (factor of 2) enhancements of both H(+) and D(+) transfer channels, such that its effect on the total cross section at collision energies below ～2 eV is comparable to those from the two stretch modes, even though the bend excitation energy is much smaller. For collision energies above ～2 eV, the vibrational effects become approximately proportional to the vibrational energy, though still much larger than the effects of equivalent addition of collision energy. Measurements of the product recoil velocity distributions show that reaction is direct at all collision energies, with roughly half the products in a sharp peak corresponding to stripping dynamics and half with a broad and approximately isotropic recoil velocity distribution. Despite the large effects of vibrational excitation on reactivity, the effects on recoil dynamics are small, indicating that vibrational excitation does not cause qualitative changes in the reaction mechanism or in the distribution of reactive impact parameters.     

Coriolis-coupled wave packet dynamics of H + HLi reaction

We investigated the effect of Coriolis coupling (CC) on the initial state-selected dynamics of H+HLi reaction by a time-dependent wave packet (WP) approach. Exact quantum scattering calculations were obtained by a WP propagation method based on the Chebyshev polynomial scheme and ab initio potential energy surface of the reacting system. Partial wave contributions up to the total angular momentum J=30 were found to be necessary for the scattering of HLi in its vibrational and rotational ground state up to a collision energy approximately 0.75 eV. For each J value, the projection quantum number K was varied from 0 to min (J, K(max)), with K(max)=8 until J=20 and K(max)=4 for further higher J values. This is because further higher values of K do not have much effect on the dynamics and also because one wishes to maintain the large computational overhead for each calculation within the affordable limit. The initial state-selected integral reaction cross sections and thermal rate constants were calculated by summing up the contributions from all partial waves. These were compared with our previous results on the title system, obtained within the centrifugal sudden and J-shifting approximations, to demonstrate the impact of CC on the dynamics of this system.     

Reaction of C2H2(+) (n · ν2, m · ν5) with NO2: reaction on the singlet and triplet surfaces

Integral cross sections and product recoil velocity distributions were measured for reaction of C(2)H(2)(+) with NO(2), in which the C(2)H(2)(+) reactant was prepared in its ground state, and with mode-selective excitation in the cis-bend (2ν(5)) and CC stretch (n · ν(2), n = 1, 2). Because both reactants have one unpaired electron, collisions can occur with either singlet or triplet coupling of these unpaired electrons, and the contributions are separated based on distinct recoil dynamics. For singlet coupling, reaction efficiency is near unity, with significant branching to charge transfer (NO(2)(+)), O(-) transfer (NO(+)), and O transfer (C(2)H(2)O(+)) products. For triplet coupling, reaction efficiency varies between 13% and 19%, depending on collision energy. The only significant triplet channel is NO(+) + triplet ketene, generated predominantly by O(-) transfer, with a possible contribution from dissociative charge transfer at high collision energies. NO(2)(+) formation (charge transfer) can only occur on the singlet surface, and appears to be mediated by a weakly bound complex at low energies. O transfer (C(2)H(2)O(+)) also appears to be dominated by reaction on the singlet surface, but is quite inefficient, suggesting a bottleneck limiting coupling to this product from the singlet reaction coordinate. The dominant channel is O(-) transfer, producing NO(+), with roughly equal contributions from reaction on singlet and triplet surfaces. The effects of C(2)H(2)(+) vibration are modest, but mode specific. For all three product channels (i.e., charge, O(-), and O transfer), excitation of the CC stretch fundamental (ν(2)) has little effect, 2 · ν(2) excitation results in ～50% reduction in reactivity, and excitation of the cis-bend overtone (2 · ν(5)) results in ～50% enhancement. The fact that all channels have similar mode dependence suggests that the rate-limiting step, where vibrational excitation has its effect, is early on the reaction coordinate, and branching to the individual product channels occurs later.     

Distorted wave calculations of vibrational excitation in CO2+He and CO2+Ar collisions

Cross sections for the vibrational excitation and relaxation of CO₂ through collisions with He and Ar atoms have been computed using several approximate methods. They all employ the infinite-order sudden approximation for the rotational motion. The vibrational motion is treated using both exact close-coupling and distorted wave techniques. The latter approximate method permits the extension of the calculation to much higher collision energies. The validity of the distorted wave approximation is examined and is shown to be particularly good for the dominant inelastic processes in He+CO₂, leading to errors of the order of 7%. These become progressively greater for smaller cross sections. The excitation cross sections are reported for several vibrational transitions over an extended energy range up to 1.36 eV.     

A pulsed-field ionization photoelectron secondary ion coincidence study of the H2+ (X,upsilon+=0-15,N+=1)+He proton transfer reaction

The endothermic proton transfer reaction, H2+(upsilon+)+He-->HeH+ + H(DeltaE=0.806 eV), is investigated over a broad range of reactant vibrational levels using high-resolution vacuum ultraviolet to prepare reactant ions either through excitation of autoionization resonances, or using the pulsed-field ionization-photoelectron-secondary ion coincidence (PFI-PESICO) approach. In the former case, the translational energy dependence of the integral reaction cross sections are measured for upsilon+=0-3 with high signal-to-noise using the guided-ion beam technique. PFI-PESICO cross sections are reported for upsilon+=1-15 and upsilon+=0-12 at center-of-mass collision energies of 0.6 and 3.1 eV, respectively. All ion reactant states selected by the PFI-PESICO scheme are in the N+=1 rotational level. The experimental cross sections are complemented with quasiclassical trajectory (QCT) calculations performed on the ab initio potential energy surface provided by Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. The QCT cross sections are significantly lower than the experimental results near threshold, consistent with important contributions due to resonances observed in quantum scattering studies. At total energies above 2 eV, the QCT calculations are in excellent agreement with the present results. PFI-PESICO time-of-flight (TOF) measurements are also reported for upsilon+=3 and 4 at a collision energy of 0.6 eV. The velocity inverted TOF spectra are consistent with the prevalence of a spectator-stripping mechanism.     

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.     

Quantum dynamics study of Li + HF reaction

We report rigorous quantum dynamics studies of the Li + HF reaction using the time-dependent wavepacket approach. The dynamics study is carried out on a recent ab initio potential energy surface, and state-selected reaction probabilities and cross sections are calculated up to 0.4 eV of collision energy. Many long-lived resonances (as long as 10 ps) at low collision energies (below 0.1 eV) are uncovered from the dynamics calculation. These long-lived resonances play a dominant role in the title reaction at low collision energies (below 0.1 eV). At higher energies, the direct reaction process becomes very important. The reaction probabilities from even rotational states exhibit a different energy dependence than those from odd rotational states. Our calculated integral cross section exhibits a broad maximum near the collision energy of 0.26 eV with small oscillations superimposed on the broad envelope which is reminiscent of the underlying resonance structures in reaction probabilities. The energy dependence of the present CS cross section is qualitatively different from the simple J-shifting approximation, in which a monotonic increase of cross section with collision energy was obtained. Received: 8 January 1997 / Accepted: 14 January 1997     

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.     

The exchange reaction H2 (ν = 1) + H = H + H2

The classical trajectory method is applied to calculate the total cross section for the exchange reaction H₂(ν = 1) + H = H + H₂. The vibrational excitation is shown to influence efficiently the threshold value. Partial reaction rate-constants calculated on the basis of the cross sections obtained are in good agreement with those measured in H-maser experiments.     

Seven-dimensional quantum dynamics study of the O(3P)+CH4 reaction

The initial state selected time-dependent wave packet calculations have been carried out to study the title reaction with seven degrees of freedom included by restricting the nonreacting CH(3) group under C(3V) symmetry and the CH bond length in the group. Total reaction probabilities as well as integral cross sections were calculated for the ground and four vibrationally excited reagent states. Our calculation shows that the reactivity is very small for the reaction for collision energy up to 1.0 eV for all the initial states. Initial vibration excitation of CH(4), in particular, the CH stretch excitation, enhances the reactivity, but only part of the excitation energy deposited can be used to reduce the reaction threshold. The rate constant for the ground initial state agrees rather well with that from a recent quasiclassical trajectory study and is larger than that from the semirigid vibrating rotor target calculations, in particular, in the low temperature region. On the other hand, the thermal rate constant calculated from the integral cross sections for these five vibrational states is about a factor of 20 smaller than that from the multiconfiguration time-dependent Hartree study.     

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.     

Full-dimensional quantum dynamics study of the H + H2O and H + HOD exchange reactions

The initial state-selected time-dependent wave packet approach is employed to study the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H, HOH' + D exchange reactions with both OH bonds in the H(2)O reactant and OH(D) bond in the HOD reactant treated as reactive bonds. The total reaction probabilities for different partial waves, as well as the integral cross sections, which are the exact CC (coupled-channel) results, are first obtained in this study for the H(2)O(HOD) reactant initially in the ground rovibrational state. Because of the shallow C(3v) minimum along the reaction path, the reaction probabilities for the three reactions present several resonance peaks, with one dominant resonance peak just above the threshold. The cross sections for the H' + HOD → HOH' + D reaction are substantially smaller than those for the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H reactions, indicating that the H'/H exchange reactions are much more favored. In the CC calculations, the resonance peaks in the reaction probabilities diminish quickly with the increase in total angular momenta J, resulting in the existence of a clear step-like feature just above the threshold in the cross sections for the title reactions, which manifests the signature of shape resonances in these reactions. In the CS calculations, the resonance peaks on reaction probabilities persist in many partial waves, and thus the resonance structures can no longer survive the partial-wave summation and are washed out completely in the CS cross sections for the title reactions.