Mg II h&k spectral line properties computed using 3D radiative transfer in an enhanced network region simulated with the MURaM-ChE code

技術的概要：強化ネットワーク領域における 3 次元放射輸送を用いて計算された Mg ii h&k スペクトル線の特性

問題と背景
太陽彩層は、光球とコロナの間の複雑で不均質な界面であり、非局所熱平衡（NLTE）条件下でスペクトル線の形成が起こる。Mg ii h&k などの強い共鳴線は、部分周波数再分配（PRD）効果を示す。3 次元放射磁気流体力学（rMHD）シミュレーションは現実的な大気モデルを提供するが、これらのモデルからのスペクトル線の合成は、往々にして 1.5 次元放射輸送（RT）近似に依存している。1.5 次元 RT では、各垂直カラムが独立した平面平行大気として扱われ、水平放射輸送および水平速度場が無視される。公開された Bifrost シミュレーションのスナップショットを用いた先行研究では、1.5 次元と 3 次元 RT の結果間に不一致が示されたが、より動的で高解像度なモデルにおけるこれらの効果の大きさは、まだ完全に定量化されていない。具体的には、強化ネットワーク（EN）領域における Mg ii h&k 線の形成の詳細と、そのような動的環境における確立された大気診断の堅牢性は、さらなる検証を必要とする。

手法
著者らは、MURaM-ChE コード（MURaM の彩層拡張版）で生成された EN 領域の自己整合的な 3 次元 rMHD シミュレーションを利用した。このモデルは、水平方向 23.4 km、垂直方向 20 km の解像度を持ち、双極性磁気構造を特徴とし、NLTE エネルギー輸送および非平衡水素電離を組み込んでいる。

本研究では、2 つの主要な合成アプローチを採用した：

1. 3 次元放射輸送： Multi3D コードを用いてスペクトルを計算し、完全な 3 次元幾何学において PRD 処理を伴いながら、放射輸送方程式（RTE）と統計的平衡方程式を同時に解いた。
2. 1.5 次元放射輸送： 比較のために、同一の大気スナップショットから、平面平行近似（1.5 次元 RT）を用いて Multi3D 内でスペクトルを計算し、水平放射結合および水平速度を無視した。

合成スペクトルは、Interface Region Imaging Spectrograph（IRIS）の機器条件に一致するように劣化させ、空間的およびスペクトル的な畳み込みを施した。これらは以下のものと比較された：

• 円盤中心付近の双極性構造の IRIS 観測。
• 公開された Bifrost スナップショットからの先行結果。
• 角度の離散化と水平速度の影響を分離した制御実験。

主要な貢献と結果

• スペクトル線プロファイルへの 3 次元 RT の影響： 3 次元 RT で計算された空間平均された Mg ii h&k スペクトルは、1.5 次元 RT による計算に比べて、コア強度とその分布において IRIS 観測と著しくよく一致する。しかし、1.5 次元アプローチは、ピーク強度（k2v および k2r）と中央極小（k3）を著しく過大評価する。MURaM-ChE モデルでは、1.5 次元と 3 次元 RT の不一致は、公開された Bifrost スナップショットよりも顕著である。線コアにおいて 1.5 次元強度は 3 次元結果より最大 72% 高く、一方、以前の Bifrost 研究ではその差は約 12% であった。
• 水平速度の役割： 制御実験により、完全な 3 次元 RT には自然に含まれるが標準的な 1.5 次元 RT には存在しない水平速度が、線プロファイルを著しく広げ、ピーク強度を減少させることが明らかになった。平面平行計算に水平速度を含めることで、結果は完全な 3 次元解に近づき、MURaM-ChE モデルにおける動的な水平速度場が、1.5 次元と 3 次元の結果間の大きな不一致の主要な駆動因子であることを示唆している。
• 形成高度と速度の相関： 本研究は、k2v（青方）および k2r（赤方）のピークが異なる大気高度で形成され、異なる流れ構造を追跡することを確認した：
  • k2v： より低い大気高度（より高い密度、約 6.8 kK）における上昇流で優先的に形成される。
  • k2r： より高い大気高度（より低い密度、約 8 kK）における下降流で優先的に形成される。
  • k3： k3 極小のドップラーシフトは、その形成高度における垂直速度と tight な相関（ピアソン係数 $R_P = 0.94$）を示し、速度診断としての有用性を確認したが、散らばりはより動的な MURaM-ChE 大気のため、Bifrost スナップショットよりも大きい。
• 温度診断： ピーク輝度温度と局所大気温度の相関は、MURaM-ChE モデルでは弱いことが判明した（k2v で $R_P \approx 0.12$、k2r で $0.06$）。約 5.25 kK 以上の強度に対して相関は存在するものの、源関数はしばしば局所ガス温度から脱結合しており、この極めて動的なモデルにおいて Mg ii h&k ピークを直接の温度診断として使用する信頼性は制限される。
• 統計的分布： 3 次元 RT 計算におけるピーク輝度温度およびピーク分離の分布は、1.5 次元の結果よりも IRIS 観測とより密に一致する。しかし、3 次元 RT シミュレーションは、観測に比べてピーク分離を依然として過小評価しており、これはシミュレーションの垂直運動が実太陽よりも低いこと、あるいはシミュレーションボックスのサイズがネットワーク要素と静穏太陽の間の相互作用を制限している可能性を示唆している。

意義と結論
本論文は、特に強化ネットワークのような動的な彩層領域において、Mg ii h&k 線の正確な前方モデリングには 3 次元放射輸送が不可欠であると結論づけている。本研究は、水平放射輸送および水平速度を無視することが、線コア強度の著しい過大評価と線プロファイル形状の誤った表現につながることを実証している。

著者らは、MURaM-ChE モデルが中程度の解像度であるにもかかわらず、公開された Bifrost スナップショットよりもより動的な大気を生成し、より複雑なスペクトルプロファイルおよび 1.5 次元と 3 次元 RT の間のより大きな不一致をもたらしていると強調している。モデルは観測された線幅および一般的な強度分布を成功裡に再現しているが、ピーク分離および強度比における残存する不一致は、太陽彩層の質量負荷およびダイナミクスを完全に捉えるために、シミュレーション解像度および水平ドメインサイズのさらなる改善が必要であることを示唆している。この研究は、速度構造を制約するための Mg ii h&k 線の診断的潜在能力を強化するが、3 次元 RT によって明らかにされた複雑な形成物理を考慮することなく、極めて動的で NLTE 支配的な環境において、それらを直接の温度診断として使用することには警告を発している。