强激光场驱动Ar 原子电离中的隧穿延时*

赵猛 全威† 肖智磊 许松坡 王志强 王明辉 成思进 吴文卓 王艳兰 赖炫扬 柳晓军‡

1) (中国科学院精密测量科学与技术创新研究院,波谱与原子分子物理国家重点实验室,武汉 430071)

2) (中国科学院大学物理科学学院,北京 100049)

3) (湖北大学物理系,武汉 430062)

阿秒钟实验方案是当前研究原子分子体系的价电子在强激光场中隧穿延时问题的有效手段.基于阿秒钟方案,本文实验研究了Ar 原子在800 nm 椭圆偏振激光场中的光电子动量分布随激光光强的演化规律.理论上采用包含库仑场效应、非绝热效应、Stark 效应、多电子屏蔽和极化效应的半经典模型对Ar 原子的强场电离动力学进行了模拟.通过对比实验测量和数值模拟结果发现,在本文所研究的光强范围内,Ar 原子的价电子在800 nm 椭圆偏振激光场中隧穿延时上限为10 as.进一步分析表明,阿秒钟方案中,考虑多电子屏蔽效应对得到的隧穿延时影响最小,而考虑非绝热效应的影响最大.

1 引言

动能小于势垒的粒子以量子隧穿的方式从势垒中通过是否需要有限时间(隧穿延时问题)? 尽管该问题已经被争论了90 年[1],但它仍然是当前的研究热点.随着超快强激光和谱仪技术的发展,研究者可以在实验室中,以超快的时间分辨能力,研究电子在激光场和库仑场联合势垒中的超快隧穿过程.最近十几年的进展催生了一系列新的实验技术和方案,如阿秒条纹相机[2]、高次谐波结合泵浦探测[3]、阿秒瞬态吸收[4]和阿秒钟[5−8]等.

近年来,阿秒钟实验方案受到了广泛关注[9,10],采用该方案仅使用飞秒激光脉冲即可实现阿秒量级的时间分辨.阿秒钟的基本工作原理为,在近圆偏振激光电场的偏振平面上,特定激光周期内不同时刻隧穿的光电子会出射到不同的方向,因此,光电子的末态动量角度投影可反映其隧穿时刻.对于800 nm 近圆偏振激光场而言,光电子动量分布中出射角度偏转1°对应时间变化为7.4 as.

基于阿秒钟实验方案,研究者从实验和理论两个方面对隧穿延时问题开展了大量研究,然而研究结论却不一致.Eckle 等[5]使用阿秒钟方案研究了He 原子的隧穿延时,发现在实验误差允许范围内,He 原子价电子在强激光场中的隧穿过程是瞬时发生的.Torlina 等[11]数值模拟了H 原子在椭圆偏振激光场中的电离过程,计算结果再次表明H 原子价电子的隧穿过程不需要时间.Sainadh 等[9]制备出H 原子样品,并基于阿秒钟方案开展实验研究,结合数值求解三维含时薛定谔方程的结果,证明了H 原子的隧穿过程是瞬时的.Quan 等[12]针对H2分子开展实验研究,通过与多种理论方法对比,给出电子隧穿延时上限为10 as.然而,与以上研究结论相反,也有部分研究的阿秒钟实验结果表明电子隧穿需要有限时间.Landsman 等[13]分别采用速度影像谱仪和冷靶反冲离子动量谱仪研究了He 原子在较大光强范围内的隧穿延时,通过与相应的理论结果[14]对比,发现实验结果与Lamor时间定义和费曼路径积分方法得到的隧穿延时接近.Camus 等[15]实验研究了Ar 和Kr 在相同条件下的光电子末态动量分布,并考察了二者的差别.由于理论计算中必须包含隧穿延时和初始纵向速度分布才能较好地重复实验结果,该工作推断隧穿过程必然伴随着确定的延时.

