Stellar Interpretation of Meteoritic Data and PLotting for Everyone (SIMPLE): Isotope Mixing Lines for Six Sets of Core-Collapse Supernova Models

技術要約：SIMPLE – 全ての人に向けた、隕石データの恒星学的解釈とプロットのためのツール

問題提起
バルク隕石とその包含物（インクルージョン）は、多くの化学元素において同位体変異性を示しており、これは太陽系形成以前の恒星や超新星における核合成イベントのシグネチャーである。これらの変異を特定することは太陽の誕生環境に関する知見を提供するが、その正確な天文学的起源については依然として議論の余地がある。この問題を解決する上での主要な課題は、隕石の同位体データと恒星核合成モデルによる予測との比較である。特に、前駆体の寿命が短く、プレソーラー塵の生成に重要な役割を果たすため、コア崩壊型超新星（CCSN）が特に関心の対象となっている。しかし、CCSNモデルの予測を隕石データと比較することは、内部正規化によって導入される非線形効果、データ処理の複雑さ、および多くの現代的なCCSNモデルにおいて詳細な組成プロファイル（放出質量に対する関数としての組成）へのアクセスが困難であることによって困難を極めている。歴史的には、少数の古いモデル（例：Meyer et al. 1995; Rauscher et al. 2002）のみがこのような解析に使用されてきたが、これでは20年間のCCSN研究の進展を捉えきれていない。さらに、化学的な質量依存の分別と放射性崩壊を区別するには、モル分率と質量分率の厳密な扱い、および不安定同位体の特定の処理が必要となる。

手法
これらの課題に対処するため、著者らは、CCSN噴出物の同位体組成を隕石データと比較するために設計されたオープンソースのPythonツールキットであるSIMPLE（Stellar Interpretation of Meteoritic Data and Plotting for Everyone）を開発した。その手法は以下の通りである：

1. データの統合： 本ツールキットは現在、6つの異なる公開セット（Ra02, Pi16, Si18, Ri18, LC18, La22）からなる18のCCSNモデルを取り込んでおり、これらは太陽金属量における初期質量15, 20, 25 $M_\odot$ をカバーしている。これらのモデルは、様々な前駆体コード（Kepler, MESA, GENEC, FRANEC）および爆発メカニズム（ピストン駆動、セドフ・ブラストウェーブ、キネティック・ボム、ニュートリノ駆動）を利用している。
2. データ処理： コードは、放出質量（mass fraction）の関数としての同位体存在量（質量分率）を読み取る。これは、比率計算における質量数（$A_i$）因子による人工的な異常を防ぎ、標準的な太陽/地球の存在量との整合性を確保するために、モル分率（$Y_i$）へと変換される。
3. 放射性崩壊の取り扱い： ユーザーは、放射性核を省略するか、あるいはそれらの全存在量を安定した娘核に加えるかを選択できる。コードは、CCSN噴出物の短いタイムスケール（爆発後数時間）と、ダスト形成（数ヶ月/数年）の間の差異を考慮しており、短寿命の親核（例：$^{60}$Co）を含めることができる一方で、$^{60}$Feのようなより長寿命の親核に対する柔軟性も提供している。
4. 正規化と希釈： 隕石データと比較するために、コードは混合線を計算する。コードは、噴出物を太陽組成で希釈し、参照同位体比を用いて内部正規化を行う。コードは、正規化のための標準的な指数関数則（$\epsilon$ 値の計算）を実装しており、代替案として線形近似法（Dauphas et al. 2004）も含んでいる。また、核合成が最小限の領域における偽の結果を避けるため、最大のオフセットを $\pm 1$ $\epsilon$ 単位にスケールするための希釈係数を自動的に決定する。
5. 可視化： 本ツールは、放出質量に対する同位体比のプロットを生成し、Meyer et al. (1995) の分類に従って、噴出物を主要な燃焼生成物（例：Si, O/Si, O/C）に基づくゾーンに分割する。また、混合線を特定するために、同位体比間の相関（例：$\epsilon^{64}\text{Ni}$ 対 $\epsilon^{62}\text{Ni}$）をプロットする。

主な貢献

• ソフトウェアのリリース： 多様なCCSNモデルの出力を標準化する、コラボレーティブなオープンソースのPythonパッケージ（GitHubおよびPyPIで入手可能）としてのSIMPLEのリリース。
• 包括的なデータセット： 2002年から2022年までのモデルを網羅する6つの独立した研究グループからの18のモデルを統合し、CCSN核合成における系統的な不確実性とロバストな特徴に関するより広い視野を提供する。
• 手法の厳密性： 質量分率とモル分率の一貫した変換、および同位体比計算における人工的な誤差を防ぐための放射性崩壊鎖の柔軟な取り扱いの実装。
• ケーススタディへの適用： 全18モデルを用いたニッケル（Ni）同位体の詳細な分析により、これまでの知見を再現し拡張する本ツールの能力を実証。

結果
Ni同位体をケーススタディとして、著者らは全6セットの20 $M_\odot$ モデルを分析した：

• ゾーンのシグネチャー： 分析により、異なる噴出物ゾーンにおける明確な同位体シグネチャーが明らかになった。Siに富むゾーンは、$^{58}\text{Ni}/^{61}\text{Ni}$ に内部正規化された $\epsilon^{64}\text{Ni}$ 対 $\epsilon^{62}\text{Ni}$ の傾きが、一貫して約3.0となる。この結果は、Siゾーンが噴出される限り、異なる前駆体質量およびモデルセットにわたってロバストである。
• 隕石との比較： バルク隕石データに見られる観測された傾き（$3.003 \pm 0.166$）は、CCSN噴出物の最も内側のSiに富む領域に由来する物質と適合している。
• モデルの変動性： Siゾーンは一貫した結果をもたらすが、他のゾーン（例：O/Si, O/C, He/C）は、特定の爆発メカニズム、C-Oシェル合体の発生、および中性子捕獲過程（例：$^{22}\text{Ne}(\alpha,n)^{25}\text{Mg}$ 反応）の活性化に応じて、傾きに大きな変動を示す。
• 放射性の寄与： 本ツールは、放射性親核（具体的には $^{60}$Co による $^{60}$Fe）が安定な $^{60}\text{Ni}$ の存在量に与える寄与を定量化することに成功し、崩壊のタイミングに関する仮定の違いが最終的な同位体組成にどのように影響するかを明らかにしている。

意義
本論文は、幅広い恒星モデルと比較するための標準化された柔軟なプラットフォームを提供することにより、SIMPLEが隕石データの解釈を促進すると主張している。バルク隕石のNi異常がCCSN噴出物のSiに富む領域の物質と適合することを確認することで、本研究は、そのような超新星が初期太陽系の不均質性に寄与したという仮説を強化している。著者らは、本ツールが恒星核合成および隕石研究の発展を促進するためのコラボレーティブなプラットフォームとして設計されていることを強調している。また、現在のバージョンは鉄コア崩壊型超新星に焦点を当てているものの、フレームワークは将来的なモデル（3Dシミュレーション、電子捕獲型超新星、および他の恒星源（例：AGB星、ノーヴァ）を含む）や、異なる噴出層の混合機能を取り込めるように構築されている。