Quantum stereodynamics of the F+OH(v, j) reactive collision on the 13A″ state

Quantum and quasi-classical calculations for the S⁺ + H₂(v,j) → SH⁺(v',j') + H reactive collisions

State-to-state cross-sections for the S(+) + H2(v,j) → SH(+)(v',j') + H endothermic reaction are obtained using quantum wave packet (WP) and quasi-classical (QCT) methods for different initial ro-vibrational H2(v,j) over a wide range of translation energies. The final state distribution as a function of the initial quantum number is obtained and discussed. Additionally, the effect of the internal excitation of H2 on the reactivity is carefully studied. It appears that energy transfer among modes is very inefficient that vibrational energy is the most favorable for the reaction, and rotational excitation significantly enhances the reactivity when vibrational energy is sufficient to reach the product. Special attention is also paid to an unusual discrepancy between classical and quantum dynamics for low rotational levels while agreement improves with rotational excitation of H2. An interesting resonant behaviour found in WP calculations is also discussed and associated with the existence of roaming classical trajectories that enhance the reactivity of the title reaction. Finally, a comparison with the experimental results of Stowe et al. for S(+) + HD and S(+) + D2 reactions exhibits a reasonably good agreement with those results.

Product rotational angular momentum polarization of H+FCl (v=0-5; j=0, 3, 6, 9) → HF+Cl and HCl+F at Erel=0.5-20 kcal mol(-1)

The rotational angular momentum polarizations of product molecules of the title reactions on the ground potential energy surface 1 (2)A' of DHTSN [Deskevic et al. J Chem Phys 2006, 124, 224303] have been studied using the quasi-classical trajectory method. Reaction dynamic results of the HF product channel comparing with another channel of HCl with 100,000 trajectories can be accurately resolved. We show the value of the polar p(ϑr) in the range of 0° ≤ ϑr ≤ 180(°), azimuthal p(φr) in the range of 0° ≤ φr ≤ 360(°), and dihedral p(ϑr, φr) in the ranges of 0(°) ≤ ϑr ≤ 180(°) and 0(°) ≤ φr ≤ 360(°); the angular distributions of the product molecules HF and HCl at relative Erel = 0.5, 1, 2, 5, 10, 15, and 20 kcal mol(-1); and four polarization-dependent differential cross sections (PDDCSs) of HF and HCl at Erel = 0.5, 1, 2, 5, 10, and 15 kcal mol(-1). p(φr) distributions at v = 0-5, and j = 0, 3, 6, 9 at every Erel are plotted cylindrically together. The stereo dynamic transformation reaction dependent upon the rovibrational states of the reactant molecule FCl and its relative translational energies around 0.5-5 kcal mol(-1) can be significantly differentiated. Translational and rovibrational enhancements of the title reactions on both early barrier potential energy surfaces have been shown in great detail and clarified. Reaction mechanisms of forward and backward scattering of the product molecules HF and HCl, respectively, have been obtained. Graphical Abstract H + FCl → either HF + Cl (left) or HCl + F (right) is moving along a trajectory on the respective PES.

Wave packet calculations on nonadiabatic effects for the O(3P)+HF(1Σ+) reaction under hyperthermal conditions

We present wave packet calculations of total and state-to-state reaction probabilities and integral cross sections for the nonadiabatic dynamics of the O((3)P)+HF → F((2)P)+OH((2)Π) reaction at hyperthermal collision energies ranging from 1.2 to 2.4 eV. The validity of the centrifugal sudden approximation is discussed for the title reaction and a comprehensive investigation of the influence of nonadiabatic effects on the dynamics of this reactive system at high (hyperthermal) collision energies is presented. In general, nonadiabatic effects are negligible for averaged observables, such as total reaction probabilities and integral cross sections, but they are clearly observed in detailed observables such as rotationally state-resolved reaction probabilities. A critical discussion of nonadiabatic effects on the dynamics of the title reaction is carried out by comparing with the reverse reaction and the characteristics of the adiabatic and diabatic potential energy surfaces involved.

EFFECT OF THE REACTANT ALIGNMENT ON STEREODYNAMICS FOR THE REACTION F + OH

Quasi-classical trajectory (QCT) calculations are performed for the reaction F + OH → HF + O based on the adiabatic potential-energy surface (PES) of the 13A″ triplet state. The reaction probability as a function of incident angle has been presented. The differential cross sections (DCSs), the distribution of angle between k and j′, P(θr) and the distribution of dihedral angle denoting k – k′ – j′ correlation, P(ϕr) have also been calculated at the different incident angles in the center-of-mass (CM) frame, respectively. The effects of reactant alignment on stereodynamics of the reaction are firstly revealed, which provides the theoretical foundation for the related experiment and enriches the theories of the stereodynamics.

STATE-TO-STATE REACTION DYNAMICS OF THE D(H) + FO → OD(OH) + F REACTIONS ON THE LOWEST TWO TRIPLET STATES

Quasiclassical trajectory (QCT) calculations have been performed to investigate the reaction dynamics of the title reactions at a state-to-state level. The recently developed adiabatic potential energy surfaces (PESs) of the two triplet electronic states of 13A″ and 13A′ (Gómez-Carrasco S. et al., J Chem Phys121:4605, 2004; J Chem Phys123:114310, 2005) are employed in the present QCT calculations. Product vibrational and rotational state distributions have been calculated at three collision energies of 0.5, 0.7 and 0.9 eV. The product vibrational state distributions are found to be Gaussian distributions. Both the vibrational and rotational state distributions depend clearly on the isotope substitution, electronic PESs and collision energies. The observed phenomena are probably attributed to the competition of direct and indirect reactions.

