激光幅度和光束腰半径对交叉碰撞中紧密聚焦的圆偏振场中的阿秒脉冲列的影响

《Frontiers in Physics》：Influence of laser amplitude and beam waist radius on attosecond pulse trains in a cross-collision tightly focused circularly polarized field

　　摘要
非线性汤姆逊散射（RNTS）产生的辐射已成为生成高度定向和宽带阿秒发射的重要机制。然而，激光场结构对电子动力学及其产生的辐射特性的影响尚未得到系统的研究。在这项工作中，我们采用了一个高度聚焦的圆偏振高斯激光模型，该模型包含了高阶场分量，以研究在激光振幅a0（归一化激光振幅a0a0）和光束腰半径b0（归一化光束腰半径b0=ω0/λ0）共同变化下的高能电子的轨迹演化和辐射特性。研究结果表明，a0主要控制纵向加速和光谱展宽，而b0则调节相互作用长度，并显著影响辐射的方向性。这些参数共同决定了发射的阿秒脉冲的空间峰位置、光束准直度和结构特征。本研究揭示了激光场结构在RNTS中的调节作用，并为优化高定向性和宽带阿秒光源提供了理论指导。

1 引言
随着超强激光技术的快速发展，极端电磁场下的激光-物质相互作用已成为现代强场物理学的中心课题。此类相互作用广泛应用于材料加工和医疗治疗[1, 2]，以及原子和强场物理学的基础研究[3, 4]。特别是啁啾脉冲放大（CPA）技术的出现，使得激光峰值强度大幅提高，同时减少了脉冲持续时间并增强了可实现的脉冲能量[5–8]。这种超强和超短激光脉冲的出现，使得激光-电子相互作用进入了相对论性和非线性效应起决定性作用的领域。因此，涉及强激光场中自由电子的汤姆逊散射进入了相对论性非线性领域，通常称为RNTS。由于其根本重要性和作为紧凑辐射机制的潜力，RNTS一直吸引着理论和实验研究的持续关注。
早在1998年，Umstadters就首次观察到了由高强度激光场与自由电子相互作用产生的非线性汤姆逊散射，包括谐波生成和辐射角分布的测量[9]。随后，Frühling等人结合实验和理论分析，研究了激光偏振（如 circular 和 elliptical 偏振）对散射辐射的偏振特性和产生的谐波结构的影响[10]。在经典理论框架内，Boca和Oprea 对90°交叉碰撞配置进行了深入研究，比较了不同碰撞几何形状对电子轨迹和相应辐射光谱的影响，同时还探讨了载波-包络相位（CEP）在短脉冲场中的作用[11, 12]。此外，Taira等人提出，在圆偏振激光照射下产生的高阶γ射线谐波具有光学涡旋特性，展示了圆偏振调节散射辐射结构独特的能力[13]。
以往的实验和理论研究，特别是基于严格经典处理的方法，显著推进了对RNTS的理解[14–16]。然而，大多数这些研究采用了平面波模型或弱聚焦激光场，这些模型无法描述高度聚焦的圆偏振光束所固有的高阶场分量和纵向效应。此外，现有研究通常分别考察激光或电子参数的影响，而多个激光参数的耦合效应尚未得到充分探索。此外，尽管在几项理论分析中考虑了交叉碰撞配置，但高能电子与高度聚焦的圆偏振激光脉冲相互作用产生的辐射的结构特性尚未得到系统研究。大多数以往的工作主要关注频率域中的谐波光谱和角分布，而负责相干辐射调制和阿秒脉冲形成的时域机制仍不清楚。

