[論文レビュー] Relativistic particle transport in extragalactic jets: I. Coupling MHD and kinetic theory

Fabien Casse, Alexandre Marcowith|ArXiv.org|Mar 26, 2003

Astrophysics and Cosmic Phenomena参考文献 35被引用数 10

ひとこと要約

本稿では、銀河間ジェット内の相対論的粒子輸送をシミュレートするための新規な2.5次元MHD-SDE連成手法を提示する。衝撃加速、空間的拡散、およびシンクロtron損失を統合する。その結果、$η_T \approx 0.2-0.3$ の乱流磁場が粒子エネルギーの増大を最大にすることが判明し、100 TeVに達する非累乗法則的スペクトルが得られる。これは初期のスペクトルカットオフにより、X線放射に対するシンクロtronモデルに疑問を呈する。

ABSTRACT

Multidimensional magneto-hydrodynamical (MHD) simulations coupled with stochastic differential equations (SDEs) adapted to test particle acceleration and transport in complex astrophysical flows are presented. The numerical scheme allows the investigation of shock acceleration, adiabatic and radiative losses as well as diffusive spatial transport in various diffusion regimes. The applicability of SDEs to astrophysics is first discussed in regards to the different regimes and the MHD code spatial resolution. The procedure is then applied to 2.5D MHD-SDE simulations of kilo-parsec scale extragalactic jets. The ability of SDE to reproduce analytical solutions of the diffusion-convection equation for electrons is tested through the incorporation of an increasing number of effects: shock acceleration, spatially dependent diffusion coefficients and synchrotron losses. The SDEs prove to be efficient in various shock configuration occurring in the inner jet during the development of the Kelvin-Helmholtz instability. The particle acceleration in snapshots of strong single and multiple shock acceleration including realistic spatial transport is treated. In chaotic magnetic diffusion regime, turbulence levels $η_T=/(B^2+)$ around $0.2-0.3$ are found to be the most efficient to enable particles to reach the highest energies. The spectrum, extending from 100 MeV to few TeV (or even 100 TeV for fast flows), does not exhibit a power-law shape due to transverse momentum dependent escapes. Out of this range, the confinement is not so efficient and the spectrum cut-off above few hundreds of GeV, questioning the Chandra observations of X-ray knots as being synchrotron radiation. The extension to full time dependent simulations to X-ray extragalactic jets is discussed.

研究の動機と目的

• 相対論的粒子輸送をモデル化するため、MHDシミュレーションと確率微分方程式（SDE）を連成する堅牢な数値フレームワークの開発。
• 複雑で時間的に変化するジェット構造における、衝撃加速、空間的拡散、および放射損失（特にシンクロtron損失）の相乗作用の解明。
• 乱流レベル（$\eta_T$）の変動が、銀河間ジェット内の粒子エネルギースペクトルおよび加速効率に与える影響の評価。
• SDEが、現実的な天体物理学的条件下で、拡散・対流方程式の解析的解を正確に再現できるかの検証。
• スペクトルカットオフおよび閉じ込め効果を考慮した場合、シンクロtronモデルがkpcスケールのジェットにおける観測されたX線放射を説明できるかどうかの評価。

提案手法

• ケルビン＝ヘルムホルツ（KH）不安定に起因する乱流を示す銀河間ジェットの2.5次元MHDシミュレーションを用い、乱流の発達と衝撃の形成を捉える。
• 粒子輸送をモデル化するため、確率微分方程式（SDE）を用い、衝撃加速、断熱損失、およびシンクロtron放射損失を組み込む。
• 空間的に変化する拡散係数を扱えるようにSDEを適合させ、拡散・対流方程式の解析的解との整合性を検証する。
• 特に$\eta_T = \langle \delta B^2 \rangle / (B^2 + \langle \delta B^2 \rangle)$ のような混沌とした磁場領域における、横方向の運動量拡散に起因する粒子の脱出を考慮する。
• MHD流れのスナップショットを用いて、キロパーセクスケールのジェットに適した衝撃構造および乱流レベルを焦点に、シミュレーションフレームワークを適用する。
• テスト粒子条件を仮定し、相対論的粒子がMHD流れおよび磁場構造に逆作用を及ぼすことを無視する。

実験結果

リサーチクエスチョン

• RQ1 乱流磁場強度（$\eta_T$）は、銀河間ジェット内の相対論的電子が達成可能な最大エネルギーにどのように影響を与えるか？
• RQ2 空間的に変化する拡散係数および横方向の粒子脱出は、エネルギースペクトルが累乗法則的形状からどれほど歪められるか？
• RQ3 SDEに基づく運動論的モデルは、複雑で多次元的な天体物理学的流れにおいて、拡散・対流方程式の解析的解を正確に再現できるか？
• RQ4 シミュレーションで観測されたスペクトルカットオフを考慮すると、シンクロtron放射モデルは銀河間ジェットからのX線放射を説明できるのか？
• RQ5 複数の衝撃構造および乱流拡散領域は、全体の粒子加速効率およびスペクトルの硬さにどのように影響を与えるか？

主な発見

• $\eta_T \approx 0.2-0.3$ の乱流レベルが、粒子加速に対して最適であることが判明。これは効率的な拡散輸送に起因し、電子が最高エネルギーに到達可能となる。
• 電子エネルギースペクトルは100 MeVから数TeVにわたり、高速な流れの場合は100 TeVに達するケースも存在するが、横方向の運動量依存的脱出に起因し、凸型の非累乗法則的形状を示す。
• シンクロtron損失により、数百GeV以上でスペクトルカットオフが生じる。これは、Chandra観測における観測されたX線放射をシンクロtronモデルが完全に説明できない可能性を示唆する。
• 混沌とした磁場拡散領域では、有効な衝撃圧縮比が曲がった衝撃の平均圧縮比でよく近似され、特定の条件下では平面衝撃の等価物を用いることが妥当であることが裏付けられる。
• MHD空間分解能によって制限される$D_{\text{min}}$ よりも$D \geq D_{\text{min}}$ の条件下では、SDE手法が拡散・対流方程式の解析的解を正確に再現できた。
• 複数の衝撃構造は、GeVエネルギー領域でより硬いスペクトルをもたらす。これは、衝撃の多重性が高エネルギー粒子生成を強化することを示唆する。