A global potential energy surface for H2S+(X 4A′′) and quasi-classical trajectory study of the S+(4S) + H2(X1Σ+g) reaction

A novel global potential energy surface for H₂S⁺(X?⁴A″) based on accurate ab initio calculations is presented. Energies are calculated at the multi-reference configuration interaction level with Davidson correction using aug-cc-pVQZ basis set plus core-polarisation high-exponent d functions. A grid of 4552 points is used for the least-square fitting procedure in the frame of a many-body expansion. The topographical features of the new potential energy surface are here discussed in detail. Such a surface is then employed for dynamic studies of the S(⁴S) + H₂(X?¹Σ⁺_g) →SH⁺(X?³Σ^?) + H(²S) reaction using the quasi-classical trajectory method. State specific trajectories are calculated, for both ground and ro-vibrationally excited initial states of H₂(X?¹Σ⁺_g). Corrections to the zero point energy leakage of the classical calculations are also presented. Calculated reaction cross sections and rate constants are here reported and compared with available literature.     

Mechanism analysis of reaction S~+(~2D)+H_2(X~1Σ_g~+) → SH~+(X~3Σ~-) + H(~2S) based on the quantum state-to-state dynamics

We present a state-to-state dynamical calculation on the reaction S~++ H_2→ SH~+ +H based on an accurate ~X2 A~″ potential surface. Some reaction properties, such as reaction probability, integral cross sections, product distribution, etc.,are found to be those with characteristics of an indirect reaction. The oscillating structures appearing in reaction probability versus collision energy are considered to be the consequence of the deep potential well in the reaction. The comparison of the present total integral cross sections with the previous quasi-classical trajectory results shows that the quantum effect is more important at low collision energies. In addition, the quantum number inversion in the rotational distribution of the product is regarded as the result of the heavy–light–light mass combination, which is not effective for the vibrational excitation. For the collision energies considered, the product differential cross sections of the title reaction are mainly concentrated in the forward and backward regions, which suggests that there is a long-life intermediate complex in the reaction process.     

Ab initio dynamics of the He + H+ 2 → HeH+ +H reaction: a new potential energy surface and quantum mechanical cross-sections

The reaction He + H⁺ ₂(v,j = 0) → HeH⁺(v′ = 0, j′) for v = 0, 1,2 and 3 and for scattering energies near the threshold (0.95–1.15 eV) has been studied by calculating ab initio points at MRCI level and ‘exact’ integral quantum reactive cross-sections. More than 1400 nuclear geometries have been chosen to cover the most important regions for the dynamics, an extended set of points being taken directly on a hyperspherical coordinate grid. A many-body expansion with a large number of terms permits an accurate analytical representation of the potential energy surface with a root-mean-square deviation <12meV. The hyperquantization algorithm has been extended to obtain quantum mechanical integral cross-sections which are compared with previous calculations and with experimental results.     

A new global potential energy surface of the ground state of SiH2+ (X2A1) system and dynamics calculations of the Si+ + H2 (v0 = 2, j0 = 0) → SiH+ + H reaction

A global potential energy surface (PES) of the ground state of SiH$_{2}^{+}$ system is built by using neural network method based on 18223 ab initio points. The topographic properties of PES are presented and compared with previous theoretical and experimental studies. The results indicate that the spectroscopic parameters obtained from the new PES are in good agreement with the experimental data. In order to further verify the validity of the new PES, a test dynamics calculation of the Si$^{+} +$ H$_{2}$ ($v_0 = 2, j_{0} = 0$) $\to $ H $+$ SiH$^{+}$ reaction has been carried out by using the time-dependent wave packet method. The integral cross sections and rate constants are computed for the title reaction. The reasonable dynamical behavior indicates that the newly constructed PES is suitable for relevant dynamics investigations.     

State-to-state quantum dynamics of the N(~4S)+H_2(X~1Σ~+) → NH(X~3Σ~-)+H(~2S) reaction and its reaction mechanism analysis

