Quasi-classical trajectory investigation on the stereodynamics of Li ＋ DF（v=1-6,j=0）→LiF＋D reaction

In this paper, the stereodynamics of Li ＋ DF → Li F ＋ D reaction is investigated by the quasi-classical trajectory（QCT）method on the ^2A＇ potential energy surface（PES） at a relatively low collision energy of 8.76 kcal/mol. The scalar properties of the title reaction such as reaction probability and cross section are studied with vibrational quantum number of v = 1–6. The product angular distributions P（θr） and P（φr） are presented in the same vibrational level range. Moreover, two polarization-dependent generalized differential cross sections（PDDCSs）, i.e., the PDDCS00 and PDDCS22＋are calculated as well. These stereodynamical results demonstrate sensitive behaviors to the vibrational quantum numbers.     

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.     

Quasi-classical trajectory study of the isotope effect on the stereodynamics in the reaction H（2S） + CH（X2Π;u= 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 11A 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.     

Influence of reagent vibration on the stereodynamics of the Li ＋ HF → LiF ＋ H reaction

We investigate the influence of reagent vibration on the stereodynamics of the title reaction by the quasi-classical trajectory on the Aguado-Paniagua2-potential energy surface developed by Aguado et al.(J.Chem.Phys.1997 106 1013).The cross sections and reaction probability as functions of the reagent vibration are calculated in the centre-ofmass frame.The product angular distributions of p(θr),p(φr),and p(θr,φr),which reflect the vector correlation,are also presented and discussed.The results indicate that the vector properties are sensitively affected by the vibrational excitation.     

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.     

Effects of a reagent＇s rotational and vibrational excitations on reactionO（3p）＋ H2（v=0, 3, j = 0, 3, 5, 7, 9, 12, 15） → OH ＋ H

To investigate the effect of a reagent’s rotational and vibrational excitations on the stereo-dynamics of the reaction product, the title reaction is theoretically simulated using the quasi-classical trajectory (QCT) method on the 3 A and 3 A potential energy surfaces (PESs). The reaction cross section is considered as the only scalar property in this work at four different collision energies. Furthermore the vector properties including two polarization-dependent differential cross sections (PDDCSs), the angular distributions of product’ rotational momentum are discussed at one fixed collision energy. Effects of reagents’ rotational excitation on the reaction do exist regularly.     

Quasi-classical trajectory study of the isotope effect on the stereodynamics in the reaction H（2S） ＋ CH（X2∏; v=o,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 11At 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.     

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

Theoretical prediction of the optimal conditions for observing the stereodynamical vector properties of the C(~3P) + OH(X~2Π) → CO(X~1 Σ~+) + H(~2S) reaction

The best optimal initial reactant state and collision energy for observing the stereodynamical vector properties of the title reaction in the ground electronic state X2A potential energy surface(PES) [Zanchet et al. 2006 J. Phys. Chem. A 110 12017] are theoretically predicted using the quasi-classical trajectory(QCT) method for the first time. The calculated results reveal that the smallest value of the rotational quantum number j, larger vibrational quantum number v, and the lower strength of collision energy should be selected for offering the most obvious picture about the stereodynamical vector properties. Polarization-dependent differential cross sections and the angular momentum alignment distribution, P(θr) and P(Φr) in the center-of-mass frame, are obtained to gain an insight into the alignment and orientation of the product molecules. The rotational angular momentum vector j of CO is aligned to be perpendicular to reagent relative velocity k. The product polarizations align along the y axis, pointing to the positive direction of the y axis. A new method is developed to investigate massive reactions with various initial states and to further study the vector properties of the fundamental reactions in detail.     

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.     

A typical slow reaction H（2S） ＋ S2（X3∑g） → SH（X2П） ＋S（3P） ona new surface： Quantum dynamics calculations

Quantum dynamics calculations for the title reaction H(2S) + S2(X3-Σg) → SH(X2Π) + S(3P) are performed by using a globally accurate double many-body expansion potential energy surface [J. Phys. Chem. A 115 5274(2011)].The Chebyshev real wave packet propagation method is employed to obtain the dynamical information, such as reaction probability, initial state-specified integral cross section, and thermal rate constant. It is found not only that there is a reaction threshold near 0.7 eV in both reaction probabilities and integral cross section curves, but also that both the probability and cross section increase firstly and then decrease as the collision energy increases. The existence of the resonance structure in both the probability and cross section curves is ascribed to the deep potential well. The calculation of the rate constant reveals that the reaction occurring on the potential energy surface of the ground-state HS2is slow to take place.     

Dynamics of the reaction of O with H2 and its isotopic variants in different rotational excited states

Scalar properties and vector correlations of the reactions of O+H 2 →OH+H, O+HD→OH+D, O+DH→OD+H, and O+D 2 →OD+D at collision energies of 25 and 34.6 kcal/mole have been studied via the quasi-classical-trajectory (QCT) method based on a BMS1 potential energy surface (PES). The generalized polarization-dependent differential cross section and the distributions of the dihedral angle at the collision energy of 34.6 kacl/mole are presented. The calculated results indicate that both the reagent rotational angular momentum and the mass factor have a significant influence on the scalar properties and vector correlations of the title reactions.     

Influence of reagent vibration on the stereodynamics of the Li + HF → LiF + H reaction

