The Effects of Cosmic Ray Protons on Galactic Nonthermal Filaments

技術要約：銀河中心の非熱的フィラメントに対する宇宙線陽子の影響

問題提起
銀河中心（GC）には、シンクロトロン放射によって照らされた、長さ数十パーセクに及ぶ非熱的フィラメント（NTF）の集団が存在する。その物理的な起源は未解決である。主に2つの注入メカニズムが提案されている。(1) NTFはパルサー風星雲からのジェットによって供給され、レプトン（電子・陽電子）が支配的な粒子群となる。(2) NTFは星間衝撃波から加速された粒子によって供給され、陽子が支配的な粒子群（陽子とレプトンの比率は約100:1）となる。銀河系の宇宙線集団は、NTFの観測結果と一致する冪乗則のエネルギー分布を示しているが、注入源とフィラメントの形態との間の具体的な因果関係は不明である。さらに、中央分子領域（CMZ）内の極方向磁場トポロジーは複雑であり、多くのフィラメントが非線形であることは、星間乱流の役割を示唆している。本研究は、銀河中心の様々な条件下における宇宙線の伝播、加熱、およびプラズマ流をシミュレーションすることにより、レプトン支配シナリオと陽子支配シナリオを区別することを目的としている。

手法
著者らは、放射および衝突損失項を組み込むように修正した磁気流体力学（MHD）コード Athena++ を用い、フィラメントを表す1次元フラックスチューブ内での宇宙線伝播をシミュレートした。モデルは、フィラメントの中点への宇宙線エネルギーの注入によって摂動を受けた一様な背景媒質をモデル化している。

• 物理プロセス： モデルには、拡散（アルフベン波による散乱）とストリーミング（アルフベン速度による対流）が含まれる。レプトンの損失メカニズムには、シンクロトロン放射、逆コンプトン散乱、クーロン損失、および制動放射が含まれる。二次生成については、一次陽子が周囲の陽子と相互作用して二次電子および陽電子を生成する過程をモデル化している。
• パラメータ： 本研究では、主要なパラメータである磁場強度（$B$）、プラズマ密度（$n_p$）、宇宙線拡散係数（$\kappa$）、および総注入エネルギー（$E_{\pm}$）を変化させている。
• シミュレーション・セット： 各パラメータ・セットについて、レプトンのみの場合と、陽子・レプトン混合の場合（陽子とレプトンの比率を考慮して総注入率を100倍にする）の2回ずつシミュレーションを行った。基準となるパラメータは、$B = 200\,\mu\text{G}$、$n_p = 100\,\text{cm}^{-3}$、$\kappa = 1 \times 10^{26}\,\text{cm}^2\,\text{s}^{-1}$ とした。
• 解析： 著者らは、空間エネルギー密度プロファイル、誘起ガス流、加熱率、および結果として得られるシンクロトロン輝度とスペクトル指数を、時間経過（最大40 kyrまで）とともに解析した。

主な結果

1. 注入メカニズムの識別不能性： シミュレーションの結果、陽子支配の場合とレプトン支配の場合の間で、観測可能な差異はほとんどないことが明らかになった。陽子の注入は、より高い総エネルギー密度によって、より強いガス流と加熱を引き起こすが、これらの効果は検出可能な閾値を下回っている。検討された最低密度においても、誘起速度（$\sim 2.5\,\text{km s}^{-1}$）や加熱率は、放射冷却や典型的な熱線幅と比較して無視できる程度である。
2. 伝播の制約： NTFの観測される平坦で均一な表面輝度は、典型的な銀河系の値（$\kappa \sim 10^{28}\,\text{cm}^2\,\text{s}^{-1}$）よりも著しく低い拡散係数を示唆している。著者らは、観測された輝度プロファイルを再現するには、より低い拡散係数（$\kappa \sim 10^{26}\,\text{cm}^2\,\text{s}^{-1}$）が必要であることを発見した。
3. フィラメントの長さと年齢： シミュレーションは、伝播を遮断するメカニズムがない限り、フィラメントは無限に拡大し続けることを示している。したがって、NTFの長さはその年齢に依存する。放射および衝突による損失は、二次生成を組み合わせても定常状態の長さを生み出さない。なぜなら、二次レプトンは注入サイトから離れた場所で生成され、一次粒子と同様に拡散するためである。
4. スペクトル指数の変動： 本研究では、観測されるシンクロトロン・スペクトル指数（$-2$から$+0.5$ の範囲）の変動メカニズムを調査した。
   • エネルギー依存拡散： 拡散係数がエネルギーに応じてスケールする場合（$\kappa \propto \gamma^q$）、スペクトル指数はフィラメントに沿って変化し、一般に注入サイトから離れるにつれて硬くなる。
   • シンクロトロン老化： 高エネルギー粒子はエネルギーをより早く失うため、時間の経過とともにスペクトルが急峻になる（ステープニング）。
   • 二次生成： 後期においてスペクトルを急峻にする。
   • これらのメカニズムを検討したものの、著者らは、これらがいずれも、実際のNTFで観察されるギザギザとした非対称なスペクトル指数の変動を完全には説明できないことを指摘している。

意義および結論
本研究の主要な貢献は、大規模な宇宙線陽子の存在が、現在の1次元MHDシミュレーションに基づいた現在の観測技術では、パルサー駆動型（レプトン）か衝撃波駆動型（陽子）かの起源を区別できるような、形態、流れ、または加熱における明確で観測可能な特徴を生み出さないことを示した点にある。

したがって、著者らは、これら2つの標準的な注入モデルを区別できないこと、およびNTFに見られる構造的な不完全性（キンク、分裂、非線形な整列）を考慮すると、第三の形成メカニズムを検討する必要があると主張している。本論文は、NTFが銀河中心の乱流における間欠的な構造から生じている可能性を提案している。このシナリオでは、不均質なMHD乱流が磁場を増幅し、渦を薄いフィラメントへと引き延ばすことで、特定の局所的な宇宙線注入源を必要とせずに、観察される構造を作り出す。著者らは、NTFの起源をより深く理解するために、今後の研究はこれらの乱流効果のシミュレーションに焦点を当てるべきであると結論付けている。