[论文解读] Modeling of photon and pair production due to quantum electrodynamics effects in particle-in-cell simulation

Wenqi Wang, Z. M. Sheng+2 more|arXiv (Cornell University)|Aug 23, 2016

Quantum Electrodynamics and Casimir Effect被引用 3

一句话总结

本文提出KLAPS，一种通过引入量子电动力学（QED）效应增强的粒子-网格（PIC）代码，用于模拟超强激光-等离子体相互作用中通过康普顿散射产生的光子生成以及通过布雷特-惠勒过程产生的电子-正电子对。该方法采用两种等价的QED公式计算辐射速率，并结合蒙特卡罗事件生成器，模拟结果与理论预测及先前研究的基准数据高度一致，验证了其在QED主导等离子体模拟中的可靠性。

ABSTRACT

We develop the particle-in-cell (PIC) code KLAPS to include the photon generation via the Compton scattering and electron-positron creation via the Breit-Wheeler process due to quantum electrodynamics (QED) effects. We compare two sets of existing formulas for the photon generation and different Monte Carlo algorithms. Then we benchmark the PIC simulation results.

研究动机与目标

• 开发一种能够模拟超强激光-等离子体相互作用中量子电动力学（QED）效应的粒子-网格（PIC）代码。
• 实现基于康普顿散射的光子生成和基于布雷特-惠勒过程的粒子对生成的精确模型。
• 将模拟结果与理论公式及文献中已有的QED-PIC结果进行基准测试。
• 验证不同蒙特卡罗算法在不同时间步长和粒子密度下的事件生成准确性与统计收敛性。
• 实现对下一代10-PW和100-PW激光装置中预期出现的QED主导区域的可靠模拟。

提出的方法

• KLAPS代码采用两种等价的光子发射速率公式来引入QED效应：一种基于Elkina和Nerush（公式1），另一种基于Erber和Kirk（公式4），两者均基于弱场和准稳态近似推导得出。
• 光子生成通过蒙特卡罗事件生成器建模，总发射速率决定每个时间步长内光子发射的概率，能量分布通过累积分布函数的逆变换采样获得。
• 对于粒子对生成，采用相同的理论框架，使用粒子对生成速率公式（公式6），并假设电子与正电子具有对称的发射行为。
• 实现了三种不同的蒙特卡罗事件生成方法并进行比较：(I) 使用均匀随机数的拒绝采样，(II) 使用指数分布的累积概率法，(III) 改进的累积条件法。
• 模拟采用128×192的网格，每个单元一个电子，时间步长自适应调整以确保$W \times Dt \leq 0.2$，保证数值稳定性。
• 基准测试基于辐射和粒子对生成速率的解析解，以及Elkina等人发表的级联模拟结果。

实验结果

研究问题

• RQ1 在高QED参数区域（$\chi_e \gtrsim 1$）下，PIC代码对康普顿散射过程中的光子发射速率模拟精度如何？
• RQ2 现有两种QED辐射速率理论公式（Elkina-Nerush与Erber-Kirk）中，哪一种与模拟结果更吻合？
• RQ3 不同蒙特卡罗事件生成算法如何影响光子与粒子对生成速率的统计收敛性与准确性？
• RQ4 QED-PIC模拟高能过程时，需要多高的时间步长分辨率以保证精度？
• RQ5 所开发的QED-PIC代码能否重现强磁场中电子-正电子级联发展的既定基准结果？

主要发现

• 在不同$\chi_e$值下，两种光子发射速率理论公式（公式1与公式4）的计算结果几乎完全一致，验证了其在模拟参数范围内的等价性。
• 经典极限公式（公式5）高估了高能光子的发射速率，证实了在强场区域必须使用完整的QED表达式。
• 所有三种蒙特卡罗事件生成方法（I、II、III）产生的结果在统计上等价，表明该实现对不同采样策略具有鲁棒性。
• 光子生成速率的模拟结果在不同$\chi_e$值范围内与理论预测高度一致，最佳一致性出现在$Dt = 0.2/W$时。
• 当时间步长增加至$Dt = 6/W$时，模拟结果与理论值显著偏离，凸显了小时间步长对精度的重要性。
• 该代码成功重现了Elkina等人发表的基准级联模拟，经过4000次运行后，能量高于100 MeV的粒子对数量与已发表数据在统计不确定性范围内一致。