Quasi-classical trajectory study of the O(1D) + HF reaction dynamics on 11A′ potential energy surface

The quasi-classical trajectory (QCT) calculations for the title reaction were carried out using the recently developed, accurate potential energy surface (PES) of the [Formula: see text] singlet state of the OHF system The integral cross section and the product rotational alignment factor [Formula: see text] were calculated as a function of collision energy. In addition, I discovered the effect of isotopic substitution on stereodynamics for the title reaction, and the influence of the rotation excitation of the reagent on stereodynamics is also presented. Both the scalar and vector properties of the reaction O(1D) + HF → OH + F(2P) are studied in this paper. It was found that the reaction is mainly controlled by an indirect reaction mechanism, and that the deep noncollinear insertion HOF well has a great impact on the dynamics of the reaction. The conclusions drawn in this paper will draw from references to similar reactions, and provide a theoretical foundation for related experiments.

Product rotational angular momentum polarization in the H+FCl(v=0-5, j=0, 3, 6, 9)→HF+Cl reaction

The product alignment and orientation of the title reaction on the ground potential energy surface of 1 (2)A' have been studied using the quasi-classical trajectory method. The calculations were carried out for case (a) at collision energies of 0.5-20 kcal mol(-1) with the initially rovibrational state of the reagent FCl molecule being at the v = 0 and j = 0 level to especially reveal in detail the dependence of the product integral cross section on collision energy. Further calculations at the collision energy of 15 kcal mol(-1) for case (b) at v = 0-5, and j = 0, and (c) at v = 0, and j = 3, 6, 9 initial states were carried out to reveal the effect of initially vibrational and rotational excitations on stereodynamics, respectively. Possessing final relative velocity k' (defined as a vector in the xz-plane), product alignment perpendicular to the reagent relative velocity vector k (defined as z- or parallel to the z-axis), for case (a) is found to be weaker at all collision energies, for case (b) is found to be vibrationally enhanced by the reactant molecule FCl, but for case (c), rather insensitive to initially rotational excitation. The rotational vector of product molecular orientation pointing to either negative or positive direction of the y-axis in the center of mass frame, e.g. origin of the coordinate system, is enhanced by collision energies regarding to 0.5-20 kcal mol(-1), while it becomes weaker at higher vibrational (v = 0-5) or rotational (j = 0, 3, 6, 9) excitation levels. Effects of collision energies and of rotational excitation at these collision energies, with 15 kcal mol(-1) as an example on the calculated PDDCSs are also shown and discussed. Detailed plots P(φ(r)) in the range of 0 ≤φ(r)≤ 360(o), and P(θ(r), φ(r)) in the ranges of 0 ≤θ(r)≤ 180° and 0 ≤φ(r)≤ 360° at collision energies 0.5-20 kcal mol(-1) have been presented. Overall, results of PDDCSs of the product alignment and product orientation at these collision energies in the title reaction are not very strongly distinguishable.

Quantum mechanics and quasiclassical study of the H/D+FO --> OH/OD+F, HF/DF+O reactions: chemical stereodynamics

The time-dependent quantum wave packet and the quasi-classical trajectory (QCT) calculations for the title reactions are carried out using three recent-developed accurate potential energy surfaces of the 1(1)A', 1(3)A', and 1(3)A'' states. The two commonly used polarization-dependent differential cross sections, dsigma(00)/domega(t), dsigma(20)/domega(t), with omega(t) being the polar coordinates of the product velocity omega', and the three angular distributions, P(theta(r)), P(Phi(r)), and P(theta(r),Phi(r)), with theta(r), Phi(r) being the polar angles of the product angular momentum, are generated in the center-of-mass frame using the QCT method to gain insight into the alignment and the orientation of the product molecules. Influences of the potential energy surface, the collision energy, and the isotope mass on the stereodynamics are shown and discussed. Validity of the QCT calculation has been examined and proved in the comparison with the quantum wave packet calculation.

F+OH reactive collisions on new excited A″3 and A′3 potential-energy surfaces

Global three-dimensional adiabatic potential-energy surfaces for the excited 2(3)A" and 1(3)A' triplet states of OHF are obtained to study the F(2P)+OH(2pi)-->O(3P)+HF(1sigma+) reaction. Highly accurate ab initio calculations are obtained for the two excited electronic states and fitted to analytical functions with small deviations. The reaction dynamics is studied using a wave-packet treatment within a centrifugal sudden approach, which is justified by the linear transition state of the two electronic states studied. The reaction efficiency presents a marked preference for perpendicular orientation of the initial relative velocity vector and the angular momentum of the OH reagent, consistent in the body-fixed frame used with an initial collinear geometry which facilitates the access to the transition state. It is also found that the reaction cross section presents a rather high threshold so that, in an adiabatic picture, the two excited triplet states do not contribute to the rate constant at room temperature. Thus, only the lowest triplet state leads to reaction under these conditions and the simulated rate constants are too low as compared with the experimental ones. Such disagreement is likely to be due to nonadiabatic transitions occurring at the conical intersections near the transition state for this reaction.