严以律，周 希，任山令，刘 辉，田友伟*

(1.南京邮电大学 贝尔英才学院，南京 210023；2.南京邮电大学 通信与信息工程学院，南京 210023；3.南京邮电大学 理学院，南京 210023)

引 言

自20世纪末啁啾脉冲放大(chirped pulse amplification,CPA) 技术问世以来，已有众多学者在强激光与物质相互作用领域[1-3]，尤其是超快激光物理方向[4-6]开展了相关研究。相对论强度激光束中的光子经与高能电子相互作用产生汤姆逊散射，而由于多普勒频移的存在，该过程是非线性的。利用非线性汤姆逊散射效应制成的激光器是X射线[7]与γ射线[8-9]的产生途径之一，也是单色激光[10-11]和超短阿秒(1as=10-18s)激光[12-13]的重要辐射源，这在生物医药[14-16]、原子物理与核物理[15,17]等学科中观测超快过程有着诸多应用。

应用场景的广泛促使研究者对相对论情况下非线性汤姆逊散射进行理论分析。LEE等人利用计算机模拟了起始位于坐标原点且静止的单个电子在不同偏振态的强激光作用下汤姆逊散射的空间特征，发现了频谱的调制结构[18]。ZHENG等人在LEE的研究基础之上探究了电子初速度不为0时电子的辐射能量及功率角分布[19]。LI等人则以入射激光强度为变量，提出在线偏振激光下空间辐射具有双叶性[20]。WANG等人以圆偏振紧聚焦激光与初始静止于原点的电子相互作用为研究对象，在3维坐标系下总结了角辐射与电子轨迹间的联系[21]。

本文中探讨了初始状态静止的单电子与线偏振紧聚焦激光相互作用的运动轨迹，并在前人的研究基础上研究了电子初始位置对高能电子空间辐射的影响。模拟结果表明，电子初始位置对全空间最大能量辐射方向及相应能量值均有较大的影响，随着电子沿z轴正向移动，最大能量辐射方向的极角逐步减小并趋于一常数，而最大辐射能量具有极大值，这说明通过合理设置电子的初始位置可以获得强度尽可能高的辐射。

1 电子与激光作用模型

笛卡尔坐标系下激光脉冲的相位φ具有如下关系式[22]:

φ=η+φR-φG+φ0

(1)

于是高斯激光脉冲的矢势可用下式表示[22]：

(2)

式中,al=exp(-η2/L2-ρ2/w2)，ρ2=x2+y2,δ为偏振参量，对于线偏振光，δ=0；x和y分别表示x,y方向单位矢量。应注意的是，上述模型中的时空坐标已分别被ω0-1和k0-1归一化，其中ω0和k0分别是激光的频率和波数。

直角坐标系中光场的矢势分量可以表示为[22]:

(3)

(4)

式中,θ=π-arctanφG。

线偏振高斯激光脉冲与单电子相互作用的示意图如图1所示。图中φ与θ分别代表电子与激光相互作用发出辐射的方位角与极角。这里，假设激光脉冲沿+z轴传播，初始时刻电子静止于z轴。

Fig.1 Schematic diagram of the interaction between the electron and the incident laser pulse

电子在电磁场中的运动可以用拉格朗日方程和电子的能量方程描述[22]：

(5)

(6)

将(3)式和(4)式代入(5)式和(6)式，得到如下方程组[22]：

(7)

式中，ux,uy,uz分别为电子在相应坐标方向上的速度分量，ax，ay，az分别表示矢势a在x,y,z这3个坐标轴上的分量。

通过求解以上偏微分方程组，可以确定电子在与激光相互作用过程中运动和能量的变化过程。

由电动力学相关知识，以相对论速度运动的电子会向空间发出电磁辐射。单位立体角内的辐射总能量可以表示为[21]：

(8)

t=t′+R0-n·r

(9)

式中,R0是观察点与电子和激光相互作用点之间的距离，r表示电子的位矢。

2 数值模拟结果

