[论文解读] Diffusion of relativistic charged particles and field lines in isotropic turbulence: II. Analytical models

Marco Kuhlen, Vo Hong Minh Phan+1 more|arXiv (Cornell University)|Nov 10, 2022

Solar and Space Plasma Dynamics被引用 4

一句话总结

本论文通过引入一个三阶段解析模型——沿磁力线的初始粒子输运、磁力线扩散以及粒子去相关——解决了各向同性湍流中垂直宇宙射线输运领域长期存在的理论模型与模拟结果之间的差异。该研究指出，磁力线输运的亚扩散行为是垂直扩散系数中非标准刚度标度的根源，且模拟结果与模型预测在所有刚度范围内（包括标准理论失效的低刚度区域）表现出极佳的一致性。

ABSTRACT

The transport of high-energy particles in the presence of small-scale, turbulent magnetic fields is a long-standing issue in astrophysics. Analytical theories on transport perpendicular to the large-scale magnetic field disagree with numerical simulations at rigidities where the particles' gyroradii are slightly smaller than the correlation length of turbulence. At the same time, extending the numerical simulations to lower rigidities has proven computationally prohibitive. We present an analytical model for the perpendicular transport, based on (1) initial particle transport along field lines, (2) the transport of field lines and (3) the eventual decorrelation of particles from field lines. Transport parallel to the large-scale field is governed by pitch-angle scattering and so for times larger than the inverse pitch-angle diffusion coefficient, particles spatially diffuse in the parallel direction. Our results suggest that perpendicular diffusion occurs when particles have displaced in the perpendicular direction by a few correlation lengths of turbulence. We have tested the analytical theory by running a large suite of test particle simulations at unprecedentedly low rigidities, making extensive use of graphical processing units (GPUs). Our numerical results exhibit a non-standard rigidity-dependence for the perpendicular diffusion coefficient at intermediate rigidities. At the lowest rigidities, the standard rigidity-dependence is recovered. The simulated diffusion coefficients are nicely reproduced by our analytical model. We have traced the non-standard rigidity-dependence to a subdiffusive phase in the field line transport. Our study confirms our understanding of the escape of cosmic rays from the Galactic halo and its rigidity-dependence. [abridged]

研究动机与目标

• 解决在中等刚度下，解析模型与数值模拟在垂直宇宙射线输运中长期存在的不一致问题。
• 探究在模拟中观察到但标准理论未能预测的垂直扩散系数中非标准刚度标度的物理起源。
• 开发一个能准确复现宽刚度范围内模拟结果（包括低刚度区域）的解析模型。
• 阐明磁力线输运与宇宙射线垂直扩散之间的物理解释联系，特别是亚扩散磁力线运动的作用。
• 为湍流天体等离子体中的宇宙射线输运提供一个物理解释一致的建模框架，适用于日球层和星际介质。

提出的方法

• 利用GPU加速计算进行大规模测试粒子模拟，以达到前所未有的低刚度范围（r_g/L_c ≲ 10⁻³），突破以往模拟的极限。
• 将垂直扩散系数建模为三阶段过程：(1) 沿磁力线的初始弹道输运，(2) 受湍流磁场结构控制的磁力线扩散，(3) 粒子最终与磁力线去相关。
• 采用半解析的常微分方程（ODE）方法建模随时间演化的磁力线扩散系数，在Corrsin假设和各向同性Kolmogorov湍流条件下求解σ_x²(z)与d_FL(z)的耦合方程。
• 对湍流磁场采用Kolmogorov k⁻⁵/³功率谱，外尺度为L_max，关联长度L_c ≈ L_max/5。
• 应用磁静近似，假设粒子速度超过阿尔芬速度，并对湍流磁场实现与粒子初始方向进行系综平均。
• 通过比较不同湍流强度（η = 0.2, 0.5）和刚度下的运行扩散系数与渐近扩散系数，验证解析模型与模拟结果的一致性。

实验结果

研究问题

• RQ1 为何标准解析模型无法复现数值模拟在中等刚度下观察到的刚度依赖的垂直扩散系数？
• RQ2 在各向同性湍流中，垂直扩散系数非标准刚度标度的物理机制是什么？
• RQ3 磁力线输运如何随时间演化？其是否表现出可解释异常粒子扩散的亚扩散行为？
• RQ4 通过将过程分解为磁力线扩散与粒子去相关，能否准确建模粒子垂直扩散系数？
• RQ5 基于磁力线动力学的解析模型是否能在整个刚度范围内（包括标准理论适用的低刚度区域）准确复现模拟结果？

主要发现

• 在中等刚度下（r_g/L_c ≲ 1），模拟得到的垂直扩散系数表现出非标准刚度依赖性，λ_⊥/λ_∥先上升后在高刚度下急剧下降。
• 在最低刚度下（r_g/L_c ≲ 10⁻³），恢复了标准刚度标度（λ_⊥ ∝ r_g²），与理论预期一致。
• 非标准刚度依赖性源于磁力线输运中的亚扩散阶段，其中d_FL(z)最初增速超过线性，随后出现过冲与振荡，最终趋于恒定渐近值。
• 基于磁力线ODE与粒子去相关的解析模型，在所有刚度与湍流水平下均能高精度复现模拟得到的扩散系数。
• ODE模型成功捕捉了运行磁力线扩散系数中的振荡过冲行为，该现象源于积分中的干涉项（cos(k cosθ z)），并准确再现了从弹道到扩散的磁力线运动转变。
• 通过施加微小的重标度因子（η=0.2时为0.8，η=0.5时为0.56），模拟与解析模型的一致性进一步提升，表明尽管振荡相位与周期存在微小差异，模型仍具稳健性。