导致研究结论不一致的一个重要原因是,阿秒钟方案的结果可能受到诸多物理效应的影响,如库仑场效应、非绝热效应、Stark 效应、多电子屏蔽和极化效应.目前并不清楚这些效应将如何影响阿秒钟方案的结果.针对这一研究现状,基于阿秒钟方案,实验测量了一系列光强条件下Ar 原子在椭圆偏振激光场中电离产生的光电子动量分布.理论上使用包含库仑场效应、非绝热效应、Stark 效应、多电子屏蔽和极化效应的半经典模型,数值计算了对应的光电子动量分布.通过对比实验测量和数值模拟结果,提取出了Ar 原子价电子的隧穿延时.研究结果表明,在本文研究的光强范围内,800 nm 近圆偏振激光场中Ar 原子的价电子隧穿延时上限为10 as.数值计算和分析表明,阿秒钟方案中,考虑多电子屏蔽效应对得到的隧穿延时影响最小,而考虑非绝热效应的影响最显著.

2 实验方案

本文的实验测量采用冷靶反冲离子动量成像谱仪(COLTRIMS)完成.商用钛宝石飞秒激光器(Femtopower Compact PRO CE-Phase HP/HR)产生中心波长800 nm,重复频率5 kHz,脉宽30 fs,单脉冲能量0.8 mJ 的线偏振激光脉冲.实验中采用消色差半波片和格兰泰勒棱镜的组合来调节激光光强,并进一步利用宽带四分之一波片得到椭偏率为0.7,激光偏振长轴平行于z轴的椭圆偏振激光场.进入COLTRIMS 谱仪作用腔的激光脉冲被焦距为75 mm 的凹面镜反射,并聚焦于超声Ar原子束上,与其相互作用后产生光电子和离子.在谱仪内部均匀直流电场(E=25 V/cm) 和磁场(B=9 Gs,1 Gs=10–4T)的作用下,光电子和离子飞向各自带有延迟线阳极的微通道板探测器,并被其探测.光电子和离子的三维动量分布可由探测器测量的飞行时间和碰撞位置得到.本文中,激光光强由近圆偏振激光场中原子电离产生的电子动量分布标定[16].

3 理论方法

为揭示不同物理效应对阿秒钟实验方案研究结果的影响,使用半经典理论模型[17−19]对原子在激光场中的电离过程进行数值模拟.半经典理论方法的优势在于计算量小,且可通过研究特定电子轨道揭示直观物理图像.在本工作中,采用半经典理论方法,可方便地比较多种物理效应对实验结果的影响.该模型目前已被广泛应用于强场电离研究[20−24].除非专门说明,本文默认采用原子单位制(a.u.).

在半经典理论计算中,电子通过隧穿跃迁到连续态[25,26]后的运动使用牛顿方程来描述[27−32]:

其中E(t)=(Ex(t),0,Ez(t))是椭圆偏振激光电场:

E0是激光电场振幅,ε是椭偏率,ω是激光角频率,a(t)是包络函数.

式中,T为激光光学振荡周期.对于只考虑库仑场效应的情况,电子感受到的有效势为

其中,Zeff=为有效核电荷数,r表示隧穿电子和离子实之间的距离.

如果认为电子在激光场和库仑场联合势垒中的隧穿过程是绝热的,可采用准静态近似[25]来处理电子的隧穿过程.在抛物坐标系中,隧穿电子的初始位置可通过解薛定谔方程来获得[26]:

而每条轨道的权重,即电离率为

如果认为电子在激光场和库仑场联合势垒中的隧穿过程是非绝热的,则应考虑非绝热效应对隧穿电子的隧穿出口、电子初始速度和电离率的影响[33,34].在采用鞍点近似的S-Matrix 理论中[35,36],鞍点解ts=t0+iti对应电子轨道的贡献最重要,这里t0是上文提到的隧穿时刻,ti为电子在势垒下演化的虚时间.电子末态动量p与初始动量p0的关系可以表示为p=p0-A(t).实验室坐标系下的初始动量p0可由旋转坐标系下平行于瞬时激光电场的分量p//和垂直于瞬时激光电场的分量p⊥来表示,px0=p//sinθ+p⊥cosθ,pz0=p//cosθ-p⊥sinθ,θ为偏振平面内瞬时激光电场与z轴的夹角,而p//可由下式得到:

当t=t0时,r0(t) 的实部就是隧穿出口位置:

