相对论皮秒激光在低密度等离子体中直接加速的电子束的横向分布特征研究*

张晓辉 董克攻 华剑飞 朱斌 谭放 吴玉迟 鲁巍 谷渝秋3)†

1) (清华大学工程物理系, 北京 100084)

2) (中国工程物理研究院激光聚变研究中心, 等离子体物理重点实验室, 绵阳 621000)

3) (深圳技术大学工程物理学院, 深圳 518118)

相对论皮秒激光与低密度等离子体作用可以通过“激光直接加速”机制获得超有质动力定标率的高能电子, 且电荷量可以达到百nC级, 在伽马射线产生、正电子产生等方面具有重要应用.然而激光直接加速电子束相比激光尾场加速电子束具有更大的发散角, 同时实验观测的横向束分布也不均匀, 但是其中的物理机制研究较少.本文通过二维粒子模拟证明, 相对论皮秒激光在低密度等离子体中驱动的激光直接加速中, 高能电子束会在激光偏振方向分叉, 而且电子能量越高这种现象越明显.文章通过细致的理论分析解释了这种高能电子横向分布产生“分叉”结构的内在原因.在激光直接加速的过程中, 电子在纵向获得加速的时候, 它在激光偏振方向(横向) betatron振荡的动能也会随之增加, 当电子的能量足够高时, 二者呈线性关系, 因此高能电子的横向速度的振幅近似相等, 这种相等的振幅最终导致了高能电子束在激光偏振方向的分叉.

1 引 言

在过去的40多年里, 激光等离子体加速技术获得了快速的发展[1,2], 使得廉价的、桌面型的加速器成为一种可能.在2004年的三个里程碑式的工作[3−5]之后, 基于飞秒 (femtosecond, fs)激光的激光尾场加速机制 (laser wakefield acceleration,LWFA)逐渐成为激光等离子体加速中的主流.LWFA可以获得更高能的电子, 同时具有良好的束流品质(如能散、发散角等), 然而电子束的电荷量受限于空泡尺寸很难提升, 同时对激光参数与等离子体参数的匹配有很高要求[6].利用皮秒(picosecond, ps)激光在低密度等离子体中加速电子[7,8], 则对等离子体参数与激光参数的匹配没有太多要求, 可以工作在较宽等离子体参数范围, 特别是在较高的等离子体密度条件下.虽然电子能量很难被加速到数百MeV以上[8], 但是可以在数十MeV的能量范围获得远高于LWFA机制的电荷量[9,10](可以达到数十nC甚至数百nC以上), 因而在韧致辐射[11]、正电子[12]、光核中子产生[13,14]等方面具有独特的优势和潜在应用.

相对论ps激光在等离子体里传播时, 电子会被激光的有质动力朝横向排开, 从而形成一条长长的等离子体通道[15].由于电荷分离场和纵向电流,通道内部存在径向的准静态电场Er和角向的准静态磁场Bθ, 这些准静态场对电子施加横向的聚焦力[8,16].电子在准静态场的作用下做横向betatron振荡.由于通道内还同时存在着向前传播的激光场, 这种横向振荡就会通过v×B的机制引起纵向的加速或者减速.在这个过程中电子直接从激光的横向电场获得能量, 因此这种加速机制被称作激光直接加速机制 (direct laser acceleration, DLA)[17].如果满足共振条件电子就可以持续的被激光场加速, 其中代表betatron振荡的频率,代表电子感受到的多普勒频移之后的激光频率,ωp代表等离子体频率,ω0代表激光频率,γ代表电子的相对论因子,c代表真空中的光速,vx代表电子的纵向速度(激光传播方向),vph代表激光的相速度.

在LWFA中, 等离子体波内同样存在横向的准静态聚焦场.如果激光脉宽长到足以覆盖被注入的电子, 电子就可能受到LWFA-DLA的共同作用[18−20].实验表明, 无论是 ps激光驱动的 DLA 加速还是fs激光驱动的LWFA-DLA混合加速, 激光偏振方向的电子散角都会更大一些[11,21].此外,在LWFA-DLA混合加速的实验中, 人们发现在垂直于激光偏振方向测量电子能谱的时候, 电子能谱上会出现叉子状的结构[18,19].尽管粒子模拟(particle in cell, PIC)的结果已经证明这种分叉结构是由DLA引起的[18], 但是DLA究竟如何导致电子束的分叉, 目前仍然没有一个很好的物理解释和定量分析.本文通过二维PIC模拟表明, 在ps激光驱动的DLA中电子的能谱也会在激光偏振方向呈现叉子状分布.我们通过进一步的理论分析获得了产生这种横向分布的原因:在DLA的过程中, 电子在激光偏振方向的能量会随电子能量的增加而成比例地增加, 而垂直于偏振方向的能量则随着电子能量的增加而减小, 从而导致电子在激光偏振方向有更大的散角; 偏振方向能量的增加还会导致高能电子在偏振方向的速度具有相同的振幅, 这种相同的振幅进而导致了电子能谱的分叉结构.本文第2部分为PIC数值模拟结果, 第3部分为解析模型分析及讨论, 最后为结论.

