The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

技術要約：三次元相対論的リコネクションにおける粒子加速への逆コンプトン損失の影響

問題提起
相対論的磁気リコネクションは、コンパクト天体のマグネットスフィアや相対論的ジェットなどの天体物理学的環境において、磁気エネルギーを散逸させ粒子を加速させる主要なメカニズムである。近年の第一原理粒子・イン・セル（PIC）シミュレーションは、3Dリコネクションにおける粒子加速のモデル（Zhang et al. 2023）を確立しており、粒子がアップストリーム領域で活発に加速される「自由フェーズ（free phase）」と、粒子が顕著なエネルギー獲得なしにダウンストリームのフラックスロープ内に閉じ込められる「トラップフェーズ（trapped phase）」を区別している。このモデルは、放射冷却がない場合、自由粒子に対しては普遍的なべき乗則スペクトル（$dN/d\gamma \propto \gamma^{-1}$）を、トラップ粒子に対してはより急峻なスペクトル（$dN/d\gamma \propto \gamma^{-2}$）を予測する。

しかし、放射冷却、特に逆コンプトン（IC）損失が存在する場合のこの加速メカニズムの挙動については、まだ十分に解明されていない。冷却損失がリコネクション率、層の構造、あるいは自由フェーズとトラップフェーズの根本的な物理を変化させるのかどうかは不明であり、特に放射反作用制限ローレンツ因子（$\gamma_{\text{rad}}$）が磁化率（$\sigma$）を上回る「弱冷却」領域においては特に不明である。

手法
著者らは、IC冷却損失を伴う相対論的リコネクションを調査するために、TRISTAN-MPコードを用いて一連の3D PICシミュレーションを実施した。シミュレーションでは、ガイドフィールド（$B_g = 0.1 B_0$）を持つハリス電流シート構成を用い、電流方向に周期境界条件を設定した。

主なシミュレーションパラメータは以下の通りである：

• プラズマ組成： 電子・陽電子ペアとし、放射冷却は電子にのみ適用した。陽電子は冷却させないことで、電子・イオンプラズマにおける重いイオンの挙動を模倣し、層の力学的バランスへの変化を最小限に抑え、Zhang et al. (2023) の非冷却モデルとの直接比較を可能にした。
• 磁化率： 固定された磁化率 $\sigma = 10$。
• 冷却領域： 本研究は弱冷却領域（$\gamma_{\text{rad}} > \sigma$）に焦点を当てた。放射反作用制限ローレンツ因子について、$\gamma_{\text{rad}} = 30, 43,$ および $60$の3つの特定の値をテストした。これらは、それぞれ約$3.5, 7.2,$および$14.2$の冷却ブレイク（$\gamma_{\text{cool}}$）に対応する。
• 解析： 著者らは、混合係数の基準を用いて「自由（アップストリーム）」粒子と「トラップ（ダウンストリーム）」粒子を区別するために、大規模な粒子サンプル（1実行あたり50万個の電子）を追跡した。加速率、脱出時間、およびエネルギースペクトルを分析し、それらを冷却を含むように拡張された Zhang et al. (2023) のモデルから導出された解析的予測と比較した。

主な貢献と結果

1. リコネクション動力学の堅牢性： シミュレーションは、電子の冷却損失がリコネクション率や電流層のマクロな構造を著しく変化させないことを示している。定常状態のリコネクション率は $\eta_{\text{rec}} \approx 0.05$ に落ち着き、これは非冷却の場合と一致する。より強い冷却は層をわずかに薄くし、プラズモイドを小さくするものの、テストされた冷却強度全般にわたって基本的な層の物理は変化しない。

2. 自由フェーズ物理の検証： 本研究は、加速の「自由フェーズ」が冷却損失に対して鈍感であることを確認した。自由電子は一貫して最適な速度（$\dot{\gamma}_{\text{acc}} \sim \eta_{\text{rec}}\omega_c$）で加速され、非冷却の場合と同一の $dN_{\text{free}}/d\gamma \propto \gamma^{-1}$ のべき乗則スペクトルを示す。このスペクトルのカットオフは、システムサイズではなく $\gamma_{\text{rad}}$ によって決定される。

3. トラップフェーズ・スペクトルの変調： 重要な知見は、冷却が存在する場合のトラップ粒子のスペクトルの変調である。

   • 非冷却の場合、トラップスペクトルは $dN/d\gamma \propto \gamma^{-2}$ となる。
   • 弱冷却領域では、冷却ブレイクを超えるエネルギー（$\gamma_{\text{cool}} < \gamma < \gamma_{\text{rad}}$）において、トラップされた電子のスペクトルは $dN/d\gamma \propto \gamma^{-3}$ へと急峻になる。
   • この急峻化は、Zhang et al. (2023) のモデルの根本的な仮定、すなわち「トラップフェーズ中には有意なエネルギー獲得が起こらない」という仮定を裏付けている。この急峻化は、トラップされた粒子が移流時間と同等またはそれよりも短いタイムスケールで冷却されるため、非冷却の場合と比較して高エネルギー領域の粒子数が減少することに起因する。

4. 粒子の運動学： シミュレーションにより、自由粒子はリコネクション電場にほぼ沿って移動し（$\beta_z \approx -1$）、迅速な加速を促進する一方で、トラップ粒子は系統的なエネルギー獲得の欠如と一致する、ゼロを中心とした対称的な $\beta_z$ 分布を示すことが検証された。

意義と天体物理学的含意
本論文は、弱放射冷却が存在する状況下での3D相対論的リコネクションに関する Zhang et al. (2023) の解析モデルを検証している。主要な含意は、リコネクション駆動の天体源（ガンマ線バーストやブレーザー・ジェットなど）の放射モデルにおいて、**二領域構造（two-zone structure）**を採用しなければならないということである：

• 領域1（自由）： スペクトル $dN/d\gamma \propto \gamma^{-1}$ を持つ、急速に加速される粒子。
• 領域2（トラップ）： 冷却ブレイク以上でスペクトルが $\gamma^{-2}$ から $\gamma^{-3}$ へと急峻化する、受動的に冷却される粒子。

著者らは、放出に関与する粒子の大部分はトラップされているものの、最大ローレンツ因子（$\gamma_{\text{rad}}$）付近では自由粒子がより豊富になることを指摘している。自由粒子はリコネクション電場に沿って優先的にビームされるため、その放射は高速な時間変動などのユニークな観測的シグネチャーをもたらす可能性がある。

限界と注意点
著者らは、自らの結果が $\gamma_{\text{rad}} > \sigma$ である弱冷却領域に適用されることを明示している。このモデルは、$\gamma_{\text{rad}} < \sigma$ となる強冷却シナリオには適用できない。さらに、シミュレーションにおける $\gamma_{\text{cool}}$ と $\gamma_{\text{rad}}$ の間のダイナミックレンジは、現実の天体物理学的システムよりも小さい。今後の課題として、これらの知見を現実的な質量比を持つ電子・陽子プラズマや、有意なガイドフィールドを持つリコネクションへと拡張することが示唆されている。