在计算中,隧穿时刻t0和初始电子横向速度p⊥通过随机数给出,ti可通过求解鞍点方程[35,36]得到,而电子初始纵向速度p//,电离概率Ω以及隧穿出口x0,y0,z0分别由公式(7),(8),(9)给出,从而将非绝热效应对隧穿过程的影响体现在电子轨道的初始条件.

为考虑多电子屏蔽效应的影响，这里采用Hartree-Fock-Slater 模型势作为一价Ar 离子的束缚势.该束缚势可写成如下形式[37,38]:

其中Φ(r)为多电子屏蔽效应的影响.在本文中,(11)式的参数取值为Z=18,b=5.4,c=1,d=3.682 时可准确模拟一价Ar 离子的束缚势.如果进一步考虑离子实多电子极化效应的影响,则束缚势可表示为[31,39]

在本文中,为考虑Stark 效应的影响,原子的电离势可表示为[40]

其中Ip0是不考虑外场情况下的电离势,αN是Ar 原子极化率,αI是一价Ar 离子极化率.

为更准确地对比实验测量和理论计算结果,在计算中考虑了聚焦平均效应[41].作用区光强分布为

其中I0为激光场的峰值光强,r为待研究原子到激光束中心轴的距离,z为沿激光传播方向的坐标,R(z)=R0[1+(z/zR)2]1/2,这里瑞利长度zR=,λ是激光波长,R0是激光焦点处的束腰半径.谱仪中激光束瑞利长度0.9 mm,而作用腔中超声束半径1.0 mm.考虑到超声束半径与瑞利长度接近,可认为在作用区内光强分布与三维高斯分布接近.考虑三维聚焦平均时,单位体积内电离率由Ω(I,p)表示,则.焦点附近的光束体积可以表示为V=,这里r(I,z)可由上文中的I(r,z)得到.最终电离产量分布为

4 结果分析

实验测量得到的Ar 原子在800 nm 椭圆偏振光偏振平面中的典型光电子二维动量分布如图1(a)所示,其对应的角分布如图1(b)所示.激光峰值光强为1.5×1014W/cm2,椭偏率为0.7.图1(c)和图1(d)是相同条件下的半经典理论计算结果.图1(a)中的插图标识出激光偏振长轴和短轴的方向.简单推导可知[42,43],如果忽略库仑场的影响,且电子在激光电场旋转到z方向(激光偏振长轴方向)时发生隧穿,则电子的最终速度将沿着x方向(垂直于长轴方向).然而,由图1(a),(b)可见,光电子角分布的最大值并未准确分布在短轴(x)方向.进一步观察可知,动量分布和角分布也并非理想高斯分布.图1(c),(d)中的半经典计算结果也有类似特征.根据文献[17]的分析可知,这一现象的来源是库仑场效应.长程库仑势一方面会导致整个光电子二维动量分布发生转动,另外一方面也会使光电子二维动量分布和角分布发生畸变.

本文采用文献[17]提出的方法得到激光峰值电场对应的电子发射方向,并将该方向与激光场短轴方向的夹角定义为α.图1(a),(c)标示出α的正方向.图1(b)展示了从图1(a)中的实验数据得到的光电子角分布.对于本文中的多周期超快激光场而言,激光电场的光学振荡可看成振幅不变的正弦或余弦振荡.对于此类振荡而言,通过选择激光电场为零的两个相邻时间点之间的半个光学周期,可使得激光电场峰值附近的电场形式前后对称.进一步由激光电场与原子电离率具有一一对应的关系可知,该激光电场峰值对应的电子最终出射角度可将该时间范围内的电子产量二等分.文献[17]的数值计算表明这一结论可推广到常用多周期激光场整个脉冲的情况.因此,可通过寻找图1(b)的角分布中将电子产量二等分的角度来确定α.本文采用该方法从实验测量和数值计算数据中得到的α如图1(b)和图1(d)中向上箭头所示.图1(b)和图1(d)中β是相应光电子产量最大值所对应的角度,也是前人阿秒钟实验中所使用的偏转角.相较β,α更准确地给出了在激光峰值电场发生隧穿的光电子轨道对应的角度.

