网站地图
下载:
-
强激光场与原子、分子相互作用可以产生许多非线性光学现象。高次谐波(high-order harmonic generation, HHG)作为其中一种非常重要的现象在探测原子、分子内电子运动以及产生孤立阿秒脉冲中起到非常重要的作用并被广泛研究[1-3]。
目前，高次谐波的辐射过程可由半经典的三步模型[4]来描述。例如：束缚电子首先由隧道电离进入连续态(第1步：电离过程)；随后电子在激光场的作用下加速并远离核子(第2步：加速过程)；最后电子在激光场驱动下反向加速与核子发生碰撞并辐射高次谐波(第3步：回碰过程)。对于原子体系来说，三步模型可以很好地解释谐波辐射过程。但对于分子体系，由于体系多出的核自由度，分子中电子的电离以及回碰与原子相比有很大的不同，进而导致一些分子谐波特有的现象。例如：多光子电离[5]、电荷共振增强电离[6]、离解态电离[7]在谐波辐射中的贡献；多通道谐波辐射[8]；谐波红移[9]；谐波蓝移[10]；谐波辐射的空间分布[11]等。
最近，随着远紫外(extreme ultraviolet, XUV)光源的发展，利用XUV光源调控谐波过程得到了很多关注。例如：作者[12]利用超短XUV光源驱动He原子增强谐波强度，并获得1个脉宽为43as的孤立阿秒脉冲。ZHANG等人[13]利用XUV光源与双色激光场的组合场驱动He+离子获得了一个47as脉冲。随后，这种利用超短XUV光源增强谐波辐射强度的方案被TAKAHASHI等人[14]实验证实。但是上述方法中都是用XUV光源驱动原子体系，并且都是用于控制三步模型中的电离过程，即利用XUV光源增大电离几率，进而增强谐波辐射强度。
因此，鉴于上述原因，本文中研究了H2+在红外激光与远紫外激光场驱动下谐波辐射的过程，并且通过经典与量子的分析给出了由于远紫外光源的加入所导致的多重谐波截止能量以及谐波频移的原因。若无说明，本文中采用原子单位(atomic units，a.u.)。
-
$ \begin{array}{l} {\rm{i}}\frac{{\partial \psi (z, R, t)}}{{\partial t}} = \left[ {} \right. - \frac{1}{{{m_{\rm{p}}}}}\frac{{{\partial ^2}}}{{\partial {R^2}}} - \frac{{2{m_{\rm{p}}} + 1}}{{4{m_{\rm{p}}}}}\frac{{{\partial ^2}}}{{\partial {z^2}}} + \\ V\left( {z, R} \right) + \left( {1 + \frac{1}{{2{m_{\rm{p}}} + 1}}} \right)zE\left( t \right)\left. {} \right]\psi (z, R, t) \end{array} $
(1) 式中，mp，R，z分别为核质量、核间距与电子坐标; t是时间; $ V\left( {z, R} \right) = 1/R - 1/\sqrt {{{\left( {z - R/2} \right)}^2} + 1} - 1/\sqrt {{{\left( {z + R/2} \right)}^2}} + 1$为H2+势能; $\psi \left( {z, R, t} \right) $为电子波函数。激光场形式为：
$ \begin{array}{l} E\left( t \right) = {E_1}{\rm{exp}}\left[ { - 4{\rm{ln}}\left( 2 \right){{\left( {t/{\tau _1}} \right)}^2}} \right]{\rm{cos}}({\omega _1}t) + \\ {E_{{\rm{XUV}}}}{\rm{exp}}\{ - 4{\rm{ln}}\left( 2 \right){\left[ {\left( {t - {t_{\rm{d}}}} \right)/{\tau _{{\rm{XUV}}}}} \right]^4}\} \times \\ {\rm{cos}}\left[ {{\omega _{{\rm{XUV}}}}\left( {t - {t_{\rm{d}}}} \right)} \right] \end{array} $
(2) 式中，E1，ω1，τ1为红外(infrared, IR)场振幅、频率和半峰全宽; EXUV，ωXUV，τXUV为XUV光源的振幅、频率和半峰全宽; td为两束激光场的延迟时间。
高次谐波频谱图可表示为：
$ S\left( \omega \right) = |\frac{1}{{\sqrt {2{\rm{ \mathsf{ π} }}} }}\smallint a\left( t \right){{\rm{e}}^{{\rm{ - i}}{\omega _1}t}}{\rm{d}}t|{^2} $
(3) 式中, $a\left( t \right) = - \left\langle {\psi \left( {z, R, t} \right)|\frac{{\partial V(z, R)}}{{\partial z}} + E\left( t \right)|\psi \left( {z, R, t} \right)} \right\rangle $为偶极加速度[21]，ω为谐波频次。由于谐波辐射能E=hω，h为普朗克常量。因此在本文中，通过分析谐波阶次ω/ω1来讨论截止能量。
-
图 1中给出H2+在单色10fs/800nm IR场以及该IR场与6fs/30nm XUV光源的组合场分别驱动下谐波辐射的特点。IR场场强为I1=4.0×1014W/cm2，XUV光源的场强为IXUV=1.0×1014W/cm2以及IXUV=1.0×1015W/cm2。两束激光场延迟时间为零。由图可知，在单色IR场驱动下，谐波光谱呈现经典的截止能量，即Ip+3.17Up(标记为“1”), 其中Ip是H2+的电离势, Up=I/4ω12为自由电子的有质动力势。当XUV光源的加入后，谐波光谱呈现2个明显特点：(1)第26阶谐波(H26)的强度被明显增强；(2)谐波截止能量明显延伸，形成能量间隔为XUV光子能量(ωXUV)的多重谐波截止能量。并且随着XUV光源强度增大，可以获得更多的截止能量, 但是其强度会逐渐减小。例如：当IXUV=1.0×1014W/cm2，谐波截止能量只延伸1个ωXUV，标记为“2”；当IXUV=1.0×1015W/cm2，谐波截止能量可以延伸2个ωXUV，标记为“3”。对于H26的增强现象很好理解，这是由于选用的XUV光源能量与H26能量相近形成了共振增强电离现象。对于多重谐波截止能量的延伸，作者认为是由于电子在回碰过程中吸收额外的XUV光子所形成的。但是由三步模型可知，谐波截止能量的延伸可以发生在电离-加速-回碰中的任何一步，因此，为了证明上述观点，本文中做了如下分析。
Figure 1. Harmonic spectra from H2+ driven by single IR field and the combined field (IR+XUV)
图 2a~图 2c中给出H2+在单色IR场以及上述组合场下的电离几率。T表示800nm激光场光学周期。由图可知，XUV光源的加入对H2+电离影响很小。因此，多重谐波截止能量延伸不可能发生在电离过程。图 3a~图 3c是电子波包随时间的演化图。从图中可知，在单色IR场下(见图 3a)，电子加速最远可以运动到35a.u.，这与经典结果xmax=E1/ω12是一致的。随着XUV光源的加入可见(见图 3b和图 3c)，电子加速的最远距离几乎保持不变，这说明多重谐波截止能量的延伸不是在加速过程中产生的。由此可见，多重谐波截止能量的延伸只能产生于回碰过程。
Figure 2. Ionization probabilities of H2+
Figure 3. Time-dependent electronic wave functions
为了进一步证明此观点以及理解分子谐波辐射过程，图 4中给出了上述激光场的波形。H2+在激光场驱动下辐射谐波的时频分析图[22]，谐波光谱随核间距R的变化关系以及核间距R随时间的演化图。在单色IR场下(见图 4a~图 4d)，由三步模型可知，电子电离可以发生在激光场振幅附近，即第1个~第7个电离时刻i1~7附近；随后电子在激光场中加速，并在激光场反向时与母核发生回碰辐射谐波，即第1个~第7个回碰时刻r1~7附近(见图 4a)。因此在谐波辐射时频分析图中可以观测到7个主要的谐波辐射过程，即第1个~第7个辐射能量峰P1~7(见图 4b)。为了更好地说明谐波辐射过程，基于牛顿方程的经典电子回碰能$ K\left( t \right) = \frac{1}{2}{[\int {} E(t){\rm{d}}t]^2}$也在图 4b中给出。可见，经典结果与量子结果符合得很好。分析图 4c和图 4d可知，谐波辐射主要发生在R=2.0a.u.到R=15.0a.u.区间。具体来说:(1)当t < -2.25T时(T表示800nm激光场周期)，由于激光强度较弱，H2+不能被充分电离，因此导致在较小核间距下的1sσg态和2pσu态之间的拉比共振[23]，这也是当R < 3.0a.u.时只有低阶谐波产生的原因; (2)随着激光在-2.25T < t < 0.75T时达到峰值，电离几率以及核运动逐渐增大，进而多光子共振电离在4.0a.u. < R < 6.0a.u.区间起主要作用，这是导致高阶谐波在4.0a.u. < R < 6.0a.u.产生的原因; (3)当0.75T < t < 2.25T时，虽然激光强度逐渐减小，但是由于随着R的增大，H2+的电离能逐渐减小，因此导致6.0a.u. < R < 10.0a.u.区间的电荷共振增强电离以及谐波辐射过程; (4)随着激光继续与H2+作用，虽然核间距可以延伸到更远的距离，例如：R>10.0a.u.，即离解态电离区间，但是由于激光强度迅速减小，因此只导致低阶谐波在此区间产生。对于组合场的情况(见图 4e~图 4h和图 4i~图 4l)，分析图 4e和图 4i可知，第2个~第6个电离时刻i2~6与回碰时刻r2~6被XUV光源覆盖，因此当加速电子返回母核回碰的过程中，由于XUV光源的作用，电子有几率在回碰过程中再吸收XUV光子，进而增大了回碰能，同时导致辐射出更高能量的谐波。为了证明上述分析过程，图 4f和图 4j中分别给出XUV光源强度为IXUV=1.0×1014W/cm2和IXUV=1.0×1015W/cm2时，组合场驱动H2+辐射谐波的时频分析图以及解经典牛顿方程得出的经典回碰能，其中，图 4f中两条曲线由上到下分别表示回碰能为K(t)以及K(t)+ωXUV；图 4j中3条曲线由上到下分别表示回碰能为K(t)，K(t)+ωXUV以及K(t)+2ωXUV。由图可知，当XUV光强较弱时(见图 4f)，电子在r2~6时刻回碰时会吸收1个XUV光子，导致P2~6的谐波截止能量以ωXUV的间隔延伸了1个XUV光子的能量。随着XUV光源强度增强(见图 4j)，电子在r2和r6时刻回碰时吸收XUV光子的几率增大，但依然只能吸收1个XUV光子，因此导致P2和P6的谐波截止能量只以ωXUV的间隔延伸了1个XUV光子的能量。但当电子在r3~5时刻回碰时，由于电离几率的增大以及较强的XUV光强，电子会吸收2个XUV光子，因此导致P3~5的谐波截止能量光谱呈现出Ip+3.17Up+2ωXUV的多重谐波截止能量延伸。但是由于每多吸收一个XUV光子的几率都会下降，因此导致多重谐波截止能量的强度呈阶梯形状的下降趋势。分析图 4g、图 4h和图 4k、图 4l可知，谐波辐射的贡献依然来自4个部分，即，Rabi振荡(R < 3.0a.u.)、多光子共振电离(4.0a.u. < R < 6.0a.u.)、电荷共振增强电离(6.0a.u. < R < 10.0a.u.)以及离解态电离(R>10.0a.u.)，并且多重谐波截止能量主要出现在多光子共振电离区间。
Figure 4. Time-dependent laser profile(a, e, i), harmonic order(b, f, j), R-dependent harmonic order (c, g, k) and time-dependent nuclear motion
图 5a中给出了改变XUV光源延迟时间后对多重谐波截止能量的影响。其中XUV光强为IXUV=1.0×1015W/cm2。计算结果表明，当延迟时间在-0.5T≤td≤0.5T时，谐波光谱可以呈现出多重谐波截止能量。但当td>0.5T时，谐波光谱的多重截止能量会逐渐消失。这里只给出td=±0.5T和td=±1.0T时的谐波谱图。图 5b中给出了XUV光源延迟时间对谐波频移的影响。由图可知，当td < 0或td>0时，谐波光谱呈现蓝移或者红移的现象(例如H11和H13)。
Figure 5. The effect of delay time
为了解释上述现象，图 6是上述激光场的波形图以及谐波辐射的时频分析图。由图可知，当延迟时间在td=-0.5T(见图 6c和图 6d)或td=0.5T时(见图 6e和图 6f)，电子回碰的r1~5时刻或r2~7时刻依然可以被XUV光源的主体覆盖，因此电子在回碰时依然可以吸收XUV光子，进而导致多重谐波截止能量的延伸。但当延迟时间在td=-1.0T(见图 6a和图 6b)或td=1.0T(见图 6g和图 6h)时，XUV光源的主体只覆盖激光包络的上升区间(r1，r2)或下降区间(r5~7)，因此电子在回碰过程r1，r2或r5~7时可以吸收XUV光子，进而导致谐波辐射过程P1，P2或P5~7的能量延伸。但是由于其延伸能量小于或只略大于经典的谐波截止能量(Ip+3.17Up)，因此在谐波光谱中多重谐波截止能量的现象消失了。由三步模型可知，在激光场的上升区间(dI(t)/dt>0)，由于激光强度的持续增强，后电离的电子会获得更多的能量，进而导致谐波光谱的蓝移；相反在激光场的下降区间(dI(t) /dt < 0)，由于激光强度的持续减小，后电离的电子将获得较低的能量，进而导致谐波光谱的红移[9-10]。分析谐波辐射的时频分析图可知，当td < 0时，激光上升区间的谐波辐射强度增强，因此导致谐波光谱的蓝移；当td>0时，激光下降区间的谐波辐射强度增强，因此导致谐波光谱的红移。
Figure 6. a, c, e, g—laser profiles of single IR field and XUV pulses under different delay time b, d, f, h—time-frequency analyses of different harmonics
-
理论研究了H2+在红外激光与远紫外激光场驱动下谐波辐射的特点。结果表明，谐波辐射的贡献分别来源于拉比振荡、多光子共振电离、电荷共振增强电离以及离解态电离。随着远紫外光的加入，谐波光谱呈现能量间隔为远紫外光子能量的多重谐波截止结构。适当调节远紫外光与红外激光场的延迟时间，谐波光谱呈现蓝移和红移的现象。
非常感谢中国科学院大连化学物理研究所韩克利研究员所提供的计算资源。
红外激光与远紫外激光场驱动H2+辐射谐波
H2+ radiation harmonics driven by infrared and extreme ultraviolet field
-
摘要: 为了了解H2+谐波辐射的过程，采用数值求解非玻恩-奥本海默近似薛定谔方程的方法，理论研究了H2+在10fs/800nm红外激光与6fs/30nm远紫外激光驱动下谐波辐射的特点。结果表明，谐波辐射的贡献主要来源于拉比振荡、多光子共振电离、电荷共振增强电离以及离解态电离；随着远紫外光的加入，谐波光谱呈现能量间隔为远紫外光子能量的多重谐波截止结构；当远紫外光与红外激光的延迟时间大于零或小于零时，谐波光谱呈现红移和蓝移的现象。该研究对理解分子谐波辐射过程是有帮助的。Abstract: In order to understand H2+ harmonic emission process, by using numerical solution of non-Bohn-Oppenheimer time-dependent Schrödinger equation, high-order harmonic generation (HHG) from H2+ driven by 10fs/800nm infrared (IR) field and 6fs/30nm extreme ultraviolet (XUV) pulse was theoretically investigated. The results show that, the contributions of HHG are mainly from Rabi-type oscillation, multi-photon resonance ionization, charge-resonance-enhanced ionization and dissociative ionization. With the introduction of XUV pulse, multi-harmonic cutoff extension separated by XUV photon energy can be found on HHG spectra. When the delay time between IR field and XUV pulse is larger than zero or smaller than zero, the red-shift and blue-shift of harmonics spectra can be obtained. The investigation is helpful to understand the molecular high-order harmonic generation.
-
Figure 1. Harmonic spectra from H2+ driven by single IR field and the combined field (IR+XUV)
Figure 2. Ionization probabilities of H2+
a—single IR field b—the combined field (IR+XUV) with IXUV=1.0×1014W/cm2 c—the combined field (IR+XUV) with IXUV=1.0×1015W/cm2
Figure 3. Time-dependent electronic wave functions
a—single IR field b—the combined field (IR+XUV) with IXUV=1.0×1014W/cm2 c—the combined field (IR+XUV) with IXUV=1.0×1015W/cm2
Figure 4. Time-dependent laser profile(a, e, i), harmonic order(b, f, j), R-dependent harmonic order (c, g, k) and time-dependent nuclear motion
a~d—single IR field e~h—the combined field (IR+XUV) with IXUV=1.0×1014W/cm2 i~l—the combined field (IR+XUV) with IXUV=1.0×1015 W/cm2
Figure 5. The effect of delay time
a—on multi-harmonic cutoff extension b—on frequency modulation of harmonics
-
[1] WANG Ch, LIU H L, TIAN J Sh, et al. Analysis of intrinsic atomic phase in process of extreme ultraviolet attosecond pulse generation[J]. Laser Technology, 2012, 36(3):342-345(in Chinese). [2] LIU H, FENG L Q. Mid-infrared field phase measurement and attosecond pulse generation[J]. Laser Technology, 2017, 41(2):151-158(in Chinese). [3] GOULIELMAKIS E, SCHULTZE M, HOFSTETTER M, et al. Single-cycle nonlinear optics[J]. Science, 2008, 320(5883):1614-1617. doi: 10.1126/science.1157846 [4] CORKUM P B. Plasma perspective on strong field multiphoton ionization[J]. Physical Review Letters, 1993, 71(13):1994-1997. doi: 10.1103/PhysRevLett.71.1994 [5] STAUDTE A, COCKE C L, PRIOR M H, et al. Observation of a nearly isotropic high-energy Coulomb explosion group in the fragmentation of D2 by short laser pulses[J]. Physical Review, 2002, A65(2):020703. [6] ZUO T, CHELKOWSKI S, BANDRAUK A D. Harmonic generation by the H2+ molecular ion in intense laser fields[J]. Physical Review, 1993, A48(5):3837-3844. [7] CHELKOWSKI S, FOISY C, BANDRAUK A D. Electron-nuclear dynamics of multiphoton H2+ dissociative ionization in intense laser fields[J]. Physical Review, 1998, A57(2):1176-1185. [8] BIAN X B, BANDRAUK A D. Multichannel molecular high-order harmonic generation from asymmetric diatomic molecules[J]. Physical Review Letters, 2010, 105(9):093903. doi: 10.1103/PhysRevLett.105.093903 [9] BIAN X B, BANDRAUK A D. Probing nuclear motion by frequency modulation of molecular high-order harmonic generation[J]. Physical Review Letters, 2014, 113(19):193901. doi: 10.1103/PhysRevLett.113.193901 [10] LI M Z, JIA G R, BIAN X B. Alignment dependent ultrafast electron-nuclear dynamics in molecular high-order harmonic generation[J]. The Journal of Chemical Physics, 2017, 146(8):084305. doi: 10.1063/1.4976973 [11] ZHANG J, GE X L, WANG T, et al. Spatial distribution on high-order-harmonic generation of an molecule in intense laser fields[J]. Physical Review, 2015, A92(1):013418. [12] FENG L Q, LIU H. Ultrashort ultraviolet source assisted enhancement of the attosecond pulse intensity[J]. Journal of Atomic and Molecular Physics, 2015, 32(4):616-620(in Chinese). [13] ZHANG G T, LIU X S. Generation of an extreme ultraviolet supercontinuum and isolated sub-50as pulse in a two-colour laser field[J]. Journal of Physics, 2009, B42(12):125603. [14] TAKAHASHI E J, KANAI T, ISHIKAWA K L, et al. Dramatic enhancement of high-order harmonic generation[J]. Physical Review Letters, 2007, 99(5):053904. doi: 10.1103/PhysRevLett.99.053904 [15] LU R F, ZHANG P Y, HAN K L. Attosecond-resolution quantum dynamics calculations for atoms and molecules in strong laser fields[J]. Physical Review, 2008, E77(6):066701. [16] HU J, HAN K L, HE G Z. Correlation quantum dynamics between an electron and D2+ molecule with attosecond resolution[J]. Physical Review Letters, 2005, 95(12):123001. doi: 10.1103/PhysRevLett.95.123001 [17] FENG L Q, CHU T S. Generation of an isolated sub-40as pulse using two-color laser pulses:Combined chirp effects[J]. Physical Review, 2011, A84(5):053853. [18] FENG L Q. Molecular harmonic extension and enhancement from H2+ ions in the presence of spatially inhomogeneous fields[J]. Physical Review, 2015, A92(5):053832. [19] LIU Sh Sh, MIAO X Y. Enhancement of high-order harmonic emission by using a coherent superposition in a two-color laser field[J]. Journal of Atomic and Molecular Physics, 2012, 29(5):881-885(in Chinese). [20] DU H N, MIAO X Y. Theoretical exploration of intensity effects on high-order harmonic generation from H2+ molecules[J]. Journal of Atomic and Molecular Physics, 2013, 30(1):115-120(in Chinese). [21] BURNETT K, REED V C, COOPER J, et al. Calculation of the background emitted during high-harmonic generation[J]. Physical Review, 1992, A45(5):3347-3349. [22] ANTOINE P, PIRAUX B, MAQUET A. Time profile of harmonics generated by a single atom in a strong electromagnetic field[J]. Physical Review, 1995, A51(3):R1750-R1753. [23] SILVA R E F, CATOIRE F, RIVIERE P, et al. Correlated electron and nuclear dynamics in strong field photoionization of H2+[J]. Physical Review Letters, 2013, 110(11):113001. doi: 10.1103/PhysRevLett.110.113001 -
点击查看大图
图(6)
计量
- 文章访问数: 6560
- HTML全文浏览量: 4878
- PDF下载量: 272
- 被引次数: 0
