Differential and integral cross sections for the H + O2 --> OH + O combustion reaction

We present accurate differential and integral cross sections for the H + O2 --> OH + O reaction obtained on a newly developed ab initio potential energy surface using time-independent and time-dependent quantum mechanical methods. The product angular distributions near the reaction threshold show pronounced forward and backward peaks, reflecting the complex-forming mechanism. However, the asymmetry of these peaks suggests certain nonstatistical behaviors, presumably due to some relatively short-lived resonances. The integral cross section increases monotonically with the collision energy above a reaction threshold.     

Fully Coriolis-coupled quantum studies of the H + O2 (upsilon i = 0-2, j i = 0,1) --> OH + O reaction on an accurate potential energy surface: integral cross sections and rate constants

We present accurate quantum calculations of the integral cross section and rate constant for the H + O2 --> OH + O combustion reaction on a recently developed ab initio potential energy surface using parallelized time-dependent and Chebyshev wavepacket methods. Partial wave contributions up to J = 70 were computed with full Coriolis coupling, which enabled us to obtain the initial state-specified integral cross sections up to 2.0 eV of the collision energy and thermal rate constants up to 3000 K. The integral cross sections show a large reaction threshold due to the quantum endothermicity of the reaction, and they monotonically increase with the collision energy. As a result, the temperature dependence of the rate constant is of the Arrhenius type. In addition, it was found that reactivity is enhanced by reactant vibrational excitation. The calculated thermal rate constant shows a significant improvement over that obtained on the DMBE IV potential, but it still underestimates the experimental consensus.     

Searching for resonances in the reaction Cl+CH4-->HCl+CH3: quantum versus quasiclassical dynamics and comparison with experiments