图1 (a) 实验测量的Ar 原子光电子动量分布,激光波长800 nm,光强1.5×1014 W/cm2,椭偏率0.7,插图标识出激光场偏振面的长轴和短轴方向;(b) 黑色方框表示从(a)中提取的Ar 原子光电子角度分布;(c) 半经典理论计算的光电子动量分布,激光参数与(a)相同;(d) 黑色圆圈表示从(c)中提取的Ar 原子光电子角度分布.(b)(d)中的红色曲线代表采用最小二乘法对实验测量和数值模拟数据的拟合结果Fig.1.(a) The measured photoelectron momentum distributions of Ar subject to the laser electric field with a wavelength of 800 nm,an intensity of 1.5×1014 W/cm2,and an ellipticity of 0.7,the inset shows the directions of the major and minor axes of the polarization ellipse of the laser field;(b) the black boxes indicate the photoelectron angular distribution extracted from the data in (a);(c) the photoelectron momentum distributions calculated by the semiclassical model with the laser parameters identical to those of(a);(d) the black circles represent the photoelectron angular distribution extracted from the data in (c).The red curves in (c) and(d) show the results of the least square fittings of experimental data and numerical calculations,respectively.

为系统研究椭圆偏振激光场中Ar 原子的价电子隧穿延时问题,实验测量了一系列光强条件下Ar 原子的光电子动量分布.使用上文描述的方法,提取了在激光峰值电场发生隧穿的电子轨道对应的偏转角度α,如图2(a)中黑色实心方框所示.测量结果表明随着激光光强升高,α逐渐减小.该结果与前人针对惰性气体原子和小分子体系的阿秒钟方案实验测量结果一致.采用第3 节中理论方法部分介绍的半经典模型,分别计算了仅包含库仑场效应(CP 海军蓝空心三角形)、包含库仑场效应和多电子屏蔽效应(Screen.+CP 蓝色空心四边形)、包含库仑场效应和非绝热效应(Nonad.+CP 橄榄绿空心五边形)、包含库仑场效应和多电子极化效应(Polar.+CP 洋红色空心五角星)、包含库仑场效应和Stark 效应(Stark+CP 深黄色空心六边形)、以及同时考虑以上提到的所有物理效应(All红色空心圆形)的情况下α随光强的演化.从图2(a)可见,在本文所研究的光强范围内,与仅考虑库仑场效应的情况相比,同时考虑多电子屏蔽效应和库仑场效应，以及同时考虑Stark 效应和库仑场效应两种情况下,α变化较小.而同时考虑多电子极化效应和库仑场效应,以及同时考虑非绝热效应和库仑场效应两种情况下,α显著变小,且这些效应的影响随光强增大而增强.在本文所考察的几种物理效应中,相对于仅有库仑场效应的情况而言，对α变化影响最小的是多电子屏蔽效应,其次是Stark 效应,再次是多电子极化效应,而影响最大的是非绝热效应.阿秒钟实验方案中,需要通过比较实验测量与半经典理论计算得到的α来提取隧穿延时.因此,对于半经典理论计算结果影响越大的物理效应,也将对最终提取的电子隧穿延时结果影响越大.由以上讨论可知,阿秒钟方案中,考虑多电子屏蔽效应对得到的隧穿延时影响最小,而考虑非绝热效应的影响最大.

为从实验测量和理论计算结果中提取隧穿延时,需使用角度α与隧穿时刻的对应关系.对椭圆偏振激光场而言,激光电场矢量旋转角 Δθ与隧穿时刻 Δt之间的关系[44]为

使用(15)式,并假定 Δθ=α,可从图2(a)中的α得到电子隧穿发生的时刻,结果如图2(b)所示.图2(b)中的黑色实心方块为从图2(a)中的实验测量α得到的对应电子隧穿时刻.如该图所示,实验得到的电子隧穿时刻随光强增大而呈缓慢减小的趋势.图2(b)中黑色空心圆(without FA)为忽略聚焦平均效应的半经典理论计算结果,黑色空心四边形(FA)为考虑三维聚焦平均效应后的半经典理论计算结果.考虑聚焦平均效应之后,半经典理论计算结果与实验数据更加接近.图2(b)中,将二者的偏差用洋红色实心圆表示(图2(b)右侧纵轴).由该结果可知,在本文所研究的光强范围内,强激光场中Ar 原子的价电子隧穿延时接近或小于实验测量误差,并可知其上限为10 as.

