下载:
-
自1941年Douglas和Herzberg[1]在星际介质中发现CH+离子的高丰度含量后, CH+离子受到了大量的实验[2-6]和理论[7-29]研究. CH+离子是高活性的, 在宇宙环境中很容易与电子、氢原子、氢分子等发生反应, 从而被迅速破坏掉. 因此, CH+离子在宇宙中的高丰度就变得难以解释, 这也是现代天体物理学中的一个难题. 为解释其在宇宙中高丰度的现象, 在接下来的几十年中, 人们对CH+离子的形成和破坏机制进行了广泛的研究. 在实验方面, 1997年, Hierl等[2]采用高温流动余热装置测量了C+(2P)和H2分子及其同位素D2分子反应在400—1300 K内的热速率常数. 实验结果表明, 反应物分子的转动能和平动能对反应的影响类似, 而反应物分子的振动激发对反应的影响远大于平动能和转动能. 与反应物的振动基态v = 0比较, v = 1的振动激发, 使H2分子速率常数增大约1000倍, 使D2分子速率常数增大约6000倍. 2011年, Plasil等[3]采用温度可变的22极离子阱和冷喷氢原子束, 研究了CH+离子与H原子的碰撞反应. 他们报道了该反应在10—1000 K温度范围内的速率常数, 发现当温度低于100 K时, 速率常数会出现一个迅速的下降. 这说明非转动的CH+离子在低温情况下很难与H原子发生碰撞反应.
在理论方面, 大量的理论研究工作者采用构建高精度势能面并进行动力学计算的方式来获得指定电子态的速率常数. 2005年, Stoecklin和Halvick[7]采用多项式拟合的方法构建了
$ {\mathrm{C}\mathrm{H}}_{2}^{+} $ 体系的基态势能面, 并基于该势能面采用RIOSA-NIP方法计算了H + CH+反应的反应概率、积分截面和速率常数. 2007年, Halvick等[8]基于该势能面采用准经典轨线方法(quasiclassical trajectory method, QCT)和相空间理论重新研究了H + CH+反应, 并报道了该反应的速率常数. 2013年, Zanchet等[9]基于该势能面采用含时量子波包方法研究了C+ + H2反应, 并报道了该反应量子态分辨的速率常数. 2011年, Warmbier和Schneider[10]报道了${\rm CH}_{2}^{+} $ 体系的一个新势能面. 他们采用AVTZ和AVQZ基组和多参考组态相互作用方法(multireference configuration interaction method, MRCI)进行了从头算计算. 此外, 基于该势能面还采用了QCT方法和ABC非含时的量子方法对H + CH+反应进行了动力学计算, 并报道了积分截面和速率常数等动力学信息. 2015年, Li等[11]采用MRCI方法和AV6Z基组进行了高精度的从头算计算, 并采用多项式拟合的方法构建了CH${}_{2}^{+} $ 体系高精度的势能面. 该势能面得到广泛的应用. Guo等[12], Sundaram等[13-15]和 Wu等[16]采用该势能面对H + CH+反应及其逆反应进行了动力学计算并报道了积分截面、速率常数等动力学信息. 最新的势能面是由Guo等[17]构建的, 在从头算计算中, 他们把基组外推到了完全基组极限, 同样采用多项式拟合的方法构建了新的势能面.综上可知, 研究者们对
$ {\rm CH}_{2}^{+} $ 体系进行了大量的理论研究. 大部分势能面都是采用多项式拟合的方式构建的, 而神经网络拟合方法与传统的拟合方法相比具有精度高、运算快的优势. 本文的目标就是构建更高精度的势能面并在其基础上进行动力学计算. -
本文采用Molpro程序包进行从头算计算, 并且对所有的构型都采用了Cs对称性. 在从头算计算中, 采用了一种明显依赖于电子-电子距离的内收缩多参考组态相互作用方法(MRCI-F12)[30]. 此外, 对C+离子和H原子采用了ACVQZ相关一致基组[31,32]. 对于所有构型, 首先执行包含4个等权重电子态(3个2A', 1个2A'')的态平均完全活化空间自洽场(state average-complete active space self-consistent field, SA-CASSCF)[33,34]计算, 然后利用SA-CASSCF计算提供的参考波函数进行了MRCI-F12计算. 为了补偿高阶截断误差, 计算中考虑了Davidson修正. 在计算中采用雅克比坐标来排布空间格点(RQ, RHH, α). 其中, RQ是C+离子到H2分子的中心距离, 在0.6—25 Bohr的范围内放置55个非均匀格点; RHH是H2分子的键长, 在0.6—25 Bohr的范围内放置56个非均匀格点; α是RQ和RHH之间的夹角, 其变化范围是1°到89°, 间隔为15°. 在拟合过程中, 删掉了能量高于基态最低能量10 eV的从头算能量点. 最终, 在拟合过程中采用了18222个从头算能量点.
-
置换不变多项式神经网络(permutive invariant polynomial-neuronal network, PIP-NN)[35,36]方法已经广泛地应用于势能面构建, 例如H + HBr[37]和K + H2反应[38]等. 在此, 简要地介绍一下PIP-NN方法. 众所周知, 神经网络方法由3个部分组成, 即输入端、 隐藏层、 输出端. 对于输入端, 为了避免在边界处一阶导数不连续, 采用PIP方法对雅克比坐标进行如下变换:
$ {G}_{1}=\dfrac{{P}_{{\mathrm{C}\mathrm{H}}_{a}^+}+{P}_{{\mathrm{C}\mathrm{H}}_{b}^+}}{2},$ $ {G}_{2}={P}_{{\mathrm{C}\mathrm{H}}_{a}^+}\cdot {P}_{{\mathrm{C}\mathrm{H}}_{b}^+}, $ $ {G}_{3}={P}_{\mathrm{H}\mathrm{H}}, $ 式中,
${P}_{{\mathrm{C}\mathrm{H}}_{a}^{+}}=1/{R}_{{\mathrm{C}\mathrm{H}}_{a}^{+}}$ ,${P}_{{\mathrm{C}\mathrm{H}}_{b}^{+}}=1/{R}_{{\mathrm{C}\mathrm{H}}_{b}^{+}}$ ,${P}_{\mathrm{H}\mathrm{H}}=1/{R}_{\mathrm{H}\mathrm{H}}$ , 其中${R}_{\mathrm{H}\mathrm{H}}$ 为氢分子之间的键长,${R}_{{\mathrm{C}\mathrm{H}}_{a}^{+}}, {R}_{{\mathrm{C}\mathrm{H}}_{b}^{+}}$ 为CH+分子的键长. 在计算中, 采用了2个隐藏层, 每个隐藏层包含15个神经元. 此外, Levenberg-Marquardt方法[39]用来训练神经网络. 神经网络的表达式为$ \begin{split} {V}_{{\mathrm{C}\mathrm{H}}_{2}^+}=\;&\sum _{k=1}^{15}{w}_{k}^{3}\cdot {f}^{2}\cdot \left[\sum _{j=1}^{15}{w}_{j}^{2}\cdot {f}^{1}\right.\\ &\left.\times \left(\sum _{i=1}^{{N}_{\mathrm{t}\mathrm{o}\mathrm{t}}}{w}_{i}^{1}\cdot {G}_{i}+{b}_{j}^{1}\right)+{b}_{k}^{2}\right]+{b}_{1}^{3},\end{split} $ 其中,
$ {N}_{\mathrm{t}\mathrm{o}\mathrm{t}} $ 是从头算能量点的个数$ ;{w}_{i}^{1} $ ,$ {w}_{j}^{2} $ 和$ {w}_{k}^{3} $ 是连接权重$ ;{b}_{1}^{3} $ ,$ {b}_{k}^{2} $ 和$ {b}_{j}^{1} $ 是相位$ ;{f}^{1} $ 和$ {f}^{2} $ 是训练函数, 其数学表达式为$ f\left(x\right)=\frac{{\rm e}^{x}-{\rm e}^{-x}}{{\rm e}^{x}+{\rm e}^{-x}}. $ 过拟合问题一直是神经网络方法中的一个棘手问题, 为了避免过拟合的发生, 将输入数据随机分为3个部分, 即训练部分(90%)、 测试部分(5%)和验证部分(5%). 采用均方根误差(root mean square error, RMSE)来测试势能面的精确度, 它的表达式为
$ {y}_{\mathrm{R}\mathrm{M}\mathrm{S}\mathrm{E}}=\sqrt{\frac{1}{N}\sum _{i=1}^{N}{\left({y}_{i}^{\mathrm{o}\mathrm{u}\mathrm{t}\mathrm{p}\mathrm{u}\mathrm{t}}-{y}_{i}^{ab}\right)}^{2}}, $ 其中, N是总的能量点个数,
$ {y}_{i}^{\mathrm{o}\mathrm{u}\mathrm{t}\mathrm{p}\mathrm{u}\mathrm{t}} $ 是拟合值,$ {y}_{i}^{ab} $ 是从头算能量. 对于本文的$ {\mathrm{C}\mathrm{H}}_{2}^{+} $ 体系, 拟合RMSE仅为1.27 meV. -
测试势能面的精确程度的一个有效途径就是将计算得到的光谱常数与实验进行比对. 当第3个原子远离另外两个原子中心距离为30 Bohr时, 改变另外两个原子之间的间距, 就可以在势能面上得到CH+和H2分子的势能曲线. 将该曲线用最小二乘法拟合得到光谱常数, 与实验结果和其他理论结果一并列入到了表1中.
re /Bohr ωe /cm–1 ωexe /cm–1 βe /cm–1 αe /cm–1 De /eV CH+(${X}^{1}{{\Sigma } }^{+}$) 本文结果 2.136 2860.31 59.32 14.217 0.501 4.257 理论[11] 2.136 2853.03 58.52 14.201 0.489 4.252 理论[25] 2.136 2851.0 58.1 14.199 0.489 4.244 理论[26] 2.144 2849.03 66.45 14.094 0.490 理论[17] 2.135 2861.95 59.63 14.311 0.447 4.257 实验[40] 2.137 2857.56 59.32 14.178 0.495 4.26 H2(${X}^{1}{{\Sigma } }_{\mathrm{g} }^{+}$) 本文结果 1.401 4407.63 139.43 60.85 3.012 4.751 理论[17] 1.401 4404.61 126.64 60.861 2.233 4.749 理论[27] 1.401 4403.60 126.60 60.864 2.232 4.748 理论[28] 1.403 4395.22 126.12 60.735 2.221 4.748 理论[29] 1.401 4389.66 121.56 60.826 3.162 4.711 实验[41] 1.401 4401.21 121.33 60.853 3.062 4.746 表 1 CH+和H2分子的光谱常数
Table 1. Spectroscopic constants of the CH+ and H2 molecules.
如表1所列, 本文拟合得到的CH+和H2分子的光谱数据与实验值十分吻合. 对于CH+离子, 平衡核间距re、谐振频率ωe、离解能De与实验值的差距仅为0.001 Bohr, 2.75 cm–1, 3 meV. 对于H2分子, 平衡核间距re与实验结果一致, 谐振频率ωe和离解能De与实验值的差距为6.42 cm–1和5 meV. 实验与理论的微小差距说明新构建的势能面很好地描述了当第3个原子远离另外两个原子时的势能曲线.
图1和图2展示了这个新构建的
$ {\rm CH}_{2}^{+} $ 势能面的主要地形特征. 由图1可知, 在研究涉及的空间范围内等势线正确并光滑, 其等势线对应反应过程中一些主要的驻点. 图1(a)描述了键角∠HC+H = 136.68°时,$ {R}_{{\mathrm{C}\mathrm{H}}^{+}} $ 键伸缩的等势线图. 很明显, 图1(a)中有一个很深的势阱, 对应整个势能面的最低能量点, 对应的结构为$ {R}_{{\mathrm{C}\mathrm{H}}^{+}} $ = 2.08 Bohr,$ {R}_{\mathrm{H}\mathrm{H}} $ = 3.866 Bohr, 相比C+ + H2通道的能量, 势阱的深度约为4.30 eV. 图1(b)描述了C+离子垂直接近H2分子中心点过程的等势线. 可看到两个驻点分别对应于局域最小值和整体最小值. 其中局域最小值的构型为$ {R}_{{\mathrm{C}\mathrm{H}}^{+}} $ = 2.602 Bohr,$ {R}_{\mathrm{H}\mathrm{H}} $ = 1.695 Bohr. 图2(a)给出了当CH+离子处于平衡核间距的构型下, H原子围绕CH+离子运动的等势线图. 在C+离子一侧有一个明显的深势阱. 当H原子在∠HQC+ = 0°—60°范围接近CH+离子时, 很容易被该势阱捕获, 从而形成长寿命的中间产物. 图2(b)给出了当H2分子处于平衡键长时, C+离子围绕H2分子的等势线图. 在90°附近有一个明显的深势阱. 可以预见: 当C+离子接近H2分子时, 很可能被该势阱吸引, 从而进入T型结构; 对于C+ + H2 → CH+ + H反应, 插入反应机理对该反应起着重要作用.
图 1 (a)键角在136.68°时, 化学键伸缩的等势线图(等势线的起点为–9.1 eV, 间隔为0.4 eV); (b) C+离子以C2v对称性接近H2分子中心的等势线图(等势线的起点为–9.1 eV, 间隔为0.37 eV)
Figure 1. (a) Contour plot for chemical bond stretching, in which the angle is fixed at 136.68° (Contours starting at –9.1 eV and equally spaced by 0.4 eV); (b) contour plot for the C+ ion approach to the midpoint of H2 molecule in the C2v symmetry (Contours starting at –9.1 eV and equally spaced by 0.37 eV).
图 2 (a) 当
$ {R}_{{\mathrm{C}\mathrm{H}}^{+}}=2.136 $ Bohr时, H原子绕CH+离子运动的等势线; (b) 当$ {R}_{\mathrm{H}\mathrm{H}}=1.401 $ Bohr时, C+离子围绕H2分子运动的等势线.Figure 2. (a) Contour plot for the H atom moves around CH+ ion at the bond distance
$ {R}_{{\mathrm{C}\mathrm{H}}^{+}}=2.136 $ Bohr; (b) contour plot for C+ ion moves around H2 molecule at its equilibrium geometry$ {R}_{\mathrm{H}\mathrm{H}}=1.401 $ Bohr.图3给出了CH
${}_2^+ $ 体系基态势能面在不同角度下的最小能量路径, 为了与之前的理论结果进行比较, Li等[11]的结果也列入了图中. 如图3所示, H + CH+通道的能量设置为0. 很明显, C+ + H2 → CH+ + H反应是一个吸热反应, 在不考虑零点能的情况下, 其吸热能约为0.5 eV, 这与Li等[11]报道的结果十分接近. 此外, 在小角度时, 能量路径上存在两个势阱, 分别位于势垒的两侧. 然而随着角度的增加, 势垒消失, 势阱变为一个, 并且势阱的深度持续变浅. 角度为180°时, 在能量路径上出现了一个微小势垒, 该势垒的高度几乎与H + CH+通道齐平. 此外, Li等[11]的最小能量路径与本文的最小能量路径基本一致, 一些细微的差别可能是因为使用不同的拟合方法导致的. -
为了进一步验证势能面的准确性, 基于新构建的势能面, 采用QCT[42]对C+ + H2反应以0.05 eV的步长, 在0.4—1.0 eV的碰撞能范围内进行了动力学计算. 当入射C+原子与H2分子的质心的距离小于8 Å时, 轨线开始计算, 并且对于每一个能量随机采样1 × 105条轨线. 在轨线的计算中, 为了确保数值的稳定性, 采用的积分步长为0.02 fs. 在正式开始计算之前, 对于每一个能量点, 需要对最大碰撞参数bmax进行调节. 在本文的计算中, 采取测试的轨线数目为5 × 103条, 在测试过程中逐渐增大bmax直至没有反应轨线出现. 在积分截面和微分截面的计算中, 碰撞参数b在0和最大碰撞参数bmax之间随机取值. 当产物分子之间的距离大于8 Å时, 轨线终止计算. 积分截面和微分截面可以通过下列公式得到:
$ \sigma =\frac{{R}_{\mathrm{t}\mathrm{r}\mathrm{a}}}{{N}_{\mathrm{t}\mathrm{r}\mathrm{a}}}\cdot \pi \cdot {b}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{2},$ $\begin{split} & \\ &\frac{\mathrm{d}\sigma }{\mathrm{d}\varOmega }=\frac{{R}_{\mathrm{t}\mathrm{r}\mathrm{a}}}{{N}_{\mathrm{t}\mathrm{r}\mathrm{a}}}\cdot \frac{{b}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{2}}{2\mathrm{s}\mathrm{i}\mathrm{n}\left( {\mathrm{\pi }{\varDelta }_{\theta }}/{180}\right)}\cdot \frac{{N}_{\theta }}{\mathrm{\pi }}, \end{split} $ 其中, Rtra是反应轨线数目, Ntra是总轨线数目, Δθ是角度步长, Nθ是角方向格点数目.
图4给出了C+ + H2 → H + CH+反应的积分截面. 为了便于与之前的理论结果进行比较, Guo等[12]和Herráez-Aguilar等[18]的结果也列入到了图中. 如图4所示, 积分截面的大小随着碰撞能的增加而增大, Guo等[12]和Herráez-Aguilar等[18]的结果与本文结果具有相同的趋势. 尽管如此, 本文的结果与Guo等[12]和Herráez-Aguilar等[18]的结果之间还存在着明显的不同. Herráez-Aguilar等[18]的结果是基于Warmbier和Schneider[10]的势能面, 采用QCT得到的. 与本文的结果比较, 在研究的碰撞能范围内, 本文的结果明显大于Herráez-Aguilar等[18]的结果, 这是因为在计算中采用了不同的势能面. Guo等[12]的结果是基于Li等[11]的势能面, 采用含时量子波包方法计算得到的. 然而在他们的计算中, 对角动量J在z轴方向上的投影进行了截断(Kblock = 20). 对于固定的J值, 完全的科里奥利耦合计算在z轴方向上的投影应为J + 1, 这表明Guo等[12]的结果在J ≥ 20时就可能得不到收敛的动力学结果. 因此, Guo等[12]的结果很可能在低能处是收敛的, 而对于较高的碰撞能, 他们的结果很可能不收敛. 本文采用的QCT方法不能考虑低能处的量子效应, 因此, 在低能处, 本文的结果与Guo等[12]的结果有一些差距. 在高能处, 量子效应变得不那么明显, 量子结果应与QCT结果取得类似的结果. 然而, Guo等[12]的结果与本文的QCT结果依然存在着较大的差距, 这可能是因为Guo等[12]截断了J在z轴方向上的投影, 从而导致量子结果不收敛.
图5给出了C+ + H2 → H + CH+反应在不同碰撞能下的微分截面. 从图5不难看出, 对于所有能量点, 散射信号都集中在0°和180°附近, 并且前后向对称. 这说明, 在反应过程中插入反应机理占据主导地位. 反应路径上的深势阱(图3)支持长寿命的束缚态和准束缚态, 从而使得产物可以到达任何可接触到的通道.
图 5 C+ + H2反应若干碰撞能下的微分截面
Figure 5. Differential cross sections of the C+ + H2 reaction at several collision energies.
-
本文采用PIP-NN构建了
$ {\rm CH}_2^+ $ 体系基态势能面, 使用MRCI-F12方法结合ACVQZ基组, 对$ {\rm CH}_2^+ $ 体系进行了从头算计算, 详细描述了势能面的特征, 并与之前的理论报道进行了详细的比较. 结果表明, 新构建的势能面足够精确, 可以用于$ {\rm CH}_2^+ $ 体系一切类型的动力学计算. 为了进一步验证新构建势能面的准确性, 基于新构建的势能面, 采用QCT在0.4—1.0 eV的碰撞能范围内对C+ + H2 → H + CH+反应进行了计算. 报道了总积分截面和微分截面, 并与之前的理论结果进行了比较. 对于积分截面, 本文的结果大体上与之前的理论结果保持一致. 微分截面的结果表明插入反应机理在反应过程中占据主导地位.
-
采用一种明显依赖于电子-电子距离的多参考组态相互作用方法, 结合ACVQZ基组进行了
${\rm CH}_2^+ $ 体系的从头算计算. 采用置换不变多项式神经网络方法对从头算得到的18222个单点能进行了拟合. 此外, 将势能面的特征与可得到的实验结果和理论结果进行了比较. 结果表明新构建的势能面比以往的势能面更加精确. 基于新构建的势能面, 采用准经典轨线方法进行了C+ + H2反应的动力学计算, 报道了积分截面和微分截面等动力学信息, 并与之前的理论结果进行了比较. 动力学结果表明在反应过程中插入反应机理占据主导地位.-
关键词:
- $ {\rm CH}_2^+ $体系 /
- 势能面 /
- 积分截面 /
- 微分截面
The multi-reference interaction method is explicitly dependent on the electron-electron distance, and ACVQZ basis set is used in the ab initio calculation. The potential energy surface (PES) is fitted by using the permutation invariant polynomial neural network method based on 18222 ab initio points. In addition, the topographical features of the PES are compared with available theoretical and experimental data. The results indicate that the present PES is more accurate and can be applied to any type of dynamic study. In order to validate the PES, the dynamic study of the C+ + H2 → H + CH+ reaction is carried out by using the quasi-classical trajectory method in a collision energy range of 0.4–1.0 eV. The integral cross sections and differential cross sections are calculated and compared with previous theoretical studies. For the integral cross section, the present results are, in general, in good agreement with previous theoretical studies, both of which increase with collision energy increasing. The forward and backward symmetric differential cross sections indicate that the “complex-forming” mechanism plays a dominant role in the reaction.-
Keywords:
- $ {\rm CH}_2^+ $ system /
- potential energy surface /
- integral cross section /
- differential cross section
[1] Douglas A E, Herzberg G 1941 Astrophys. J. 94 381
Google Scholar
[2] Hierl P M, Morris R A, Viggiano A A 1997 J. Chem. Phys. 106 10145
Google Scholar
[3] Plasil R, Mehner T, Dohnal P, Kotrik T, Glosik J, Gerlich D 2011 Astrophys. J. 737 60
Google Scholar
[4] Armentrout P B 2000 Int. J. Mass Spect. 200 219
Google Scholar
[5] Luca A, Borodi G, Gerlich D 2006 Photonic, Electronic and Atomic Collisions (Singapore: World Scientific)
[6] Federer W, Villinger H, Howorka F, Lindinger W, Tosi P, Bassi D, Ferguson E 1984 Phys. Rev. Lett. 52 2084
Google Scholar
[7] Stoecklin T, Halvick P 2005 Phys. Chem. Chem. Phys. 7 2446
Google Scholar
[8] Halvick P, Stoecklin T, Larrégaray P, Bonnet L 2007 Phys. Chem. Chem. Phys. 9 582
Google Scholar
[9] Zanchet A, Godard B, Bulut N, Roncero O, Halvick P, Cernicharo J 2013 Astrophys. J. 766 80
Google Scholar
[10] Warmbier R, Schneider R 2011 Phys. Chem. Chem. Phys. 13 10285
Google Scholar
[11] Li Y Q, Zhang P Y, Han K L 2015 J. Chem. Phys. 142 124302
Google Scholar
[12] Guo J, Zhang A J, Zhou Y, Liu J Y, Jia J F, Wu H S 2017 Chem. Phys. Lett. 689 121
Google Scholar
[13] Sundaram P, Manivannan V, Padamanadan R 2017 Phys. Chem. Chem. Phys. 19 20172
Google Scholar
[14] Sundaram P, Padamanadan R 2018 J. Chem. Phys. 148 164306
Google Scholar
[15] Sundaram P, Padamanadan R 2020 J. Phys. B-At. Mol. Opt. 53 105201
Google Scholar
[16] Wu H, Duan Z, Chen G 2020 Chem. Phys. Lett. 755 137783
Google Scholar
[17] Guo L, Ma H Y, Zhang L L, Song Y Z, Li Y Q 2018 RSC Adv. 8 13635
Google Scholar
[18] Herráez-Aguilar D, Jambrina P G, et al. 2014 Phys. Chem. Chem. Phys. 16 24800
Google Scholar
[19] Werfelli G, Halvick P, Honvault P, Kerkeni B, Stoecklin T 2015 J. Chem. Phys. 143 114304
Google Scholar
[20] Bovino S, Grassi T, Gianturco F A 2015 J. Phys. Chem. A 119 11973
Google Scholar
[21] Faure A, Halvick P, Stoecklin T, et al. 2017 Mon. Not. R. Astron. Soc. 469 612
Google Scholar
[22] Guo L, Yang Y F, Fan X X, Ma F C, Li Y Q 2017 Commun. Theor. Phys. 67 549
Google Scholar
[23] Gerlich D, Horning S 1992 Chem. Rev. 92 1509
Google Scholar
[24] Gerlich D, Borodi G, Luca A, Mogo C, Smith M 2011 Z. Phys. Chem. 225 475
Google Scholar
[25] Biglari Z, Shayesteh A, Maghari A 2014 Comput. Theor. Chem. 1047 22
Google Scholar
[26] Reddy R R, Nazeer A Y, Rama G K, Baba B D 2004 J. Quant. Spectrosc. Radiat. Transfer 85 105
Google Scholar
[27] Varandas A J C 1996 J. Chem. Phys. 105 3524
Google Scholar
[28] Song Y Z, Zhang Y, Zhang L L, Gao S B, Meng Q T 2015 Chin. Phys. B 24 063101
Google Scholar
[29] Yang C L, Huang Y J, Zhang X, Han K L 2003 J. Mol. Struct.: Theochem. 625 289
Google Scholar
[30] May A J, Valeev E F, Polly R, Manby F R 2005 Phys. Chem. Chem. Phys. 7 2710
Google Scholar
[31] Dunning T H 1989 J. Chem. Phys. 90 1007
Google Scholar
[32] Kendall R A, Dunning T H, Harrison R J 1992 J. Chem. Phys. 96 6796
Google Scholar
[33] Werner H J, Knowles P J 1985 J. Chem. Phys. 82 5053
Google Scholar
[34] Knowles P J, Werner H J 1985 Chem. Phys. Lett. 115 259
Google Scholar
[35] Jiang B, Li J, Guo H 2016 Int. Rev. Phys. Chem. 35 479
Google Scholar
[36] Jiang B, Guo H 2013 J. Chem. Phys. 139 054112
Google Scholar
[37] Li W T, He D, Sun Z G 2019 J. Chem. Phys. 151 185102
Google Scholar
[38] Li W T, Wang X M, Zhao H L, He D 2020 Phys. Chem. Chem. Phys. 22 16203
Google Scholar
[39] Hagan M T, Menhaj M B 1994 IEEE Trans. Neural Netw. Learn. Syst. 5 989
Google Scholar
[40] Hakalla R, Kępa R, Szajna W, Zachwieja M 2006 Eur. Phys. J. D 38 481
Google Scholar
[41] Huber K P, Herzberg G 1979 Molecular Spectra and Molecular Structure IV: Constants of Diatomic Molecules (New York: Van Nostrand Reinhold)
[42] Hase W L 1998 Classical Trajectory Simulations: Initial Conditions, A Chapter in Encyclopedia of Computational Chemistry (Vol. 1) (New York: Wiley) pp 399–402
-
图 1 (a)键角在136.68°时, 化学键伸缩的等势线图(等势线的起点为–9.1 eV, 间隔为0.4 eV); (b) C+离子以C2v对称性接近H2分子中心的等势线图(等势线的起点为–9.1 eV, 间隔为0.37 eV)
Fig. 1. (a) Contour plot for chemical bond stretching, in which the angle is fixed at 136.68° (Contours starting at –9.1 eV and equally spaced by 0.4 eV); (b) contour plot for the C+ ion approach to the midpoint of H2 molecule in the C2v symmetry (Contours starting at –9.1 eV and equally spaced by 0.37 eV).
图 2 (a) 当
$ {R}_{{\mathrm{C}\mathrm{H}}^{+}}=2.136 $ Bohr时, H原子绕CH+离子运动的等势线; (b) 当$ {R}_{\mathrm{H}\mathrm{H}}=1.401 $ Bohr时, C+离子围绕H2分子运动的等势线.Fig. 2. (a) Contour plot for the H atom moves around CH+ ion at the bond distance
$ {R}_{{\mathrm{C}\mathrm{H}}^{+}}=2.136 $ Bohr; (b) contour plot for C+ ion moves around H2 molecule at its equilibrium geometry$ {R}_{\mathrm{H}\mathrm{H}}=1.401 $ Bohr.
图 5 C+ + H2反应若干碰撞能下的微分截面
Fig. 5. Differential cross sections of the C+ + H2 reaction at several collision energies.
表 1 CH+和H2分子的光谱常数
Table 1. Spectroscopic constants of the CH+ and H2 molecules.
re /Bohr ωe /cm–1 ωexe /cm–1 βe /cm–1 αe /cm–1 De /eV CH+(${X}^{1}{{\Sigma } }^{+}$) 本文结果 2.136 2860.31 59.32 14.217 0.501 4.257 理论[11] 2.136 2853.03 58.52 14.201 0.489 4.252 理论[25] 2.136 2851.0 58.1 14.199 0.489 4.244 理论[26] 2.144 2849.03 66.45 14.094 0.490 理论[17] 2.135 2861.95 59.63 14.311 0.447 4.257 实验[40] 2.137 2857.56 59.32 14.178 0.495 4.26 H2(${X}^{1}{{\Sigma } }_{\mathrm{g} }^{+}$) 本文结果 1.401 4407.63 139.43 60.85 3.012 4.751 理论[17] 1.401 4404.61 126.64 60.861 2.233 4.749 理论[27] 1.401 4403.60 126.60 60.864 2.232 4.748 理论[28] 1.403 4395.22 126.12 60.735 2.221 4.748 理论[29] 1.401 4389.66 121.56 60.826 3.162 4.711 实验[41] 1.401 4401.21 121.33 60.853 3.062 4.746 
下载: 导出CSV
-
[1] Douglas A E, Herzberg G 1941 Astrophys. J. 94 381
Google Scholar
[2] Hierl P M, Morris R A, Viggiano A A 1997 J. Chem. Phys. 106 10145
Google Scholar
[3] Plasil R, Mehner T, Dohnal P, Kotrik T, Glosik J, Gerlich D 2011 Astrophys. J. 737 60
Google Scholar
[4] Armentrout P B 2000 Int. J. Mass Spect. 200 219
Google Scholar
[5] Luca A, Borodi G, Gerlich D 2006 Photonic, Electronic and Atomic Collisions (Singapore: World Scientific)
[6] Federer W, Villinger H, Howorka F, Lindinger W, Tosi P, Bassi D, Ferguson E 1984 Phys. Rev. Lett. 52 2084
Google Scholar
[7] Stoecklin T, Halvick P 2005 Phys. Chem. Chem. Phys. 7 2446
Google Scholar
[8] Halvick P, Stoecklin T, Larrégaray P, Bonnet L 2007 Phys. Chem. Chem. Phys. 9 582
Google Scholar
[9] Zanchet A, Godard B, Bulut N, Roncero O, Halvick P, Cernicharo J 2013 Astrophys. J. 766 80
Google Scholar
[10] Warmbier R, Schneider R 2011 Phys. Chem. Chem. Phys. 13 10285
Google Scholar
[11] Li Y Q, Zhang P Y, Han K L 2015 J. Chem. Phys. 142 124302
Google Scholar
[12] Guo J, Zhang A J, Zhou Y, Liu J Y, Jia J F, Wu H S 2017 Chem. Phys. Lett. 689 121
Google Scholar
[13] Sundaram P, Manivannan V, Padamanadan R 2017 Phys. Chem. Chem. Phys. 19 20172
Google Scholar
[14] Sundaram P, Padamanadan R 2018 J. Chem. Phys. 148 164306
Google Scholar
[15] Sundaram P, Padamanadan R 2020 J. Phys. B-At. Mol. Opt. 53 105201
Google Scholar
[16] Wu H, Duan Z, Chen G 2020 Chem. Phys. Lett. 755 137783
Google Scholar
[17] Guo L, Ma H Y, Zhang L L, Song Y Z, Li Y Q 2018 RSC Adv. 8 13635
Google Scholar
[18] Herráez-Aguilar D, Jambrina P G, et al. 2014 Phys. Chem. Chem. Phys. 16 24800
Google Scholar
[19] Werfelli G, Halvick P, Honvault P, Kerkeni B, Stoecklin T 2015 J. Chem. Phys. 143 114304
Google Scholar
[20] Bovino S, Grassi T, Gianturco F A 2015 J. Phys. Chem. A 119 11973
Google Scholar
[21] Faure A, Halvick P, Stoecklin T, et al. 2017 Mon. Not. R. Astron. Soc. 469 612
Google Scholar
[22] Guo L, Yang Y F, Fan X X, Ma F C, Li Y Q 2017 Commun. Theor. Phys. 67 549
Google Scholar
[23] Gerlich D, Horning S 1992 Chem. Rev. 92 1509
Google Scholar
[24] Gerlich D, Borodi G, Luca A, Mogo C, Smith M 2011 Z. Phys. Chem. 225 475
Google Scholar
[25] Biglari Z, Shayesteh A, Maghari A 2014 Comput. Theor. Chem. 1047 22
Google Scholar
[26] Reddy R R, Nazeer A Y, Rama G K, Baba B D 2004 J. Quant. Spectrosc. Radiat. Transfer 85 105
Google Scholar
[27] Varandas A J C 1996 J. Chem. Phys. 105 3524
Google Scholar
[28] Song Y Z, Zhang Y, Zhang L L, Gao S B, Meng Q T 2015 Chin. Phys. B 24 063101
Google Scholar
[29] Yang C L, Huang Y J, Zhang X, Han K L 2003 J. Mol. Struct.: Theochem. 625 289
Google Scholar
[30] May A J, Valeev E F, Polly R, Manby F R 2005 Phys. Chem. Chem. Phys. 7 2710
Google Scholar
[31] Dunning T H 1989 J. Chem. Phys. 90 1007
Google Scholar
[32] Kendall R A, Dunning T H, Harrison R J 1992 J. Chem. Phys. 96 6796
Google Scholar
[33] Werner H J, Knowles P J 1985 J. Chem. Phys. 82 5053
Google Scholar
[34] Knowles P J, Werner H J 1985 Chem. Phys. Lett. 115 259
Google Scholar
[35] Jiang B, Li J, Guo H 2016 Int. Rev. Phys. Chem. 35 479
Google Scholar
[36] Jiang B, Guo H 2013 J. Chem. Phys. 139 054112
Google Scholar
[37] Li W T, He D, Sun Z G 2019 J. Chem. Phys. 151 185102
Google Scholar
[38] Li W T, Wang X M, Zhao H L, He D 2020 Phys. Chem. Chem. Phys. 22 16203
Google Scholar
[39] Hagan M T, Menhaj M B 1994 IEEE Trans. Neural Netw. Learn. Syst. 5 989
Google Scholar
[40] Hakalla R, Kępa R, Szajna W, Zachwieja M 2006 Eur. Phys. J. D 38 481
Google Scholar
[41] Huber K P, Herzberg G 1979 Molecular Spectra and Molecular Structure IV: Constants of Diatomic Molecules (New York: Van Nostrand Reinhold)
[42] Hase W L 1998 Classical Trajectory Simulations: Initial Conditions, A Chapter in Encyclopedia of Computational Chemistry (Vol. 1) (New York: Wiley) pp 399–402
-
[1] 赵文丽, 王永刚, 张路路, 岳大光, 孟庆田. 基于新CH2( \begin{document}${\tilde {\bf{X}}{}^3}{\bf{A''}}$\end{document} \begin{document}${\bf C}{\bf({}^3}{\bf{P})} + {\bf{H}_2(}{\bf X^1}\Sigma _{\bf g}^ + {\bf )} $\end{document} \begin{document}$ \to {\bf H({}^2}{\bf S}) + {\bf CH}{(\bf{}^2}\Pi ) $\end{document} [2] 袁美玲, 李文涛. O++H2 → OH++H反应的动力学研究. 物理学报, 2019, 68(8): 083401. doi: 10.7498/aps.68.20182141 [3] 刘丽娟, 颉录有, 陈展斌, 蒋军, 董晨钟. 镁原子碰撞激发微分截面和Stokes参数的理论研究. 物理学报, 2012, 61(10): 103102. doi: 10.7498/aps.61.103102 [4] 沈光先, 汪荣凯, 令狐荣锋, 周勋, 杨向东. He-HD (HT, DT) 非对称碰撞体系振转势能面及微分散射截面的理论计算. 物理学报, 2012, 61(21): 213101. doi: 10.7498/aps.61.213101 [5] 李勇军, 冯灏, 孙卫国, 曾阳阳, 王小炼, 李会东, 樊群超. 基于严格交换势的低能电子与H2分子碰撞振动激发散射截面的研究. 物理学报, 2011, 60(4): 043401. doi: 10.7498/aps.60.043401 [6] 李劲, 令狐荣锋, 司冠杰, 杨向东. 低能He原子与Li2分子碰撞散射截面理论计算. 物理学报, 2010, 59(8): 5424-5428. doi: 10.7498/aps.59.5424 [7] 王悦, 董德智, 李伟艳, 凤尔银, 崔执凤. He-Na2体系低温下的冷碰撞研究. 物理学报, 2009, 58(10): 6913-6919. doi: 10.7498/aps.58.6913 [8] 王斌, 冯灏, 孙卫国, 曾阳阳, 戴伟. 低能电子与氢分子碰撞的振动激发积分散射截面的研究. 物理学报, 2009, 58(10): 6932-6937. doi: 10.7498/aps.58.6932 [9] 余春日, 汪荣凯, 张杰, 杨向东. He同位素原子与HBr分子碰撞的微分截面. 物理学报, 2009, 58(1): 229-233. doi: 10.7498/aps.58.229 [10] 汪荣凯, 沈光先, 杨向东. He-BH碰撞体系微分截面的理论计算. 物理学报, 2009, 58(8): 5335-5341. doi: 10.7498/aps.58.5335 [11] 潘 宇, 王凯俊, 方祯云, 汪先友, 彭庆军. 精确计算n-n重正化链图传播下n+n→２π0反应截面. 物理学报, 2008, 57(8): 4817-4825. doi: 10.7498/aps.57.4817 [12] 王 平, 李芳昱, 何晓宇. 电磁场中光子-轴子的转化微分截面. 物理学报, 2008, 57(9): 5442-5447. doi: 10.7498/aps.57.5442 [13] 汪荣凯, 沈光先, 宋晓书, 令狐荣锋, 杨向东. He同位素对He-NO碰撞体系微分截面的影响. 物理学报, 2008, 57(7): 4138-4142. doi: 10.7498/aps.57.4138 [14] 余春日, 凤尔银, 程新路, 杨向东. He-HI复合物势能面及微分散射截面的理论研究. 物理学报, 2007, 56(8): 4441-4447. doi: 10.7498/aps.56.4441 [15] 施德恒, 孙金锋, 朱遵略, 杨向东, 刘玉芳, 马 恒. 中、高能电子被SO2分子散射的微分截面、动量转移截面及弹性积分截面. 物理学报, 2007, 56(8): 4435-4440. doi: 10.7498/aps.56.4435 [16] 汪荣凯, 令狐荣锋, 杨向东. He-NO碰撞体系微分截面的理论计算. 物理学报, 2007, 56(4): 2067-2072. doi: 10.7498/aps.56.2067 [17] 韩慧仙, 彭 谦, 文振翼, 王育彬. S2O分子的局域势能面和振动光谱的解析. 物理学报, 2005, 54(1): 78-84. doi: 10.7498/aps.54.78 [18] 施德恒, 孙金锋, 朱遵略, 刘玉芳, 杨向东. 中高能电子被O2及CF4分子散射的微分截面、弹性积分截面及动量转移截面. 物理学报, 2005, 54(8): 3548-3553. doi: 10.7498/aps.54.3548 [19] 王晓艳, 丁世良. 用李代数方法构造四原子分子的势能面. 物理学报, 2004, 53(2): 423-426. doi: 10.7498/aps.53.423 [20] 孙桂华, 杨向东. H+H2反应截面的全量子力学研究. 物理学报, 2002, 51(3): 506-511. doi: 10.7498/aps.51.506
计量
- 文章访问数: 573
- PDF下载量: 24
- 被引次数: 0
