The Reaction X + Cl2→XCl + Cl (X = Mu,H, D). II. Comparison of experimental data with theoretical results derived from a new potential energy surface

We consider experimental implications for the Mu + Cl₂, H + Cl₂, and D + Cl₂ reactions of the extended London—Eyring—Polanyi—Sato (LEPS) potential energy surface derived from experimental data in paper I. In the present calculations, it is necessary to make additional implicit and explicit assumptions concerning the three-dimensional (3D) nature of the potential surface, since the inversion procedure of paper I yields information only on the collinear (1D) part of the surface. We have performed accurate 1D quantum calculations of reaction probabilities, which are then transformed into 3D by an information theoretic 1D → 3D transformation incorporating a constraint to allow for angular momentum transfer effects in light+heavy—heavy atom reactions. This procedure implicitly accounts for the 3D nature of the potential surface. The calculated vibrational and vibrotational product distributions are in good agreement with those determined in thermal chemiluminescence experiments. The Sato parameters for the 1D surface also define a full 3D surface. This is used as an approximation to the true surface, and its properties are explored in 3D quasiclassical trajectory calculations. Comparison is made for the H and D reactions with available chemiluminescence, molecular beam and kinetic experimental data for differential and total reaction cross sections, energy disposal, rate coefficients and Arrhenius parameters. Some kinetic isotope effects in the Mu, H, and D reactions are discussed using vibrationally adiabatic theory. Comparison is also made with results from other calculations in the literature for the H + Cl₂ and D + Cl₂ reactions.     

The effect of reagent translational energy in the reaction O(D2) + H2 → OH(Π) + H

Quasiclassical trajectory calculations have been performed to determine the effect of reactant collision energy on product state distributions in the reaction O(¹D) + H₂ → OH(²Π) + H. The product vibrational distribution becomes more excited as the collision energy is increased. This is not due to an increase in the cross section for collinear abstraction. A detailed analysis has shown that strong O---H₂ repulsion, which occurs during the insertion of the O into the H---H bond, converts the kinetic energy of the reacting system to vibrational motion of the intermediate.     

Quantum stereodynamics of Li + HF reactive collisions: the role of reactants polarization on the differential cross section