2 理论与公式
首先需要指出，在以下推导中，所有空间和时间坐标均以归一化单位表示，其中空间坐标按 k?10=λ0/2π 缩放，时间坐标按 ω?10=λ0/2πc 缩放。这里，λ0=1μm 表示激光波长，c=2.998×108m/s 表示光速。
在以下各部分中，所有对 b0 的引用均指归一化光束腰半径 b0=ω0/λ0。因此，在这项工作中，我们系统研究了激光振幅 a0 和光束腰半径 b0 的共同变化如何影响高能电子与高度聚焦的圆偏振激光场之间交叉碰撞中的电子轨迹、辐射方向性和产生的阿秒脉冲的特性。我们旨在阐明 a0 的增加如何增强纵向加速，使电子更接近近轴区域，而 b0 的增加如何延长相互作用范围，从而进一步积累纵向能量。这两个参数的共同作用产生了可调的横向振荡幅度，从而塑造了发射辐射的空间分布。在空间功率-角分布中，a0 主要决定了辐射强度和主要峰方向，而在高强度条件下，b0 引入了峰位置和光束准直度的协同调制。适当选择 a0 和 b0 的组合可以显著提高辐射方向性。时域和频率域分析进一步揭示了明显的相干调制结构：a0 控制光谱展宽和高频分量的产生，而 b0 则精细化了光谱特征和能量分布。

3 结果与讨论
在这项研究中，我们研究了高能单个电子与高度聚焦的圆偏振高斯激光脉冲交叉碰撞产生的非线性汤姆逊散射的辐射特性。所有空间坐标均归一化为激光波长 λ。电子沿着正的XX轴传播，最初位于x=-15λ的位置，确保在与激光脉冲相互作用之前有足够的自由传播距离。电子的横向和纵向初始坐标分别设为y=0和z=0，使得电子与激光场中心的相互作用发生在焦点(0, 0, 0)。电子的初始能量设为γ0=10，意味着它携带的能量大约是静止能量的九倍（相当于约4.6 MeV的动能）。圆偏振激光脉冲沿着正的ZZ方向传播，脉冲长度为L=3λ，对应于时间持续时间为τ=L/c。激光的振幅a0和束腰b0分别在6-10和1-10的范围内变化，相关的瑞利长度为zR=πb02/λ，确保电子在主要相互作用期间保持在焦点附近。所有参数都在经典电动力学的有效范围内；因此，本工作中采用的数值模拟方法是适当且可靠的[28, 29]。激光是圆偏振的（δ=1），并且工作在高强度范围内，能够驱动明显的横向振动并使电子获得显著的能量增益。

