Correlating Particle Acceleration Rates with Plasma Conditions in Colliding Wind Binaries

技術要約：衝突風連星における粒子加速率とプラズマ条件の相関

問題提起
超新星残骸（SNe）は、従来、$\sim$1 TeVまでの銀河宇宙線の主要な加速源であると考えられてきたが、LHAASO実験による最近の観測では、超新星とは関連のない星形成領域（SFRs）から広域の超高エネルギー（VHE）放射が特定されている。衝突風連星（CWBs）は、これらの粒子をTeVおよびPeVエネルギーまで加速できる潜在的な光源として提案されている。しかし、CWBsの複雑で乱流的な環境において、この加速を駆動する具体的な物理メカニズムは依然として不明である。先行研究では、CWBsが非熱的放射を生み出すことは確立されているが、古典的な拡散衝撃波加速（DSA）と、乱流や磁気駆動プロセス（確率的加速、磁気リコネクション、ミラー効果など）による寄与を十分に分離できていない。本研究は、どのプラズマ条件および物理量が粒子加速率と最も強く相関するかを特定し、CWBsにおける支配的な加速メカニズムを明らかにすることを目的としている。

手法
著者らは、$\eta$ Carinae連星の近星点通過時をモデルとしたCWBシステムの高解像度3次元磁気流体力学（MHD）シミュレーションを実施した。シミュレーションの設定は以下の通りである：

• MHDモデリング: システムは、AMUNコードを用いて、異なる風パラメータ（速度はそれぞれ500 km s$^{-1}$および2500 km s$^{-1}$）を持つ主星と副星からなる理想圧縮性MHD方程式を用いて発展させた。シミュレーションでは、有効ポリトロピック指数（$\Gamma_{\text{eff}}$）を断熱的（5/3）から強力な冷却状態（1.1）まで変化させることで、異なるプラズマ圧縮性レジームを探索した。
• 解像度: 適応格子細分化法（AMR）を採用し、衝撃波、乱流、および磁気リコネクション構造を正確に捉えるために、最大で$960 \times 1920 \times 1920$セルの有効解像度を達成した。
• テスト粒子積分: MHD発展の後、$10^5$個のテスト粒子（計算効率のために主にFe$^+$を用い、結果はH$^+$に一般化）を乱流領域に注入した。粒子の軌道は、シミュレーションされた速度場および磁場に結合した相対論的運動方程式を用いて積分された。
• 相関分析: 本研究では、局所的な粒子加速率（$dE/dt$）を、シミュレーションのスナップショットから抽出された4つの特定のプラズマ量、すなわち磁場強度（$|B|$）、電流密度（$|J|$）、渦度（$|\nabla \times \mathbf{u}|$）、および速度ダイバージェンス（$|\nabla \cdot \mathbf{u}|$）と相関させることに焦点を当てた。統計解析は、加速率が100 eV/sを超えた領域のみを対象として、2次元確率密度関数（2D-PDF）およびピアソンの相関係数を用いて行われた。

主な結果

1. 加速メカニズム: 解析の結果、粒子加速は古典的なDSAよりも、乱流および磁気駆動プロセスによって支配されていることが明らかになった。

   • 磁場と渦度: 加速率と、磁場強度（$|B|$）および渦度（$|\nabla \times \mathbf{u}|$）との間に強い正の相関が見られた。ピアソンの相関係数は、これらの変数に対して0.46から0.55の範囲であった。回帰勾配は超線形（最大1.79）であり、より強い磁場と高い乱流がエネルギー獲得を著しく増強することを示している。
   • 速度ダイバージェンス: 対照的に、加速率と速度ダイバージェンス（衝撃圧縮のプロキシ）との相関は弱から存在しない程度であり（ピアソン係数0.18–0.33、負またはゼロに近い勾配）、明確な平面状の衝撃波が加速において果たす役割は限定的であることを示唆している。
   • 電流密度: 電流密度は中程度の相関（0.44–0.51）を示したが、回帰勾配はより低かった。著者らはこれを、抵抗効果が理想MHD極限において無視できるため、直接的な駆動因子ではなく、磁場の複雑性の間接的なトレーサーであると解釈している。

2. エネルギー分布:

   • 粒子は数時間から数日のうちに熱的分布から相対論的冪乗則へと進化し、数百TeVのエネルギーに達した。
   • 最大エネルギーおよびエネルギー・スペクトルの「膝（knee）」の位置は、ヒルスの閉じ込め基準と一致して、表面磁場強度（$B_*$）に依存する。
   • ポリトロピック指数（$\Gamma_{\text{eff}}$）はスペクトルの形状に大きな影響を与えた。冷却がより強いモデル（より低い$\Gamma_{\text{eff}}$）は、乱流と磁場の複雑さが増すことで、より硬いスペクトルと高い加速効率を生み出した。断熱モデル（$\Gamma_{\text{eff}} = 5/3$）は、より急な勾配と低い加速効率を示した。

3. 空間分布: 最高レベルの加速が発生する領域は、最強の衝撃波面（速度ダイバージェンス）と厳密には一致せず、むしろ磁場が増幅され複雑化している乱流の下流領域と一致していた。

意義と主張
本論文は、定常状態の仮定を超えて、粒子加速効率と局所的なプラズマ条件との間の決定的な結びつきを提供すると主張している。主な意義は、乱流と磁場の複雑さが粒子加速の主要な駆動因子であり、古典的な拡散衝撃波加速を凌駕していることを示した点にある。

著者らは以下のように結論付けている：

• CWBsは、適度な磁化（$>1$ G）があれば、ハドロン粒子を数百TeVまたはPeVエネルギーまで加速できる。
• 加速プロセスは本質的に確率的かつ磁気駆動型であり、単純な衝撃波交差ではなく、乱流構造（アルフベン波、磁気ミラー、リコネクション）との相互作用に依存している。
• これらの知見は、将来の高感度ガンマ線観測（LHASTやCTAOなど）の解釈に直接的な影響を与える。すなわち、CWBsが宇宙線源として検出可能かどうかは、冷却効率や軌道位相によって変調される、乱流のレベルや磁場の複雑さに大きく依存することを示唆している。

本研究は新しい実験セットアップを提案するものではなく、既存および将来の観測データを解釈するための理論的枠組みを提供するものであり、一部のCWBsにおける強力なガンマ線検出の欠如は、加速能力の欠如ではなく、乱流や磁化のレベルが低いことによるものである可能性を強調している。