The stereodynamic properties of the F+HO(v,j) → HF+O reaction on~1 A' and ~3A' potential energy surfaces by quasi-classical trajectory calculations:Initial excitation effect(v=1-3, j=0 and v= 0, j=1-3)

The stereodynamic properties of the F ＋ HO （v, j） reaction are explored by quasi-classical trajectory （QCT） calculations performed on the 1At and 3At potential energy surfaces （PESs）. Based on the polarization-dependent differential cross sections （PDDCSs） and the angular distributions of the product angular momentum with the reactant at different values of initial v or j, the results show that the product scattering and product polarization have strong links with initial vibrationalrotational numbers of v and j. The significant manifestation of the normal DCSs is that the forward scattering gradually becomes predominant with the initial vibrational excitation increasing, and the scattering angle of the HF product taking place on the 3At potential energy surface is found to be more sensitive to the initial value of v. The product orientation and alignment are strongly dependent on the initial rovibrational excitation effect. With enhancement in the initial rovibrational excitation effect, there is an overall decrease in the product orientation as well as in the product alignment either perpendicular to the reagent relative velocity vector k or along the direction of the y axis, for which the initial rotational excitation effect is much more noticeable than the vibrational excitation effect. Moreover, the initial rovibrational excitation effect on the product polarization is more pronounced for the 3At potential energy surface than for the 1At potential energy surface.     

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.     

Product polarization and mechanism of Li + HF(v=0, j=0) → LiF(v', j') + H collision reaction

A state-to-state dynamics analysis for the Li+HF(v=0,j=0)→LiF(v',j')+H collision reaction has been performed through quasiclassical trajectory(QCT)calculations.It is found that the differential cross section(DCS)of the LiF products from the title reaction is preferentially backward scattering for v=0,yet forward scattering for v=1 and 2.For v=3,the DCS exhibits forward,backward,and sideways scatterings.The variation of the internuclear distances and angles along the propagation time reveals that more than 99.08%of reaction trajectories undergo the direct reaction mechanism.The values of the polarization parameters a{1}1and a{2}0demonstrate that the product rotational angular moment j' is not only aligned perpendicular to the reagent relative velocity vector,but also oriented along the negative y axis.These product polarization results agree well with the recent quantum mechanical studies.The mechanism of these results was proposed and discussed in detail.     

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.     

C+OH(v=0-3,j=0-6)→CO+H反应的准经典轨线研究

运用准经典轨线方法,在碰撞能为1.0eV时,研究了反应物OH分子的振转激发对C+OH反应的立体动力学性质的影响.详细地讨论了该反应在不同反应物振转态下的矢量性质.计算表明,OH分子的振转激发对C+OH反应的矢量性质非常敏感,这种现象在对该反应的标量性质的研究中是不存在的.     

Theoretical study of stereodynamics for the D'+DS(ν = 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 3A' 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π/σ)(ds₀₀/dω_t), (2π/σ)(ds₂₀/dω_t), (2π/σ)(ds₂₂₊/dω_t), (2π/σ)(ds_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.     

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.     

C+OH(v=0—3,j=0—3)→CO+H反应的准经典轨线研究

运用准经典轨线方法,在碰撞能为1.0eV时,研究了反应物OH分子的振转激发对C+OH反应的立体动力学性质的影响.详细地讨论了该反应在不同反应物振转态下的矢量性质.计算表明,OH分子的振转激发对C+OH反应的矢量性质非常敏感,这种现象在对该反应的标量性质的研究中是不存在的.     

低碰撞能下Li+HF(v=0－3，j=0－40)→LiF+H反应的立体动力学研究

基于2003年势能面,运用准经典轨线法(QCT)研究Li+HF→LiF+H反应立体动力学.探究较低碰撞能(1.15 kcal·mol-1-5.00 kcal·mol-1)下碰撞能、振转激发对极化微分反应截面(PDDCSs)和三矢量相关的P(θr,r)分布函数的影响,将积分散射截面与已有的理论及实验结果比较.结果显示,在较低碰撞能下碰撞能、振转激发对极化微分散射截面和三矢量相关的P(θr,r)分布函数有影响,但振转激发对极化微分反应截面和P(θr,r)分布的影响更大,碰撞能的增加使产物转动角动量后向散射的极化强度增大.在计算的能量范围内积分散射截面与其它的理论及实验结果符合较好.     

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.     

Non-adiabatic quantum dynamical studies of Na(3p) + HD(v = 1, j5) = 0)→NaH/NaD + D/H reaction

