本文へ移動
ISEGORIABenjamin Haire
English

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 万年になります。

参考文献: MIT OpenCourseWare · 8.901 Astrophysics I

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 と、同じ質量が中心へますます強く詰め込まれるからです。

参考文献: MIT OpenCourseWare · 8.901 Astrophysics I

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}\) です。

参考文献: MIT OpenCourseWare · 8.901 Astrophysics I

次に読む

線形代数:変換の幾何フーリエ解析:波で関数を作る力学系:安定性とカオス最適化:最良の選択の幾何確率:不確実性から学ぶ群論:対称性を代数にする代数的位相幾何:穴を見つける数値解析:計算が誤解を招くとき数論:整数のパターングラフ理論:経路・木・ネットワーク偏微分方程式:動く場情報理論:不確実性と符号変分法:経路と原理古典力学:運動と力電磁気学:場と誘導光学:光線・波・色熱力学:エネルギー・仕事・エントロピー量子力学:振幅とスピン複素解析:写像・留数・調和場流体力学:流れ・圧力・渦度統計と推測:データの中の信号特殊相対性理論:空間・時間・光微分幾何:曲率と形統計力学:ミクロ状態と温度測度論:大きさ・近似・収束マルコフ連鎖:遷移・定常性・吸収微分形式:循環・カール・引き戻し一般相対性理論:曲率・時計・光関数解析:ノルム・射影・作用素プラズマ物理:遮蔽・軌道・波リー群とリー代数:連続的な対称性ハミルトン力学:相空間とその幾何確率過程:ブラウン運動とノイズ固体物理:結晶の中の波制御理論:フィードバック・極・安定性論理と計算可能性:何が計算できるか双曲幾何:平行線が増える世界剛体力学:回転・転倒・歳差楕円曲線:足し算ができる幾何原子物理:軌道とスペクトル四元数:掛け算としての回転押し出しと引き戻し:写像を通した積分結び目理論:絡まりを見分けるフラクタル幾何:整数の間の次元量子情報:もつれとその限界宇宙論:膨張する宇宙セル・オートマトン:局所規則から生まれる計算複雑系:単純な多数の部品から生まれる秩序ガロア理論:方程式の対称性磁性とイジング模型:整列から生まれる秩序振動子と脱進機:時計はどうやって時を刻むか歯車と機構:運動を正確に伝える水晶振動子:時を刻む結晶スピーカー:ムービングコイル型ドライバーフィルターとクロスオーバー:音をドライバーに振り分ける室内音響:部屋はスピーカーの一部ウェーブレット:信号を拡大して見るレーザー物理:自らを複製する光表現論:行列として作用する群半導体物理:バンド・ドーピング・接合組合せ論:組み立てて数える音響学:弦、板、動く音源スイッチング電源:刻んで変換するゲーム理論:相手も選ぶときの戦略原子核物理:結合、崩壊、トンネル効果ピストンエンジン:クランクからサイクルまでエルゴード理論:時間平均が空間平均に等しくなるとき連続体力学:固体はどう伸び、曲がり、座屈するか電動機:電流をトルクに変える重力波:空間の幾何のさざ波ロボットアーム:運動学とヤコビ行列リーマン面:多価関数のすみか射影幾何:平行線が交わる場所構造力学:橋はどう荷重を支えるか

数学百科事典に戻る