搜索

x

留言板

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

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

HD+分子的强场光解离动力学及其量子调控的理论研究

姚洪斌 蒋相站 曹长虹 李文亮

引用本文:
Citation:

HD+分子的强场光解离动力学及其量子调控的理论研究

姚洪斌, 蒋相站, 曹长虹, 李文亮

Theoretical study of dissociation dynamics of HD+ and its quantum control with an intense laser field

Yao Hong-Bin, Jiang Xiang-Zhan, Cao Chang-Hong, Li Wen-Liang
PDF
HTML
导出引用
  • 利用精确求解原子核与电子耦合运动的三维含时量子波包法, 理论研究了HD+分子在强激光场中的光解离动力学, 并给出了量子调控HD+分子光解离通道的理论方案. 通过分析HD+分子在不同的初始振动态和激光场强度下的光解离动力学过程及其解离核动能谱, 得出了HD+分子的光解离机理及其随激光场强度的变化规律. 研究结果表明, 利用激光场的强度可以实现HD+分子光解离通道的量子调控. 当激光场强度I1 = 4.0 × 1013 W/cm2时, HD+分子的光解离主要是通过净单光子吸收解离和净双光子吸收解离; 当激光场强度增大到I2 = 2.0 × 1014 W/cm2时, 直接双光子吸收解离取代了净单光子吸收解离, 净双光子吸收解离的比重也下降了.
    The dissociation dynamics of HD+ molecule in an intense field is investigated by using an accurate three-dimensional time-dependent wave packet approach. When the 790-nm laser pulse interacts with HD+ molecule, the lowest electronic 1sσ and 2pσ states are coupled. Due to the existence of the permanent electric dipole moment, the transitions in HD+ molecule involve the direct absorption of an odd and even number of photons, thereby opening different pathways for dissociation. The model of the photon-dressed states is presented to analyze the possible dissociation pathways of HD+ molecule. The laser-induced dissociation of HD+ molecule is mainly composed of the four pathways: the direct one-photon absorption, the net two-photon absorption, the direct two-photon absorption, and the direct two-photon absorption. To reveal the dissociation mechanism of HD+ molecule, the kinetic energy resolved spectra are calculated at the given laser intensities. It is found that the dissociation pathways are strongly dependent on laser intensity, especially for the net one-photon absorption dissociation and direct two-photon absorption dissociation. With further research, the dissociation pathways of HD+ are controlled by regulating the intensity of laser pulse. At a laser intensity of 4.0 × 1013 W/cm2, the kinetic energy resolved spectrum for the vibrational state ν = 3 includes the contributions from the net two-photon absorption dissociation and the direct two-photon absorption dissociation. For the vibrational state ν = 6, HD+ molecule is preferentially dissociated via the net one-photon absorption. However, the dissociation mechanism of HD+ molecule at the vibrational states ν = 3 and ν = 6 have significant changes as the laser intensity increases to 2.0 × 1014 W/cm2. For the vibrational state ν = 3, the branching ratio between the dissociation pathway of the net two-photon absorption and that of the direct two-photon absorption has a dramatic change with the increase of laser intensity. Compared with the kinetic energy resolved spectra at laser energy of 4.0 × 1013 W/cm2, the height of the dissociation peak from the net two-photon absorption decreases, and that of the direct two-photon absorption increases at laser intensity of 2.0 × 1014 W/cm2. For the vibrational state ν = 6, the dissociation process of the net one-photon absorption almost disappears at laser intensity of 2.0 × 1014 W/cm2, and it is replaced by the dissociation pathway of the direct two-photon absorption.
      通信作者: 李文亮, wenliangli@vip.126.com
    • 基金项目: 新疆维吾尔自治区高等学校科研计划项目(批准号: XJEDU2016I051, XJEDU2017M042, XJEDU2018I020)、国家自然科学基金(批准号: 11764041, 51462034)、新疆维吾尔自治区天山青年计划(批准号: 2017Q033, 2017Q034)和新疆维吾尔自治区自然科学基金(批准号: 2019D01A32)资助的课题.
      Corresponding author: Li Wen-Liang, wenliangli@vip.126.com
    • Funds: Project supported by the Scientific Research Program of the Higher Education Institution of Xinjiang, China (Grant Nos. XJEDU2016I051, XJEDU2017M042, XJEDU2018I020), the National Natural Science Foundation of China (Grant Nos. 11764041, 51462034), the Tianshan Youth Program of Xinjiang, China (Grant Nos. 2017Q033, 2017Q034), and the Natural Science Foundation of Xinjiang, China (Grant No. 2019D01A32).
    [1]

    秦朝朝, 黄燕, 彭玉峰 2017 物理学报 66 193301Google Scholar

    Qin C C, Huang Y, Peng Y F 2017 Acta Phys. Sin. 66 193301Google Scholar

    [2]

    Kling M F, Siedschlag C, Verhoef A J, Khan J I, Schultze M, Uphues T, Ni Y, Uiberacker M, Drescher M, Krausz F, Vrakking M J J 2006 Science 312 246

    [3]

    Sun Z P, Yao H B, Wang C Y, Zhao, W K, Yang C L 2019 Laser Phys. Lett. 16 016001Google Scholar

    [4]

    Chang Z C, Li C M, Guo W, Yao H B 2018 Chin. Phys. B 27 053301Google Scholar

    [5]

    Yao H B, Guo W, Hoffmann M R, Han K L 2014 Phys. Rev. A 90 063418Google Scholar

    [6]

    刘灿东, 贾正茂, 郑颖辉, 葛晓春, 曾志男, 李儒新 2016 物理学报 65 223206Google Scholar

    Liu C D, Jia Z M, Zheng Y H, Ge X C, Zeng Z N, Li R X 2016 Acta Phys. Sin. 65 223206Google Scholar

    [7]

    姚洪斌, 李文亮, 张季, 彭敏 2014 物理学报 63 178201Google Scholar

    Yao H B, Li W L, Zhang J, Peng M 2014 Acta Phys. Sin. 63 178201Google Scholar

    [8]

    Yao H B, Zhen Y J 2011 Phys. Chem. Chem. Phys. 13 8900Google Scholar

    [9]

    Bucksbaum P H, Zavriyev A, Muller H G, Schumacher D W 1990 Phys. Rev. Lett. 64 1883Google Scholar

    [10]

    Frasinski L J, Posthumus J H, Plumridge J, Codling K, Taday P F, Langley A J 1999 Phys. Rev. Lett. 83 3625Google Scholar

    [11]

    Jolicard G, Atabek O 1992 Phys. Rev. A 46 5845Google Scholar

    [12]

    Frasinski L J, Plumridge J, Posthumus J H, Codling K, Taday P F, Divall E J, Langley A J 2001 Phys. Rev. Lett. 86 2541Google Scholar

    [13]

    Seideman T, Ivanov M Y, Corkum P B 1995 Phys. Rev. Lett. 75 2819Google Scholar

    [14]

    Posthumus J H 2004 Rep. Prog. Phys. 67 623Google Scholar

    [15]

    Orr P A, Williams I D, Greenwood J B, Turcu I C E, Bryan W A, Pedregosa-Gutierrez J, Walter C W 2007 Phys. Rev. Lett. 98 163001Google Scholar

    [16]

    Kiess A, Pavičić D, Hänsch T W, Figger H 2008 Phys. Rev. A 77 053401Google Scholar

    [17]

    McKenna J, Sayler A M, Gaire B, Johnson N G, Zohrabi M, Carnes K D, Esry B D, Ben-Itzhak I 2009 J. Phys. B: At. Mol. Opt. Phys. 42 121003Google Scholar

    [18]

    McKenna J, Sayler A M, Gaire B, Johnson N G, Parke E, Carnes K D, Esry B D, Ben-Itzhak I 2009 Phys. Rev. A 80 023421Google Scholar

    [19]

    Liu Z T, Yuan K J, Shu C C, Hu W H, Cong S L 2010 J. Phys. B: At. Mol. Opt. Phys. 43 055601Google Scholar

    [20]

    He H X, Lu R F, Zhang P Y, Guo Y H, Han K L, He G Z 2011 Phys. Rev. A 84 033418Google Scholar

    [21]

    He H X, Lu R F, Zhang P Y, Han K L, He G Z 2012 J. Chem. Phys. 136 024311Google Scholar

    [22]

    Lu R F, Zhang P Y, Han K L 2008 Phys. Rev. E 77 066701Google Scholar

    [23]

    Hu J, Han K L, He G Z 2005 Phys. Rev. Lett. 95 123001

    [24]

    Feuerstein B, Thumm U 2003 Phys. Rev. A 67 043405Google Scholar

    [25]

    姚洪斌, 张季, 彭敏, 李文亮 2014 物理学报 63 198202Google Scholar

    Yao H B, Zhang J, Peng M, Li W L 2014 Acta Phys. Sin. 63 198202Google Scholar

    [26]

    Yao H B, Zhao G J 2014 J. Phys. Chem. A 118 9173Google Scholar

    [27]

    McKenna J, Sayler A M, Anis F, Gaire B, Johnson N G, Parke E, Hua J J, Mashiko H, Nakamura C M, Moon E, Chang Z, Carnes K D, Esry B D, Ben-Itzhak I 2008 Phys. Rev. Lett. 100 133001Google Scholar

  • 图 1  HD+分子在1sσ和2pσ光子缀饰态势能曲线上的解离示意图[15,25]

    Fig. 1.  The schematic diagram for the dissociation of HD+ molecule at the dressed potentials of electronic 1sσ and 2pσ states.

    图 2  激光场强度I1 = 4.0 × 1013 W/cm2时, HD+分子在初始振动态ν = 3 (a)和ν = 6 (b)上的光解离核动能谱 (激光场的波长λ = 790 nm, 脉冲宽度τ = 40 fs)

    Fig. 2.  At the laser intensity of I1 = 4.0 × 1013 W/cm2, the kinetic energy resolved distributions of dissociation for the vibrational states ν = 3 (a) and ν = 6 (b) of HD+ molecule. The laser wavelength is 790 nm and the pulse duration is 40 fs, respectively.

    图 3  激光场强度I2 = 2.0 × 1014 W/cm2时, HD+分子在初始振动态ν = 3 (a)和ν = 6 (b)上的光解离核动能谱 (激光场的波长λ = 790 nm, 脉冲宽度τ = 40 fs)

    Fig. 3.  At the laser intensity of I2 = 2.0 × 1014 W/cm2, the kinetic energy resolved distributions of dissociation for the vibrational states ν = 3 (a) and ν = 6 (b) of HD+ molecule. The laser wavelength is 790 nm and the pulse duration is 40 fs, respectively.

    图 4  当激光场强度I1 = 4.0 × 1013 W/cm2 (蓝线)和I2 = 2.0 × 1014 W/cm2 (红线)时, HD+分子在初始振动态ν = 6上光解离通道(激光场的波长λ = 790 nm, 脉冲宽度τ = 40 fs)

    Fig. 4.  The related pathways of dissociation for the vibrational state ν = 6 of HD+ molecule at the laser intensities of 4.0 × 1013 W/cm2 (blue line) and 2.0 × 1014 W/cm2 (red line). The laser wavelength is 790 nm, the pulse duration is 40 fs.

  • [1]

    秦朝朝, 黄燕, 彭玉峰 2017 物理学报 66 193301Google Scholar

    Qin C C, Huang Y, Peng Y F 2017 Acta Phys. Sin. 66 193301Google Scholar

    [2]

    Kling M F, Siedschlag C, Verhoef A J, Khan J I, Schultze M, Uphues T, Ni Y, Uiberacker M, Drescher M, Krausz F, Vrakking M J J 2006 Science 312 246

    [3]

    Sun Z P, Yao H B, Wang C Y, Zhao, W K, Yang C L 2019 Laser Phys. Lett. 16 016001Google Scholar

    [4]

    Chang Z C, Li C M, Guo W, Yao H B 2018 Chin. Phys. B 27 053301Google Scholar

    [5]

    Yao H B, Guo W, Hoffmann M R, Han K L 2014 Phys. Rev. A 90 063418Google Scholar

    [6]

    刘灿东, 贾正茂, 郑颖辉, 葛晓春, 曾志男, 李儒新 2016 物理学报 65 223206Google Scholar

    Liu C D, Jia Z M, Zheng Y H, Ge X C, Zeng Z N, Li R X 2016 Acta Phys. Sin. 65 223206Google Scholar

    [7]

    姚洪斌, 李文亮, 张季, 彭敏 2014 物理学报 63 178201Google Scholar

    Yao H B, Li W L, Zhang J, Peng M 2014 Acta Phys. Sin. 63 178201Google Scholar

    [8]

    Yao H B, Zhen Y J 2011 Phys. Chem. Chem. Phys. 13 8900Google Scholar

    [9]

    Bucksbaum P H, Zavriyev A, Muller H G, Schumacher D W 1990 Phys. Rev. Lett. 64 1883Google Scholar

    [10]

    Frasinski L J, Posthumus J H, Plumridge J, Codling K, Taday P F, Langley A J 1999 Phys. Rev. Lett. 83 3625Google Scholar

    [11]

    Jolicard G, Atabek O 1992 Phys. Rev. A 46 5845Google Scholar

    [12]

    Frasinski L J, Plumridge J, Posthumus J H, Codling K, Taday P F, Divall E J, Langley A J 2001 Phys. Rev. Lett. 86 2541Google Scholar

    [13]

    Seideman T, Ivanov M Y, Corkum P B 1995 Phys. Rev. Lett. 75 2819Google Scholar

    [14]

    Posthumus J H 2004 Rep. Prog. Phys. 67 623Google Scholar

    [15]

    Orr P A, Williams I D, Greenwood J B, Turcu I C E, Bryan W A, Pedregosa-Gutierrez J, Walter C W 2007 Phys. Rev. Lett. 98 163001Google Scholar

    [16]

    Kiess A, Pavičić D, Hänsch T W, Figger H 2008 Phys. Rev. A 77 053401Google Scholar

    [17]

    McKenna J, Sayler A M, Gaire B, Johnson N G, Zohrabi M, Carnes K D, Esry B D, Ben-Itzhak I 2009 J. Phys. B: At. Mol. Opt. Phys. 42 121003Google Scholar

    [18]

    McKenna J, Sayler A M, Gaire B, Johnson N G, Parke E, Carnes K D, Esry B D, Ben-Itzhak I 2009 Phys. Rev. A 80 023421Google Scholar

    [19]

    Liu Z T, Yuan K J, Shu C C, Hu W H, Cong S L 2010 J. Phys. B: At. Mol. Opt. Phys. 43 055601Google Scholar

    [20]

    He H X, Lu R F, Zhang P Y, Guo Y H, Han K L, He G Z 2011 Phys. Rev. A 84 033418Google Scholar

    [21]

    He H X, Lu R F, Zhang P Y, Han K L, He G Z 2012 J. Chem. Phys. 136 024311Google Scholar

    [22]

    Lu R F, Zhang P Y, Han K L 2008 Phys. Rev. E 77 066701Google Scholar

    [23]

    Hu J, Han K L, He G Z 2005 Phys. Rev. Lett. 95 123001

    [24]

    Feuerstein B, Thumm U 2003 Phys. Rev. A 67 043405Google Scholar

    [25]

    姚洪斌, 张季, 彭敏, 李文亮 2014 物理学报 63 198202Google Scholar

    Yao H B, Zhang J, Peng M, Li W L 2014 Acta Phys. Sin. 63 198202Google Scholar

    [26]

    Yao H B, Zhao G J 2014 J. Phys. Chem. A 118 9173Google Scholar

    [27]

    McKenna J, Sayler A M, Anis F, Gaire B, Johnson N G, Parke E, Hua J J, Mashiko H, Nakamura C M, Moon E, Chang Z, Carnes K D, Esry B D, Ben-Itzhak I 2008 Phys. Rev. Lett. 100 133001Google Scholar

  • [1] 何鑫, 李鑫焱, 李景辉, 张振华. Fe原子吸附的锑烯/WS2异质结的磁电子性质及调控效应. 物理学报, 2022, 0(0): . doi: 10.7498/aps.71.20220949
    [2] 赵嘉琳, 程开, 于雪克, 赵纪军, 苏艳. 几种典型含能材料光激发解离的含时密度泛函理论研究. 物理学报, 2021, 70(20): 203301. doi: 10.7498/aps.70.20210670
    [3] 马东飞, 侯文清, 徐春华, 赵春雨, 马建兵, 黄星榞, 贾棋, 马璐, 刘聪, 李明, 陆颖. 脂质体包裹荧光受体方法研究α-突触核蛋白在磷脂膜上的结构和动态特征. 物理学报, 2020, 69(3): 038701. doi: 10.7498/aps.69.20191607
    [4] 谢武, 沈斌, 张勇军, 郭春煜, 许嘉诚, 路欣, 袁辉球. 重费米子材料与物理. 物理学报, 2019, 68(17): 177101. doi: 10.7498/aps.68.20190801
    [5] 李文涛, 于文涛, 姚明海. 采用量子含时波包方法研究H/D+Li2LiH/LiD+Li反应. 物理学报, 2018, 67(10): 103401. doi: 10.7498/aps.67.20180324
    [6] 张斯淇, 陆景彬, 刘晓静, 刘继平, 李宏, 梁禺, 张晓茹, 刘晗, 吴向尧, 郭义庆. 运用理想光子禁带模型实现对激发态原子系统演化的调控. 物理学报, 2018, 67(9): 094205. doi: 10.7498/aps.67.20172050
    [7] 王文彬, 朱银燕, 殷立峰, 沈健. 复杂氧化物中电子相分离的量子调控. 物理学报, 2018, 67(22): 227502. doi: 10.7498/aps.67.20182007
    [8] 陈建宏, 郑小平, 张正荣, 吴学勇. 氢负离子在少周期激光场中解离时的干涉效应. 物理学报, 2016, 65(8): 083202. doi: 10.7498/aps.65.083202
    [9] 郑小丰, 樊群超, 孙卫国, 范志祥, 张燚, 付佳, 李博. 用差分收敛法研究NaLi分子部分电子态的完全振动能谱. 物理学报, 2015, 64(20): 203301. doi: 10.7498/aps.64.203301
    [10] 冯小静, 郭玮, 路兴强, 姚洪斌, 李月华. 三态K2分子飞秒含时光电子能谱的理论研究. 物理学报, 2015, 64(14): 143303. doi: 10.7498/aps.64.143303
    [11] 杨增强, 张力达. 红外激光载波包络相位对氦原子的极紫外光(XUV)吸收谱的量子调控研究. 物理学报, 2015, 64(13): 133203. doi: 10.7498/aps.64.133203
    [12] 姚洪斌, 李文亮, 张季, 彭敏. K2分子在强激光场下的量子调控:缀饰态选择性分布. 物理学报, 2014, 63(17): 178201. doi: 10.7498/aps.63.178201
    [13] 姚洪斌, 张季, 彭敏, 李文亮. H2+在强激光场中的解离及其量子调控的理论研究. 物理学报, 2014, 63(19): 198202. doi: 10.7498/aps.63.198202
    [14] 徐天宇, 何峰. H2+在阿秒以及双色飞秒激光脉冲中解离时电子位置的相干控制. 物理学报, 2013, 62(6): 068201. doi: 10.7498/aps.62.068201
    [15] 黄仙山, 刘海莲. 运用动态腔环境实现对原子自发辐射过程的调控. 物理学报, 2011, 60(3): 034205. doi: 10.7498/aps.60.034205
    [16] 黄仙山, 刘海莲, 羊亚平, 石云龙. 运用动态Lorentz库实现对激发态原子动力学特性的调控. 物理学报, 2011, 60(2): 024205. doi: 10.7498/aps.60.024205
    [17] 吴 勇, 刘 玲, 王建国. O3+与H2碰撞中非解离电荷转移过程的全量子计算. 物理学报, 2008, 57(2): 947-956. doi: 10.7498/aps.57.947
    [18] 郑敦胜, 郭锡坤. 高激发振动态氰化氢分子的解离研究. 物理学报, 2004, 53(10): 3347-3352. doi: 10.7498/aps.53.3347
    [19] 胡正发, 王振亚, 孔祥蕾, 张先燚, 李海洋, 周士康, 王娟, 武国华, 盛六四, 张允武. 甲胺分子的同步辐射光电离解离质谱. 物理学报, 2002, 51(2): 235-239. doi: 10.7498/aps.51.235
    [20] 向天翔. 碰撞诱导离解(CID)(Ⅰ)——含时量子理论. 物理学报, 1990, 39(5): 718-725. doi: 10.7498/aps.39.718
