Variational Formulation of a Hybrid Kinetic and Gyrokinetic Model for Astrophysical and Laboratory Plasmas

技術要約：縮退した運動論的プラズマの変分定式化

問題提起
太陽風の乱流における運動論的スケールでのエネルギー散逸は、宇宙物理学において依然として未解決の重要な問題である。磁気流体力学（MHD）は慣性範囲を効果的に記述するが、イオンおよび電子スケールにおける熱エネルギーへの遷移には、波と粒子の相互作用や不可逆的な熱力学的プロセスを含む複雑な運動論的効果が関与している。標準的な完全運動論的モデルは、多くの天体物理学的応用において計算コストが極めて高く、一方で従来のジャイロ運動論的モデルは、しばしば低周波数領域（$\omega \ll \Omega_i$）に制限されており、高周波現象や宇宙乱流に関連する特定のパラメータ選択を十分に捉えられない場合がある。したがって、流体/ハイブリッド近似と完全運動論的扱いの間のギャップを埋める、特に電子の運動論的効果を保持しつつ、アクセス可能な周波数範囲を拡張できる縮退モデルが必要とされている。

手法
著者らは、完全運動論的なイオンとジャイロ運動論的な電子を、一貫した数学的形式を用いて結合させたラグランジュ・モデルを提示している。導出は、ハミルトニアン形式のガイド・センター運動に適用された微分幾何学、変分法、および高次のリエ・変換摂動法に基づいている。

1. 座標変換： モデルは、電子のジャイロ位相依存性を排除するために、ジャイロ運動論的な座標縮退を用いている。これには、背景磁場が摂動に対して支配的である場合に磁気モーメントの保存を保証するように、電位および磁気ポテンシャルに関する1次の摂動展開が含まれる。この変換は、変動を除去する新しいハミルトニアンを生成するためにリエ変換を利用しており、最終的にジャイロセンター座標系が得られる。
2. 変分原理： 前提となるハイブリッド・アプローチがプルバック／プッシュフォワード法を使用しているのに対し、本研究では変分法を用いて場の方程式を導出している。全作用 $S$ は、ジャイロ運動論的電子の粒子ラグランジアン ($L^p_{gy}$) と完全運動論的イオンの粒子ラグランジアン ($L^p_{fk,i}$)、および電磁場ラグランジアン ($L^f$) から構成される。
3. 場の方程式： 作用をスカラーポテンシャル $\phi_1$ およびベクトルポテンシャル $A_1$ に関して変分することにより、著者らはポアソン方程式と平行方向のアンプエールの法則を導出している。この導出の重要な特徴は、不均質な多様体の整合的な取り扱いであり、場の方程式は、ジャイロ平均化および座標変換を介して、物理座標空間（イオン）と縮退したジャイロ運動論的座標空間（電子）を結合させている。
4. 数値実装： 得られた方程式は「Super Simple Vlasov (ssV)」コードに実装された。このソルバーは、Vlasov方程式を解くために、5次単調性・正値保持（SLMP5）スキームを用いたセミ・ラグランジュ法を使用している。時間発展にはStrang分割を採用し、周期境界条件を持つスラブ幾何学で動作する。

主要な貢献

• ハイブリッド運動論的・ジャイロ運動論的フレームワーク： 本論文は、イオンを完全運動論的に、電子をジャイロ運動論的に扱う変分ハイブリッドモデルを導入している。このアプローチは、完全運動論的シミュレーションと比較して計算コストを抑えつつ、選択された電子の運動論的効果を保持することを目的としている。
• 変分による導出： 場の方程式は単一の作用原理から一貫して導出されており、これは非変分的な導出と比較して、より優れた保存特性（エネルギー、運動量、シンプレクティック容量）を保証すると著者らは主張している。
• 拡張された周波数範囲： ポアソン方程式に「マクスウェル型ラプラシアン」項（電場の発散と電荷密度の不均衡から導出）を組み込むことで、モデルは標準的なジャイロ運動論的オーダーを超えて拡張されている。これにより、標準的なジャイロ運動論的扱いでは到達できない、ローハイブリッド波やランダウ周波数に接近するモードなどの高周波枝を捉えることが可能になる。
• 一貫した結合： 本モデルは、電磁場を通じて、異なる位相空間次元（イオンの6次元 vs 電子の5次元）を持つ種を厳密に結合させる手法を提供する。

結果
著者らは、線形分散解析と非線形シミュレーションを通じてモデルを検証した。

1. 線形波解析： モデルは、完全運動論的ソルバーであるDSHARKおよびハイブリッド流体・運動論的ソルバーであるHYDROSと比較検証された。
   • イオン音波 (IAW): モデルはIAWの線形分散関係とランダウ減衰率を正常に再現しており、解析解およびDSHARKの結果と一致している。マクスウェル型ラプラシアンの導入により、高波数（$k$）におけるモデルの挙動が改善された。
   • イオンベルンスタイン波 (IBW): モデルはIBW枝および高次IBWモードを捉えており、DSHARKとの一致を示している。これは、モデルが完全運動論的な線形特徴を再現できる能力を持つことを裏付けている。
   • 高周波波 ($\omega \gg \Omega_{ci}$): モデルは、ローハイブリッド波を含む高周波現象を正常に記述している。伝播角のスキャンにより、縮退ドリフト運動論的極限で見られる偽の周波数爆発を回避し、準平行伝播において電子プラズマ周波数に接近することを示した。
2. 非線形ランダウ減衰： ssVソルバーは、非線形ランダウ減衰のベンチマークに対してテストされた。
   • シミュレーションは、線形減衰フェーズ（$\gamma_1 \approx -0.289$）とその後の非線形成長フェーズ（$\gamma_2 \approx 0.078$）を正確に再現した。
   • システムは優れた全エネルギー保存を示し、長時間のシミュレーションにおいても永続的なドリフトは見られなかった。
   • 位相空間解析により、トラッピングおよびフィラメント構造の形成が示され、より細かい速度分解能を用いることで、より鋭い勾配と数値拡散の低減が得られた。

意義と主張
本論文は、本研究が、完全運動論的かつジャイロ運動論的である電磁非線形ハイブリッドモデルを構築するための第一歩であると主張している。主な意義は、以下の特性を持つ縮退モデルを提供することにある：

• 太陽風におけるエネルギー散逸の理解に必要な、本質的な電子の運動論的効果を保持する。
• 標準的なジャイロ運動論的扱いを超えてアクセス可能な周波数範囲を拡張し、イオンサイクロトロン範囲の散逸に関連する物理を捉える可能性がある。
• 保存則を自然に保存する、幾何学的に一貫したフレームワーク（接触幾何学およびシンプレクティック幾何学を介して）を提供する。

著者らは、このモデルは完全な高周波電磁フレームワークではないものの、特定の極限においてランダウ周波数に接近する波の振る舞いを捉えるのに十分な物理を保持していると控えめに述べている。彼らは、次のステップとして、運動論的アルヴェン波（KAW）とIBWの結合を研究するための完全な電磁分散関係ソルバーの実装を進める予定であり、それが太陽風の乱流散逸メカニズムをさらに解明することになると強調している。現在の研究は、この縮退アプローチの理論的基礎と線形妥当性を確立するものである。