ISEGORIA / 数学百科事典
恒星天体物理:星の輝きと死
星は、自分の重さに対して圧力で支えられたガスの球で、星のほとんどの性質はこのつり合いから決まります。明るさを表面温度に対してプロットすると、星は斜めの帯に集まり、そこでは質量だけで色・大きさ・光度・寿命が決まります。内部では一本の方程式が中心から表面までの密度・圧力・温度の変化を決めます。燃料が尽きると量子力学が主役になります。これ以上詰められない電子の球が白色矮星を支え、太陽の約 1.4 倍の質量を超えるとそれすら支えきれなくなります。
前提知識: 熱力学と原子核物理
予想してから操作し、各実験の練習問題で理解を確かめてください。練習問題は三つの難易度で毎回新しく生成されます。グラフは数式の図示であり、証明ではありません。
1. ヘルツシュプルング・ラッセル図
星の性質のほとんどは二つの数で表せます。放つ光の量(光度 L)と、表面の温度(有効温度 T)です。一方を他方に対してプロットすると、星はばらばらには散らばりません。大半は斜めの帯、主系列に乗ります。そこでは星は中心で水素を核融合していて、質量だけで他のすべてが決まります。重い星ほど高温で、明るく、大きく、寿命が短いのです。ステファン・ボルツマンの法則により、図のどの点にも半径が決まるので、巨星と白色矮星はすぐに区別できます。金色の点を主系列に沿ってドラッグすると、太陽と並んだ星の色と大きさが変わり、スペクトルはウィーンの法則に従って短い波長へ動きます。「任意の星」に切り替えると、点を自由に置いて、その点が要求する半径を読み取れます。
計算例. 太陽の 10 倍の質量の主系列星では \(L=1.4\,M^{3.5}=4400\,L_\odot\)、\(R=M^{0.6}=4.0\,R_\odot\) なので、\(T_{\rm eff}=5772\,{\rm K}\,(L/R^2)^{1/4}=23{,}600\) K の青白い B 型星で、スペクトルのピークは \(\lambda_{\max}=2.898\times10^{-3}/23600=123\) nm の紫外線にあります。寿命はおよそ \(10\ {\rm Gyr}\times10/4400=23\) Myr で、燃料は 10 倍あるのに太陽の約 400 分の 1 です。
注意点. ここに描いた質量光度関係と質量半径関係は経験的な近似です(四つの質量範囲で \(L\propto M^{2.3}\)、\(M^{4}\)、\(M^{3.5}\)、\(M\)、また \(R\) は太陽質量以下で \(M^{0.9}\)、以上で \(M^{0.6}\))。精度は数十 % で、実際の星は組成と年齢によってばらつき、年を取るにつれて出発点から右上へずれていきます。図中の既知の星は概算値です。色は黒体の色で、実際の恒星のスペクトルはそれにおおよそしか従いません。寿命は水素の 10% が 0.7% の効率で燃えるものとして計算しています。
太陽の 10 倍の質量の星は、燃料が 10 倍しかないのに、なぜ寿命が太陽の数百分の 1 なのですか?
燃料をはるかに速く燃やすからです。光度は太陽の約 4400 倍(\(L\propto M^{3.5}\))で、10 倍の燃料を 440 倍の速さで使い、寿命は約 1/440、およそ 2300 万年になります。
2. 流体静力学的平衡とポリトロープ
星はなぜつぶれもせず、吹き飛びもしないのでしょうか。どの半径でも、圧力の勾配がその外側のガスの重さを支えなければなりません。\(dP/dr=-Gm\rho/r^2\) です。圧力が密度のべきに従うとき、\(P=K\rho^{1+1/n}\)、構造の方程式は一本の無次元の方程式、レーン・エムデン方程式にまとまり、その解 θ(ξ) が内部全体、つまり密度・圧力・温度・内側の質量を決めます。指数 n はガスのやわらかさを表します。星を切り開いて、核が見えるようにしてあります。n が大きいほど同じ質量がより小さく密な中心に詰まり、中心温度はビリアル則 \(T_c\propto\mu M/R\) から決まります。収縮する星が熱くなる理由です。探針を赤道に沿ってドラッグして、その場の状態を読み取ってください。
計算例. 太陽を \(n=3\) のポリトロープとみなすと(\(M=1\,M_\odot\)、\(R=1\,R_\odot\)、\(\mu=0.6\))、\(\xi_1=6.897\)、\(\rho_c/\bar\rho=54.18\) なので \(\rho_c=54.18\times1.41=76\) g/cm³ です。中心圧力は \(P_c=11.05\,GM^2/R^4=1.2\times10^{16}\) Pa、中心温度は \(T_c=11.8\) MK です。質量の半分は \(r=0.283\,R\) の内側にあり、重力エネルギーは \(|\Omega|=1.5\,GM^2/R\) です。
注意点. ポリトロープは構造のモデルであって、構造の計算ではありません。実際の太陽の \(\rho_c\approx150\) g/cm³、\(T_c\approx15.7\) MK で、組成と不透明度が半径とともに変わるので、一つの n では表せません。ここではガスを一定の平均分子量 μ の理想気体とし(電離した太陽組成で 0.6、純水素で 0.5、純ヘリウムで 4/3)、\(T\propto P/\rho=\theta\) としています。放射圧と縮退は無視しています。n ≥ 5 では有限の半径で密度がゼロになる解がないので、n は 4.5 までにしています。
n を 0 から増やすと、中心密度と平均密度の比はどうなりますか?
限りなく増えます。n = 0(一様な球)で 1、n = 3/2 で 5.99、n = 3 で 54.2 と、同じ質量が中心へますます強く詰め込まれるからです。
3. 縮退圧とチャンドラセカール限界
白色矮星は燃料を使い果たし、熱以外のもので支えられている星です。電子がパウリの原理により、温度によらず、ますます高い運動量の状態へ押し込まれるほど密に詰められています。この縮退圧は密度だけで決まります。電子が遅いうちは \(P\propto\rho^{5/3}\) で、重い矮星ほど小さくなります。さらに押しつぶすと電子は光速に近づき、圧力は \(\rho^{4/3}\) へやわらかくなって、支えられる質量に上限が生じます。これがチャンドラセカール限界で、炭素と酸素なら太陽質量の約 1.46 倍、基本定数と電子一個あたりの核子数だけで決まります。点を質量半径曲線に沿ってドラッグすると、地球と並んだ矮星が縮み、密度分布が n = 3/2 のポリトロープから n = 3 のものへ移っていきます。
計算例. \(\mu_e=2\)、中心密度 \(10^6\) g/cm³ のとき、電子の運動量は \(x_c=p_F/m_ec=0.80\) で、相対論的になり始めた程度です。流体静力学的平衡を外側へ積分すると \(M=0.40\,M_\odot\)、\(R=1.72\,R_\oplus=11{,}000\) km になります。\(3\times10^7\) g/cm³ では質量が \(1.0\,M_\odot\) に達し、半径は \(0.88\,R_\oplus\) まで縮み、シリウス B とほぼ同じです。限界は炭素・酸素で \(M_{\rm Ch}=5.83/4=1.46\,M_\odot\)、鉄(\(\mu_e=2.15\))で \(1.26\,M_\odot\) です。
注意点. これは冷たく、自転せず、ニュートン重力で、電子が相互作用しない、チャンドラセカール(1935)の星です。実際の白色矮星は温かいのでわずかに大きく(観測されるシリウス B は曲線より約 5% 大きく、点は概算の観測値で描いています)、さらにクーロン引力、一般相対性理論、逆ベータ崩壊(10⁹ から 10¹⁰ g/cm³ 付近で効き始めます)により、曲線が M_Ch に届く前に実際の星は不安定になります。したがって曲線の右端は理想化した極限で、実際の星が到達する状態ではありません。純水素(μₑ = 1)は思考実験です。
白色矮星に質量を加えると、なぜ小さくなるのですか?
質量が増えると支える重さが増えるので、中心圧力が上がらなければなりません。縮退圧は密度でしか上がらない(\(P\propto\rho^{5/3}\))ので、星はより密に押しつぶされ、質量が大きいのに半径は小さくなります。\(R\propto M^{-1/3}\) です。