Non-adiabatic dynamical calculations are carried out for the Na(3 p)+HD(ν = 1, j = 0)→NaH/NaD+D/H reaction on the diabatic potential energy surfaces of Wang et al.(Sci. Rep. 2018, 8, 17960) by using the time-dependent wave packet method. The state-to-state integral cross sections and differential cross sections of two reaction channels(NaH/NaD+D/H)are calculated for collision energy up to 0.4 eV. The cross section branching ratio indicates that the dominant reaction channel changes from NaD+H to NaH+D when the collision energy is larger than 0.227 eV. The products from two reaction channels both prefer to form in vibrationally cold but rotationally hot states, and they both tend to forward scattering.     

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.     

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.     

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.     

低碰撞能下H+HF(v=3,j=0)反应的量子动力学研究

报道了H原子和振动激发的HF(v=3,j=0)分子在低碰撞能下的量子反应动力学研究.计算结果表明:在低碰撞能下,散射主要以非反应过程为主;振动激发有利于反应的进行;在小于10-4eV碰撞能下,非反应非弹性散射截面和反应截面的比值约为3.在反应截面上发现了Feshbach共振现象,证实是由反应通道上紧邻反应势垒的H…HF(v=3,j=1～3)的范德瓦尔斯聚合体存在的准束缚态所形成.     

碰撞能对S(1D) +H2(v=0, j=0)→SH+H反应立体动力学性质的影响

利用1A′态的势能面[ Ho et al., J. Chem. Phys. 116, 4124 (2002)],采用准经典轨线方法研究了在不同碰撞能条件下,S(1D) +H2(v=0, j=0)→SH+H反应的立体动力学性质. 通过计算得到了描述反应物速度矢量k与产物的转动角动量矢量j′这两个矢量相关的分布函数P(r)、描述反应物速度矢量k、产物速度矢量k′与产物的转动角动量矢量j′这三个矢量相关的二面角分布函数P(r)以及描述反应产物角动量极化的分布函数P(r,r).计算结果表明产物的转动角动量矢量j′在空间具有明显的定向和取向效应,并且产物的转动角动量具有强烈的极化. 另外,计算结果还表明这些立体动力学性质对碰撞能非常敏感.     

碰撞能对S(1D) +H2(v=0, j=0)→SH+H反应立体动力学性质的影响

利用1A′态的势能面[ Ho et al., J. Chem. Phys. 116, 4124 (2002)],采用准经典轨线方法研究了在不同碰撞能条件下,S(1D) +H2(v=0, j=0)→SH+H反应的立体动力学性质. 通过计算得到了描述反应物速度矢量k与产物的转动角动量矢量j′这两个矢量相关的分布函数P(r)、描述反应物速度矢量k、产物速度矢量k′与产物的转动角动量矢量j′这三个矢量相关的二面角分布函数P(r)以及描述反应产物角动量极化的分布函数P(r,r).计算结果表明产物的转动角动量矢量j′在空间具有明显的定向和取向效应,并且产物的转动角动量具有强烈的极化. 另外,计算结果还表明这些立体动力学性质对碰撞能非常敏感.     

A full-dimensional analytical potential energy surface for the F+CH_4→HF+CH_3 reaction

A full-dimensional analytical potential energy surface(APES) for the F + CH4 →HF + CH3 reaction is developed based on 7127 ab initio energy points at the unrestricted coupled-cluster with single,double,and perturbative triple excitations.The correlation-consistent polarized triple-split valence basis set is used.The APES is represented with a many-body expansion containing 239 parameters determined by the least square fitting method.The two-body terms of the APES are fitted by potential energy curves with multi-reference configuration interaction,which can describe the diatomic molecules(CH,H2,HF,and CF) accurately.It is found that the APES can reproduce the geometry and vibrational frequencies of the saddle point better than those available in the literature.The rate constants based on the present APES support the experimental results of Moore et al.[Int.J.Chem.Kin.26,813(1994)].The analytical first-order derivation of energy is also provided,making the present APES convenient and efficient for investigating the title reaction with quasiclassical trajectory calculations.