Quantum state-to-state dynamics of the N(4S) + H-2(X1+Σ) → NH(X3Σ) + H(2S) reaction is reported in an accurate novel potential energy surface constructed by Zhai et al.(2011 J. Chem. Phys. 135 104314). The time-dependent wave packet method, which is implemented on graphics processing units, is used to calculate the differential cross sections. The influences of the collision energy on the product state-resolved integral cross sections and total differential cross sections are calculated and discussed. It is found that the products NH are predominated by the backward scattering due to the small impact parameter collisions, with only minor components being forward and sideways scattered, and have an inverted rotational distribution and no inversion in vibrational distributions; both rebound and stripping mechanisms exist in the case of high collision energies.     

Quasi-classical trajectory study of the isotope effect on the stereodynamics in the reaction H(2S)+CH(X2Π; v=0, j=1)→C(1D)+H2(X1Σg+)

The isotope effect on the stereodynamic properties in the title reaction is investigated by a quasi-classical trajectory (QCT) method on the 1¹A' potential energy surface at a collision energy of 23.06 kcal/mol. The angular distributions P(θ_r), P(ø_r), P(θ_r, ø_r), and the polarization-dependent generalized differential cross sections are calculated, which demonstrate the observable influences on the rotational polarization of the product by the isotopic substitution of H with D.     

Study of the H＋HS reaction on a newly built potential energy surface using the quasi-classical trajectory method

<正>The quasi-classical trajectory（QCT） method is used to study the H+HS reaction on a newly built potential energy surface（PES） of the triplet state of H₂S（³A″） in a collision energy range of 0-60 kcal/mol.Both scalar properties, such as the reaction probability and the integral cross section（ICS）,and the vector properties,such as the angular distribution between the relative velocity vector of the reactant and that of the product,etc.,are investigated using the QCT method.It is found that the ICSs obtained by the QCT method and the quantum mechanical（QM） method accord well with each other.In addition,the distribution for the product vibrational states is cold,while that for the product rotational states is hot for both reaction channels in the whole energy range studied here.     

State-to-state dynamics of F(~2P) + HO(~2Π) → O(~3P) + HF(~1Σ~+)reaction on 1~3A〞 potential energy surface

State-to-state time-dependent quantum dynamics calculations are carried out to study F(~2P) + HO(~2Π) → O(~3P) +HF(~1Σ~+) reaction on 1~3A〞 ground potential energy surface(PES). The vibrationally resolved reaction probabilities and the total integral cross section agree well with the previous results. Due to the heavy–light–heavy(HLH) system and the large exoergicity, the obvious vibrational inversion is found in a state-resolved integral cross section. The total differential cross section is found to be forward–backward scattering biased with strong oscillations at energy lower than a threshold of 0.10 eV, which is the indication of the indirect complex-forming mechanism. When the collision energy increases to greater than 0.10 eV, the angular distribution of the product becomes a strong forward scattering, and almost all the products are distributed at θ_t = 0°. This forward-peaked distribution can be attributed to the larger J partial waves and the property of the F atom itself, which make this reaction a direct abstraction process. The state-resolved differential cross sections are basically forward-backward symmetric for v' = 0, 1, and 2 at a collision energy of 0.07 eV; for a collision energy of 0.30 eV,it changes from backward/sideward scattering to forward peaked as v increasing from 0 to 3. These results indicate that the contribution of differential cross sections with more highly vibrational excited states to the total differential cross sections is principal, which further verifies the vibrational inversion in the products.     

A theoretical study of energy transfer in the quenching of O(1 D 2) by CO(1Σ+)

Energy transfer processes occurring during the collisional quenching of O(¹ D) by CO(¹Σ⁺) are studied using a classical collision complex model together with potentials previously derived for the C(³ P) + O₂(³Σ _g ^-) reaction. Room temperature quenching rate constants, electronic-vibrational transfer efficiencies and product vibrational state distributions are in good agreement with experiment. The calculated and experimental temperature dependences of the electronic-vibrational transfer efficiencies and vibrational populations, however, disagree. The effect of using vibrationally excited CO as a quenching partner is studied and shown to result in a lowering of the quenching rate constant by a factor of 4 at room temperature. Enhancement of initial translational energy by the equivalent of a vibrational quantum of energy leads to an even larger decrease in the rate constant. This difference between vibrational and translational energy enhancement is interpreted in terms of an increased centrifugal barrier in the latter case.     

