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.     

Stereodynamics study of reactions N(~2D)+HD→NH+D and ND+H

The product polarizations of the title reactions are investigated by employing the quasi-classical trajectory (QCT) method. The four generalized polarization-dependent differential cross-sections (PDDCSs) $({2\pi } / \sigma )(\d\sigma _{00} / \d\omega _t )$, $({2\pi } / \sigma )(\d\sigma _{20} / \d\omega _t )$, $({2\pi } / \sigma )(\d\sigma _{22 + } / \d\omega _t )$, and $({2\pi } / \sigma )(\d\sigma _{21 - } / \d\omega _t )$ are calculated in the centre-of-mass frame. The distribution of the angle between ${{\bm k}}$ and ${{\bm j^\prime }}$, $P(\theta _r )$, the distribution of the dihedral angle denoting ${{\bm k}}${--}$\bm k^\prime $--$\bm j^\prime $ correlation, $P(\phi _r )$, as well as the angular distribution of product rotational vectors in the form of polar plots $P(\theta _r ,\phi _r )$ are calculated. The isotope effect is also revealed and primarily attributed to the difference in mass factor between the two title reactions.     

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.     

Time-dependent quantum dynamics study for reaction of D+CH4→CH3+HD

The semirigid vibrating rotor target (SVRT) model has been applied to the study of the reaction of D+CH_4→CH_3+HD using a time-dependent wave packet method. The energy dependence of the calculated reaction probability shows oscillatory structures similar to those observed in the abstraction reaction of H+H_2, H+CH_4 etc. We have also studied the influence of rotational and vibrational excitation of the reacting molecule (CH_4) on reaction probability. The excitation of the H－CH_3 stretching vibration gives significant enhancement of reaction probability, which rises significantly with the enhancement of rotational quantum number j. Finally, we have compared the cross section and the rate constant of the D+CH_4 system with that of the H+CH_4 system.     

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.     

Isotopic effects of the N(2D) + H2 → NH + H reaction: a quantum time-dependent wave packet investigation

Based on the potential energy surface reported by Li and co-workers (J. Comput. Chem. 34 1686–1696 (2013)), the dynamics calculations of N(²D)?+?H₂(v ₀?=?0, j ₀?=?0) reaction and its isotopic variants HD and D₂ are studied using time-dependent wave packet method in the collision energy range of 0.01–1.0?eV. Dynamics properties such as reaction probability, differential cross section, and integral cross section are studied at state-to-state level of theory. Present values are compared with available theoretical and experimental results. The results indicate that the integral cross sections of N(²D)?+?D₂ reaction are in general good agreement with the experimental data at collision energy below 0.15?eV. The rotational state-resolved integral cross sections of N(²D)?+?H₂/HD/D₂ reactions are compared with experimental values for the first time, with the obtained values being in good agreement with the experimental data.     

State-to-state dynamics of reactions H+DH'(v=0,j=0)→HH'(v',j')+D/HD(v',j')+H'with time-dependent quantum wave packet method

State-to-state time-dependent quantum dynamics calculations have been carried out to study H+DH'→HH'+D/HD+H' reactions on BKMP2 surface.The total integral cross sections of both reactions are in good agreement with earlier theoretical and experimental results,moreover the rotational state-resolved reaction cross sections of H+DH'→HH‘+D at collision energy Ec=0.5 eV are closer to the experimental values than the ones calculated by Chao et al [J.Chem.Phys.117 8341(2002)],which proves the higher precision of the quantum calculation in this work.In addition,the state-to-state dynamics of H+DH'→ HD'+H reaction channel have been discussed in detail,and the differences of the micro-mechanism of the two reaction channels have been revealed and analyzed clearly.     

Exact quantum dynamics study of the H(2S)+SiH+(X1Σ+) reaction on a new potential energy surface of SiH2+(X2A1)

Based on a new global potential energy surface of SiH₂⁺(X²A₁), the exact quantum dynamical calculation for the H(²S)+SiH⁺(X¹Σ⁺)→H₂+Si⁺ reaction has been carried out by using the Chebyshev wave packet method. The initial state specified (ν_i=0, j_i=0) probabilities, integral cross sections (ICS) and thermal rate constants of the title reaction are calculated. All partial wave contributions up to J=90 are calculated in exact quantum calculation including the full Coriolis coupling (CC) effect. The dynamical behaviors of probabilities, ICSs and rate constants are found to be in accord with an exothermic reaction without potential barrier. By comparing the probabilities of CC with the corresponding centrifugal sudden (CS) approximation ones, it can be concluded that neglecting CC effect will decrease the collision time, increase the amplitude of oscillation and lead to overestimation or underestimation of the reaction probability. For ICSs and rate constants, it is found that the deviation of CC and CS ICSs is small in the most of collision energy range except for the range of 0 eV-0.05 eV, while the deviation of both rate constants is considerable in the temperature range of 16 K-1000 K.     

