Counter-streaming heat-flux closure for electron-only collisionless magnetic reconnection

技術要約：電子のみの衝突のない磁気リコネクションにおけるカウンター・ストリーミング熱流閉鎖

問題提起
乱流プラズマにおける磁気リコネクション（MR）は、リコネクション電流シートが特性イオンスケールよりも小さくなり、イオンがダイナミクスからデカップリングされる「電子のみ」のレジームをますます示すようになっている。運動論的な粒子・衝突・粒子（PIC）シミュレーションは、電子のみのMRが、超イオン・アルヴェン速度のジェットなどの特異な現象を示すことを明らかにしているが、従来の流体モデルはこれらのダイナミクスを捉えることができない。ハメット・パーキンス（Hammett-Perkins）閉鎖（準マックスウェル分布に依存）や Le et al. の圧力閉鎖（高$\beta$レジーム向けに導出）といった既存の流体閉鎖は、電子のみのMRで観察される強く非マックスウェル的でマルチモーダルな電子分布を再現することが根本的に不可能である。具体的には、既存のモデルは、異なる流体速度を持つコヒーレントな電子サブ集団が同じ空間的位置を占める「電子のカウンター・ストリーミング（逆ストリーミング）」を説明できない。この限界により、地球の磁気尾部のような環境の理解に不可欠な、電子のみのMRの効率的な流体モデリングが妨げられている。

手法
著者らは、電子のみの衝突のないMRにおける電子熱流の物理的起源を調査するために、完全運動論的PICシミュレーション（OSIRISコードを使用）を用いている。磁気尾部を代表するパラメータ（$\beta_{Reconn.} = 0.1$、 $B_{Guide}/B_{Reconn.} = 0.1$）を用いて、2D-3Vの自由力学的なダブル・ハリス電流シートを初期化している。電流シートの中面の上側または下側のどちらかに初期位置に基づき電子をタグ付けすることで、著者らは全電子分布関数 $f$ を2つのサブ集団（$f_1$ と $f_2$）に分解している。

この分解を通じて、著者らは全熱圧および熱流テンソルの解析的な表現を導出している。彼らは、全熱流がサブ集団内の非断熱的挙動からではなく、2つのサブ集団間の相対的なドリフトと温度差から生じる「カウンター・ストリーミング」項（$q_{CS}$）によって支配されていることを示している。

この知見に基づき、著者らは新しい2電子流体（2-e-Fluid）閉鎖を開発した。このモデルでは：

1. 電子集団は、電流シートの反対側に初期化された2つの異なる流体に分割される。
2. 各サブ集団は近似的に断熱的（$q_{1,ijk} \approx q_{2,ijk} \approx 0$）として扱われる。
3. 全熱流は、2つの流体の相対運動によって閉鎖される（$q_{ijk} \approx q_{CS,ijk}$）。これにより、3次モーメント方程式の必要性が排除される。

著者らは、この閉鎖をGkeyll流体コードに実装し、10個のモーメント（密度、速度、および圧力テンソル成分）を持つ2つの電子流体を進化させている。比較のために、同じ条件下で従来の単一電子流体（1-e-Fluid）シミュレーションも実行している。

主な結果

• 熱流の起源: PIC解析により、電子のみのMRにおける電子熱流は、電流シートの反対側から来たサブ集団間のカウンター・ストリーミングにほぼ完全に由来することが確認された。全分布は高度に非マックスウェル的であるが、各サブ集団はリコネクション領域を通じて近似的に断熱的な性質を保持している。
• 流体モデルの性能: 2-e-Fluidモデルは、以下のものを含む、PICシミュレーションで観察される主要な運動論的特徴を正常に再現している：
  • リコネクション電流密度 ($J_z$): 密度減少による非物理的な密度欠乏（これにより $J_z$ が抑制される）に苦しむ1-e-Fluidモデルとは異なり、2-e-Fluidモデルはサブ集団が中面を横切ることを許容することで、電流シートを満たした状態に保つ。
  • 熱圧 ($p_{xx}$): 1-e-Fluidモデルは、熱流チャネルの欠如により、アウトフローにおいて偽の圧力蓄積を生じさせる。2-e-Fluidモデルは、カウンター・ストリーミング熱流を捉え、熱エネルギーをアウトフローから輸送することで、圧力の蓄積を防ぐ。
  • 熱流構造: 2-e-Fluidモデルは、全熱流の大部分を構成するカウンター・ストリーミング熱流（$q_{CS,xxx}$）の空間構造と符号を再現する。
• 有効領域: $\beta_{Reconn.}$ および $B_{Guide}/B_{Reconn.}$ に関するパラメータスキャンにより、この閉鎖はこれらの値が低い場合（地球の磁気尾部に相当する領域）に最も正確であることが明らかになった。精度は、初期電子熱速度と初期バルク流速度の比（$v_{th}/v_{fl}$）によって制御される。この比が増大したり、サブ集団が磁気島内で循環し始めたりして単純なカウンター・ストリーミングの仮定が崩れたりすると、モデルの精度は低下する。

意義と主張
本論文は、カウンター・ストリーミングを電子熱流の物理的起源として特定することにより、電子のみの衝突のないMRのモデリングにおける根本的な限界を解決したと主張している。複雑な非マックスウェル分布が2つの断熱的なサブ集団に分解できることを示すことで、著者らは、これまで運動論的シミュレーションでのみアクセス可能であった物理を捉える、計算効率の高い流体閉鎖を確立した。

著者らは、この2-e-Fluidアプローチが、従来の単一流体モデルでは「根本的に再現できない」電子のみの衝突のないMRのモデリングを可能にすると主張している。彼らは、この手法が、熱圧に影響を与える熱流の動的な重要性を捉えることを強調している。本論文は、3D電子のみのMRおよび乱流リコネクションへの潜在的な拡張を示唆しているが、現在の閉鎖は低$\beta$・低ガイド磁場レジームにおいて最も堅牢であり、複雑な島内ダイナミクスを解明するには追加の流体が必要になる可能性があることを控えめに述べている。本研究は、不可欠な運動論的物理であるカウンター・ストリーミングを犠牲にすることなく、効率的な電子のみのリコネクションの流体モデリングへの道筋を提供するものである。