ExTraSS: a Domain Decomposed 3D NLTE Radiative Transfer spectral synthesis code for nebular phase transients

3D NLTE 放射伝達スペクトル合成コード ExTraSS の技術的概要

問題提起
超新星（SNe）の星雲相のモデル化は、大きな計算上の課題を呈する。この相では、噴出物の密度が十分に低下し、構造が広範に光学的に薄くなるため、体積全体が放射に寄与する。しかし、正確なモデル化には、複雑な微物理を考慮するために非局所熱平衡（NLTE）物理が必要であり、特に水素に富む核崩壊型超新星の場合、放射伝達（RT）効果を無視することはできない。さらに、超新星の非対称な形態は、3 次元入力モデルを必要とする。

3 次元 NLTE 放射伝達における主要なボトルネックは、光励起率の生成と保存である。$\sim 10^5$ 個のセルを持つグリッドにおいて、数百万の遷移に対するこれらの率を保存するには、1 TB を超えるデータが必要となり、現代のスーパーコンピュータにおける単一計算ノードのメモリ容量を超えてしまう。以前の 3 次元 NLTE の取り組みは、1 次元に限定されていたか、光学的に薄い近似を用いていたか、あるいは全球的な放射場効果を無視する「その場（on-the-spot）」処理を採用していた。放射場と原子準位集団の間の全球的な結合を伴う完全な 3 次元 NLTE 放射伝達を実行できるコードが必要とされていた。

手法
著者らは、完全な 3 次元 NLTE 放射伝達を実行するようにアップグレードされたコード ExTraSS（EXplosive TRAnsient Spectral Simulator）を提示する。このコードは「スナップショット」モードで動作し、光学的深さが低い後期相に適した、時間とは独立した光子輸送を追跡する。

• グリッドと幾何学: ExTraSS は球座標グリッドを採用しており、内側の低速で移動する噴出物では高分解能を、外側の半径ではより大きなセルを可能にする。放射場の角度変化は、$N_{\theta, VA} \times N_{\phi, VA}$ 個のビューア（デフォルト 20 $\times$ 20）を用いて離散化される。
• $\gamma$ 線輸送: コードは、コンプトン散乱を用いたモンテカルロ法により、放射性崩壊（主に $^{56}$Ni と $^{56}$Co）を追跡する。陽電子は局所的に閉じ込められると仮定される。$\gamma$ 線からのエネルギー付与は、Spencer-Fano 形式を通じて非熱的加熱、電離、励起を駆動する。
• NLTE ソルバ: ソルバは、局所温度と準位集団を反復的に決定する。熱的および非熱的物理を結合し、$\gamma$ 線および光電離からの加熱と、束縛 - 束縛、自由 - 束縛、および 2 光子放射からの冷却とのバランスを取る。ソルバは、準位集団を収束させるために、適応的ダンピングを備えたニュートン - ラフソン法を使用する。
• 放射伝達（RT）: コードは、全球的な放射場を計算するために完全なレイトレーシング手法を実装している。レイは、セルの中心からすべての視角に向けて放射される。コードは、光電離（励起状態には水素様近似を、基底状態には Verner らのデータを使用）および光励起（線相互作用にはソボレフ近似を使用）に対する光学的深さを計算する。
• ドメイン分解: メモリ制限を克服するために、ExTraSS はドメイン分解アルゴリズムを採用する。計算ドメインは $\phi$ 基準のスライス（方位角セクター）に分割され、各ノードが特定のセクターを処理する。
  • 通信アーキテクチャ: 各ノードは、ノード間通信を処理する 1 つの「マネージャ」コアを指定し、残りの「ワーカー」コアはレイトレーシングと局所計算を実行する。
  • 最適化: システムは、ワーカーのアイドル時間を最小化するために非ブロッキング MPI 通信を使用する。データは個別にではなく、束（レイ）として転送される。マネージャコアはドメイン間のレイ転送を調整し、ワーカーがデータ待ちではなく計算に最大限の時間を費やすことを保証する。

主な貢献

1. 超新星における最初の 3 次元 NLTE 放射伝達: 著者らの知る限り、この研究は超新星放射伝達コードに対してドメイン分解を初めて実装したものであり、メモリ制約のために以前は不可能であった完全な 3 次元 NLTE 計算を可能にする。
2. メモリ最適化: 本論文は、共有メモリ、「ジャギッド配列」（既存の種/準位のみに対してメモリを割り当てる）、およびドメイン分解スキームの実装の詳細を記述しており、これによりメモリ使用量は使用されるノード数に反比例してスケーリングされる。
3. アルゴリズムの洗練: 著者らは、計算コストの高い Spencer-Fano サブルーチンの呼び出し頻度を最適化し（1、4 回目、およびその後は 10 回ごとに呼び出す 4+10 スキームを使用）、ニュートン - ラフソン準位集団ソルバのダンピングスキームを改良して、収束速度と安定性のバランスを取った。

結果と検証
著者らは、いくつかのベンチマークを通じて ExTraSS を検証した。

• 収束性と効率: Dardel スーパーコンピュータ（128 コアノードを使用）でのテストにより、負荷分散が最適化された場合、放射伝達モジュールの計算効率は $\sim 0.8$ に達することが示された。全球的な収束（全波長光度と $\gamma$ 線付与の比率）は $\sim 17$ 回反復で達成され、通常 0.96–1.0 の比率で安定する。
• SUMO（1 次元）との比較: ExTraSS は、「1 次元から 3 次元へ変換された」モデル（J18）を用いて、1 次元コード SUMO と比較された。
  • 物理的条件: 温度と電子分率はよく一致したが、ExTraSS は中性分率が低いため、内側の噴出物でわずかに高い温度を示した。
  • 率: 光電離率と光励起率は一般的に互換性があったが、ExTraSS は最も内側の領域でわずかに高い率を示した。
  • スペクトル: 合成スペクトルは良好な一致を示したが、Fe II 放射と H$\alpha$ 強度にわずかな不一致があり、これは含まれる原子準位数の違い（ExTraSS はイオンあたり最大 100 準位、SUMO は数百準位）に起因すると考えられた。
• マルチコードベンチマーク（Blondin ら 2022）との比較: ExTraSS は、'toy06'モデルに対して ARTIS-nebular、CMFGEN、および SUMO とともにテストされた。
  • 物理的条件: ExTraSS は、内側の噴出物でいくつかのコードよりもわずかに低い温度と、コバルトの電離分率を計算したが、結果は全体的に集合体と一致していた。
  • スペクトル: スペクトル出力は、特に 5100–6400 Å および >8800 Å の範囲で、SUMO および他のコードとよく一致した。4600–5000 Å 周辺の Fe 放射ピークには不一致が認められ、ExTraSS はより低い光度を予測した。これは、準位集団の切り捨てまたは電離の違いに起因する可能性がある。

重要性
本論文は、ExTraSS が現在、完全に機能するドメイン分解された 3 次元 NLTE 放射伝達スペクトル合成コードであると主張している。その重要性は、3 次元 NLTE 物理の莫大なメモリ要件を処理する能力にあり、星雲相における非対称な超新星爆発の高精度なモデル化への道を開く。現在、核崩壊型超新星モデル（特に Type II）に適用されているが、著者らは、このコードが将来のアップデートで追加の物理成分（時間依存性、塵、または分子など）が組み込まれれば、多様な入力モデルおよび爆発的遷移現象に適用可能になるように設計されていると述べている。この研究は、詳細なスペクトルシグネチャを通じて爆発特性、噴出物質量、および組成をより厳密に制約するための基盤を確立する。