Failed ejection and oscillations of a current-carrying filament balanced by gravity

問題点
太陽フィラメント（プロミネンス）は、太陽コロナの中に浮遊する高密度で低温のプラズマ構造であり、しばしば太陽噴出やコロナ質量放出（CME）に関連している。多くのモデルはフィラメントの安定性と噴出を記述しているが、既存の数値シミュレーションにおける重大な限界は、重力を無視したゼロ・プラズマ・ベータの仮定である。現実的なコロナ条件には、有限のプラズマ圧と重力による成層が含まれており、これらはフィラメントの安定性、質量負荷、および熱的構造を理解する上で極めて重要である。さらに、フィラメントの噴出に関する観測データは、噴出プロセスが停止し、フィラメントが落下する「失敗した」イベントを頻繁に示しており、この現象を完全に捉えるには、重力と現実的な電流構成を考慮できるモデルが必要である。

手法
著者らは、Solov'ev (2010) によって提案された重力バランスモデルを用いて、フィラメントの不安定化後の進化を調査している。本研究では、以下の二つのアプローチを採用している：

1. 解析的研究： 著者らは、Solov'evモデルに基づき、電流密度、プラズマ密度、およびそれらの空間分布に関する解析的な式を導出している。このモデルでは、電流を運ぶ高密度のフィラメントが、フィラメント電流とその下層フォトスフェア内の鏡像電流との間の反発的なローレンツ力によって重力に対して支えられている。
2. 数値シミュレーション： Lare2d 数値コードを用いて、2次元磁気流体力学（MHD）モデルを構築している。このモデルは、重力と電気抵抗を明示的に含んだデカルト座標系において、抵抗性・圧縮性MHD方程式を解くものである。シミュレーション領域（$x \in [-10, 10]$ Mm, $y \in [0, 40]$ Mm）は、$2000 \times 4000$ の格子によって解像されている。初期平衡状態は、磁場 $B_0 = 10^{-3}$ T、質量密度 $\rho_0 \sim 1.3 \times 10^{-9}$ kg m$^{-3}$、温度 $T \sim 13000$ K によって特徴付けられる。フォトスフェアとコロナの間の遷移領域をシミュレートするために、修正された温度プロファイルが使用されている。

システムは、(a) フィラメント内の電流の増加（磁場強度 $B_0$ の増加による）および (b) フィラメントの質量密度の減少という、二つの異なるシナリオを通じて不安定化される。系の進化を追跡することで、噴出、停止、および落下ダイナミクスを観察する。

主な貢献

• 解析的式の導出： 本論文は、全電流およびフィラメントの重量を含むモデルパラメータの関数として、電流密度、プラズマ密度、およびそれらの空間分布に関する解析的な公式を提供している。
• 重力バランスMHDモデル： 重力的効果を無視するゼロ・ベータ・モデルの限界に対処するため、反発的な磁気力と重力のバランスを明示的に組み込んだ2次元MHDモデルを実装している。
• 失敗した噴出の特性評価： 著者らは、不安定化したフィラメントが噴出し、特定の高度で停止し、その後に落下して振動サイクルに入るという「失敗した噴出」現象を、体系的にシミュレートし、特性評価を行っている。
• 定量的依存性： 最大噴出高度は、不安定化させる電流の大きさと、フィラメントのプラズマ希釈の度合いの両方の関数として決定される。

結果

• 振動挙動： 両方の不安定化レジーム（電流増加および密度減少）において、フィラメントは上方へ噴出し、ある高度で停止し、その後落下する。このサイクルは約600秒の周期で繰り返される。
• 噴出速度： 最大噴出速度は、不安定化のメカニズムによって異なる。電流増加シナリオでは約80 km s$^{-1}$、密度減少シナリオでは約40 km s$^{-1}$ である。
• 電流シートの形成： 噴出したフィラメントの下部では、磁気リコネクションが発生する電流シートが形成される。このリコネクションは、上部と下部の電流を結ぶ磁力線を分断する。
• パラメータ依存性：
  • 電流（磁場強度）を増加させると、噴出高度が高くなる。
  • フィラメント密度を減少させると、これも噴出高度が高くなり、密度比（$\rho/\rho_0$）が低いほど、より高い高度に達する。
  • 密度減少シナリオでは、電流増加シナリオと比較して、フィラメントが最大高度に到達する時間が早く、そこに留まる時間も長い（密度減少は約200秒、電流増加は約100秒）が、振動周期は同程度である。
• 非対称性： 電流増加の場合、フィラメント下の密度構造は最終的に非対称性を発達させ、フィラメントの垂直軸からの変位を引き起こす。密度減少の場合、フィラメントは垂直軸に沿って振動しながら安定状態へと向かう。
• 検証： 垂直方向の領域を拡張した比較シミュレーションにより、元の数値ドメインのサイズが磁場構成や到達する最大高度に影響しないことが確認された。

意義
本論文は、重力と有限のプラズマ・ベータを組み込んだモデルが、噴出前の平衡状態と、フィラメントが安定性を失う特定の条件を理解するために不可欠であることを示している。観測される「失敗した噴出」現象とその後の振動挙動を再現することにより、本研究は、フィラメントの形成と不安定化に関する観測データを解釈するための理論的枠組みを提供する。これらの結果は、質量負荷と熱的構造が、不安定性の閾値および噴出プロセスの非線形ダイナミクスに影響を与える重要な要因であることを強調しており、理想的な真空モデルに代わる、より現実的な選択肢を提示している。