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.     

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.     

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.     

Theoretical study of stereodynamics for the D'+ DS(v=0,j=0) → D'D + S abstraction reaction

Quasiclassical trajectory (QCT) calculations have been performed for the abstraction reaction, D'+ DS(v = 0, j = 0) → D'D + S on a new LZHH potential energy surface (PES) of the adiabatic 3 A electronic state [Lü et al. 2012 J. Chem. Phys. 136 094308]. The collision energy effect on the integral cross section and product polarization are studied over a wide collision energy range from 0.1 to 2.0 eV. The cross sections calculated by the QCT procedure are in good accordance with previous quantum wave packet results. The three angular distribution functions, P(θr), P(φr), and P(θr,φr), together with the four commonly used polarization-dependent differential cross sections ((2π/σ)(dσ00/dωt), (2π/σ)(dσ20/dωt), (2π/σ)(dσ22+/dωt), (2π/σ)(dσ21/dωt)) are obtained to gain insight into the chemical stereodynamics of the title reaction. Influences of the collision energy on the product polarization are exhibited and discussed.     

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.     

The characteristics of O＋HD （v = 0, j = 0） reaction dynamics

The analytical potential energy function of HDO is constructed at first using the many-body expansion method.The reaction dynamics of O+HD(v = 0,j = 0) in five product channels are all studied by quasi-classical trajectory(QCT) method.The results show that the long-lived complex compound HDO is the dominant product at low collision energy.With increasing collision energy,O+HD → OH+D and O+HD → OD+H exchange reactions will occur with remarkable characteristics,such as near threshold energies,different reaction probabilities,and different reaction cross sections,implying the isotopic effect between H and D.With further increasing collision energy(e.g.,up to 502.08 kJ/mol),O+HD → O+H+D will occur and induce the complete dissociation into single O,H,and D atoms.     

Theoretical prediction of energy dependence for D+BrO→DBr+O reaction:The rate constant and product rotational polarization

A quasi-classical trajectory(QCT) calculation is used to investigate the vector and scalar properties of the D + Br O → DBr + O reaction based on an ab initio potential energy surface(X1A state) with collision energy ranging from 0.1 kcal/mol to 6 kcal/mol. The reaction probability, the cross section, and the rate constant are studied. The probability and the cross section show decreasing behaviors as the collision energy increases. The distribution of the rate constant indicates that the reaction favorably occurs in a relatively low-temperature region(T 100 K). Meanwhile, three product angular distributions P(θr), P(φr), and P(θr, φr) are presented, which reflect the positive effect on the rotational angular momentum j' polarization of the DBr product molecule. In addition, two of the polarization-dependent generalized differential cross sections(PDDCSs), PDDCS00 and PDDCS20, are computed as well. Our results demonstrate that both vector and scalar properties have strong energy dependence.     

A theoretical study of the stereodynamics on the abstraction reactions H/D ＋ HS/DS

The quasi-classical trajectory(QCT) method is employed to calculate the stereodynamics of the abstraction reactions H/D+HS/DS based on an accurate potential energy surface [L S J,Zhang P Y,Han K L and He G Z 2012 J.Chem.Phys.136 094308].The reaction cross sections of the title reaction are computed,and the vector correlations for different collision energies and different initial vibrational states are presented.The influences of the collision energy and reagent vibration on the product polarization are studied,and the product polarizations of the title reactions are found to be distinctly different,which arises from the different mass factors,collision energies,and reagent vibrational states.     

Quantum stereodynamics study for the reaction F + HD

This paper studies the quantum stereodynamics of the F + HD(υ = 0, j = 0) → HD + F/HF + D reaction at the collision energies of 0.52 and 0.87~kcal/mol. The quantum scattering calculations, based on Stark-Werner potential energy surfaces, show that the differential cross sections for the HF(υ' = 2) + D and DF(υ' = 3) + H channels are consistent with the recent theoretical results. Furthermore, the product rotational angular momentum orientation and alignment have been determined for some selected rovibrational states of the HF + D and DF + H channels.     

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).     

Stereodynamics of the O（3P） with H2（D2） （ν=0,j=0） reaction

Quasi-classical trajectory theory is used to study the reaction of O(3 P) with H 2 (D 2) based on the ground 3 A″ potential energy surface (PES).The reaction cross section of the reaction O+H 2 →OH+H is in excellent agreement with the previous result.Vector correlations,product rotational alignment parameters P 2 (j · k) and several polarizeddependent differential cross sections are further calculated for the reaction.The product polarization distribution exhibits different characteristics that can be ascribed to different motion paths on the PES,arising from various collision energies or mass factors.     

Effects of the vibrational and rotational excitation of reagent on the stereodynamics of the reaction S（3P） + H2→ SH + H

Quasiclassical trajectory （QCT） calculations are first carried out to study the stereodynamics of the S （3p） ＋ H2 → SH ＋ H reaction based on the ab initio 13Atr potential energy surface （PES） （Lii etal. 2012 J. Chem. Phys. 136 094308）. The QCT-calculated reaction probabilities and cross sections for the S ＋ H2 （v = 0, j = 0） reaction are in good agreement with the previous quantum mechanics （QM） results. The vector properties including the alignment, orientation, and polarization- dependent differential cross sections （PDDCSs） of the product SH are presented at a collision energy of 1.8 eV. The effects of the vibrational and rotational excitations of reagent on the stereodynamics are also investigated and discussed in the present work. The calculated QCT results indicate that the vibrational and rotational excitations of reagent play an important role in determining the stereodynamic properties of the title reaction.     

