优化组合激光场驱动原子产生高次谐波及单个超短阿秒脉冲理论研究*

汉琳 苗淑莉 李鹏程†

1) (汕头大学理学院物理系,先进光学与光子学研究中心,汕头 515063)

2) (兰州大学物理科学与技术学院,兰州理论物理中心,甘肃省理论物理重点实验室,兰州 730000)

高次谐波是获得阿秒紫外光源最主要的方法之一,是强场超快领域研究的热点问题,具有非常广泛的应用前景.本文围绕如何产生超连续高次谐波平台及单个超短阿秒脉冲面临的问题,概述了这方面研究的进展,并从理论上展示了一种有效可行的方案,即将强激光场中的含时薛定谔方程与非约束优化算法相结合,以扩展谐波平台最宽为目标函数,分别优化双色和三色组合激光场并驱动氦原子产生超连续高次谐波谱.优化后的双色组合激光场驱动氦原子产生的超连续谐波谱平台达到了100 阶,叠加获得了最短25 as 的单个阿秒脉冲;优化后的三色组合激光场驱动氦原子产生的超连续谐波谱平台宽度达到了170 阶,叠加获得最短17 as的单个阿秒脉冲,同时谐波转换效率也有所提高.为了给实验提供切实可行的参考,本文以优化的双色组合激光场情况为例,基于同时求解含时薛定谔方程和麦克斯韦方程,进一步考虑了介质宏观演化效应对单原子层次产生阿秒脉冲的影响,发现利用远场轴外量子通道的空间选择性可以获得更短的单个阿秒脉冲.

1 引言

激光与物质相互作用可以产生许多非线性光学现象,如多光子电离[1–7](multiphoton ionization,MPI)、阈上电离[8–12](above-threshold ionization,ATI)、隧穿电离[13–19](tunneling ionization,TI)、高 次谐波[20–24](high-order harmonic generation,HHG)等,其中高次谐波辐射是产生阿秒相干远紫外辐射源最有前途的方法之一.阿秒脉冲可以实时探测物质中电子动力学过程,使得研究物质结构的层次从分子扩展到原子内部,为理解和控制微观世界的动力学过程提供了强有力的工具[25].

几十年来,高次谐波及阿秒脉冲产生的研究引起了人们极大的兴趣.1987 年,McPherson 等[26]利用波长248 nm 的KrF 激光与稀有气体Ne 原子相互作用,首次在实验上观测到了高达17 阶的高次谐波辐射.1988 年,Ferray 等[27]用1064 nm的Nd:YAG 激光与稀有气体Xe 原子相互作用,观察到了高达33 阶的高次谐波.1993 年,L’Huilier和Balcou[28]用1053 nm 的Nd:Glass 激光分别与稀有气体Xe,Ar,Ne 原子相互作用,对应观测到了29 阶、57 阶以及135 阶高次谐波辐射.同年,Macklin 等[29]利用800 nm 的红外激光与Ne 原子气体相互作用,观察到了109 阶高次谐波.高次谐波特征谱拥有独特的频谱结构[30],低阶谐波随着谐波阶次增加呈指数型衰减,紧接着出现了一个平台区,谐波强度随谐波阶次的增加无明显变化,最后在特定阶次出现一个谐波强度快速下降的截止.

Corkum[31]和Kulander 等[32]分别提出了半经典的“三步模型”,很好地解释了强场高次谐波的产生过程.首先,库仑场在外加激光场作用下被压低形成势垒,电子可以隧穿通过势垒,紧接着在激光场作用下加速获得能量,当激光场反向时电子开始减速,部分电子反向加速回到母核附近与母核碰撞复合,释放出多余能量,并辐射出高次谐波.电子复合辐射出的光子能量是电子动能和电离势能的总和,最大瞬时动能约等于 3.17Up,Up为有质动力能,所以高次谐波截止位置的能量[31]ECutoff=Ip+3.17Up,其中Ip代表原子的电离能.Lewenstein等[33,34]基于量子理论提出了一种强场近似模型(strong-field approximation,SFA),该模型计算简单,可以解释大多数强激光场与原子分子相互作用产生的非线性光学现象,但是该方法没有考虑库仑势的影响,并且忽略了原子激发态对谐波的贡献,无法精确描述原子分子在强激光场中强烈的非线性动力学过程.数值求解含时薛定谔方程(timedependent Schrödinger equation,TDSE)可以避免强场近似理论的不足,但是求解过程计算量大,非常耗时.Kulander[35]和Krause 等[36]提出了简化的单电子(single-active-electron,SAE)模型理论,极大地减少了计算量,可以基于TDSE 更精确地模拟强激光场与原子分子相互作用产生的非线性光学现象.对于多电子体系,以最简单的两电子He 原子为例,完全求解TDSE 需要在超级计算机上进行,很难实现.1984 年,Runge 和Gross[37]发展了含时密度泛函理论,用含时的电子密度代替含时电子波函数,体系的维数不再随电子数目增加,大大减少了计算量,使得处理强激光场中复杂多电子体系的模拟成为可能.