Quasi-classical trajectory study of H+LiH(v=0,1,2,j=0)→Li+H_2 reaction on a new global potential energy surface

Quasi-classical trajectory(QCT) calculations are reported for the H+LiH(v = 0–2, j = 0)→Li+H_2 reaction on a new ground electronic state global potential energy surface(PES) of the LiH_2 system. Reaction probability and integral cross sections(ICSs) are calculated for collision energies in the range of 0 eV–0.5 eV. Reasonable agreement is found in the comparison between present results and previous available theoretical results. We carried out statistical analyses with all the trajectories and found two main distinct reaction mechanisms in the collision process, in which the stripping mechanism(i.e., without roaming process) is dominated over the collision energy range. The polarization dependent differential cross sections(PDDCSs) indicate that forward scattering dominates the reaction due to the dominated mechanism. Furthermore,the reactant vibration leads to a reduction of the reactivity because of the barrierless and attractive features of PES and mass combination of the system.     

Direct mapping of insertion reaction dynamics: S(1D)+H2→SH+H

The Doppler-selected time-of-flight method was applied to map out the differential cross sections of the title reaction at two collision energies in a crossed-beam experiment. Roughly symmetric and highly forward–backward peaking angular distributions were observed at both energies. Vibrational structures of the SH product were resolved from the velocity measurements of the counter-product H-atom. Most of the angle-integrated observables can readily be understood on statistical grounds, which suggests that statistics plays the dominant role in determining the outcomes of this prototypical insertion reaction. In terms of more detailed angle-specific reaction attributes, significant discrepancies from statistical considerations were revealed, indicative of some hidden dynamics being buried under the statistical factor. Received: 26 February 2000 / Revised version: 20 April 2000 / Published online: 13 September 2000     

State-to-state quantum dynamics studies of the H + LiH+ → Li+ + H2 reaction

ABSTRACT

Quantum dynamical calculations of the H?+?LiH⁺?→?Li⁺?+?H₂ reaction were performed based on the potential energy surface (PES) reported by Dong et al. (RSC Adv. 7, 7008 (2017)) using the time-dependent quantum wave packet method in collision energy range from 0.01 to 1.0?eV. Dynamics properties such as reaction probability, integral cross section, differential cross section (DCS), and thermal rate constant of the H?+?LiH⁺?→?Li⁺?+?H₂ reaction were reported at the state-to-state level of theory and compared with available theoretical calculations. The results indicated that present values are in good agreement with results obtained from the quasi-classical trajectory method. However, large differences can be found between present values and previous quantum results. This can be attributed to the different PESs used in the calculation and the CS approximation was adopted in previous theoretical studies. In addition, the ‘rebound’ reaction mechanism was proposed in previous theoretical studies in a high collision energy range. However, the DCS scattering signals calculated in the present work indicated that complex-forming and direct abstract reaction mechanisms are dominant in low and high collision energies, respectively.     

Vibrational and rotational excitation effects of the N(2D) + D2(X1Σg +) → ND(X3Σ+) + D(2S) reaction

The effects of the rovibrational excitation of reactants in the N(²D) + D₂(X¹Σ_g⁺) → ND(X³Σ⁺) + D(²S) reaction are calculated in a collision energy range from the threshold to 1.0 eV using the time-dependent wave packet approach and a second-order split operator. The reaction probability, integral cross-section, differential cross-section and rate constant of the title reaction are calculated. The integral cross-section and rate constant of the initial states v = 0, j = 0, 1, are in good agreement with experimental data available in the literature. The rotational excitation of the D₂ molecule has little effect on reaction probability, integral cross-section and the rate constant, but it increased the sideways and forward scattering signals. The vibrational excitation of the D₂ molecule reduced the threshold and broke up the forward–backward symmetry of the differential cross-section; it also increased the forward scattering signals. This may be because the vibrational excitation of the D₂ molecule reduced the lifetime of the intermediate complex.     