图2 (a)实验测量与半经典理论计算得到的α 随激光光强的依赖曲线;黑色实心方块(Exp.)为实验测量结果.海军蓝空心三角形(CP)为仅包含库仑场效应的计算结果、蓝色空心四边形(Screen.+CP)为包含库仑场效应和多电子屏蔽效应的计算结果、橄榄绿空心五边形(Nonad.+CP)为包含库仑场效应和非绝热效应的计算结果、洋红色空心五角星(Polar.+CP)为包含库仑场效应和多电子极化效应的计算结果、深黄色空心六边形(Stark+CP)为包含库仑场效应和Stark 效应的计算结果,红色空心圆形(All)为以上所有效应同时考虑的情况.(b) 由(a)中α 得到的电子隧穿时刻随光强的依赖曲线.黑色实心方块(Exp.)为实验结果,黑色空心圆(Without FA)为单一光强计算结果,黑色空心四边形(FA)为考虑聚焦平均效应的结果.采用阿秒钟方案得到的隧穿延时随激光光强的变化由洋红色实心圆(Tunneling delay time)表示(右侧纵轴)Fig.2.(a) The measured and calculated intensity dependence of α.The filled black boxes (Exp.) show the experimental results.The open navy triangles (CP),open blue diamonds (Screen.+CP),open olive pentagons (Nonad.+CP),open magenta pentacles(Polar.+CP),open dark yellow hexagons (Stark+CP),and open red circles (All) indicate the semiclassical calculations where only the influence of Coulomb potential,the influence of both Coulomb potential and the multi-electron screening effect,the influence of both Coulomb potential and the non-adiabatic effect,the influence of both Coulomb potential and multi-electron polarization effect,the influence of both Coulomb potential and the Stark effect,and all the physical effects have been considred,respectively.(b) The laser intensity dependence of tunneling instant obtained from the measured and calculated α in (a).The filled black boxes (Exp.)represent the experimental data.The open black circles (Without FA) indicate the calculation without focusing average.The open black diamonds (FA) represent the calculation where the focusing average has been considered.The filled magenta circles (Tunneling delay time) indicate the intensity dependence of the tunneling delay time which is obtained based on the attoclock scheme.

为进一步深入理解阿秒钟实验中相关物理效应的影响,分别计算了在激光峰值电场发生隧穿的典型电子轨道在考虑多种物理效应影响的结果.图3 对比了仅考虑库仑场效应和同时考虑库仑场效应及非绝热效应两种情况下的典型电子轨道.

图3(a)展示了在激光峰值电场发生隧穿的典型电子轨道在隧穿出口附近的空间演化.计算结果表明,考虑相同物理效应的情况下,激光光强越高,隧穿出口距离库仑场中心越近(红色粗点线和比红色细实线更接近z=0,蓝色粗点划线比蓝色细划线更接近z=0).与仅有库仑场效应的计算结果(红色细实线和红色粗点线)相比，非绝热效应(蓝色细划线和蓝色粗点划线)会导致隧穿出口离库仑场中心更近,这一影响在低光强条件下更加明显.如图3(c)所示,pz(沿着激光峰值电场方向的动量)随激光场振荡的幅度更小(蓝色细划线位于红色细实线上方,红色粗点线与蓝色粗点划线几乎重合),表明非绝热效应所带来的隧穿出口差别使得低光强条件下电子在隧穿出口感受到的库仑场影响更大.如图3(b)所示,非绝热效应会导致电子轨道在x方向上(垂直于激光峰值电场方向)有明显的初始速度,且光强越高,这个初始速度越大(在0 时刻,蓝色粗点划线位于蓝色细划线下方).电子初始横向速度会导致最终光电子能量升高,相应地减弱库仑场的影响.另一方面,考虑到非绝热效应会导致隧穿出口离库仑场中心更近,从而增强库仑场影响,因此理解非绝热效应的总体影响还需要进一步考察光电子最终的动量和位置.