3.1 电子轨迹
在本工作中建立的三维笛卡尔坐标系中，XX轴和YY轴被定义为横向方向，而ZZ轴被定义为纵向方向。对于所有的a0和b0组合，高能电子在激光场中既表现出横向振动也表现出纵向漂移，如图2所示。纵向漂移是由ponderomotive力引起的，而电子的横向振动是由激光的振荡电场驱动的。
图2：在不同a0-b0参数设置下，与激光脉冲交叉相互作用后产生的电子轨迹，标记为(a–i)。
关于横向振动，通过对轨迹的检查和坐标数据的仔细比较，可以发现沿XX轴和YY轴的运动之间存在明显的不对称性：XX方向的振动幅度明显大于YY方向。在理想的圆偏振激光场中，电子应该在横向yx平面上进行幅度相等的对称旋转。然而，数值模拟显示沿XX轴的振动比沿YY轴的振动更强。一个原因是电子并没有准确地从光轴开始，因此XX方向的样本场强高于YY方向，导致XX方向的振动更为明显。另一个原因是与光束的空间结构和聚焦有关：在激光振幅a0和束腰b0的共同变化下，沿XX方向的振动幅度可以估计为x∝(a0λ)/(2π)，表明其对于a0有线性敏感度，而对b0的依赖性较弱。相比之下，沿YY方向的分量主要由与包络梯度相关的ponderomotive力驱动；在轴附近，以瑞利长度zR=πb02/λ作为特征相互作用尺度，位移可以估计为y∝(π2a02b03)/(8λ2)。因此，沿YY方向的振动主要受束腰b0的控制，较小的b0会导致更大的振动幅度。
对于纵向漂移，即电子沿ZZ轴的运动，无论是绘制的轨迹还是详细的坐标数据都表明，其位移随时间单调增加。这种行为是由两个因素引起的：首先，电子从激光的电场和磁场的共同作用中获得净向前动量；其次，激光包络沿ZZ方向的强度梯度通过ponderomotive力产生额外的向前推动。这两种机制的联合效应导致了整体上单调的纵向漂移。当a0和b0同时变化时，沿ZZ轴的漂移主要归因于在横向振动期间通过vx×By力获得的净向前动量，其周期平均效应可以由ponderomotive力描述。因此，纵向漂移的幅度与激光振幅的平方成正比，可以近似表示为pz∝a02（归一化到mc）。与横向振动不同，纵向漂移对b0的依赖主要通过场的空间非均匀性体现出来：对于较大的b0，激光沿传播方向几乎是均匀的，漂移呈现出相对平滑、单调的增加；对于较小的b0，增强的场梯度导致向前推力的相位依赖性变化，从而导致漂移轨迹的压缩和变形，以及更复杂的整体形状。根据上述分析，a0和b0的共同变化对纵向漂移的影响可以总结如下：漂移对a0的依赖性很强，而对b0的依赖性则体现在漂移模式的复杂性增加上。
因此，同时变化激光振幅a0和束腰b0会在激光场中引入电子的纵向加速度和横向振动的新行为。a0的增加显著增强了纵向加速度，使电子更容易在靠近轴的区域获得高能量；同时，较大的b0延长了激光与电子之间的相互作用区域，允许纵向加速度在更长的距离上积累，并进一步增强了靠近轴的运动。此外，横向振动表现出明显的方向依赖性：沿XX轴，振动幅度仅受b0的轻微影响，因此耦合变化对位移的影响有限；然而，沿YY轴，a0和b0的共同增加同时增强了包络梯度和场驱动力，导致更明显的位移变化。a0和b0的耦合控制不仅放大了纵向加速度，还引入了不对称的横向振动。这种不对称性直接体现在辐射的空间功率-角度分布上，这将在下一节中详细分析。
为了澄清a0和b0的不同作用，纵向加速度主要来源于周期平均的v⊥×B力。由于横向动量与a0成正比，纵向动量遵循pz∝a02，表明a0直接决定了相对论性纵向加速度的强度。相比之下，束腰参数b0控制了瑞利长度zR=πb02，这决定了电子与聚焦激光场之间的有效相互作用距离。因此，虽然a0控制了加速度强度，但b0主要影响相互作用长度。

