[論文レビュー] 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

ひとこと要約

本研究は、イオン化された銀河間媒体（ISM）と中性雲の間を通過する宇宙線（CR）輸送を、初めて自己無撞着なMHD-PICシミュレーションで行った。その結果、自己生成されたアルベール波による波-粒子散乱が、流体的挙動を引き起こすことが示された。シミュレーションは、キネティックシミュレーションで得られた散乱率が準線形理論およびフォッカー・プランク理論の予測と一致する場合に、CRの流体的輸送が正確に記述されることを確認した。これにより、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から中性雲への宇宙線輸送を調査すること。
• 埋め込まれた雲によって破壊される空間的対称性と強い空間勾配を伴う領域において、CR流体的輸送の有効性を検証すること。
• 準線形理論およびフォッカー・プランク理論に基づく散乱率と照らし合わせ、キネティックシミュレーションと比較することで、輸送係数を補正すること。
• 特に低イオン化、弱い導電性を示す環境（例：分子雲）において、宇宙線の流体近似が破綻するか、維持される条件を特定すること。
• 将来的な銀河規模のシステムにおけるCR輸送モデリングのため、補正済みの流体係数を用いた第一原理的基盤を確立すること。

提案手法

• イオン-中性子間の摩擦力によって中性雲領域のアルベール波を減衰させる、一次元的周期的MHD-PICシミュレーションを実施する。
• 宇宙線ストリーミング不安定性（CRSI）による自己整合的な波生成を伴う、粒子法（PIC）を用いて宇宙線のダイナミクスをモデル化する。
• 波-粒子散乱の存在下での空間勾配と時間発展を抽出するために、CR分布関数およびそのモーメント（エネルギー密度およびフラックス）を計算する。
• シミュレーションにおけるCRモーメントの空間的・時間的発展から、波-粒子散乱率を準線形理論、フォッカー・プランク理論、および直接的に算出する。
• 異なる理論的枠組みにおける散乱率を比較し、CR流体的輸送との整合性を評価する。
• 雲のサイズやイオン-中性子衝突率を変化させることで、異なる物理的状態における結果の頑健性を検証する。

実験結果

リサーチクエスチョン

• RQ1 中性雲に起因するイオン-中性子間の減衰が、宇宙線分布関数の空間的構造と並進対称性をどのように破壊するか。
• RQ2 自己生成されたアルベール波が存在する状況下で、キネティックシミュレーションによるCR輸送が、流体的CR流体力学とどの程度一致するか。
• RQ3 シミュレーションのモーメントの時間的・空間的発展から導かれる波-粒子散乱率は、準線形理論およびフォッカー・プランク理論の予測と整合的か。
• RQ4 低イオン化、弱い導電性を示す領域（例：分子雲）において、宇宙線の流体近似がどの条件下で有効に保たれるか。
• RQ5 第一原理的キネティックシミュレーションから得られた輸送係数は、天体物理学的モデルにおけるマクロなCR流体力学の補正に信頼性を持って利用可能か。

主な発見

• 波-粒子散乱による拡散的輸送が駆動するため、イオン化されたISM領域で宇宙線エネルギー密度に負の空間勾配が形成される。
• 散乱が抑制される中性雲領域では、強い散乱が存在しない状況下でフラックスが減少するため、宇宙線エネルギー密度に正の空間勾配が現れる。
• シミュレーションにおけるCRモーメントの空間的・時間的発展から算出された散乱率は、平均自由行程がシミュレーションボックスサイズに比べて短い領域では、準線形理論およびフォッカー・プランク理論の予測と良好に一致する。
• 3つの異なる散乱率の計算結果が一致することは、特にISM領域内に複数の平均自由行程が収容可能な運動量領域の宇宙線に対して、CR流体力学の第一原理的妥当性を裏付ける。
• 雲のサイズやイオン-中性子衝突率を変化させても結果が頑健であるため、研究対象の条件下では流体近似の有効性がこれらのパラメータに敏感でないことが示された。
• 本研究は、自己生成されたアルベール波による散乱を用いた、CR流体力学における拡散係数の第一原理的補正の道筋を確立した。これにより、多相性ISMにおけるCR輸送のより高精度なモデリングが可能になる。