2 皮秒激光与低密度等离子体相互作用的PIC模拟

为理解激光直接加速的物理过程和高能电子束的特性, 我们利用二维PIC程序OSIRIS[22]模拟了ps激光在低密度等离子体中驱动的激光直接加速的过程.模拟选用全离化的He等离子体, 其密度分布由 500 µm 的上升沿, 2 mm 的平台和500 µm 的下降沿组成, 最大电子密度为 2 × 1019cm–3.激光偏振方向在模拟平面内, 聚焦在等离子体上升沿的尽头.激光波长为 1053 nm, 脉宽为 0.8 ps, 焦斑为 7.5 µm, 激光无量纲化的矢量势a0= 3.模拟窗以光速移动, 大小为 700 µm × 300 µm, 被分为22000 × 1884 个网格, 对应的单位激光波长内的网格数分别别为33个和6.28个.我们在每个网格里面放2个宏电子和1个宏He2+粒子.模拟时间步尽量接近库朗条件(Courant condition)以避免数值噪声.

图1(a)和图1(b)分别为 PIC模拟中t=时刻的激光强度包络和电子的密度分布.我们可以看到, 激光仅在其头部驱动出等离子体波, 同时在这个区域激光本身也被等离子体波调制, 遭受到严重的不稳定性.随后, 大部分电子被激光排开形成一条长长的等离子体通道, 而激光则被很好地约束在这条通道里.由于通道里面大部分低能电子都被排开了, 因此通道内没有明显的等离子体波.图1(c)展示了通道内的聚焦场Es=Eys−cBzs.在二维模拟中, 径向的准静态电场Er退化为y方向的准静态电场Eys, 而角向的准静态磁场Bθ则退化为z方向的准静态磁场Bzs.在轴线附近Es与横坐标y近似为线性关系.因为通道内同时存在聚焦场与激光场, 因此电子可以通过DLA获得加速.同时, 由于通道内不存在等离子体波, 电子经历的是纯净的DLA, 因此这个区域非常适合研究DLA产生电子束的特征.

图1 PIC 模拟中 t=5965时刻的 (a)激光强度包络,(b)电子的电荷密度分布与(c)通道内的聚焦场Es=Eys−cBzs, 模拟中, 等离子体密度为 2 × 1019 cm–3, 激光脉宽为0.8 ps, a0= 3Fig.1.(a) The laser envelope; (b) electron density; (c)channel focusing force Es=Eys−cBzs at t=5965into the simulation, in which the plasma density is 2 ×1019 cm–3 and the laser have a duration 0.8 ps with a0= 3.

图2 t=5965 时刻电子在相空间的分布 (a)电子在能量-vy相空间的分布, 白色虚线是电子横向速度振幅的理论值, 右侧的黑色实线代表着能量大于60 MeV的电子的vy的分布, 为了更好地展示, 其计数值做了归一化处理; (b)能量在60—70 MeV之间的电子在y-py相空间的分布Fig.2.Electron phase space at t=5965 :(a) Energy -vy phase space, the white dashed lines denote the amplitude of vy from analytical solution,the black solid line denotes the vy distribution of electrons above 60 MeV, the counts are normalized to achieve a better illustration; (b) the y-py phase space of electrons within energy range from 60 MeV to 70 MeV.

3 理论分析与讨论

在DLA中, 电子主要受到两种场的作用:一个是向前传播的激光场, 另外一个是通道内的准静态电场Er和准静态磁场Bθ.为了推导电子的横向运动, 我们将激光场简化为无限长均匀平面波

其中η(>c) 代表激光电场与磁场之间的比值, 激光沿+x方向传播,y为激光偏振方向.为了简化模型, 我们将准静态场简化为线性的准静态电场Es=κr, 其中Es=Er−cBθ代表Er和Bθ的复合作用,r=y+z代表电子在横向的位置, k 是聚焦场强度的系数.那么电子在横向所受的来自准静态场的力可以表示为Fy= −eκy,Fz= −eκz, 进而电子动力学方程可以写为

根据图1可知, 准静态场远小于激光场, 因此等式(6)可以简化为

我们定义电子在y,z方向的总能量分别为

其中等式右边第一项为电子在该方向的动能, 后一项为电子在该方向的势能.对ϵy和ϵz分别求导, 可以得到

将(4)和(5)式代入 (10)和(11)式, 上面两式可以简化为

而对高能电子而言,

其中 ∆η=η−c, 可以当作常数来处理.那么等式(12)可以进一步化简为

当电子能量比较高的时候, 假设z方向的能量被完全阻尼掉,vz= 0, 那么电子的运动就变成了二维的, 此时等式(15)可以简化为

对上式积分就可得到

为了验证等式(17)是否成立, 我们在PIC模拟中追踪了100个获得DLA加速的电子的轨迹.图3展示了这些电子在这段时间内的能量g与横向能量ϵy的变化关系.这段时间通道内的准静态场相对稳定, 且电子的能量比较高, 因此比较适合与我们的解析分析进行对比.从图3可以看出, 尽管局部有些偏差, 但整体上电子横向能量ϵy在模拟中随着电子能量的近乎线性地增加.图中红色实线是根据等式(17)拟合得到的结果, ∆η=5.48×10−3c,b=1.11m0c2.这说明我们对电子横向能量的分析是可信的.