图3 半经典理论计算的在激光峰值电场发生隧穿的电子轨道 (a),(d)隧穿刚发生时和激光脉冲快结束时电子轨道的空间演化;(b),(e)隧穿刚发生时和激光脉冲快结束时,沿x 方向电子动量(px)随时间的演化;(c),(f)隧穿刚发生时和激光脉冲快结束时,沿z 方向电子动量(pz)随时间的演化,其中红色细实线和红色粗点线表示仅考虑库仑场效应的轨道,蓝色细划线和蓝色粗点划线表示同时考虑库仑场和非绝热效应的轨道,红色细实线和蓝色细划线表示激光光强为 1.5×1014 W/cm2 的计算结果,红色粗点线和蓝色粗点划线表示激光光强为 2.3×1014 W/cm2 的计算结果Fig.3.The typical trajectories calculated by the semiclassical model for the photoelectrons tunneling from the peak of laser field:(a),(d) The spatial evolution of electron trajectories around the tunneling instant and the end of the laser pulse;(b),(e) the temporal evolution of photoelectron momenta along the x direction (px) around the tunneling instant and the end of the laser pulse;(c),(f) the temporal evolution of photoelectron momenta along z direction (pz) around the tunneling instant and the end of the laser pulse.The thin red solid and thick red dotted lines represent the trajectories calculated by the semiclassical model where only the influence of Coulomb potential is considered.The thin blue dashed and thick blue dot-dashed lines represent the trajectories calculated by the semiclassical model where the influence of both Coulomb potential and the nonadiabatic effect are considered.The thin red solid and thin blue dashed lines indicate the trajectories which are calculated at 1.5×1014 W/cm2.The thick red dotted and thick blue dot-dashed lines indicate the trajectories which are calculated at 2.3×1014 W/cm2.

图3(d)—(f)为激光脉冲快结束时典型电子轨道随时间演化情况.如图3(d)所示,在激光场结束时,两种光强条件下,同时考虑非绝热效应和库仑场效应的电子轨道到达的位置离库仑场中心更远,表明非绝热效应总体上使得库仑场的影响减弱,且在x方向上的相对差别比z方向上相对差别更大.由于z方向的相对差别主要来自于非绝热效应对隧穿出口的影响(如图3(a)所示),而x方向的相对差别主要来自于非绝热效应对电子初始横向速度的影响(如图3(b)所示),图3(d)的结果说明电子初始横向速度的变化所产生的影响更大.比较图3(e)和图3(f)的结果可知,尽管随着光强增大,px和pz都在增大,但px的相对值增加更快.结合图1 中α的定义可知,这一影响将导致α随光强减小,从而解释了图2(a)中非绝热效应导致α随光强减小的现象.

