Computation of cross sections for the F+H2(v=0,j=0) → FH(v′j)+H reaction by the hyperspherical method

Summary A quantum mechanical calculation of cross sections for the reaction F+H₂(v=0,j=0) FH(vj)+H has been performed on the T5A semiempirical potential surface using hyperspherical coordinates. State-to-state integral and differential cross sections converge rapidly with the number of components of the total angular momentum projection onto the axis of least inertia. Thev=3 differential cross section has a forward peak whose magnitude increases with energy whereas thev=2 differential cross section has a backward maximum, in qualitative agreement with cross-beam experiments. Thev=2 andv=3 rotational distributions are in rather good agreement with experiment, but not the vibrational branching ratios.     

Quasi-classical trajectory study of the Ne + H2(+) → NeH(+) + H reaction based on global potential energy surface

Using the multireference configuration interaction method with a Davidson correction and a large orbital basis set (aug-cc-pVQZ), we obtain an energy grid that includes 32 038 points for the construction of a new analytical potential energy surface (APES) for the Ne + H(2)(+) → NeH(+) + H reaction. The APES is represented as a many-body expansion containing 142 parameters, which are fitted from 31?000 ab initio energies using an adaptive nonlinear least-squares algorithm. The geometric characteristics of the reported APES and the one presented here are also compared. On the basis of the APES we obtained, reaction cross sections are computed by means of quasi-classical trajectory (QCT) calculations and compared with the experimental and theoretical data in the literature.     

The polarization dependent differential cross sections of the reactions:H+LiH~+(v=0,j=0)→H_2+Li~+ and H~++ LiH(v=0,j=0)→H_2~++Li

<正>Quasi-classical trajectory(QCT) calculations have been carried out to study the generalized polarization dependent differentialcross sections(PDDCSs) for the reactions H + LiH~+(v = 0,j = 0)→H_2 + Li~+ and H~+ + LiH(v = 0,j = 0)→H_2~+ + Li occurring onthe two lowest-lying electronic states of the LiH_2~+ system,using the ab initio potential energy surfaces(PESs) of Martinazzo et al.[3].Four PDDCSs,i.e.,(2π/σ)(dσ_(00)/dω_t),(2π/σ)(dσ_(20)/dω_t),(2π/σ)(dσ_(22+)/dω_t),(2π/σ)(dσ_(21-)/dω_t) have been discussed in detail.     

Nonadiabatic quantum dynamics of C(1D)+H2→CH+H: coupled-channel calculations including Renner-Teller and Coriolis terms

The Renner-Teller (RT) coupled-channel dynamics for the C((1)D)+H(2)(X(1)Σ(g) (+))→CH(X(2)Π)+H((2)S) reaction has been investigated for the first time, considering the first two singlet states ??(1)A' and b(1)A' of CH(2) dissociating into the products and RT couplings, evaluated through the ab initio matrix elements of the electronic angular momentum. We have obtained initial-state-resolved probabilities, cross sections and thermal rate constants via the real wavepacket method for both coupled electronic states. In contrast to the N((2)D)+H(2)(X(1)Σ(g)(+)) system, RT effects tend to reduce probabilities, cross sections, and rate constants in the low energy range compared to Born-Oppenheimer (BO) ones, due to the presence of a repulsive RT barrier in the effective potentials and to long-lived resonances. Furthermore, contrary to BO results, the rate constants have a positive temperature dependence in the 100-400 K range. The two-state RT rate constant at 300 K, lower than the BO one, remains inside the error bars of the experimental value.     

Differential reaction cross sections in the bending corrected rotating non-linear model: Li + HF → LiF + H

The 3D reaction Li + HF → LiF + H is studied using the bending corrected rotating non-linear model (BCRNM). Differential cross sections are computed and compared to experimental results at two energies, 0.38 and 0.63 eV. Results are found to be in good agreement with experiment at the lower energy, which is near the reaction threshold.     

Communication: state-to-state quantum dynamics study of the OH + CO → H + CO2 reaction in full dimensions (J = 0)

A full dimensional state-to-state quantum dynamics study is carried out for the prototypical complex-formation OH + CO → H + CO(2) reaction in the ground rovibrational initial state on the Lakin-Troya-Schatz-Harding potential energy surface by using the reactant-product decoupling method. With three heavy atoms and deep wells on the reaction path, the reaction represents a huge challenge for accurate quantum dynamics study. This state-to-state calculation is the first such a study on a four-atom reaction other than the H(2) + OH ? H(2)O + H and its isotope analogies. The product CO(2) vibrational and rotational state distributions, and product energy partitioning information are presented for ground initial rovibrational state with the total angular momentum J = 0.     

Ab initio potential energy surface and quantum dynamics for the H + CH4 → H2 + CH3 reaction