Reaction mechanism of D + ND→N + D_2 and its state-to-state quantum dynamics

The quantum state-to-state calculations of the D + ND→N + D_2 reaction are performed on a potential energy surface of 4 A' state. The state-resolved integral and differential cross sections and product state distributions are calculated and discussed. It is found that the rotational distribution, rather than the vibrational distribution, of the product has an obvious inversion. Due to the fact that it is a small-impact-parameter collision, its product D_2 is mainly dominated by rebound mechanism, which can lead to backward scattering at low collision energy. As the collision energy increases, the forward scattering and sideward scattering begin to appear. In addition, the backward collision is also found to happen at high collision energy, through which we can know that both the rebound mechanism and stripping mechanism exist at high collision energy.     

Theoretical study of stereodynamics for the reaction O（3P）+D2（v= 0,j=0） → OD+D and isotope effect

Quasi-classical trajectory (QCT) calculations have been performed to study the product polarization behaviours in the reaction O(3P) + D2 (v = 0, j = 0) → OD + D. By running trajectories on the 3A and 3A potential energy surfaces (PESs), vector correlations such as the distributions of the polarization-dependent differential cross sections (PDDCSs), the angular distributions of P (θr) and P (φr) are presented. Isotope effect is discussed in this work by a comprehensive comparison with the reaction O(3P) + H2 (v = 0, j = 0) → H + H. Common characteristics as well as differences are discussed in product alignment and orientation for the two reactions. The isotope mass effect differs on the two potential energy surfaces: the isotope mass effect has stronger influence on P (θr) and PDDCSs of the 3A PES while the opposite on P (φr) of the 3A potential energy surface.     

Stereodynamics of the O（^3p） with H： （D2） （v = O,j =O） reaction

Quasi-classical trajectory theory is used to study the reaction of O（3p） with H2 （D2） based on the ground 3A″ potential energy surface （PES）. The reaction cross section of the reaction O＋H2→＋OH＋H is in excellent agreement with the previous result. Vector correlations, product rotational alignment parameters （P2（j′. k）） and several polarizeddependent differential cross sections are further calculated for the reaction. The product polarization distribution exhibits different characteristics that can be ascribed to different motion paths on the PES, arising from various collision energies or mass factors.     

Isotope effect on the stereodynamics for the collision reaction H＋LiF（v = 0, j = 0） → HF＋Li

Stereodynamics for the reaction H+LiF(v=0, j=0) → HF+Li and its isotopic variants on the ground-state (1 2 A′) potential energy surface (PES) are studied by employing the quasi-classical trajectory (QCT) method. At a collision energy of 1.0 eV, product rotational angular momentum distributions P (θr), P (φr), and P (θr ,φr), are calculated in the center-of-mass (CM) frame. The results demonstrate that the product rotational angular momentum j′ is not only aligned along the direction perpendicular to the reagent relative velocity vector k, but also oriented along the negative y axis. The four generalized polarization-dependent differential cross sections (PDDCSs) are also computed. The PDDCS 00 distribution shows a preferential forward scattering for the product angular distribution in each of the three isotopic reactions, which indicates that the title collision reaction is a direct reaction mechanism. The isotope effect on the stereodynamics is revealed and discussed in detail.     

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.     

The effect of collision energy on the stereodynamics of the reaction H(~2S)+NH(X~3∑~-,v = 0,j = 0)→ N(~4S)+H_2

The quasi-classical trajectory(QCT) is calculated to study the stereodynamics properties of the title reaction H(2S)+NH(X3∑-) →N(4S)+H2 on the ground state 4A' potential energy surface(PES) constructed by Zhai and Han [2011 J.Chem.Phys.135 104314].The calculated QCT reaction probabilities and cross sections are in good agreement with the previous theoretical results.The effects of the collision energy on the k-k' distribution and the product polarization of H2 are studied in detail.It is found that the scattering direction of the product is strongly dependent on the collision energy.With the increase in the collision energy,the scattering directions of the products change from backward scattering to forward scattering.The distribution of P(θr) is strongly dependent on the collision energy below the lower collision energy(about 11.53 kcal/mol).In addition,the P(φr) distribution dramatically changes as the collision energy increases.The calculated QCT results indicate that the collision energy plays an important role in determining the stereodynamics of the title reaction.     

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.     

Quasi-classical trajectory study of collision energy effect on the stereodynamics of H + Br O → O + HBr reaction

Quasi-classical trajectory(QCT) studies on the stereodynamics of H + Br O → O + HBr reaction have been performed on the X1A′state of ab initio potential energy surface by Peterson [Peterson K A 2000 J. Chem. Phys. 113 4598] in a collision energy range from 0 kcal/mol to 6 kcal/mol. Two of the polarization-dependent generalized differential cross sections(PDDCSs),(2π /σ)( dσ00/ dωt)(PDDCS00) and(2π /σ)( dσ20/ dωt)(PDDCS20) are considered. The rotational polarizations of these products show sensitive behaviors to the calculated collision energy range. Furthermore, in order to gain more knowledge about vector correlations, the product angular distribution, P(θr), and the dihedral angle, P(φr),are calculated, and the results indicate that both the rotational alignment and orientation of the product are enhanced as collision energy increases.