当前,高次谐波的产生研究主要集中在两个方面: 一是增强高次谐波的发射效率;二是扩展高次谐波的平台.高次谐波的发射效率受电子复合概率和空间传播相位匹配等因素的影响[38].高次谐波平台截止频率正比于Iλ2,其中I代表激光场强度,λ代表激光场波长.通过增加驱动激光场的强度,可以增加电子电离概率,实现电子复合概率的增加,也可以达到扩展高次谐波平台的目的.1997年,Spielmann 等[39]和Chang 等[40]利用钛蓝宝石激光器产生中心波长为800 nm 的短波长高强度激光,并驱动He 原子产生了接近500 eV 的高次谐波谱.相比其他中性稀有气体原子,He 原子的饱和强度最高,可承受的激光强度最高,但是He 原子产生的高次谐波效率比其他气体原子低几个数量级[41].由于受限于电子电离饱和效应,不可能无限制地增加驱动激光的强度来扩展谐波平台,因此增加驱动激光的波长就成为另一个必然的选择[42–53].虽然利用长波长驱动可以增大谐波截止频率,但是会出现谐波效率降低的问题.已有研究表明[54],谐波效率受波长影响,转换效率以波长标定为λ-5.5±0.5.为了克服这些困难,人们提出了多色组合激光场[55–72]驱动原子分子的方案,目的是获得兼具高发射效率和超宽平台的谐波辐射.实验研究方面,Shelton等[73]通过控制两路同步飞秒激光产生了短阿秒脉冲;Hassan 等[74]和Wirth 等[75]利用飞秒钛宝石激光多路参量放大的方案,开展了三路飞秒激光相干合成的实验研究;Hassan 等[76]进一步采用四路可见光波段的飞秒激光调控合成激光的电场波形得到了380 as 脉冲.理论研究方面,Takahashi等[59]通过混合多周期双色激光场产生连续高次谐波谱,获得强的单个阿秒脉冲.Feng 和Chu[77]通过控制三个激光场的激光参数,研究了He 原子在三色场中产生超短阿秒脉冲的最佳条件,表明通过合成两色或三色组合激光场,优化激光场波形,不仅可以增加谐波的截止频率,还可以在不增加激光总功率的情况下增强谐波强度到两个数量级,从而产生更强更短的单个阿秒脉冲.

实际上,理论模拟高次谐波及阿秒脉冲的产生还需考虑宏观介质演化效应的影响,求解宏观介质演化效应需结合含时薛定谔方程和麦克斯韦方程,计算工作量巨大.Jin 等[78]基于强场近似的重散射模型(quantitative rescattering theory,QRS)结合麦克斯韦方程求解两色和三色组合激光场,研究了宏观相位匹配效应对Ne 原子气体产生高次谐波的影响,该方法计算量相对较小,还可以解释大部分实验现象.为了更精确描述强激光场驱动原子分子产生高次谐波及单个阿秒脉冲,更深层次理解背后的物理机制,同时求解了含时薛定谔方程和麦克斯韦方程[79–82],首次将非约束优化算法(NEW unconstrained optimization algorithm,NEWUOA)[83]和含时薛定谔方程相结合,以扩展谐波平台最宽为目的构造目标函数,分别优化双色和三色组合激光场,并驱动He 原子产生超宽平台高次谐波及单个阿秒脉冲,同时基于求解麦克斯韦方程考虑了宏观演化效应的影响.研究结果表明,尽管宏观演化效应对单原子层次获得的阿秒脉冲产生有影响,但利用远场轴外量子轨道的空间选择性可以获得更短的单个阿秒脉冲,本文的研究为实验上产生单个超短阿秒脉冲提供了一种实际可行的新思路.