3.2 电子辐射的空间功率-角度分布
在本节中，我们关注电子与高强度激光脉冲发生交叉碰撞时产生的辐射的空间功率-角度分布，以及与之相关的指标，如Zmax、θpeak和?peak。空间功率-角度分布揭示了脉冲辐射的能量如何在不同角度方向上分布，从而提供了关于其定向性和聚焦性的关键见解。这里，ZZ表示归一化的角度辐射功率密度，Zmax表示其最大值，表征了辐射发射的相对强度和角度集中度。此外，a0和b0的共同变化对辐射的方向性和强度都有不同的影响。这些发现对于优化辐射收集和设计高亮度辐射源具有实际意义。
图3显示了九种参数组合下的空间功率-角度分布的二维投影图，其中a0=6,8,10，b0=3,6,9。从这些图中可以明显看出，随着a0和b0的增加，辐射的方向性和强度都发生了显著变化。整体图案表现出一个特征性的“马蹄形”结构，辐射能量主要集中在θ方向，并且相对于?轴大致对称。随着a0和b0的同时增加，辐射峰值逐渐向较小的θ角度移动，亮区扩大，辐射结构变得更加明显。虽然不同b0组合下的辐射分布在?方向上的对称性基本保持一致，但随着a0的增加，出现了更加集中的结构，导致辐射方向更加集中。总体而言，辐射保持了良好的方位角对称性，尽管在某些样本中可能会发生轻微的偏差。
图3：电子与紧密聚焦的激光脉冲交叉碰撞后产生的归一化NITS辐射功率的角分布，对于几种a0-b0参数设置，标记为(a–i)。
在两个参数同时变化的情况下，a0（激光振幅）的增加对应于场强的增加。这导致电子在与激光的交叉碰撞过程中获得更高的相对论能量。此外，更强的激光场增强了电子的非线性振动和纵向漂移，从而增加了高阶谐波辐射的贡献。束腰b0的增加延长了瑞利长度zR=πb02/λ，这减弱了横向强度梯度。结果，横向ponderomotive力Fp∝??a2减小，电子的偏转被抑制了。a20与b20的比率也是影响横向发散程度的另一个重要因素。对于选定的三种情况，a0=6, b0=3、a0=8, b0=6和a0=10, b0=9，计算出的比率分别约为4、1.78和1.23，显示出明显的减弱趋势。这表明随着a0和b0的共同增加，横向发散有效地减小，确保了辐射方向更加一致。
我们从所有样本分布中提取峰值功率Zmax，并构建了一个二维热图来说明峰值功率如何随着a0和b0的共同变化而变化。这里，a0反映了激光强度，而b0表示激光束腰半径。颜色尺度和标记大小编码了不同参数组合的Zmax的大小。所有Zmax的值都通过全局最大值进行了归一化，因此颜色条的范围是从0到1。
随着a0的增加，激光强度上升，电子可以被驱动发出更强的辐射峰值，从而显著增强了Zmax。b0的增加效果相对于a0来说较弱，只有当a0已经较大时才变得显著。较小的束腰b0提供了更强的聚焦，但限制了相互作用体积，从而限制了可实现的峰值功率。相反，即使b0较大，较小的a0也无法产生显著的峰值增强；然而，在高a0下，较大的b0允许系统提供更强的峰值功率。这表明两个参数之间存在协同放大效应：较大的a0确保了强烈的单电子辐射，而较大的b0提供了电子与激光之间更长且更均匀的相互作用区域，共同导致了峰值功率的增强。这些趋势在图4所示的二维峰值功率图中进行了总结。
图4：在a0-b0参数空间中的峰值辐射功率Zmax的二维图，颜色和标记大小表示相对峰值功率水平。
对于每个a0值，我们在b0=1,5,10时选择最大的极角θ和方位角Φ，并将它们作为a0的函数绘制出来，水平轴表示a0，垂直轴表示θpeak或?peak。使用线图和散点图来可视化趋势。在散点图中，为了防止具有相同值的点重叠，在代码中向?peak添加了小的随机抖动，从而使得观察和分析更加清晰。提取的θpeak和?peak值分别在图5a和图5b中标出。
图5(a)：在束腰值b0=1,5,10的情况下，峰值极角θpeak作为激光振幅a0的函数，说明了在不同聚焦条件下的主导发射方向的变化。(b) 对于相同的参数集，相应的方位角峰值偏差Φpeak相对于Φ=0的散点分布，标记中应用了轻微的抖动以减少重叠并提高可区分性。
对于θpeak曲线，当b0=1时，情况接近72°，并且随着a0的增加仅表现出轻微的变化。相比之下，对于b0=5和10，θpeak随着a0的增加而显著变化：更强的激光场增强了电子的相对论加速度，使得辐射更加向前定向，因此减少了θpeak。同时，b0的增加减弱了聚焦效果并减少了横向场梯度，降低了电子的横向ponderomotive偏转，导致辐射角度变窄。当a0a0和b0b0同时增加时，横向力场会减弱，而纵向加速度会增强，这种综合效应使得θpeakθpeak更快地收敛到较小的角度。在?peak??peak散布图中，纵轴表示与?=0?=0的偏差。对于a0=1a0=1和5的情况，?peak??peak接近0°，表明辐射在??方向上基本保持对称。然而，当a0=10a0=10时，所有三个数据集都表现出明显的偏差。在如此高的a0a0值下，电子运动进入了强非线性区域，在该区域中横向振荡不再完全对称，导致辐射模式在??方向上发生系统性偏移。在这种情况下，横向运动中出现了谐波和混合频率成分；它们在xx和yy方向上的振幅不同，打破了轨迹在??平面上的镜像对称性。结果，电子轨迹包含非对称成分，发出的辐射干涉在?≠0?≠0处产生一个峰值，表现为可观察到的方位角偏移。第3.3节分析的时间域和频率域辐射是沿着空间功率-角度分布的峰值发射方向进行的，其中θobs = 75°和φobs = 10°。这确保了时间和频谱信号与主导辐射方向相对应。空间功率-角度分布的分析阐明了辐射功率是如何分布在θθ和??上的。ZmaxZmax热图显示，激光振幅a0a0的变化在决定峰值辐射强度方面起主导作用，而光束腰厚b0b0的影响则较为温和，尽管在高a0a0值下这两个参数表现出协同放大效应。θpeakθpeak和?peak?peak的行为进一步揭示了a0a0和b0b0的变化如何影响辐射的方向性和聚焦性。通过适当选择a0a0和b0b0的组合，可以有效改善脉冲辐射的方向性和聚焦性。这些发现对于提高辐射捕获系统中的能量利用效率具有指导意义：通过跟踪θpeakθpeak的收缩，可以相应调整探测器位置以提高收集效率。在设计高亮度X/γ射线源时，调节激光强度和光束腰厚可以优化辐射的方向性和收集效率。

