搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

O++H2 → OH++H反应的动力学研究

袁美玲 李文涛

引用本文:
Citation:

O++H2 → OH++H反应的动力学研究

袁美玲, 李文涛

Dynamics studies of O+ + H2→ OH+ + H reaction

Yuan Mei-Ling, Li Wen-Tao
PDF
HTML
导出引用
  • 采用含时量子波包方法结合二阶分裂算符传播子对初始态为(v = 0, j = 0)的O+ + H2 → OH+ + H反应体系在0.01—1.00 eV的碰撞能范围内进行了态分辨理论水平上的动力学计算. 对反应概率、积分截面、微分截面以及固定初始态的热速率常数等动力学信息进行了计算并与文献报道的实验和理论结果进行了比较. 结果表明本文的理论结果与实验结果十分符合. 从微分截面的散射信息可知, 在低碰撞能范围内, 插入反应机制在反应中占据主导地位. 随着碰撞能的增加, 反应机制逐渐由插入机制变为抽取反应机制.
    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.
      通信作者: 李文涛, wtlee1982@163.com
    • 基金项目: 辽宁省教育厅青年项目(批准号: LQ2017001)资助的课题.
      Corresponding author: Li Wen-Tao, wtlee1982@163.com
    • Funds: Project supported by the Youth Fund of Education Department of Liaoning Province, China (Grant No. LQ2017001).
    [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 3346Google Scholar

    [4]

    Wang B B, Han Y C, Gao W, Cong S L 2017 Phys. Chem. Chem. Phys. 19 22926Google Scholar

    [5]

    Fehsenfeld F C, Schmeltekopf A L, Ferguson E E 1967 J. Chem. Phys. 46 2802Google Scholar

    [6]

    Kim J K, Theard L P, Huntress W T 1975 J. Chem. Phys. 62 45Google Scholar

    [7]

    Federer W, Villinger H, Howorka F, Lindinger W, Tosi P, Bassi D, Ferguson E 1984 Phys. Rev. Lett. 52 2084Google Scholar

    [8]

    Smith D, Adams N G, Miller T M 1978 J. Chem. Phys. 69 308Google Scholar

    [9]

    Burley J, Ervin K M, Armentrout P 1987 Int. J. Mass Spectrom. Ion Process 80 153Google Scholar

    [10]

    Sunderlin L, Armentrout P 1990 Chem. Phys. Lett. 167 188Google Scholar

    [11]

    Flesch G D, Ng C Y 1991 J. Chem. Phys. 94 2372Google Scholar

    [12]

    Li X, Huang Y L, Flesch G D, Ng C Y 1997 J. Chem. Phys. 106 564Google Scholar

    [13]

    Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 58 5373Google Scholar

    [14]

    Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 59 6380Google Scholar

    [15]

    Ng C Y 2002 J. Phys. Chem. A 106 5953Google Scholar

    [16]

    Martínez R, Millán J, González M 2004 J. Chem. Phys. 120 4705Google Scholar

    [17]

    Martínez R, Sierra J D, González M 2005 J. Chem. Phys. 123 174312Google Scholar

    [18]

    Martínez R, Lucas J M, Giménez X, Aguilar A, González M 2006 J. Chem. Phys. 124 144301Google Scholar

    [19]

    Martínez R, Sierra J D, Gray S K, González M 2006 J. Chem. Phys. 125 164305Google Scholar

    [20]

    Kłos J, Bulut N, Akpinar S 2012 Chem. Phys. Lett. 532 22Google Scholar

    [21]

    Xu W, Li W, Lv S, Zhai H, Duan Z, Zhang P 2012 J. Phys. Chem. A 116 10882Google 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 33Google 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 11951Google 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 1039Google Scholar

    [25]

    段志欣, 邱明辉, 姚翠霞 2014 物理学报 63 063402Google Scholar

    Duan Z X, Qiu M H, Yao C X 2014 Acta Phys. Sin. 63 063402Google Scholar

    [26]

    李文涛, 于文涛, 姚明海 2018 物理学报 67 103401Google Scholar

    Li W T, Yu W T, Yao M H 2018 Acta Phys. Sin. 67 103401Google Scholar

    [27]

    张静, 魏巍, 高守宝, 孟庆田 2015 物理学报 64 063101Google Scholar

    Zhang J, Wei W, Gao S B, Meng Q T 2015 Acta Phys. Sin. 64 063101Google Scholar

    [28]

    Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 214303Google Scholar

    [29]

    Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 064104Google Scholar

    [30]

    Li W T, Chen M D, Sun Z G 2015 Chin. J. Chem. Phys. 28 415Google Scholar

    [31]

    Fleck J A, Morris J R, Feit M D 1976 Appl. Phys. 10 129

  • 图 1  新旧两个势能面的长程相互作用势比较

    Fig. 1.  Comparison of the long range interaction potential for the old and new potential energy surfaces

    图 2  90°, 120°, 150°和180°的最小能量路径

    Fig. 2.  The minimum energy paths at 90°, 120°, 150°, and 180° angles

    图 3  对于几个选定的总角动量J在能量范围0—1.0 eV的反应概率

    Fig. 3.  The total reaction probabilities for several selected total angular momentum J in the collision energy from 0 to 1.0 eV

    图 4  总角动量J = 0反应概率的振动分辨

    Fig. 4.  The vibrational resolved reaction probabilities of the total angular momentum J = 0

    图 5  O+ + H2(v = 0, j = 0) → OH+ (v' = 0, j' ) + H反应在总角动量J = 0时转动分辨的反应概率

    Fig. 5.  Rotationally resolved reaction probabilities for the O+ + H2 (v = 0, j = 0) → OH+(v' = 0, j' ) + H reaction at total angular momentum J = 0

    图 6  O+ + H2反应的积分截面以及文献[22]的理论结果和文献[8]的实验结果

    Fig. 6.  The integral cross section of the O+ + H2 reaction as well as the theoretical and experimental results obtained from Ref. [22] and Ref. [8], respectively

    图 7  O+ + H2反应若干能量点的微分截面 (a)低能部分; (b) 高能部分

    Fig. 7.  The differential cross sections of O+ + H2 reaction for several collision energies: (a) The low energy; (b) the high energy part

    图 8  O+ + H2反应的速率常数和文献[5]的实验结果

    Fig. 8.  The rate constant of the O+ + H2 reaction and the experimental values obtained from Ref. [5] in the temperature range of 200 to 1000 K

    表 1  计算中使用的参数(除了特殊声明, 均采用原子单位a.u.)

    Table 1.  Parameters used in the calculation (The atomic unit is used in the calculation unless otherwise stated)

    格点范围和大小$R \in \left[ {0.1,16} \right],{N_R} = 279,N_R^{{\rm{int}}} = 159$
    $r \in \left[ {0.1,15} \right],{N_r} = 279,N_r^{{\rm{asy}}} = 159$
    初始波包Rc = 11.0, ${k_0} = \sqrt {2{E_0}{\mu_R}} $
    ${\varDelta _R} = 0.2 {\rm{exp}}\left[ { - \dfrac{{{{\left( {R - {R_{\rm{c}}}} \right)}^2}}}{{2\varDelta _R^2}}} \right]{\rm{cos}}\left( {{k_0}R} \right)$;

    其中E0 = 0.5 eV
    总传播时间 30000
    最大总角
    动量J
    70
    下载: 导出CSV
  • [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 3346Google Scholar

    [4]

    Wang B B, Han Y C, Gao W, Cong S L 2017 Phys. Chem. Chem. Phys. 19 22926Google Scholar

    [5]

    Fehsenfeld F C, Schmeltekopf A L, Ferguson E E 1967 J. Chem. Phys. 46 2802Google Scholar

    [6]

    Kim J K, Theard L P, Huntress W T 1975 J. Chem. Phys. 62 45Google Scholar

    [7]

    Federer W, Villinger H, Howorka F, Lindinger W, Tosi P, Bassi D, Ferguson E 1984 Phys. Rev. Lett. 52 2084Google Scholar

    [8]

    Smith D, Adams N G, Miller T M 1978 J. Chem. Phys. 69 308Google Scholar

    [9]

    Burley J, Ervin K M, Armentrout P 1987 Int. J. Mass Spectrom. Ion Process 80 153Google Scholar

    [10]

    Sunderlin L, Armentrout P 1990 Chem. Phys. Lett. 167 188Google Scholar

    [11]

    Flesch G D, Ng C Y 1991 J. Chem. Phys. 94 2372Google Scholar

    [12]

    Li X, Huang Y L, Flesch G D, Ng C Y 1997 J. Chem. Phys. 106 564Google Scholar

    [13]

    Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 58 5373Google Scholar

    [14]

    Gillen K T, Mahan B H, Winn J S 1973 J. Chem. Phys. 59 6380Google Scholar

    [15]

    Ng C Y 2002 J. Phys. Chem. A 106 5953Google Scholar

    [16]

    Martínez R, Millán J, González M 2004 J. Chem. Phys. 120 4705Google Scholar

    [17]

    Martínez R, Sierra J D, González M 2005 J. Chem. Phys. 123 174312Google Scholar

    [18]

    Martínez R, Lucas J M, Giménez X, Aguilar A, González M 2006 J. Chem. Phys. 124 144301Google Scholar

    [19]

    Martínez R, Sierra J D, Gray S K, González M 2006 J. Chem. Phys. 125 164305Google Scholar

    [20]

    Kłos J, Bulut N, Akpinar S 2012 Chem. Phys. Lett. 532 22Google Scholar

    [21]

    Xu W, Li W, Lv S, Zhai H, Duan Z, Zhang P 2012 J. Phys. Chem. A 116 10882Google 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 33Google 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 11951Google 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 1039Google Scholar

    [25]

    段志欣, 邱明辉, 姚翠霞 2014 物理学报 63 063402Google Scholar

    Duan Z X, Qiu M H, Yao C X 2014 Acta Phys. Sin. 63 063402Google Scholar

    [26]

    李文涛, 于文涛, 姚明海 2018 物理学报 67 103401Google Scholar

    Li W T, Yu W T, Yao M H 2018 Acta Phys. Sin. 67 103401Google Scholar

    [27]

    张静, 魏巍, 高守宝, 孟庆田 2015 物理学报 64 063101Google Scholar

    Zhang J, Wei W, Gao S B, Meng Q T 2015 Acta Phys. Sin. 64 063101Google Scholar

    [28]

    Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 214303Google Scholar

    [29]

    Zhao B, Sun Z G, Guo H 2016 J. Chem. Phys. 144 064104Google Scholar

    [30]

    Li W T, Chen M D, Sun Z G 2015 Chin. J. Chem. Phys. 28 415Google Scholar

    [31]

    Fleck J A, Morris J R, Feit M D 1976 Appl. Phys. 10 129

  • [1] 李文涛, 袁美玲, 王杰敏. C++H2反应的动力学研究: 基于一个新构建的势能面. 物理学报, 2022, 71(9): 093402. doi: 10.7498/aps.71.20212241
    [2] 袁方园, 朱子亮. D + DBr反应的态-态动力学研究. 物理学报, 2020, 69(11): 113401. doi: 10.7498/aps.69.20200321
    [3] 李文涛, 于文涛, 姚明海. 采用量子含时波包方法研究H/D+Li2LiH/LiD+Li反应. 物理学报, 2018, 67(10): 103401. doi: 10.7498/aps.67.20180324
    [4] 张静, 魏巍, 高守宝, 孟庆田. H+Li2: 一个典型的释能反应体系及其含时动力学研究. 物理学报, 2015, 64(6): 063101. doi: 10.7498/aps.64.063101
    [5] 段志欣, 邱明辉, 姚翠霞. 采用量子波包方法和准经典轨线方法研究S(3P)+HD反应. 物理学报, 2014, 63(6): 063402. doi: 10.7498/aps.63.063402
    [6] 刘丽娟, 颉录有, 陈展斌, 蒋军, 董晨钟. 镁原子碰撞激发微分截面和Stokes参数的理论研究. 物理学报, 2012, 61(10): 103102. doi: 10.7498/aps.61.103102
    [7] 沈光先, 汪荣凯, 令狐荣锋, 周勋, 杨向东. He-HD (HT, DT) 非对称碰撞体系振转势能面及微分散射截面的理论计算. 物理学报, 2012, 61(21): 213101. doi: 10.7498/aps.61.213101
    [8] 李勇军, 冯灏, 孙卫国, 曾阳阳, 王小炼, 李会东, 樊群超. 基于严格交换势的低能电子与H2分子碰撞振动激发散射截面的研究. 物理学报, 2011, 60(4): 043401. doi: 10.7498/aps.60.043401
    [9] 李劲, 令狐荣锋, 司冠杰, 杨向东. 低能He原子与Li2分子碰撞散射截面理论计算. 物理学报, 2010, 59(8): 5424-5428. doi: 10.7498/aps.59.5424
    [10] 王悦, 董德智, 李伟艳, 凤尔银, 崔执凤. He-Na2体系低温下的冷碰撞研究. 物理学报, 2009, 58(10): 6913-6919. doi: 10.7498/aps.58.6913
    [11] 王斌, 冯灏, 孙卫国, 曾阳阳, 戴伟. 低能电子与氢分子碰撞的振动激发积分散射截面的研究. 物理学报, 2009, 58(10): 6932-6937. doi: 10.7498/aps.58.6932
    [12] 余春日, 汪荣凯, 张杰, 杨向东. He同位素原子与HBr分子碰撞的微分截面. 物理学报, 2009, 58(1): 229-233. doi: 10.7498/aps.58.229
    [13] 汪荣凯, 沈光先, 杨向东. He-BH碰撞体系微分截面的理论计算. 物理学报, 2009, 58(8): 5335-5341. doi: 10.7498/aps.58.5335
    [14] 王 平, 李芳昱, 何晓宇. 电磁场中光子-轴子的转化微分截面. 物理学报, 2008, 57(9): 5442-5447. doi: 10.7498/aps.57.5442
    [15] 潘 宇, 王凯俊, 方祯云, 汪先友, 彭庆军. 精确计算n-n重正化链图传播下n+n→2π0反应截面. 物理学报, 2008, 57(8): 4817-4825. doi: 10.7498/aps.57.4817
    [16] 汪荣凯, 沈光先, 宋晓书, 令狐荣锋, 杨向东. He同位素对He-NO碰撞体系微分截面的影响. 物理学报, 2008, 57(7): 4138-4142. doi: 10.7498/aps.57.4138
    [17] 施德恒, 孙金锋, 朱遵略, 杨向东, 刘玉芳, 马 恒. 中、高能电子被SO2分子散射的微分截面、动量转移截面及弹性积分截面. 物理学报, 2007, 56(8): 4435-4440. doi: 10.7498/aps.56.4435
    [18] 汪荣凯, 令狐荣锋, 杨向东. He-NO碰撞体系微分截面的理论计算. 物理学报, 2007, 56(4): 2067-2072. doi: 10.7498/aps.56.2067
    [19] 施德恒, 孙金锋, 朱遵略, 刘玉芳, 杨向东. 中高能电子被O2及CF4分子散射的微分截面、弹性积分截面及动量转移截面. 物理学报, 2005, 54(8): 3548-3553. doi: 10.7498/aps.54.3548
    [20] 白丽华, 张庆刚, 刘新国. D+CD4→CD3+D2反应的四维量子散射计算. 物理学报, 2003, 52(11): 2774-2780. doi: 10.7498/aps.52.2774
计量
  • 文章访问数:  3264
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-05
  • 修回日期:  2019-02-20
  • 上网日期:  2019-04-01
  • 刊出日期:  2019-04-20

O++H2 → OH++H反应的动力学研究

  • 1. 枣庄学院光电工程学院, 枣庄 277160
  • 2. 渤海大学数理系, 锦州 121000
  • 3. 中国科学院大连化学物理研究所, 分子反应动力学国家重点实验室, 大连 116023
  • 通信作者: 李文涛, wtlee1982@163.com
    基金项目: 辽宁省教育厅青年项目(批准号: LQ2017001)资助的课题.

摘要: 采用含时量子波包方法结合二阶分裂算符传播子对初始态为(v = 0, j = 0)的O+ + H2 → OH+ + H反应体系在0.01—1.00 eV的碰撞能范围内进行了态分辨理论水平上的动力学计算. 对反应概率、积分截面、微分截面以及固定初始态的热速率常数等动力学信息进行了计算并与文献报道的实验和理论结果进行了比较. 结果表明本文的理论结果与实验结果十分符合. 从微分截面的散射信息可知, 在低碰撞能范围内, 插入反应机制在反应中占据主导地位. 随着碰撞能的增加, 反应机制逐渐由插入机制变为抽取反应机制.

English Abstract

    • 离子与分子类型反应的研究在星际化学、电荷转移、燃烧化学以及化学键的弱相互作用等领域都有着重要的意义[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小节给出结论.

      图  1  新旧两个势能面的长程相互作用势比较

      Figure 1.  Comparison of the long range interaction potential for the old and new potential energy surfaces

    • 作为一个强有力的计算工具, TDWP方法广泛应用于原子与双原子散射反应[25-27]以及包含多原子的反应[28, 29]中. 关于TDWP方法的详细信息可以参考文献[30], 本文仅简单介绍TDWP方法. 对于一个给定的总角动量J, O+ + H2体系的哈密顿量可以写为

      $ \hat H = - \frac{{{\hbar ^2}}}{{2{{\mu} _R}}}\frac{{{\partial ^2}}}{{\partial {R^2}}} - \frac{{{\hbar ^2}}}{{2{{\mu} _r}}}\frac{{{\partial ^2}}}{{\partial {r^2}}} + \frac{{{{(\hat J - \hat j)}^2}}}{{2{{\mu} _R}{R^2}}} + \frac{{{{\hat j}^2}}}{{2{{\mu} _r}{r^2}}} + V, $

      其中${{\mu} _R}$${{\mu} _r}$分别为Rr方向的约化质量, J是总角动量, j是H2分子的转动角动量. 含时波函数是由体坐标下的平动、振动和转动基矢构成$\left\{ {u_n^v(R){\varphi _v}(r)y_{jK}^{JMp}(\hat R,\hat r)} \right\}$,

      $\begin{split} &\psi _{{v_0}{j_0}{K_0}}^{JMp}(R,r,t) \\[5pt ]= \;& \sum\limits_{n,v,j,K} {F_{nvjK,{v_0}{j_0}{K_0}}^{JMp}} (t)u_n^v(R){\varphi _v}(r) \\ & \times y_{jK}^{JMp}(\hat R,\hat r),\end{split}$

      其中下标n为平动基, p为体系的宇称, (v0, j0, K0)是体系的初始振-转态, M是总角动量J在空间坐标系下z轴的投影.

      在采用TDWP方法计算过程中, 二阶分裂算符[31]用来传播初始波包. 如果$\psi \left( {R,r,t} \right)$t时刻的波包, 那么经过一个时间步, 波包可以写为

      $\begin{split} &\psi (R,r,t) \\ =\; &{{\rm{e}}^{ - {\rm{i}}V{\varDelta _t}/(2\hbar) }}{{\rm{e}}^{ - {\rm{i}}{T_j}{\varDelta _t}/(2\hbar) }}{{\rm{e}}^{ - {\rm{i}}({T_r} + {T_R}){\varDelta _t}/\hbar }} \\ &\times {{\rm{e}}^{ - {\rm{i}}{T_j}{\varDelta _t}/(2\hbar) }}{{\rm{e}}^{ - {\rm{i}}V{\varDelta _t}/(2\hbar) }} \times \psi (R,r,t + {\varDelta _t}),\end{split}$

      其中V是体系的势能, ${T_j} = [{(J - j)^2}/(2{{\mu} _R}{R^2})] + $${j^2}/(2{{\mu} _r}{r^2})$是角动量算符, ${T_r} = - \left[ {{\hbar ^2}/(2{{\mu} _r})} \right]\left( {{\partial ^2}/{\partial ^2}r} \right)$${T_R} = - \left[ {{\hbar ^2}/(2{{\mu} _R)}} \right]\left( {{\partial ^2}/{\partial ^2}R} \right)$分别为rR方向的径向动能算符.

      最终, 通过传播初始波包和抽取态-态散射矩阵, 可以得到体系的反应概率、积分截面、微分截面和速率常数:

      $ P_{{v_0}{j_0}{k_0}}^J(E) = \frac{\hbar }{{{{\mu} _R}}}{\mathop{\rm Im}\nolimits} \left[ {\left\langle {\psi (E)\left| {\delta (r - {r_0})\frac{\partial }{{\partial r}}} \right|\psi (E)} \right\rangle } \right]\!, $

      $\begin{split} {\sigma _{{\upsilon _v}{j_v} \leftarrow {\upsilon _0}{j_0}}} =\; & \frac{{\text{π}}}{{(2{j_0} \!+\! 1)k_{{\upsilon _0}{j_0}}^2}}\!\sum\limits_{{K_v}} {\sum\limits_{{K_0}} {\sum\limits_J {(2J \!+\! 1)} } }\\ & \times|S_{{\upsilon _v}{j_v}{K_v} \leftarrow {\upsilon _0}{j_0}{K_0}}^J{|^2},\end{split}$

      $\begin{split} &\frac{{{\rm{d}}{\sigma _{{\upsilon _v}{j_v} \leftarrow {\upsilon _0}{j_0}}}(\vartheta ,E)}}{{{\rm{d}}\varOmega }} \\=\; & \frac{1}{{2{j_0} + 1}}\sum\limits_{{K_v}} {\sum\limits_{{K_0}} {|\frac{1}{{2{\rm{i}}k_{{\upsilon _0}{j_0}}^2}}} }\\ &\times\sum\limits_J {(2J + 1)d_{{K_v}{K_0}}^J} (\vartheta )S_{{\upsilon _v}{j_v}{K_v} \leftarrow {\upsilon _0}{j_0}{K_0}}^J{|^2},\end{split}$

      $\begin{split} {k_{{v_0}{j_0}}}(T) = & {\left( {\frac{{8kT}}{{{\text{π}}{\mu} }}} \right)^{1/2}}{\left( {{k_{\rm{B}}}T} \right)^{ - 2}} \\ & \times \int_0^\infty {E{{\rm{e}}^{ - E/{(k_{\rm{B}}}T)}}} {\sigma _{{v_0}{j_0}}}(E){\rm{d}}E,\end{split}$

      其中ψ(E)可以通过对含时波函ψ(t)进行快速傅里叶变换得到, ${{S}}_{{\upsilon _v}{j_v}{K_v} \leftarrow {\upsilon _0}{j_0}{K_0}}^J$是散射矩阵, $\vartheta $为散射角, kB为玻尔兹曼常数.

    • 图2给出了O+ + H2反应势能面在90°, 120°, 150°和180°时的最小能量路径, 其中180°对应线性的O+-H-H结构也是体系的最小能量反应路径. 由图2可知, O+ + H2反应是一个放热反应, 其放热能约为0.54 eV(考虑了零点能). 随着碰撞能的增加, O+ + H2反应可以在偏离最小反应路径的位置进行, 在反应路径上多出一个势垒和一个势阱, 并且随着角度的减少, 势垒高度不断增加. 进一步了解势能面的特性, 读者可以参考文献[24].

      图  2  90°, 120°, 150°和180°的最小能量路径

      Figure 2.  The minimum energy paths at 90°, 120°, 150°, and 180° angles

      为了获得收敛的动力学结果, 在进行动力学计算之前需要对诸多动力学参数进行调试. 本文计算了总角动量J = 0的反应概率, 得到了最优的动力学参数, 并将这些参数用于J大于0的计算中. 具体的参数见表1

      格点范围和大小$R \in \left[ {0.1,16} \right],{N_R} = 279,N_R^{{\rm{int}}} = 159$
      $r \in \left[ {0.1,15} \right],{N_r} = 279,N_r^{{\rm{asy}}} = 159$
      初始波包Rc = 11.0, ${k_0} = \sqrt {2{E_0}{\mu_R}} $
      ${\varDelta _R} = 0.2 {\rm{exp}}\left[ { - \dfrac{{{{\left( {R - {R_{\rm{c}}}} \right)}^2}}}{{2\varDelta _R^2}}} \right]{\rm{cos}}\left( {{k_0}R} \right)$;

      其中E0 = 0.5 eV
      总传播时间 30000
      最大总角
      动量J
      70

      表 1  计算中使用的参数(除了特殊声明, 均采用原子单位a.u.)

      Table 1.  Parameters used in the calculation (The atomic unit is used in the calculation unless otherwise stated)

    • 图3给出了若干总角动量J的反应概率. 由图可知, 在反应概率上有很多共振峰, 并且这些共振峰随着能量和总角动量J的增加而变宽, 振幅变小. 这是由于在最小能量路径上有一个支持很多束缚态和准束缚态的深势阱. 随着碰撞能和J的增加, 中间络合物O+-H-H的寿命变短, 因此反应概率上的共振峰变得不那么明显了. 此外, 还可以发现, 对于J < 40, 随着碰撞能的增加其反应概率逐渐降低. 这是一个典型的放热反应特征, 这是因为随着碰撞能以及离心势的增加O+离子的速率越来越快, H2分子来不及与O+离子形成OH键, O+离子就跑开了, 因此反应概率会越来越低. 而对于J > 60, 其反应概率随着碰撞能的增加逐渐升高. 结合图1, 本文认为这是因为随着碰撞能和离心势的增加, O+ + H2反应已经不局限在最小反应路径上, 反应可以在其他角度的能量路径上进行, 在反应路径上将会有势垒的存在, 这导致反应概率随着碰撞能的增加而增加.

      图  3  对于几个选定的总角动量J在能量范围0—1.0 eV的反应概率

      Figure 3.  The total reaction probabilities for several selected total angular momentum J in the collision energy from 0 to 1.0 eV

      图4给出了J = 0反应概率的振动分辨. 由图可知, 产物OH+ (v′ = 0, 1)的振动态在低能时就是打开的, 随着碰撞能的增加, 其他振动激发态陆续打开, 这表明随着碰撞能的增加有更多的能量转化为体系的内能. 从图4中不难看出, v′ = 0和v′ = 1所占据的反应概率的比重一直是最高的. 随着产物振动态的增加, 其所占的比例越低. 这表明在O+ + H2→ OH+ + H反应过程中振动激发的效率不高, 这可能是因为OH+产物振动态之间的能隙间隔过宽所致.

      图  4  总角动量J = 0反应概率的振动分辨

      Figure 4.  The vibrational resolved reaction probabilities of the total angular momentum J = 0

      图5给出了J = 0反应概率的转动分辨. 图5展示了所有共振结构的特征, 如上所述, 这是由反应路径上的深势阱所致. 产物OH+ (v′= 0)的转动通道在j′= 11之前都是打开的. 在图5中仅列出j′从0到5的反应概率. 对比J = 0反应概率的振动分辨, 发现OH+的转动激发更容易实现. 图5有一个奇怪的现象, 产物OH+(v′= 0, j′ = 1)的激发效率远低于其他转动态. 这可能是过渡态产物的特性导致的.

      图  5  O+ + H2(v = 0, j = 0) → OH+ (v' = 0, j' ) + H反应在总角动量J = 0时转动分辨的反应概率

      Figure 5.  Rotationally resolved reaction probabilities for the O+ + H2 (v = 0, j = 0) → OH+(v' = 0, j' ) + H reaction at total angular momentum J = 0

    • 图6给出了O+ + H2反应体系在0.01—1.0 eV能量范围内的总积分截面. 为了与理论和实验文件进行比较, 由文献[22]获得理论结果和文献[8]得到的实验结果也在图6中给出. 由图6可知, 总积分截面随着碰撞能的增加而降低. 比较文献[22]的理论结果与本文的结果, 当碰撞能小于0.1 eV时, 文献[22]的结果高于本文结果, 当碰撞能大于0.4 eV时, 本文的结果略高于文献[22]的结果. 通过与实验结果的比较可知, 本文结果与实验结果符合得更好, 而文献[22]的结果在低碰撞能范围内略高于实验结果.

      图  6  O+ + H2反应的积分截面以及文献[22]的理论结果和文献[8]的实验结果

      Figure 6.  The integral cross section of the O+ + H2 reaction as well as the theoretical and experimental results obtained from Ref. [22] and Ref. [8], respectively

    • 图7给出了O+ + H2反应若干能量点的微分截面. 由图7可知, 微分散射的信号主要集中在前向(0°)和后向(180°)的方向上. 边缘散射的信号分布很少, 并且随着碰撞能的升高逐渐降低. 从图7左图相应的低能部分可以看出, 前向散射的信号和后向散射的信号都非常明显, 但是偏向于前向散射. 在低碰撞能时, 反应主要是在最低能量路径上进行, 反应路径上的深势阱所产生的中间络合物具有很长的寿命. 因此, 产物既可以向前又可以向后散射, complex-forming反应机制占据主导地位. 从图7右图可以看出, 随着碰撞能的增加, 反应已经不局限于最小能量路径上, 反应路径上出现了一个势垒两个势阱, 前向散射明显高于后向散射. 这说明O+离子碰撞H2分子直接形成OH键并带走一个H原子, 没有经过过多的振荡, 表明反应机制逐渐由间接的插入机制转变为直接的抽取机制.

      图  7  O+ + H2反应若干能量点的微分截面 (a)低能部分; (b) 高能部分

      Figure 7.  The differential cross sections of O+ + H2 reaction for several collision energies: (a) The low energy; (b) the high energy part

    • 图8给出了O+ + H2反应在200—1000 K温度范围内的速率常数. 从图8可以看出, 随着温度的增加, 速率常数以约1.25的因子增加. 通过与实验比较发现, 本文的理论结果明显偏低. 速率常数可以通过积分截面获得, 由图6可以看出本文的积分截面与实验结果十分吻合, 所以预期本文的速率常数结果将会与实验结果有比较好的符合. 然而通过图8可以看出, 两者有较大的差别. 本文认为导致这个差距的原因主要源于和实验数据进行对比需要包含反应物H2分子的转动激发态的影响, 而本文的计算为初始态为(v = 0, j = 0)的动力学计算, 没有包含转动激发的影响.

      图  8  O+ + H2反应的速率常数和文献[5]的实验结果

      Figure 8.  The rate constant of the O+ + H2 reaction and the experimental values obtained from Ref. [5] in the temperature range of 200 to 1000 K

    • 本文采用含时量子波包方法和二阶分裂算符传播子, 并基于新势能面进行O+ + H2 → OH+ + H反应的态-态动力学计算. 对反应概率、积分截面、微分截面和速率常数等动力学信息进行了计算, 并与之前的理论和实验结果进行了比较. 结果表明本文得到的积分截面与实验符合得很好, 而本文的速率常数值低于实验值. 从微分截面得到的信息可知, 在反应过程中的低能部分complex-forming反应机制占据主导地位. 随着碰撞能的升高, 反应机制由间接变为直接, 在高能部分抽取反应机制占据主导地位.

参考文献 (31)

目录

    /

    返回文章
    返回