2 理论方法

在偶极近似和长度规范下,强激光场中原子的含时薛定谔方程为(原子单位):

其中,H0表示无外场时原子的哈密顿量,U(r,t) 表示原子与激光场的相互作用项,分别表示如下:

式中,Vl是He 原子模型势,是球谐函数,f(t) 是激光场包络,E0是激光场电场振幅,ω0是激光场频率,ϕ是相位.为了精确计算He 原子的谐波谱,构造了角动量依赖的模型势[84],表示如下:

式中,α是 He+偶极极化率;rc为有效核半径;N=1;S,A1,B1,A2,B2分别为常数,Wn为截断函数,表示如下:

表1 给出了He 原子角动量依赖模型势的具体参数.

表1 He 原子模型势参数(原子单位制)[65]Table 1.Model potential parameters for He (in a.u.)[65].

利用伪谱法(generalized pseudospectral method,GPS)[83]精确数值求解了含时薛定谔方程(1).伪谱法数值求解有两方面优点: 首先,伪谱法在空间采用非均匀离散,在核附近离散取点较密,而离核较远取点较稀疏,与传统的等间距网格离散方法相比,不但计算精度高,而且计算效率也成倍提高;其次,伪谱法把空间径向坐标从 (0,rmax] 转换到[-1,1],有效地解决了库仑奇点的问题.

采用分裂算符方法处理含时演化问题,含时演化波函数表示如下:

一旦得到随时间变化的波函数ψ(r,t),就可以计算加速度形式下诱导偶极矩的期望值:

通过对加速度形式诱导偶极矩做傅里叶变换,可以得到单原子层次的高次谐波谱:

宏观介质中的传播效应需要同时考虑激光场和谐波场的演化,由麦克斯韦方程描述.对于激光场的演化,三维的麦克斯韦方程表示如下:

式中,El(r′,z′,t)是激光脉冲的电场分量,Jabs(t) 是由于介质的电离作用而引起的吸收项,ηeff是气体的有效折射率,Jabs(t) 和ηeff分别表示为

其中,γ(t)是电离率,Ip是电离势,η0是考虑中性原子的折射和吸收效应,η2描述的是依赖瞬时激光强度I(t)的非线性光学克尔效应,ωp为等离子体频率,可以通过自由电子密度ne(r′,z′,t) 来计算,表示为

在圆柱坐标中,r'是径向距离,z'是轴向传播距离.

对于谐波场演化,三维的麦克斯韦方程表示如下:

其中,极化场D(r′,z′,t) 可以分为线性分量Dl(r′,z′,t)和非线性分量Dnl(r′,z′,t).线性分量包括高次谐波的色散和吸收效应,非线性分量和介质有关.在本文的计算中,气体压强较低,色散和吸收效应贡献很小[85],所以线性分量可以忽略.非线性分量Dnl(r′,z′,t)可以表示为

其中,n0是中性原子密度,dA(r′,z′,t) 是单原子加速度形式偶极矩.三维的麦克斯韦方程(9)和方程(13)在频域空间通过Crank-Nicolson 方法数值求解,获得宏观高次谐波谱,表示如下:

相应的阿秒脉冲强度由下式得出:

为了更好地分析高次谐波谱的发射特性,利用小波变换[86]对时间域谐波谱Eh(t′) 进行时频分析,

式中,W[ω(t′-t)] 是Morlet 小波变换的窗函数:

3 结果与讨论

3.1 双色组合激光场优化及其驱动原子产生高次谐波和单个阿秒脉冲

首先利用NEWUOA 算法,基于优化控制理论(optimal control theory,OCT)对双色组合激光场进行优化.优化过程中,固定第一色激光场波长和强度以及第二色激光场强度,优化第二色激光场波长和两束激光场之间的延迟时间,以扩展谐波平台最宽为目标函数,表示如下:

其中,

式中,ωL和ωH分别是优化目标函数对应积分频率的下限和上限.

计算中采用的双色组合激光场为

其中,E1(t)为第一色激光场,E2(t) 为第二色激光场,分别表示如下:

式中,A1,A2分别为双束激光场的电场振幅;ω1,ω2分别为两束脉冲的角频率;τ1,τ2为两束脉冲的持续时间;td为两束脉冲间的延迟时间.