为比较不同物理效应的影响,在图4 中给出了,在考虑多个物理效应的情况下,激光峰值电场对应典型电子轨道隧穿时刻的电子初始位置、电子初始动量以及脉冲结束时电子最终位置和最终动量.由图4(a),(c)可知,除了非绝热效应,其它物理效应对电子轨道初始条件的影响主要体现在位于z方向隧穿出口的差别.相应地,不同物理效应对应的典型轨道在脉冲结束时的最终位置(图4(b))和最终动量(图4(d))差别也主要体现在z方向上.除了多电子极化效应的影响以外,差别的大小顺序正好与隧穿出口一致.由公式(12)可知,多电子极化效应对应的束缚势不仅与电子所处的位置有关,而且与激光振荡电场的大小有关.而其他效应,比如: 库仑场效应和多电子屏蔽效应,所对应的束缚势仅仅与电子所处的位置有关.所以,随着激光电场不断振荡,多电子极化效应对应的束缚势变化更加剧烈.并导致对应的电子典型轨道受束缚势的影响,与仅考虑库仑场影响的情况相比,在z方向和x方向均变大,因此多电子极化效应的影响比除了非绝热效应以外的其他效应影响更大.不仅如此,由于在z方向激光峰值电场更大,所以多电子极化效应会导致典型电子轨道在z方向变化更大.为了清楚地表现这一影响,注意到同时考虑库仑场效应和Stark 效应,以及同时考虑库仑场效应和多电子极化效应,两种情况下隧穿出口非常接近(图4(a)中深黄色六边形与洋红色五角星的位置非常接近),可通过对比这两种情况下最终电子位置和动量,定性揭示多电子极化效应在轨道演化过程中的影响.由图4(b)和图4(d)可见,与同时考虑库仑场效应和Stark 效应的轨道(深黄色六边形)相比,同时考虑库仑场效应和多电子极化效应的轨道(洋红色五角星)在z方向电子最终位置和动量的偏移量均大于在x方向的偏移量.这一结果与上文关于多电子极化效应对电子轨道演化过程影响的分析一致.图4(d)给出了在考虑多种物理效应的情况下,激光峰值对应电子轨道在脉冲结束时的动量大小.该图中,采用一条直线将不同物理效应对应的电子最终动量与原点相连,标示出电子最终出射方向.与仅考虑库仑场效应的情况(海军蓝三角形)相比,偏角越大,说明相应的物理效应影响越大.由图4(d)可知,非绝热效应(橄榄绿五边形)的影响最大,其次是多电子极化效应(洋红色五角星),而影响最小的是多电子屏蔽效应(蓝色四边形).

图4 (a),(c)半经典理论计算得到的激光峰值电场对应电子轨道在隧穿时刻(initial)电子的位置和动量;(b),(d)半经典理论计算得到的激光峰值电场对应电子轨道在脉冲结束时(final)电子的位置和动量;(d)中连接原点和对应符号的实线标示出电子的最终出射方向,海军蓝三角形(CP)为仅包含库仑场效应的计算结果、蓝色四边形(Screen.+CP)为包含库仑场效应和多电子屏蔽效应的计算结果、橄榄绿五边形(Nonad.+CP)为包含库仑场效应和非绝热效应的计算结果、洋红色五角星(Polar.+CP)为包含库仑场效应和多电子极化效应的计算结果、深黄色六边形(Stark+CP)为包含库仑场效应和Stark 效应的计算结果,实心(Low)和空心(High)分别表示激光光强为1.5×1014W/cm2 和2.3×1014 W/cm2 的相应计算结果Fig.4.(a),(c) The calculated positions and the momenta of the photoelectrons tunneling from the peak of laser field at the tunneling instant (initial);(b),(d) the calculated positions and the momenta of the photoelectrons tunneling from the peak of laser field at the end of the laser pulse (final);in (d),the solid lines connecting the origin and the symbols indicate the corresponding emission angles of the photoelectrons at the end of the laser pulse.The navy triangles (CP),blue diamonds (Screen.+CP),olive pentagons(Nonad.+CP),magenta pentacles (Polar.+CP),and dark yellow hexagons (Stark+CP) indicate the semiclassical calculations where only the influence of Coulomb potential,the influence of both Coulomb potential and the multi-electron screening effect,the influence of both Coulomb potential and the non-adiabatic effect,the influence of both Coulomb potential and multi-electron polarization effect,and the influence of both Coulomb potential and the Stark effect have been considred,respectively.The filled (Low) and open (High) Symbols indicate the results calculated at 1.5×1014 W/cm2 and 2.3×1014 W/cm2,respectively.

5 结论

实验测量了Ar 原子在椭圆偏振激光场中的光电子动量分布随激光光强的演化规律.理论上使用包含库仑场效应、非绝热效应、Stark 效应、多电子屏蔽和极化效应的半经典理论方法对Ar 原子的强场电离动力学进行了模拟,得到了Ar 原子电离产生的光电子动量角分布.基于阿秒钟方案,通过对比实验测量和理论计算结果,发现在本文所研究的光强范围内,Ar 原子的价电子在椭圆偏振强激光场中的隧穿延时小于10 as.理论上研究了考虑库仑场效应、非绝热效应、Stark 效应、多电子屏蔽和极化效应对阿秒钟方案结果的影响,并结合对典型轨道的分析阐明了其物理机制.