カシオペアAを例に超新星爆発のエネルギー収支と宇宙線伝搬について簡単にまとめてみる

はじめに

この記事では、カシオペアAを例に超新星爆発のエネルギー収支と宇宙線伝搬について簡単にまとめてみます。

超新星残骸 カシオペアA（Cassiopeia A, Cas A） は、

• 若い（年齢 ≈ 340 yr）
• 明るい
• X線・γ線・電波・赤外の multi-messenger 的データが豊富
• Hadronic 宇宙線加速の候補でもある

という理由から、超新星爆発の エネルギー収支 と 宇宙線伝搬 を理解する上で最適な天体の一つです。

この記事では、Cas A を例にしながら

1. 宇宙線のエネルギーは銀河規模でどのくらい必要？
2. 超新星爆発ではどの物理がどれだけのエネルギーを担う？
3. なぜニュートリノが「爆発エネルギーの99%」を持ち出すのか？
4. 親星の初期質量・イジェクタ質量・中性子星質量の関係は？
5. Cas A からの宇宙線電子・陽子は地球に届く？
6. “直接”届くとはどういう意味か？

を一気に整理します。

1. 銀河宇宙線のエネルギー収支：超新星 100 年に 1 回で賄える？

銀河円盤の宇宙線のエネルギー密度はだいたい

1 ~ {\rm eV/cm}^3 ≈ 1.6\times10^{-12} ~{\rm erg}/cm^3。

これを銀河の有効体積に掛けると、

E_{\rm CR,tot} \sim 10^{55}~{\rm erg}

宇宙線の銀河内滞在時間はおよそ

\tau_{\rm res} \sim 10^{7}~{\rm yr}

なので、銀河が宇宙線を維持するために必要な“注入パワー”は、

L_{\rm CR} \sim \frac{10^{55}}{10^{7}} \sim 10^{48}~{\rm erg/yr}
\approx 3\times10^{40}~{\rm erg/s}

一方で、典型的な超新星の運動エネルギー(補足参照)は

E_{\rm SN} \sim 10^{51}\ {\rm erg}

超新星爆発が 100 年に 1 回 起きるとすると：

\frac{10^{51}}{100} = 10^{49}~{\rm erg/yr}

その 10% が宇宙線に変換されると仮定すると：

10^{49} \times 0.1 = 10^{48}~{\rm erg/yr} \approx 3\times10^{40}~{\rm erg/s}

→ 銀河宇宙線に必要なパワーとピッタリ一致
（これが “SNR が銀河宇宙線の主供給源” と言われる理由）

"ピッタリ一致"の意味は、超新星爆発の全エネルギーの99%はニュートリノが持ち運ぶので、超新星爆発のエネルギーの1%(噴出物の運動エネルギー)が100年に一回の頻度で供給されると宇宙線のエネルギーが説明可能、という意味。

2. 超新星爆発のエネルギー収支の基本：なぜニュートリノが 99%？

コア崩壊型超新星では、最終的に中心に 中性子星 (NS) が残ります。

中性子星の質量と半径が決まれば、その“重力束縛エネルギー”は

E_{\rm bind} \sim \frac{3}{5}\frac{GM^{2}}{R}
\sim (2\text{–}4)\times10^{53}\ {\rm erg}

これは 太陽の 100 万年の光度の総量を 10 秒で放つような天文学的エネルギーですが、
このうち 99%以上がニュートリノとして放出されます。

補足：なぜ $E_{\rm bind} \sim \frac{3}{5}\frac{GM^2}{R}$ なのか？

一様密度の球（質量 $M$、半径 $R$） を、無限遠から少しずつ集めて作ると考える。

半径 $r$ まで出来上がったときの内部質量

M(r) = M\left(\frac{r}{R}\right)^3

その外側から薄い殻 (dm)（面密度一定）を持ってくるときの仕事は

dE = - \frac{G M(r)~dm}{r}.

ここで

dm = 4\pi r^2 \rho~dr, \quad
\rho = \frac{3M}{4\pi R^3}.

これを代入すると

dE = - \frac{G~ M(r)~ 4\pi r^2 \rho}{r}dr
= - 4\pi G \rho M(r) r~dr
= - 4\pi G \rho M\left(\frac{r}{R}\right)^3 r~dr.

さらに $\rho$ を代入して積分すると

E_{\rm bind}
= \int_0^R dE
= - \int_0^R \frac{3GM^2}{R^3} \left(\frac{r^4}{R^3}\right) dr
= - \frac{3GM^2}{R^6} \int_0^R r^4 dr
= - \frac{3GM^2}{R^6} \left[\frac{R^5}{5}\right]
= -\frac{3}{5}\frac{GM^2}{R}.

符号は「束縛エネルギー（負のポテンシャル）」を意味するので、
絶対値として

E_{\rm bind} \sim \frac{3}{5}\frac{GM^2}{R}

となる。

※ 中性子星は一様密度ではないが、オーダー見積もりとしてこの式がよく使われる。

理由は：

• 中心温度 $T \sim 10^{11}$ K
• 電子捕獲・逆β崩壊・pair annihilation (対消滅) が激しく発生
• 光は完全に閉じ込められるが、ニュートリノは外へ逃げられる

つまり、中性子星が“束縛される”ことで失うエネルギー = ニュートリノの総放射エネルギーになります。

実際、SN1987A の観測されたニュートリノエネルギーは

