Semiclassical glory analyses in the time domain for the H + D2(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction

We make the first application of semiclassical (SC) techniques to the plane-wavepacket formulation of time-domain (T-domain) scattering. The angular scattering of the state-to-state reaction, H + D(2)(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D, is analysed, where v and j are vibrational and rotational quantum numbers, respectively. It is proved that the forward-angle scattering in the T-domain, which arises from a delayed mechanism, is an example of a glory. The SC techniques used in the T-domain are: An integral transitional approximation, a semiclassical transitional approximation, a uniform semiclassical approximation (USA), a primitive semiclassical approximation and a classical semiclassical approximation. Nearside-farside (NF) scattering theory is also employed, both partial wave and SC, since a NF analysis provides valuable insights into oscillatory structures present in the full scattering pattern. In addition, we incorporate techniques into the SC theory called "one linear fit" and "two linear fits", which allow the derivative of the quantum deflection function, Θ?(')(J), to be estimated when Θ?J exhibits undulations as a function of J, the total angular momentum variable. The input to our SC analyses is numerical scattering (S) matrix data, calculated from accurate quantum collisional calculations for the Boothroyd-Keogh-Martin-Peterson potential energy surface No. 2, in the energy domain (E-domain), from which accurate S matrix elements in the T-domain are generated. In the E-domain, we introduce a new technique, called "T-to-E domain SC analysis." It half-Fourier transforms the E-domain accurate quantum scattering amplitude to the T-domain, where we carry out a SC analysis; this is followed by an inverse half-Fourier transform of the T-domain SC scattering amplitude back to the E-domain. We demonstrate that T-to-E USA differential cross sections (DCSs) agree well with exact quantum DCSs at forward angles, for energies where a direct USA analysis in the E-domain fails.     

The polarization dependent differential cross sections of the reactions:H+LiH~+(v=0,j=0)→H_2+Li~+ and H~++ LiH(v=0,j=0)→H_2~++Li

<正>Quasi-classical trajectory(QCT) calculations have been carried out to study the generalized polarization dependent differentialcross sections(PDDCSs) for the reactions H + LiH~+(v = 0,j = 0)→H_2 + Li~+ and H~+ + LiH(v = 0,j = 0)→H_2~+ + Li occurring onthe two lowest-lying electronic states of the LiH_2~+ system,using the ab initio potential energy surfaces(PESs) of Martinazzo et al.[3].Four PDDCSs,i.e.,(2π/σ)(dσ_(00)/dω_t),(2π/σ)(dσ_(20)/dω_t),(2π/σ)(dσ_(22+)/dω_t),(2π/σ)(dσ_(21-)/dω_t) have been discussed in detail.     

Angular scattering using parameterized S matrix elements for the H + D2(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction: an example of Heisenberg's S matrix programme

A neglected topic in the theory of reactive scattering is the use of parameterized scattering (S) matrix elements to calculate differential cross sections (DCSs). We construct four simple parameterizations, whose moduli are smooth step-functions and whose phases are quadratic functions of the total angular momentum quantum number. Application is made to forward glory scattering in the DCS of the H + D(2)(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction at a translational energy of 1.81 eV, where v and j are vibrational and rotational quantum numbers respectively. The parameterized S matrix elements can reproduce the forward scattering for centre-of-mass reactive scattering angles up to 30° and can identify the total angular momenta (equivalently, impact parameters) that contribute to the glory. The theoretical techniques employed to analyze structure in the DCS include: nearside-farside theory, local angular momentum theory--in both cases incorporating resummations of the partial wave series representation of the scattering amplitude--and the uniform semiclassical theory of forward glory scattering. Our approach is an example of Heisenberg's S matrix programme, in which no potential energy surface is used. Our calculations for the DCS using the four parameterized S matrix elements are counterexamples to the following universal statements often found in the chemical physics literature: "every molecular scattering investigation needs detailed information about the interaction potential," and "an accurate potential energy surface is an essential element in carrying out simulations of a chemical reaction". Both these statements are false.     

准经典轨线研究Li+HF(v=0,j=O)→LiF+H反应动力学

基于Aguado等人拟合的APW势能面(PES),运用准经典轨线(QCT)方法,对反应Li+HF(v=0,j=0)→LiF+H的动力学性质进行了计算.主要研究了不同碰撞能条件下的反应截面、转动取向、产物散射角分布和竞争反应模式等.结果表明,该反应存在直接提取型和间接插入型两种反应模式,在低能量下反应以间接插入反应模式为主,能量大于200 meV时则以直接提取反应为主.     

Product rotational angular momentum polarization in the H+FCl(v=0-5, j=0, 3, 6, 9)→HF+Cl reaction

The product alignment and orientation of the title reaction on the ground potential energy surface of 1 (2)A' have been studied using the quasi-classical trajectory method. The calculations were carried out for case (a) at collision energies of 0.5-20 kcal mol(-1) with the initially rovibrational state of the reagent FCl molecule being at the v = 0 and j = 0 level to especially reveal in detail the dependence of the product integral cross section on collision energy. Further calculations at the collision energy of 15 kcal mol(-1) for case (b) at v = 0-5, and j = 0, and (c) at v = 0, and j = 3, 6, 9 initial states were carried out to reveal the effect of initially vibrational and rotational excitations on stereodynamics, respectively. Possessing final relative velocity k' (defined as a vector in the xz-plane), product alignment perpendicular to the reagent relative velocity vector k (defined as z- or parallel to the z-axis), for case (a) is found to be weaker at all collision energies, for case (b) is found to be vibrationally enhanced by the reactant molecule FCl, but for case (c), rather insensitive to initially rotational excitation. The rotational vector of product molecular orientation pointing to either negative or positive direction of the y-axis in the center of mass frame, e.g. origin of the coordinate system, is enhanced by collision energies regarding to 0.5-20 kcal mol(-1), while it becomes weaker at higher vibrational (v = 0-5) or rotational (j = 0, 3, 6, 9) excitation levels. Effects of collision energies and of rotational excitation at these collision energies, with 15 kcal mol(-1) as an example on the calculated PDDCSs are also shown and discussed. Detailed plots P(φ(r)) in the range of 0 ≤φ(r)≤ 360(o), and P(θ(r), φ(r)) in the ranges of 0 ≤θ(r)≤ 180° and 0 ≤φ(r)≤ 360° at collision energies 0.5-20 kcal mol(-1) have been presented. Overall, results of PDDCSs of the product alignment and product orientation at these collision energies in the title reaction are not very strongly distinguishable.     

H+HF(v,j)→H+HF(v′,j′)的传能动力学理论研究

本文运用非含时量子动力学方法研究了H+HF(v=1,j)→H+HF(v'=0,j')传能过程在295~500 K的振动弛豫速率常数.在此温度范围内,所有转动分辨的振动弛豫速率常数随着温度升高而单调递增,速率常数最大的末态转动量子数随着初态转动量子数的增加而增加.在室温下,振动态分辨的振动弛豫速率常数与实验值符合较好.同时,我们也计算了H+HF(v=1,j)→H+HF(v'=1,j')纯转动传能过程在500 K的速率常数,发现它们整体上比振动弛豫速率常数大了几个数量级,并且△j=-1的速率常数一般大于△j=-1的速率常数.     

Quasi-Classical Trajectory Study on O++DH (v=0, j=0)→OD++H Reaction at Different Collision Energy

The quasi-classical trajectory calculations O⁺+DH (v=0, j=0)→OD⁺+H reactions on the RODRIGO potential energy surface have been carried out to study the isotope effect onstereo-dynamics at the collision energies of 1.0, 1.5, 2.0, and 2.5 eV. The distributions of dihedral angle P(φr) and the distributions of P(θr) are discussed. Furthermore, the angular distributions of the product rotational vectors in the form of polar plot in θr and φr are calculated. The differential cross section shows interesting phenomenon that the reaction is dominated by the direct reaction mechanism. Reaction probability and reaction cross section are also calculated. The calculations indicate that the stereo-dynamics properties of the title reactions are sensitive to the collision energy.     

三维H+H2(v,j)→H2(v′,j′)+H反应中复合态生成及产物转动态分布的量子散射理论研究

用排列通道线性组合-散射波函数(LCAC-SW,linear combination of arrangement channelsscattering wavefunction)量子反应散射方法计算了H+H2(v,j)→H2(v′,j′)+H三维态-态反应几率,分析了反应体系的复合态生成(或能量共振结构),并由产物的转动态分布解释了能量共振的起源来自于平动态-内态之间的干涉效应.     

H+H2(v=0,1)→H2(v′)+H反应截面的电子非绝热三维量子散射近似计算

本文把电子非绝热一维量子散射反应几率和三维量子散射反应截面的近似公式结合起来, 对于反应物分子(H_2)不同的量子振动态(v=0, 1) 分别计算了H+H_2(v=0)→H_2(v′=0, 1)+H和H+H_2(v=1)→H_2(v′=0, 1)+H的平均反应截面σ_0和σ_1, 并同文献上用电子绝热理论计算的结果作了比较, 表明对这类中性原予-分子反应碰撞的过程, 特别是当反应物分子处于振动激发态时, 电子非绝热效应是存在的。     

I+HI(v=0)→IH(v′=0)+I体系散射共振态的理论研究

在 QCISD ( T) / / MP2水平下 ,分别采用 6-3 1 1 ++G* * 基组和 SDD基组对重 -轻 -重反应 I+HI(ν=0 )→ IH(ν′=0 ) +I中的 H和 I的偏分势能面进行了 ab initio计算 ,指认出在 0～ 0 .5 8e V碰撞能范围内所产生的 6个散射共振态为 Feshbach共振 ,并与文献报道的量子散射理论计算与高分辨阈值光分离光谱实验结果进行了比较 .     

H+O2(n0, j0)→HO+O及C+H2(n0, j0)→CH+H体系选态反应截面的过渡态理论计算

本文用微正则过渡态理论计算了H+O_2(n_0,j_0)→HO+O和C+H_2(n_0, j_0)→CH+H在ab initio势能面上的选态反应截面σ_(n_0,j_0); E.分析了势能面性质对反应截面的影响。计算结果表明, 在指定反应物分子的振动态n_0、转动态j_0时, 两个反应体系的反应截面随相对平动能的增加先是增加后是减小(j_0=1, n_0=0除外); 在给定相对平动能和反应物分子的转动态j_0时, 随反应物分子的振动量子数n_0的增加, 两个体系的选态反应截面均有较显著的增加, 在指定相对平动能和反应物分子的振动态n_0时, H+O_2体系的选态反应截面随j_0的变化较为复杂, 而C+H_2体系则比较简单(j_0=1除外)。对于H+O_2反应体系, 本文得到的反应截面与实验结果及准经典轨迹理论的计算结果符合得很好。     

Theoretical Studies of O(1D)+HD (v=0, j=0, 1, 2, 3)→OD(H)+H(D) Reaction

The quasi-classical trajectory calculation for the reaction O(¹D)+HD is carried out based on the Dobbyn and Knowles potential energy surface. In this work, the reaction cross section and product branching ratio are obtained. The product branching ratio OD/OH was discussed. The calculated results show that the cross-section decreases thoroughly with the increasing of the collision energy from 4.6 kJ/mol to 46.0 kJ/mol. The average branching ratio decrease with the increase of rotational quantum number of reactant HD.     

Collision induced absorption in the a1Δ(v = 2) ← X3Σg(-)(v = 0) band of molecular oxygen

Using cavity ring-down spectroscopy we measured the collision induced absorption spectrum associated with the a(1)Δ(v = 2) ←X(3)Σ(g)(-)(v = 0) band of oxygen near 922 nm both in pure oxygen and in mixtures of oxygen and nitrogen. For pure oxygen, we report for this band an integrated absorption of (1.56 - 0.04/+0.40) × 10(-5) cm(-2) amg(-2). We find that collisions between oxygen and nitrogen do not result in any measurable CIA signal. At 1 bar of oxygen, this collision induced transition is much stronger than the allowed magnetic dipole and electric quadrupole transitions.     

准经典轨线研究Li+HF(ν=0, j=0)→LiF+H反应动力学

基于Aguado等人拟合的APW势能面(PES), 运用准经典轨线(QCT)方法, 对反应Li+HF(ν=0, j=0)→LiF+H的动力学性质进行了计算. 主要研究了不同碰撞能条件下的反应截面、转动取向、产物散射角分布和竞争反应模式等. 结果表明, 该反应存在直接提取型和间接插入型两种反应模式, 在低能量下反应以间接插入反应模式为主, 能量大于200 meV时则以直接提取反应为主.     

Kinetics and dynamics study of the H + CCl4?��?HCl(v��, j��)?+?CCl3 reaction

Based on a previous potential energy surface describing the H?+?CCl₄ reaction, a new analytical surface named PES-2010 was developed modifying both the functional form to give it more flexibility, and the calibration process in which exclusively theoretical information was used. Thus, the surface is completely symmetric with respect to the permutation of the four methane chlorine atoms, and no experimental information is used in the process. For the kinetics, the thermal rate constants were calculated using variational transition-state theory with semiclassical transmission coefficients over a wide temperature range, 300?C2,500?K. The theoretical results reproduce the experimental variation with temperature. The influence of the tunneling factor is small, since the abstraction reaction involves the motion of a heavy particle (a chlorine atom) that cannot easily tunnel through the reaction barrier. The coupling between the reaction coordinate and the vibrational modes shows qualitatively that the HCl stretching mode in the products appears vibrationally excited. The dynamics study was performed using quasi-classical trajectory calculations, including corrections to avoid the zero-point energy problem. First, we found that the HCl(????, j??) product mostly appears with small rotational energy and vibrational population inversion. Second, the state-specific scattering distributions show backward scattering, which becomes more noticeable as the HCl(????) vibrational state increases. Unfortunately, no experimental dynamics data are available for the title reaction, but the comparison with the kinematically similar and well-studied H?+?Cl₂ reaction shows good agreement, indicative of similar mechanisms. These kinetics and dynamics results seem to indicate that the potential energy surface is adequate to describe this reaction, and the reasonable agreement with experiment lends further confidence to this new surface.     

Scattering Resonance States and Partial Potential Energy Surface of Reaction I+HI(v=0)→IH(v′=0)+I

The partial potential energy surface of the I HI→IH I reaction involving the translational and vibrational motions has been constructed at the QCISD(T)//MP4SDQ level with the pseudo potential method that is helpful to interpreting the scattering resonance states. The lifetimes of the scattering resonance states in the title reaction obtained from the partial potential energy surface are about 90-120 fs, which agrees with the result of high-resolved threshold photodetachment spectroscopy of anion IHI- measured by Neumark.     

Scattering Resonance States and Partial Potential Energy Surface of Reaction I+HI(v=0)→IH(v′=0) +I

Introduction Scattering resonance states in chemical reactionsattract much interest of both theorists and experimental-ists because the rate of a chemical reaction, state distri-bution and space distribution of the products are con-trolled by those states[1]. Scattering resonance states,also called quasi-bound states, are formed in theprocesses of energy transitions, such as collision andray radiation processes. Those states are very closelyrelated to the structure and spectroscopy of the tran…     

反应物振动激发对反应O++DH(V=0,1,2,3,j=0)→OD++H 影响的准经典轨线研究

运用准经典轨线方法,基于RODRIGO势能面,对碰撞能为192.97 kJ/mol时,O++DH反应的立体动力学性质进行了理论研究,对k-j',k-k'两矢量相关和k-k'-j'三矢量相关的分布函数、反应几率、反应截面,以及产物转动取向参数进行了详细的讨论.结果表明,反应物的振动量子数对o++DH (v=0,1,2,3...     

Determination of the absolute scattering cross section for the reaction D + H2(v=1) → HD+H at 0.33 eV

The reaction D+H₂(v=1) has been investigated in a crossed molecular beam experiment at the most probable collision energy of E₀=0.33 eV. Angular and time-of-flight distributions have been measured and the total absolute cross section has been determined to be σ_r(j_i=0, v = 1, E_c.m. = 0.33 eV) = 1.14 ± 0.50 Ų. This value, as well as the distributions, are in good agreement with the results of quasiclassical trajectory calculations (QCT) and the reactive infinite-order sudden approximation (RIOSA).