Radiation magnetohydrodynamics modeling of an impulsively driven chromospheric jet in the solar atmosphere

技術的概要：衝動的駆動による彩層ジェットの数値モデル化における放射磁気流体力学

問題提起
太陽ジェット、特にスピキュールおよびマクロスピキュールは、複数の波長で観測される動的なプラズマ噴出現象である。磁気リコネクションや衝動的な擾乱（圧力パルスや速度パルスなど）が引き金となるメカニズムとして広く受け入れられているが、これらのジェット進化における放射過程の役割は依然として調査の対象となっている。過去の数値研究では、理想的または非理想的な磁気流体力学（MHD）を用いて、しばしば経験的な放射冷却関数と結合させた衝動的駆動ジェットをモデル化してきた。しかし、周波数平均および角度平均の放射輸送方程式の系と完全に結合された非理想的 MHD 方程式を用いてこれらのジェットをシミュレーションする研究には欠落がある。具体的には、放射エネルギー密度と放射フラックス成分の両方を独立に時間発展させ、光学薄な太陽コロナにおけるジェット形態および散逸への具体的な影響を決定するシミュレーションが存在しない。

手法
本研究では、PLUTO コード（v4.4-patch3）を用いた 2 次元数値シミュレーションにより、衝動的に駆動される彩層ジェットをモデル化した。物理的枠組みは以下の通りである：

• 支配方程式： 系は、灰色近似の下での 2 モーメント放射流体力学スキームと結合された非理想的 MHD 方程式を解く。これには、プラズマ密度、運動量、全エネルギー、磁場の進化に加え、放射エネルギー密度（$E_r$）および放射フラックス（$F_r$）の進化が含まれる。
• 放射輸送： 放射モーメント方程式を閉じるために M1 閉じ込みが用いられ、拡散と自由飛行の両極限を記述可能である。放射 - 物質相互作用項は、吸収と散乱の光深を考慮し、$\beta = v/c$ の 1 次までのローレンツ変換を通じて導出される。
• 太陽大気モデル： 初期状態は、コロナまで拡張された半経験的 C7 モデルに基づく静水圧平衡である。プラズマ組成は、平均分子量 $\mu = 0.6$（完全電離した水素とヘリウム）を仮定する。
• 物理パラメータ： シミュレーションでは、光学薄な太陽コロナに特徴的な、低いプランク平均吸収光深（$\kappa = 1 \times 10^{-6}$ cm$^2$ g$^{-1}$）およびロースランド平均散乱光深（$\sigma = 0.34$ cm$^2$ g$^{-1}$）を用いる。熱伝導はフラックス制限アプローチでモデル化され、経験的なコロナ加熱関数が背景大気を維持する。
• 摂動： ジェット形成は、$y = 1.75$ Mm における垂直速度成分（$v_y$）の局所的なガウスパルス（振幅 100 km s$^{-1}$）によって開始される。
• 数値スキーム： 方程式は 2 次 TVD ランゲ・クッタ時間積分器を用いて解かれる。MHD 方程式には HLLD リーマンソルバが、放射輸送には HLLC ソルバが用いられる。剛性の高い放射 - 物質相互作用項は、陰陽（IMEX）ランゲ・クッタスキームで処理され、計算の実行可能性を確保するため、$\hat{c} = 0.1c$ の減速光近似（RSLA）が採用される。

主要な結果
本シミュレーションは、マクロスピキュールに似たコリメートされたプラズマ構造の形成と進化を成功裡に再現した：

• ジェット形態と力学： 衝動的パルスは衝撃波を生成し、プラズマを上方へ押し上げ、約 8 Mm の高さに達するジェットを形成する。ジェットは、衝撃波前面で最大垂直速度 $\sim$60 km s$^{-1}$、ジェット本体内部で $\sim$30 km s$^{-1}$ を示す。構造は後期段階で二重スレッド形態を示し、これは太陽スピキュールに共通する特徴である。
• 熱的進化： ジェット頭部は $\sim$2 $\times$ 10$^5$ K に達し、ジェット本体は約 2 $\times$ 10$^5$ K で安定化する。衝撃領域は冷却前に一時的に $\sim$9 $\times$ 10$^5$ K に達する。
• 放射効果： 自己無撞着な放射輸送の導入により、光学薄なコロナにおいて放射が効率的な冷却機構として機能することが明らかになった。放射エネルギー密度は線形に増加し、放射フラックスは減少するが、これは光学薄な放射損失と一致する。
• 散逸メカニズム： 決定的なことに、放射輸送項はジェットの散逸に寄与する。プラズマと放射場の結合は、コロナ内で負の流速（下方運動）を誘起し、ジェットを拡散させて下部大気へ散逸させる。放射変数はジェットの初期形状や形態を決定するものではないが、その進化と寿命を著しく支配する。
• 背景の安定性： ジェット摂動がない場合、背景太陽大気は熱源項（熱伝導、加熱、放射）に対して全体的に安定しており、わずかな温度変動のみが観測される。

意義と主張
本論文は、放射エネルギー密度とフラックスを経験的冷却関数に依存せず、2 モーメント輸送方程式を通じて自己無撞着に時間発展させる、衝動的駆動彩層ジェットの数値シミュレーションとして初めてのものであると主張する。主な意義は以下の点にある：

1. 放射冷却の重要性： 光深が低い光学薄領域であっても、放射輸送の明示的な扱いが、マクロスピキュール様現象の高さ、持続時間、エネルギー含有量を制限する支配的な冷却機構として機能する。
2. 散逸力学： 放射項は、ジェットの急速な減衰と散逸を引き起こす下方流を誘起する。この過程は経験的損失関数を用いた過去の研究と整合的であるが、現在は放射輸送の枠組み内で第一原理から導出されたものである。
3. 方法論的進展： 本研究は、非理想的 MHD と結合された M1 閉じ込みが太陽ジェットモデル化に有効であることを検証した。放射は初期ジェット形態を変化させないが、正しい熱力学的進化と散逸時間スケールを捉えるためには不可欠であることを示した。

著者らは、シミュレートされたジェットが観測されたマクロスピキュールの主要な動的特徴（衝撃波形成や二重スレッド構造など）を再現している一方で、その特定のパラメータ（高さ、速度、温度）は観測範囲の下限に位置しており、モデル化された事象は比較的低エネルギーの発生であることを示唆していると結論づけている。この研究は、このような太陽大気現象の散逸を正確にモデル化するためには、放射輸送の自己無撞着な扱いが必要であると強調している。