超新星光度曲線シミュレーター

超新星光度曲線シミュレーター

恒星が一生の最後に起こす爆発「超新星 (Supernova)」。Type Ia から II-P・IIn・Ib/Ic まで、放出質量・Ni-56 質量・爆発エネルギー・距離を変えながら、Arnett 解析モデルで光度曲線とピーク絶対等級・Phillips ΔM15 をリアルタイム計算します。

パラメータ設定

超新星タイプ

爆発機構と典型的なエジェクタ組成

放出質量 M_ej

M☉

爆発で吹き飛ぶ全質量。Ia≈1.4、II-P≈10〜20

Ni-56 質量 M_Ni

M☉

爆発で合成される ⁵⁶Ni の量。ピーク光度の主決定因子

爆発エネルギー E_kin

foe

運動エネルギー (1 foe = 10⁵¹ erg)。Ia≈1.5、II≈1.0

放出速度 v_ej

km/s

エジェクタ表面の膨張速度。光球の特性速度

距離 d

Mpc

超新星までの距離。1 Mpc ≈ 326 万光年

観測バンド

B：青色、V：可視（標準）、Bolometric：全波長合計

計算結果

—

拡散時間 τ_diff (day)

—

ピーク光度 (× 10⁹ L☉)

—

絶対等級 M_V (mag)

—

見かけ等級 m_V (mag)

—

ピーク到達時間 (day)

—

Phillips ΔM15 (mag)

—

超新星爆発 — エジェクタ・衝撃波・観測者

中心の Ni-56 領域（白色）と外側の衝撃波面（赤外殻）が膨張し、γ 線が光子拡散を経て観測者に届きます。リング外周の輝度バーが光度曲線の現在位置。

光度曲線 — L(t) [対数スケール]

SN タイプ別ピーク絶対等級 M_V の比較

理論・主要公式