当电子运动到轴线上时, 也就是y= 0的时,其中下标A代表电子在轴线上的时刻, 那么我们就可以得到

进而电子在轴线上时的横向速度vyA可以表示为

y= 0的时刻往往也是其横向速度最大的时候.等式(19)表明, 在DLA中电子的横向速度的振幅仅由电子的能量以及初始条件所决定.当电子能量比较高的时候, 等式右边的后两项趋近于0, 就可以进一步化简为

vyA趋近于一个常数, 也就是说, 当电子的能量足够高以后, 它的横向速度的振幅是恒定的, 这个振幅主要由激光电场和磁场的比值决定.等式(20)表明, 在DLA加速中, 高能电子在激光偏振方向具有相同的最大速度.

接下来讨论相同的横向最大速度如何导致了电子束的分叉.我们先考察单个电子的横向速度的分布.在DLA中, 电子的横向振荡近似为简谐振荡.在此假设下, 电子的横向速度和加速度可以表示为:

其中vm代表vy的振幅, 根据前面的分析我们知道vm≈vyA.电子的速度和加速度满足关系

如果我们在足够长的一段时间内(远大于电子振荡的周期)统计这个电子vy的概率分布就会得到

根据等式 (24), 当vy=±vm时,趋近于无穷大,也就是说,vy在它的极值处出现的概率最高.等式(24)描述的是单个电子在一段时间内横向速度的概率分布, 对于一团电子束而言, 如果它们横向最大速度相同, 而且在各自振动的相位上均匀分布, 那么在某个时刻对这些电子的横向速度做统计, 它们的分布就与单个电子在一段时间内横向速度的概率分布相同.前面我们已经知道, 在DLA中高能电子横向速度具有相同的振幅, 因此它们横向速度的分布也可以用等式(24)描述.根据等式(24), 电子束会聚集在 ±vm附近, 这点也在图2(a)中得到了验证.图2(a)中的白色虚线是根据等式(19)得到的电子横向速度的振幅, 其中Dh和b的值是从图3拟合得到的.可以看到, 电子出现概率最高的地方正是横向速度的最大值.尽管在图3的拟合中方差比较大, 但是在图2(a)中电子最大值却与理论值符合得非常好.

图3 在 t=4965到 t=6965这段时 间 内 100 个被追踪的电子的能量g以及横向能量 ϵy 的变化, 图中红色实线是根据等式(17)拟合得到的结果Fig.3.The transverse energy ϵy of 100 electrons as a function of g from t=4965 to t=6965.The red dashed line is the fitted result according to Eq.(17).

在实际情况中,vy的振幅不可能完全相同.假设vy的振幅围绕vm有一定的展宽s, 我们在图4中通过数值计算展示了不同展宽下vy的分布.图4(a)展示的是三种不同的振幅展宽, 图1(b)展示了对应展宽下vy的分布.我们可以看出, 当s比较小的时候,呈现两边高中间低的近中空结构, 峰值的位置出现在 ±vm处.随着s的增大, ±vm处的峰逐渐展宽, 峰值逐渐降低,趋向于均匀分布.图2(b)中的黑线具有与图4(b)中的蓝线相似的结构, 并且波峰与波谷的比值更大, 这说明在PIC模拟中,vy的振幅的展宽很小.以上PIC模拟、理论分析以及数值计算的结果足以说明, 我们对电子分叉的解释是正确的.虽然我们的数值模拟和解析分析是针对皮秒激光直接加速开展的, 但是经过简单的修正, 也可以适用于飞秒激光直接加速的情况.

图4 (a)电子速度振幅 vm 的三种不同展宽; (b) 三种展宽下对应的 vy 的分布Fig.4.(a) Three distributions of vm ; (b) the corresponding distributions of vy.

4 结 论

我们通过二维PIC模拟证明了相对论皮秒激光在低密度等离子体中驱动的激光直接加速中, 电子能谱会在激光偏振方向分叉.我们通过理论分析解释了电子束分叉的成因, 并对分叉的程度进行了定量的分析.在激光直接加速的过程中, 电子在纵向获得加速的时候, 它在激光偏振方向(横向)betatron振荡的动能也会成随之增加.当电子的能量足够高时, 二者呈线性关系, 因此高能电子的横向速度的振幅近似相等.这种近似相等的振幅导致了电子束在激光偏振方向的分叉.我们的推导还解释了在DLA中为什么电子的散角总在激光偏振方向更大.这些解释虽然是针对皮秒激光直接加速研究的结果, 但是经过简单修正后同样适用于飞秒激光直接加速的情形.电子横向速度的定量推导具有重要意义, 可以为很多场合进行指导.比如在LWFADLA混合加速中, 如果已知电子的横向速度, 就可以通过调控等离子体密度来实现锁相加速, 或者在DLA中通过调控Dh来控制betatron辐射的散角等.