-
20世纪80年代, 人们在实验中观测到了最早的高次谐波谱[1,2]. 随后, 高次谐波的产生(HHG)过程就成为强场物理领域的研究热点之一[3-6]. 该过程由激光与物质相互作用时产生的极端非线性效应导致, 其物理机制可用半经典的“三步模型”[7,8]直观地描述, 即电离过程、加速过程和复合过程. 谐波谱中最大的能量通常被称为截止能量, 可以简单地用截止公式
$ 1.32{I_{\text{p}}}{\text{ + 3}}{\text{.17}}{U_{\text{p}}} $ [9]进行估算, 其中Ip和Up分别表示电离能和有质动力势. 虽然关于高次谐波的研究已经从原子分子[7,10,11]拓展到了固体[12-15]和液体[16,17]等体系, 但是原子分子强场下的超快电子动力学研究依然十分重要, 许多超快物理现象还没有得到充分理解, 其对分子结构和分子势场的依赖仍有待深入研究. 其中相对于非极性分子, 极性分子因其在轴向上电子分布的不均匀性(例如图1(a)), 展现出了更加丰富的超快物理现象, 如空间两侧辐射不对称的极性分子阿秒脉冲产生[18]以及各向异性的极性分子轨道层析成像[19]. 特别是极性分子形状共振(shape resonance)的不对称性对于光电离时间延迟有着重要的影响[20]. 因此, 对于极性分子HHG的研究, 有着十分重要的物理意义.图 1 (a) 激光场与极性分子CO的几何关系示意图, 其中激光偏振(红色双箭头)与分子轴(沿Z轴)的夹角为
$ \theta $ ; (b) 当激光偏振沿分子轴方向(即$ \theta = {0^ \circ } $ )时, 在单色场(红色粗实线)的作用下, CO分子每半个激光周期产生的阿秒脉冲(蓝色和绿色细实线)Figure 1. (a) The geometric relation of the laser field and polar molecule CO, where
$ \theta $ is the angle between the molecular (Z) axis and the driving laser polarization (red double arrow); (b) attosecond bursts (blue and green thin solid lines) generated every half cycle of a monochromatic field (red thick line) from polar molecule CO when$ \theta = {0^ \circ } $ .本文研究了极性分子CO在线偏振激光场作用下HHG过程. 第2节简单介绍理论方法和计算细节; 第3节主要从C和O两侧高次谐波辐射的强度、频率范围以及形状共振三个角度, 分析了极性分子HHG过程中的不对称性. 最后, 总结全文的研究内容.
-
本文使用三维含时Hartree-Fock (3D-TDHF)方法来研究CO分子的高次谐波谱. 依托DMTDHF程序包[21]对上述问题展开研究, 该程序包基于Born-Oppenheimer近似下的定核近似, 已经被成功应用于双原子的多电子动力学研究[22-24]. 此外, 鉴于该程序包的开源性, 可以根据自己的个性化需求, 对其加以修改和完善, 进而开发出新功能. 选取CO分子的基态作为含时演化的初始态. 此时分子的总能量为–112.7909 a.u. (1 a.u.= 27.2114 eV), 与文献[25]符合较好. 分子的核间距
$ R $ 取实验测量的平均值2.14 a.u. (1 a.u. = 0.0529 nm), 令分子轴沿Z轴, 且C核和O核分别位于Z轴的$ {{ \pm R} \mathord{\left/ {\vphantom {{ \pm R} 2}} \right. } 2} $ 处(如图1(a)所示). 所使用的激光场可表示为$ E{\text{(}}t{\text{)}} = f{\text{(}}t{\text{)}}{E_0}{\text{sin(}}{\omega _0}t{\text{)}} \text{, } $ 其中激光的角频率
$ {\omega _0} $ = 0.05695 a.u. (1.55 eV), 对应的波长为800 nm;$ {E_0} = 0.0755 $ a.u.为电场的峰值强度, 对应的光强为$ {I_0} = 2 \times {10^{14}} $ W/cm2. 此外, 脉冲包络$ f{\text{(}}t{\text{)}} $ 为Ramped包络, 其中上升区和下降区分别为2个光学周期(T ), 中间的平台区占4个光学周期.根据高次谐波产生的三步模型[7,8]可知, 对于平台区的谐波, 除截止位置外, 其他谐波均有两条路径(长路径和短路径)的电子波包参与高次谐波的辐射. 虽然两条路径下电子波包间的相互干涉对高次谐波有着十分重要的影响, 但短路径的电子波包对实验上测量的谐波数据起主要贡献[26,27]. 因此, 接下来主要关注短路径电子波包贡献的高次谐波谱. 目前有两种常见的方法可用来过滤掉长路径的电子波包, 即空间过滤法和时间过滤法. 对于复杂的多电子分子, 在时域上区分长路径和短路径是非常困难的. 此外, 若所选择的计算空间越小, 则进行3D-TDHF方法所需的计算量越小. 基于以上因素, 本文选择空间过滤法筛选出CO分子短路径电子波包贡献的谐波谱.
在每一步的含时演化结束后, 分子波函数均乘以如下形式的吸收函数,
$ \begin{split}& M(r) = \\ \;&\left\{\begin{aligned} &1, & r\leqslant {r}_{\text{mask}}, \;\;\;\qquad\\ &{\rm{cos}}^{\tfrac{1}{8}}\left(\frac{\pi (r-{r}_{\text{mask}})}{2({r}_{\text{max}}-{r}_{\text{mask}})}\right), & {r}_{\text{mask}} < r\leqslant {r}_{\text{max}}. \end{aligned} \right.\end{split} $ 通过电子轨迹的经典分析, 可以得到短路径电子能够到达的最远距离为L = 26.6 a.u. 在L附近调节吸收参数rmask和rmax, 直至长路径电子波包的贡献被完全吸收. 最后, 通过对加速度形式下的含时偶极矩做傅里叶变换, 然后取其模的平方得到高次谐波谱[10,28], 其中含时偶极矩为包含了所有电子与外场相互作用下的总偶极矩.
-
研究了CO分子在线偏振激光场下产生的高次谐波谱. 当激光偏振沿分子轴方向时, 即
$ \theta = {0^ \circ } $ , CO分子产生的谐波谱及其对应的时频分析谱如图2所示. 仔细观察图2(a), 发现整个谐波谱同时包含奇次和偶次谐波. 由于在沿CO分子轴方向上, 一个光学周期中相邻半周期产生的阿秒脉冲是不对称的(如图1(b)所示), 这一对称性的破缺导致偶次谐波的产生. 谐波谱的截止位置在H36阶附近, 与截止公式估算的数值基本一致. 由于分子的电离能约等于最高占据轨道(HOMO)电子的束缚能, 故此处取$ {I_{\text{p}}} \approx {E_{5\sigma }} = 0.5549{\text{ a}}{\text{.u}}{\text{. }} $ (15.10 eV). 图2(a)中的红色实线为谐波的均线, 结果显示在H12.6阶 (19.5 eV)和H18阶 (27.9 eV)谐波附近有两个明显的峰值. 其中H12.6阶谐波(绿色虚线)峰值能量正好等于$ {I_{\text{p}}} + E \cdot R $ , 代表着从一个核附近电离出去的电子返回到距离较远的另一个核附近; H18阶谐波附近的峰值(玫红色虚线)反映的是CO分子的形状共振[29,30]. 形状共振是一种出现在分子中的常见现象—在光电离过程, 电离的电子在吸引势与排斥势组合作用下暂时被捕获而形成准束缚态的现象[31]. 通常, 形状共振所对应的能级寿命在百阿秒量级.图 2 当激光偏振沿分子轴方向时, CO分子产生的高次谐波谱(a)及其对应的时频分析(b)
Figure 2. High harmonic spectrum (a) and the corresponding time-frequency analysis (b) of CO molecules, when the laser polarization is along the molecular axis.
显然, 通过整个谐波谱很难获得极性分子在激光相邻半周期产生的阿秒脉冲的辐射差异. 为了获取极性分子HHG过程中的时域信息, 利用Morlet小波变化[32]得到了CO分子的时频分析谱(图2(b)).通过时频分析谱可以清晰地观察到长路径电子波包的吸收效果. 当吸收参数取值为
$ {r_{{\text{mask}}}} = 17{\text{ a}}{\text{.u}}{\text{. }} $ 和$ {r_{{\text{max}}}} = 40{\text{ a}}{\text{.u}}. $ 时, 长路径电子波包对高次谐波的贡献基本被完全吸收. 可以看到相邻半周期的阿秒脉冲辐射也存在着明显的差异, 例如, (3.5T, 4.0T)周期和(4.0T, 4.5T)周期. 通常情况下, 对于复杂的多电子体系, 多电子效应对于体系的电子动力学将产生重要的影响[11,19,23,33]. 显然, 由于相邻半个周期的阿秒脉冲均受到内壳层电子的影响, 因此, 通过多电子理论来分析相邻半周期的辐射差异十分困难. 为简化问题, 以隧穿电离时刻出射电子的方向来对电子进行区分. 当电场由C指向O时, 电离电子的出射方向为由O指向C, 此时我们标记这个电子为C侧电离的电子. 通常, 从C侧电离出去的电子也会从C侧返回到母离子附近, 与母离子复合并释放出高次谐波. 同理, 当电场反向的时候, 标记电子为从O侧出射的电子.将上述对电离电子的分类方法同所用的激光场(图2(b)中红色虚线所示)相结合, 可以得到不同时刻辐射的阿秒脉冲是由哪侧电离的电子贡献而来. 换句话说, 我们可获得从不同方向返回的电子所辐射的阿秒脉冲在时域上的分布情况. 例如, (3.5T, 4.0T)半周期内, 电场
$ E(t) $ 为正值, 电场的方向由O指向C, 此时由C侧电离出去的电子在电场的作用下被拉回母离子附近, 与母离子复合并释放出阿秒脉冲, 故此时间段内辐射的阿秒脉冲主要由C侧出射的电子贡献而来. 同理, (4.0T, 4.5T)半周期内所辐射的阿秒脉冲主要来自于O侧电离电子的贡献. 显然C和O两侧辐射的阿秒脉冲存在明显的差异. 一方面, 在(3.5T, 4.0T)半周期内, 由C侧电子贡献的阿秒脉冲强度明显要高于(4.0T, 4.5T)半周期由O侧电子贡献的结果. 这是由于C和O两侧的电离率存在的巨大差异导致[22].另一方面, 来自C侧的谐波辐射的频率范围也要明显宽于O侧的结果. 此处, 对于平台区的高次谐波, H7阶到H17阶的谐波主要由C侧电离的电子贡献产生, 其余的谐波(即从H18阶到H36阶)则是由C和O两端的电子共同参与的结果. 两端电离的电子对不同阶谐波的贡献差异体现了多电子动力学的复杂性. 最后, 对于图2(a)谐波谱中出现的两个峰值, 在时频分析结果中也有明显的体现. 对于H12.6阶谐波所在的峰值, 此时可以明确是由C侧电离的电子复合到了较远的原子核O附近形成的; 而H18阶的谐波峰值则是由C和O两侧的能量相近的形状共振所导致.为了进一步明确C和O两侧出射电子对谐波辐射强度的贡献差异, 我们研究了不同取向角
$ \theta $ 下CO分子产生的高次谐波谱及其对应的时频分析谱, 如图3所示. 其中, 左栏为高次谐波谱, 右栏为对应的时频分析结果. 注: 除$ \theta $ 外, 其余的计算参数同图2一致. 当$ \theta \ne {0^ \circ } $ 时, 谐波辐射并不是线偏振的. 本文重点关注谐波辐射的主要分量, 即谐波沿激光偏振方向的分量. 已知电子的超快动力学受到分子取向的强烈调制, 故对于平台区的谐波而言, 来自两端的阿秒辐射同样受到$ \theta $ 的重要影响. 连同图2(b)一起观察, 可以发现当改变$ \theta $ 时, 来自两端的谐波强度和辐射贡献范围, 均发生了明显的变化.当$ \theta = {0^ \circ } $ 时, 相邻半周期的差异达到最大; 随着取向角$ \theta $ 的增加, 差异越来越小; 直至$ \theta = {90^ \circ } $ 时, 两端差异完全消失. 通过时频分析结果可知, 相邻半周期的阿秒脉冲完全一样, 与谐波谱中仅包含奇次谐波相互验证, 如图3(c)所示. 这一结果反映出沿激光偏振方向的轨道波函数完全对称导致来自两端的辐射差异完全消失.图 3 不同取向角
$ \theta $ 下, CO分子产生的高次谐波谱及其对应的时频分析 (a)$ \theta = {30^ \circ } $ ; (b)$ \theta = {60^ \circ } $ ; (c)$ \theta = {90^ \circ } $ . 除$ \theta $ 外, 其他计算参数同图2保持一致, 其中$ \theta $ 为分子轴与激光偏振的夹角, 如图1(a)所示Figure 3. High harmonic spectra and corresponding time-frequency analyses at (a)
$ \theta = {30^ \circ } $ $ \theta $ , (b)$ \theta = {60^ \circ } $ and (c)$ \theta = {90^ \circ } $ . Note all parameters remain the same with Fig. 2; except for$ \theta $ , where$ \theta $ is the angle between the molecular axis and the driving laser polarization as shown in Fig.1 (a). -
为了验证3.1节的结论, 下面通过谐波强度对极性分子的不对称性进行数值分析. 已知分子总的谐波谱是由C和O两侧辐射的阿秒脉冲相互干涉而来的, 为获取C和O两侧的辐射强度差异, 应分别得到两侧的谐波辐射强度. 首先, 根据谐波的时频分析谱, 可以获取任意一端谐波辐射的时间范围. 然后, 利用这个时间信息, 调整合适的时间窗[34], 再同含时偶极矩相结合便可以滤出仅来自一端的谐波谱. 用IC(或IO)来标记仅来自C端(或O端)电离出去的电子返回母离子后辐射的谐波强度, 其中对谐波的积分范围为从H7阶到H50阶. 最终, 计算来自C和O两端的谐波强度比值IC/IO, 结果如图4所示. 整体来看, IC/IO随着取向角的增加, 越来越小, 表示极性分子的极化特性随着取向角的增加越来越不明显. 在θ = 0º附近, 分子的极性最大. 当θ = 90º时, 该比值为1, 分子的极性消失. 这与之前的预测一致.
图 4 IC/IO随着取向角θ的变化关系, 其中IC(或IO)代表从C侧(或O侧)电离电子所辐射的谐波强度. 谐波强度的积分范围为H7阶到H50阶
Figure 4. IC/IO as a function of the
$ \theta $ . IC (or IO) denotes the harmonic intensity emitted by the electrons ionized from C (or O) side with an integration range from H7th to H50th上面关于极性分子差异的理论分析, 在实验上并不能实现. 如何在实验上获取分子的极性信息呢? 在沿激光偏振的方向上, 如果轨道波函数是对称的, 那么谐波谱中仅包含奇次谐波. 当对称性被破坏时, 偶次谐波就出现了, 因此奇偶次谐波的相对强度在一定程度上可以反映出分子体系的对称性. 因此, 下面研究平台区偶次和奇次谐波的强度比I2n/I2n+1随着取向角θ的变化关系, 如图5所示. 结果显示, 当θ = 0º时, I2n/I2n+1达到最大值; 当θ = 90º时, I2n/I2n+1达到最小值. 整体来看, 平台区的偶次和奇次谐波的强度比I2n/I2n+1随着取向角θ的增加呈递减的趋势. 图4和图5均表明, CO分子的极性差异随着激光偏振与分子轴的夹角的变化存在一定的角分布结构.
图 5 平台区的偶次和奇次谐波的强度比I2n/I2n+1随着取向角θ的变化关系, 其中I2n(或I2n+1)表示对H7到H50阶谐波中偶(或奇)次谐波的强度积分结果
Figure 5. The intensity ratio of even and odd order harmonics I2n/I2n+1 versus the alignment angle θ, where I2n (or I2n+1) represents the integrated intensity of even (or odd) order harmonics with a frequency range from H7th to H50th.
-
3.1节对于图2的描述中提到, 谐波谱第H18阶附近有个明显的峰值, 表示形状共振的位置. 从时频分析结果中明显看到沿C和O不同方向上的形状共振强度存在明显差异. 下面通过强场近似(SFA)理论来分析形状共振的不对称性.
$ S(\omega ) \propto N{\omega ^4}{\left| {a\left[ {k(\omega )} \right]d(\omega )} \right|^2} , $ 其中, N为电离速率;
$ a\left[ {k(\omega )} \right] $ 表示返回电子波包的振幅;$ d(\omega ) = \left\langle {{\varPsi _0}\left| r \right|k(\omega )} \right\rangle $ 表示基态$ {\varPsi _0} $ 到连续态$ k(\omega ) = \sqrt {2\omega } $ 的跃迁矩阵元. 光电离截面表示使原子分子中占据电子的壳层产生空穴的概率. 光电离截面$ \sigma (\omega ) $ 与跃迁矩阵元$ d(\omega ) $ 满足如下关系式:$ \sigma (\omega ) \propto {d^2}(\omega ) . $ 对于, 任意一侧(C或O)的谐波辐射均满足(3)式, 故有
$ {S_{{\text{C/O}}}}(\omega ) \propto {N_{{\text{C/O}}}}{\omega ^4}{\left| {a\left[ {k(\omega )} \right]{d_{{\text{C/O}}}}(\omega )} \right|^2} . $ 根据上述表达式, 可以得到来自C侧辐射的谐波强度与来自O侧的比值:
$ \frac{{{S_{\text{C}}}(\omega )}}{{{S_{\text{O}}}(\omega )}} = \frac{{{N_{\text{C}}}}}{{{N_{\text{O}}}}}\frac{{{\sigma _{\text{C}}}(\omega )}}{{{\sigma _{\text{O}}}(\omega )}} . $ 因此, 如果已知C和O两侧的谐波强度比和电离比,就可以得到光电离截面的比值
$ \dfrac{{{\sigma _{\text{C}}}(\omega )}}{{{\sigma _{\text{O}}}(\omega )}} $ . 下面通过3种不同的处理方法来计算截面比值, 并进行对比. 方法一, 我们在仅考虑HOMO轨道的前提下, 使用ePloyScat程序[37,38]计算所得的CO分子的光电离截面, 结果如图6所示. 位于27.5 eV附近的形状共振所对应的光电离截面比值为$ \dfrac{{{\sigma _{\text{C}}}(\omega )}}{{{\sigma _{\text{O}}}(\omega )}} = 3.02 $ . 方法二, 利用前面提到的通过选取恰当的时间窗函数可挑选出C侧或O侧单独贡献的谐波谱. 当θ = 0º时, 计算得到H18阶C和O两侧的谐波辐射强度比为$ \dfrac{{{S_{\text{C}}}(\omega )}}{{{S_{\text{O}}}(\omega )}} = 15.92 $ . 此外, 电离速率随着时间的变化关系如图7所示. 图7中黑色实线为归一化的电子占据数, 通过求解电子占据数对时间的导数可得分子的电离速率(绿色实线). 蓝色箭头所示区域为激光的平台区所对应的电离速率区间. 在同一个光学周期内, 电离速率的幅值较大的半周期对应为C侧电离的时间区间, 反之则为O侧电离的区间. 然后分别对处于平台区的C和O两侧的电离率求均值, 可得两侧电离率的比值为$ \dfrac{{{N_{\text{C}}}}}{{{N_{\text{O}}}}} = 2.71 $ . 利用(6)式, 最终可得$ \dfrac{{{\sigma _{\text{C}}}(\omega )}}{{{\sigma _{\text{O}}}(\omega )}} = 5.87 $ . 方法三, 由于直接在实验上分离出C或O侧的谐波贡献十分困难, 因此亟需另辟蹊径重新得到两侧的谐波强度比.图 6 沿O(黑色实线)和C(红色虚线)两侧的最外层轨道的光电离截面, 该结果由ePloyScat程序计算所得
Figure 6. Photoionization cross sections of the HOMO from O side (black solid line) and C side (red dotted line), which are simulated by ePloyScat.
图 7 归一化的电子占据数(黑色实线)和电离速率(绿色实线)随着时间的变化关系. 红色虚线为激光场, 蓝色箭头所示的区域为激光脉冲的平台区
Figure 7. Time profile of normalized population (black solid line) and ionization rates (green solid line) of electrons. The red dotted line denotes the used laser field. The time range marked by blue arrows is the plateau region of the used laser field.
已知, C和O两侧辐射的谐波强度同总的谐波强度之间满足如下关系:
$ \begin{split}&\frac{{{{\left| {\sqrt {{S_{\text{C}}}} + \sqrt {{S_{\text{O}}}} } \right|}^2}}}{{{{\left| {\sqrt {S{}_{\text{C}}} - \sqrt {{S_{\text{O}}}} } \right|}^2}}}\\ =\;& \frac{{2S(2n{\omega _0})}}{{S[(2n - 1){\omega _0}] + S[(2n + 1){\omega _0}]}} = {\lambda ^2} \text{, }\end{split} $ 其中对于总谐波谱的某阶谐波强度, 在实验上是可以获取的, 故(7)式中λ就会变成已知量. 对于H18阶处的形状共振位置的分析, 利用
$ \dfrac{{2 S({\text{H}}18)}}{{S[{\text{H}}17] + S[{\text{H}}19]}} $ , 可得$ \lambda = 1.73 $ . 那么, 通过$ \dfrac{{{{\left| {\sqrt {{S_{\text{C}}}} + \sqrt {{S_{\text{O}}}} } \right|}^2}}}{{{{\left| {\sqrt {{S_{\text{C}}}} - \sqrt {{S_{\text{O}}}} } \right|}^2}}} = {\lambda ^2} $ , 可得$\dfrac{{{S_{\text{C}}}}}{{{S_{\text{O}}}}} = {\left( {\dfrac{{1 + \lambda }}{{\lambda - 1}}} \right)^2} = 14.59$ . 再结合方法二中得到的两侧电离率的比值$\dfrac{{{N_{\text{C}}}}}{{{N_{\text{O}}}}}$ , 最终可得$ \dfrac{{{\sigma _{\text{C}}}(\omega )}}{{{\sigma _{\text{O}}}(\omega )}} = 5.38 $ .上面通过两种方法(方法二和方法三)对TDHF结果进行了分析, 所得的截面比值十分接近, 但与基于单轨道假设的ePloyScat的计算结果存在较大差异. 这反映出内壳层电子的影响增大了极性分子谐波辐射的不对称性, 因此多电子效应对于极性分子的高次谐波辐射有着十分重要的影响.
-
本文使用三维含时Hartree-Fock方法研究了极性分子CO在强激光场的作用下产生高次谐波的过程. 当激光偏振沿分子轴方向时, CO分子的谐波辐射在相邻半周期间存在明显的不对称性. 对于平台区的低阶谐波, 主要是C侧的贡献, 而对于高阶谐波, 则是由C和O两侧共同参与而来. 并且, C侧与O侧的谐波辐射强度比值随着分子取向角
$ \theta $ 的增加逐渐变小. 此外, C和O两侧的形状共振也存在明显的差异, 通过强场近似理论得到C和O两侧的光电离截面比值在5.5附近, 大于只考虑HOMO轨道下应用ePloyScat程序计算所得的比值3. 由此可见多电子效应增强了极性分子的不对称性. 本文的研究工作探索了高次谐波在超快电子动力学方面的应用, 加深了对极性分子的高次谐波辐射过程中不对称性的理解, 为后续探究复杂体系的多电子动力学问题提供了有益的启示.
-
极性分子由于其本身的不对称性, 在强激光场的作用下展现出丰富且复杂的电子动力学现象. 本文利用三维含时Hartree-Fock方法研究了极性分子CO的高次谐波产生过程. 通过高次谐波谱和时频分析结果可知,当激光偏振沿分子轴方向时, 来自C和O两侧的电子对高次谐波的产生具有不同的贡献. 对于平台区较低阶的谐波, 仅C侧电离的电子参与谐波的产生. 而对于较高阶的谐波, C和O两侧的电子共同参与谐波的辐射. 并且, 随着激光偏振与分子轴的夹角θ逐渐增大, C和O两侧电子对谐波强度贡献的差异越来越小. 在高次谐波谱中能量28 eV附近发现了明显的形状共振峰, 随后通过强场近似理论解析了其不对称性. 本文的工作有助于推动高次谐波谱在追踪电子超快动力学和研究极性分子结构方面的应用.
-
关键词:
- 高次谐波 /
- 极性分子 /
- 含时Hartree-Fock /
- 形状共振
Compared with nonpolar molecules, owing to the inherent asymmetry, polar molecules exhibit rich and very complex electronic dynamics under the interaction with strong laser fields. In this work, high-order harmonic generation (HHG) of polar molecules CO is investigated by using the three-dimensional time-dependent Hartree-Fock (3D-TDHF) theory, with all electrons active. Through the high harmonic spectra and time-frequency analyses, it is found that when the laser field polarizes along the molecular axis, the ionized electrons from the two sides (C side and O side) contribute differently to the harmonic radiation. On the one hand, the harmonic intensity from the C side is greater than that from the O side, which is caused by the ionization rate. On the other hand, for the lower-order (7th–17th order) harmonics of plateau region, only the electrons from the C side participate in the HHG. However, for its higher part (18th–36th order), the electrons from both C side and O side contribute to high harmonics simultaneously. Moreover, the difference between contributions from two sides is related to the alignment angle θ between the laser polarization and the molecular axis, and it reaches a maximum value around θ = 0º and a minimum value around θ = 90º. There are two strong resonances around harmonic order H12.6 (19.5 eV) and H18 (27.9 eV) in the harmonic spectra when θ = 0º. The first resonance around H12.6 reveals that part of electrons ionized from the C side recombine to the vicinity of the further O nucleus. Near the second resonance around H18, there appears a shape resonance. Nevertheless, the shape resonances from the C and O sides are disparate. Based on the strong-field approximation theory, the ratio between photoionization cross sections from C and O sides around the shape resonance is calculated. The ratio is about 5.5 from 3D-TDHF, which is greater than the result of 3 simulated by ePloyScat, where only HOMO is considered. This discrepancy reveals that multi-electron effects enhance the asymmetry of polar molecules. This work provides an in-depth insight into the asymmetry in HHG of polar molecules, which benefits the generation of isolated attosecond pulse . It also promotes the application of high harmonic spectra in tracking the ultrafast dynamics of electrons.-
Keywords:
- high-order harmonic generation /
- polar molecule /
- time-dependent Hartree–Fock /
- shape resonance
[1] McPherson A, Gibson G, Jara H, Johann U, Luk T S, McIntyre I A, Boyer K, Rhodes C K 1987 J. Opt. Soc. Am. B: Opt. Phys. 4 595
Google Scholar
[2] Ferray M, L'Huillier A, Li X F, Lompre L A, Mainfray G, Manus C 1988 J. Phys. B: At. Mol. Opt. Phys. 21 L31
Google Scholar
[3] Brabec T, Krausz F 2000 Rev. Mod. Phys. 72 545
Google Scholar
[4] Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163
Google Scholar
[5] Pazourek R, Nagele S, Burgdörfer J 2015 Rev. Mod. Phys. 87 765
Google Scholar
[6] Mourou G 2019 Rev. Mod. Phys. 91 030501
Google Scholar
[7] Krause J L, Schafer K J, Kulander K C 1992 Phys. Rev. Lett. 68 3535
Google Scholar
[8] Corkum P B 1993 Phys. Rev. Lett. 71 1994
Google Scholar
[9] Lewenstein M, Balcou P, Ivanov M Y, L’Huillier A, Corkum P B 1994 Phys. Lett. A 49 2117
Google Scholar
[10] Yang Y, Liu L, Zhao J, Tu Y, Liu J, Zhao Z 2021 J. Phys. B:At. Mol. Opt. Phys. 54 144009
Google Scholar
[11] Huang Y, Zhao J, Shu Z, Zhu Y, Liu J, Dong W, Wang X, Lü Z, Zhang D, Yuan J, Chen J, Zhao Z 2021 Ultrafast Science 2021 9837107
Google Scholar
[12] Li L, Zhang Y, Lan P, Huang T, Zhu X, Zhai C, Yang K, He L, Zhang Q, Cao W, Lu P 2021 Phys. Rev. Lett. 126 187401
Google Scholar
[13] Ghimire S, Reis D A 2019 Nat. Phys. 15 10
Google Scholar
[14] Liu L, Zhao J, Dong W, Liu J, Huang Y, Zhao Z 2017 Phys. Lett. A 96 053403
Google Scholar
[15] Liu L, Zhao J, Yuan J M, Zhao Z X 2019 Chin. Phys. B 28 114205
Google Scholar
[16] Luu T T, Yin Z, Jain A, Gaumnitz T, Pertot Y, Ma J, Wörner H J 2018 Nat. Commun. 9 3723
Google Scholar
[17] Zeng A W, Bian X B 2020 Phys. Rev. Lett. 124 203901
Google Scholar
[18] Hu H, Li N, Liu P, Li R, Xu Z 2017 Phys. Rev. Lett. 119 173201
Google Scholar
[19] Ren Z, Yang Y, Zhu Y, Zan X, Zhao J, Zhao Z 2021 J. Phys. B: At. Mol. Opt. Phys. 54 185601
Google Scholar
[20] Vos J, Cattaneo L, Patchkovskii S, Zimmermann T, Cirelli C, Lucchini M, Kheifets A, Landsman A S, Keller U 2018 Science 360 1326
Google Scholar
[21] Zhang B, Yuan J, Zhao Z 2015 Comput. Phys. Commun. 194 84
Google Scholar
[22] Zhang B, Yuan J, Zhao Z 2013 Phys. Rev. Lett. 111 163001
Google Scholar
[23] Zhang B, Yuan J, Zhao Z 2014 Phys. Rev. A 90 035402
Google Scholar
[24] Zhang B, Zhao J, Zhao Z 2018 Chin. Phys. Lett. 35 043201
Google Scholar
[25] Kobus J 1993 Chem. Phys. Lett. 202 7
Google Scholar
[26] Dong W, Hu H, Zhao Z 2020 Opt. Express 28 22490
Google Scholar
[27] Zhang B, Lein M 2019 Phys. Rev. A 100 043401
Google Scholar
[28] Le C T, Vu D D, Ngo C, Le V H 2019 Phys. Rev. A 100 053418
Google Scholar
[29] Kraus P M, Baykusheva D, Wörner H J 2014 Phys. Rev. Lett. 113 023001
Google Scholar
[30] Frumker E, Kajumba N, Bertrand J B, Wörner H J, Hebeisen C T, Hockett P, Spanner M, Patchkovskii S, Paulus G G, Villeneuve D M, Naumov A, Corkum P B 2012 Phys. Rev. Lett. 109 233904
Google Scholar
[31] 徐克尊 2012 高等原子分子物理学 第3版 (合肥: 中国科学技术大学出版社) 第51—52页
Xu K Z 2012 Advanced Atomic and Molecular Physics (3rd Ed.) (Hefei: University of Science and Technology of China Press) pp51—52 (in Chinese)
[32] Antoine P, Piraux B, Maquet A 1995 Phys. Lett. A 51 R1750
[33] Shu Z, Liang H, Wang Y, Hu S, Chen S, Xu H, Ma R, Ding D, Chen J 2022 Phys. Rev. Lett. 128 183202
Google Scholar
[34] 刘璐, 赵晶, 袁建民, 赵增秀 2017 中国科学: 物理学 力学 天文学 047 033006
Google Scholar
Liu L, Zhao J, Yuan J M, Zhao Z X 2017 Sci Sin-Phys. Mech. Astron. 047 033006
Google Scholar
[35] Chen Y J, Fu L B, Liu J 2013 Phys. Rev. Lett. 111 073902
Google Scholar
[36] Huang Y, Meng C, Wang X, Lü Z, Zhang D, Chen W, Zhao J, Yuan J, Zhao Z 2015 Phys. Rev. Lett. 115 123002
Google Scholar
[37] Gianturco F A, Lucchese R R, Sanna N 1994 J. Chem. Phys. 100 6464
Google Scholar
[38] Natalense A P P, Lucchese R R 1999 J. Chem. Phys. 111 5344
Google Scholar
-
图 1 (a) 激光场与极性分子CO的几何关系示意图, 其中激光偏振(红色双箭头)与分子轴(沿Z轴)的夹角为
$ \theta $ ; (b) 当激光偏振沿分子轴方向(即$ \theta = {0^ \circ } $ )时, 在单色场(红色粗实线)的作用下, CO分子每半个激光周期产生的阿秒脉冲(蓝色和绿色细实线)Fig. 1. (a) The geometric relation of the laser field and polar molecule CO, where
$ \theta $ is the angle between the molecular (Z) axis and the driving laser polarization (red double arrow); (b) attosecond bursts (blue and green thin solid lines) generated every half cycle of a monochromatic field (red thick line) from polar molecule CO when$ \theta = {0^ \circ } $ .图 3 不同取向角
$ \theta $ 下, CO分子产生的高次谐波谱及其对应的时频分析 (a)$ \theta = {30^ \circ } $ ; (b)$ \theta = {60^ \circ } $ ; (c)$ \theta = {90^ \circ } $ . 除$ \theta $ 外, 其他计算参数同图2保持一致, 其中$ \theta $ 为分子轴与激光偏振的夹角, 如图1(a)所示Fig. 3. High harmonic spectra and corresponding time-frequency analyses at (a)
$ \theta = {30^ \circ } $ $ \theta $ , (b)$ \theta = {60^ \circ } $ and (c)$ \theta = {90^ \circ } $ . Note all parameters remain the same with Fig. 2; except for$ \theta $ , where$ \theta $ is the angle between the molecular axis and the driving laser polarization as shown in Fig.1 (a).图 5 平台区的偶次和奇次谐波的强度比I2n/I2n+1随着取向角θ的变化关系, 其中I2n(或I2n+1)表示对H7到H50阶谐波中偶(或奇)次谐波的强度积分结果
Fig. 5. The intensity ratio of even and odd order harmonics I2n/I2n+1 versus the alignment angle θ, where I2n (or I2n+1) represents the integrated intensity of even (or odd) order harmonics with a frequency range from H7th to H50th.
图 7 归一化的电子占据数(黑色实线)和电离速率(绿色实线)随着时间的变化关系. 红色虚线为激光场, 蓝色箭头所示的区域为激光脉冲的平台区
Fig. 7. Time profile of normalized population (black solid line) and ionization rates (green solid line) of electrons. The red dotted line denotes the used laser field. The time range marked by blue arrows is the plateau region of the used laser field.
-
[1] McPherson A, Gibson G, Jara H, Johann U, Luk T S, McIntyre I A, Boyer K, Rhodes C K 1987 J. Opt. Soc. Am. B: Opt. Phys. 4 595
Google Scholar
[2] Ferray M, L'Huillier A, Li X F, Lompre L A, Mainfray G, Manus C 1988 J. Phys. B: At. Mol. Opt. Phys. 21 L31
Google Scholar
[3] Brabec T, Krausz F 2000 Rev. Mod. Phys. 72 545
Google Scholar
[4] Krausz F, Ivanov M 2009 Rev. Mod. Phys. 81 163
Google Scholar
[5] Pazourek R, Nagele S, Burgdörfer J 2015 Rev. Mod. Phys. 87 765
Google Scholar
[6] Mourou G 2019 Rev. Mod. Phys. 91 030501
Google Scholar
[7] Krause J L, Schafer K J, Kulander K C 1992 Phys. Rev. Lett. 68 3535
Google Scholar
[8] Corkum P B 1993 Phys. Rev. Lett. 71 1994
Google Scholar
[9] Lewenstein M, Balcou P, Ivanov M Y, L’Huillier A, Corkum P B 1994 Phys. Lett. A 49 2117
Google Scholar
[10] Yang Y, Liu L, Zhao J, Tu Y, Liu J, Zhao Z 2021 J. Phys. B:At. Mol. Opt. Phys. 54 144009
Google Scholar
[11] Huang Y, Zhao J, Shu Z, Zhu Y, Liu J, Dong W, Wang X, Lü Z, Zhang D, Yuan J, Chen J, Zhao Z 2021 Ultrafast Science 2021 9837107
Google Scholar
[12] Li L, Zhang Y, Lan P, Huang T, Zhu X, Zhai C, Yang K, He L, Zhang Q, Cao W, Lu P 2021 Phys. Rev. Lett. 126 187401
Google Scholar
[13] Ghimire S, Reis D A 2019 Nat. Phys. 15 10
Google Scholar
[14] Liu L, Zhao J, Dong W, Liu J, Huang Y, Zhao Z 2017 Phys. Lett. A 96 053403
Google Scholar
[15] Liu L, Zhao J, Yuan J M, Zhao Z X 2019 Chin. Phys. B 28 114205
Google Scholar
[16] Luu T T, Yin Z, Jain A, Gaumnitz T, Pertot Y, Ma J, Wörner H J 2018 Nat. Commun. 9 3723
Google Scholar
[17] Zeng A W, Bian X B 2020 Phys. Rev. Lett. 124 203901
Google Scholar
[18] Hu H, Li N, Liu P, Li R, Xu Z 2017 Phys. Rev. Lett. 119 173201
Google Scholar
[19] Ren Z, Yang Y, Zhu Y, Zan X, Zhao J, Zhao Z 2021 J. Phys. B: At. Mol. Opt. Phys. 54 185601
Google Scholar
[20] Vos J, Cattaneo L, Patchkovskii S, Zimmermann T, Cirelli C, Lucchini M, Kheifets A, Landsman A S, Keller U 2018 Science 360 1326
Google Scholar
[21] Zhang B, Yuan J, Zhao Z 2015 Comput. Phys. Commun. 194 84
Google Scholar
[22] Zhang B, Yuan J, Zhao Z 2013 Phys. Rev. Lett. 111 163001
Google Scholar
[23] Zhang B, Yuan J, Zhao Z 2014 Phys. Rev. A 90 035402
Google Scholar
[24] Zhang B, Zhao J, Zhao Z 2018 Chin. Phys. Lett. 35 043201
Google Scholar
[25] Kobus J 1993 Chem. Phys. Lett. 202 7
Google Scholar
[26] Dong W, Hu H, Zhao Z 2020 Opt. Express 28 22490
Google Scholar
[27] Zhang B, Lein M 2019 Phys. Rev. A 100 043401
Google Scholar
[28] Le C T, Vu D D, Ngo C, Le V H 2019 Phys. Rev. A 100 053418
Google Scholar
[29] Kraus P M, Baykusheva D, Wörner H J 2014 Phys. Rev. Lett. 113 023001
Google Scholar
[30] Frumker E, Kajumba N, Bertrand J B, Wörner H J, Hebeisen C T, Hockett P, Spanner M, Patchkovskii S, Paulus G G, Villeneuve D M, Naumov A, Corkum P B 2012 Phys. Rev. Lett. 109 233904
Google Scholar
[31] 徐克尊 2012 高等原子分子物理学 第3版 (合肥: 中国科学技术大学出版社) 第51—52页
Xu K Z 2012 Advanced Atomic and Molecular Physics (3rd Ed.) (Hefei: University of Science and Technology of China Press) pp51—52 (in Chinese)
[32] Antoine P, Piraux B, Maquet A 1995 Phys. Lett. A 51 R1750
[33] Shu Z, Liang H, Wang Y, Hu S, Chen S, Xu H, Ma R, Ding D, Chen J 2022 Phys. Rev. Lett. 128 183202
Google Scholar
[34] 刘璐, 赵晶, 袁建民, 赵增秀 2017 中国科学: 物理学 力学 天文学 047 033006
Google Scholar
Liu L, Zhao J, Yuan J M, Zhao Z X 2017 Sci Sin-Phys. Mech. Astron. 047 033006
Google Scholar
[35] Chen Y J, Fu L B, Liu J 2013 Phys. Rev. Lett. 111 073902
Google Scholar
[36] Huang Y, Meng C, Wang X, Lü Z, Zhang D, Chen W, Zhao J, Yuan J, Zhao Z 2015 Phys. Rev. Lett. 115 123002
Google Scholar
[37] Gianturco F A, Lucchese R R, Sanna N 1994 J. Chem. Phys. 100 6464
Google Scholar
[38] Natalense A P P, Lucchese R R 1999 J. Chem. Phys. 111 5344
Google Scholar
计量
- 文章访问数: 627
- PDF下载量: 63
- 被引次数: 0