E_\nu \sim 3\times10^{53}\ {\rm erg}

であり、まさに中性子星の束縛エネルギーと一致します。

3. 親星の初期質量・イジェクタ・残った中性子星の関係

Cas A の場合：

量	値の典型	注釈
初期質量 $M_{\rm ZAMS}$	15–25 M☉	IIb SN progenitor
イジェクタ質量 $M_{\rm ej}$	2–4 M☉	X線・IR 観測から推定
中性子星の重力質量 $M_{\rm grav}$	1.3–1.5 M☉	観測値より
中性子星のバリオン質量 $M_{\rm baryon}$	1.5–1.7 M☉	EOS-dependent (状態方程式依存)
差分（束縛エネルギー）	0.1–0.3 M☉	→ ニュートリノへ

超新星で重要なのは、“中心に最終的に落ち込むバリオン質量”がニュートリノエネルギーを決める（初期質量やイジェクタ質量はほぼ関係ない）という点です。

4. 宇宙線電子は Cas A から地球に届くのか？

ほぼ届きません。理由は 2 つあります。

(1) ラーモア半径が小さすぎる（直進不可）

銀河磁場 $B\sim3\ \mu{\rm G}$ とすると電子のラーモア半径は：

エネルギー	ラーモア半径
10 GeV	~0.7 AU
1 TeV	~70 AU
1 PeV	0.3 pc

距離 3.4 kpc に対して 1千万分の1〜1万分の1 しかなく、
完全に拡散輸送になります。

荷電粒子のラーモア半径は

r_{\rm L} = \frac{p c}{e B}

で、ここで

• $p$：運動量
• $e$：電荷（電子の電荷の大きさ）
• $B$：磁場の強さ

です。

（補足）電子のラーモア半径を実際に計算してみる

0. 横方向の運動と「横質量」の歴史的導出と特殊相対論的円運動

磁場中で荷電粒子が円運動するとき、
ローレンツ力は常に速度に垂直である：

\mathbf{F} \perp \mathbf{v}.

このとき重要なのは 速度の大きさ $|\mathbf{v}|$ が時間微分で不変 な点である。

実際、

\frac{d}{dt}(v^2)
= \frac{d}{dt}(\mathbf{v}\cdot\mathbf{v})
= 2~\mathbf{v}\cdot\frac{d\mathbf{v}}{dt}

であり、右辺の $\mathbf{v}\cdot \frac{d\mathbf{v}}{dt}=0$（加速度が速度に垂直）より

\frac{d}{dt}(v^2)=0
\quad\Rightarrow\quad
v=\mathrm{const.}

となる。
したがって $\gamma = (1-v^2/c^2)^{-1/2}$ も一定で、時間微分の外に出せる。
（ここが歴史的導出の「肝」）

これにより 4-運動量の変化は

\frac{d}{dt}(\gamma m\mathbf{v})
= \gamma m \frac{d\mathbf{v}}{dt},

と簡単化され、横方向の慣性として

m_\perp = \gamma m

が定義される（古典的教科書での“横質量”）。

すると円運動の力学は

F_\perp
= m_\perp \frac{v^2}{r}
= \gamma m \frac{v^2}{r}.

一方、ローレンツ力の大きさは

F_\perp = e v B.

両者を等置すると

\gamma m \frac{v^2}{r} = e v B
\quad\Rightarrow\quad
r = \frac{\gamma m v}{eB}.

ここで運動量 $p=\gamma m v$ を用いれば、

\boxed{
r = \frac{p}{eB}
}

が得られる。

相対論的ラーモア半径の最終形

歴史的導出は上記のように「横質量」を使うが、
最終結果は現代相対論の運動方程式

\frac{d\mathbf{p}}{dt} = e \mathbf{v}\times\mathbf{B}

から直接得られる：

r_{\rm L} = \frac{p}{eB}.

超相対論的粒子では $p \simeq E/c$ なので、

\boxed{
r_{\rm L} \simeq \frac{E}{e B c}
}

となり、宇宙線加速・伝搬の議論ではこの式が標準的に使われる。

（例：10 GeV 電子、$B=3~\mu\mathrm{G}$ なら $r_{\rm L}\sim 0.7~\mathrm{AU}$）

1. まず「超相対論的」かを確認する

電子の静止質量エネルギーは

m_e c^2 = 511~{\rm keV} = 5.11\times 10^5~{\rm eV}.

例えば 10 GeV 電子なら

E = 10~{\rm GeV} = 10^{10}~{\rm eV}.

なのでローレンツ因子は

\gamma = \frac{E}{m_e c^2} 
       \approx \frac{10^{10}}{5\times 10^5}
       \sim 2\times 10^4.

$\gamma \gg 1$ なので 完全に超相対論的。
したがって運動量は

p \simeq \frac{E}{c}