We investigate the influence of reagent vibration on the stereodynamics of the title reaction by the quasi-classical trajectory on the Aguado-Paniagua2-potential energy surface developed by Aguado et al. (J. Chem. Phys. 1997 106 1013). The cross sections and reaction probability as functions of the reagent vibration are calculated in the centre-of-mass frame. The product angular distributions of p(θ_r), p(φ_r), and p(θ_r, φ_r), which reflect the vector correlation, are also presented and discussed. The results indicate that the vector properties are sensitively affected by the vibrational excitation.     

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.     

Theoretical prediction of the optimal conditions for observing the stereodynamical vector properties of the C（3P）＋OH （X2∏）→CO（X1S＋）＋H（2S） reaction

The best optimal initial reactant state and collision energy for observing the stereodynamical vector properties of the title reaction in the ground electronic state X2A’ potential energy surface （PES）[Zanchet et al. 2006 J. Phys. Chem. A 110 12017] are theoretically predicted using the quasi-classical trajectory （QCT） method for the first time. The calculated results reveal that the smallest value of the rotational quantum number j, larger vibrational quantum number v, and the lower strength of collision energy should be selected for offering the most obvious picture about the stereodynamical vector properties. Polarization-dependent differential cross sections and the angular momentum alignment distribution, P（θr） and P（Φr） in the center-of-mass frame, are obtained to gain an insight into the alignment and orientation of the product molecules. The rotational angular momentum vector j’ of CO is aligned to be perpendicular to reagent relative velocity k. The product polarizations align along the y axis, pointing to the positive direction of the y axis. A new method is developed to investigate massive reactions with various initial states and to further study the vector properties of the fundamental reactions in detail.     

Stereodynamics study of the exchange reaction O(~3P)+CH_4 →H+OCH_3

A new London-Eyring-Polanyi-Sato potential energy surface is employed in this work to study the stereo properties of the O(3P)+CH4 →H+CH3O reaction in its rovibrationally ground state using the quasiclassical trajectory method(QCT).Our calculations are performed at a range of collision energies,Ec=1.5eV～3.5eV,and the excitation function obtained by the QCT method accords well with the experimental data.The product rotational polarization is calculated,and the product shows a strong rotational polarization in the centre-of-mass coordinate system.The orientation of the product rotational angular momenta is sensitive to the increase in collision energy,and the alignment of the product rotational angular momenta shows some of the properties of the heavy heavy-light mass combination reactions.In the isotopic substituted reaction study,when the H atoms in methane are replaced by D atoms,the rotational polarization is obviously reduced.The polarization-dependent differential cross section is also studied by this QCT calculation to provide detailed information about the rotational alignment and orientation of the product.     

The reagent vibrational excitation effect on the stereodynamics of the reaction O(~1D)+HBr→OH+Br

Calculations on the dynamics of the reaction O( 1 D) + HBr → OH + Br are performed on the ab initio potential energy surfaces (PESs) of the ground state given by Peterson Peterson K A J. Chem. Phys. 113 4598 (2000)using the quasiclassical trajectory (QCT) method. The product distribution of the dihedral angle, P (φ r ), and that of the angle between and , P (θ r ), are presented in three dimensions. Moreover, we also investigate the reagent vibrational excitation effects on the two polarization-dependent generalized differential cross sections (PDDCS), PDDCS 00 and PDDCS 20 , in the center-of-mass frame. The results indicate that the vector properties are sensitive to the reagent vibrational quantum number.     

Theoretical Study of Isotopic Effect of Oxygen Atom on the Stereodynamics for the O（^3P） ＋ D2 → OD ＋ D Reaction

Quasi-classical trajectory theory is used to study the isotope effect of oxygen atoms on the vector correlations in the O（^3P） ＋ D reaction at a collision energy of 25kcal/mol using accurate potential energy surface of the 3A＇ triplet state. The distributions of p（θr） and the distribution of dihedral angel p（φr） as well as p（θr,φr） are calculated. Moreover, four polarization-dependent generalized differential cross sections （PDDCSs） of product are presented in the center-of-mass frame. The results indicate that the polarization of the product presents different characters for the isotope effect of oxygen atoms. Isotopic substitute can cause obviously different effects on the four PDDCSs.     

Effect of isotope on state-to-state dynamics for reactive collision reactions ■ and ■ in ground state 1~2 A’’ and first excited 1~2 A’ potential energy surfaces

We carry out quantum scattering dynamics and quasi-classical trajectory(QCT) calculations for the O+H_2~+ reactive collision in the ground(1~2 A").nd first excited(1~2 A') potential energy surface.We calculate the reaction probabilities of O+H_2~+(v=0,j=0)→OH~++H and O+H_2~+(v=0,j=0)→OH+H~+reaction for total angular momentum J=0.The results calculated by QCT are consistent with those from quantum mechanical wave packet.Using the QCT method,we generate in the center-of-mass frame the product state-resolved integral cross-sections(ICSs);two commonly used generalized polarization-dependent differential cross-sections(PDDCS s),(2π/σ)(dσ_(00)/dω_t),(2π/σ)(dσ_(20)/dω_t));and three angular distributions of the product rotational vectors,p(θ_r),P(φ_r),and p(θ_r,φ_r).We discuss the influence on the scalar and vector properties of the potential energy surface,the collision energy,and the isotope mass.Since there are deep potential wells in these two potential energy surfaces,their kinetic characteristics are similar to each other and the isotopic effect is not obvious.However,the well depths and configurations of the two potential energy surfaces are different,so the effects of isotopic substitution on the integral cross-section and the rotational polarization of product are different.