$$L(t) = \int_0^t [\dot E_{Ni} + \dot E_{Co}]\, e^{-(t-t')^2/2\tau_d^2}\,dt',\quad \tau_d = \sqrt{\frac{\kappa M_{ej}}{c\,v}}$$

Arnett (1982) の解析モデル。τ_d は光子拡散時間、κ は不透明度、M_ej は放出質量、v は膨張速度、Ni→Co→Fe 連鎖崩壊で発生する γ 線が放出物質中で熱化して光となる。

$$L_{\rm peak} \approx 2\times10^{43}\,\frac{M_{Ni}}{M_\odot}\;[{\rm erg/s}], \qquad M_V = 4.74 - 2.5\log_{10}(L/L_\odot)$$

Arnett 則による簡略ピーク光度と、それを絶対 V 等級に換算する式。M_Ni = 0.6 M☉ なら M_V ≈ -19.0 となり SN Ia の標準光源値に一致。

$$m - M = 5\log_{10}(d/10\,{\rm pc}), \qquad \Delta M_{15} \simeq 0.7 - 0.5\,(M_{Ni} - 0.6)$$

距離指数（見かけ等級と絶対等級の関係）と Phillips (1993) 関係。ΔM15 はピークから 15 日後までの B バンド減光量で、SN Ia の標準化に用いられる。

超新星光度曲線シミュレーターとは

🙋

「超新星」って、星が爆発する現象ですよね。タイプ Ia とか II とか聞きますが、何がそんなに違うんですか？

🎓

爆発の「きっかけ」が根本的に違うんだ。Type Ia は太陽くらいの質量の白色矮星が、伴星から物質を吸い込んでチャンドラセカール限界 (≈1.4 M☉) に達した瞬間、芯の炭素が一気に核融合して星全体を吹き飛ばす「熱核暴走」。一方 Type II・Ib/Ic は質量 8 M☉ 以上の重い星が、鉄の芯で核融合エネルギーを取り出せなくなり、自重で内側に潰れる「重力崩壊」。崩壊した中心は中性子星かブラックホールになり、跳ね返りの衝撃波で外層を吹き飛ばすんだ。同じ「超新星」でも、片や白色矮星の最後、片や巨星の最後、と全然違う星の死に方なんだよ。

🙋

それなのに、明るさは Type Ia がほぼ一定で「標準光源」になるって聞きました。なぜ揃うんですか？

🎓

そこが SN Ia の魔法のところ。爆発のスタート条件 ── チャンドラセカール限界の白色矮星 ── がほぼ同じだから、合成される放射性 Ni-56 の量がだいたい 0.4〜0.8 M☉ に揃うんだ。光度曲線のピーク付近は、その Ni-56 が Co-56 → Fe-56 と崩壊して γ 線を出し、それが熱化して光に変わるエネルギーで決まる。Arnett の経験則で L_peak ≈ 2×10^43 erg/s × M_Ni/M☉ となるから、左で M_Ni=0.6 にすると約 3×10^9 L☉、絶対 V 等級にして -19.0 等。これが宇宙のどこで起きてもほぼ同じ値だから、見かけの明るさを測れば距離が分かる。これがハッブル定数や宇宙の加速膨張（1998 年 Riess・Perlmutter・Schmidt → 2011 年ノーベル賞）の決定に使われた仕組みだよ。

🙋

面白い！じゃあ「Arnett 則」って何ですか？理論メモにも τ_d = √(κM/cv) って出てきますね。

🎓

Arnett (1982) が打ち立てた解析モデルで、超新星の光は「Ni-56 → Co-56 → Fe-56 の崩壊エネルギー」を「光子拡散時間 τ_d」で時間的にぼかしたもの、と考える。τ_d は光子が散乱されながら表面に脱出する時間で、不透明度 κ・放出質量 M・膨張速度 v に依存する。SN Ia の典型値（κ≈0.1 m²/kg, M=1.4 M☉, v=10⁴ km/s）で τ_d ≈ 20 日。これがそのままピーク到達時間（rise time）に対応するんだ。τ_d が大きいと光がじわじわ拡散するからピークが遅く・低く・幅広に、τ_d が小さいと早く・高く・鋭くなる。II-P 型のように水素を保ったまま潰れる赤色超巨星は、放出質量 10〜20 M☉ で τ_d が長く、しかも水素再結合で光球が一定温度に止まる「プラトー」が 100 日続くのが特徴だよ。

🙋

Phillips 関係 ΔM15 って何ですか？標準光源が「標準化される」ってどういうことです？

🎓

Phillips (1993) が発見した経験則で、SN Ia は「ピークの絶対等級」と「ピークから 15 日後までの B バンド減光量 ΔM15」が強い相関を持つ。明るい SN Ia ほどゆっくり暗くなり、暗い SN Ia ほど早く暗くなる、という法則。これで個々の SN Ia を ±0.15 mag、距離にして ±7% の精度で揃えられる。Riess らはこの ΔM15 補正を使って高赤方偏移 SN Ia を見ることで、宇宙が加速膨張しているのを突き止めたんだ。現代の LSST/Rubin Observatory・JWST・Zwicky Transient Facility (ZTF) では毎晩数百個の SN を検出して、Phillips 関係をさらに精密化しているよ。ちなみに観測史上の有名な超新星は SN 1987A（大マゼラン雲 LMC、ニュートリノ初検出）、SN 2011fe（風車銀河 M101、最も観測された Ia）、SN 2014J（M82、近傍 Ia）あたりだね。

よくある質問

熱核暴走型（Type Ia）と重力崩壊型（Type II, Ib/Ic）のいずれも、ピーク付近の光は爆発で合成された Ni-56 が Co-56 → Fe-56 と崩壊する際に放出される γ 線を、放出物質中で熱化したものです。Arnett の経験則ではピーク光度 L_peak ≈ 2×10^43 erg/s × (M_Ni / M_sun) と書け、Ni-56 質量 0.6 M☉ のとき L_peak ≈ 3.14×10^9 L_sun、絶対 V 等級にして約 -19.0 等になります。これが SN Ia を宇宙論距離指標として使える物理的根拠です。

放出ガス（エジェクタ）の中で光子が散乱されながら表面に脱出するまでの時間スケールで、Arnett (1982) は τ_diff = √(κ M_ej / (c v)) と与えました。κ は不透明度、M_ej は放出質量、v は膨張速度です。SN Ia の典型値（M_ej=1.4 M☉, v=10⁴ km/s, κ≈0.1 m²/kg）で τ_diff ≈ 20 日となり、これがピーク到達時間とほぼ一致します。II-P 型のように M_ej が大きい超新星はピークが遅れ、剥がれ星 Ib/Ic は M_ej が小さいので早く明るくなります。

Phillips (1993) が発見した SN Ia の経験則で、ピーク等級 M_peak と「ピークから 15 日後までの B バンド減光量 ΔM15」が強い相関を持つというものです。明るい SN Ia ほど減光が遅く、暗い SN Ia ほど早く減光します。この関係で個々の SN Ia のピーク等級を ±0.15 mag 程度の精度で標準化でき、Riess・Perlmutter・Schmidt の宇宙加速膨張発見（2011 年ノーベル物理学賞）の鍵になりました。本ツールでは ΔM15 ≈ 0.7 - 0.5·(M_Ni - 0.6) という簡略式で表示します。

はい、距離が分かって初めて意味が出ます。見かけ等級 m と絶対等級 M、距離 d (pc) の間には m - M = 5·log10(d/10pc) の距離指数（distance modulus）が成り立ちます。SN Ia は M_V ≈ -19.0 がほぼ一定なので、観測で m_V を測れば d が逆算でき、これが「標準光源」として宇宙論的距離測定に使われる仕組みです。10 Mpc で約 11 等、100 Mpc で約 16 等、1 Gpc で約 21 等となり、ハッブル望遠鏡や Rubin/LSST、JWST の検出可能領域に直結します。

実世界での応用

宇宙論距離測定とハッブル定数：SN Ia の絶対等級がほぼ一定（M_V ≈ -19.0）であることを利用して、近傍銀河から赤方偏移 z ≈ 1 を超える遠方銀河までの距離を測ります。SH0ES プロジェクトはセフェイド変光星と SN Ia を組み合わせて H_0 ≈ 73 km/s/Mpc を導き、Planck 衛星の CMB 由来値 H_0 ≈ 67 と対立して「ハッブル・テンション」と呼ばれる現代天文学の最大論争を生んでいます。SN Ia 1 つの精度がそのまま宇宙論パラメータの不確実性に直結する場面です。

大規模時間領域サーベイ（LSST/ZTF/Rubin）：Vera C. Rubin Observatory (LSST) は 2025〜2035 年の 10 年運用で毎晩 2,000 件以上の超新星候補を検出します。ZTF は既に 5,000 件超の SN を分類済み。これら大量データから光度曲線パラメータ（rise time, ΔM15, ピーク等級）を統計処理することで、宇宙加速膨張のさらなる精密測定、暗黒エネルギー状態方程式 w の制限、変則的な SN サブクラス（SN Iax, super-Chandra など）の発見が進みます。

元素合成と銀河化学進化：宇宙の鉄の約半分は SN Ia 由来、酸素・マグネシウムなどの α 元素は主に重力崩壊型 SN 由来とされます。本ツールで M_Ni を変えるとピーク光度が変わりますが、それは宇宙に放出される鉄（Ni-56 → Fe-56）の量を直接示します。銀河の星種族の年齢・金属量分布は、SN Ia と core-collapse SN の発生率の歴史によって決まり、銀河形成シミュレーション（IllustrisTNG, EAGLE）の基本入力になります。

マルチメッセンジャー天文学：SN 1987A は光・ニュートリノ・電波・X 線で多波長観測された最初の超新星で、Kamiokande・IMB が捉えた 24 個のニュートリノが、重力崩壊型 SN の核ニュートリノ放出を初めて実証しました。次に銀河系内で爆発する SN では、SK-Gd・KamLAND・IceCube・JUNO によって 10^4 個級のニュートリノが検出される見込みで、コア崩壊物理の精査が期待されています。

よくある誤解と注意点

まず最大の誤解が、「SN Ia のピーク光度は完全に一定」という思い込み。実際は ±0.4 mag 程度のばらつきがあり、Phillips 関係（ΔM15 補正）・色補正（SALT2 など）を適用して初めて ±0.15 mag に揃います。さらに「サブ-Chandra」や「super-Chandra」と呼ばれる非標準 SN Ia があり、これらは標準化を逸脱します。宇宙論解析ではこれらを慎重に除外しないと H_0 や w の値がバイアスを受けます。本ツールの簡略 Phillips 式 (ΔM15 ≈ 0.7 - 0.5(M_Ni-0.6)) はあくまで概念理解用で、実観測の SALT2 / SNooPy などのフィッターはもっと複雑です。

次に、「Arnett モデルは万能ではない」こと。Arnett (1982) は「光球が球対称・均質・光学的に厚い」「ピーク付近では光子拡散時間 τ_d で崩壊エネルギーが時間的に平均化される」という仮定を置きます。これは Type Ia のピーク付近では良い近似ですが、Type II-P のプラトー期（水素再結合）、Type IIn の CSM 相互作用（外側の星周物質と衝突して光る）、Type Ib/Ic で観測される非対称性（GRB 関連）には適用できません。本ツールでは Ia 以外も同じ枠組みで扱っているので、II-P/IIn/Ibc は概算と理解してください。

最後に、「見かけ等級は距離だけで決まらない」こと。実際の観測ではホスト銀河・天の川銀河内の星間塵による減光（extinction、E(B-V)）が加わります。SN 2014J は M82 のダストレーンで起きたため、補正前は約 1.6 mag 暗く見えました。さらに高赤方偏移 SN では K-correction（観測バンドと静止系バンドのズレ）、宇宙論 K-correction、銀河間ダストなど多段の補正が必要です。本ツールの距離指数 m - M = 5log(d/10pc) は最もシンプルな「真空中の幾何学的減衰」のみを扱っていることに注意してください。