A new full-dimensional potential energy surface for the title reaction has been constructed using the modified Shepard interpolation scheme. Energies and derivatives were calculated using the UCCSD(T) method with aug-cc-pVTZ and 6-311++G(3df,2pd) basis sets, respectively. A total number of 30,000 data points were selected from a huge number of molecular configurations sampled by trajectory method. Quantum dynamical calculations showed that the potential energy surface is well converged for the number of data points for collision energy up to 2.5 eV. Total reaction probabilities and integral cross sections were calculated on the present surface, as well as on the ZBB3 and EG-2008 surfaces for the title reaction. Satisfactory agreements were achieved between the present and the ZBB3 potential energy surfaces, indicating we are approaching the final stage to obtain a global potential energy surface of quantitative accuracy for this benchmark polyatomic system. Our calculations also showed that the EG-2008 surface is less accurate than the present and ZBB3 surfaces, particularly in high energy region.     

The reaction Ne+H 2 + (v=0, 1, 2, 3, 4)→NeH++H: 3D potential energy surface and quasiclassical trajectory calculations

A three-dimensional potential energy surface for the endoergic reaction Ne+H ₂ ⁺ →NeH⁺+H in the² A′ ground state of the system NeH ₂ ⁺ has been calculated by quantum chemical ab initio methods (CEPA approximation). The calculated points on this surface were fitted to an analytic ansatz in terms of an extended LEPS functional form augmented by a correction function. The latter was expanded in polynomials in inverse powers of the internuclear distances. This analytic form was used for quasiclassical trajectory calculations of reactive cross sections. In agreement with experimental investigations a strong vibrational enhancement is observed, i.e. the reaction is markedly favored if the necessary reaction energy is supplied as vibrational energy of H ₂ ⁺ rather than as relative translational energy. Other properties of the reaction dynamics such as the backward to forward scattering ratio, the lifetime of the collision complex NeH ₂ ⁺ , and final rotational and vibrational state distributions are also discussed on the basis of the quasiclassical trajectory calculations.     

