超新星光度曲线模拟器

参数设置

超新星类型

爆炸机制和典型喷出物组成

喷出质量 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 (天)

—

峰值光度 (× 10⁹ L☉)

—

绝对星等 M_V (mag)

—

见光星等 m_V (mag)

—

峰值到达时间 (天)

—

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 的标准化。

超新星光度曲线模拟器介绍

🙋

「超新星」就是星星爆炸对吧？听说有 Type Ia、Type 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 mag。这个值在宇宙任何地方都差不多，所以只要测见光亮度，就能算出距离。这是发现宇宙加速膨胀（1998 年 Riess、Perlmutter、Schmidt → 2011 年诺贝尔奖）的核心工具。

🙋

那「Arnett 法则」是什么？理论备忘录里还有 τ_d = √(κM/cv)，这些参数什么意思？

🎓

Arnett (1982) 建立的分析模型，核心思想是「超新星的光来自 Ni-56→Co-56→Fe-56 衰变释放的能量，但光子在膨胀的喷出物中散射延迟，延迟时间就是光子扩散时间 τ_d」。τ_d 等于 √(κM/cv)，其中 κ 是不透明度、M 是喷出质量、v 是膨胀速度。SN Ia 典型值（κ≈0.1 m²/kg，M=1.4 M☉，v=10⁴ km/s）给出 τ_d ≈ 20 天。这个时间就对应从爆炸到峰值的上升时间。τ_d 大了光扩散慢，峰值就迟、低、宽；τ_d 小了峰值就早、高、尖。II-P 型的红色超巨星喷出质量 10～20 M☉，所以 τ_d 很长，峰值会拖到数周后；剥离星 Ib/Ic 喷出质量小，增亮很快。

🙋

Phillips 关系 ΔM15 又是什么？说是「标准化」，怎么理解？

🎓

Phillips (1993) 发现的经验规律，SN Ia 有个很强的规律：「峰值星等」和「从峰值到第 15 天的 B 波段减光量 ΔM15」成相关。亮的 SN Ia 减光缓慢，暗的减光快速。这个关系让每个 SN Ia 都能标准化到 ±0.15 mag、距离精度 ±7%。Riess 他们就用这个 ΔM15 修正去观测远处的 SN Ia，从而发现了宇宙在加速膨胀。现在的 LSST/Rubin、JWST、ZTF 每晚检测数百个 SN，继续精化 Phillips 关系。著名的超新星有 SN 1987A（大麦哲伦云，首次中微子检测）、SN 2011fe（M101 风车星系，观测最全面）、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 mag。这是 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 型因喷出质量大，峰值延迟；剥离星 Ib/Ic 因喷出质量小，增亮较快。

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 的精度直接影响宇宙学参数的不确定性。

大规模时间巡天（LSST/ZTF/Rubin）：Vera C. Rubin 天文台（LSST）在 2025～2035 年的 10 年运行中每晚将检测 2000+ 个超新星候选体。ZTF 已归档 5000+ 超新星。通过这些大数据统计处理光度曲线参数（上升时间、ΔM15、峰值星等），可进一步精化宇宙加速膨胀测量、限制暗能量状态方程 w、发现异常超新星亚类（SN Iax、超钱德极限等）。

元素合成与星系化学演化：宇宙铁元素约一半来自 SN Ia，氧、镁等 α 元素主要来自重力崩壊型 SN。本工具改变 M_Ni 时峰值光度改变，直接反映喷向宇宙的铁（Ni-56 → Fe-56）量。星系星族的年龄、金属量分布由 SN Ia 和重力崩壊 SN 的历史发生率决定，是星系形成模拟（IllustrisTNG、EAGLE）的基本输入。

多信使天文学：SN 1987A 首次在光、中微子、电波、X 射线多波段被观测，其中 Kamiokande、IMB 捕获的 24 个中微子首次证实重力崩壊型超新星的核中微子释放。下一次银河系内超新星爆炸时，SK-Gd、KamLAND、IceCube、JUNO 有望检测 10⁴ 级中微子，深入研究核心崩壊物理。

常见误解与注意事项

最大的误解是认为「SN Ia 峰值光度完全恒定」。实际上散布约 ±0.4 mag，需要用 Phillips 关系（ΔM15 修正）、色修正（SALT2 等）才能标准化到 ±0.15 mag。还有「亚钱德拉极限」和「超钱德拉极限」等非标准 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 的「环星物质碰撞」（外层星周物质发光）、Type Ib/Ic 的观测非对称性（GRB 关联）都不适用。本工具对 Ia 外的类型也套用同一框架，应视为粗略估算。

最后，「见光星等不仅决定于距离」。实际观测中还有银河系和宿主星系的星际尘埃消光（extinction，E(B-V)），会额外减光。SN 2014J 在 M82 的尘埃线中爆发，修正前暗 1.6 mag。高红移超新星还要加 K-correction（观测波段与静止系波段的转换）、宇宙 K-correction、星系间尘埃等多层修正。本工具的距离模数 m - M = 5log(d/10pc) 仅是「真空中的几何衰减」，需注意这一点。