第3.3节中分析的时间域和频率域光谱也受到激光振幅a0a0和光束腰厚b0b0不同组合的影响。在所有情况下，随时间变化的辐射功率都是沿着第3.2节中定义的观测方向（θobs = 75°, φobs = 10°）进行评估的。对于每一组参数，横轴代表时间t，纵轴表示单位立体角的瞬时辐射功率。这些时间光谱允许直接视觉比较光束腰厚b0b0和激光振幅a0a0的变化如何影响辐射的峰值强度。从时间光谱可以看出，a0a0主要影响整体辐射强度，但对峰值形状的影响较小，而增加b0b0不仅增强了辐射强度，还修改了发射的时间结构。随着b0b0的增加，辐射从单峰结构演变为三峰结构，如图6所示。这种转变源于在较大的光束腰厚中，电子在其运动过程中多次满足静止相位条件——即电子的速度和加速度暂时与激光相位对齐——从而产生多个辐射峰值。中间峰值对应于在焦点中心附近发出的辐射。在不同b0b0值下，左侧和右侧峰值的强度交替变化，是因为电子在进入和离开焦点区域时遇到了略有不同的激光相位、加速度方向和传播路径。这些差异导致到达观测方向的辐射波之间的相位延迟不同，从而产生了观察到的不对称性。此外，对于b0=6b0=6的情况，观察到辐射强度的异常减少，这是由于三个发射段之间的相对相位差异接近ππ，导致部分破坏性干涉而不是相干增强。