借助MATLAB软件对高能电子与线偏振高斯激光脉冲的作用过程进行计算模拟，得到电子在z轴具有不同初始位置时的运动轨迹及相应的空间辐射特性。所使用的激光脉冲为沿+z轴传播且极化方向沿x轴的线偏振光，考虑到当前实验室所能达到的条件，选取激光的归一化激光振幅a0=6(对应激光峰值光强I=4.97×1019W/cm2)、脉宽L=3λ0(10fs)、束腰半径w0=3λ0(3μm)，以及电子的归一化初始能量γ0=1(初速度为0)。电子初始位于z轴且位置区间为[-10λ0,10λ0]。

Fig.2 Motion trajectory of the electron and its spatial distribution of energy radiation when the electron has the initial position (0,0,-10λ0)

调整电子初始位置至(0,0,-7λ0)得到图3。此时电子轨迹的振荡部分和直线运动部分在z轴上的投影区域都有所展宽，且振荡幅度最大值较初始位置在(0,0,-10λ0)处的有所增大，其余性质与图2所示无异。空间能量辐射的分布特征几乎保持一致，能量仍是集中在关于z轴对称的两支上，此时的最大能量辐射方向(θmax,φmax)=(23.5°,180°)，对应的辐射能量为2.76×107，较前有所增加。

Fig.3 Motion trajectory of the electron and its spatial distribution of energy radiation when the electron has the initial position (0,0,-7λ0)

进一步将电子初始位置置于(0,0,-4λ0)得到图4。此时电子轨迹的振荡部分和直线运动部分在z轴上的投影区域进一步展宽，振荡幅度最大值也进一步略微增大。空间能量辐射分布仍旧几乎保持不变，此时最大能量辐射方向(θmax,φmax)=(22.5°,180°)，在该方向上的辐射能量为2.58×107，较前又有所减小。

Fig.4 Motion trajectory of the electron and its spatial distribution of energy radiation when the electron has the initial position (0,0,-4λ0)

由上可见，在电子初始位置由-10λ0逐渐向+z轴移动的过程中，最大能量辐射方向的方位角φmax始终为180°，而极角θmax则有逐渐减小的趋势。此外，最大辐射能量呈现出先增大后减小的规律。为进一步探究电子初始位置对空间能量辐射的影响，以λ0为间隔绘制了在[-10λ0,10λ0]区间内最大能量辐射方向及该方向上辐射能量与电子初始位置的关系，如图5所示。由于在该范围内φmax始终保持180°不变，故未在图中画出其变化趋势。

Fig.5 Relationship between the maximum radiated energy as well as its corresponding angle and the initial position of the electron

由图5可知，θmax随电子初始位置沿+z轴移动逐渐减小，这表明具有最大空间辐射能量的一支逐渐靠近z轴，并在z=2λ0处开始稳定在20.5°，此后保持该值不变。而最大辐射能量随z轴坐标的增加呈现先增后减的变化规律，且在z=-7λ0处达到最大值2.76×107。此外，还注意到即使θmax不同，但最大辐射能量在区间[-10λ0,-4λ0]上关于z=-7λ0具有一定对称性。

综上所述，初始状态静止的具有不同初始位置的电子与高能激光相互作用的轨迹都呈现沿+z轴振荡行进，最后沿直线运动的特征，且当电子初始位置位于(0,0,-7λ0)时，在最大能量辐射方向(θmax,φmax)=(23.5°,180°)处具有最大辐射能量2.76×107。

3 结 论

依据拉格朗日方程构建了高能电子与紧聚焦高斯激光脉冲的相互作用模型，利用MATLAB在单电子情形下通过对运动轨迹及辐射特性的模拟，探究了电子初始位置对其空间能量辐射的影响。结果表明，高能电子经与强激光的相互作用，其在xOz平面内先朝+z轴沿之字形运动，最后做直线运动。此外，模拟结果还显示出电子初始位置对空间最大辐射能量及其辐射方向均有较大影响，随着电子初始位置向+z轴移动，空间最大辐射能量具有先增大后减小的特征，存在一个极大值，且最大能量辐射方向的极角将逐渐减小并最终稳定在20.5°，而方位角始终保持180°不变。空间最大辐射能量在电子初始位置(0,0,-7λ0)、方向(θ,φ)=(23.5°,180°)处取得。