反应N(2D)+H2→NH+H和N(2D)+D2→ND+D的立体动力学理论研究.

基于Ho等人的精确势能面(J. Chem. Phys. 119, 3063(2003))研究,运用准经典轨线方法计算了21.3 kJ/mol碰撞能下反应N(²D)+H₂→NH+H和N(²D)+D₂→ND+D的产物与反应物之间的矢量相关．发现两个反应的产物角分布都是前向和后向呈现峰值分布,产物的转动角动量矢量j′不仅是取向的,而且是在y轴负方向上定向的．两个反应显示出的同位素效应主要归因于同位素质量的差别．     

Quasiclassical trajectory theoretical study on the chemical stereodynamics of the O(~1D)+H_2→OH+H reaction and its isotopic variants (HD,D_2)

The quasiclassical trajectory (QCT) method is used to study stereodynamic information about the reaction O ( 1 D) + H 2 →OH + H on the DK (Dobbyn and Knowles) (1 1 A' ) ab initio potential energy surface (PES). A wide scale of collision energy (E c ) from 0.05 eV to 0.5 eV is considered in the dynamic calculations. To reveal the rovibrational excitation effect, calculations at a collision energy of 0.52 eV are carried out for the v = 0 ～ 5, j = 0 and v = 0, j = 0 ～ 15 initial states. The two popularly used polarization-dependent differential cross sections (PDDCSs), dσ 00 /dω t (0, 0) and dσ 20 /dω t (2, 0), and two angular distributions, P(θ r ) and P( r ) are calculated to obtain an insight into the alignment and the orientation of the product molecules. From the calculations, we can obtain that the alignment of the OH product is weaker at high collision energy and becomes stronger with the increase of initial vibrational level, and it is almost insensitive to the initially rotational excitation. Influences of the mass values of isotopes (HD, D 2 ) on the stereodynamics are also shown and discussed. Comparisons between available theoretical results and experimental results are made and discussed.     

Influence of rotational excitation and collision energy on the stereo dynamics of the reaction： N（4S）＋H2 （v = 0, j = 0, 2, 5, 10） →NH（X3∑-）＋H

The N+H₂ reaction has attracted a great deal of attention from both the experimental and the theoretical community, and most of the attention has been paid to the first excited state N(²D) atoms in collisions with hydrogen molecules and the scalar properties of the reaction. In this paper, we study the stereo dynamical properties and calculate the reaction cross sections of the N(⁴S) + H₂ (v=0, j=0, 2, 5, 10) → NH(X³Σ^-) + H using the quasi-classical trajectory (QCT) method on an accurate NH₂ potential energy surface (PES) reported by Poveda and Varandas [Poveda L A and Varandas A J C 2005 Phys. Chem. Chem. Phys. 7 2867], in a collision energy range of 25 kcal·mol^-1-140 kcal·mol^-1. Results indicate that the reactant rotational excitation and initial collision energy both have a considerable influence on the distributions of the k-j′ correlation, the k-k′-j′ correlation and k-k′ correlation. The differential cross section is found to be sensitive to collision energy.     

Vector correlations study of the reaction N(~2D) + H_2(X~1Σ_g~+) →NH(a~1?) + H(~2S) with different collision energies and reagent vibration excitations

Vector correlations of the reaction N(2D)+ H2(X1Σ+g) → NH(a1?)+ H(2S) are studied based on a recent DMBESEC PES for the first excited state of NH2[J. Phys. Chem. A 114 9644(2010)] by using a quasi-classical trajectory method.The effects of collision energy and the reagent initial vibrational excitation on cross section and product polarization are investigated for v = 0–5 and j = 0 states in a wide collision energy range(10–50 kcal/mol). The integral cross section could be increased by H2 vibration excitation remarkably based on the DMBE-SEC PES. The different phenomena of differential cross sections with different collision energies and reagent vibration excitations are explained. Particularly,the NH molecules are scattered mainly in the backward hemisphere at low vibration quantum number and evolve from backward to forward direction with increasing vibration quantum number, which could be explained by the fact that the vibrational excitation enlarges the H–H distance in the entrance channel, thus enhancing the probability of collision between N atom and H atom. A further study on product polarization demonstrates that the collision energy and vibrational excitation of the reagent remarkably influence the distributions of P(θr), P(φr), and P(θr, φr).