[论文解读] Diffusive Shock Acceleration by Multiple Weak Shocks

Hyesung Kang|arXiv (Cornell University)|Jun 16, 2021

Astrophysics and Cosmic Phenomena参考文献 16被引用 4

一句话总结

本文研究了星系团介质（ICM）中多个弱激波对宇宙射线（CR）的扩散性激波加速（DSA），表明通过 Mach 数 M ≤ 3 的激波反复再加速，可使宇宙射线质子能谱变平，并提升 DSA 效率。关键结果为：由于累积再加速作用，最终 CR 能谱指数趋近于 p⁻³，且与单个激波的压缩比无关，经过三次激波作用后效率显著提高。

ABSTRACT

The intracluster medium (ICM) is expected to experience on average about three passages of weak shocks with low sonic Mach numbers, $M\lesssim 3$, during the formation of galaxy clusters. Both protons and electrons could be accelerated to become high energy cosmic rays (CRs) at such ICM shocks via diffusive shock acceleration (DSA). We examine the effects of DSA by multiple shocks on the spectrum of accelerated CRs by including {\it in situ} injection/acceleration at each shock, followed by repeated re-acceleration at successive shocks in the test-particle regime. For simplicity, the accelerated particles are assumed to undergo adiabatic decompression without energy loss and escape from the system, before they encounter subsequent shocks. We show that in general the CR spectrum is flattened by multiple shock passages, compared to a single episode of DSA, and that the acceleration efficiency increases with successive shock passages. However, the decompression due to the expansion of shocks into the cluster outskirts may reduce the amplification and flattening of the CR spectrum by multiple shock passages. The final CR spectrum behind the last shock is determined by the accumulated effects of repeated re-acceleration by all previous shocks, but it is relatively insensitive to the ordering of the shock Mach numbers. Thus multiple passages of shocks may cause the slope of the CR spectrum to deviate from the canonical DSA power-law slope of the current shock.

研究动机与目标

• 研究星系团介质（ICM）中多次激波通过如何通过反复再加速影响宇宙射线（CR）质子能谱。
• 评估激波 Mach 数变化及绝热膨胀稀释对测试粒子模型中 CR 能谱演化的影响。
• 确定多个激波的累积效应是否可解释射电遗迹激波中 X 射线测得的 Mach 数（M_X）与射电测得的 Mach 数（M_radio）之间的差异。
• 基于 Liouville 定理与先前的离子体模拟，将 DSA 的半解析模型扩展至包含每次激波的膨胀与注入过程。

提出的方法

• 采用半解析方法模拟多次激波通过后的 CR 质子能谱，假设处于测试粒子模型，并假设注入与动量无关。
• 应用 Liouville 定理以考虑激波后区域的绝热膨胀稀释，动量缩放关系为 p′ = Rp，其中 R = (D/r)¹ᐟ³。
• 在每次激波处建模超热粒子的注入，动量满足 p ≥ 3.8p_th，使用各向异性的注入动量 p_inj,k 与归一化因子 f_o,k。
• 同时考虑相同激波（M = 3）与不同 Mach 数序列（如 M = 1.7, 2.0, 3.0）的情形，以检验激波顺序对最终能谱的影响。
• 计算 DSA 效率 η 与 CR 能量比 E_CR,k / E_th,k，以评估非线性效应与反作用潜力。
• 对比 D = 1 与 D = 0.5（即无膨胀与有膨胀）的情形，以分离膨胀对能谱形状的影响。

实验结果

研究问题

• RQ1 在 M ≤ 3 的多个弱激波通过后，CR 质子能谱的幂律指数如何演化？
• RQ2 不同 Mach 数激波的通过顺序是否影响最终 CR 能谱或其幂律指数？
• RQ3 激波之间的绝热膨胀在多大程度上削弱了能谱的变平与放大效应？
• RQ4 多个激波的累积再加速是否可解释射电遗迹中 X 射线 Mach 数（M_X）与射电 Mach 数（M_radio）之间的观测差异？
• RQ5 在 ICM 中，DSA 效率（η）与 CR 能量增益如何随激波通过次数变化？

主要发现

• 在 D = 1 条件下，经过三次相同激波（M = 3）作用后，最终 CR 质子能谱变平至 f(p) ∝ p⁻³，且与单个激波的压缩比无关，与先前理论预期一致。
• DSA 效率 η 从第一次激波后的 7.1×10⁻³ 提升至三次激波后的 0.13，表明通过反复加速实现了显著的 CR 能量增益。
• CR 能量比 E_CR,k / E_th,k 从第一次的 9.6×10⁻³ 上升至第三次的 0.17，证实测试粒子假设仍成立（η < 0.1），适用于 M ≤ 3 的激波。
• 最终 CR 能谱几乎与不同 Mach 数激波的通过顺序无关，且在 p/m_p c ≫ 1 区域，总幂律指数 q_sum(p) 几乎保持恒定。
• 激波后区域的膨胀稀释削弱了能谱的变平与放大效应，D = 0.5 情形下产生的 CR 能量低于 D = 1 情形，原因在于预激波密度降低与动量缩放减弱。
• 多个激波的累积效应主导了最终能谱的指数，其值可能与最近一次激波的 DSA 指数 q_DSA 显著不同，提示 M_radio 可能反映的是激波历史的累积效应，而非单次激波的特性。