[论文解读] Accurate numerical simulation of radiation reaction effects in strong electromagnetic fields

N. Elkina, A. M. Fedotov+2 more|arXiv (Cornell University)|Jan 30, 2014

Superconducting Materials and Applications参考文献 5被引用 3

一句话总结

本文提出了一种基于隐式龙格-库塔-尼斯特伦（RKN）方法求解协变兰道-利夫希茨方程的高精度数值求解器，用于强电磁场中的辐射反作用问题。结果表明，四阶、六阶和八阶隐式配点RKN方法在保持质量壳条件和能量守恒方面显著优于显式方法，从而实现了对强激光场中电子动力学的精确长期模拟，应用结果显示由于辐射反作用，粒子穿透焦点区域的程度降低。

ABSTRACT

The Landau-Lifshitz equation provides an efficient way to account for the effects of radiation reaction without acquiring the non-physical solutions typical for the Lorentz-Abraham-Dirac equation. We solve the Landau-Lifshitz equation in its covariant four-vector form in order to control both the energy and momentum of radiating particle. Our study reveals that implicit time-symmetric collocation methods of the Runge-Kutta-Nyström type are superior in both accuracy and better maintaining the mass-shell condition than their explicit counterparts. We carry out an extensive study of numerical accuracy by comparing the analytical and numerical solutions of the Landau-Lifshitz equation. Finally, we present the results of simulation of particles scattering by a focused laser pulse. Due to radiation reaction, particles are less capable for penetration into the focal region, as compared to the case of radiation reaction neglected. Our results are important for designing the forthcoming experiments with high intensity laser fields.

研究动机与目标

• 开发一种在强电磁场中模拟辐射反作用的数值稳定且高精度的方法。
• 克服洛伦兹-亚伯拉罕-狄拉克方程的局限性，如发散解和因果性破坏问题。
• 通过采用隐式、时间对称、基于配点的RKN方法，改进显式积分器，以实现二次不变量的精确守恒。
• 在恒定磁场和平面波场中，通过与解析解对比验证该方法的准确性。
• 将求解器应用于实际场景，如聚焦激光脉冲中的电子散射。

提出的方法

• 将兰道-利夫希茨方程以协变四维矢量形式表述，以确保能量与动量守恒的一致性。
• 利用高斯求积节点推导出四阶、六阶和八阶的隐式配点高斯-勒让德龙格-库塔-尼斯特伦（RKN）方法。
• 通过配点构造方法，确保时间对称性以及对闵可夫斯基范数等二次不变量的精确守恒。
• 将数值格式与恒定磁场中粒子运动及圆偏振平面波中的解析解进行对比验证。
• 建议采用自适应时间步长控制，以减轻强场梯度区域的精度损失。
• 将该方法应用于聚焦激光束中相对论性电子散射的模拟，对比考虑和不考虑辐射反作用的结果。

实验结果

研究问题

• RQ1 隐式、时间对称的RKN方法是否能在强场中提供优于显式方法的精度与守恒特性？
• RQ2 在保持质量壳条件和能量守恒方面，隐式RKN方法与显式方法相比表现如何？
• RQ3 辐射反作用对聚焦激光脉冲中电子散射动力学有何影响？
• RQ4 在强梯度场区域，数值精度与碰撞参数之间存在何种关联？
• RQ5 采用隐式RKN积分的兰道-利夫希茨方程能否准确模拟出与半经典蒙特卡罗模拟相当的辐射反作用效应？

主要发现

• 四阶、六阶和八阶隐式配点RKN方法在求解兰道-利夫希茨方程时，相比显式方法展现出显著更高的精度。
• 这些隐式方法在长时间积分过程中能更有效地保持质量壳条件 $u_{\nu}u^{\nu} = c^2$。
• 在聚焦激光脉冲中电子散射的模拟中，辐射反作用导致粒子穿透焦点区域的深度低于无辐射反作用的情况。
• 在强场梯度区域，数值精度会下降，尤其对较大碰撞参数的情况更为明显，表明需要采用自适应时间步长。
• 兰道-利夫希茨求解器的结果与半经典蒙特卡罗模拟中硬光子发射的特征具有定性一致性。
• 采用高阶隐式RKN方法可在保持高精度的同时使用更大的时间步长，尽管单步计算成本较高，但整体计算效率得到提升。