という近似でよいです（厳密には $p = \gamma m_e v \simeq \gamma m_e c) ですが、(\gamma m_e c^2 = E$ なので同じ）。

2. 10 GeV 電子のラーモア半径を cgs 単位で計算

単位を全部 cgs にそろえます。

• 電子のエネルギー

  E = 10~{\rm GeV}
  = 10^{10}~{\rm eV}
  \simeq 1.6\times 10^{-12}\ {\rm erg/eV} \times 10^{10}
  \simeq 1.6\times 10^{-2}\ {\rm erg}.
  

• したがって運動量

  p \simeq \frac{E}{c}
  \simeq \frac{1.6\times 10^{-2}}{3\times 10^{10}}
  \simeq 5.3\times 10^{-13}\ {\rm g~cm/s}.
  

• 電子電荷（cgs・esu）：

  e \simeq 4.8\times 10^{-10}\ {\rm statC}.
  

• 銀河磁場：

  B = 3~\mu{\rm G} = 3\times 10^{-6}~{\rm G}.
  

これをラーモア半径の式に代入すると

r_{\rm L} 
  = \frac{p c}{e B}
  \simeq \frac{(5.3\times 10^{-13})(3\times 10^{10})}
               {(4.8\times 10^{-10})(3\times 10^{-6})}.

上の段（分子）：

(5.3\times 10^{-13})(3\times 10^{10})
\simeq 1.6\times 10^{-2}\ {\rm g\,cm^2/s^2}
= 1.6\times 10^{-2}\ {\rm erg}.

下の段（分母）：

(4.8\times 10^{-10})(3\times 10^{-6})
\simeq 1.44\times 10^{-15}.

したがって

r_{\rm L}
 \simeq \frac{1.6\times 10^{-2}}{1.44\times 10^{-15}}
 \sim 1.1\times 10^{13}\ {\rm cm}.

1 天文単位（AU）は

1~{\rm AU} \simeq 1.5\times 10^{13}\ {\rm cm}

なので、

r_{\rm L}(10~{\rm GeV}, 3~\mu{\rm G})
 \sim 0.7~{\rm AU}

というオーダーが出ます（表の値と一致）。

3. 一般式と 1 TeV / 1 PeV への拡張

10 GeV を基準にすると、ラーモア半径は

r_{\rm L} \simeq 1.1\times 10^{13}\ {\rm cm}\ 
\left(\frac{E}{10~{\rm GeV}}\right)
\left(\frac{B}{3~\mu{\rm G}}\right)^{-1}.

• 1 TeV = 100 × 10 GeV のとき：

  r_{\rm L} \sim 100\times 1.1\times 10^{13} {\rm cm}
             \sim 1.1\times 10^{15}{\rm cm}.
  

  1 AU = 1.5×10¹³ cm なので、

  r_{\rm L} \sim \frac{1.1\times 10^{15}}{1.5\times 10^{13}} {\rm AU}
             \sim 70~{\rm AU}.
  

• 1 PeV = 10⁵ × 10 GeV のとき：

  r_{\rm L} \sim 10^{5}\times 1.1\times 10^{13}{\rm cm}
             \sim 1.1\times 10^{18}{\rm cm}.
  

  1 pc ≃ 3.1×10¹⁸ cm なので、

  r_{\rm L} \sim \frac{1.1\times 10^{18}}{3.1\times 10^{18}} {\rm pc}
             \sim 0.3~{\rm pc}.
  

これで、表の

エネルギー	ラーモア半径
10 GeV	~0.7 AU
1 TeV	~70 AU
1 PeV	~0.3 pc

が、511 keV の電子質量を踏まえた相対論的近似 $p\simeq E/c$ からちゃんと出ていることが確認できます。

(2) 冷却時間が短すぎる（そもそも生き残れない）

電子は、磁場でシンクロトロン放射を、光子場（CMB など）で逆コンプトン (IC) 放射を起こすため、超相対論的電子は非常に速く冷えるという特徴があります。

● 冷却時間の前提（超相対論的電子）

超相対論的電子 $γ ≫ 1$のエネルギー損失率は

-\frac{dE}{dt}
= \frac{4}{3}\sigma_{\rm T} c~ \gamma^2 \left(U_B + U_{\rm rad}\right),

ここで

• $\sigma_{\rm T}$：トムソン断面積
• $U_B = B^2/(8\pi)$：磁場エネルギー密度
• $U_{\rm rad}$：CMB などの放射場のエネルギー密度

です。

冷却時間は

t_{\rm loss}
= \frac{E}{|dE/dt|}
\simeq \frac{3 m_e c}{4 \sigma_{\rm T} \gamma (U_B+U_{\rm rad})}.

シンクロトロン・IC の冷却は γ に反比例し、エネルギーに反比例します。

● 数値化するとこうなる

銀河磁場 $B \sim 3~\mu\mathrm{G}$、CMB のエネルギー密度 $U_{\rm CMB} \simeq 0.26~ \mathrm{eV/cm^3}$ を入れた代表式はよく

t_{\rm loss} \approx 3\times10^{5}\ {\rm yr}\ \left(\frac{E}{1~{\rm GeV}}\right)^{-1}

と書かれます。

これから、

• 1 TeV 電子：

  t_{\rm loss} \sim 300~\mathrm{yr}
  

• 10 TeV 電子：

  t_{\rm loss} \sim 30~\mathrm{yr}
  

Cas A の年齢は 約 340 年。
→ 1 TeV 以上の電子は、加速されてもほぼ全部冷え切ってしまう。

⇒ 結論：電子はそもそも“生き残れない”ので地球に届くはずがない

距離 3.4 kpc の問題以前に、電子は 1,000 年単位の旅をするだけの寿命がない。

では陽子はどうか？（質量依存性に着目）

ここも重要で、シンクロトロン損失は電子のような軽い粒子にだけ強く効くという点です。

シンクロトロン損失の一般式は

P_{\rm syn} \propto \frac{q^4}{m^2} \gamma^2 B^2.

電荷 $q$ は同じでも、質量 $m$ が違うと損失率が大きく変わる。

陽子は電子の

\frac{m_p}{m_e} \approx 1836

倍重いので、

\frac{P_{\rm syn,p}}{P_{\rm syn,e}}
\simeq \left(\frac{m_e}{m_p}\right)^2
\approx (1/1836)^2
\sim 3\times 10^{-7}.

つまり、

● 陽子の放射冷却は電子の 1,000,000 分の 1 以下。

したがって、陽子の冷却時間は

t_{\rm loss,p} \sim 10^6~ t_{\rm loss,e}

となり、宇宙線陽子は銀河内で“ほぼ冷えない”。

⇒ 陽子は 10⁵〜10⁷ 年スケールでエネルギーを保持できる

（だから銀河宇宙線の主成分であり、電子よりはるかに長距離を伝搬できる）

まとめ（冷却の観点）

粒子	冷却機構	冷却時間	SNR → 地球に届くか？
電子	synch + IC	E⁻¹ で急速に短くなる。TeV では 10–300 年	届かない（寿命が短すぎる）
陽子	radiative loss は電子の 10⁻⁶〜10⁻⁷	10⁵–10⁷ 年（ほぼ冷えない）	到達可能だが、340 年の Cas A からはまだ無理

5. 宇宙線陽子・イオンなら届くのか？

電子と違って、

• Hadronic CR（陽子・イオン）は放射冷却が極めて弱い
• TeV〜PeV 領域でも数百年〜数十万年スケールでエネルギーを保てる

ため、将来的には銀河宇宙線の一部として地球に届く可能性は十分にある。

しかし、Cas A までの距離は 3.4 kpc、年齢は 340 年。

この若さでは、宇宙線が拡散で到達できる距離は

\ell_{\rm diff} \sim \sqrt{2 D(E) t}

となり、TeV〜PeV でも 10 〜 1000 pc 程度にとどまる。
（後述したB/C 比から推定される銀河平均の拡散係数に基づく）

→ 3.4 kpc には全く届かない。

（補足）陽子のラーモア半径は陽子電子質量比「約2000倍」でスケール

ラーモア半径は

r_{\rm L} = \frac{p}{eB}

なので、同じ速度・同じガンマ因子なら 運動量 $p=\gamma m v$ に比例する。

陽子質量は

m_p \simeq 938\ {\rm MeV}/c^2,
\quad
m_e \simeq 0.511\ {\rm MeV}/c^2,

したがって

\frac{m_p}{m_e} \approx 1836.

つまり、同じローレンツ因子 $\gamma$ の電子と陽子を比べると陽子のラーモア半径は電子の 約 1800 倍 になる。

この質量スケールを使えば、電子で得た値を陽子に読み替えることができる：

• 10 TeV 電子 → $r_L \sim 70\ \mathrm{AU}$
• 10 TeV 陽子 → $r_L \sim 70 \times 1800 \sim 1.3\times10^5\ \mathrm{AU} \approx 0.6\ \mathrm{pc}$

同様に PeV では：

• 1 PeV 陽子 → $r_L \sim 1\ \mathrm{pc}$
• 10 PeV 陽子 → $r_L \sim 10\ \mathrm{pc}$

これで、先に示した値
（1 PeV → 1 pc、10 PeV → 10 pc）
が自然に理解できる。

（結論）

陽子は電子より遥かに大きなラーモア半径を持つが、それでも

• 銀河磁場に対しては依然として“曲げられる”領域
• 直進して方向情報を保持するには r_L ≪ 3.4 kpc

であり、Cas A から地球へ「ビーム状に」届くことはない。

来るとしても 完全に拡散輸送され、“どこから来たか分からない銀河宇宙線の一部”として将来的に寄与するだけである。

6. では、カシオペアAの宇宙線はどう観測する？

地球には来ないが、SNR 内部では

• 電子：シンクロトロン (radio–X)
• 電子：IC(CMB) → TeV γ線
• 陽子：p–p 衝突 → π⁰ 崩壊 → GeV–TeV γ線

で電磁放射を出します。

Fermi-LAT や MAGIC は Cas A の γ線スペクトルから Hadronic 成分を示唆しています。

つまり、宇宙線そのものは地球に届かないが、加速現場の“放射シグナル”は観測できる

7. 最後に：Cas A と「超新星と宇宙線」の全体像

まとめると：

• 銀河の宇宙線パワー（10⁴⁰–10⁴¹ erg/s）は「100年に一度、10⁵¹erg の超新星」で丁度よい
• 超新星のエネルギーの 99% は中性子星の重力束縛エネルギーとしてニュートリノで放出される
• 親星の初期質量とイジェクタ質量は「外層の問題」。
• ニュートリノエネルギーを決めるのは最終的な 中性子星の質量と半径
• 宇宙線電子は距離・冷却の両面で ほぼ確実に地球に届かない
• 宇宙線陽子は将来的には届くが、Cas A は若すぎてまだ到達していない
• しかし Cas A は X線・γ線で「宇宙線加速現場」を明瞭に見せてくれる
  → SNR が Galactic CR を供給するという標準モデルの最重要例

おわりに

Cas A は、

• ニュートリノ主導爆発エネルギー
• 中性子星の質量–半径–EOS
• 宇宙線加速と放射
• 銀河宇宙線のエネルギー収支
• 宇宙線の伝搬（拡散・冷却）

といった様々なテーマを一度に学べる、天体物理学にとって絶好の天体です。

◆ 補足1：宇宙線の拡散係数はどう推定される？

宇宙線の拡散係数 $D(E)$ は、主に 銀河内での二次生成元素の比（B/C 比） から推定されます。

宇宙線が銀河磁場を拡散しながら伝搬する間に、一次 cosmic-ray（C, N, O など）が ISM と衝突して生成される二次 cosmic-ray（B, Li, Be など）が作られます。

● もし拡散係数が小さい（＝宇宙線が長く銀河に閉じ込められる）

→ 二次生成が増える → B/C 比が高くなる

● 逆に拡散係数が大きい（＝すぐ銀河外へ逃げる）

→ 二次生成が起きにくい → B/C 比が低くなる

この B/C 比（一次 CR：C, O など／二次 CR：B, Be, Li）のエネルギー依存性をAMS-02、ACE/CRIS、Voyager などの観測データと、銀河宇宙線伝搬モデル（拡散方程式）
（例：Strong & Moskalenko 1998, Maurin et al. 2001）を比較してフィットすると、銀河平均の拡散係数としては

D(E) \sim 10^{28}
\left(\frac{E}{3~\mathrm{GeV}}\right)^{\delta}
\ {\rm cm^{2}/s},
\quad 
\delta \sim 0.3\text{–}0.6

という経験則が得られることが知られています。

● 重要な点

この推定は「銀河平均」の値であり、SNR 周辺や spiral-arm 内部では局所的に 1–2 桁小さい（turbulence が強い） という議論もあります。

Cas A のような若い SNR 周辺は、通常より拡散が遅い（D が小さい）可能性が高く、電子・陽子が“長く閉じ込められる”ことで、シンクロトロン・π⁰ 崩壊 γ線が強く見えると解釈されています。

◆ 補足2：なぜニュートリノは外に逃げられるのか？

超新星中心部は、

\rho \sim 10^{14-15}\ {\rm g/cm^3}

という核物質密度で、光子は完全にトラップされますが、
ニュートリノだけは弱相互作用によって透過できるという特徴があります。

● 断面積の違い

• 電磁相互作用の断面積：巨大 → 即トラップ
• ニュートリノの弱相互作用の断面積：

  \sigma_\nu \sim 10^{-42}\ \mathrm{cm^2}
  

  → 10⁵ km の物質を通してもほぼ進める小ささ

● “逃げられる”と言っても自由に飛べるわけではない

中心の密度が高いため、ニュートリノの mean free path（自由行程）は数 m〜数十 m 程度に縮むことがあります。それでも光子の "10⁻¹² cm レベルの自由行程" と比べて 圧倒的に長い。

その結果：

• コア中心 → ニュートリノは拡散しながら脱出（diffusive escape）
• 外層では密度が急低下 → ほとんど相互作用せず自由に外へ飛ぶ

これにより、コア崩壊の熱と束縛エネルギーが 10 秒スケールでニュートリノとして放出される。

◆ 補足3：“超新星爆発”のどのタイミングでニュートリノが出て、ejecta が出るのか？

タイムラインでまとめると分かりやすいです。

★ Step 1：Core collapse（ミリ秒〜数 ms）

Feコアが電子捕獲を起こし、重力崩壊。この時点で 大量の電子ニュートリノが burst（neutronization burst） として出る。

★ Step 2：Bounce（10 ms スケール）

コアが核密度に達すると反発して "bounce"。
衝撃波（shock wave）が誕生するが、数 ms で失速する（stalled shock）。

★ Step 3：Neutrino heating（10–300 ms）

proto-NS からの ニュートリノ加熱によって衝撃波が“再活性化”される（neutrino-driven explosion）。

この時点が 実質的な「爆発の開始」 に対応します。

• ニュートリノは常に中心から供給
• ejecta は shock の復活後に外層を押し出す

★ Step 4：Shock breakout（数秒〜）

衝撃波が外層に伝わり、stellar envelope を吹き飛ばす。
ここで ejecta（外層）が一気に外へ噴出する。

★ Step 5：Neutrino cooling phase（~10 秒）

proto-NS が冷えて縮む過程で、中性子星の重力束縛エネルギー（2–4×10⁵³ erg）が
10秒間でニュートリノとして放射される。

つまり：

• ニュートリノ放射は collapse の瞬間から始まり、爆発前にピーク
• ejecta は shock が revived したあと、数秒〜数時間で外に放射状に展開

という構造になっています。

◆ 補足4：なぜ Cas A は RSG 〜 WR の progenitor だと考えられているのか？

Cas A の progenitor（親星）が Red Supergiant (RSG) か Wolf-Rayet (WR) と考えられる理由はいくつもあります。

(1) H 層がほとんど失われていた（Type IIb）

Cas A の光回折エコー観測から、

「爆発時には水素殻が極めて薄い Type IIb だった」

ことがわかっています（Krause et al. 2008）。

Type IIb は、

• 大部分の H envelope を失っている
• RSG → WR → SN という evolution の典型例

です。

(2) CNO 処理された物質が ejecta に見られる

X線で観測される N/O 比の異常は、“CNO サイクルで形成された外層” があらかじめ取り除かれていたことを示します。RSG や WR 風による強い質量損失があれば自然に説明できます。

(3) 周囲の QSF (quasi-stationary flocculi) が dense-RSG 風に一致

Cas A 周囲の “QSF” と呼ばれる凝縮塊は、

• 高密度
• 低速 (100 km/s 程度)
• 重元素の組成が RSG 風と整合

という特徴を持つため、爆発前にRSG風があったという有力証拠になっています。

(4) 初期質量（15–25 M☉）の理論的進化経路が RSG → WR を示唆

この質量帯の星は強い質量損失を受けると：

• RSG で外層を失う
• さらに stripping が進むと WR になる
• 最終的に水素殻をほぼ失って IIb を起こす

Cas A の観測と整合的です。

(5) ejecta の質量（~3 M☉）が「薄い envelope」を示唆

大きな RSG envelope を残したまま爆発すれば ejecta は 10–15 M☉ になりますが、Cas A では明らかに小さい。→ 外層がすでに剥がされていた（RSG/WR 風 or binary mass transfer） と判断されます。

Red Supergiant（RSG）と Wolf–Rayet（WR）とは？

Red Supergiant（RSG）と Wolf–Rayet（WR）は、どちらも大質量星の進化段階ですが、
物理状態・大きさ・表面温度・風の性質・質量損失率がまったく異なる“対照的な天体”です。以下に、本質的な違いを整理します。

◆ RSG（赤色超巨星）とは？

● 主な特徴

• 半径がとてつもなく大きい（数百〜1000 R☉）
• 表面温度が低い（T ≈ 3000–4000 K） → 赤く見える
• 水素 envelope をまだ保っている段階
• 対流が非常に発達した巨大な envelope
• 質量損失はあるが、速度は遅い風（10–20 km/s）

● 典型例

• ベテルギウス
• アンタレス

● 質量損失率

\dot{M}_{\rm RSG} \sim 10^{-6} \sim 10^{-4} \ M_\odot/{\rm yr}

遅いが“量は多い”風。
星の周囲に密度の高い envelope を作る。

◆ WR（Wolf–Rayet 星）とは？

● 主な特徴

• 極端な高温（T ≈ 50,000–150,000 K）
• 半径が非常に小さい（~1–10 R☉）
• 外層の水素がほぼ失われ、
  He, C, N, O がむき出しの“大量に剥ぎ取られたコア”
• 核融合生成物が表面に露出している

● 質量損失率

\dot{M}_{\rm WR} \sim 10^{-5} \sim 10^{-4} \ M_\odot/{\rm yr}

風速度は 1000–3000 km/s の高速風。
→ 周囲の星間物質（ISM）を大きく吹き飛ばし、広大な bubble を形成する。

● なぜ外層がない？

• 自身の強風で envelope を失った
• または binary mass transfer で剥ぎ取られた

◆ RSG → WR の進化（大質量星の典型ルート）

初期質量が 15–25 M☉ 程度の星は、

1. 主系列（H 燃焼）
2. RSG（He 核＋厚い H envelope）
   　→ 強い質量損失
3. WR（H envelope がほぼ失われ、He や CNO が露出）
4. 超新星爆発（IIb / Ib / Ic）

という進化をたどる可能性が高い。

• H envelope が多少残れば IIb
• H が完全に取れれば Ib
• He も失えば Ic

Cas A が Type IIb だったのは“H をほとんど失ったが、わずかに残っていた”という状況を示します。

→ これは RSG または WR を経た星で自然に説明できる。

◆ RSG と WR の違いをまとめると

項目	RSG	WR
半径	超巨大（100–1000 R☉）	小さい（1–10 R☉）
温度	低温（3000–4000 K）	高温（50,000 K 以上）
外層	厚い H envelope	水素ほぼ無し、He/C/N/O が露出
風	遅いが密度が高い	高速・低密度だが質量損失率が大きい
スペクトル	分子吸収線が多い	強い広がった発光線（HeII, CIII, NIVなど）
典型的 SN 型	II-P, II-L	IIb, Ib, Ic

(6) エネルギーの“どれが何なのか”をもう一度明確化

超新星（特にコア崩壊型）の場合

1. ニュートリノが持ち去るエネルギー

   E_\nu \sim (2–4)\times10^{53}~{\rm erg}
   

   → 中性子星の重力束縛エネルギー（＝コア崩壊エネルギー）

2. ejecta の運動エネルギー（爆発エネルギー）

   E_{\rm SN} \sim 10^{51}~{\rm erg}
   

   → shock revival の“残りカス”で、ニュートリノの ~1％未満

3. 光学的輻射（SN light curve）が放射するエネルギー

   E_{\rm rad} \sim 10^{49}\ {\rm erg}
   

   → さらにその 1/100

のような“完全に階層構造の違う 3 つのエネルギー”が存在します。

◆ 図式にするとさらに明確

(1) コア崩壊で解放された重力束縛エネルギー 〜 3×10^53 erg
         ↓ 99%以上
       ニュートリノがすべて持ち去る
         ↓ 1%未満
(2) 爆発の運動エネルギー（ejecta）〜 10^51 erg
         ↓ 1%
(3) 光度（電磁波）として出るエネルギー 〜 10^49 erg

つまり、宇宙線の議論で登場する $10^{51}$ erg は「爆発の力」、つまり ejecta の 運動エネルギーのことです。この量は、

• 外層が秒速数千 km で吹き飛ぶ
• Cas A の shock radius が現在 2.5 pc まで広がっている
• Sedov–Taylor 期に入る前の自由膨張の動力源

などを説明するための “爆発エネルギー（kinetic energy of ejecta）” です。
これはニュートリノとは別筋の話です。

◆ では、ニュートリノはどれくらいのエネルギーか？

中性子星の重力束縛エネルギー：

E_\nu \approx 0.1\text{–}0.3\ M_\odot c^2
\approx (2\text{–}5)\times10^{53}{\rm erg}

SN1987A の観測からも

E_\nu \approx 3\times10^{53}{\rm erg}

と推定され、この値は 爆発エネルギー(ejecta が持つ運動エネルギー)の ~300 倍 です。

◆ まとめ

物理量	大きさ	意味
$E_\nu$	~10⁵³ erg	中性子星形成の重力束縛エネルギー（ニュートリノが持ち去る）
$E_{\rm SN}$	~10⁵¹ erg	ejecta の運動エネルギー（爆発の力）
放射	~10⁴⁹ erg	SN(超新星) の光として出る

→ 10⁵¹ erg はニュートリノではなく、爆発の運動エネルギーのこと。

X線観測やTESの開発者向けのまとめ

最後に、X線観測やTES開発に関わる人向けを念頭に、もう少し広い視点で超新星残骸研究を整理してみます。

X線観測で粒子加速を調べることの位置づけ

1. 研究目的

• 超新星残骸の X線シンクロトロン放射を通じて、どこで・どのくらいの効率で電子が加速されているかを定量的に知ること。
• ショック速度・磁場強度・シンクロトロン冷却のバランスを測ることで、「SNR が銀河宇宙線の主要な加速源になりうるか？」 を検証すること。
• Cas A のような若い SNR を例に、爆発エネルギーのどのくらいが粒子加速に回っているかを、エネルギー収支の文脈で理解すること。

2. 得られた成果とその確実性

• X線シンクロトロンの急峻なフィラメント構造や時間変動から、少なくとも数十 TeV 級の電子が衝撃波近傍で加速されている ことはかなり確実であろう。
• 電子の最大エネルギーと冷却時間から、星間磁場(B ~ μG クラス)に比して、磁場増幅（B ~ 数百 μG クラス）が起きている ことも強く示唆されており、衝撃波での効率的な Fermi 加速の場としての妥当性が支持されている。
• γ線（Fermi-LAT, MAGIC など）との組み合わせにより、少なくとも一部の SNR（Cas A を含む）で Hadronic 成分（陽子・イオン）の加速が起きている可能性も高いと考えられている。

3. 研究の弱点・限界

• X線は基本的に 電子シンクロトロンに敏感であり、陽子・イオンのスペクトルや最大エネルギーを直接は測れない（Hadronic 成分は γ線側でしか見えない）。
• スペクトルだけを見ると、磁場強度と電子スペクトル形状の縮退（degeneracy） が残りやすく、「磁場 B が強いのか、電子が少ないのか」の切り分けが難しい。
• 3次元構造を 2次元投影で見ているため、幾何学的な投影効果・混合成分（thermal vs non-thermal）の分離にも不確実性が残る。
• 拡散係数や乱流スペクトルなど、宇宙線伝搬のパラメータは依然としてモデル依存 であり、「Cas A で測った値＝銀河平均」とは限らない。

4. 将来的な方向性

• 高分解能・高感度の X線分光像（XRISM の次世代、Athena など）により、熱成分と非熱成分を空間的・スペクトル的により精密に分離し、衝撃波構造と粒子加速効率の対応を明確にする。
• X線だけでなく、MeV–GeV–TeV γ線・電波・光赤外とのマルチ波長解析で、電子・陽子・環境（CSE/ISM）のパラメータを総合的に俯瞰する方向が重要。
• 理論側では、MHD + PIC シミュレーションと実観測（Cas A のような若い SNR）の空間構造・時間変動を 1:1 で比較できるレベルまでモデルを高度化し、「どのタイプの SNR が、どのエネルギー帯の銀河宇宙線に主に寄与しているか」を統計的に評価することが今後の大きな目標。
• 特に Cas A のような「若くて明るい残骸」は、時間変動（数年スケール）を追える貴重な“実験室” なので、継続モニタリングと次世代X線ミッションの組み合わせが鍵になる。

おまけ：X線TES・可視光TESによる「粒子加速現場の超精密観測」へ

将来的には、超伝導遷移端検出器（TES）による高分解能面分光も、粒子加速研究と相性のよいツールになっていきます。

• X線TESマイクロカロリメータ
  XRISM/Resolve に代表されるように、数 eV 以下のエネルギー分解能で

  • ショック近傍のドップラーシフト・ライン幅から「乱流・バルク流」を直接測る
  • イオン温度と電子温度の非平衡や、元素ごとの加熱効率を見分ける
    といった解析が可能になります。これは、「加速されていない側」のプラズマの状態を精密に測ることで、粒子加速に使われたエネルギーの“残り”を逆算するアプローチとも言えます。

• 可視光TESによる高精度フォトン計測
  可視〜近赤外領域でも TES を用いれば、

  • 1 光子ごとのエネルギーと到来時刻を同時に記録
  • ショックに対応するラインプロファイルのテイル成分や微細構造
  • 高時間分解能の変動（パルサー・ショック相互作用など）
    を同時に追いかけることができます。

こうした TES 系の計測は、

• 「加速された粒子（宇宙線）」そのものではなく
• 「加速を支えているプラズマと放射場の状態」を、エネルギー・速度・時間の 3 つの軸で高精度にマッピングする

という意味で、SNR における粒子加速の“環境側の制約”を劇的に強くするポテンシャルを持っています。Cas A のような若い超新星残骸を、TES でどこまで「解像」できるかは、今後の楽しみなテーマのひとつだと思います。

(おまけのおまけ) “若いSNRにTESは役立つか？”という指摘に対して

最後に、X線・可視光 TES（Transition-Edge Sensor）について、専門家からよく指摘されるポイントと、それに対する技術的な回答をまとめておきます。TES に期待されている将来性と限界を、公正な形で整理するための“おまけ”です。

1. 若い SNR の線幅は広すぎて TES が活かせないのでは？

Cas A のような若い SNR では、v ≃ 3000–6000 km/s に達する Doppler broadening により、Fe-K などの線幅が 数十〜数百 eV まで広がります。

そのため、

「TES の eV 級分解能は無意味では？」

という指摘があります。これは半分くらいは正しいツッコミです。以下の 3 点を考えると、TES の強みは別のところにあります。

（A）TES 分光の価値は「線を細く見る」ことではなく、「線の広がりの原因を空間的に分解できる」こと

若い SNR の線が広い理由は、単一の速度場ではなく、

• ショック前後の bulk 流
• ejecta knots
• reverse shock の異なる ion populations
• turbulence
  など、複数の速度成分が空間的に重なっているためです。

TES（位置分解型 TES, or IFU 的アプローチ）の強みは、空間的に細かく切り分けて、混ざった速度コンポーネントを分離できることです。(IFU＝Integral Field Unit：天文分野で「面分光（空間×スペクトルを同時取得）」を意味する計測手法）

分散型分光では slit に入った複数成分が混ざって1 本の“太い線” になるのに対し、TES は

• shock rim
• Fe-rich knots
• turbulence-dominated patches

といった局所構造を別々のスペクトルとして取れます。

（B）線幅が広くても、line 比 や イオン温度には TES 分解能が効く

広い線でも、相対位置が近い

• He-like triplet（r, i, f ライン）
• Kβ/Kα、衛星線（satellite lines）
• charge-state 分布

のような指標は eV 級の分解能で確度よく分離可能です。若い SNR はむしろ ionization nonequilibrium（NEI）やion–electron 非平衡が強く、TES 分光が効く典型領域です。(broadening の程度にも依存する話ですが)

（C）複雑な速度場を“成分ごと”に抽出できるのは TES の独自能力

• shock heating
• ion temperature の非平衡
• 乱流の空間分布
• 逆衝撃波での再加熱

などの情報は、空間分解能＋eV 分解能で初めて本格的に取得できるものです。

つまり、

「線が広いからTESが不要」ではなく、「広がりの原因を分離できるからTESが必要」

という捉え方も可能です。

2. 可視光TESは分散型分光器に勝てないのでは？

もう一つの典型的ツッコミ：

「可視光では分散型分光器（R ≈ 10⁴–10⁵）の方が圧勝。TES はスペック的に勝ち目がないのでは？」

これは 分散型と TES の用途が異なることを知らないと生じる誤解です。

（A）可視光 TES の強みは「分散分光の代替」ではなく「新しい計測モード」

可視光 TES は、

• 1 光子ごとのエネルギー（0.数eV 精度）
• 1 光子ごとの到来時刻（< ~μs）
• 極めて高い量子効率（QE）

を同時に実現します。

つまり TES が得意なのは、

• 高分解能スペクトル ではなく
• 高 S/N・高時間分解能・高 photon throughput の“1 photon spectroscopy”

です。

（B）分散型が苦手な領域で圧倒的に強い

例えば：

• 微弱な forbidden lines の時間変動を単光子で追う
• transient（重力波源・爆発天体など）の光子を漏れなくスペクトル化
• background が極端に低い領域で光子統計を最大化
• IFU 的な構造を作れば、空間×時間×スペクトルの 3D データを取得

といった「分散型が構造的に苦手な世界」で、TES は特徴が生きてきます。

（C）“勝ち筋”は分解能ではなく「取得できる情報量の次元」

TES のコア強みは、

• Space（空間）：pixel TES or IFU-like 配列で空間情報
• Energy（スペクトル）：単光子エネルギー測定
• Time（時間）：μs precision
• Flux（絶対光子数）：高 QE・無分散でロスが少ない

という 4D 情報の同時取得にあります。

それぞれの次元が“浅く広く”ではなく“深く高効率”で得られるため、分散分光とは 競合ではなく補完関係です。

TESの現実的なメリットのまとめ

• 若い SNR の線幅が広いことは事実だが、TES の価値は線を細く見ることだけではなく、混ざった速度成分を空間的に分離できること。

• 可視光 TES は分散型分光の競合ではなく、単光子エネルギー分光＋時間分解能＋高 QE という新しい観測パラダイムを提供する。

• TES 系計測は、粒子加速の“環境側”の物理（乱流・速度場・電離状態）を高精度で制約する手段であり、将来観測の重要な選択肢である。