[论文解读] Cosmic Ray Transport in the Ionized and Neutral ISM: MHD-PIC Simulations and Effective Fluid Treatments

Christopher J. Bambic, Xue‐Ning Bai+1 more|arXiv (Cornell University)|Feb 23, 2021

Astrophysics and Star Formation Studies参考文献 1被引用 4

一句话总结

本研究首次开展了宇宙射线（CR）在电离星际介质（ISM）与中性云之间传输的自洽MHD-PIC模拟，证明由自激产生的阿尔芬波引起的波粒散射可导致宇宙射线表现出流体行为。模拟结果表明，当动力学模拟所得到的散射速率与拟线性理论和福克-普朗克理论的预测值一致时，宇宙射线流体动力学能准确描述宇宙射线的输运，从而验证了在ISM中使用流体近似方法的合理性。

ABSTRACT

Cosmic rays (CRs) have critical impacts in the multiphase interstellar medium (ISM), driving dynamical motions in low-density plasma and modifying the ionization state, temperature, and chemical composition of higher-density atomic and molecular gas. We present a study of CR propagation between the ionized ISM and a neutral cloud. Using one-dimensional magnetohydrodynamic particle-in-cell simulations which include ion-neutral drag to damp Alfv$\acute{ ext{e}}$n waves in the cloud, we self-consistently evolve the kinetic physics of CRs and fluid dynamics of the multiphase gas. By introducing the cloud in our periodic domain, our simulations break translational symmetry and allow the emergence of spatial structure in the CR distribution function. A negative spatial gradient forms across the fully-ionized ISM region while a positive gradient forms across the neutral cloud. We connect our results with CR hydrodynamics formulations by computing the wave-particle scattering rates as predicted by quasilinear, fluid, and Fokker-Planck theory. For momenta where the mean free path is short relative to the box size, we find excellent agreement among all scattering rates. By exploring different cloud sizes and ion-neutral collision rates, we show that our results are robust. Our work provides a first-principles verification of CR hydrodynamics when particles stream down their pressure gradient, and opens a pathway toward comprehensive calibrations of transport coefficients from self-generated Alfv$\acute{ ext{e}}$n wave scattering with CRs.

研究动机与目标

• 研究宇宙射线在电离ISM向中性云过渡区域的输运行为，其中离子-中性子阻尼会改变波介导的散射过程。
• 检验在存在强空间梯度且由于嵌入云团导致平动对称性破缺的条件下，宇宙射线流体动力学的有效性。
• 通过将动力学模拟与诸如拟线性与福克-普朗克散射速率等流体理论进行比较，校准输运系数。
• 确定在低电离、中性环境中，宇宙射线流体近似在何种条件下失效或仍成立。
• 为未来在星系尺度系统中使用校准后的流体系数进行宇宙射线输运建模，建立第一性原理基础。

提出的方法

• 执行一维、周期性MHD-PIC模拟，引入离子-中性子阻尼以抑制中性云区域的阿尔芬波。
• 通过宇宙射线漂移不稳定性（CRSI）实现自洽波生成，利用粒子-网格方法模拟宇宙射线动力学。
• 计算宇宙射线分布函数及其矩（能量密度与通量），以提取在波粒散射存在下的空间梯度与时间演化特征。
• 利用拟线性理论、福克-普朗克理论以及模拟中CR矩的时空演化，计算波粒散射速率。
• 在不同理论框架之间比较散射速率，评估其与宇宙射线流体动力学的一致性。
• 改变云团尺寸与离子-中性子碰撞率，以检验结果在不同物理参数下的鲁棒性。

实验结果

研究问题

• RQ1 中性云中存在离子-中性子阻尼时，如何破坏平动对称性并改变宇宙射线分布函数的空间结构？
• RQ2 在存在自激阿尔芬波的情况下，宇宙射线输运的动力学模拟与基于流体的宇宙射线流体动力学在多大程度上达成一致？
• RQ3 从模拟矩中推导出的波粒散射速率是否与拟线性和福克-普朗克理论的预测一致？
• RQ4 在低电离、弱导电区域（如分子云）中，宇宙射线流体近似在何种条件下仍保持有效？
• RQ5 能否可靠地使用基于第一性原理的动力学模拟所获得的输运系数，来校准天体物理模型中的宏观宇宙射线流体动力学？

主要发现

• 在电离ISM区域，由于波粒散射驱动的扩散传播，宇宙射线能量密度形成负的空间梯度。
• 在中性云中，由于散射被抑制且通量随空间梯度减小，宇宙射线能量密度出现正的空间梯度。
• 基于CR矩时空演化的散射速率计算结果，与拟线性理论和福克-普朗克理论在平均自由程相对于模拟盒尺寸较短的区域预测值高度一致。
• 三种散射速率计算结果的一致性，为电离ISM中宇宙射线流体动力学提供了第一性原理验证，尤其适用于平均自由程可容纳于ISM区域多个量级的宇宙射线动量。
• 结果在不同云团尺寸与离子-中性子碰撞率下均保持稳健，表明在所研究条件下，流体近似有效性不依赖于这些参数。
• 本研究为利用自激阿尔芬波散射，通过第一性原理方法校准宇宙射线流体动力学中的扩散系数，建立了可行路径，从而可实现对多相ISM中宇宙射线输运更精确的建模。