-
利用强场近似理论开展了具有较长脉宽的偏振控制脉冲与氦原子相互作用产生高次谐波和阿秒脉冲发射的理论研究. 研究发现, 当具有10 fs脉冲宽度的偏振控制脉冲被用作驱动脉冲时, 只要恰当地调整两束反向旋转圆偏振脉冲峰值之间的时间延迟和强度比, 即使不附加二次谐波脉冲, 仍然可以得到效率较高且规则分布的高次谐波平台结构, 傅里叶变换后得到了175 as的孤立短脉冲. 该方案一方面通过调整两束脉冲峰值之间时间延迟突破了传统偏振控制方案中要求偏振门宽度为半个光学周期的限制, 另一方面通过调整两束脉冲峰值之间的强度比避免了偏振门前端多个光学周期电场引起气体介质电离不利于谐波相位匹配的弊端.Isolated attosecond pulses make it possible to study and control the ultrafast electron processes in atoms and molecules. High order harmonic generation (HHG) is the most promising way to generate such pulses, which is benefited from the broad plateau structure of the typical HHG spectrum. In previous HHG studies on the polarization gating pulse with longer pulse duration, one needs to dramatically increase the separation in time between the two counter-rotating circularly-polarized pulses to generate the nearly-linear half-cycle " polarization gate”. This leads to a low harmonic conversion efficiency because the field outside the polarization gate is much stronger than inside the polarization gate. In this paper, by using Lewenstein model, we theoretically simulate the high-order harmonic generation from helium atom subjected to the polarization gating pulse with 10 fs pulse duration. It is found that high-order harmonic spectra each with a higher efficiency and regular structure can still be obtained by reasonably adjusting the delay time ratio and the amplitude ratio of electric fields between the two counter-rotating pulses. Further, a single 175 as pulse in the time domain is obtained by Fourier transforming the 80th order harmonics into the 172nd order harmonic without compensating for the harmonic chirp. This scheme has two main advantages. First, the adjustment of the polarization gate width from half optical cycle into nearly one cycle ensures higher intensity of the synthesized electric field inside the polarization gate. Second, the suitable adjustment of the amplitude ratio between two electric fields ensures the low ionization probability before the polarization gate, and thus further fulfills the harmonic phase matching condition in the process of the propagation.
1. 引 言
孤立阿秒脉冲可用以跟踪和控制原子、分子及固体中电子的运动过程而备受人们关注[1,2]. 高次谐波发射谱因具有独特的超连续平台结构成为目前实现阿秒脉冲的唯一有效手段[3-5].
高次谐波发射是指强激光脉冲与原子、分子、团簇以及固体等介质相互作用时产生的相干辐射波[6]. 该过程可由半经典的“三步” 模型[7]给出合理解释: 当激光场辐照原子介质时, 处在束缚态的电子隧穿由原子势和激光电场形成的势垒, 电子电离后在振荡的电场中运动, 当电场反向时, 部分电离电子回到母离子附近, 与基态复合, 放出高能光子, 且光子的能量为
$h\upsilon = {I_{\rm{p}}} + {E_{{\rm{kin}}}}$ , 其中Ip是原子的电离势, Ekin为电子与原子基态复合时具有的动能. 之前的研究结果表明, 所有高次谐波辐射谱均呈现如下特征行为: 最初几次谐波, 随着谐波次数增加发射效率快速下降; 接着, 在达到某次谐波后出现一个强度几乎不随谐波次数变化的区域(称为平台区); 最后, 在平台区末端的某一阶次谐波附近, 谐波强度又快速下降, 出现截止. 由于谐波谱呈现平台区以及平台区谐波的独特优点, 使它成为获得极紫外及软X 射线波段阿秒脉冲相干辐射源的重要手段[8].当前基于高次谐波发射获得孤立阿秒脉冲的研究主要有以下几种技术: 太赫兹场方案[9]、偏振控制方案[4,10,11]、少周期方案[12]、电离控制方案[13]、阿秒灯塔方案[6]以及最近提出的相对论等离子体镜方案[14]. 迄今为止, 瑞士苏黎世联邦理工学院Gaumnitz等[15]已实现了最短43 as孤立短脉冲的产生.
偏振控制方案主要是利用高次谐波强度对激光偏振状态敏感的特性获得孤立阿秒脉冲产生. 基于偏振脉冲的椭偏率为0.2时的高次谐波强度比线偏振光情形减少一个数量级, Corkum 等[10]提出控制驱动脉冲光场在半个光学周期的时间窗口里是线性偏振, 而其余时间都是椭圆偏振的椭偏率随时间变化的激光场产生高次谐波, 最终获得了单一阿秒脉冲. 事实上, 该驱动光场可由具有一定时间延迟的两束反旋圆偏振光合成. 实验上, 通常由全波片和零级四分之一波片实现偏振控制方案中的驱动光场. 当全波片光轴与入射线偏振脉冲光轴夹角为45°时, 入射线偏振脉冲通过该波片后分成两束垂直偏振且强度相等的线偏振脉冲, 两束脉冲之间时间延迟由全波片厚度和折射率控制. 若零级四分之一波片光轴与全波片光轴夹角为45°时, 具有一定时间延迟的两束等强度线偏振脉冲转换为两束等强度反旋圆偏振脉冲, 若合理地调整两束脉冲峰值之间的时间延迟, 可在交叠区产生具有线性偏振特点的半个光学周期的偏振门.
相比于阿秒脉冲产生的其他方案, 偏振控制方案的优势是在实验上较易操作, 且谐波辐射仅发生在偏振门内的线性半个光学周期内, 因此可得到在整个平台区及截止位置附近都连续的谐波谱, 满足阿秒脉冲产生对超宽光谱带宽的要求. 然而, 该方案中为了获得较强的孤立阿秒脉冲, 需采用短脉宽的驱动脉冲获得半个光学周期的偏振门. Chang[16]利用脉宽5 fs的驱动脉冲及5 fs的时间延迟在理论上得到了58 as的孤立短脉冲. 然而, 相比于5 fs超短脉冲的产生, 具有较长脉冲宽度的飞秒脉冲在实验中较易实现, 且脉冲强度也较强. 因此, 有必要开展利用较长脉宽的偏振控制脉冲方案获得孤立阿秒脉冲的研究. 当前, 对于具有较长脉冲宽度的偏振控制脉冲, 为了确保偏振门的宽度为线性半个光学周期, 两束反向旋转圆偏振脉冲峰值之间的时间延迟应大幅增加, 这导致了偏振门外电场强度远高于偏振门内的场强, 因此得到的谐波谱效率很低. 为了克服这个困难, Zhao等[17]在2012年提出了双光学控制方案, 即在脉宽为7 fs偏振控制脉冲基础上附加一束二次谐波脉冲, 从而将偏振门的宽度从半个周期释放到一个光学周期, 最终得到了67 as的孤立短脉冲. 然而, 二次谐波脉冲的加入增加了实验操作的困难.
本文说明了当采用具有10 fs脉冲宽度的偏振控制脉冲时, 只要恰当地调整两束反向旋转圆偏振脉冲之间的时间延迟和强度比, 即使不附加二次谐波脉冲, 仍然可以得到有效孤立阿秒脉冲的产生. 调整两束脉冲峰值之间时间延迟的目的在于将偏振门宽度从半个光学周期调整到接近一个周期, 偏振门宽度释放使得偏振门内合成脉冲电场强度增加到与偏振门外的场强接近一致, 克服了激光场能量大量损失导致谐波效率低的弊端. 另外, 在通常偏振控制方案研究中, 偏振门内用以产生谐波的电场宽度为半个光学周期, 其中前1/4周期电场负责谐波产生过程中原子的电离, 后1/4周期电场负责电离电子与母离子的复合. 当偏振门的宽度调整到接近一个光学周期时, 其中前半个光学周期电场负责谐波产生过程中原子的电离, 然而研究发现, 这半个光学周期的前1/4周期电场仅对40阶次附近谐波产生有贡献, 平台及截止谐波的产生仍然来自于后1/4周期电场, 因此通过调整两束脉冲峰值之间的时间延迟有望得到规则平台且效率较高的谐波谱. 两束反旋圆偏振脉冲峰值电场之间强度比的调整不仅可以使偏振门前端电场峰值强度低于后端电场强度, 而且调整后偏振门前端的光学周期数也有所减少, 这将避免具有较长脉宽驱动脉冲因偏振门前端多个光学周期电场引起气体介质电离不利于谐波相位匹配的弊端, 同时也确保了偏振门内有足够电子电离实现有效高次谐波发射.
2. 理论方法
高次谐波产生包括单原子响应和谐波在介质中传播两个过程. 单原子响应是指单个原子在激光电场作用下的高次谐波发射, 而介质中所有原子总的高次谐波发射可以通过传播方程给出. 考虑传播过程只在增强阿秒脉冲和高次谐波发射强度及消除长轨迹等方面存在影响, 因此本文研究忽略传播过程. 本文利用Lewenstein强场近似理论模型[18-20]数值模拟了10 fs偏振控制脉冲辐照下氦原子的高次谐波发射, 以此探究用较长脉宽的偏振控制方案获得孤立阿秒脉冲的研究. 由于激光脉冲的椭偏率随时间变化, 因此需要分别计算x, y两个不同方向上的偶极矩, 沿y方向的谐波谱强度远低于沿x方向的谐波谱强度, 这一点已经在之前的研究中得到验证. 因此只需计算出x方向上的偶极矩, 其具体表达式如下[21](除非另有说明, 否则本文均采用原子单位):
$\begin{split} x\left( t \right) \approx \, &{\rm{i}}\int_{ - \infty }^t {\rm{d}}t'{\left[ {\frac{{\text{π}}}{{\varepsilon + {\rm{i}}\left( {t - t'} \right)/2}}} \right]^{3/2}}\\ & \cdot {d_x}^*\left[ {{{p}_{{\rm{st}}}}\left( {t',t} \right) - {A}\left( t \right)} \right],\\ & \cdot\exp \left[ { - {\rm{i}}{S_{{\rm{st}}}}\left( {{{p}_{{\rm{st}}}},t',t} \right)} \right]\\ & \times {d}\left[ {{{p}_{{\rm{st}}}}\left( {t',t} \right) - {A}\left( {t'} \right)} \right] \\ & \cdot {E}\left( {t'} \right)g\left( t \right) + {\rm{c}}.{\rm{c}}., \end{split}$
(1) 式中ε是一个正的小数, E(t)是脉冲电场, A(t)是它的相关矢势. 基态布居数
$g\left( t \right) = {\rm{exp}}\left[ - \int_{ - \infty }^t {\omega _{{\rm{ADK}}}}\left( {t''} \right){\rm{d}}{t^{''}}\right],$
其中, 电离速率
${\omega _{{\rm{ADK}}}}\left( {t''} \right)$ 由Ammosov-Delone-Krainov (ADK)理论计算[22].电子的准经典作用表示为
$\begin{split}{S_{{\rm{st}}}}\left( {{{p}_{{\rm{st}}}},t',t} \right) =\, & \left( {t - t'} \right){I_{\rm{p}}} - \frac{1}{2}p_{{\rm{st}}}^2\left( {t',t} \right)\left( {t - t'} \right) \\ &+ \frac{1}{2}\int_{t'}^t {A^2}({t{''}}){\rm{d}}t'',\end{split}$
(2) 这里Ip是氦原子的电离势, pst是与电子定态相位相关的正则动量, 其表示为
${{p}_{{\rm{st}}}}\left( {t',t} \right) = \frac{1}{{t - t'}}\int_{t'}^t {A}(t''){\rm{d}}t''.$
(3) 最后, 我们给出了(1)式中基态和连续态之间的无场偶极跃迁矩阵元
$\begin{split} &{d_x}\left[ {{{p}_{{\rm{st}}}}\left( {t',t} \right) - {A}\left( t \right)} \right] = {\rm{i}}\frac{{{2^{7/2}}}}{{\text{π}}}{\left( {2{I_{\rm{p}}}} \right)^{5/4}}\\ \cdot &\frac{{{p_{{\rm{st}},x}}\left( {t',t} \right) \!-\! {A_x}\left( {{t}} \right)}}{{{{\left\{ {{{\left[ {{p_{{\rm{st}},x}}\left( {t',t} \right) \!-\! {A_x}\left( {{t}} \right)} \right]}^2} \!+\! {{\left[ {{p_{{\rm{st}},y}}\left( {t',t} \right) \!-\! {A_y}\left( {{t}} \right)} \right]}^2} \!+\! 2{I_{\rm{p}}}} \right\}}^3}}}.\end{split}$
(4) 同样, 我们也可以给出
$\begin{split} & {d}\left[ {{{p}_{{\rm{st}}}}\left( {t',t} \right) - {A}\left( {t'} \right)} \right]{E}\left( {t'} \right) = {\rm{i}}\frac{{{2^{7/2}}}}{{\text{π}}}{\left( {2{I_{\rm{p}}}} \right)^{5/4}} \\ \cdot & \frac{{\left[ {{p_{{\rm{st}},x}}\!\left( {t',t} \right) \!-\!\! {A_x}(t')} \right]\!{E_x}(t') \!+\! \left[ {{p_{{\rm{st}},y}} (t',t)\! -\! \!{A_y} (t')} \right]\!{E_y} (t')}}{{{{\left\{ \!{{{\left[ {{p_{{\rm{st}},x}}\!\left( {t',t} \right)\! -\! {A_x}\left( {t'} \right)} \right]}^2} \!+\! {{\left[ {{p_{{\rm{st}},y}}\!\left( {t',t} \right) \!-\! {A_y}\left( {t'} \right)} \right]}^2}\! +\! 2{I_{\rm{p}}}} \right\}}^3}}}. \end{split}$
(5) 通过傅里叶变换偶极矩, 可以得到来自单个原子的谐波谱. 需要指出的是, 本文考虑了高于30次的谐波辐射, 30次谐波的光子能量为46 eV, 明显大于氦原子的电离势. 根据文中采用的激光脉冲参数可计算出偏振门中心的 Keldysh参数γ = 0.54 < 1, 可见 Lewenstein模型在这些给定条件下有效.
通过叠加谐波频谱上超连续分布谐波, 可得到超短阿秒脉冲的产生
$I\left( t \right) = {\left| {\mathop \sum \nolimits_q {a_q}{{\rm{e}}^{{\rm{i}}q\omega t}}} \right|^2},$
(6) 其中
${a_q} = \int a \left( t \right){{\rm{e}}^{ - {\rm{i}}q\omega t}}{\rm{d}}t$ .3. 结果与分析
在我们的数值模拟中, 由左旋和右旋圆偏振高斯脉冲叠加形成的椭偏率随时间变化的偏振控制脉冲电场形式为
$\hat E\left( t \right) = {E_{{\rm{drive}}}}\left( t \right)\hat x + {E_{{\rm{gate}}}}\left( t \right)\hat y$ , 其中驱动脉冲和控制脉冲电场形式依次为$\begin{split}{E_{{\rm{drive}}}} =\, & \frac{{{E_0}}}{{\sqrt 2 }}\left({\rm{sin}}\theta \exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t + {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}} \right]}^2}} \right\}\right. \\ &+\left. {\rm{cos}}\theta \exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t - {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}}\right]}^2}} \right\}\right)\\ & \cdot{\rm{cos}}\left( {\omega t + \varphi }\right),\end{split}$
(7) $\begin{split} {E_{{\rm{gate}}}} = \, &\frac{{{E_0}}}{{\sqrt 2 }}\left({\rm{sin}}\theta \exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t + {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}} \right]}^2}} \right\}\right.\\ & - \left.{\rm{cos}}\theta \exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t - {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}} \right]}^2}} \right\}\right)\\ & \cdot{\rm{sin}}\left( {\omega t + \varphi } \right),\end{split}$
(8) 这里E0和ω分别是入射脉冲激光的电场强度和载波频率, Td是两束脉冲峰值之间的时间延迟, τp是脉冲宽度, φ是载波包络相位. 角度θ是全波片光轴与入射线偏振脉冲光轴之间的夹角, 当θ = 45°时, 偏振控制方案中两束反旋圆偏振脉冲的峰值电场强度比为1, 我们称其为对称偏振控制方案, 本文中研究偏振门宽度效应时采用了此方案. 当θ偏离45°时, 两束反旋圆偏振脉冲的峰值电场强度不再相等, 此情形称为不对称偏振控制方案. 在不对称偏振控制方案中, 若θ < 45°, 右旋圆偏振光的峰值场强度低于左旋圆偏振光的峰值场强度, 即偏振门前端电场强度低于后端电场强度, 本文中考虑两束脉冲峰值之间强度比效应时采用了不对称偏振控制方案.
在对称偏振控制方案中, 合成脉冲电场的含时椭偏率为
$\xi \left( t \right) = \frac{{\left| {\exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t + {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}} \right]}^2}} \right\} - \exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t - {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}} \right]}^2}} \right\}} \right|}}{{\exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t + {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}} \right]}^2}} \right\} + \exp \left\{ { - 2{\rm{ln}}2{{\left[ {\left( {t - {T_{\rm d}}/2} \right)/{\tau _{\rm{p}}}} \right]}^2}} \right\}}}.$
(9) 对于高于21次的谐波, 当椭偏率从0增加到0.2时, 谐波强度下降超过一个数量级. 因此, 阿秒脉冲仅在t = tc (tc为完全线偏振时刻)附近且椭偏率ξ ≤ 0.2的时间范围内产生, 该时间范围称为偏振门宽度, 其表达式为
${\text{δ}}{t_{\rm{G}}} = \frac{{{\xi _{{\rm{th}}}}}}{{{\rm{ln}}2}}\frac{{\tau _{\rm{p}}^2}}{{{T_{\rm{d}}}}},$
(10) 式中阈值椭偏率ξth取值小于0.2. 从(10)式可以看到偏振门宽度与两束圆偏振脉冲峰值之间的时间延迟Td及入射脉冲的脉宽τp有关. 通常情况, 为了获得孤立阿秒脉冲, 偏振门宽度应该短至T0/2 (T0是入射脉冲光场的一个光学周期), 因此可通过增加两束脉冲峰值之间的时间延迟以及采用短脉宽超短脉冲两种方法减少偏振门宽度.
实验中人们通常选取脉宽为5 fs的钛宝石脉冲激光作为驱动脉冲, 为了确保偏振门的宽度为半个光学周期, 两束脉冲峰值之间的时间延迟应取为5.5 fs, 此时Td
$\approx $ τp, 偏振门内合成电场的峰值振幅$E\left( {{t_{\rm{c}}}} \right) = \sqrt 2 {E_0}$ , 很显然, 偏振门内合成脉冲电场强度高于入射驱动脉冲的电场强度, 图1(a)和图1(b)分别展示了偏振控制脉冲总电场及驱动脉冲电场随时间变化的曲线图, 图1(b)中的阴影部分是偏振门的位置. 当选取脉宽为10 fs的偏振控制脉冲辐照时, 为了确保偏振门的宽度为半个光学周期, 两束脉冲峰值之间的时间延迟应取为22.5 fs, 此时Td$\gg $ τp, 偏振门内合成电场的峰值振幅为$E\left( {{t_{\rm{c}}}} \right) = {E_0}{\rm{exp}} \Big[ - \dfrac{{{\rm{ln}}2}}{2}{{\Big( {\dfrac{{{T_{\rm{d}}}}}{{{\tau _{\rm{p}}}}}} \Big)}^2} \Big]$ , 明显低于每一束脉冲的峰值场振幅, 图1(c)和图1(d)分别展示了该情形下偏振控制脉冲总电场及驱动脉冲电场随时间的变化, 图1(d)中的阴影部分是偏振门的位置. 可以看到, 偏振门外的场强远远高于偏振门内的场强, 激光场能量损失将导致得到的谐波谱转换效率很低. 然而, 若将两束反旋圆偏振脉冲峰值之间的时间延迟减小到15 fs, 偏振门的宽度将从半个光学周期调整到接近一个光学周期, 如图1(e)和图1(f)所示, 此时偏振门内合成脉冲电场强度与偏振门外的场强接近一致, 显然偏振门外激光电场能量损失导致谐波谱转换效率低的问题将得到缓解.图2中的红色实线给出了脉宽为10 fs且偏振门宽度接近一个光学周期时的偏振控制脉冲与氦原子相互作用得到的高次谐波发射谱. 数值模拟过程中入射钛宝石驱动脉冲的载波频率ω为0.057, 脉冲的载波包络相位
$\varphi = {\text{π}}/2$ . 两束脉冲峰值之间的时间延迟Td取为15 fs, 对应的偏振门宽度为0.82 T0. 当输入脉冲的峰值场振幅为0.35 a.u. 时, 偏振门内的峰值场幅则为0.18 a.u.. 从图2可以看到, 谐波频谱的截止(频率)位置达到158次谐波, 且整个平台区及截止位置都呈现了规则分布的特点. 为了比较, 我们也给出了脉宽为10 fs且偏振门宽度为半个光学周期时的偏振控制脉冲与氦原子相互作用得到的高次谐波发射谱, 如图2中的黑色实线所示. 计算中脉冲的载波包络相位仍取为$\varphi = {\text{π}}/2$ , 两束脉冲峰值之间的时间延迟Td取为22.5 fs, 对应的偏振门宽度为0.5 T0. 当输入脉冲的峰值场振幅为0.35 a.u.时, 偏振门内的峰值场幅则为0.09 a.u.. 从图2可以看到, 谐波频谱的截止频率为55次, 而且具有规则结构的平台宽度仅为25次, 特别注意的是谐波平台的效率低了2个数量级. 如果调整偏振门内峰值场振幅为0.18 a.u., 此时需要入射驱动脉冲的峰值场振幅达到0.84 a.u.. 在该参数下, 偏振门到达前原子中的电子已全部电离, 不能得到高次谐波发射谱. 可见, 若想利用具有较长脉冲宽度的偏振控制脉冲方案获得有效孤立阿秒脉冲的产生, 应合理地调整两束反向旋转圆偏振脉冲峰值之间的时间延迟.为了解释偏振门宽度接近一个光学周期时谐波平台及截止位置呈规则分布的特点, 利用半经典三步模型理论计算了偏振门内不同时刻电离电子的复合时刻及该时刻电子具有的动能, 如图3所示. 图中黑色和红色实线分别给出了谐波随着电离时刻以及复合时刻的变化曲线. 根据半经典理论, 阿秒脉冲辐射在驱动激光场的每个光学周期内发生两次. 如果将偏振门宽度调整为一个光学周期, 由于两次谐波辐射之间的干涉效应, 应该获得不规则的高次谐波频谱. 但从图3可以看出, 若将偏振门中心tc之前半个光学周期分成前后两个1/4周期, 那么前1/4周期(区域I)电离的电子仅对40次左右的谐波有贡献, 平台区和截止区谐波产生仍然主要来源于后1/4周期电离电子的贡献(区域II), 因此我们将偏振门宽度增加到0.82个光学周期时, 仍然能够获得仅由单一长短轨迹构成的规则平台结构. 这里长轨迹是指在电场中先电离后复合的电子运动轨迹, 短轨迹指在电场中后电离而先复合的电子运动轨迹.
偏振门宽度不同导致谐波谱平台转换效率差别可以通过原子电离概率得到合理解释. 对于脉宽为10 fs的偏振控制脉冲, 偏振门宽度增加减小了两束脉冲峰值之间的时间延迟, 进而使偏振门内的合成电场强度大幅增加, 减少了电场在偏振门外的损失, 正如图1(f)所展示的(对比图1(d)). 很显然, 当该合成脉冲辐照氦原子时, 原子的电离概率将发生明显改变, 正如图4(b)中的红色实线所展示. 为了方便分析, 图中也展示了驱动脉冲电场强度(黑色曲线)以及偏振控制脉冲的椭偏率(蓝色曲线)随时间的变化曲线. 可以看到, 在偏振门内第二个1/4周期(图4(b)中的II区)开始时原子的电离概率为37.69%, 在偏振门中心原子的电离概率为38.92%, 偏振门内用以产生谐波的有效原子电离概率为1.23%, 远超偏振门为半个光学周期时的有效原子电离几率0.004%, 如图4(a)中的红色实线所示. 因此, 偏振门宽度从半个光学周期调整到接近一个光学周期后高次谐波发射谱的效率提高了2个数量级.
综上所述, 当采用具有较长脉宽的偏振控制脉冲方案时, 通过合理调整偏振门宽度可以得到在整个平台区和截止位置都规则分布的高次谐波发射谱, 然而, 该方案中, 我们发现在偏振门到达之前, 原子的电离概率已经高达36.69%, 这部分电离不仅不利于谐波产生过程中的相位匹配, 而且由于原子基态布居大幅衰退影响偏振门内高次谐波产生的效率. 接下来尝试调整两束反向旋转圆偏振脉冲峰值之间强度比使偏振门前端的电场强度小于偏振门后端电场强度(即前面所述的不对称偏振控制方案), 避免具有较长脉宽驱动脉冲因偏振门前端多个光学周期电场引起气体介质电离衰退原子基态布居的弊端, 从而确保偏振门内有足够电子电离实现较高效的高次谐波发射.
图5(a)展示了全波片光轴与入射线偏振脉冲光轴之间偏振角θ为20°时的不对称偏振控制脉冲电场随时间变化曲线图, 可以看到偏振门前端的峰值电场强度明显小于偏振门后端的峰值电场强度, 而且, 偏振门中心的位置已经从0前移到–2.7 fs, 明显减小了偏振门前端光学周期数. 图5(b)展示了该不对称偏振控制脉冲与氦原子相互作用得到的高次谐波发射谱. 计算中选择的脉冲激光参数同对称偏振控制脉冲一致. 偏振门内的电场峰值强度为0.18 a.u., 两束脉冲峰值之间的时间延迟为15 fs, 两束脉冲的脉宽为10 fs. 为了比较, 我们也给出了对称偏振控制脉冲与氦原子相互作用得到的高次谐波发射谱, 如图5(b)中的红色实线所示. 可以看到, 利用不对称偏振控制脉冲得到的谐波谱强度明显高于对称偏振控制脉冲情形, 且谐波谱截止位置从156次延展到了172次. 下面从原子电离概率分布说明谐波谱效率提高的原因. 图4(c)给出了不对称偏振控制脉冲作用下原子的电离概率随时间变化曲线图, 可以看到在偏振门的第二个1/4周期(图4(c)中的II区)开始时原子的电离概率为3.37%, 偏振门中心原子的电离概率为4.75%, 有效原子电离概率达到1.38%, 高于对称偏振控制方案中的有效电离概率1.23%. 更重要的是, 在不对称偏振控制方案中, 原子的低电离概率有利于谐波在介质传播过程中满足相位匹配条件.
为了说明图5(b)中谐波谱平台及截止位置呈现规则分布特点的原因, 在图6中利用小波变换方法计算了不对称偏振控制脉冲辐照下高次谐波发射的时频分析图像. 图中平台上的谐波主要由长短两个量子轨迹构成, 长轨迹对应于–2 fs附近时刻的发射峰, 短轨迹对应于–1 fs附近时刻的发射峰, 两个轨迹干涉后得到了规则分布的谐波谱平台结构. 但我们也发现在0.5 fs位置附近有50—75次谐波的辐射, 这部分谐波主要来自于在偏振门后半个光学周期电离的电子回到离核较远的地方复合得到的, 其低效率使得谐波谱上50—75次谐波有轻微的调整, 通过放大图5(b)中的黑色实线可以看到这一现象.
图7(a)和图7(b)分别展示了不对称和对称偏振控制方案中叠加谐波谱上呈规则分布谐波生成的阿秒脉冲时域图. 偏振角为20°时, 通过叠加80—172次谐波, 得到了一个阿秒脉冲链, 如图7(a)所示. 第一个阿秒脉冲a对应着谐波发射的短轨迹, 其半高全宽为175 as, 第二个阿秒脉冲b对应着谐波发射的长轨迹, 其半高全宽为177 as, 最后一个弱阿秒脉冲c对应着离核较远电子与核复合产生的辐射. 而对称偏振控制脉冲情形, 通过叠加80—156次谐波, 同样得到了一个阿秒脉冲链, 如图7(b)所示. 第一个阿秒脉冲a′对应着谐波发射的短轨迹, 其半高全宽为179 as, 第二个阿秒脉冲b′对应着谐波发射的长轨迹, 其半高全宽为179 as, 最后一个弱阿秒脉冲仍然来自于离核较远电子与核复合产生的辐射. 在高次谐波产生的数值模拟中, 长轨迹b和b′的贡献可通过求解传播方程被抑制, 因此最终得到的仍然是单一阿秒脉冲a和a′. 更重要的是, 当采用不对称偏振控制脉冲情形, 单一阿秒脉冲的强度提高了两倍, 这与上述原子电离概率的分析结果是一致的.
4. 结 论
本文利用三维强场近似理论说明了脉宽为10 fs的两束不对称反旋圆偏振控制脉冲可用以获得强度较高孤立阿秒脉冲的产生. 不同于通常偏振控制脉冲方案, 为了获得有效的高次谐波发射, 我们一方面将两束脉冲峰值之间的时间延迟从半个光学周期调整到接近一个光学周期, 从而确保了偏振门内合成脉冲电场具有较高的强度, 另一方面调整了两束脉冲的峰值强度比, 确保了偏振门前端原子低电离概率, 这使得谐波产生过程中的相位匹配条件容易得到满足. 相比于双光学控制方案, 该方案具有实验操作相对容易的优势.
[1] Schultze M, Fiess M, Karpowicz N, Gagnon J, Korbman M, Hofstetter M, Neppl S, Cavalieri A L, Komninos Y, Mercouris T, Nicolaides C A, Pazourek R, Nagele S, Feist J, Burgdörfer J, Azzeer A M, Ernstorfer R, Kienberger R, Kleineberg U, Goulielmakis E, Krausz F, Yakovlev V S 2010 Science 328 1658Google Scholar
[2] Drescher M, Hentschel M, Kienberger R, Uiberacker M, Yakovlev V, Scrinzi A, Westerwalbesloh T, Kleineberg U, Heinzmann U, Krausz F 2002 Nature 419 803Google Scholar
[3] Sansone G, Benedetti E, Calegari F, Vozzi C, Avaldi L, Flammini R, Poletto L, Villoresi P, Altucci C, Velotta R, Stagira S, De Silvestri S, Nisoli M 2006 Science 314 443Google Scholar
[4] Goulielmakis E, Schultze M, Hofstetter M, Yakovlev V S, Gagnon J, Uiberacker M, Aquila A L, Gullikson E M, Attwood D T, Kienberger R, Krausz F, Kleineberg U 2008 Science 320 1614Google Scholar
[5] Vincenti H, Quéré F 2012 Phys. Rev. Lett. 108 113904Google Scholar
[6] Chen J G, Yang Y J, Chen J, Wang B B 2015 Phys. Rev. A 91 043403Google Scholar
[7] Corkum P B 1993 Phys. Rev. Lett. 71 1994Google Scholar
[8] Agostini P, Dimauro L F 2004 Rep. Prog. Phys. 67 813Google Scholar
[9] Yuan K J, Bandrauk A D 2013 Phys. Rev. Lett. 110 023003Google Scholar
[10] Corkum P B, Burnett N H, Ivanov M Y 1994 Opt. Lett. 19 1870Google Scholar
[11] Li J, Ren X, Yin Y, Zhao K, Chew A, Cheng Y, Cunningham E, Wang Y, Hu S, Wu Y, Chini M, Chang Z 2017 Nat. Commun. 8 794Google Scholar
[12] Ferrari F, Calegari F, Lucchini M, Vozzi C, Stagira S, Sansone G, Nisoli M 2010 Nat. Photonics 4 875Google Scholar
[13] Lan P, Lu P, Cao W, Li Y H, Wang X L 2007 Phys. Rev. A 76 21801Google Scholar
[14] Kormin D, Borot A, Ma G J, Dallari W, Bergues B, Aladi M, Földes I, Veisz L 2018 Nat. Commun. 9 4992Google Scholar
[15] Gaumnitz T, Jain A, Pertot Y, Huppert M, Jordan I, Ardana-Lamas F, Wörner H J 2017 Opt. Express 25 27506Google Scholar
[16] Chang Z 2005 Phys. Rev. A 71 023813Google Scholar
[17] Zhao K, Zhang Q, Chini M, Wu Y, Wang X W, Chang Z H 2012 Opt. Lett. 37 3891Google Scholar
[18] Keldysh L V 1965 Sov. Phys. JETP 20 1307
[19] Faisal F H M 1973 J. Phys. B 6 L89Google Scholar
[20] Reiss H R 1980 Phys. Rev. A 22 1786Google Scholar
[21] Lewenstein M, Salières P, L’Huillier A 1995 Phys. Rev. A 52 4747Google Scholar
[22] Ammosov M V, Delone N B, Krainov V 1986 Proc. SPIE 664 1191Google Scholar
-
图 1 偏振控制脉冲总电场(红色曲线)、控制场(绿色曲线)及驱动场(蓝色曲线)随时间变化三维图 (a) 5 fs脉宽和5 fs时间延迟; (c) 10 fs脉宽和22.5 fs时间延迟; (e) 10 fs脉宽和15 fs时间延迟; (b), (d), (f)显示了与(a), (c), (e)图相对应的驱动脉冲电场随时间变化曲线图(阴影部分是偏振门)
Fig. 1. Three dimensional diagrams for the total electric field (red), gating field (green), and driving field (blue) in polarization gating pulse as a function of time: (a) 5 fs pulse width and 5 fs time delay; (c) 10 fs pulse width and 22.5 fs time delay; (e) 10 fs pulse width and 15 fs time delay. Panels (b), (d), (f) correspond to the driving electric field versus time in panels (a), (c), (e), respectively (shaded portion is polarization gate).
图 2 脉宽为10 fs的偏振控制脉冲与氦原子相互作用得到的高次谐波发射光谱, 其中黑线表示δtG = T0/2的偏振控制脉冲, 红线表示δtG = 0.82T0的偏振控制脉冲
Fig. 2. High order harmonic generation from helium atom in a polarization gating pulse with 10 fs pulse width. The black curve is from the polarization gate width
${\text{δ}}{t_{\rm{G}}} = \dfrac{{{T_0}}}{2}$ , the red curve is from the polarization gate width δtG = 0.82T0.图 4 偏振门内原子电离概率(红线)、驱动脉冲电场(黑线)、椭偏率(蓝线)随时间变化曲线 (a)对称偏振控制方案δtG = 0.5T0; (b)对称偏振控制方案δtG = 0.82T0; (c)不对称偏振控制方案δtG = 0.82T0
Fig. 4. Atomic ionization probability (red line), electric field of driving pulse (black line), ellipticity (blue line) in the polarization gate as function of time: (a) Symmetric polarization gating scheme δtG = 0.5T0; (b) symmetric polarization gating scheme δtG = 0.82T0; (c) asymmetric polarization gating scheme δtG = 0.82T0.
图 5 (a)不对称偏振控制脉冲总电场(红色曲线)、控制场(绿色曲线)及驱动场(蓝色曲线)随时间变化三维图; (b)高次谐波发射谱
Fig. 5. (a) Three-dimensional diagrams for the total electric field (red curve), gating field (green curve) and driving field (blue curve) in the asymmetric polarization gating pulse as a function of time; (b) high harmonic emission spectra.
-
[1] Schultze M, Fiess M, Karpowicz N, Gagnon J, Korbman M, Hofstetter M, Neppl S, Cavalieri A L, Komninos Y, Mercouris T, Nicolaides C A, Pazourek R, Nagele S, Feist J, Burgdörfer J, Azzeer A M, Ernstorfer R, Kienberger R, Kleineberg U, Goulielmakis E, Krausz F, Yakovlev V S 2010 Science 328 1658Google Scholar
[2] Drescher M, Hentschel M, Kienberger R, Uiberacker M, Yakovlev V, Scrinzi A, Westerwalbesloh T, Kleineberg U, Heinzmann U, Krausz F 2002 Nature 419 803Google Scholar
[3] Sansone G, Benedetti E, Calegari F, Vozzi C, Avaldi L, Flammini R, Poletto L, Villoresi P, Altucci C, Velotta R, Stagira S, De Silvestri S, Nisoli M 2006 Science 314 443Google Scholar
[4] Goulielmakis E, Schultze M, Hofstetter M, Yakovlev V S, Gagnon J, Uiberacker M, Aquila A L, Gullikson E M, Attwood D T, Kienberger R, Krausz F, Kleineberg U 2008 Science 320 1614Google Scholar
[5] Vincenti H, Quéré F 2012 Phys. Rev. Lett. 108 113904Google Scholar
[6] Chen J G, Yang Y J, Chen J, Wang B B 2015 Phys. Rev. A 91 043403Google Scholar
[7] Corkum P B 1993 Phys. Rev. Lett. 71 1994Google Scholar
[8] Agostini P, Dimauro L F 2004 Rep. Prog. Phys. 67 813Google Scholar
[9] Yuan K J, Bandrauk A D 2013 Phys. Rev. Lett. 110 023003Google Scholar
[10] Corkum P B, Burnett N H, Ivanov M Y 1994 Opt. Lett. 19 1870Google Scholar
[11] Li J, Ren X, Yin Y, Zhao K, Chew A, Cheng Y, Cunningham E, Wang Y, Hu S, Wu Y, Chini M, Chang Z 2017 Nat. Commun. 8 794Google Scholar
[12] Ferrari F, Calegari F, Lucchini M, Vozzi C, Stagira S, Sansone G, Nisoli M 2010 Nat. Photonics 4 875Google Scholar
[13] Lan P, Lu P, Cao W, Li Y H, Wang X L 2007 Phys. Rev. A 76 21801Google Scholar
[14] Kormin D, Borot A, Ma G J, Dallari W, Bergues B, Aladi M, Földes I, Veisz L 2018 Nat. Commun. 9 4992Google Scholar
[15] Gaumnitz T, Jain A, Pertot Y, Huppert M, Jordan I, Ardana-Lamas F, Wörner H J 2017 Opt. Express 25 27506Google Scholar
[16] Chang Z 2005 Phys. Rev. A 71 023813Google Scholar
[17] Zhao K, Zhang Q, Chini M, Wu Y, Wang X W, Chang Z H 2012 Opt. Lett. 37 3891Google Scholar
[18] Keldysh L V 1965 Sov. Phys. JETP 20 1307
[19] Faisal F H M 1973 J. Phys. B 6 L89Google Scholar
[20] Reiss H R 1980 Phys. Rev. A 22 1786Google Scholar
[21] Lewenstein M, Salières P, L’Huillier A 1995 Phys. Rev. A 52 4747Google Scholar
[22] Ammosov M V, Delone N B, Krainov V 1986 Proc. SPIE 664 1191Google Scholar
计量
- 文章访问数: 9337
- PDF下载量: 100
- 被引次数: 0