[论文解读] Modeling of radiative and quantum electrodynamics effects in PIC simulations of ultra-relativistic laser-plasma interaction

Mathieu Lobet, E. d’Humières+4 more|arXiv (Cornell University)|Nov 5, 2013

Laser-Plasma Interactions and Diagnostics参考文献 15被引用 51

一句话总结

本文将辐射和量子电动力学（QED）效应——特别是非线性康普顿散射和布雷特-惠勒对产生——实现于PIC代码Calder中，用于模拟超相对论性激光-等离子体相互作用。该方法结合了低χ区域的经典辐射阻尼与高χ过程的随机蒙特卡罗采样，验证结果表明模拟与解析解高度一致，证明了在超过10²³ W/cm²强度下对对产生级联和能量损失动力学的精确建模。

ABSTRACT

Next generation of ultra-intense laser facilities will lead to novel physical conditions ruled by collective and quantum electrodynamics effects, such as synchrotron-like emission of high-energy photons and electron-positron pair generation. In view of the future experiments performed in this regime, the latter processes have been implemented into the particle-in-cell code CALDER.

研究动机与目标

• 在粒子-网格（PIC）模拟中对超相对论性激光-等离子体相互作用中的辐射和量子电动力学（QED）效应进行建模。
• 在强激光场中实现高能光子发射（通过非线性康普顿散射）和$e^{+}e^{-}$对产生（通过布雷特-惠勒过程）的精确模拟。
• 通过在电子在恒定磁场中辐射的基准配置中与解析解对比，验证所实现模型的准确性。
• 通过提供一个物理解释准确的模拟框架，为下一代拍瓦和艾瓦级激光装置的未来实验计划提供支持。

提出的方法

• 扩展PIC代码Calder，引入用于光子和对产生生成的QED模型，使用超粒子表示电子和光子群体。
• 在经典区域，通过Sokolov方案建模辐射阻尼，辐射功率由量子发射率（公式2）计算得出。
• 在量子区域，采用蒙特卡罗算法，基于发射概率的光学深度演化与累积分布函数（公式6–7）采样光子发射。
• 对产生通过类似蒙特卡罗方案应用于超光子实现，通过聚并方案保持能量与动量守恒。
• 模拟中跟踪量子参数$\chi_e$、$\chi_\gamma$和$\chi_\pm$，以区分经典与量子区域。
• 采用聚并方案，在光子与对产生过程中平均保持能量与动量守恒，从而提高数值效率。

实验结果

研究问题

• RQ1 如何在超相对论性激光-等离子体相互作用的PIC模拟中精确建模辐射和QED效应，如非线性康普顿散射与布雷特-惠勒对产生？
• RQ2 光子发射与对产生中，经典与量子区域之间的过渡行为如何？该行为在模拟中如何被捕捉？
• RQ3 在高χ环境中，辐射损失与对级联如何随时间演化？它们对电子与光子能量谱有何影响？
• RQ4 在基准配置中，QED过程的蒙特卡罗实现能在多大程度上复现解析解？

主要发现

• PIC模拟结果与Fokker-Planck型方程（9）和（10）的解析解高度一致，验证了QED过程实现的正确性。
• 在量子区域（$\chi_e \simeq 5$），最初发射的光子能量接近初始电子能量（$\gamma_\gamma \sim 10^3$），导致电子发生显著能量损失。
• 初始阶段平均光子量子参数$\chi_\gamma \sim 1$，有利于通过布雷特-惠勒机制实现高效对产生。
• 观察到对级联现象：新生成的粒子通过辐射产生更多对，形成致密的对等离子体。
• 当平均光子能量降至对产生阈值以下（$\chi_\gamma < 10^{-1}$）时，对产生率趋于稳定，辐射损失成为主导的冷却机制。
• 量子区域中QED过程的随机性使电子能量谱展宽，而半经典区域中的经典辐射阻尼则倾向于使其变窄。