O + NH反应在3A’’和1A’’势能面上的量子含时波包动力学研究

对于O + NH反应,在3A’’和1A’’ 势能面（Guadagnini R, Schatz G C, Walch S P. Global potential energy surface for the lowest 1A’, 3A’’, and 1A’’ states of HNO [J]. J. Chem. Phys., 1995,10:774）上,我们运用coupled state or centrifugal sudden (CS)近似和close coupling or Coriolis coupled (CC)方法进行了量子动力学计算。通过比较两种方法得到的总的反应几率,我们发现对于两个势能面上的标题反应,CS近似是失效的。我们还讨论了用CS和CC方法得到的速率常数,并进行了结果比较。     

采用量子波包方法和准经典轨线方法研究S(~3P)+HD反应

利用量子波包方法和准经典轨线方法在一个新的3A′′势能面上研究了S(3P)+HD→SD+H和SH+D反应的动力学性质.计算得到了两个反应通道在碰撞能为0.8—2.2 eV范围内的J=0反应几率、积分反应截面、内同位素因子和产物转动取向因子.这些结果揭示了S(3P)+HD反应非常明显的内同位素效应.通过对势能面和反应轨线的分析,我们发现了SD+H反应通道一个新的反应机制.S(3P)+HD反应的内同位素效应可以利用新发现的反应机制和反应的质量组合来进行解释.     

O+NH反应在3A″和1A″势能面上的量子含时波包动力学研究

对于O+NH反应,在3A″和1A″势能面(Guadagnini R,Schatz G C,Walch S P.Global potential energy surface for the lowest 1A′,3A″,and 1A″states of HNO [J].J.Chem.Phys.,1995,10:774)上,我们运用coupled state or centrifugal sudden (CS)近似和close coupling or Coriolis coupled (CC)方法进行了量子动力学计算.通过比较两种方法得到的总的反应几率,我们发现对于两个势能面上的标题反应,CS近似是失效的.我们还讨论了用CS和CC方法得到的速率常数,并进行了结果比较.     

Cs(6D5/2)+H2→CsH[X1∑+(v,J)]+H转动分辨总截面的测量