Semiclassical glory analyses in the time domain for the H + D2(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction

We make the first application of semiclassical (SC) techniques to the plane-wavepacket formulation of time-domain (T-domain) scattering. The angular scattering of the state-to-state reaction, H + D(2)(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D, is analysed, where v and j are vibrational and rotational quantum numbers, respectively. It is proved that the forward-angle scattering in the T-domain, which arises from a delayed mechanism, is an example of a glory. The SC techniques used in the T-domain are: An integral transitional approximation, a semiclassical transitional approximation, a uniform semiclassical approximation (USA), a primitive semiclassical approximation and a classical semiclassical approximation. Nearside-farside (NF) scattering theory is also employed, both partial wave and SC, since a NF analysis provides valuable insights into oscillatory structures present in the full scattering pattern. In addition, we incorporate techniques into the SC theory called "one linear fit" and "two linear fits", which allow the derivative of the quantum deflection function, Θ?(')(J), to be estimated when Θ?J exhibits undulations as a function of J, the total angular momentum variable. The input to our SC analyses is numerical scattering (S) matrix data, calculated from accurate quantum collisional calculations for the Boothroyd-Keogh-Martin-Peterson potential energy surface No. 2, in the energy domain (E-domain), from which accurate S matrix elements in the T-domain are generated. In the E-domain, we introduce a new technique, called "T-to-E domain SC analysis." It half-Fourier transforms the E-domain accurate quantum scattering amplitude to the T-domain, where we carry out a SC analysis; this is followed by an inverse half-Fourier transform of the T-domain SC scattering amplitude back to the E-domain. We demonstrate that T-to-E USA differential cross sections (DCSs) agree well with exact quantum DCSs at forward angles, for energies where a direct USA analysis in the E-domain fails.     

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.     

Determination of the absolute scattering cross section for the reaction D + H2(v=1) → HD+H at 0.33 eV

The reaction D+H₂(v=1) has been investigated in a crossed molecular beam experiment at the most probable collision energy of E₀=0.33 eV. Angular and time-of-flight distributions have been measured and the total absolute cross section has been determined to be σ_r(j_i=0, v = 1, E_c.m. = 0.33 eV) = 1.14 ± 0.50 Ų. This value, as well as the distributions, are in good agreement with the results of quasiclassical trajectory calculations (QCT) and the reactive infinite-order sudden approximation (RIOSA).     

准经典轨线研究Li+HF(v=0,j=O)→LiF+H反应动力学

基于Aguado等人拟合的APW势能面(PES),运用准经典轨线(QCT)方法,对反应Li+HF(v=0,j=0)→LiF+H的动力学性质进行了计算.主要研究了不同碰撞能条件下的反应截面、转动取向、产物散射角分布和竞争反应模式等.结果表明,该反应存在直接提取型和间接插入型两种反应模式,在低能量下反应以间接插入反应模式为主,能量大于200 meV时则以直接提取反应为主.     

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.     

On the fine-structure dependence of the Cs(7p, 2P12,320)+H2→ CsH+H reaction cross section

The experimentally observed effect of the Cs¹ fine-structure level on the title reaction is reproduced by a rough (collinear, one dimensional) dynamical treatment using accurate long-distance entrance channel potentials, as obtained from an effective Hamiltonian method, and hemiquantal dynamics in the collision energy range 0.03 to 0.13 eV. The system is shown to behave adiabatically at large distances, to undergo a transition from adiabatic to diabatic behaviour with respect to spin-orbit coupling in the 11–8.5 au Cs-H distance range, and to transfer adiabatically the whole ²Σ⁺ weight of the electronic wavefunction into theionic channel in the harpooning region.     

Kinetics and dynamics of the NH3 + H → NH2 + H2 reaction using transition state methods, quasi-classical trajectories, and quantum-mechanical scattering

On a recent analytical potential energy surface developed by two of the authors, an exhaustive kinetics study, using variational transition state theory with multidimensional tunneling effect, and dynamics study, using both quasi-classical trajectory and full-dimensional quantum scattering methods, was carried out to understand the reactivity of the NH(3) + H → NH(2) + H(2) gas-phase reaction. Initial state-selected time-dependent wave packet calculations using a full-dimensional model were performed, where the total reaction probabilities were calculated for the initial ground vibrational state and for four excited vibrational states of ammonia. Thermal rate constants were calculated for the temperature range 200-2000 K using the three methods and compared with available experimental data. We found that (a) the total reaction probabilities are very small, (b) the symmetric and asymmetric N-H stretch excitations enhance the reactivity, (c) the quantum-mechanical calculated thermal rate constants are about one order of magnitude smaller than the transition state theory results, which reproduce the experimental evidence, and (d) quasi-classical trajectory calculations, which were performed with the main goal of analyzing the influence of the zero-point energy problem on the final dynamics results, reproduce the quantum scattering calculations on the same surface.     

Time-dependent wave packet theory for state-to-state differential cross sections of four-atom reactions in full dimensions: application to the HD + OH → H2O + D reaction

Time-dependent wave packet method has been developed to calculate differential cross section for four-atom reactions in full dimension, utilizing an improved version of reactant-product-decoupling scheme. Differential cross sections for the title reaction were calculated for collision energy up to 0.4 eV. It is found that the differential cross sections for the reaction are all peaked in the backward direction. The majority of H(2)O is produced in the first stretch excited state, with a large fraction of available energy for the reaction going into H(2)O internal motion. As compared in a previous report by Xiao et al. [Science 333, 440 (2011)], the differential cross section at E(c) = 0.3 eV and the differential cross section at the backward direction as a function of collision energy agree with experiment very well, indicating it is possible now to calculate complete dynamical information for some simple four-atom reactions, as have been done for three-atom reactions in the past decades.     

Differential cross sections for H + D2 → HD(v' = 2, j' = 0,3,6,9) + D at center-of-mass collision energies of 1.25, 1.61, and 1.97 eV

We have measured differential cross sections (DCSs) for the reaction H + D(2) → HD(v' = 2,j' = 0,3,6,9) + D at center-of-mass collision energies E(coll) of 1.25, 1.61, and 1.97 eV using the photoloc technique. The DCSs show a strong dependence on the product rotational quantum number. For the HD(v' = 2,j' = 0) product, the DCS is bimodal but becomes oscillatory as the collision energy is increased. For the other product states, they are dominated by a single peak, which shifts from back to sideward scattering as j' increases, and they are in general less sensitive to changes in the collision energy. The experimental results are compared to quantum mechanical calculations and show good, but not fully quantitative agreement.     

Angular scattering using parameterized S matrix elements for the H + D2(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction: an example of Heisenberg's S matrix programme

A neglected topic in the theory of reactive scattering is the use of parameterized scattering (S) matrix elements to calculate differential cross sections (DCSs). We construct four simple parameterizations, whose moduli are smooth step-functions and whose phases are quadratic functions of the total angular momentum quantum number. Application is made to forward glory scattering in the DCS of the H + D(2)(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction at a translational energy of 1.81 eV, where v and j are vibrational and rotational quantum numbers respectively. The parameterized S matrix elements can reproduce the forward scattering for centre-of-mass reactive scattering angles up to 30° and can identify the total angular momenta (equivalently, impact parameters) that contribute to the glory. The theoretical techniques employed to analyze structure in the DCS include: nearside-farside theory, local angular momentum theory--in both cases incorporating resummations of the partial wave series representation of the scattering amplitude--and the uniform semiclassical theory of forward glory scattering. Our approach is an example of Heisenberg's S matrix programme, in which no potential energy surface is used. Our calculations for the DCS using the four parameterized S matrix elements are counterexamples to the following universal statements often found in the chemical physics literature: "every molecular scattering investigation needs detailed information about the interaction potential," and "an accurate potential energy surface is an essential element in carrying out simulations of a chemical reaction". Both these statements are false.