Influence of ro-vibrational and isotope effects on the dynamics of the C(3 P)+ OD(X 2Π) → CO(X 1 Σ+) + D(2 S) reaction

The C(³ P)+OD(X ²Π) reaction has been studied by means of quantum mechanical real wave packet (RWP) and quasiclassical trajectory (QCT) methodologies on the ground potential energy surface of Zanchet et al. [J. Phys. Chem. A 110, 12017 (2006)]. Initial state selected total reaction probabilities at J?=?0 total angular momentum have been calculated for a wide range of collision energies. Product state-resolved integral cross-sections at selected collision energies and excitation functions have been determined from the RWP calculations using the J-shifting approximation and from QCT calculations. State-specific and thermal rate coefficients have been calculated using both methodologies up to 500 K. The effect of reagent rotational excitation on the dynamics for the C(³ P)+OH(X ²Π) and C(³ P)+OD(X ²Π) reactions has been investigated and interesting discrepancies between the QCT and RWP results have been found. The RWP results are found to be in an overall good agreement with the corresponding QCT results, although the QCT integral cross-section and rate coefficients are slightly smaller than those obtained from the RWP calculations.     

Novel potential energy surface-based quantum dynamics of ion–molecule reaction O~++ D_2→OD~++ D

According to a novel electronic ground-state potential energy surface of H_2O~+(X~4 A~'),we calculate the reaction probabilities and the integral cross section for the titled reaction O+~+ D_2→OD~++ D by the Chebyshev wave packet propagation method.The reaction probabilities in a collision-energy range of 0.0 e V–1.0 e V show an oscillatory structure for the O~++ D_2 reaction due to the existence of the potential well.Compared with the results of Mart′?nez et al.,the present integral cross section is large,which is in line with experimental data.     

Quasi-classical trajectory approach to the O（1D）＋HBr ）OH＋Br reaction stereo-dynamics on X1A＇ potential energy surface

The vector properties of reaction O(¹D)+HBr→ OH+Br on the potential energy surface (PES) of X¹A′ ground singlet state are studied by using the quasi-classical trajectory (QCT) theory. The polarization-dependent differential cross sections (PDDCSs), the average rotational alignment factor 2(j′· k)>, as well as the distributions reflecting vector correlations are also computed. The analysis of the results shows that the alignment and the orientation distribution of the rotation angular momentum vector of product molecule OH is influenced by both the effect of heavy-light-heavy (HLH) type mass combination and the deep well of PES.     

Time-dependent quantum wave packet and quasiclassical trajectory studies of the Au(2S) + H2(X1∑+g) → AuH(X1∑+g) + H(2S) reaction

The time-dependent quantum wave packet (TDWP) and quasiclassical trajectory calculations (QCT) are carried out for the Au(²S) + H₂(X¹∑⁺_g) → AuH(X¹∑⁺_g) + H(²S) reaction on a global potential energy surface. The reaction probabilities at a series of J values, integral cross sections (ICSs) and differential cross sections of the title reaction are calculated by the TDWP method. For reaction probabilities, there are a mass of sharp oscillations at low collision energy, which can be attributed to resonances supported by the potential well. Due to the endothermicity of the title reaction, the total ICS shows a threshold about 1.53 eV. In order to further investigate the reactive mechanism, the lifetime of complex is calculated by QCT method. At the low collision energy, most intermediate complexes are long lived, which implies that the reaction is governed by indirect reactive mechanism. With the collision energy increasing, the direct reactive mechanism occupies the dominant position. Due to the change of the reactive mechanism, the angular distribution shifts toward the forward direction with collision energy increasing. The isotopic variant, Au + D₂→AuD + D reaction, is also calculated by TDWP method. The calculated reaction probabilities and ICSs show that the isotope effect reduces the reactivity.