[论文解读] Anisotropic transport and early dynamical impact of Cosmic Rays around Supernova remnants

Philipp Girichidis, Thorsten Naab+2 more|arXiv (Cornell University)|Jun 18, 2014

Solar and Space Plasma Dynamics被引用 4

一句话总结

本文提出了一种在FLASH磁流体动力学代码中对宇宙射线（CRs）的新型实现方法，将CRs建模为具有能量依赖性各向异性扩散（沿磁场方向与垂直于磁场方向）、平流输运及绝热损失的独立流体。研究发现，CRs引致的加速度强度远超热速度，最高可达热值的100倍，主要发生在垂直于磁场方向；在80 pc尺度上，持续约150 kyr的时间内，1 GeV能量处的峰值注入谱比幂律谱产生更显著的动力学影响。

ABSTRACT

We present a novel implementation of cosmic rays (CR) in the magneto-hydrodynamic code FLASH. CRs are described as separate fluids with different energies. CR advection, energy dependent anisotropic diffusion with respect to the magnetic field and adiabatic losses to follow the evolution of spectra are taken into account. We present a first study of the transport and immediate (~150 kyr) dynamical impact of CRs on the turbulent magnetised interstellar medium around supernova remnants on scales up to 80 pc. CR diffusion quickly leads to an efficient acceleration of low-density gas (mainly perpendicular to the magnetic field) with accelerations up to two orders of magnitude above the thermal values. Peaked (at 1 GeV) CR injection spectra have a stronger impact on the dynamics than power-law spectra. For realistic magnetic field configurations low energy CRs (with smaller diffusion coefficients) distribute anisotropically with large spatial variations of a factor of ten and more. Adiabatic losses can change the local spectra perceptibly but do not have an integral effect on the dynamics at the spatial and temporal scales considered here. We discuss the potential global impact of CRs and anisotropic transport on the dynamical structure of the ISM and also detail the limitations of the model.

研究动机与目标

• 在磁流体动力学（MHD）代码中开发并实现一种物理解释一致的宇宙射线（CRs）建模方法，作为独立流体，考虑能量依赖的输运特性。
• 研究宇宙射线在约150 kyr内的早期动力学影响，作用于超新星遗迹（SNRs）周围尺度达80 pc的湍流磁化星际介质（ISM）。
• 评估CR扩散各向异性、谱形（峰值谱与幂律谱）以及绝热损失对CR压力梯度和气体加速度的影响。
• 评估CR能量输运的数值鲁棒性，特别是绝热损失的可逆性以及谱演化过程中总能量的守恒性。

提出的方法

• 将CRs作为具有能量依赖输运特性的独立流体处理，采用多能段能量分布以解析谱演化过程。
• 通过流体方程建模CR平流输运，并引入与粒子能量相关的平行和垂直于磁场方向的各向异性扩散系数。
• 采用基于通量的方法在能段之间实现绝热能量损失，确保能量守恒达到机器精度。
• 使用有限体积MHD求解器（FLASH），结合自洽的磁场演化，模拟超新星遗迹膨胀及CR反馈在湍流ISM中的作用。
• 对CR谱采用分段常数分布函数，并在压缩流条件下测试绝热损失的可逆性。
• 开展受控的数值测试，包括单格点超新星爆炸和正弦压缩模拟，以验证能量守恒与谱演化特性。

实验结果

研究问题

• RQ1 各向异性CR扩散（沿磁场方向较强，垂直方向较弱）如何影响CR与周围ISM之间的动力学耦合？
• RQ2 在超新星爆发后约150 kyr的早期阶段，1 GeV处的峰值谱与幂律谱在气体加速度方面有何相对动力学影响？
• RQ3 绝热损失在多大程度上改变局部CR谱，并影响80 pc尺度上CR能量向ISM的整体转移？
• RQ4 谱分辨率和能量通量可逆性方面的数值限制如何影响CR能量输运与谱演化的准确性？
• RQ5 CR压力梯度在超新星遗迹激波附近低密度气体区域中，如何驱动早期非热加速度？

主要发现

• 各向异性CR扩散导致气体加速度达到热速度的两倍以上，主要发生在垂直于磁力线的方向。
• 以1 GeV为中心的峰值CR注入谱相比幂律谱，对ISM产生显著更强的动力学影响，原因在于其更强的压力梯度。
• 低能CR表现出高度各向异性的空间分布，其变化幅度可达十倍以上，主要由较小的扩散系数驱动。
• 绝热损失会明显改变局部CR谱，但在所研究的空间和时间尺度上，对CR在ISM中的整体动力学影响影响甚微。
• 数值实现中，总CR能量守恒达到机器精度，但当压缩流速超过约10⁻¹² s⁻¹时，绝热损失的可逆性开始失效，该速率接近超新星遗迹激波中的物理观测范围。
• 单格点超新星注入模拟表明，若未解析内部激波结构，则不会出现谱硬化现象，凸显了准确解析激波结构对CR谱演化模拟的重要性。