图1(a)给出了优化前和优化后的双色组合激光场.初始激光场波长分别设定为800 nm 和400 nm,强度分别设定为 3.0×1014W/cm2和0.8×1014W/cm2,半高全宽均为8 fs 的高斯包络激光脉冲.数值求解含时薛定谔方程的过程中,结合NEWUOA 算法,以构造的目标函数优化迭代计算,优化得到第二色激光场波长为2754.9 nm,两束脉冲之间的延迟时间为 0.022T1(T1= 2π/ω1).可以看出优化后的激光场波形相比优化前变化很大.图 1(b)给出了相应的利用优化前和优化后的双色组合激光场驱动He 原子产生的高次谐波谱.优化后的谐波平台相比优化前扩展了130 多阶,谐波从70 阶到截止位置附近变得光滑连续,即优化后的激光场驱动He 原子产生了超连续平台,但是谐波强度有所下降,这是因为优化后的组合激光场有效波长增大.事实上,谐波转换效率降低这一不利因素可以通过宏观相位匹配来增强[87].

图1 (a) 优化前和优化后的双色组合激光场;(b)优化前和优化后的双色组合激光场驱动He 原子产生的高次谐波谱Fig.1.(a) Initial laser field and the optimized two-color laser field;(b) corresponding HHG power spectra.

为了进一步理解谐波的发射特性,利用小波变换对优化前和优化后的高次谐波谱进行了时频分析.图 2(a)为优化前高次谐波谱的小波时频分析结果,很明显谐波的发射是周期性变化的,在70 阶左右发射效率迅速降低,这与图 1(b)对应的谐波截止位置符合.在一个光学周期内,谐波的发射主要是由长轨道电子和短轨道电子贡献产生.根据半经典“三步模型”,长轨道电子为先电离后返回,短轨道电子为后电离先返回.图 2(b)为优化后高次谐波谱的小波时频分析结果,对比图 2(a)可以看出,谐波平台得到了极大的扩展,优化后的谐波在180 阶之后效率迅速降低,与图 1(b)对应的谐波截止位置一致,另外,谐波的发射时间主要集中在8.25 光学周期,短轨道电子贡献被抑制,长轨道电子贡献增强,这种一个量子轨道被抑制,另外一个量子轨道被增强的现象有利于超连续结构的产生.

图2 (a) 优化前高次谐波时频分析图;(b)优化后高次谐波时频分析图Fig.2.Wavelet time-frequency of the HHG spectra for the cases of (a) initial and (b) optimized two-color laser fields.

为了进一步阐明长轨道电子和短轨道电子对每一阶次谐波发射的贡献,图 3(a)给出了从70 阶到170 阶给定阶次谐波发射强度随时间的变化,在8.0 光学周期附近谐波的发射被抑制,而在8.2到8.3 光学周期发射比较集中,几乎是同时的,这与图 2(b)中短轨道电子贡献被抑制,长轨道电子贡献被增强是一致的.已有的研究结果表明[88],如果谐波发射相位是锁定的,就具有超连续平台结构,叠加这些谐波可以产生更强更短的阿秒脉冲.为了证实我们的判断,我们叠加了70 阶到170 阶的谐波谱,获得了能带宽度为155 eV 的25 as 单个阿秒脉冲,如图 3(b)所示.为了方便比较优化前与优化后合成的阿秒脉冲,图 3(b)中优化后的阿秒在强度上乘了7 倍,与优化前获得的最短单个阿秒脉冲相比,单个阿秒脉冲形状经过优化更加规则光滑,同时阿秒脉宽从130 as 缩短到了25 as.

图3 (a) 优化后给定阶次高次谐波强度随时间变化关系;(b)对比优化前后单个阿秒脉冲的产生Fig.3.(a) Dipole time profiles of harmonics from the 70th to the 170th harmonic order;(b) attosecond pulse generation for the cases of the initial and optimized two-color laser fields.

3.2 三色组合激光场优化及其驱动原子产生高次谐波和单个阿秒脉冲

为了得到更短的单个阿秒脉冲,接下来利用NEWUOA 算法对三色组合激光场进行优化.计算中,固定三色激光场波长和强度,优化第一色激光场与第二色激光场之间的延迟时间和第一色激光场与第三色激光场之间的延迟时间,以扩展谐波平台最宽为目标函数,构造目标函数如下:

