CIGaRS I: Combined simulation-based inference from type Ia supernovae and host photometry

技術的サマリー：CIGaRS I – タイプ Ia 超新星と宿主銀河の光度測定からの統合シミュレーションベース推論

問題提起
タイプ Ia 超新星（SNe Ia）は主要な宇宙論的プローブとして機能するが、その有用性は、明るさと宿主銀河の性質（例えば、恒星質量、金属量、年齢、星形成率）との間の経験的相関によって制限されている。現在の宇宙論解析では、通常、まず「質量ステップ」のような宿主恒星質量に基づくアドホックな方法を用いて超新星の観測量を補正し、その後に宇宙論的フィッティングを行うという「二段階」アプローチが採用されている。この分離は、エディントンバイアス（大質量の宿主への過剰な重み付け）や、塵の消光や選択バイアスなどの外的効果から本質的な物理的依存関係を解きほぐせないことなど、統計的バイアスを導入する。さらに、分光赤方偏移への依存は、LSST（Legacy Survey of Space and Time）のような将来の大規模サーベイにおいて、分光フォローアップが得られる超新星と宿主がごく一部しか存在しないため、サンプルサイズを制限している。

手法：CIGaRS フレームワーク
著者は、CIGaRS（Combined Inference and Galaxy-related Standardization）を提示する。これは、純粋な光度測定データからの同時推論を行うように設計された統合的なベイズ階層モデルである。このフレームワークは、以下の 3 つの中核コンポーネントを統合している：

1. 統合前方シミュレータ：このモデルは、銀河進化、星形成履歴、化学的増殖のためのProspector-βと、SNe Ia の光度曲線要約のためのSimple-BayeSNを組み合わせる。これらは以下の要素を介して連結される：

   • 遅延時間分布（DTD）： progenitor（前身星）の年齢の関数としての SNe Ia 発生率を支配するべき乗則パラメータ化。
   • 物理的宿主 - 超新星の結合： progenitor の金属量（$Z^*$）と年齢（$t^*$）に対する SNe Ia の本質的等級の明示的なパラメータ化された依存関係、および残余の「質量ステップ」（$\Delta M$）。
   • 塵と選択効果：シミュレータは宿主銀河と SNe 両方の塵消光をモデル化し、宇宙論的距離モジュールを適用し、現実的な観測選択関数（LSST の wide-fast-deep サーベイ条件を模倣）を組み込む。

2. シミュレーションベース推論（SBI）：超新星あたり数十個の潜在パラメータという膨大な数の潜在パラメータに対して、扱いにくい尤度関数を構築する代わりに、著者は切断された周辺化ニューラル比推定（TMNRE）を採用する。条件付き deep set++ ニューラルネットワークは、パラメータと観測量のペアからなるシミュレーションデータ上で訓練され、グローバル/集団パラメータとオブジェクト固有の潜在変数の事後分布を同時に近似する。

3. 入力データ：推論ネットワークは、宿主銀河の広帯域光度測定（$ugrizy$）と、SNe Ia の要約統計量（ピーク明るさ、ストレッチ、色、および関連する不確かさ）のみを利用する。分光赤方偏移や遷移現象の分類データは仮定されないが、シミュレーションにおいては完全な遷移現象の分類が仮定される。

主要な貢献

• 統合モデリング：本論文は、標準化と宇宙論的解析を分離することなく、宇宙論的パラメータ、DTD、宿主 - 超新星の物理的相関、およびオブジェクト固有の性質（赤方偏移、質量、金属量）を共同で推論する自己整合的なフレームワークを導入する。
• 混同の解離：このモデルは、SNe Ia の本質的な金属量依存性と年齢依存性を区別する能力を実証する。これらは、宿主恒星質量分布において（滑らかなステップ対線形トレンドとして）異なって現れる。
• スケーラビリティ：ニューラル SBI を活用することで、このフレームワークは階層モデルの高次元パラメータ空間を処理し、従来のマルコフ連鎖モンテカルロ法の計算上の扱いにくさなしに、約 1,600 から約 16,000 のオブジェクトへ効果的に拡張する。

結果
現在の（DES）および将来の（LSST）サーベイ規模を代表するモックデータセットを用いて：

• パラメータ回復：このフレームワークは、DTD 傾き（$b$）を$\pm 0.05$以内に回復し、金属量（$\gamma_{Z^*}$）および年齢（$\gamma_{t^*}$）の相関係数を解離することに成功した。
• 光度赤方偏移：このモデルは、中央値のばらつきが$\sim 0.01$の光度赤方偏移を推論し、標準的な LSST の光度測定みの推定（$\sim 0.04$）を大幅に上回り、より広い UV - 光学 - IR カバレッジを持つ解析の結果と一致する。この精度は、Prospector-βの事前分布と共同宇宙論的推論を通じて質量 - 年齢 - 赤方偏移の混同を破ることで達成される。
• 宇宙論的制約：SNe と宿主の光度測定の共同解析は、分光赤方偏移を持つオブジェクトの小部分に依存する解析と比較して、宇宙論的制約（$\Omega_{m0}$–$\Omega_{\Lambda0}$平面内）を約 4 倍改善する。
• 質量ステップ対物理的相関：シミュレーションは、金属量相関が$10^{10} M_\odot$にわたって観測される「質量ステップ」を模倣できることを確認し、年齢相関は明確な線形トレンドを生み出すことを示した。これは、従来の質量ステップが、背後にある物理的駆動力の代理指標である可能性を示唆している。

意義
著者は、CIGaRS が LSST 時代における「光度測定データの全能力を解き放つ」と主張している。アドホックな補正や二段階解析から離れることで、このフレームワークは、選択、塵、赤方偏移、宿主の性質といったすべての不確かさを最終的な宇宙論的制約へと伝播させる、厳密なエンドツーエンドのシミュレーションベースパイプラインを提供する。このアプローチは、将来のサーベイから期待される大規模で主に光度測定データから堅牢な科学的洞察を抽出するためのスケーラブルな解決策として提示されており、SNe Ia の形成メカニズムとその環境依存性に関する長年の論争を解決する可能性がある。著者は、現在のモデルが特定の仮定（例えば、固定された塵の法則や恒星の初期質量関数）によって制約されていることに言及しつつも、SBI フレームワークは、将来のイテレーションにおいてそのような「既知の未知」を nuisance 変数として体系的に含めることを可能にすると指摘している。