计量
  • 文章访问数:  3509
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-21
  • 修回日期:  2019-06-03
  • 上网日期:  2019-09-01
  • 刊出日期:  2019-09-05

HD+分子的强场光解离动力学及其量子调控的理论研究

  • 新疆工程学院, 新能源材料研究重点实验室, 乌鲁木齐 830091
  • 通信作者: 李文亮, wenliangli@vip.126.com
    基金项目: 新疆维吾尔自治区高等学校科研计划项目(批准号: XJEDU2016I051, XJEDU2017M042, XJEDU2018I020)、国家自然科学基金(批准号: 11764041, 51462034)、新疆维吾尔自治区天山青年计划(批准号: 2017Q033, 2017Q034)和新疆维吾尔自治区自然科学基金(批准号: 2019D01A32)资助的课题.

摘要: 利用精确求解原子核与电子耦合运动的三维含时量子波包法, 理论研究了HD+分子在强激光场中的光解离动力学, 并给出了量子调控HD+分子光解离通道的理论方案. 通过分析HD+分子在不同的初始振动态和激光场强度下的光解离动力学过程及其解离核动能谱, 得出了HD+分子的光解离机理及其随激光场强度的变化规律. 研究结果表明, 利用激光场的强度可以实现HD+分子光解离通道的量子调控. 当激光场强度I1 = 4.0 × 1013 W/cm2时, HD+分子的光解离主要是通过净单光子吸收解离和净双光子吸收解离; 当激光场强度增大到I2 = 2.0 × 1014 W/cm2时, 直接双光子吸收解离取代了净单光子吸收解离, 净双光子吸收解离的比重也下降了.