图6显示了不同a0a0和b0b0组合的时间分辨辐射功率，标记为(a–i)。主曲线显示了单位立体角的瞬时辐射功率，而插图提供了峰值的放大视图，以揭示精细的时间结构。通过放大时间光谱，我们观察到主峰和侧峰都是由两个接近的波形逐渐合并形成的。沿不同电子轨迹发射的辐射带有固有的时间延迟，它们在峰值位置的汇聚反映了相干增强，从而形成了一个尖锐的高强度峰值。为了量化这一主发射事件的持续时间，我们测量了其半高全宽（FWHM），得到大约0.160.16阿秒的值。这个脉冲宽度处于阿秒范围内，证实我们的配置能够产生具有特征性阿秒级时间持续时间的超短脉冲。时间光谱通常显示出明显的多峰结构，表明辐射过程具有强烈的相干性。这种时间调制对应于频谱干涉结构，并可能伴随着频域中的带宽重新分布。因此，我们继续分析和频域光谱，以研究不同频率成分中的辐射能量分布。从时间光谱中，我们已经识别出峰位置和峰值强度的参数依赖性变化。为了进一步揭示与这些时间行为相关的频率组成和非线性展宽机制，我们沿着空间功率-角度分布中辐射最大的方向分析光谱。在每个光谱中，横轴代表归一化频率ω/ω0ω/ω0，纵轴表示单位立体角和单位频率的辐射能量密度d2I/dωdΩ，描述了不同频率成分沿观测方向的强度分布。为了清晰起见，我们将光谱分为两组，A组和B组。A组使用原始的纵轴单位来突出辐射能量密度的绝对差异，而B组采用归一化、降采样和平滑处理来强调光谱形状的变化。在激光振幅a0=8a0=8固定的情况下，对于光束腰厚b0=4,6,b0=4,6和1010获得的三个光谱显示出高度一致的整体形状。因此，只有对应于a0=8,b0=6a0=8,b0=6的原始光谱被选为主文中的代表性例子（图7）。该图展示了保留物理纵轴单位的完整能量密度分布，说明在这一参数范围内的辐射光谱具有特征性的单峰结构。

图7表示了a0=8,b0=6a0=8,b0=6的情况的代表性辐射光谱。图显示了保留物理单位的完整能量密度分布，并展示了在这一参数范围内观察到的典型单峰结构。为了量化光束腰厚的影响，从所有三个光谱中提取了高频截止频率和光谱质心，并在表1中进行了总结。结果显示，在所有情况下，截止频率都达到了计算扫描的上限（ω/ω0=8000ω/ω0=8000），表明当a0=8a0=8时，光束腰厚对最大辐射频率的影响较小。然而，随着b0b0的增加，光谱质心略有上升，表明较大的光束腰厚使电子能够在有效相互作用区域内停留更长时间，导致能量向高频侧弱分布。

表1列出了在固定激光振幅a0=8a0=8的情况下，三种光束腰厚b0=4,6,b0=4,6和1010的参数组合。所有情况下的截止频率都达到了扫描范围的上限，而光谱质心随着b0b0的增加而略有上升。在固定光束腰厚b0=3b0=3的情况下，a0=6,8,a0=6,8和1010的三个光谱分别显示在图8中。这些光谱经过降采样和移动平均平滑处理，以保持主要光谱特征的同时提高不同a0a0值之间的可比性。随着激光振幅的增加，主光谱峰值向更高频率移动，伴随着高频尾部的显著扩展。这种行为表明，更强的激光场诱导了电子更强烈的非线性加速，从而产生了更宽的高频辐射带。观察到的展宽趋势与时间域光谱中的主要峰值增强和辐射持续时间的延长一致，证实了a0a0是控制光谱展宽的主要参数。

总结来说，这项工作基于一个涉及圆偏振激光脉冲和高能电子的交叉碰撞模型，我们系统地研究了激光振幅a0a0和光束腰厚b0b0对电子动力学和辐射特性的协同影响。通过控制这些激光场参数，我们对RNTS的机制有了更深入的了解。利用数值模拟，我们研究了电子动力学、发射辐射的空间分布以及发射的时间-频谱特性。除了传统的分析之外，a0a0和b0b0的耦合变化还使我们能够对电子轨迹的空间形态、能量的演化以及辐射的方向性和聚焦特性进行新的观察。此外，这种协调调节被证明可以产生具有阿秒级时间结构的高频辐射信号。本研究的结果为设计高亮度、方向可控的非线性汤姆逊散射源提供了理论指导，并为开发超快X射线和阿秒辐射源提供了潜在价值。尽管如此，当前模型依赖于单电子近似和理想的高斯场；未来的研究应该结合光束能量扩散、聚焦缺陷和辐射-反应效应，以获得更接近实验条件的预测。总体而言，激光振幅a0a0和光束腰厚b0b0的协调控制为理解和优化强场条件下的相对论辐射过程提供了新的视角，并为设计可调的超快辐射源奠定了基础。