搜索

x

留言板

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

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

对称分子${\text{H}}_{\text{2}}^{\text{ + }}$在强短波激光场中高次谐波椭偏率性质的研究

于术娟 刘竹琴 李雁鹏

引用本文:
Citation:

对称分子${\text{H}}_{\text{2}}^{\text{ + }}$在强短波激光场中高次谐波椭偏率性质的研究

于术娟, 刘竹琴, 李雁鹏

Ellipticity properties of symmetric molecules $ {\text{H}}_{\text{2}}^ + $ in strong and short-wavelength laser fields

Yu Shu-Juan, Liu Zhu-Qin, Li Yan-Peng
PDF
HTML
导出引用
  • 通过数值计算, 研究了强短波($ 400—600{\text{ nm}} $)激光场中${\text{H}}_2^+$分子高次谐波辐射的椭偏率性质. 研究表明, ${\text{H}}_2^+$分子在不同激光强度、不同激光波长、不同核间距及不同取向角下高次谐波的椭偏率性质是不同的; 特别是在两中心干涉区, 激发态在高次谐波产生中起着重要作用, 但在不同取向角下, 激发态对谐波椭偏率的影响不同; 分析表明, 这些不同的影响源于沿平行和垂直于激光偏振方向辐射的高次谐波的相对产量, 以及激发态对平行和垂直谐波产量的影响; 此外, 椭偏率的测量检验了强场近似和平面波近似在强短波激光场中是成立的, 并对强短波激光场中分子动力学做了更充分的研究.
    Ellipticity properties of high-order harmonic generation (HHG) from symmetric molecules $ {\text{H}}_{\text{2}}^ + $ in strong and short wavelength (less than 800nm) laser fields are numerically investigated. In this study, the ellipticity of harmonic is compared with the corresponding harmonic spectrum and dipole, and the calculation results are analyzed and the results obtained at different laser intensities, different laser wavelengths, different internuclear distances and different orientation angles are compared with each other. Our numerical simulations show that the influences of laser intensity, laser wavelength, internuclear distance and orientation angle on the ellipticity of harmonic are different. Especially in a two-center interference region, the excited state plays an important role in the HHG, but the effects of the excited state on the ellipticity of harmonic are different at different orientation angles. Further analysis shows that these different effects are due to the influence of the excited state on the harmonic yield. Using the numerical scheme, it is determined that in the two-center interference region, the excited state plays an important role in the parallel harmonic spectrum, while the effects of the excited state on the perpendicular harmonics at different angles are all very small, which results in different phase differences between the accurate harmonic spectrum and the harmonic spectrum only returning to the ground state. Overall, the relative yields of the accurate perpendicular harmonics are lower (higher) than those of the accurate parallel harmonics, but the intensities of the perpendicular harmonics, which only return to the ground state, are comparable to (or farther away from) those of the parallel harmonics which are only to return to ground state in the two-center interference regions. Therefore, the small (large) intensity ratio between the accurate perpendicular harmonic and accurate parallel harmonic can be attributed to the contributions of the excited state to harmonics. Then we can conclude that the harmonic spectra that only go back to the ground state show high (small) ellipticity, whereas the accurate harmonic spectra show small (high) ellipticity, resulting in a strong angle dependence of the influence of the excited state on the ellipticity of harmonic. In addition, in the high-order harmonic plateau region, the relative yields of harmonics can be well predicted by the corresponding dipoles, indicating the applicability of tunneling pictures and plane wave approximation in the strong and short-wave laser fields. When the ellipticity of harmonic occurs in the interference region due to the two-center characteristics of the symmetric potential, the results show that the polarization measurement can also be used to detect the structures of symmetric molecules and track the dynamic behaviors of excited states.
      通信作者: 李雁鹏, liyanpeng1107@163.com
    • 基金项目: 国家自然科学基金青年科学基金(批准号: 11904310)和延安市专项基金(批准号: 2019-03)资助的课题.
      Corresponding author: Li Yan-Peng, liyanpeng1107@163.com
    • Funds: Project supported by the Young Scientists Fund of National Natural Science Foundation of China (Grant No. 11904310), and the Special Foundation of Yan'an City, China (Grant No. 2019-03).
    [1]

    L’Huillier A, Schafer K J, Kulander K C 1991 J. Phys. B 24 3315Google Scholar

    [2]

    Itatani J, Levesque J, Zeidler Hiromichi Niikura D, Pepin H, Kieffer J C, Corkum P B, Villeneuve D M 2004 Nature 432 867Google Scholar

    [3]

    Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163Google Scholar

    [4]

    Le V H, Le A T, Xie R H, Lin C D 2007 Phys. Rev. A 76 013414Google Scholar

    [5]

    Vozzi C, Negro M, Calegari F, Sansone G, Nisoli M, Silvestri S De, Stagira S 2011 Nat. Phys. 7 822Google Scholar

    [6]

    Smirnova O, Mairesse Y, Patchkovskii S, Dudovich N, Villeneuve D, Corkum P, Ivanov M Yu 2009 Nature 460 972Google Scholar

    [7]

    Shafifir D, Soifer H, Bruner B D, Dagan M, Mairesse Y, Patchkovskii S, Ivanov M Yu, Smirnova O, Dudovich N 2012 Nature 485 343Google Scholar

    [8]

    Corkum P B 1993 Phys. Rev. Lett. 71 1994Google Scholar

    [9]

    Lewenstein M, Balcou Ph, Ivanov M Yu, L'Huillier Anne, Corkum P B 1994 Phys. Rev. A 49 2117Google Scholar

    [10]

    Keldysh L V 1965 Sov. Phys. JETP 20 1307

    [11]

    Kim V V, Ganeev R A, Boltaev G S, Alnaser A S 2020 J. Phys. B 53 155405Google Scholar

    [12]

    Dong F L, Tian Y Q, Yu S J, Wang S, Yang S P, Chen Y J 2015 Opt. Express 23 18106Google Scholar

    [13]

    Zhou X, Lock R, Wagner N, Li W, Kapteyn H C, Murnane M M 2009 Phys. Rev. Lett. 102 073902Google Scholar

    [14]

    Levesque J, Mairesse Y, Dudovich N, epin H P, Kieffer J C, Corkum P B, Villeneuve D M 2007 Phys. Rev. Lett. 99 243001Google Scholar

    [15]

    Ramakrishna S, Sherratt P A J, Dutoi A D, Seideman T 2010 Phys. Rev. A 81 021802(RGoogle Scholar

    [16]

    Strelkov V V, Gonoskov A A, Gonoskov I A, Ryabikin M Yu 2011 Phys. Rev. Lett. 107 043902Google Scholar

    [17]

    Akagi H, Otobe T, Staudte A, Shiner A, Turner F, Dörner R, Villeneuve D M, Corkum P B 2009 Science 325 1364Google Scholar

    [18]

    Bian X B, Bandrauk A D 2010 Phys. Rev. Lett. 105 093903Google Scholar

    [19]

    Chen Y J, Zhang B 2011 Phys. Rev. A 84 053402Google Scholar

    [20]

    Qin M, Zhu X, Zhang Q, Hong W, Lu P 2011 Opt. Express 19 25084Google Scholar

    [21]

    Du H, Luo L, Wang X, Hu B 2012 Phys. Rev. A 86 013846Google Scholar

    [22]

    Yu S J, Zhang B, Li Y P, Yang S P, Chen Y J 2014 Phys. Rev. A 90 053844Google Scholar

    [23]

    Barreau L, Veyrinas K, Gruson V, Weber S J, Auguste T, Hergott J F, Lepetit F, Carré B, Houver J -C, Dowek D, Salières P 2018 Nat. Commun. 9 4727Google Scholar

    [24]

    Sun F J, Chen C, Li W Y, Liu X, Li W, Chen Y J 2021 Phys. Rev. A 103 053108Google Scholar

    [25]

    Habibović D, Milošević D B 2020 Photonics 7 110Google Scholar

    [26]

    Huo X X, Xing Y H, Qi T, Sun Y, Li Bo, Zhang J, Liu X S 2021 Phys. Rev. A 103 053116Google Scholar

    [27]

    Chen Y J 2011 Phys. Rev. A 84 043423Google Scholar

    [28]

    Chen Y J, Zhang B 2012 Phys. Rev. A 86 023415Google Scholar

    [29]

    Bian X B, Bandrauk A D 2012 Phys. Rev. A 86 053417Google Scholar

    [30]

    Miao X Y, Zhang C P 2014 Phys. Rev. A 89 033410Google Scholar

    [31]

    Han Y C 2013 Phys. Rev. A 87 043404Google Scholar

    [32]

    Feit M D, Jr Fleck J A, Steiger A 1982 J. Comput. Phys. 47 412Google Scholar

    [33]

    Shiner A D, Trallero-Herrero C, Kajumba N, Bandulet H C, Comtois D, Légaré F, Giguère M, Kieffffer J C, Corkum P B, Villeneuve D M 2009 Phys. Rev. Lett. 103 073902Google Scholar

    [34]

    Yu S J, Li W Y, Li Y P, Chen Y J 2017 Phys. Rev. A 96 013432Google Scholar

    [35]

    Chen Y J, Hu B 2009 J. Chem. Phys. 131 244109Google Scholar

    [36]

    Chen Y J, Chen J, Liu J 2006 Phys. Rev. A 74 063405Google Scholar

  • 图 1  ${\text{H}}_2^ + $分子在激光强度$ I = 4 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $及波长$ \lambda = 600{\text{nm}} $时谐波的椭偏率、谐波谱及偶极子的比较 (a), (d), (g) 核间距为$ R = 1.7{\text{ a}}{\text{.u}}{\text{. }} $; (b), (e), (h) 核间距为$ R = 2{\text{ a}}{\text{.u}}{\text{. }} $; (c), (f), (i) 核间距为$ R = 2.3{\text{ a}}{\text{.u}}{\text{. }} $

    Fig. 1.  Comparison of the harmonic ellipticity, harmonic spectra and corresponding dipole for ${\text{H}}_2^ + $ molecules at the laser intensity $ I = 4 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $and wavelength $ \lambda = 600{\text{ nm}} $: (a), (d), (g) $ R = 1.7{\text{ a}}. {\text{u}}. $; (b), (e), (h) $ R = 2{\text{ a}}. {\text{u}}. $; (c), (f), (i) $ R = 2.3{\text{ a}}. {\text{u}}. $.

    图 2  与图1相同, 但激光强度$I = 5 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$和波长$\lambda = 400\;{\text{nm}}$

    Fig. 2.  Same as Fig.1, but the laser intensity $ I = 5 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $ and wavelength $ \lambda = 400{\text{ nm}} $.

    图 3  与图1相同, 但激光强度$I = 5 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$和波长$\lambda = 500\;{\text{nm}}$

    Fig. 3.  Same as Fig.1, but the laser intensity $I = 5 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$ and wavelength $\lambda = 500\;{\text{nm}}$.

    图 4  ${\text{H}}_2^ + $分子$I = 4 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$, $ \lambda = 600{\text{ nm}} $R=2 a.u.时谐波谱、偶极子和谐波椭偏率的比较 (a), (d), (g) $ \theta = {10^\circ} $; (b), (e), (h) $ \theta = {30^\circ} $; (c), (f), (i) $ \theta = {50^\circ} $

    Fig. 4.  Comparison of the harmonic ellipticity, harmonic spectra and corresponding dipole for ${\text{H}}_2^ + $ molecules with R=2 a.u. at the laser intensity $ I = 4 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $and wavelength $ \lambda = 600{\text{ nm}} $: (a), (d), (g) $ \theta = {10^\circ} $; (b), (e), (h) $ \theta = {30^\circ} $; (c), (f), (i) $ \theta = {50^\circ} $.

    图 5  ${\text{H}}_2^ + $分子$I = 4 \times {10^{14}}\;{\text{ W}}/{\text{c}}{{\text{m}}^2}$, $\lambda = 600\;{\text{nm}}$, R=2 a.u.及$ \theta = {10^\circ} $时谐波谱、相位差和谐波椭偏率的比较 (a) 精确平行与跃迁回基态的平行谐波谱(基态平行谐波)比较; (b) 精确平行与精确垂直谐波谱的比较; (c) 精确与基态垂直谐波谱的比较; (d) 基态平行与基态垂直谐波谱的比较; (e) 精确与基态谐波相位差比较; (f) 精确与基态谐波椭偏率比较

    Fig. 5.  The comparison of the harmonic spectra, the phase differences $ {\delta _1}\left( d \right) $ and the ellipticity of harmonic for ${\text{H}}_2^ + $ molecules with R=2 a.u. and $ \theta = {10^\circ} $ at the laser intensity $I = 4 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$and wavelength $\lambda = 600\;{\text{ nm}}$: (a) The comparison of accurate parallel harmonic and the parallel harmonic that transition back to ground state (ground-state parallel harmonic); (b) the comparison of accurate parallel and perpendicular harmonic; (c) comparison between accurate and ground state perpendicular harmonic spectrum; (d) comparison between ground-state parallel and ground-state perpendicular harmonic; (e) the comparison of the phase difference of accurate and ground-state harmonic; (f) comparison between accurate and ground-state harmonic ellipticity.

    图 6  与图5相同, 但$ \theta = {30^\circ} $

    Fig. 6.  Same as Fig. 5, but $ \theta = {30^\circ} $.

    图 7  与图5相同, 但$ \theta = {50^\circ} $

    Fig. 7.  Same as Fig.5, but $ \theta = {50^\circ} $.

  • [1]

    L’Huillier A, Schafer K J, Kulander K C 1991 J. Phys. B 24 3315Google Scholar

    [2]

    Itatani J, Levesque J, Zeidler Hiromichi Niikura D, Pepin H, Kieffer J C, Corkum P B, Villeneuve D M 2004 Nature 432 867Google Scholar

    [3]

    Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163Google Scholar

    [4]

    Le V H, Le A T, Xie R H, Lin C D 2007 Phys. Rev. A 76 013414Google Scholar

    [5]

    Vozzi C, Negro M, Calegari F, Sansone G, Nisoli M, Silvestri S De, Stagira S 2011 Nat. Phys. 7 822Google Scholar

    [6]

    Smirnova O, Mairesse Y, Patchkovskii S, Dudovich N, Villeneuve D, Corkum P, Ivanov M Yu 2009 Nature 460 972Google Scholar

    [7]

    Shafifir D, Soifer H, Bruner B D, Dagan M, Mairesse Y, Patchkovskii S, Ivanov M Yu, Smirnova O, Dudovich N 2012 Nature 485 343Google Scholar

    [8]

    Corkum P B 1993 Phys. Rev. Lett. 71 1994Google Scholar

    [9]

    Lewenstein M, Balcou Ph, Ivanov M Yu, L'Huillier Anne, Corkum P B 1994 Phys. Rev. A 49 2117Google Scholar

    [10]

    Keldysh L V 1965 Sov. Phys. JETP 20 1307

    [11]

    Kim V V, Ganeev R A, Boltaev G S, Alnaser A S 2020 J. Phys. B 53 155405Google Scholar

    [12]

    Dong F L, Tian Y Q, Yu S J, Wang S, Yang S P, Chen Y J 2015 Opt. Express 23 18106Google Scholar

    [13]

    Zhou X, Lock R, Wagner N, Li W, Kapteyn H C, Murnane M M 2009 Phys. Rev. Lett. 102 073902Google Scholar

    [14]

    Levesque J, Mairesse Y, Dudovich N, epin H P, Kieffer J C, Corkum P B, Villeneuve D M 2007 Phys. Rev. Lett. 99 243001Google Scholar

    [15]

    Ramakrishna S, Sherratt P A J, Dutoi A D, Seideman T 2010 Phys. Rev. A 81 021802(RGoogle Scholar

    [16]

    Strelkov V V, Gonoskov A A, Gonoskov I A, Ryabikin M Yu 2011 Phys. Rev. Lett. 107 043902Google Scholar

    [17]

    Akagi H, Otobe T, Staudte A, Shiner A, Turner F, Dörner R, Villeneuve D M, Corkum P B 2009 Science 325 1364Google Scholar

    [18]

    Bian X B, Bandrauk A D 2010 Phys. Rev. Lett. 105 093903Google Scholar

    [19]

    Chen Y J, Zhang B 2011 Phys. Rev. A 84 053402Google Scholar

    [20]

    Qin M, Zhu X, Zhang Q, Hong W, Lu P 2011 Opt. Express 19 25084Google Scholar

    [21]

    Du H, Luo L, Wang X, Hu B 2012 Phys. Rev. A 86 013846Google Scholar

    [22]

    Yu S J, Zhang B, Li Y P, Yang S P, Chen Y J 2014 Phys. Rev. A 90 053844Google Scholar

    [23]

    Barreau L, Veyrinas K, Gruson V, Weber S J, Auguste T, Hergott J F, Lepetit F, Carré B, Houver J -C, Dowek D, Salières P 2018 Nat. Commun. 9 4727Google Scholar

    [24]

    Sun F J, Chen C, Li W Y, Liu X, Li W, Chen Y J 2021 Phys. Rev. A 103 053108Google Scholar

    [25]

    Habibović D, Milošević D B 2020 Photonics 7 110Google Scholar

    [26]

    Huo X X, Xing Y H, Qi T, Sun Y, Li Bo, Zhang J, Liu X S 2021 Phys. Rev. A 103 053116Google Scholar

    [27]

    Chen Y J 2011 Phys. Rev. A 84 043423Google Scholar

    [28]

    Chen Y J, Zhang B 2012 Phys. Rev. A 86 023415Google Scholar

    [29]

    Bian X B, Bandrauk A D 2012 Phys. Rev. A 86 053417Google Scholar

    [30]

    Miao X Y, Zhang C P 2014 Phys. Rev. A 89 033410Google Scholar

    [31]

    Han Y C 2013 Phys. Rev. A 87 043404Google Scholar

    [32]

    Feit M D, Jr Fleck J A, Steiger A 1982 J. Comput. Phys. 47 412Google Scholar

    [33]

    Shiner A D, Trallero-Herrero C, Kajumba N, Bandulet H C, Comtois D, Légaré F, Giguère M, Kieffffer J C, Corkum P B, Villeneuve D M 2009 Phys. Rev. Lett. 103 073902Google Scholar

    [34]

    Yu S J, Li W Y, Li Y P, Chen Y J 2017 Phys. Rev. A 96 013432Google Scholar

    [35]

    Chen Y J, Hu B 2009 J. Chem. Phys. 131 244109Google Scholar

    [36]

    Chen Y J, Chen J, Liu J 2006 Phys. Rev. A 74 063405Google Scholar

  • [1] 汉琳, 苗淑莉, 李鹏程. 优化组合激光场驱动原子产生高次谐波及单个超短阿秒脉冲理论研究. 物理学报, 2022, 71(23): 233204. doi: 10.7498/aps.71.20221298
    [2] 魏博宁, 焦志宏, 周效信. 非对称波形激光驱动的氢原子高次谐波频移及控制. 物理学报, 2022, 71(7): 073201. doi: 10.7498/aps.71.20212146
    [3] 徐新荣, 仲丛林, 张铱, 刘峰, 王少义, 谭放, 张玉雪, 周维民, 乔宾. 强激光等离子体相互作用驱动高次谐波与阿秒辐射研究进展. 物理学报, 2021, 70(8): 084206. doi: 10.7498/aps.70.20210339
    [4] 蔡怀鹏, 高健, 李博原, 刘峰, 陈黎明, 远晓辉, 陈民, 盛政明, 张杰. 相对论圆偏振激光与固体靶作用产生高次谐波. 物理学报, 2018, 67(21): 214205. doi: 10.7498/aps.67.20181574
    [5] 罗香怡, 刘海凤, 贲帅, 刘学深. 非均匀激光场中氢分子离子高次谐波的增强. 物理学报, 2016, 65(12): 123201. doi: 10.7498/aps.65.123201
    [6] 管仲, 李伟, 王国利, 周效信. 激光驱动晶体发射高次谐波的特性研究. 物理学报, 2016, 65(6): 063201. doi: 10.7498/aps.65.063201
    [7] 李雁鹏, 于术娟, 陈彦军. 不同取向角下CO2分子波长依赖的垂直谐波效率. 物理学报, 2015, 64(18): 183102. doi: 10.7498/aps.64.183102
    [8] 汪礼锋, 贺新奎, 滕浩, 运晨霞, 张伟, 魏志义. 5fs驱动激光脉冲的高次谐波选择性优化. 物理学报, 2014, 63(22): 224103. doi: 10.7498/aps.63.224103
    [9] 李小刚, 李芳, 何志聪. 双色场驱动下高次谐波的径向量子轨道干涉. 物理学报, 2013, 62(8): 087201. doi: 10.7498/aps.62.087201
    [10] 卢发铭, 夏元钦, 张盛, 陈德应. 飞秒强激光脉冲驱动Ne高次谐波蓝移产生相干可调谐极紫外光实验研究. 物理学报, 2013, 62(2): 024212. doi: 10.7498/aps.62.024212
    [11] 葛德彪, 魏兵. 互易定理计算分层半空间上方任意取向偶极子的远区场. 物理学报, 2012, 61(5): 050301. doi: 10.7498/aps.61.050301
    [12] 曹卫军, 成春芝, 周效信. 原子在双色组合场中产生高次谐波的转换效率与激光波长的关系. 物理学报, 2011, 60(5): 054210. doi: 10.7498/aps.60.054210
    [13] 崔磊, 王小娟, 王帆, 曾祥华. 脉冲激光偏振方向对氧分子高次谐波的影响——基于含时密度泛函理论的模拟. 物理学报, 2010, 59(1): 317-321. doi: 10.7498/aps.59.317
    [14] 李会山, 李鹏程, 周效信. 强激光场中模型氢原子的势函数对产生高次谐波强度的影响. 物理学报, 2009, 58(11): 7633-7639. doi: 10.7498/aps.58.7633
    [15] 张春丽, 祁月盈, 刘学深, 丁培柱. 双色激光场中高次谐波转化效率的提高. 物理学报, 2007, 56(2): 774-780. doi: 10.7498/aps.56.774
    [16] 赵松峰, 周效信, 金 成. 强激光场中模型氢原子和真实氢原子的高次谐波与电离特性研究. 物理学报, 2006, 55(8): 4078-4085. doi: 10.7498/aps.55.4078
    [17] 顾 斌, 崔 磊, 曾祥华, 张丰收. 超强飞秒激光脉冲作用下氢分子的高次谐波行为——基于含时密度泛函理论的模拟. 物理学报, 2006, 55(6): 2972-2976. doi: 10.7498/aps.55.2972
    [18] 崔 磊, 顾 斌, 滕玉永, 胡永金, 赵 江, 曾祥华. 脉冲激光偏振方向对氮分子高次谐波的影响--基于含时密度泛函理论的模拟. 物理学报, 2006, 55(9): 4691-4694. doi: 10.7498/aps.55.4691
    [19] 张秋菊, 盛政明, 张 杰. 超短脉冲强激光与固体靶作用产生的高次谐波红移. 物理学报, 2004, 53(7): 2180-2183. doi: 10.7498/aps.53.2180
    [20] 李 昆, 张 杰, 余 玮. 激光与固体靶作用产生高次谐波的振荡镜面模型. 物理学报, 2003, 52(6): 1412-1417. doi: 10.7498/aps.52.1412
计量
  • 文章访问数:  641
  • PDF下载量:  35
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-10-12
  • 修回日期:  2022-11-12
  • 上网日期:  2022-12-02
  • 刊出日期:  2023-02-20

对称分子${\text{H}}_{\text{2}}^{\text{ + }}$在强短波激光场中高次谐波椭偏率性质的研究

  • 延安大学物理与电子信息学院, 延安 716000
  • 通信作者: 李雁鹏, liyanpeng1107@163.com
    基金项目: 国家自然科学基金青年科学基金(批准号: 11904310)和延安市专项基金(批准号: 2019-03)资助的课题.

摘要: 通过数值计算, 研究了强短波($ 400—600{\text{ nm}} $)激光场中${\text{H}}_2^+$分子高次谐波辐射的椭偏率性质. 研究表明, ${\text{H}}_2^+$分子在不同激光强度、不同激光波长、不同核间距及不同取向角下高次谐波的椭偏率性质是不同的; 特别是在两中心干涉区, 激发态在高次谐波产生中起着重要作用, 但在不同取向角下, 激发态对谐波椭偏率的影响不同; 分析表明, 这些不同的影响源于沿平行和垂直于激光偏振方向辐射的高次谐波的相对产量, 以及激发态对平行和垂直谐波产量的影响; 此外, 椭偏率的测量检验了强场近似和平面波近似在强短波激光场中是成立的, 并对强短波激光场中分子动力学做了更充分的研究.

English Abstract

    • 当原子或分子暴露于强激光场中时, 库仑势被电场严重扭曲. 因此, 电子可通过隧穿从原子核中脱离, 脱离的电子在激光场中进行加速运动, 当激光场改变方向时, 电子再次返回原子核, 进而被原子核俘获, 并伴随着高能光子的释放, 这个过程称为高次谐波产生. 由于高次谐波可用作阿秒脉冲光源, 近年来, 学者们越来越关注高次谐波辐射的研究[13], 并提供了一种独特的方式来探测物质在超快时间尺度上的结构[4,5]和电子动力学[6,7].

      高次谐波产生的物理过程可通过三步模型[8]和列文斯坦模型[9]来理解, 模型中有两个重要的假设, 即强场近似. 首先, 只有基态对系统的演化起重要作用; 第二, 当电子脱离原子核时, 忽略库伦势的影响. 这些假设适用于具有Keldysh[10]参数$ \gamma = \sqrt {{I_{\text{p}}}/\left( {2{U_{\text{p}}}} \right)} \ll 1 $强的低频激光场, 其中, ${U_{\text{p}}} = {E^2}/(4\omega _0^2)$是有质动能, ${I_{\text{p}}}$是原子的电离能, $ E $$ {\omega _0} $分别是激光的振幅和频率, 因此, 在低频激光场中隧穿图像和平面波近似都是成立的.

      以往学者们主要关注沿激光偏振方向辐射的高次谐波功率谱, 最近, 准直对称分子在低频激光场中谐波的椭偏率问题成为学者们关注的焦点, 例如$ {\text{H}}_{\text{2}}^ + $[11,12]$ {{\text{N}}_{\text{2}}} $[13-15]; 根据实验和理论研究结果, 在线偏振激光场作用下, 准直分子可以产生谐波的椭圆偏振态, 而谐波椭偏率的起源是多电子效应[13]、库伦效应[15]、两中心干涉效应[11]以及电子运动的量子不确定性[16]. 随后研究了不对称分子谐波的偏振性质[1722], 理论研究表明, 由于最高占据分子轨道的不对称性, 在$ {\text{CO}} $的高次谐波产生中可观察到大的椭偏率[20,21]; 由于不同高次谐波产生路径间的相互作用, 在$ {\text{He}}{{\text{H}}^{{\text{2}} + }} $的谐波产生中也可观察到大的椭偏率[22]; 同时圆[23]、椭圆[24]和正交[25,26]偏振激光脉冲驱动高次谐波产生的椭偏率一直也是实验和理论上的研究热点. 以前主要研究强的低频激光场中谐波的椭偏率问题, 但在强的高频激光场中(Keldysh参数$ \gamma $接近或略大于1), 谐波的椭偏率未被研究过, 由于高次谐波的椭偏率不仅与平行谐波及垂直于激光偏振方向发射的谐波(即垂直谐波)密切相关, 还与高次谐波的相位密切相关, 为研究高频激光场中分子高次谐波的产生提供了深入的见解; 另外, 高频激光场中电子的动力学非常复杂, 激发态效应非常强[27], 激发态对谐波发射的影响已被研究过[28-30]; 且激发态在平衡核间距及大核间距附近对低频激光场中高次谐波谱上干涉最小的影响也被研究过[31]; 目前的研究表明, 当Keldysh参数$ \gamma \approx 1 $时, 隧穿图像仍然适用[11], 但在Keldysh参数$ \gamma \gt 1 $的情况下, 隧穿图像是否成立, 平面波近似是否成立, 激发态对谐波椭偏率有什么影响? 为了回答上述问题, 研究高频场中谐波的椭偏率问题成为必然.

      通过数值求解二维含时Schrödinger方程, 研究了高频激光场中对称分子${\text{H}}_{\text{2}}^{\text{ + }}$高次谐波的偏振特性; 结果表明: ${\text{H}}_{\text{2}}^{\text{ + }}$分子谐波的偏振特性与激光强度、激光波长、核间距及取向角有关; 利用数值方案, 确定了在两中心干涉区, 激发态在对称分子辐射的平行谐波谱中起重要作用, 对于不同角度的垂直谐波, 激发态的影响非常小, 激发态对平行谐波和垂直谐波的不同影响导致精确垂直和精确平行谐波谱的强度比与仅跃迁回基态的垂直和平行谐波谱的强度比不同, 且精确谐波谱与仅跃迁回基态的谐波谱的相位差不同, 进而导致激发态对谐波椭偏率的影响具有强烈的角依赖性. 此外, 偏振的目的是对高频激光场中强场近似的检验, 尤其是对平面波近似的检验.

    • 原子单位下, 二维模型$~ {\text{H}}_2^ + ~$分子的哈密顿为$H(t) = {{{{\boldsymbol{p}}^2}}}/{2} + V({\boldsymbol{r}}) + {\boldsymbol{r}} \cdot {\boldsymbol{E}}(t)$. 式中$V({\boldsymbol{r}}) = -Z/\sqrt {\xi + {\boldsymbol{r}}_1^2} - Z/\sqrt {\xi + {\boldsymbol{r}}_2^2}$为软心势, $ {{\boldsymbol{r}}_1} = {\boldsymbol{r}} - {{\boldsymbol{R}}_1} $, ${{\boldsymbol{r}}_2} = {\boldsymbol{r}} - {{\boldsymbol{R}}_2}$, 其中$ {{\boldsymbol{R}}_1} $$ {{\boldsymbol{R}}_2} $分别是两个核的坐标$ (R/2, 0) $$ (-R/2, 0) $. $ Z $是有效电荷, 为避免库仑奇异, 软心势参数$\xi = 0.5$. 激光的电场分量${\boldsymbol{E}}(t) = {\boldsymbol{e}}f(t)E \sin{\omega _0}t$, $ f(t) $为包络函数, $ E $$ {\omega _0} $分别为激光场的振幅和频率, $ {\boldsymbol{e}} $是激光电场方向的单位向量. ${\text{H}}_2^ + $分子的电离能$ {{I} _{\text{p}}} = 1.11{\text{ }}{\kern 1 pt} {\text{a}}{\text{.u}}{\text{. }} $. 数值计算中, 采用10个周期的激光脉冲, 其中前三个周期线性上升, 后三个周期线性下降, 剩下四个周期的激光强度是常数. 含时薛定谔方程采用谱方法求解[32], 空间格子的大小$ {L_x} \times {L_y} = 204.8{\text{ a}}. {\text{u}}. \times 204.8{\kern 1 pt} {\text{ a}}. {\text{u}}. $, 空间步长为0.4 a.u.. 为消除边界处对波包的反射作用, 采用滤波函数$ F(x) = {\text{co}}{{\text{s}}^{1/8}}\left[ {{{{\text{π }}(\left| x \right| - {x_0})} \mathord{\left/ {\vphantom {{{\text{π }}(\left| x \right| - {x_0})} {({L_x} - 2{x_0})}}} \right. } {({L_x} - 2{x_0})}}} \right] $对波函数从边界至${x_0} = {{{L_x}} \mathord{\left/ {\vphantom {{{L_x}} 8}} \right. } 8}$处进行滤波, 对y轴的处理与x轴相似.

      靶分子被激光照射时产生的精确平行谐波和精确垂直谐波产量[33]的计算式为

      $ \begin{gathered} {F_{\parallel ( \bot )}}(\omega ) = \int {\left\langle {\psi (t)} \right|} {{\boldsymbol{e}}_{\parallel ( \bot )}} \cdot \nabla V\left| {\psi (t)} \right\rangle {{\text{e}}^{{\text{i}}\omega t}}{\text{d}}t \\ \\ \end{gathered} \text{, } $

      其中$ \psi (t) $为与哈密顿$ H(t) $对应的含时波函数. 基于强场近似理论, 与基态再结合的电子对谐波辐射起主要贡献. 因此有

      $ F_{\parallel ( \bot )}^0(\omega ) = \int {\left\langle 0 \right|{{\boldsymbol{e}}_{\parallel ( \bot )}} \cdot \nabla V\left| {\psi (t)} \right\rangle {a^ * }_0{{\text{e}}^{{\text{i}}\omega t}}{\text{d}}t} . $

      (2)式表示电子从连续态返回至基态辐射的高次谐波, 其中$ {a_0}(t) = \left\langle {0} \mathrel{\left | {\vphantom {0 {\psi (t)}}} \right. } {{\psi (t)}} \right\rangle $, 式中未考虑激发态对高次谐波产生的贡献. 此外, 垂直谐波和平行谐波的振幅率及相位差决定了谐波的椭偏率[12]:

      $ \varepsilon {\text{ = }}\sqrt {\frac{{1{\text{ + }}{\mu ^2} - \sqrt {1 + 2{\mu ^2}\cos (2\delta ) + {\mu ^4}} }}{{1{\text{ + }}{\mu ^2} + \sqrt {1 + 2{\mu ^2}\cos (2\delta ) + {\mu ^4}} }}} \text{, } $

      式中, $ \mu = {S_ \bot }/{S_{||}} $为垂直与平行谐波的强度比, $ \delta = {\phi _ \bot } - {\phi _{||}} $是相位差, 其中$ {S_{||( \bot )}} = {\left| {{F_{||( \bot )}}(\omega )} \right|^2} $, ${\phi _{||( \bot )}}(\omega ) = {\text{arg}}\left[ {{F_{||( \bot )}}(\omega )} \right]$, $ \left| \delta \right|{\text{ = }}k{\text{π }} + {\delta _1} $, 这里$ k = 0, 1 $$ \left|{\delta }_{1}\right|\in \left[0, \text{π}\right] $, 椭偏率范围为$ 0 \leqslant \varepsilon \leqslant 1 $, 线偏振、椭圆偏振及圆偏振分别对应于$ \varepsilon {\text{ = }}0 $, $0 \lt \varepsilon \lt 1$$ \varepsilon {\text{ = 1}} $. 由(3)式可知, 当平行和垂直高次谐波的相位差$ {\delta _1} $$ \pi /2 $, 且垂直谐波强度接近平行谐波强度时, 显示高的椭偏率[12].

    • 基于强场近似, 对于原子和对称分子, 由于势的对称性, 只辐射奇次谐波, 其产生路径可简单地表示为$ \left\langle {\boldsymbol{p}} \right|{{\boldsymbol{e}}_x} \cdot {\boldsymbol{r}}{\left| {\text{0}} \right\rangle _{\text{i}}}\left\langle {\text{0}} \right|{\boldsymbol{e}} \cdot {\boldsymbol{r}}{\left| {\boldsymbol{p}} \right\rangle _{\text{r}}} $, 意味着电子从相同的初始态$ \left| 0 \right\rangle $电离并与之复合, 其中初始态$ \left| 0 \right\rangle $具有确定的宇称; 在高次谐波产生中, 偶极子矩阵元$ \left\langle {\boldsymbol{p}} \right|{{\boldsymbol{e}}_x} \cdot {\boldsymbol{r}}{\left| {\text{0}} \right\rangle _{\text{i}}} $为与电子沿激光极化方向$ {{\boldsymbol{e}}_x} $的电离有关, $ \left\langle {\text{0}} \right|{\boldsymbol{e}} \cdot {\boldsymbol{r}}{\left| {\boldsymbol{p}} \right\rangle _{\text{r}}} $与电子沿$ e $方向的复合有关; 此外, 在电离过程中, 激光场导致库仑势明显畸变; 因此, 电离过程与分子结构无关, 对于谐波的产生, 只有返回过程与分子结构紧密相关的.

      对于对称分子${\text{H}}_2^ + $, ${D_{//}}\left( {\omega , \theta } \right) = \left\langle 0 \right|{{\boldsymbol{e}}_x} \cdot {\boldsymbol{r}}{\left| {\boldsymbol{p}} \right\rangle _{\text{r}}}$是连续态$ \left| {p} \right\rangle $和基态$\left| 0 \right\rangle $之间的偶极矩阵元, 假设对称分子${\text{H}}_2^ + $的基态波函数为${\psi _0}\left( {{\boldsymbol{r}}, \theta } \right){\text{ = }}\left\langle {{\boldsymbol{r}}} \mathrel{\left | {\vphantom {{\boldsymbol{r}} 0}} \right. } {0} \right\rangle = {N_{\text{f}}}\left( {{{\text{e}}^{ - \kappa {r_{\text{a}}}}} + {{\text{e}}^{ - \kappa {r_{\text{b}}}}}} \right)$, $ {r_{\text{a}}} = \left| {r - {R_1}} \right| $, $ {r_{\text{b}}} = \left| {r - {R_2}} \right| $$ \kappa {\text{ = }}\sqrt {2{I_{\text{p}}}} $, $ {N_{\text{f}}} $是归一化系数且$ {I_{\text{p}}} $是基态$ \left| 0 \right\rangle $的电离能. 连续态$ \left| {\boldsymbol{p}} \right\rangle $由平面波函数$\left| {{{\text{e}}^{{\text{i}}{{\boldsymbol{p}}_k} \cdot {\boldsymbol{r}}}}} \right\rangle $近似得到, 动量$ {p_k} $为连续态$ \left| {\boldsymbol{p}} \right\rangle $的有效动量, 其中${p_k}{\text{ = }}\left| {{{\boldsymbol{p}}_k}} \right| = {\left[ {2\left( {{I_{\text{p}}} + {E_{\text{p}}}} \right)} \right]^{1/2}}$考虑了库仑效应[34], $ {E_{\text{p}}}{\text{ = }}\omega - {I_{\text{p}}} $为连续态$ \left| {\boldsymbol{p}} \right\rangle $的能量, 对于${\text{H}}_2^ + $分子, 沿激光极化方向的偶极矩阵元计算公式为[35]

      $ {D_{{\text{||}}}}(\omega ,\theta ) \propto G(\omega ,\theta ){D_a} \text{, } $

      式中, $ G(\omega , \theta ) = 2{\cos}[{p_k}R/2\cos \theta ] $, ${D_a} = \displaystyle\int {\text{d}}r[{{\text{e}}^{ - \kappa r}}{{\boldsymbol{e}}_{//}} \cdot {\boldsymbol{r}}{\sin}({{\boldsymbol{p}}_k} \cdot {\boldsymbol{r}})]$为原子的偶极矩阵元.

      基于强场近似, 谐波辐射谱由分子结构决定[2]:

      $ {S_{//} }(\omega ,\theta ) \sim {\left| {a(\omega ){D_{//} }(\omega ,\theta )} \right|^2} \text{, } $

      式中, $a(\omega )$是连续态电子隧穿时的谱振幅, 不随分子取向角而变化, 且不依赖于连续态$ \left| {\boldsymbol{p}} \right\rangle $的能量, 因此可得沿激光偏振方向的高次谐波谱${S_{//} }(\omega , \theta )$与对应偶极矩${D_{{//}}}(\omega , \theta )$之间的定性关系, 进而可知隧穿图像成立.

    • 为研究高频激光场中${\text{H}}_2^ + $分子高次谐波偏振特性受核间距及激光参数的影响, 在以下讨论中, 主要关注三种核间距情况$ R = 1.7{\text{ a}}. {\text{u}}. $, $ R = 2{\text{ a}}. {\text{u}}. $$ R = 2.3{\text{ a}}. {\text{u}}. $, 为便于比较, 不同核间距的${\text{H}}_2^ + $分子具有相同的电离能$ {I_{\text{p}}} = 1.1{\text{ a}}{\text{.u}}{\text{. }} $, 另外, 为研究${\text{H}}_2^ + $分子谐波椭偏率对取向的依赖性, 以下讨论主要关注三个典型角度$\theta = {10^ \circ }, {\text{ }}{30^ \circ }, {\text{ }}{50^ \circ }$.

    • 图1(a)图1(b)图1(c)分别表示由(3)式计算的${\text{H}}_2^ + $分子在核间距$ R = 1.7{\text{ a}}. {\text{u}}. $, $ R = 2{\text{ a}}. {\text{u}}. $$ R = 2.3{\text{ a}}. {\text{u}}. $及不同取向角下的谐波椭偏率曲线, 随着核间距的变化, 谐波最大椭偏率所处的谐波阶的位置会改变; 图1中第二、三行图是与谐波椭偏率曲线对应的角依赖的谐波谱和偶极子曲线, 其与谐波椭偏率曲线趋势相反, 在不同取向角下, 图1(d)图1(e)图1(f)中谐波相对产量与图1(g)图1(h)图1(i)中相应偶极子曲线符合的很好, 如(5)式所述, (5)式是基于隧穿图像得到的, 该现象表示隧穿图像在高频激光场中成立, 同时(5)式中偶极子的连续态$ \left| {\boldsymbol{p}} \right\rangle $由平面波函数$\left| {{{\text{e}}^{{\text{i}}{{\boldsymbol{p}}_k} \cdot {\boldsymbol{r}}}}} \right\rangle $近似得到, 进而可得在高频激光场中平面波近似也成立. 但是谐波谱和偶极子曲线之间存在显著的差异, 图1(g)图1(h)图1(i)中的偶极子曲线展示了一个清晰的最小值, 这个最小值起源于两中心干涉效应, 但在图1(d)图1(e)图1(f)中对应的谐波谱曲线上未观察到该最小值, 文献[27]中详细讨论了高频激光场中对称分子谐波谱上由两中心干涉效应而导致的最小值消失的原因,即在短波长(高频激光场)的情况下, 由于激发态对分子高次谐波的重要贡献, 导致干涉引起的凹槽在高次谐波谱中消失. 通过对谐波椭偏率曲线及相应偶极子曲线的深入比较发现偶极子曲线上的最小值对应最大椭偏率所处的谐波阶位置, 暗示了通过谐波椭偏率测量来探测高频激光场中对称分子偶极子上的最小值是可行的, 因此在高频激光场中, 两中心干涉效应是谐波椭偏率的起源. 因为在$\theta = {50^{\circ}}$下偶极子曲线上未出现最小值, 因此谐波在$\theta = {50^{\circ}}$基本不极化, 即使在某些区域有极化, 也是因为两中心干涉效应, 谐波在$\theta = {50^{\circ}}$不极化这一现象说明在没有干涉效应的情况下(5)式依然是成立的, 即在没有干涉效应的情况下, 隧穿图像和平面波近似在高频激光场中都成立.

      图  1  ${\text{H}}_2^ + $分子在激光强度$ I = 4 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $及波长$ \lambda = 600{\text{nm}} $时谐波的椭偏率、谐波谱及偶极子的比较 (a), (d), (g) 核间距为$ R = 1.7{\text{ a}}{\text{.u}}{\text{. }} $; (b), (e), (h) 核间距为$ R = 2{\text{ a}}{\text{.u}}{\text{. }} $; (c), (f), (i) 核间距为$ R = 2.3{\text{ a}}{\text{.u}}{\text{. }} $

      Figure 1.  Comparison of the harmonic ellipticity, harmonic spectra and corresponding dipole for ${\text{H}}_2^ + $ molecules at the laser intensity $ I = 4 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $and wavelength $ \lambda = 600{\text{ nm}} $: (a), (d), (g) $ R = 1.7{\text{ a}}. {\text{u}}. $; (b), (e), (h) $ R = 2{\text{ a}}. {\text{u}}. $; (c), (f), (i) $ R = 2.3{\text{ a}}. {\text{u}}. $.

      图1类似, 图2图3给出了激光强度$ I = 5 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $、波长$ \lambda $分别为400和500 nm的谐波椭偏率曲线、谐波谱曲线和偶极子曲线. 随着激光参数的变化, 由(4)式计算的偶极子曲线上的最小值移向其他谐波阶, 如图2图3第三行所示. 同时图2图3中的谐波相对产量与相应偶极子曲线依旧一一对应,如(5)式所述, 上述现象说明随着激光参数的改变, 隧穿图像和平面波近似在高频激光场中成立, 另外, 由于图2(d)图2(e)图2(f)图3(d)图3(e)图3(f)中谐波谱曲线上的最小值未出现, 谐波的最大椭偏率仍旧对应偶极子曲线最小值所在的谐波阶, 因此, 两中心干涉效应为谐波最大椭偏率的起源. 最重要的是随着激光参数和核间距的改变, 谐波在$\theta = 50^\circ$基本不极化, 即使在某些区域有极化, 也是因为两中心干涉效应, 如图2图3所示. 因此谐波在$\theta = {50^{\circ}}$不极化这一现象说明在没有干涉效应的情况下(5)式依然是成立的, 因此椭偏率测量至少在本文计算的Keldysh参数范围内($ 1.028 \leqslant \gamma \leqslant {\text{1}}{\text{.377}} $)给出了(5)式的适用性.

      图  2  与图1相同, 但激光强度$I = 5 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$和波长$\lambda = 400\;{\text{nm}}$

      Figure 2.  Same as Fig.1, but the laser intensity $ I = 5 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $ and wavelength $ \lambda = 400{\text{ nm}} $.

      图  3  与图1相同, 但激光强度$I = 5 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$和波长$\lambda = 500\;{\text{nm}}$

      Figure 3.  Same as Fig.1, but the laser intensity $I = 5 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$ and wavelength $\lambda = 500\;{\text{nm}}$.

    • 由于激发态对分子高次谐波的重要贡献, 导致干涉效应引起的凹槽在高次谐波谱中消失, 那么激发态是如何影响谐波椭偏率的? 因此在图4中, 研究了激发态对谐波椭偏率特性的影响, 将(1)式计算的结果(实线)与(2)式计算的结果(短虚线)进行比较, 其中(1)式包含所有束缚态对高次谐波的贡献(精确谐波), (2)式只考虑电子返回到基态对高次谐波的贡献(基态谐波); 对于图4(a)$ \theta = {10^\circ} $的小角及图4(b)$ \theta = {30^\circ} $的中间角情况, 只跃迁至基态的谐波谱显示一个大的抑制区域, 其在精确谐波谱中消失. 另外, 图4(a)图4(b)中基态谐波(短虚线)上的抑制区域对应于图4(g)图4(h)中偶极子曲线上的凹槽, 如虚线箭头所示, 该偶极子曲线由(4)式计算得到; 结果表明, 激发态在干涉诱导凹槽区影响高次谐波的产生; 然而, 在不同取向角下, 激发态对谐波椭偏率的影响是不同的, 如图4(d)(f)所示; 例如, 对于图4(d)$ \theta = {10^\circ} $的情况, 基态谐波的椭偏率(短虚线)显示了一个尖锐的峰值, 该曲线由(2)式代入(3)式计算得到, 而精确谐波椭偏率(实曲线)峰值显示了更广泛的分布, 该实曲线由(1)式代入(3)式计算得到; 图4(e)$ \theta = {30^\circ} $的中间角情况则相反, 由(2)式计算的跃迁回基态的谐波谱在此角度上的椭偏率很小(短虚线), 而由(1)式计算的精确谐波(实曲线)具有显著的椭偏率. 对于$ \theta = {50^\circ} $的大角情况, 在两中心干涉区, 由(2)式计算的基态谐波谱略偏离由(1)式计算的精确谐波谱, 如图4(c)所示; 然而, 当图4(f)中的实曲线显示小的椭偏率时, 图4(f)中的短虚线则显示大的椭偏率.

      图  4  ${\text{H}}_2^ + $分子$I = 4 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$, $ \lambda = 600{\text{ nm}} $R=2 a.u.时谐波谱、偶极子和谐波椭偏率的比较 (a), (d), (g) $ \theta = {10^\circ} $; (b), (e), (h) $ \theta = {30^\circ} $; (c), (f), (i) $ \theta = {50^\circ} $

      Figure 4.  Comparison of the harmonic ellipticity, harmonic spectra and corresponding dipole for ${\text{H}}_2^ + $ molecules with R=2 a.u. at the laser intensity $ I = 4 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $and wavelength $ \lambda = 600{\text{ nm}} $: (a), (d), (g) $ \theta = {10^\circ} $; (b), (e), (h) $ \theta = {30^\circ} $; (c), (f), (i) $ \theta = {50^\circ} $.

      下面, 探讨了激发态对谐波椭偏率影响的起源. 对于图5$ \theta = {10^\circ} $的情况, 其具有与图4相同的激光和分子参数; 在图5(a)中, 在由(2)式计算的基态谐波谱(短虚线)上可观察到大的抑制区域, 该区域在由(1)式计算的精确谐波谱(实曲线)中消失; 然而对于图5(c)中垂直谐波情况, 由(2)式计算的基态谐波谱(短虚线)与由(1)式计算的精确谐波谱(实曲线)非常相似; 上述现象表明, 激发态对平行和垂直谐波的影响不同, 导致精确谐波的相位差较大, 如图5(e)所示. 然而, 当分别比较图5(b)中由(1)式计算的平行和垂直谐波谱和图5(d)中由(2)式计算的平行和垂直谐波谱时, 发现在图5(b)的整个谐波平台区域, 精确垂直谐波的相对产量低于精确平行谐波的; 在此情况下, 在两中心干涉区, 由于消除了激发态对高次谐波的贡献, 由(2)式计算的基态垂直谐波的相对产量与基态平行谐波的相当, 如图5(d)所示; 因此, 在两中心干涉区, 由(1)式计算的精确垂直谐波和精确平行谐波之间小的强度比归因于激发态对谐波的贡献, 该影响导致图5(f)中精确谐波谱显示小的椭偏率. 换言之, 即使精确谐波的相位差大于基态谐波的相位差, 但由(1)式获得的垂直和平行谐波之间的强度比非常小, 导致图5(f)中精确谐波频谱显示出小的椭偏率, 在这种情况下, 激发态抑制了谐波的椭偏率.

      图  5  ${\text{H}}_2^ + $分子$I = 4 \times {10^{14}}\;{\text{ W}}/{\text{c}}{{\text{m}}^2}$, $\lambda = 600\;{\text{nm}}$, R=2 a.u.及$ \theta = {10^\circ} $时谐波谱、相位差和谐波椭偏率的比较 (a) 精确平行与跃迁回基态的平行谐波谱(基态平行谐波)比较; (b) 精确平行与精确垂直谐波谱的比较; (c) 精确与基态垂直谐波谱的比较; (d) 基态平行与基态垂直谐波谱的比较; (e) 精确与基态谐波相位差比较; (f) 精确与基态谐波椭偏率比较

      Figure 5.  The comparison of the harmonic spectra, the phase differences $ {\delta _1}\left( d \right) $ and the ellipticity of harmonic for ${\text{H}}_2^ + $ molecules with R=2 a.u. and $ \theta = {10^\circ} $ at the laser intensity $I = 4 \times {10^{14}}\;{\text{W}}/{\text{c}}{{\text{m}}^2}$and wavelength $\lambda = 600\;{\text{ nm}}$: (a) The comparison of accurate parallel harmonic and the parallel harmonic that transition back to ground state (ground-state parallel harmonic); (b) the comparison of accurate parallel and perpendicular harmonic; (c) comparison between accurate and ground state perpendicular harmonic spectrum; (d) comparison between ground-state parallel and ground-state perpendicular harmonic; (e) the comparison of the phase difference of accurate and ground-state harmonic; (f) comparison between accurate and ground-state harmonic ellipticity.

      对于图6$ \theta = {30^\circ} $的中间角, 在平台区域, 激发态对平行和垂直高次谐波的影响与图5类似, 即激发态对平行和垂直谐波的贡献不同, 导致精确谐波产生较大的相位差. 此外, 在图6(b)的两中心干涉区, 由(1)式计算的精确垂直谐波强度与精确平行谐波强度相当, 由于图5(b)中精确垂直谐波与精确平行谐波的强度比小, 因此$ \theta = {10^\circ} $图5(f)中精确谐波的椭偏率低于$ \theta = {30^\circ} $图6(f)中精确谐波的椭偏率, 这些数值分析揭示了为什么图1(a)$ \theta = {10^\circ} $的谐波椭偏率峰值低于图1(a)$ \theta = {30^\circ} $的谐波椭偏率峰值. 且由于激发态的缺失, 由(2)式计算的图6(a)中基态谐波谱(短虚线)出现了一个较宽的抑制区域, 进而使得该区域基态平行谐波的相对产量偏离基态垂直谐波的相对产量, 如图6(d)所示. 这些影响表明, 由(1)式计算的精确垂直谐波和精确平行谐波之间大的强度比归因于激发态对谐波的贡献, 因此由(1)式计算的精确谐波谱显示高的椭偏率, 如图6(f)所示; 激发态增强了谐波的椭偏率.

      图  6  与图5相同, 但$ \theta = {30^\circ} $

      Figure 6.  Same as Fig. 5, but $ \theta = {30^\circ} $.

      应当指出, 对于图7$ \theta = {50^\circ} $的情况, 激发态在两中心干涉区对平行谐波具有较小的影响, 如图7(a)所示; 在整个平台区域, 激发态不影响垂直谐波, 见图7(c), 该结果导致在两中心干涉区, 由(2)式计算的图7(d)中垂直和平行谐波之间的强度比大于由(1)式计算的图7(b)中垂直和平行谐波的强度比, 且使得精确谱的相位差$ {\delta _1} \sim 0.16{\text{π }} $; 相应的在两中心干涉区, 由(1)式计算的谐波谱具有低椭偏率, 因此, 激发态抑制了谐波的椭偏率.

      图  7  与图5相同, 但$ \theta = {50^\circ} $

      Figure 7.  Same as Fig.5, but $ \theta = {50^\circ} $.

      为使研究更全面, 在不同取向角及其他激光参数下, 例如取向角$ \theta = {20^\circ} $, $ \theta = {40^\circ} $, $ \theta = {60^\circ} $, 激光强度$ I = 5 \times {10^{14}}{\text{ W}}/{\text{c}}{{\text{m}}^2} $及波长$ \lambda = 400{\text{ nm}} $, $ \lambda = 500{\text{ nm}} $, 检验了谐波椭偏率的行为, 上述观察到的现象基本不变, 即在两中心干涉区, 激发态对平行谐波有重要影响, 但垂直谐波不受激发态影响, 进而导致精确谐波和基态谐波相位差不同; 且在不同取向角, 激发态对谐波椭偏率的影响是不同的. 总的来说, 虽然由(1)式计算的精确垂直谐波的相对产量低于(高于)精确平行谐波的, 但在两中心干涉区, 由于激发态的缺失, 由(2)式计算的基态垂直谐波的强度更接近(远离)基态平行谐波的强度. 相应的, 由(1)式计算的精确垂直谐波和精确平行谐波之间小(大)的强度比归因于激发态对谐波的贡献. 因此, 可得出结论, 而由(1)式代入(3)式计算的精确谐波椭偏率具有低(高)的椭偏率, 激发态抑制(增强)了谐波的椭偏率, 大角度谐波椭偏率的起源与图7(f)中的曲线相似.

      需要强调的是在上述讨论中, 我们没有考虑较大核间距的情况. 对于高频激光场的情况, 由于激光频率较高, 核间距较大分子的基态和激发态之间的耦合可能会对谐波辐射产生重要的影响. 此时只考虑基态贡献的强场近似模型不适用, 需要在通常的强场近似模型中引入激发态贡献[36]例如R =4 a.u., 此时基态和激发态强耦合, 其贡献可比. 关于高频场中核间距大的分子谐波椭偏率的研究, 将在未来的工作中做进一步的详细研究.

    • 本文通过数值模拟研究了在高频强激光场作用下对称分子${\text{H}}_2^ + $高次谐波的椭偏率特性. 在研究中, 利用谐波椭偏率与相应谐波谱和偶极子进行比对, 并在不同激光强度之间、不同激光波长之间、不同核间距之间及不同取向角之间进行比对、分析计算结果, 结果表明: 激光强度、激光波长、核间距及取向角对谐波椭偏率的影响不同, 特别是在分子两中心干涉区, 谐波椭偏率通常较强, 在该区域, 激发态对谐波有重要影响, 但激发态对谐波椭偏率的影响在不同角度是不同的; 并分析了其原因, 即当精确垂直谐波的相对产量低于(高于)精确平行谐波的相对产量, 但在两中心干涉区, 由于激发态的缺失, 跃迁回基态的垂直谐波的强度更接近(远离)基态平行谐波的强度; 相应地, 精确垂直谐波和精确平行谐波之间小(大)的强度比归因于激发态对谐波的贡献; 因此, 精确谐波具有低(高)的椭偏率, 激发态抑制(增强)了谐波的椭偏率. 此外, 在高次谐波平台区, 谐波的相对产量可通过相应偶极子很好地预测, 说明隧穿图像和平面波近似在高频激光场中的适用性. 由于对称势的两中心特性, 高次谐波在干涉区出现高的椭偏率, 结果表明, 椭偏率测量也可用于探测对称分子的结构及追踪其激发态的动力学.

参考文献 (36)

目录

    /

    返回文章
    返回