English Abstract

    • 随着激光脉冲技术的飞速发展, 利用超快强激光脉冲场探索分子的超快动力学特性[15], 并对其实施高效量子调控已是当前强场动力学领域的重点研究课题之一[68]. 在众多的分子体系中, 氢及其同位素分子离子(H2+, D2+和HD+)通常是研究分子强场光解离动力学特性的范例, 这不仅是因为它具有便于理论研究的简单构型—两个原子核和一个电子, 更重要的是它在强激光场的作用下会呈现出许多重要的非线性光学现象. 譬如, 键软化解离[9]、键硬化解离[10]、阈上解离[11]、阈下解离[12]、电荷共振增强电离[13]等. 这些非线性光学现象所蕴含的动力学机理在物理化学、量子光学、量子信息、光化学等诸多领域都有着巨大的应用前景.

      同核双原子分子H2+和D2+与强激光场的相互作用, 由于受到对称性的限制, 在两个电子态(1sσ和2pσ态)之间的跃迁过程中所吸收或放出的光子能量必须为奇数个光子能量[14]. 异核双原子分子HD+由于存在固有跃迁偶极矩, 其光吸收过程将不再受到对称性的限制, 在两个电子态(1sσ和2pσ态)之间的跃迁过程中可以吸收或放出奇数和偶数个光子[15]. 譬如, 处在1sσ态上的H2+分子跃迁到2pσ态的过程中只能吸收一个光子和三个光子的能量, 而HD+分子在跃迁过程中不仅可以吸收一个光子和三个光子的能量, 还可以吸收两个光子的能量.

      然而, 早期的相关研究对固有跃迁偶极矩是否对HD+分子的光解离动力学产生影响, 以及HD+分子在光解离过程中是否存在双光子吸收都存在着较大的争议[1619]. 譬如, Kiess等[16]首次利用波长790 nm的激光脉冲场研究了HD+分子的光解离动力学, 并在HD+分子的光解离碎片中探测到双光子吸收解离峰. McKenna等[17]利用三维动量高分辨成像技术研究了HD+分子在波长790 nm的激光脉冲场作用下的光解离动力学, 其在激光场强度为2.0 × 1014 W/cm2时并没有观测到双光子吸收解离所对应的光解离碎片. 他们认为Kiess等[16]所探测的直接双光子吸收解离峰有可能是HD+分子在高振动态(ν > 6)上的单光子吸收而引发的键软化解离. 随后, McKenna等[18]实验研究了HD+分子在波长395 nm和790 nm的激光脉冲场中的光解离和电离动力学, 其认为HD+分子在波长790 nm、强度为3.0 × 1013 W/cm2的激光脉冲场中有可能发生净单光子吸收解离(分子在1sσ态上先吸收两个光子跃迁到2pσ态上, 随后再释放一个光子并发生解离).

      本文拟选取HD+分子为研究对象, 利用精确求解原子核与电子耦合的三维含时量子波包法理论研究HD+分子在强激光场中的光解离动力学. 通过深入分析分子在不同的初始振动态和激光场强度下的解离核动能谱, 得出HD+分子的光解离机理及其随激光场强度的变化规律. 在此基础之上, 进一步提出利用激光场强度量子调控HD+分子光解离通道的理论方案.

    • HD+分子在强激光场作用下的含时薛定谔方程可以写成(原子单位制)

      ${\rm{i}}\frac{\partial }{{\partial t}}\psi ({{R}},{{r}},t) = [{H_0} + W(t)]\psi ({{R}},{{r}},t),$

      其中H0是分子在无外场时的哈密顿量, W(t)是分子与强激光场的相互作用项, R是核间距矢量, r是电子相对于原子核质心的位置矢量. 当激光场的偏振方向平行于分子轴时, HD+分子的原子核在偏振方向上做一维运动, 电子的运动在柱坐标系下可以简化为二维运动. 在柱坐标系下, HD+分子在强激光场中的运动可以用三维的空间坐标描述: (R, ρ, z)[20,21]. 其中, R是分子的核间距, (ρ, z)是电子的坐标位置.

      HD+分子在无外场时的三维哈密顿量可以写成:

      ${H_0} = {H_{\rm{n}}} + {H_{\rm{e}}},$

      ${H_{\rm{n}}} = {T_{\rm{n}}} + {V_{{\rm{nn}}}} = - \frac{1}{{2{\mu _{\rm{n}}}}}\frac{{{\partial ^2}}}{{\partial {R^2}}} + \frac{1}{R},$

      $\begin{split} {H_{\rm{e}}} = & {T_{\rm{e}}} + {V_{{\rm{en}}}}\\ = & - \frac{1}{{2{\mu _{\rm{e}}}}}\left( {\frac{{{\partial ^2}}}{{\partial {\rho ^2}}} + \frac{1}{\rho }\frac{\partial }{{\partial \rho }} + \frac{{{\partial ^2}}}{{\partial {z^2}}}} \right)\\ & - \frac{1}{{\sqrt {{\rho ^2} + {{(z - {z_{\rm{p}}})}^2}} }} - \frac{1}{{\sqrt {{\rho ^2} + {{(z + {z_{\rm{d}}})}^2}} }}, \end{split}$

      其中${\mu _{\rm{n}}} = \dfrac{{{m_{\rm{p}}}{m_{\rm{d}}}}}{{{m_{\rm{p}}} + {m_{\rm{d}}}}}$${\mu _{\rm{e}}} = \dfrac{{{m_{\rm{p}}} + {m_{\rm{d}}}}}{{{m_{\rm{p}}} + {m_{\rm{d}}} + 1}}$是原子核和电子的折合质量, ${z_{\rm{p}}} = \dfrac{{{m_{\rm{d}}}}}{{{m_{\rm{p}}} + {m_{\rm{d}}}}}R$${z_{\rm{d}}} = \dfrac{{{m_{\rm{p}}}}}{{{m_{\rm{p}}} + {m_{\rm{d}}}}}R$是质子和氘核的相对位置. ${m_{\rm{p}}}$${m_{\rm{d}}}$分别是质子和氘核的质量.

      HD+分子与强激光场的相互作用项W(t)可以写成(长度规范下)

      $W(t) = - E(t)\left( {\frac{{{m_{\rm{p}}} + {m_{\rm{d}}} + 2}}{{{m_{\rm{p}}} + {m_{\rm{d}}} + 1}}z - \frac{{{m_{\rm{d}}} - {m_{\rm{p}}}}}{{{m_{\rm{p}}} + {m_{\rm{d}}}}}R} \right),$

      其中$E(t) = {E_0}f(t)\cos (\omega t)$. E0是激光场的峰值场强, ω是激光场的角频率, f(t)是激光场的脉冲包络, 本文采用高斯型的脉冲形式:

      $f(t) = \exp \left[ { - {{\left( {\frac{t}{\tau }} \right)}^2}} \right],$

      其中$\tau $是激光场的脉冲宽度.

      为了求解上述含时薛定谔方程, 我们在ρ方向上使用Crank-Nicholson (CN)方法, 在Rz方向使用以sine-DVR为基的二阶分裂算符方法[22,23]. 为了提高计算效率, 我们根据原子核与电子运动的时间尺度差, 优化了时间传播算符, 其表示形式为

      ${{\rm{e}}^{ - {\rm{i}}H\Delta t}} \approx {{\rm{e}}^{ - {\rm{i}}{T_R}\Delta t/2}}{\{U_{\rm{e}}^{{\rm{SPO}}}(\Delta t/N)\}^N}{{\rm{e}}^{ - {\rm{i}}{T_R}\Delta t/2}},$

      其中,

      $U_{\rm{e}}^{{\rm{SPO}}}(\Delta t) = {{\rm{e}}^{ - {\rm{i}}{T_z}\Delta t/2}}{{\rm{e}}^{ - {\rm{i}}V\Delta t}}{{\rm{e}}^{ - {\rm{i}}{T_z}\Delta t/2}}$

      是电子的时间传播算符, TRTzRz方向上的动能传播算符, V是系统的势能, N是原子核与电子的时间步长的比值.

      本文利用虚时传播的方法获取HD+分子的初始波函数, 并结合流算符和虚拟探测的方法[24]计算HD+分子的解离核动能谱. 在数值计算中, R的取值范围0—25 a.u., z的取值范围–40—40 a.u., ρ的取值范围0—15 a.u.. R, zρ的空间格点数目分别取250, 400和30, 时间步长分别取1.00, 0.05和0.05 a.u..

    • 由于固有跃迁偶极矩的存在, HD+分子的光吸收过程不再受对称性的限制, 其在1sσ和2pσ电子态之间的跃迁过程中可以吸收或放出奇数和偶数个光子. 图1给出了HD+分子在1sσ和2pσ光子缀饰态势能曲线上的解离示意图.

      图  1  HD+分子在1sσ和2pσ光子缀饰态势能曲线上的解离示意图[15,25]

      Figure 1.  The schematic diagram for the dissociation of HD+ molecule at the dressed potentials of electronic 1sσ and 2pσ states.

      图1中, $\left| {1{{s\sigma}} - 0\omega } \right\rangle $, $\left| {1{{s\sigma}} - 1\omega } \right\rangle $, $\left| {1{{s\sigma}} - 2\omega } \right\rangle $, $\left| {1{{s\sigma}} - 3\omega } \right\rangle $态是1sσ电子态在零光子、单光子、双光子、三光子缀饰下的光子缀饰态; $\left| {2{\rm{p\sigma}} - 0\omega } \right\rangle $, $\left| {2{{p\sigma}} - 1\omega } \right\rangle $, $\left| {2{{p\sigma}} - 2\omega } \right\rangle $, $\left| {2{{p\sigma}} - 3\omega } \right\rangle $态是2pσ电子态在零光子、单光子、双光子、三光子缀饰下的光子缀饰态. $\left| {1{{s\sigma}} - 0\omega } \right\rangle $态与$\left| {2{{p\sigma}} - 1\omega } \right\rangle $, $\left| {2{{p\sigma}} - 2\omega } \right\rangle $, $\left| {2{{p\sigma}} - 3\omega } \right\rangle $态在HD+分子的初始振动态ν = 9, ν = 6, ν = 3附近发生了交叉. 在强激光场的作用下, HD+分子易在交叉区域附近吸收一个光子、两个光子、三个光子的能量. 若HD+分子吸收一个光子能量后直接解离, 此过程称之为直接单光子吸收解离; 若HD+分子吸收两个光子能量后直接解离, 此过程称之为直接双光子吸收解离; 若HD+分子先吸收两个光子能量, 再释放一个光子后发生解离, 此过程称之为净单光子吸收解离; 若HD+分子先吸收三个光子能量, 再释放一个光子后发生解离, 此过程称之为净双光子吸收解离.

      图1所示, HD+分子的光解离过程主要由以下四个解离通道组成: 初始振动态ν = 9附近的分子在核间距R = 4.7 a.u.附近吸收一个光子跃迁到$\left| {2{{p\sigma}} - 1\omega } \right\rangle $态, 在$\left| {2{{p\sigma}} - 1\omega } \right\rangle $态上的分子运动到较大的核间距后通过直接单光子吸收解离, 解离后的能量约为0.1 eV; 初始振动态ν = 3附近的分子在核间距R = 3.2 a.u.附近吸收三个光子跃迁到$\left| {2{{p\sigma}} - 3\omega } \right\rangle $态, 当$\left| {2{{p\sigma}} - 3\omega } \right\rangle $态上的分子运动到较大的核间距R = 4.7 a.u.附近后又释放出一个光子, 进而回到$\left| {1{{s\sigma}} - 2\omega } \right\rangle $态并通过净双子吸收解离, 解离后的能量约为1.1 eV; 初始振动态ν = 6附近的分子在核间距R = 3.8 a.u.附近吸收两个光子跃迁到$\left| {2{{p\sigma}} - 2\omega } \right\rangle $态上, 运动到较大的核间距区域后通过直接双光子吸收解离, 解离后的能量约为1.6 eV; 此外, 处于$\left| {2{{p\sigma}} - 2\omega } \right\rangle $态上的分子还有可能运动到核间距R = 4.7 a.u.附近后又释放出一个光子, 进而回到$\left| {1{{s\sigma}} - 1\omega } \right\rangle $态并通过净单光子吸收解离, 解离后的能量约为0.3 eV. 简而言之, 在初始振动态ν = 3附近的HD+分子主要是通过净双子吸收解离, 在初始振动态ν = 9附近的HD+分子主要是通过直接单光子吸收解离, 而在初始振动态ν = 6 附近的HD+分子有两种解离通道, 其不仅可以通过直接双光子吸收解离, 还可以通过净单光子吸收解离.

      对于HD+分子, 在光解离过程中是否存在直接双光子吸收解离和净单光子吸收解离, 以及两种解离通道的竞争关系一直是早期实验研究[1618]争论的焦点, 本文重点讨论HD+分子的直接双光子吸收解离和净单光子吸收解离, 以及其解离概率随激光场强度的变化规律. 对于初始振动态ν = 9附近的HD+分子, 其光解离主要是通过直接单光子吸收解离, 对应的解离机理是键软化解离. 根据我们早期的理论研究[26], 可以得出直接单光子吸收解离过程受激光场强度的影响很小, 本文不做详细讨论. 对于初始振动态ν = 3和ν = 6附近的HD+分子, 其光解离动力学过程与激光场强度密切相关. 接下来进一步深入探讨HD+分子在初始振动态ν = 3和ν = 6上的光解离动力学及其随激光场强度的变化规律.

      为了进一步研究HD+分子在初始振动态ν = 3和ν = 6上的光解离过程及其动力学机理, 我们计算了HD+分子的光解离核动能谱. 图2给出了激光场强度I1 = 4.0 × 1013 W/cm2时, HD+分子在初始振动态ν = 3和ν = 6上的光解离核动能谱(KER). 根据图2(a)所示, HD+分子在初始振动态ν = 3上的光解离核动能谱有两个峰: 主峰的能量约为1.1 eV, 对应的解离通道为净双光子吸收解离[27]; 次峰的能量约为1.6 eV, 对应的解离通道为直接双光子吸收解离. 净双光子吸收解离和直接双光子吸收解离过程的解离机理为阈上解离, 其解离概率随着激光场强度的增强而增大[26]. 根据图2(b)所示, HD+分子在初始振动态ν = 6上的光解离核动能谱只有一个明显的解离峰, 其能量约为0.3 eV, 对应的解离通道为净单光子吸收解离.

      图  2  激光场强度I1 = 4.0 × 1013 W/cm2时, HD+分子在初始振动态ν = 3 (a)和ν = 6 (b)上的光解离核动能谱 (激光场的波长λ = 790 nm, 脉冲宽度τ = 40 fs)

      Figure 2.  At the laser intensity of I1 = 4.0 × 1013 W/cm2, the kinetic energy resolved distributions of dissociation for the vibrational states ν = 3 (a) and ν = 6 (b) of HD+ molecule. The laser wavelength is 790 nm and the pulse duration is 40 fs, respectively.

      图3给出了激光场强度增大到I2 = 2.0 × 1014 W/cm2时, HD+分子在初始振动态ν = 3和ν = 6上的光解离核动能谱. 根据图3(a)所示, HD+分子在初始振动态ν = 3上的光解离核动能谱与激光场强度I1 = 4.0 × 1013 W/cm2时相比, 能量为1.1 eV的主峰和1.6 eV的次峰的高度虽然都升高了, 但是主峰和次峰的相对高度却发生了显著的变化. 因此, 我们可以得出HD+分子在光解离过程中, 两种解离通道的比重发生了明显的改变, 主要表现在: HD+分子的净双光子吸收解离减弱了, 而直接双光子吸收解离明显增强了. 这主要是因为随着激光场强度的增加, HD+分子的净双光子吸收解离和直接双光子吸收解离的概率虽然都增大了, 但是当激光场强度I2 = 2.0 × 1014 W/cm2时, 处于净双光子吸收解离通道内的分子易发生电荷共振增强电离[26], 从而减弱了净双光子吸收解离的比重, 使得分子更有利于通过直接双光子吸收解离. 根据图3(b)所示, HD+分子在初始振动态ν = 6上的光解离核动能谱与激光场强度I1 = 4.0 × 1013 W/cm2时的相比, 其解离峰的位置发生了显著的变化. 主要表现在: 净单光子吸收解离峰(能量约为0.3 eV)几乎消失了, 光解离通道转变为直接双光子吸收解离, 对应的解离能约为1.6 eV. 这主要是因为直接双光子吸收解离引发的阈上解离概率随激光场强度的增强而增大.

      图  3  激光场强度I2 = 2.0 × 1014 W/cm2时, HD+分子在初始振动态ν = 3 (a)和ν = 6 (b)上的光解离核动能谱 (激光场的波长λ = 790 nm, 脉冲宽度τ = 40 fs)

      Figure 3.  At the laser intensity of I2 = 2.0 × 1014 W/cm2, the kinetic energy resolved distributions of dissociation for the vibrational states ν = 3 (a) and ν = 6 (b) of HD+ molecule. The laser wavelength is 790 nm and the pulse duration is 40 fs, respectively.

      综上所述, HD+分子的光解离通道与激光场强度息息相关. 当激光场强度I1 = 4.0 × 1013 W/cm2时, HD+分子在初始振动态ν = 6上的光解离主要是通过净单光子吸收解离, 如图4所示(蓝线); 当激光场强度增加到I2 = 2.0 × 1014 W/cm2时, 其光解离主要是通过直接双光子吸收解离, 如图4所示(红线). 因此, 我们可以得出: 调节激光场强度可以量子调控HD+分子的光解离通道. 当选用较弱的激光场强度时, 净单光子吸收解离是HD+分子的主要光解离通道; 当选取较强的激光场强度时, 净单光子吸收解离消失了, 直接双光子吸收解离成为HD+分子的主要光解离通道. 此外, 改变激光场强度还可以调控HD+分子的净双光子吸收解离和直接双光子吸收解离通道的分支比. 当选取较弱的激光场强度时, HD+分子的净双光子吸收解离起主导作用; 当选取较强的激光场强度时, HD+分子的直接双光子吸收解离将起主导作用.

      图  4  当激光场强度I1 = 4.0 × 1013 W/cm2 (蓝线)和I2 = 2.0 × 1014 W/cm2 (红线)时, HD+分子在初始振动态ν = 6上光解离通道(激光场的波长λ = 790 nm, 脉冲宽度τ = 40 fs)

      Figure 4.  The related pathways of dissociation for the vibrational state ν = 6 of HD+ molecule at the laser intensities of 4.0 × 1013 W/cm2 (blue line) and 2.0 × 1014 W/cm2 (red line). The laser wavelength is 790 nm, the pulse duration is 40 fs.

    • 本文利用精确求解原子核与电子耦合运动的三维含时量子波包法, 系统地研究了HD+分子在强激光场中的解离动力学及其量子调控的理论方案. HD+分子的光解离动力学过程与分子的初始振动态分布和激光场强度息息相关. 当激光场强度I1 = 4.0 × 1013 W/cm2时, HD+分子在初始振动态ν = 3上的光解离主要是通过净双光子吸收解离, 在初始振动态ν = 6上的光解离主要是通过净单光子吸收解离. 当激光场强度增加到I2 = 2.0 × 1014 W/cm2时, HD+分子在初始振动态ν = 3上的光解离除了通过净双光子吸收解离外, 还可以通过直接双光子吸收解离. 初始振动态ν = 6上的光解离却从净单光子吸收解离转变为直接双光子吸收解离.

      调节激光场强度不仅可以控制分子的净单光子吸收解离和直接双光子吸收解离, 还可以调节分子的净双光子吸收解离和直接双光子吸收解离的分支比. 这种量子调控方案在操控复杂分子的反应动力学方面有着巨大的应用前景, 譬如: 利用强激光场调控化学反应体系的光子缀饰态势能面, 使其沿着某个给定的反应通道进行化学反应, 进而实现了对产物通道的高效操控. 在后续的研究中, 我们将叠加两束超快强激光脉冲场(基准脉冲场及其二次谐波), 通过改变两束激光脉冲场的延迟时间和相对强度等参数量子调控分子的光解离动力学.

参考文献 (27)

目录

    /

    返回文章
    返回