利用泵浦-检测方法,在样品池条件下,研究了Cs(6D5/2)与H2反应碰撞传能过程.利用激光感应荧光(LIF)光谱技术,确定了CsH[X1∑+(v,J)]振转能级上的布居分布,转动态分布与热统计分布基本一致.Cs激发态原子密度由激光能量吸收得到.记录A1∑+(v',J+1)→X1∑+(v,J)的时间分辨荧光,从荧光强度的对数值给出的直线斜率确定(v',J+1)→(v,J)的自然辐射率,结合(v,J)→(v',J+1)吸收系数的测量,得到反应生成物CsH[x1∑+(v,J)]态的分子密度.由速率方程分析,给出反应截面σ(v,J),对J求和,得到σ(v)[10-16 cm2单位]分别为(0.64士0.19)(v=0)和(0.58士0.17)(v=1).     

Cs(6D5/2)+H2→CsH[ X1Σ+(v,J)]+H转动分辨总截面的测量

利用泵浦-检测方法,在样品池条件下,研究了Cs(6D5/2)与H2反应碰撞传能过程。利用激光感应荧光（LIF）光谱技术,确定了CsH[X1Σ+(v,J)]振转能级上的布居分布,转动态分布与热统计分布基本一致.Cs激发态原子密度由激光能量吸收得到.记录A1Σ+(v',J+1)→X1Σ+(v,J)的时间分辨荧光,从荧光强度的对数值给出的直线斜率确定(v',J+1)→(v,J)的自然辐射率,结合(v,J)→(v',J+1)吸收系数的测量,得到反应生成物CsH[X1Σ+(v,J)]态的分子密度.由速率方程分析,给出反应截面(v,J),对J求和,得到(v)[10-16cm2单位]分别为（0.64±0.19）（v=0）和（0.58±0.17）（v=1）.     

相对论多组态相互作用方法计算Mg~+离子同位素位移

采用相对论多组态相互作用方法研究了Mg~+离子3s~2S_(1/2)—3s~2P_(1/2)和3s~2S_(1/2)—3s~2P_(3/2)两条跃迁谱线的特殊质量位移系数和场位移因子,并计算了中子数8≤N≤20的Mg~+离子的同位素位移.计算结果与其他理论的计算值符合得比较好,与最新的实验测量结果比较,相对误差在0.13%到0.28%范围,是目前最接近Mg~+离子同位素位移实验测量的理论计算结果.该计算结果可为Mg~+离子同位素位移实验和理论研究提供参考,能够用于Mg~+离子的短寿命同位素的光谱测量实验以及利用Mg~+离子开展幻中子数N=8和N=20附近的奇异原子核特性研究等.所用的计算方法和电子激发模式也可以推广到其他核外电子数为11的多电子体系,用于开展相应的原子光谱结构计算和同位素位移的理论研究.     

N+ND反应的量子力学速率常数

在两个耦合的势能面1²A′和2²A′上对N+ND反应进行了非绝热量子动力学研究. 计算了N+ND→N₂+D反应和N′+ND→N′D+N反应在5 meV～1.0 eV碰撞能的反应几率和积分截面. 结果发现N+ND→N₂+D反应是N+ND反应的主要反应通道.另外,计算了N+ND→N₂+D 反应的速率常数.     

AuCl3催化炔酮分子内环化反应的机理

采用密度泛函理论的B3LYP方法研究炔酮发生分子内环化生成菲衍生物的反应机理.结果表明：在无催化剂和AuCl₃催化剂作用下,反应均可通过[2+2]和[6+2]途径生成产物.无催化剂时,两条途径控制步骤的能垒都较高,但是[2+2]途径的能垒比[6+2]途径的低32.01 kJ·mol^-1,故反应主要通过[2+2]途径进行.在AuCl₃催化作用下,反应的优势途径为[2+2]途径,其能垒为137.05 kJ·mol^-1.比较活化能发现,AuCl₃催化剂明显地降低了反应能垒,使反应能够顺利进行.     

气相中2Ca＋介入N2O与CO, N2O与H2反应的计算研究

本文采用密度泛函理论DFT-UB3LYP方法, 在6-311+G(2d, p)的基组下, 计算研究了气相中Ca+离子介入N2O (1∑+)和CO (1∑+) 与N2O (1∑+) 和 H2 (1∑+g)反应的微观机理. 报道了二重态势能面上各反应物、中间体和过渡态的构型特征及能量, 并用频率分析和内禀反应坐标（IRC）方法对过渡态进行了验证. 计算结果表明,金属离子参与N2O和CO与N2O和H2的反应都分两步进行, 其中Ca+离子对反应N2O (X1∑+) + CO (1∑+)生成N2 (X1∑+g) + CO2 (1∑+g)比较Fe+, Ir+, Pt+等的金属离子有良好的催化作用, 而对反应N2O (X1∑+) + H2 (1∑+g) → N2 (1∑+g) + H2O (1A1) 催化作用不是很好,N2、CaOH+和H是该反应的主要产物,与实验观测结果相符, 并通过对金属离子亲氧性（OA）的计算, 从热力学方面进一步说明主题反应的可行性.     

气相中2Ca＋介入N2O与CO, N2O与H2反应的计算研究

本文采用密度泛函理论DFT-UB3LYP方法, 在6-311+G(2d, p)的基组下, 计算研究了气相中Ca+离子介入N2O (1∑+)和CO (1∑+) 与N2O (1∑+) 和 H2 (1∑+g)反应的微观机理. 报道了二重态势能面上各反应物、中间体和过渡态的构型特征及能量, 并用频率分析和内禀反应坐标（IRC）方法对过渡态进行了验证. 计算结果表明,金属离子参与N2O和CO与N2O和H2的反应都分两步进行, 其中Ca+离子对反应N2O (X1∑+) + CO (1∑+)生成N2 (X1∑+g) + CO2 (1∑+g)比较Fe+, Ir+, Pt+等的金属离子有良好的催化作用, 而对反应N2O (X1∑+) + H2 (1∑+g) → N2 (1∑+g) + H2O (1A1) 催化作用不是很好,N2、CaOH+和H是该反应的主要产物,与实验观测结果相符, 并通过对金属离子亲氧性（OA）的计算, 从热力学方面进一步说明主题反应的可行性.     

二溴卡宾与CH3OH多通道反应的理论计算

对单、三重态CBr2与CH3OH的多通道反应进行了理论计算. 在b3lyp/6-311g*水平上优化了势能面上构型并进行了构型确认. 计算了各物种经零点能校正的CCSD(T)能量. 结果表明, 单重态CBr2既可以与O－H键插入生成1,1-二溴代甲醚P1 (CH3OCHBr2), 也可以与C－H键插入生成2, 2-二溴代乙醇P3(Br2HCCH2OH). 同时, 存在单重态CBr2与CH3OH抽提双氢生成Pabs(CH2O+CH2Br2)的反应. 采用统计热力学及Wigner校正的Erying理论研究了主反应通道的热力学及动力学性质, 并对反应通道上的关键点进行了NBO和AIM分析.     

二溴卡宾与CH3OH多通道反应的理论计算

对单、三重态CBr2与CH3OH的多通道反应进行了理论计算. 在b3lyp/6-311g*水平上优化了势能面上构型并进行了构型确认. 计算了各物种经零点能校正的CCSD(T)能量. 结果表明, 单重态CBr2既可以与O－H键插入生成1,1-二溴代甲醚P1 (CH3OCHBr2), 也可以与C－H键插入生成2, 2-二溴代乙醇P3(Br2HCCH2OH). 同时, 存在单重态CBr2与CH3OH抽提双氢生成Pabs(CH2O+CH2Br2)的反应. 采用统计热力学及Wigner校正的Erying理论研究了主反应通道的热力学及动力学性质, 并对反应通道上的关键点进行了NBO和AIM分析.     

Hamilton–Jacobi equation,reaction action surface and the emergence of the force concept in chemical reaction dynamics

A combination of Hamilton–Jacobi equation and fast marching algorithm can be used to study reaction dynamics by converting the potential energy surface to a reaction action surface. The reaction action surface has been found to be an important tool in theoretical chemistry, allowing us to provide a different force-based perspective of chemical reactions. Several properties such as reaction force, reaction force surface, reaction path force and reaction path force constant have been defined and calculated by using the reaction action surface. This paper investigates these newly defined properties in order to understand the role they play in chemical reaction with reference to a model 4-well potential energy surface.     

Stereo-dynamics of the exchange reaction Ha+LiHb→LiHa+Hb and its isotopic variants

Quasi-classical trajectory (QCT) method is used to calculate the stereo-dynamics of the exchange reaction H_a+LiH_b→LiH_a+H_b and its isotopic variants based on an accurate potential energy surface reported by Prudente et al. [Prudente F V, Marques J M C and Maniero A M 2009 Chem. Phys. Lett. 474 18]. The reactive probability of the title reaction is computed. The vector correlations and four polarization-dependent generalized differential cross sections (PDDCSs) at different collision energies are presented. The influences of the collision energy and the reagent rotation on the product polarization are studied in the present work. The results indicate that the product rotational angular momentum j' is not only aligned, but also oriented along the direction perpendicular to the scattering plane. The product polarization distributions of the title reaction and its isotopic variants exhibit distinct differences which may arise from different mass combinations.     

Stereo-dynamics of the exchange reaction HaWLiHb → LiHa＋Hb and its isotopic variants

<正>The quasi-classical trajectory(QCT) method is used to calculate the stereo-dynamics of the exchange reaction H_a+LiH_b→LiH_a+H_b and its isotopic variants based on an accurate potential energy surface reported by Prudente et al.[Prudente F V,Marques J M C and Maniero A M 2009 Chem.Phys.Lett.474 18].The reactive probability of the title reaction is computed.The vector correlations and four polarization-dependent generalized differential cross sections(PDDCSs) at different collision energies are presented.The influences of the collision energy and the reagent rotation on the product polarization are studied in the present work.The results indicate that the product rotational angular momentum j’ is not only aligned,but also oriented along the direction perpendicular to the scattering plane. The product polarization distributions of the title reaction and its isotopic variants exhibit distinct differences which may arise from different mass combinations.     

Effect of ro-vibrational excitation of NeH~+ on the stereodynamics for the reactions H+NeH~+(v=1-3,j=1,3,5) →H_2~++Ne

The stereodynamics of the abstraction reaction H + NeH+(v = 1-3,j = 1,3,5) → H2+ + Ne is studied theoretically with a quasi-classical trajectory method on a new ab initio potential energy surface [ S J,Zhang P Y,Han K L and He G Z 2012 J.Chem.Phys.132 014303].The effects of vibrational and rotational excitation of reagent molecules on the polarization of the product are investigated.The reaction cross sections,the distributions of P(θr),P(φr),and polarizationdependent differential cross sections(PDDCSs) are calculated.The obtained cross sections indicate that the title reaction is a typical barrierless atom(ion)-ion(molecule) reaction.The initial vibrational excitation and rotational excitation of reagent molecules have distinctly different influences on stereodynamics of the title reaction,and the possible reasons for the differences are presented.     

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.     

Quasi-classical trajectory study of the stereodynamics of a Ne+H2+→NeH++H reaction

We have carried out a quasi-classical trajectory calculation for the reaction of Ne + H₂⁺ (v=0, j=1) → NeH⁺ + H on the ground state (1²A') using the LZHH potential energy surface constructed by Lü et al. [Lü S J, Zhang P Y, Han K L and He G Z 2010 J. Chem. Phys. 132 014303]. Differential cross sections at many collision energies indicate that the reaction is dominated by forward-scattering. In addition, the NeH⁺ product shows rotationally hot and vibrationally cold distributions. Stereodynamical results indicate that the products are strongly polarized in the direction perpendicular to the scattering plane and that the products rotate mainly in planes parallel to the scattering plane.     

Electronic transition of C3H− in the vicinity of the lowest photodetachment threshold

The electronic spectrum of the C³H(D)^? anion has been studied near the lowest electron detachment threshold associated with the chain isomer. On the basis of the vibrational and rotational analyses, the photodetachment spectrum is assigned to the ³A″ ← [Xtilde]³A″ electronic transition of the cumulene-like C³H(D)^? anion. The spectrum is characterized by a rich vibrational structure leading to the determination of almost all fundamental frequencies. The vibrational progression involving the CCH(D) in plane bending mode (V¹ ₆) is fitted with a 1-dimensional potential energy curve along the CCH(D) angle. This fit and the rotational analysis enable the geometry of the upper state to be determined. Although this excited state is expected to have dipole bound character, its neutral counterpart has a significantly different structure. The observation of an excited triplet state near the electron detachment threshold and close to the the calculated transition state for the electron catalysed isomerization reaction leads to new insight into this process.