A complete quantum study for the state-to-state Li + HF(v,j,m) → LiF(v',j',Ω') + H reactive collisions has been performed using a wave packet method, for different initial rotational states and helicity states of the reactants. The state-to-state differential cross section has been simulated, and the polarization of products extracted. It is found that the reactivity is enhanced for nearly collinear collisions, which produces a vibrational excitation of HF, needed to overcome the late barrier. It is also found that LiF(v' = 0) products are preferentially forward scattered, while vibrationally excited LiF(v' = 1 and 2) are backward scattered. These results are interpreted with a simple reaction mechanism, based on the late character and bent geometry of the transition state, originating from a covalent/ionic crossing, which consists of two steps: the arrival at the transition state and the dissociation. In the first step, in order to get to the saddle point some HF vibrational excitation is required, which favors head-on collisions and therefore low values of m. In the second step a fast dissociation of H atom takes place, which is explained by the ionic Li(+)F(-)H character of the bent transition state: the FH(-) is repulsive making that H depart rapidly leaving a highly rotating LiF molecule. For the higher energy analyzed, where resonances slightly contribute, the orientation and alignment of product rotational states, referred to as reactants frame (with the z-axis parallel to k), are approximately constant with the scattering angle. The alignment is close to -1, showing that j' is perpendicular to k, while starting from initial states with well defined rotational orientation, as states with pure m values, the final rotational are also oriented. It is also found that when using products frame (with the z'-axis parallel to k') the rotational alignment and orientation of products varies a lot with the scattering angle just because the z' axis changes from being parallel to anti-parallel to k when varying from θ = 0 to π.     

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.     

Internal Energy State Distribution of Product CaBr in the Collision Reactions Ca+C2H5Br和Ca+nC3H7Br

The internal energy distributions of product CaBr in the collision reactions Ca+C2H5Br and Ca+nC3H7Br are studied by using the quasiclassical trajectory method. The average vibrational, rotational and translational energies and total available energies of the product CaBr molecules are calculated. The results indicate that when the collision energy is equal to 7.54 kJ/mol the energy of product CaBr is mainly the vibrational energy. As the reactant collision energy increases, the average translational and rotational energies of the product CaBr increase, the average vibrational energy decreases slightly, and the most probable vibrational state shifts to lower vibrational energy levels. The internal states of reagents have little influence on the internal energy distribution of the product. The bigger the radical group is, the higher ratio of the vibrational energy to the available energy of the product is. There exist two competitive reaction paths for the collision reactions Ca+C2H5Br and Ca+nC3H7Br, the migratory encounter and direct reaction paths. The former produces high vibrational excited state product CaBr and the latter causes C-Br bond to break. When the collision energy increases, the reactions tend to the latter path.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

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

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

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.     

High Resolution Crossed Beams Scattering Study of the F＋HD→DF＋H Reaction

Crossed beams scattering study was carried out on the F＋HD→DF＋H reaction using high- resolution H-atom Rydberg tagging time-of-flight technique. Vibrational state-resolved differential cross sections were measured, with partial rotational state resolution, at eight collision energies in the range of 2.51-5.60 kJ/mol. Experimental results indicated that the product angular distributions are predominantly backward scattered. As the collision energy increases, the backward scattered peak becomes broader gradually. Dependence of product vibration branching ratios on the collision energy was also determined. The experimental results show that the DF products are highly inverted in the vibrational state distribution and the DF （v＇=3） product is the most populated state. Furthermore, the DF （v＇=l） product has also been observed at collision energy above 3.97 kJ/mol.     

Influence of the potential energy surface on the reaction cross section: the K + HF → KF + H system

The effect of the potential energy surface on the K + HF → KF + H cross section has been studied using reasonable Sorbie-Murrell (bent saddle point) and LEPS (collinear saddle point) potential energy surfaces (PESs). Trajectory calculations for selected initial conditions (translational energies, rovibrational levels (v, J) of HF, as well as initial parallel or perpendicular alignments between the HF rotational angular momentum and the reactants relative velocity vectors) have been performed on these PESs to compare them with experiments. The Sorbie-Murrell and LEPS-4 PESs lead to steric effect ratio results quite close to the experimental ones, once the error margins are included. The resllts point towards a bent K-F-H saddle point although the PES is very isotropic. This could explain why experimental determinations lead to suggest a collinear saddle point. The K + HF → KF + H reaction exhibits an enormous vibrational enhancement of reactivity with one quantum HF vibrational excitation, even at translational energies well above the HF(v=0) threshold, where tunnelling effect contribution to reactivity can be neglected. This behaviour has not been reproduced in the trajectory calculations and no satisfactory explanation has been obtained for this fact. Nevertheless, the HF(v=1)/HF(v=0) cross section ratio at translational energies not far from the HF(v=0) threshold and the relative cross section for HF(v=0) have been satisfactorily descibed. In what regards rotation, the best theoretical results are those corresponding to the Sorbie-Murrell PES (the cross section increases with J), although important differences with experiment appear for the J = 0–3 interval at the lower translational energy values considered (0.54 and 0.77 eV).     

Cross sections of the O+ + H2 --> OH+ + H ion-molecule reaction and isotopic variants (D2, HD): quasiclassical trajectory study and comparison with experiments

A dynamics study [cross section and microscopic mechanism versus collision energy (E(T))] of the reaction O+ + H2 --> OH+ + H, which plays an important role in Earth's ionosphere and interstellar chemistry, was conducted using the quasiclassical trajectory method, employing an analytical potential energy surface (PES) recently derived by our group [R. Martinez et al., J. Chem. Phys. 120, 4705 (2004)]. Experimental excitation functions for the title reaction, as well as its isotopic variants with D2 and HD, were near-quantitatively reproduced in the calculations in the very broad collision energy range explored (E(T) = 0.01-6.0 eV). Intramolecular and intermolecular isotopic effects were also examined, yielding data in good agreement with experimental results. The reaction occurs via two microscopic mechanisms (direct and nondirect abstraction). The results were satisfactorily interpreted based on the reaction probability and the maximum impact parameter dependences with E(T), and considering the influence of the collinear [OHH]+ absolute minimum of the PES on the evolution from reactants to products. The agreement between theory and experiment suggests that the reaction mainly occurs through the lowest energy PES and nonadiabatic processes are not very important in the wide collision energy range analyzed. Hence, the PES used to describe this reaction is suitable for both kinetics and dynamics studies.     

Effect of vibrational excitation of HBr on the H+HBr abstraction and exchange reactions

The effect of vibrational excitation of HBr on the H+HBrH₂+Br and H+HBrH+HBr reactions has been investigated on the extended LEPS surface (ELEPS) constructed on the basis of quantum chemically calculated points of PES. Together with this surface the LEPS surface of Sudhakaran and Raff [1] was used for comparison at two relative translational energies. A quasiclassical trajectory method was used to study the abstraction and exchange reaction dynamics. The reactive cross section was calculated as a function of the relative collision energy and the vibrational state of HBr. The following conclusions can be drawn from the results of the study: (i) vibrational excitation v=0 v=2 more than doubles the reaction cross section, (ii) the increase in the collision energy is most effectively channelled into the product translational energy.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

Stereodynamics of the Ca + HCl → CaCl + H molecular reaction imposed by the rotational-excited states of HCl

The effects of the initial rotational-excited states of the HCl molecule on the stereodynamics properties of the Ca + HCl molecular reaction are investigated using the quasiclassical trajectory theory and the analytical potential energy surface. The orientation and alignment behaviors for the rotational angular momentum of the product, along with the generalized differential cross section (PDDCS)-dependent polarization, are calculated to explore the stereodynamics properties. The initial rotational-excited states of the HCl molecule impose a remarkable effect on the vector correlation distributions, regardless of the orientation, alignment, or PDDCS. The forward, backward, and weak sideway scatterings are found in the Ca + HCl → CaCl + H molecular reaction. The results demonstrate that the initial rotational-excited state of j = 3 results in more obvious stereodynamics effects.     

Accurate quantum mechanical calculations of differential and integral cross sections and rate constant for the O+OH reaction using an ab initio potential energy surface

The authors report accurate quantum mechanical studies of the O+OH reaction on the improved Xu-Xie-Zhang-Lin-Guo potential energy surface. The differential cross section was obtained at several energies near the reaction threshold using a time-independent method. The dominant forward and backward peaks in the angular distribution are consistent with a complex-forming mechanism, which is also confirmed by the extensive rotational excitation in the O2 product. However, the asymmetry of these peaks suggests a significant nonstatistical component. The initial state (upsilon i=0, j i=0) specified integral cross section, which was calculated up to 1.15 eV of collision energy using the Chebyshev wave packet method, shows no energy threshold and decreases with the increasing collision energy, consistent with the barrierless nature of the reaction. The resulting rate constant exhibits a negative temperature dependence for T>100 K and decays as the temperature is lowered, in qualitative agreement with available experimental data.     

Quasiclassical trajectory study of the F + CH4 reaction dynamics on a dual-level interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the F + CH4 reactive system has been constructed using the interpolation method of Collins and co-workers. The ab initio calculations have been performed using second-order M?ller-Plesset (MP2) perturbation theory to build the initial PES. Scaling all correlation (SAC) methodology has been employed to improve the ab initio calculations and to construct a dual-level PES. Using this PES, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations and internal energy distributions has been carried out for the F + CH4 and F + CD4 reactions and the theoretical results have been compared with the available experimental data.     

Quasi-Classical Trajectory Study of the Chemical Reaction Ca+CH3I

The Ca＋CH3I→CaI＋CH3 reaction system has been studied with the quasi-classical trajectory method on the extended Lond-Eyring-Polanyi-Sato（LEPS） potential energy surface. At collision energy Ecol=10.78 kJ/mol, the calculated results show that the CaI vibrational population peaks are located at v=2. The calculated cross section decreases slowly with the collision energy increasing. The angle product distributions tend toward backward scattering. The calculated （P2（J^1·K）） values deviate slightly from-0.5 and decrease with increasing collision energy. The Quasiclassical trajectory calculation（QCT） results are in reasonable agreement with experimental data. Moreover, the dynamics of the reaction has been discussed.     

Dynamics of the D+ + H2 and H+ + D2 reactions: a detailed comparison between theory and experiment

An extensive set of experimental measurements on the dynamics of the H(+) + D(2) and D(+) + H(2) ion-molecule reactions is compared with the results of quantum mechanical (QM), quasiclassical trajectory (QCT), and statistical quasiclassical trajectory (SQCT) calculations. The dynamical observables considered include specific rate coefficients as a function of the translational energy, E(T), thermal rate coefficients in the 100-500 K temperature range. In addition, kinetic energy spectra (KES) of the D(+) ions reactively scattered in H(+) + D(2) collisions are also presented for translational energies between 0.4 eV and 2.0 eV. For the two reactions, the best global agreement between experiment and theory over the whole energy range corresponds to the QCT calculations using a gaussian binning (GB) procedure, which gives more weight to trajectories whose product vibrational action is closer to the actual integer QM values. The QM calculations also perform well, although somewhat worse over the more limited range of translational energies where they are available (E(T) < 0.6 eV and E(T) < 0.2 eV for the H(+) + D(2) and D(+) + H(2) reactions, respectively). The worst agreement is obtained with the SQCT method, which is only adequate for low translational energies. The comparison between theory and experiment also suggests that the most reliable rate coefficient measurements are those obtained with the merged beams technique. It is worth noting that none of the theoretical approaches can account satisfactorily for the experimental specific rate coefficients of H(+) + D(2) for E(T)≤ 0.2 eV although there is a considerable scatter in the existing measurements. On the whole, the best agreement with the experimental laboratory KES is obtained with the simulations carried out using the state resolved differential cross sections (DCSs) calculated with the QCT-GB method, which seems to account for most of the observed features. In contrast, the simulations with the SQCT data predict kinetic energy spectra (KES) considerably cooler than those experimentally determined.     

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.     

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.     

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

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