Constraining hydrodynamic model of nearby type IIP SN 2023ixf

技術要約：近傍のII型超新星SN 2023ixfにおける流体モデルの制約

問題提起
近傍に位置するII型IIP型超新星SN 2023ixfは、膨大な観測データが利用可能であるにもかかわらず、これまでに発表された流体力学モデルは、信頼できる一連の爆発パラメータの収束に失敗している。既存の研究（例：Bersten et al., 2024; Moriya and Singh, 2024）は、推定される爆発エネルギー（$0.5–3 \times 10^{51}$ erg）や放出質量に幅広い範囲をもたらしており、現在のパラメータ導出の信頼性を損なっている。この「パラメータの縮退」は、標準的な観測量（光度曲線のプラトー輝度や光球速度など）だけでは、モデルを一意に制約するには不十分であることから生じている。さらに、従来のモデルは初期の輝度ピークを説明するために、大規模な星周（CS）シェル（$\sim 0.4–0.5 M_\odot$）を想定することが多かったが、これはすべてのデータによって必ずしも支持される特徴ではない。

手法
著者らは、ラグランジュ放射流体計算コードCRABを用い、SN 2022ixfの流体力学モデリングを再検討している。このパラメータの縮退を解消するため、本研究では、先行研究よりも広く、かつ具体的な観測的制約を組み込んでいる：

1. 初期最大放出速度： 初期のスペクトル（8〜14日目）で検出された浅いバルマー吸収を利用して、最大膨張速度（$v_{sn}$）を決定する。著者らは、不透明な放出核と不透明な星周シェル（CDS）の間の層におけるこれらの吸収の形成を、トムソン散乱とCDSの光学的厚さを考慮してモデル化している。
2. $^{56}$Ni混合範囲： JWSTのスペクトル（252日目）における[Co II] 10.521 $\mu$m放射線の青側幅から導出された、$^{56}$Ni放出体の速度範囲を組み込んでいる。
3. 星周構造： H$\alpha$プロファイル（トムソン散乱の翼部）およびX線データから推測された柱密度（$N_H$）を用いて、CSシェルの密度と範囲を制約している。
4. 放射流体力学： 前方ショックとCS媒質との相互作用をシミュレートし、放射支配型のショックから物質支配型のショックへの遷移をモデル化し、特に初期の硬X線放出の出現を扱っている。

前超新星（pre-SN）モデルは静水圧的な赤色超巨星（RSG）に基づいているが、水素外層におけるより急峻な密度勾配と、初期電離フラッシュスペクトルと矛盾しない低質量CSシェルを含むように修正されている。

主な結果
最適化された流体力学モデルは、SN 2023ixfに対して以下の具体的なパラメータを算出している：

• 爆発エネルギー： $2.8 \times 10^{51}$ erg ($\pm 0.4 \times 10^{51}$ erg)。
• 放出質量： $13.2 M_\odot$ ($\pm 1.3 M_\odot$)。
• 前超新星半径： $1540 R_\odot$ ($\pm 150 R_\odot$)。
• $^{56}$Ni質量： $0.07 M_\odot$ ($\pm 0.004 M_\odot$)。
• $^{56}$Ni混合範囲： 最大 $4000$ km s$^{-1}$ まで。
• 星周物質： 高密度の内部CSシェル（質量 $\approx 0.01 M_\odot$、半径 $\approx 5 \times 10^{14}$ cm）と希薄な外部風で構成される。内部シェルは密度パラメータ $w_{int} = 3.7 \times 10^{16}$ g cm$^{-1}$ を持ち、外部風は $w_{ext} = 1.7 \times 10^{15}$ g cm$^{-1}$ である。

このモデルは、CDSの光学的厚さとロスランド不透明度の相互作用による5日目付近の短い輝度極小を含む、ボロメトリック光度曲線を再現することに成功している。また、光球速度の進化も再現している。導出された前超新星質量 $14.8 M_\odot$（放出体 + 中性子星）は、親星の推定値である $16–17 M_\odot$ と一致している。

初期の硬X線に関しては、放射流体シミュレーションにより、H$\alpha$とX線柱密度から導出された一様なCSシェル密度は、5日目のX線光度を過大評価することが示された。より良い適合は、CSシェルが塊状（クランピー）構造を持つと仮定することで得られる。すなわち、希薄な間雲媒質（$\eta \approx 0.1$）が硬X線放出を許容し、一方で高密度の塊（クランプ）がX線吸収データと一致するような遮蔽を提供するという構造である。

意義および主張
著者らは、本研究が、爆発から硬X線放射の出現に至るまでの、前方ショックの放射支配型から断熱的な物質支配型への遷移を明示的にモデル化した、SN 2023ixfに関する初の流体力学シミュレーションであると主張している。

本研究の主要な意義は以下の通りである：

1. パラメータの縮退の解消： 初期のバルマー吸収から導出された初期最大放出速度と$^{56}$Niの速度範囲を優先することで、著者らは従来のモデルと比較して、爆発エネルギーと放出質量の不確実性を大幅に絞り込んだ。
2. CSシェルモデルの改訂： 本研究は、大規模なCSシェルは不要であり、代わりに、塊状構造を持つことを前提とすれば、低質量（$\sim 0.01 M_\odot$）のシェルで初期のスペクトル特徴とX線データを説明するのに十分であることを示している。
3. ショック相互作用の放射流体力学： 本論文は、高密度なCS環境におけるショック遷移メカニズムの詳細な記述を提示しており、これは放射流体力学の枠組みにおいてこれまで記述されていなかった現象である。

著者らは、導出された爆発エネルギー（$2.8 \times 10^{51}$ erg）が高い値であり、標準的なニュートリノ駆動型の爆発メカニズムのエネルギー限界を超える可能性があることを指摘しており、爆発物理学または親星の構造における潜在的な複雑性を示唆している。