-
In the present work, the long-range interaction potential part of potential energy surface (PES) of OH2+ system is revised and the new resulting PES apparently is more reasonable than the old one in the long-range part. Based on the new PES, the dynamics calculations of O+ +H2→ OH+ + H reaction are carried out at a state-to-state level of theory by using time-dependent quantum wave packet method with second order split operator in a collision energy range from 0.01 to 1.0 eV. The dynamic properties such as reaction probability, ro-vibrational resolved statereaction probability, integral cross section, differential cross section, and state specific rate constant are calculated and compared with available theoretical and experimental results. The results of ro-vibrational resolved state reaction probability reflect some dynamic properties such as resonances which is attributed to the deep well located on the reaction path. The vibrational resolved state reaction probability indicates that the excitation efficiency of the OH+ product is relatively low. The results of integral cross sections indicate that the present results are in better agreement with the experimental values than with previous theoretical calculations, especially in the low collision energy region. However, the state specific rate constant results underestimate the experimental values. The comparison betweenour calculations and the experimental results indicates that the contribution of the rotational excitation of H2 molecule should be included in the calculations. However, only the initial state v = 0, j = 0 is calculated in the present work. We suppose that the deviation of the present results from the experimental data is due to the fact that the rotational excitation of reactant isnot included in the present calculation. The differential cross section signals indicate that the complex-forming reaction mechanism isdominated in the case of low collision energy, but it transforms into abstract reaction mechanism as the collision energy further increases.
1. 引 言
离子与分子类型反应的研究在星际化学、电荷转移、燃烧化学以及化学键的弱相互作用等领域都有着重要的意义[1-3]. 此外, 包含离子的碰撞反应也是解释反应过程中多通道竞争现象的重要反应类型[4]. 由于电子在不同电子态之间的跃迁导致多数的离子与分子类型反应存在透热效应, 因此无论是非绝热势能面的构建还是非绝热动力学的计算都充满了挑战. 而O+ + H2反应是一个比较特殊的离子与分子反应, 在中低碰撞能下, 该反应是在绝热势能面上进行的. 这样的特性使得O+ + H2反应体系无论是在实验上[5-15]还是在理论上[16-24]都受到了广泛的关注.
实验上, 1987年Burley等[9]采用导向离子束质谱测量了基态氧离子与H2分子及其同位素反应的速率常数、分子比以及积分截面等动力学信息. 采用电荷转移探测方法、八极离子囚禁技术和2—4倍的八极装置, Li等[12]在0.02—12.0 eV的碰撞能范围内测量了O+(4S, 2D, 2P)+ H2(D2)反应的积分截面.
理论上, 2004年Martínez等[16]采用599个从头算能量点构建了OH2+体系的解析势能面. Martínez等[16-19]基于该势能面进行了一系列的动力学研究. 2005年, 采用准经典轨线方法(QCT), Martínez等[17]计算了O+与H2分子及其D2, HD同位素的反应. 紧接着采用非含时的量子机制的紧密耦合方法结合超球坐标, Martínez等[18]在0.04—1.24 eV的碰撞能范围内进一步研究了O+ + H2反应. 2006年, 采用含时量子机制的实波包方法, Martínez等[19]重新计算了O+离子与H2分子及其同位素的反应. 文中把反应概率和积分截面与之前应用非含时量子方法计算的结果进行了比较, 以此来检测科里奥利耦合效应在反应过程中的作用. 之后基于Martínez等[16]的势能面, Xu等[21], Gómez-Carrasco等[22], Bulut等[23]使用含时量子波包方法(TDWP)进行了O+ + H2反应体系的动力学计算, 报道了该反应的反应概率、积分截面和速率常数等动力学性质. 由于Martínez等[16]的势能面仅包含599个从头算能量点, 不足以精确描述具有两个势阱一个势垒的复杂反应体系, 2018年, Li等[24]采用置换不变多项式神经网络方法对63000个从头算能量点进行了拟合, 构建了新的势能面. 并采用TDWP方法在新的势能面基础上进行了动力学计算. 文中对反应概率、积分截面和微分截面等动力学信息进行了报道, 并与之前的实验与理论结果进行了比较.
之前的研究主要集中于介绍势能面的性质, 并且对于长程相互作用势的描述还不充分. 本文对该体系在核间距较远的地方增加了一些从头算点并重新拟合构建了新势能面. 图1给出了当氢分子固定在平衡位置, O+离子线性接近氢分子的图像. 由图1不难看出, 新势能面在长程相互作用势上对于老势能面有了很大的改善. 在新的势能面的基础上, 进行了O+ + H2反应的动力学计算. 本文的结构如下: 第2小节简略介绍TDWP方法; 第3小节对动力学结果进行详细的分析和讨论; 第4小节给出结论.
2. 动力学计算方法
作为一个强有力的计算工具, TDWP方法广泛应用于原子与双原子散射反应[25-27]以及包含多原子的反应[28, 29]中. 关于TDWP方法的详细信息可以参考文献[30], 本文仅简单介绍TDWP方法. 对于一个给定的总角动量J, O+ + H2体系的哈密顿量可以写为
ˆH=−ℏ22μR∂2∂R2−ℏ22μr∂2∂r2+(ˆJ−ˆj)22μRR2+ˆj22μrr2+V, (1) 其中
μR 和μr 分别为R和r方向的约化质量, J是总角动量, j是H2分子的转动角动量. 含时波函数是由体坐标下的平动、振动和转动基矢构成{uvn(R)φv(r)yJMpjK(ˆR,ˆr)} ,ψJMpv0j0K0(R,r,t)=∑n,v,j,KFJMpnvjK,v0j0K0(t)uvn(R)φv(r)×yJMpjK(ˆR,ˆr), (2) 其中下标n为平动基, p为体系的宇称, (v0, j0, K0)是体系的初始振-转态, M是总角动量J在空间坐标系下z轴的投影.
在采用TDWP方法计算过程中, 二阶分裂算符[31]用来传播初始波包. 如果
ψ(R,r,t) 是t时刻的波包, 那么经过一个时间步, 波包可以写为ψ(R,r,t)=e−iVΔt/(2ℏ)e−iTjΔt/(2ℏ)e−i(Tr+TR)Δt/ℏ×e−iTjΔt/(2ℏ)e−iVΔt/(2ℏ)×ψ(R,r,t+Δt), (3) 其中V是体系的势能,
Tj=[(J−j)2/(2μRR2)]+ j2/(2μrr2) 是角动量算符,Tr=−[ℏ2/(2μr)](∂2/∂2r) 和TR=−[ℏ2/(2μR)](∂2/∂2R) 分别为r和R方向的径向动能算符.最终, 通过传播初始波包和抽取态-态散射矩阵, 可以得到体系的反应概率、积分截面、微分截面和速率常数:
PJv0j0k0(E)=ℏμRIm[⟨ψ(E)|δ(r−r0)∂∂r|ψ(E)⟩], (4) συvjv←υ0j0=π(2j0+1)k2υ0j0∑Kv∑K0∑J(2J+1)×|SJυvjvKv←υ0j0K0|2, (5) dσυvjv←υ0j0(ϑ,E)dΩ=12j0+1∑Kv∑K0|12ik2υ0j0×∑J(2J+1)dJKvK0(ϑ)SJυvjvKv←υ0j0K0|2, (6) kv0j0(T)=(8kTπμ)1/2(kBT)−2×∫∞0Ee−E/(kBT)σv0j0(E)dE, (7) 其中ψ(E)可以通过对含时波函ψ(t)进行快速傅里叶变换得到,
SJυvjvKv←υ0j0K0 是散射矩阵,ϑ 为散射角, kB为玻尔兹曼常数.3. 结果与讨论
图2给出了O+ + H2反应势能面在90°, 120°, 150°和180°时的最小能量路径, 其中180°对应线性的O+-H-H结构也是体系的最小能量反应路径. 由图2可知, O+ + H2反应是一个放热反应, 其放热能约为0.54 eV(考虑了零点能). 随着碰撞能的增加, O+ + H2反应可以在偏离最小反应路径的位置进行, 在反应路径上多出一个势垒和一个势阱, 并且随着角度的减少, 势垒高度不断增加. 进一步了解势能面的特性, 读者可以参考文献[24].
为了获得收敛的动力学结果, 在进行动力学计算之前需要对诸多动力学参数进行调试. 本文计算了总角动量J = 0的反应概率, 得到了最优的动力学参数, 并将这些参数用于J大于0的计算中. 具体的参数见表1
表 1 计算中使用的参数(除了特殊声明, 均采用原子单位a.u.)Table 1. Parameters used in the calculation (The atomic unit is used in the calculation unless otherwise stated)格点范围和大小 R∈[0.1,16],NR=279,NintR=159 r∈[0.1,15],Nr=279,Nasyr=159 初始波包 Rc = 11.0, k0=√2E0μR ΔR=0.2exp[−(R−Rc)22Δ2R]cos(k0R); 其中E0 = 0.5 eV 总传播时间 30000 最大总角动量J 70 3.1 反应概率
图3给出了若干总角动量J的反应概率. 由图可知, 在反应概率上有很多共振峰, 并且这些共振峰随着能量和总角动量J的增加而变宽, 振幅变小. 这是由于在最小能量路径上有一个支持很多束缚态和准束缚态的深势阱. 随着碰撞能和J的增加, 中间络合物O+-H-H的寿命变短, 因此反应概率上的共振峰变得不那么明显了. 此外, 还可以发现, 对于J < 40, 随着碰撞能的增加其反应概率逐渐降低. 这是一个典型的放热反应特征, 这是因为随着碰撞能以及离心势的增加O+离子的速率越来越快, H2分子来不及与O+离子形成OH键, O+离子就跑开了, 因此反应概率会越来越低. 而对于J > 60, 其反应概率随着碰撞能的增加逐渐升高. 结合图1, 本文认为这是因为随着碰撞能和离心势的增加, O+ + H2反应已经不局限在最小反应路径上, 反应可以在其他角度的能量路径上进行, 在反应路径上将会有势垒的存在, 这导致反应概率随着碰撞能的增加而增加.
图4给出了J = 0反应概率的振动分辨. 由图可知, 产物OH+ (v′ = 0, 1)的振动态在低能时就是打开的, 随着碰撞能的增加, 其他振动激发态陆续打开, 这表明随着碰撞能的增加有更多的能量转化为体系的内能. 从图4中不难看出, v′ = 0和v′ = 1所占据的反应概率的比重一直是最高的. 随着产物振动态的增加, 其所占的比例越低. 这表明在O+ + H2→ OH+ + H反应过程中振动激发的效率不高, 这可能是因为OH+产物振动态之间的能隙间隔过宽所致.
图5给出了J = 0反应概率的转动分辨. 图5展示了所有共振结构的特征, 如上所述, 这是由反应路径上的深势阱所致. 产物OH+ (v′= 0)的转动通道在j′= 11之前都是打开的. 在图5中仅列出j′从0到5的反应概率. 对比J = 0反应概率的振动分辨, 发现OH+的转动激发更容易实现. 图5有一个奇怪的现象, 产物OH+(v′= 0, j′ = 1)的激发效率远低于其他转动态. 这可能是过渡态产物的特性导致的.
3.2 积分截面
图6给出了O+ + H2反应体系在0.01—1.0 eV能量范围内的总积分截面. 为了与理论和实验文件进行比较, 由文献[22]获得理论结果和文献[8]得到的实验结果也在图6中给出. 由图6可知, 总积分截面随着碰撞能的增加而降低. 比较文献[22]的理论结果与本文的结果, 当碰撞能小于0.1 eV时, 文献[22]的结果高于本文结果, 当碰撞能大于0.4 eV时, 本文的结果略高于文献[22]的结果. 通过与实验结果的比较可知, 本文结果与实验结果符合得更好, 而文献[22]的结果在低碰撞能范围内略高于实验结果.
3.3 微分截面
图7给出了O+ + H2反应若干能量点的微分截面. 由图7可知, 微分散射的信号主要集中在前向(0°)和后向(180°)的方向上. 边缘散射的信号分布很少, 并且随着碰撞能的升高逐渐降低. 从图7左图相应的低能部分可以看出, 前向散射的信号和后向散射的信号都非常明显, 但是偏向于前向散射. 在低碰撞能时, 反应主要是在最低能量路径上进行, 反应路径上的深势阱所产生的中间络合物具有很长的寿命. 因此, 产物既可以向前又可以向后散射, complex-forming反应机制占据主导地位. 从图7右图可以看出, 随着碰撞能的增加, 反应已经不局限于最小能量路径上, 反应路径上出现了一个势垒两个势阱, 前向散射明显高于后向散射. 这说明O+离子碰撞H2分子直接形成OH键并带走一个H原子, 没有经过过多的振荡, 表明反应机制逐渐由间接的插入机制转变为直接的抽取机制.
3.4 速率常数
图8给出了O+ + H2反应在200—1000 K温度范围内的速率常数. 从图8可以看出, 随着温度的增加, 速率常数以约1.25的因子增加. 通过与实验比较发现, 本文的理论结果明显偏低. 速率常数可以通过积分截面获得, 由图6可以看出本文的积分截面与实验结果十分吻合, 所以预期本文的速率常数结果将会与实验结果有比较好的符合. 然而通过图8可以看出, 两者有较大的差别. 本文认为导致这个差距的原因主要源于和实验数据进行对比需要包含反应物H2分子的转动激发态的影响, 而本文的计算为初始态为(v = 0, j = 0)的动力学计算, 没有包含转动激发的影响.
4. 结 论
本文采用含时量子波包方法和二阶分裂算符传播子, 并基于新势能面进行O+ + H2 → OH+ + H反应的态-态动力学计算. 对反应概率、积分截面、微分截面和速率常数等动力学信息进行了计算, 并与之前的理论和实验结果进行了比较. 结果表明本文得到的积分截面与实验符合得很好, 而本文的速率常数值低于实验值. 从微分截面得到的信息可知, 在反应过程中的低能部分complex-forming反应机制占据主导地位. 随着碰撞能的升高, 反应机制由间接变为直接, 在高能部分抽取反应机制占据主导地位.
[1] Duley W W, Williams D A 1984 Interstellar Chemistry (New York: Academic Press) p251
[2] Armentrout P 2000 Int. J. Mass Spectrom. Ion Process. 200 21933
[3] Jambrina P G, Alvarino J M, Gerlich D, Hankel M, Herrero V J, Saez-Rabanos V, Aoiz F J 2012 Phys. Chem. Chem. Phys. 14 3346
Google Scholar
[4] Wang B B, Han Y C, Gao W, Cong S L 2017 Phys. Chem. Chem. Phys. 19 22926
Google Scholar
[5] Fehsenfeld F C, Schmeltekopf A L, Ferguson E E 1967 J. Chem. Phys. 46 2802
Google Scholar
[6] Kim J K, Theard L P, Huntress W T 1975 J. Chem. Phys. 62 45
Google Scholar
[7] Federer W, Villinger H, Howorka F, Lindinger W, Tosi P, Bassi D, Ferguson E 1984 Phys. Rev. Lett. 52 2084
Google Scholar
[8] Smith D, Adams N G, Miller T M 1978 J. Chem. Phys. 69 308
Google Scholar
[9] Burley J, Ervin K M, Armentrout P 1987 Int. J. Mass Spectrom. Ion Process 80 153
Google Scholar
[10] Sunderlin L, Armentrout P 1990 Chem. Phys. Lett. 167 188
Google Scholar
[11] Flesch G D, Ng C Y 1991 J. Chem. Phys. 94 2372
Google Scholar
[12] Li X, Huang Y L, Flesch G D, Ng C Y 1997 J. Chem. Phys. 106 564
Google Scholar
[13] Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 58 5373
Google Scholar
[14] Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 59 6380
Google Scholar
[15] Ng C Y 2002 J. Phys. Chem. A 106 5953
Google Scholar
[16] Martínez R, Millán J, González M 2004 J. Chem. Phys. 120 4705
Google Scholar
[17] Martínez R, Sierra J D, González M 2005 J. Chem. Phys. 123 174312
Google Scholar
[18] Martínez R, Lucas J M, Giménez X, Aguilar A, González M 2006 J. Chem. Phys. 124 144301
Google Scholar
[19] Martínez R, Sierra J D, Gray S K, González M 2006 J. Chem. Phys. 125 164305
Google Scholar
[20] Kłos J, Bulut N, Akpinar S 2012 Chem. Phys. Lett. 532 22
Google Scholar
[21] Xu W, Li W, Lv S, Zhai H, Duan Z, Zhang P 2012 J. Phys. Chem. A 116 10882
Google Scholar
[22] Gómez-Carrasco S, Godard B, Lique F, Bulut N, Kłos J, Roncero O, Aguado A, Aoiz F J, Castillo J F, Goicoechea J R 2014 Astrophys. J. 794 33
Google Scholar
[23] Bulut N, Castillo J F, Jambrina P G, Kłos J, Roncero O, Aoiz F J, Bañares L 2015 J. Phys. Chem. A 119 11951
Google Scholar
[24] Li W T, Yuan J C, Yuan M L, Zhang Y, Yao M H, Sun Z G 2018 Phys. Chem. Chem. Phys. 20 1039
Google Scholar
[25] 段志欣, 邱明辉, 姚翠霞 2014 物理学报 63 063402
Google Scholar
Duan Z X, Qiu M H, Yao C X 2014 Acta Phys. Sin. 63 063402
Google Scholar
[26] 李文涛, 于文涛, 姚明海 2018 物理学报 67 103401
Google Scholar
Li W T, Yu W T, Yao M H 2018 Acta Phys. Sin. 67 103401
Google Scholar
[27] 张静, 魏巍, 高守宝, 孟庆田 2015 物理学报 64 063101
Google Scholar
Zhang J, Wei W, Gao S B, Meng Q T 2015 Acta Phys. Sin. 64 063101
Google Scholar
[28] Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 214303
Google Scholar
[29] Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 064104
Google Scholar
[30] Li W T, Chen M D, Sun Z G 2015 Chin. J. Chem. Phys. 28 415
Google Scholar
[31] Fleck J A, Morris J R, Feit M D 1976 Appl. Phys. 10 129
-
表 1 计算中使用的参数(除了特殊声明, 均采用原子单位a.u.)
Table 1. Parameters used in the calculation (The atomic unit is used in the calculation unless otherwise stated)
格点范围和大小 R∈[0.1,16],NR=279,NintR=159 r∈[0.1,15],Nr=279,Nasyr=159 初始波包 Rc = 11.0, k0=√2E0μR ΔR=0.2exp[−(R−Rc)22Δ2R]cos(k0R); 其中E0 = 0.5 eV 总传播时间 30000 最大总角动量J 70 -
[1] Duley W W, Williams D A 1984 Interstellar Chemistry (New York: Academic Press) p251
[2] Armentrout P 2000 Int. J. Mass Spectrom. Ion Process. 200 21933
[3] Jambrina P G, Alvarino J M, Gerlich D, Hankel M, Herrero V J, Saez-Rabanos V, Aoiz F J 2012 Phys. Chem. Chem. Phys. 14 3346
Google Scholar
[4] Wang B B, Han Y C, Gao W, Cong S L 2017 Phys. Chem. Chem. Phys. 19 22926
Google Scholar
[5] Fehsenfeld F C, Schmeltekopf A L, Ferguson E E 1967 J. Chem. Phys. 46 2802
Google Scholar
[6] Kim J K, Theard L P, Huntress W T 1975 J. Chem. Phys. 62 45
Google Scholar
[7] Federer W, Villinger H, Howorka F, Lindinger W, Tosi P, Bassi D, Ferguson E 1984 Phys. Rev. Lett. 52 2084
Google Scholar
[8] Smith D, Adams N G, Miller T M 1978 J. Chem. Phys. 69 308
Google Scholar
[9] Burley J, Ervin K M, Armentrout P 1987 Int. J. Mass Spectrom. Ion Process 80 153
Google Scholar
[10] Sunderlin L, Armentrout P 1990 Chem. Phys. Lett. 167 188
Google Scholar
[11] Flesch G D, Ng C Y 1991 J. Chem. Phys. 94 2372
Google Scholar
[12] Li X, Huang Y L, Flesch G D, Ng C Y 1997 J. Chem. Phys. 106 564
Google Scholar
[13] Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 58 5373
Google Scholar
[14] Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 59 6380
Google Scholar
[15] Ng C Y 2002 J. Phys. Chem. A 106 5953
Google Scholar
[16] Martínez R, Millán J, González M 2004 J. Chem. Phys. 120 4705
Google Scholar
[17] Martínez R, Sierra J D, González M 2005 J. Chem. Phys. 123 174312
Google Scholar
[18] Martínez R, Lucas J M, Giménez X, Aguilar A, González M 2006 J. Chem. Phys. 124 144301
Google Scholar
[19] Martínez R, Sierra J D, Gray S K, González M 2006 J. Chem. Phys. 125 164305
Google Scholar
[20] Kłos J, Bulut N, Akpinar S 2012 Chem. Phys. Lett. 532 22
Google Scholar
[21] Xu W, Li W, Lv S, Zhai H, Duan Z, Zhang P 2012 J. Phys. Chem. A 116 10882
Google Scholar
[22] Gómez-Carrasco S, Godard B, Lique F, Bulut N, Kłos J, Roncero O, Aguado A, Aoiz F J, Castillo J F, Goicoechea J R 2014 Astrophys. J. 794 33
Google Scholar
[23] Bulut N, Castillo J F, Jambrina P G, Kłos J, Roncero O, Aoiz F J, Bañares L 2015 J. Phys. Chem. A 119 11951
Google Scholar
[24] Li W T, Yuan J C, Yuan M L, Zhang Y, Yao M H, Sun Z G 2018 Phys. Chem. Chem. Phys. 20 1039
Google Scholar
[25] 段志欣, 邱明辉, 姚翠霞 2014 物理学报 63 063402
Google Scholar
Duan Z X, Qiu M H, Yao C X 2014 Acta Phys. Sin. 63 063402
Google Scholar
[26] 李文涛, 于文涛, 姚明海 2018 物理学报 67 103401
Google Scholar
Li W T, Yu W T, Yao M H 2018 Acta Phys. Sin. 67 103401
Google Scholar
[27] 张静, 魏巍, 高守宝, 孟庆田 2015 物理学报 64 063101
Google Scholar
Zhang J, Wei W, Gao S B, Meng Q T 2015 Acta Phys. Sin. 64 063101
Google Scholar
[28] Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 214303
Google Scholar
[29] Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 064104
Google Scholar
[30] Li W T, Chen M D, Sun Z G 2015 Chin. J. Chem. Phys. 28 415
Google Scholar
[31] Fleck J A, Morris J R, Feit M D 1976 Appl. Phys. 10 129
期刊类型引用(1)
1. 袁方园,朱子亮. D+DBr反应的态-态动力学研究. 物理学报. 2020(11): 160-167 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 7530
- PDF下载量: 52
- 被引次数: 2