A quantum-mechanical (QM) and quasiclassical trajectory (QCT) study was performed on the title reaction, using a pseudotriatomic ab initio based surface. Probabilities and integral cross sections present some clear peaks versus the collision energy E(col), which we assign to Feshbach resonances of the transition state, where the light H atom oscillates between the heavy Cl and CH(3) groups. For ground-state reactants, reactivity is essentially of quantum origin (QCT observables and oscillations are smaller, or much smaller, than QM ones), and the calculated integral cross section and product distributions are in reasonable agreement with the experiment. The reaction occurs through an abstraction mechanism, following both a direct and an indirect mechanism. The quasiclassical trajectory calculations show the participation of a short-lived collision complex in the microscopic reaction mechanism. Finally, QCT differential cross sections of Cl+CH(4)-->HCl (nu(')=0 and 1)+CH(3) oscillate versus E(col), whereas experimentally this only occurs for HCl (nu(')=1). This theoretical result and other oscillating properties found here could, however, be related to the existence of a Feshbach resonance for the production of HCl (nu(')=1), as suggested by experimentalists.     

An ab initio-based global potential energy surface for the SH3 system and full-dimensional state-to-state quantum dynamics study for the H2 + HS → H2S + H reaction

An accurate potential energy surface for the ground electronic state of SH₃ system has been constructed with 41,882 high level ab initio energy points and the neural network fitting method. The time-dependent wave packet method has been used to calculate the first state-to-state differential cross sections for the title reaction up to 1.2 eV in full dimensions, based on the reactant–product decoupling scheme. It is found that the majority of H₂S are produced in the ground vibrational state, with a large fraction of available energy for the reaction ending up as product translational motion. The differential cross sections at the threshold energy are dominated by a very narrow peak in the backward direction. With the increase of collision energy, the width of the angular distribution increases considerably, which is a typical feature of a direct reaction via abstract mechanism, similar to the H₂ + OH → H₂O + H reaction. © 2018 Wiley Periodicals, Inc.     

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.     

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.     

Influence of collision energy and reagent rotation on the cross sections and product polarizations of the reaction F + HCl

Quasiclassical trajectory calculations have been carried out for the F+HCl reaction in three dimensions on a recent DHSN PES of the ground 1(2)A' electronic state [M. P. Deskevich, M. Y. Hayes, K. Takahashi, R. T. Skodje, and D. J. Nesbitt, J. Chem. Phys. 124, 224303 (2006)]. The effects of the collision energy and the reagent initial rotational excitation on the cross sections and product polarization are studied for the v = 0 and j ≤ 10 states of HCl over a wide collision energy range. It has been found that either the collision energy or the HCl rotational excitation increase remarkably reaction cross sections. The QCT-calculated integral cross sections are in good agreement with previous QM results. A detailed study on product polarization for the title reaction is also performed. The calculated results show that the product rotational angular momentum j' is not only aligned, but also oriented along the direction perpendicular to the scattering plane. The orientation of the HF product rotational angular momentum vector j' depends very sensitively on the collision energy and also affected by the reagent rotation. The theoretical findings and especially the roles of the collision energy and initial rotational momentum on the product polarization are discussed and reasonably explained by the HLH mass combination, the property of the PES, as well as the reactive mechanism.     

Quasi-classical Trajectory Study of the Intramolecular Isotope Effect in the Reaction O（3p）＋H2/HD

Theoretical studies of the dynamics of the reactions O（3p）＋H2/HD（ν=0, j=0）→OH＋H have been performed with quasi-classical trajectory method （QCT） on an ab initio potential surface for the lowest triplet electronic state of H2O（aA＂）. The QCT-calculated integral cross sections are in good agreement with the earlier time-dependent quantum mechanics results. The state-resolved rotational distributions reveal that the product OH rotational distributions for O＋HD have a preference for populating highly internally excited states compared with the O＋H2 reaction. Distributions of differential cross sections show that directions of scattering are strongly dependent on the choice of quantum state. The polarization dependent generalized differential cross-sections and the distributions were calculated and a pronounced isotopic effect is revealed. The calculated results indicate that the product polarization is very sensitive to the mass factor.     

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.     

S(1D)+D2反应的立体动力学理论研究

在ab initio 势能面上利用非含时量子反应散射理论研究了S(1D)+D2反应在散射能为22.18 kJ·mol-1时的立体动力学. 计算的微分反应截面呈现弱反对称的前向后向散射, 这一结果与以前的准经典结果和实验结果相一致. 态分辨的极化依赖的微分反应截面和极化参数显示产物极化并非各向异性分布, 这意味着在反应过程中由于势能面上存在一个深势阱从而形成了长寿命的复合物.     

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.     

Quantum Dynamics of Li+HF/DF Reaction Investigated by a State-to-State Time-dependent Wave Packet Approach

Using the reactant coordinate based time-dependent wave packet method, on the APW potential energy surface, the differential and integral cross sections of the Li+DF/HF(v=0, j=0, 1) reactions were calculated over the collision energy range from the threshold to 0.25 eV. The initial state-specified reaction rate constants of the title reaction were also calculated. The results indicate that, compared with the Li+DF reaction, the product LiF of Li+HF reaction is a little more rotationally excited but essentially similar. The initial rotational excitation from j=0 to 1 has little effect on the Li+DF reaction. However, the rotational excitation of DF does result in a little more rotationally excited product LiF. The different cross section of both reactions is forward biased in the studied collision energy range, especially at relatively high collision energy. The resonances in the Li+HF reaction may be identifiable as the oscillations in the product ro-vibrational state-resolved integral cross sections and backward scattering as a function of collusion energy. For the Li+HF reaction, the rate constant is not sensitive to the temperature and almost has no change in the temperature range considered. For the Li+DF reaction, the rate constant increase by a factor of about 10 in the temperature range of 100?300 K. Brief comparison for the total reaction probabilities and integral cross section of the Li+HF reaction has been carried out between ours and the values reported previously. The agreement is good, and the difference should come from the better convergence of our present calculations.     

Three-vector correlation in statistical reactions: the role of the triatomic parity

This article presents a methodology for the determination of the k-j-k' three-vector correlation assuming a statistical model for atom-diatom reactions; k and k' are the reagent-approach and product-recoil directions, respectively, and j is the rotational angular momentum of the reagent diatomic. Although the polarization of reagent angular momentum is in most cases negligible, conservation of the triatomic parity imposes a certain polarization for some combinations involving low reagent and product rotational states. Statistical and quantum-mechanical polarization-dependent differential cross sections were calculated for the barrierless D(+) + H(2)(v = 0,j) → HD(v' = 0,j') + H(+) reaction. The agreement between the two is in most cases excellent, confirming the statistical character of the reaction at low and moderate collision energies.     

Time-dependent reactive scattering for the system H- + D2 <--> HD + D- and comparison with H- + H2 <--> H2 + H-

This work presents results of quantum mechanical calculations of reaction probabilities for the ion-neutral molecule collisions H- + D2 <--> HD + D-. Time-dependent wave packet propagations for total angular momentum J not equal to 0, including the full Coriolis coupling, are performed. The calculated state-to-state reaction probabilities using product Jacobi coordinates are compared with energy-resolved reaction probabilities calculated with the flux-operator using reactant Jacobi coordinates and with time-independent calculations. Differences between nearly converged integral cross sections and those using the J-shifting method and centrifugal sudden approximation and comparison with experimental results will be presented.     

State-to-state dynamics of H + O2 reaction, evidence for nonstatistical behavior

Converged differential and integral cross sections are reported for the H + O2 --> OH + O reaction on an improved potential energy surface of HO2(X2A') using a dynamically exact quantum wave packet method and Gaussian weighted quasi-classical trajectory method. The complex-forming mechanism is confirmed by strong forward and backward scattering peaks and by highly inverted OH rotational state distributions. Both the quantum and classical results provide strong evidence for nonstatistical behavior in this important reaction.     

Differential and integral cross sections of the N(2D)+H2-->NH+H reaction from exact quantum and quasi-classical trajectory calculations

Exact quantum mechanical state-to-state differential and integral cross sections and their energy dependence have been calculated on an accurate NH2 potential energy surface (PES), using a newly proposed Chebyshev wave packet method. The NH product is found to have a monotonically decaying vibrational distribution and an inverted rotational distribution. The product angular distributions peak in both forward and backward directions, but with a backward bias. This backward bias is more pronounced than that observed previously on a less accurate PES. Both the differential and integral cross sections oscillate mildly with collision energy, indicating the dominance of short-lived resonances. The quantum mechanical results are compared with those obtained from quasi-classical trajectories. The agreement is generally reasonable and the discrepancies can be attributed to the neglect of quantum effects such as tunneling. Detailed analysis of the trajectories revealed that the backward bias in the differential cross section stems overwhelmingly from the fast insertion component of the reaction, augmented with some flux from the abstraction channel, particularly at high collision energies.     

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.     

Differential reaction cross sections from rotationally resolved quantum scattering calculations: application to gas-phase S(N)2 reactions

Differential reaction cross sections have been computed based on previous rotationally resolved time-independent quantum-mechanical scattering calculations for the complex-forming S(N)2 reaction Cl(-) + CH(3)Br → ClCH(3) + Br(-). The results show almost isotropic cross sections for reactant molecules with high rotational quantum numbers. Backward scattering is disfavoured for reaction out of states with small rotational excitation, in particular the rovibrational ground state. This is a quantum-mechanical effect (interference of partial waves) that can partly be rationalized by simple classical arguments. In particular for higher vibrational excitations, an umbrella effect can be observed that favours the backward direction. It can be explained by the strong enhancement of the reactivity by opening a direct mechanism. The ion-dipole interaction exerts a torque onto the molecule which carries out a rotation by about 90° and then completes the reaction.     

辛准经典轨迹法研究Cl+H2反应

用辛准经典轨迹法模拟了Cl+H2反应在mBW2势能面上的动力学行为。研究了各种初始条件下的反应碰撞截面,产物的能量分配,角度分布和态分布。另外,我们还比较了反应物的三种能量形式(平动能,转动能和振动能)对反庆的有效性。