其中,

式中,ωL和ωH分别是优化目标函数对应积分频率的下限和上限.

利用三色组合激光场驱动He 原子产生高次谐波及单个阿秒脉冲,计算中采用的三色组合激光场为

其中,三束激光场分别对应为E1(t),E2(t) 和E3(t),表达式为

式中,A1,A2,A3分别为三束激光场的电场振幅;ω1,ω2,ω3分别为三束脉冲的角频率;τ1,τ2,τ3分别为三束脉冲的持续时间;td2,td3分别为第一束脉冲与第二束脉冲间的延迟时间,第一束脉冲与第三束脉冲间的延迟时间.

图4(a)给出的是优化前和优化后的三色组合激光场.初始激光场波长分别为800,2400 nm 和1600 nm,对应强度分别为I1=3.0×1014W/cm2,I2=0.8×1014W/cm2和I3=0.2×1014W/cm2,半高全宽均为15 fs 的高斯包络激光脉冲.数值求解含时薛定谔方程的过程中,结合NEWUOA 算法,以构造的目标函数优化迭代计算,优化得到第一色激光场与第二色激光场之间的延迟时间为-0.2308T1(T1=2π/ω1),第一色激光场与第三色激光场之间的延迟时间为 -0.2764T1.对比优化前后的三色组合激光场,激光场峰值强度基本没有改变,但是激光场波形被改变,这种变化使得组合激光场有效波长增大,有利于扩展高次谐波平台.图 4(b)给出了相应的优化前和优化后三色组合激光场驱动He 原子产生的高次谐波谱,优化前的高次谐波谱有三个平台,第一个平台在150 阶左右截止,第二个平台在220 阶左右截止,第三个平台在260 阶左右截止,即最大截止位置.优化后的谐波谱平台减少为两个,130 阶之前为第一平台,130—290 阶为第二平台区,而且第二个平台区变得相对规则光滑,谐波效率也有所增强.另一个明显的变化是优化后谐波的第一平台截止位置相对优化前的情况向低阶谐波移动了20 阶左右,为获得更宽谐波平台创造了条件.最终优化后有效谐波平台比优化前扩展了80 阶左右.

图4 (a) 优化前和优化后的三色组合激光场;(b)优化前和优化后的三色组合激光场驱动He 原子产生的高次谐波谱Fig.4.(a) Initial laser field and the optimized three-color laser field;(b) corresponding HHG power spectra.

如图 5(a) 利用小波变换对优化前三色组合激光场驱动He 原子产生高次谐波谱进行了时频分析,优化前高次谐波谱的发射主要由4 部分组成,分别位于14.7,16.3,17.7 和19.3 光学周期,对应的最大阶次分别为150 阶、260 阶、210 阶和75 阶.其中,平台最大峰值和次最大峰值发射都比较强,对最小峰值75 阶次以下谐波发射都有强烈的贡献,导致截止峰值75 阶次在图 4(b)中并不明显.实际上,在图 5(b)中优化后高次谐波谱的发射也主要由四部分组成,但主要的贡献来自于16.3 和17.7 光学周期的发射,对应的截止峰值分别为130阶次和290 阶次,与图 4(b) 结果一致.特别需要指出的是,在16.15 光学周期处出现了一个非常集中的发射区域,即从130 阶到280 阶谐波的发射几乎是同时的.对比优化前后的时频分析结果,发现长轨道电子贡献增强,短轨道电子贡献被抑制,这有利于单个阿秒脉冲的产生.

图5 (a) 优化前高次谐波时频分析图;(b)优化后高次谐波时频分析图Fig.5.Wavelet time-frequency of the HHG spectra for the cases of (a) initial and (b) optimized three-color laser fields.

图6(a)给出了110—200 阶给定谐波阶次发射强度随时间变化的结果.谐波发射主要集中在16.1—16.2 光学周期 (o.c.),这些阶次谐波的发射几乎是同时的,有利于单个阿秒脉冲的生成.图6(b)中叠加了110—280 阶的谐波谱,获得了能带宽度为260 eV 的17 as 单个阿秒脉冲,相比优化前90 eV 的29 as 单个阿秒脉冲,通过优化获得了小于一个原子单位时间(24 as)的单个阿秒脉冲.

