[论文解读] Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

José Ortuño-Macías, Krzysztof Nalewajko+5 more|arXiv (Cornell University)|May 4, 2022

Astrophysics and Cosmic Phenomena参考文献 73被引用 11

一句话总结

本研究利用三维动力学模拟，研究具有环向磁场的圆柱形、相对论性、磁化电子-正电子喷流中的不稳定性与粒子加速。通过混合参数（fmix）平衡环向磁场与轴向磁场及气体压强，结果表明：扭曲模（kink modes）主导不稳定性增长，粒子加速主要由快速磁能耗散阶段的垂直电场驱动，并提出适用于浅环向磁场剖面的归一化能量极限 γlim × (RBφ/R0)。

ABSTRACT

Relativistic magnetized jets, such as those from AGN, GRBs and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and non-thermal radiation. Here we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron-positron pair plasma. Generalizing recent treatments by Alves et al. (2018) and Davelaar et al. (2020), we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. (2018) is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

研究动机与目标

• 研究压强与磁场平衡在驱动相对论性、圆柱形、磁化喷流不稳定性中的作用。
• 确定轴向磁场与气体压强之间混合力平衡（即力平衡）对电流驱动与压强驱动不稳定性增长的影响。
• 研究相对论性对等离子体动力学模拟中粒子加速的效率与机制。
• 评估环向磁场剖面（幂律指数 αBφ）对不稳定性增长与最大粒子能量的影响。
• 通过探索不同 fmix 的中间构型，弥合此前在无力场（ffix=0）与 Z-箍缩（ffix=1）极限下的研究结果。

提出的方法

• 对静态、具有圆柱对称性的相对论性对等离子体喷流执行三维粒子-网格（PIC）动力学模拟。
• 采用径向环向磁场 Bφ(r) ∝ r^αBφ 的剖面，其中 −1.5 ≤ αBφ ≤ 0，并包含内、外截断。
• 通过轴向磁场 Bz(r) 与气体压强 P(r) 的组合，利用 fmix ∈ [0,1] 参数化平衡环向磁场力。
• 模拟从无力场（fmix=0）到 Z-箍缩（fmix=1）的构型，局部磁化度最高达 σhot ≈ 8。
• 追踪不稳定性演化、磁通量耗散、电流密度结构与粒子能量谱的变化。
• 采用 Zeltron PIC 代码，包含相对论性、无碰撞等离子体动力学，并在后续扩展中引入辐射冷却效应。

实验结果

研究问题

• RQ1 混合参数 fmix（控制轴向磁场与气体压强之间的平衡）如何影响相对论性喷流中不稳定性增长？
• RQ2 环向磁场剖面的陡峭程度（αBφ）在决定最大粒子能量与快速磁能耗散阶段持续时间方面起什么作用？
• RQ3 在动力学模拟中，哪种电场分量（平行或垂直于 B）主导粒子加速？
• RQ4 对于平坦 Bφ(r) 剖面（αBφ > -1）的喷流外区，是否如线性理论预测的那样出现局域化箍缩模？
• RQ5 对于非均匀磁场剖面，粒子能量极限 γmax 是否可超越 Hilla 标准 γlim = eB0R0/mc2 的一般化？

主要发现

• 对于所有 fmix < 1 的情况，m=1 扭曲模主导不稳定性增长，其线性增长时间常数 τmin 随 fmix 增加而系统性减小。
• 当 fmix = 1（Z-箍缩）时，在中等半径处（r ≈ 2R0）出现一种微弱的局域化 m=0 管道模，与 Begelman（1998）的局域色散关系一致。
• 粒子加速主要由垂直于局部磁场的电场驱动，而平行电场仅起次要作用。
• 对于陡峭的环向磁场剖面（αBφ ≤ -1），最大粒子能量达到标准 Hilla 限制 γlim = eB0R0/mc2。
• 对于平坦剖面（αBφ ≥ -0.5），最大能量按 (RBφ/R0) × γlim 缩放，其中 RBφ 是 Bφ(r) 的峰值半径，表明存在一种重新标定的约束极限。
• 驱动高效粒子加速的快速磁能耗散阶段具有有限持续时间，当扰动传播至外边界时结束，尤其在平坦剖面中更为显著。