Signatures of $^{56}$Ni Mixing and Neutron-rich Ejecta in Supernovae

技術要約：超新星における$^{56}\text{Ni}$混合と中性子過剰噴出物のシグネチャー

問題提起
超新星（SN）の光度曲線は、通常、一つのゾーンを用いた放射性崩壊モデル（例：Arnett 1980, 1982）を用いて解釈され、前駆体の性質（噴出物質量 $M_{\text{ej}}$、ニッケル質量 $M_{\text{Ni}}$、爆発エネルギー $E_k$ など）が推定される。これらのモデルは、単一の有効な拡散タイムスケールを仮定しており、放射性物質の空間分布を無視している。しかし、流体シミュレーションによれば、^{5 \text{Ni}}$ は噴出物の深部で合成された後、レイリー・テイラー不安定性、対流、および乱流によって外側へと混合される。この空間分布は、一ゾーンモデルが欠落させている極めて重要な物理変数である。その結果、高度に混合された噴出物による光曲線は、噴出物質量が少ない、あるいは追加の中央エンジン（マグネターや降着によるフォールバックなど）によるエネルギー供給が必要であると誤認される可能性がある。さらに、コラプサーにおけるrプロセス濃縮の観測可能性は、ニッケルによる加熱背景に対する中性子過剰な噴出物の配置と幾何学的形状に敏感に依存しており、これは球対称モデルではしばしば簡略化されすぎる要因である。

手法
著者らは、RedbackフレームワークのJAXベースのプラグインとして実装された、半解析的なマルチシェルモデルである snmix を導入する。一ゾーンモデルとは異なり、snmix は、速度が幾何学的に配置された $N$ 個のシェルへと、ホモロジー的に膨張する噴出物を離散化する。

• 放射性構造： 本モデルは、$^{56}\text{Ni}$ の任意の半径方向分布（トップハット型またはコア＋テール型のプロファイルにおける混合分率 $f_{\text{mix}}$ によってパラメータ化される）およびrプロセス物質を許容する。
• 輸送物理： 各シェルは独立して進化し、熱エネルギー、断熱損失、および放射拡散を追跡する。決定的なことに、モデルは、各シェルに外部に存在する累積カラム密度に基づいて、拡散およびガンマ線沈着を計算する（グローバルな平均値ではなく）。
• 不透明度とSED： モデルは、グレイな不透明度と温度依存の不透明度則を組み込んでいる。rプロセスの計算については、ランタノイド依存の不透明度と、核ネットワーク計算によって較正された解析的な加熱率を含んでいる。
• 幾何学： 主に球対称であるが、コラプサー用の二成分「ディスク・ウィンド」幾何学（極方向のSN成分と赤道方向の中性子過剰なウィンドを組み合わせ、視線角によって重み付けしたもの）へと拡張可能である。
• 推論： 著者らは、このフレームワークを用いて合成ボロメトリック光曲線を生成し、一ゾーンによるフィッティングとマルチシェルの実態を比較するために、ベイズ推論（NUTSサンプリングを使用）を行い、パラメータのバイアスを分析している。

主な貢献と結果

1. $^{56}\text{Ni}$ 混合が光曲線に与える影響：

   • $^{56}\text{Ni}$ の外側への混合は、上層の拡散カラムおよびガンマ線捕捉カラムを減少させる。これにより、固定された $M_{\text{Ni}}$, $M_{\text{ej}}$, $E_k$ に対して、より速く明るい上昇と、より急峻な放射性テイルをもたらす。
   • 速く明るい上昇は、低噴出物質量やエンジンの存在を示す排他的な証拠ではない。高度に混合された放射性モデルは、追加の動力源なしにこれらの特徴を再現できる。

2. 一ゾーン推論におけるバイアス：

   • 混合された光曲線を一ゾーンモデルでフィッティングすると、視覚的には良好なフィットが得られるが、物理パラメータにはバイアスが生じる。
   • 不透明度が固定されている場合、モデルは欠落している半径方向の構造を、より低い $M_{\text{ej}}$ とより高いニッケル含有率（$f_{\text{Ni}}$）を推論することで吸収してしまう（$M_{\text{Ni}}$ は比較的安定しており、固定不透明度の条件下では $\sim3\%$ 以内で回収される）。
   • 不透明度が変数として許容される場合、バイアスは大幅に低い有効不透明度の推論へとシフトする（完全に混合されたケースでは $\kappa_{\text{fit}}/\kappa_{\text{input}} \simeq 0.24$）。
   • これらのバイアスは、爆発エネルギーや前駆体との対応関係に関する誤った推論へと伝播する。

3. ケーススタディ：SN 1998bw：

   • SN 1998bwに snmix を適用したところ、中程度の混合分率（$f_{\text{mix}} \approx 0.59$）が得られた。一ゾーンの最大尤度推定値は事後分布と重なっているものの、$M_{\text{ej}}$ では $\sim13\%$、$E_k$ では $\sim5\%$ の差があり、ボロメトリックなデータだけでは質量、エネルギー、不透明度、および混合を一意に分解できないことを示している。

4. コラプサーのrプロセス・シグネチャー：

   • rプロセス濃縮のシグネチャーは、必ずしも晩期の近赤外線（NIR）過剰ではない。
   • 配置： 埋没したrプロセス物質は、高い不透明度によって初期の光学的な光度を抑制するが、濃縮された層が透明になった後に初めてNIR過剰が現れる。
   • 幾何学： ディスク・ウィンドモデルにおいて、中性子過剰な物質が赤道面に限定されている場合、オンアックスのGRB-SN（極方向からの視線）では、極方向のSN成分が支配的となるため、rプロセスのシグネチャーはほとんど見られない。信号は、オフアックスの観測者、または物質がより高緯度（$\gtrsim 30^\circ$）まで混合されている場合に最も強くなる。
   • 診断： 色の変化（$\Delta(r-H)$）は有用な診断指標であるが、光学的な抑制とHバンドの増光を分離する必要がある。10〜60日間の光学およびNIRの同時観測は、晩期のHバンドの点データのみよりも情報量が多い。

意義
本論文は、「放射性構造」（加熱と不透明度の半径方向分布）が、超新星の解釈において極めて重要でありながら、しばしば無視されている変数であることを主張している。

• ニッケル駆動型超新星に対して： 高速上昇を示すIc-BL型超新星やGRB-SNに対し、中央エンジンの存在を日常的に持ち出すことに異議を唱え、外側への混合が実行可能な代替案であることを示唆している。また、一ゾーンのフィッティングが「統計的に良好な」フィットを提供しながらも、物理的には誤ったパラメータを与える可能性があると警告している。
• コラプサーに対して： 単なる「NIR過剰」というヒューリスティックを超えた、rプロセス濃縮の探索戦略を洗練させている。赤道面に限定された噴出物がある場合、オンアックスのGRB-SNでの非検出は、rプロセス生成を否定するものではないことを強調している。
• 観測戦略： これらの縮退を打破するには、観測量を組み合わせることが必要であると著者らは結論付けている。低噴出物質量、高混合、あるいはエンジンパワーを区別するためには、光度、速度、色進化、スペクトル、および晩期のボロメトリックな減少を組み合わせることが不可欠である。

本研究は、snmix が簡略化された半解析的なツールではあるものの、半径方向の構造の影響を効果的に分離できることを強調しており、標準的な一ゾーン解釈に対する必要な修正を提供している。