图6 (a) 优化后给定阶次高次谐波强度随时间变化关系;(b)优化前后单个阿秒脉冲的产生Fig.6.(a) Dipole time profiles of harmonics from the 110th to the 200th harmonic order;(b) attosecond pulse generation for the cases of the initial and optimized three-color laser fields.

3.3 介质演化效应对产生高次谐波和阿秒脉冲的影响

以上讨论了单原子层次双色和三色组合激光场优化及其驱动He 原子产生高次谐波及单个阿秒脉冲,表明优化的双色和三色组合激光场驱动He 原子可以得到更宽的谐波平台和更短的单个阿秒脉冲.事实上,实验中激光场振幅在空间有一定的分布,成千上万个原子(分子)在激光场中感受到的光强是不一样的,并且发出的高次谐波辐射在空间受到介质演化效应的影响.接下来讨论介质演化效应对产生高次谐波和阿秒脉冲的影响.为了最大程度上减小宏观介质演化效应对单原子层次优化激光场波形的影响,选取了较低的气体压强和束缚较松的激光脉冲束腰.以优化的双色组合激光场驱动He 原子产生高次谐波及单个阿秒脉冲为例,结合激光场中含时薛定谔方程和介质演化效应的麦克斯韦方程,讨论宏观介质演化效应对产生高次谐波和单个阿秒脉冲的影响.宏观介质演化效应需要同时考虑激光场的演化和谐波场的演化.计算中,激光聚焦强度、中心波长、以及脉宽等参数与单原子优化后双色组合激光场情况一致.空间依赖的双色组合激光场表达式如下:

其中,E1(r′,z′) 和E2(r′,z′) 表示空间依赖的第一束和第二束激光场电场分量,f1(r′,z′,t) 和f2(r′,z′,t-td)分别表示两束激光场的高斯激光包络,ω1和ω2是对应的双色激光场的频率,ϕ1和ϕ2是相位,td是两束激光脉冲间的延迟时间,两束激光的聚焦强度分别为I1=3.0×1014W/cm2和I2=0.8×1014W/cm2.两束激光脉冲聚焦处束腰均为50 µm,He 原子气体靶位于激光聚焦位置处,气体长度为0.5 mm,气体压强为30 Torr (1 Torr≈133.322 Pa).

图7 (a)是优化后双色组合激光场进入气体靶前的空间分布,图 7(b)是相应演化后的激光场空间分布.很明显激光场场强在空间不同位置强度是不同的.为了更清楚地展示激光场演化前后的区别,图 7(c)给出了演化前与演化后的轴上激光场对比,在峰值处,演化前后强度有所不同.为了参考,图 7(d)比较了单原子层次优化后的双色组合激光场和宏观演化后轴上激光场,因为计算过程中采用的气体压强较小,激光脉冲束腰较松,所以在峰值位置处强度的差别很小,但是对比仍可以发现介质演化效应对激光场波形存在一定影响.

图7 优化的宏观双色组合激光场空间传播效应 (a)演化前激光场;(b)演化后激光场;(c)轴上演化前和演化后激光场对比;(d)驱动单原子激光场与宏观演化后轴上激光场对比Fig.7.Macroscopic propagation effects of the optimized two-color laser fields: (a) Entrance;(b) exit;(c) comparison of entrance and exit on axis;(d) comparison of the fields for single-atom case and mac-field on axis.

图8(a)是双色优化后组合激光场考虑介质宏观演化效应的高次谐波空间分布,计算中设置在激光聚焦位置后30 cm 的远场处收集谐波.从图中可以看出,谐波的空间分布是不均匀的,在轴上的分布比较强,远离轴的分布比较弱.轴上100—190 阶出现了超连续空间分布结构,与单原子情况相符合.图 8(b)给出的是远场谐波对空间求和的结果,为了比较列出了单原子结果,单原子谐波谱在强度上抬高了12 个量级.比较而言,宏观谐波谱更加光滑规则,超连续平台向较低谐波区域移动.图 8(c)给出了叠加谐波谱从60 阶到200 阶得到的阿秒脉冲空间分布.由于在轴上和轴外谐波的传播距离不同,导致在轴外的区域阿秒脉冲发射有延迟,即阿秒脉冲空间分布在轴外有明显向右的移动.由于轴外各点位置激光强度不同,谐波截止频率也不一样,用同样范围的谐波谱叠加获得阿秒脉冲波形会有很大的差异,但由于量子轨道的空间选择性,总可以在轴外合适的位置获得单个超短阿秒脉冲.

图8 (a) 远场高次谐波空间分布;(b)单原子与远场高次谐波对比;(c)远场阿秒脉冲空间分布Fig.8.(a) Spatial distribution of far-field HHG;(b) comparison of single atom and far-field HHG;(c) spatial distribution of far-field attosecond pulse.

图9(a)和图 9(b)给出的是远场轴上和轴外0.37 mm 处(图 8(c)中绿点线所示)的高次谐波时频分析.对于轴上的情况,发射主要集中在8.2到8.3 光学周期,与单原子情况一致,也就意味着在本文采用的激光场和靶结构参数下,介质演化效应对轴上谐波影响不大.对于轴外情况,不同空间位置长轨道电子和短轨道电子的贡献是不一样的,离轴越远,长轨道电子贡献越强,短轨道电子贡献越弱,谐波发射的同步性越好,有利于产生更短的单个阿秒脉冲.如图 9(c)所示,叠加轴上第66—191 阶谐波得到了28 as 的单个阿秒脉冲,与单原子结果相差不大,而轴外0.37 mm 处叠加第65—200 阶谐波得到最短为18 as 的单个阿秒脉冲.为了方便比较,轴外单个阿秒脉冲强度乘了300 倍.结果表明,利用量子轨道的空间选择性,可以在轴外获得更短的单个超短阿秒脉冲.

图9 (a) 轴上高次谐波时频分析;(b)轴外0.37 mm 处高次谐波时频分析;(c)轴上与轴外0.37 mm 处获得的单个最短阿秒脉冲Fig.9.(a) Wavelet time-frequency of the HHG spectra on axis;(b) wavelet time-frequency of the HHG spectra at 0.37 mm off axis;(c) on-axis and off-axis attosecond pulse generation.

最后,为了检验原子气体靶的位置对高次谐波超连续谱的影响,将He 原子气体靶分别放置在激光聚焦位置前和激光聚焦位置后2 mm 处,其他参数与图 7 的情况完全一致,计算了相应的轴上单个阿秒脉冲的产生,如图 10 所示.激光聚焦位置前(–2 mm)、激光聚焦位置处(0 mm)和激光聚焦位置后(2 mm)分别在轴上得到了最短25 as,28 as,28 as 的单个阿秒脉冲,阿秒结果相差不大,但激光聚焦位置之后的单个阿秒脉冲更加规则光滑,因为气体靶放置在激光聚焦位置之后高次谐波相位匹配更好,这一结论与实验结果一致[89].

图10 轴上不同气体靶位置产生的单个阿秒脉冲对比Fig.10.Comparison of single attosecond pulse for the different target position.

4 总结与展望

利用优化组合激光场驱动原子分子在强激光场中发射的高次谐波是合成超短阿秒脉冲的主要方案之一,高次谐波及阿秒脉冲的产生是强场超快过程的前沿研究领域,涌现出了许多重要研究成果,限于篇幅,不能一一提及.本文仅从理论方面研究了利用优化组合双色和三色激光场驱动原子发射的高次谐波及合成单个超短阿秒脉冲的最优方案,并展示了研究结果.将含时薛定谔方程和NEWUOA 优化算法相结合,通过优化激光场波形,极大地增加了谐波的截止频率,对于双色和三色激光场优化及其驱动He 原子分别获得了最短25 as 和17 as 的单个阿秒脉冲.同时以双色激光场优化为例,基于同时求解含时薛定谔方程和麦克斯韦方程,进一步考虑了宏观介质演化效应对高次谐波和阿秒脉冲产生的影响,发现在远场轴外通过量子通道的空间选择性可以获得更短更规则的单个18 as 超短阿秒脉冲,相关研究为实验提供了更精准的理论参考.随着激光技术的发展,实验室中已经可以产生任意波形的激光脉冲,结合优化控制理论,优化组合激光场驱动原子分子产生高次谐波及单个超短阿秒脉冲具有非常明显的优点,使获得更短单个阿秒脉冲光源成为可能,并将成为揭示新物理机制的重要手段.

感谢西北师范大学王志斌博士和陈阳阳硕士的讨论.