恆星的一生(進階):質量即命運的物理推導
從流體靜力平衡到赫羅圖,看懂恆星演化為何是一條可被計算的軌跡
如果恆星是一條河,質量就決定了它流向哪片海
你已經知道恆星「出生於分子雲、靠核融合發光、最後依質量走向白矮星、中子星或黑洞」這條主線。但這裡有個更深的問題值得我們追問:為什麼一顆恆星「一旦質量被決定」,它的一生幾乎就被注定了?我們無法回到太陽出生那一刻去調整它,但只要知道它的質量,天文學家就能相當有把握地告訴你:它會發多亮的光、活多久、何時膨脹成紅巨星、最後留下什麼殘骸。
這種「質量即命運」的鐵律,並不是觀測歸納出來的巧合,而是從幾條極簡單的物理原理——力的平衡、能量的傳輸、氣體的狀態——層層推導出來的必然。這篇進階文章不再重述各個階段「發生了什麼」,而是要帶你看見支配整條恆星生命的結構方程式(stellar structure equations) 與赫羅圖(Hertzsprung–Russell diagram, HR diagram) 這張「恆星的相空間地圖」,理解恆星演化為什麼是一條可被計算的軌跡,而不是一串隨機事件。

為什麼恆星不會塌縮也不會炸開:流體靜力平衡
恆星是一團被自身重力束縛的熾熱氣體。重力時時刻刻想把每一層氣體往中心壓,那為什麼它沒有在幾分鐘內塌成一個點?答案是流體靜力平衡(hydrostatic equilibrium):向內的重力,恰好被向外的壓力梯度抵消。
取恆星內部一個薄殼層,半徑 $r$、厚度 $dr$。它受到的向內重力與向外的壓力差必須平衡,於是得到恆星結構的第一條方程式:
$$\frac{dP}{dr} = -\frac{G\, M(r)\, \rho(r)}{r^2}$$
其中 $P$ 是壓力、$\rho$ 是密度、$M(r)$ 是半徑 $r$ 以內的累積質量、$G$ 是重力常數。負號代表壓力隨半徑增加而下降——核心壓力最高,表面趨近於零。
這條看似簡單的式子藏著驚人的力量。我們可以用它做一次「量級估算」:把微分粗略換成整顆恆星的尺度,$dP \sim P_c$(核心壓力)、$dr \sim R$(恆星半徑)、$M(r)\sim M$、$\rho \sim M/R^3$,得到
$$P_c \sim \frac{G M^2}{R^4}$$
對太陽代入數字($M_\odot \approx 2\times10^{30}\ \mathrm{kg}$、$R_\odot \approx 7\times10^{8}\ \mathrm{m}$),會算出核心壓力約 $10^{14}\ \mathrm{Pa}$ 的量級(精確的太陽模型給出約 $2.5\times10^{16}\ \mathrm{Pa}$,量級估算已捕捉到主要依賴關係)。重點不在精確值,而在那個 $M^2/R^4$ 的依賴——它告訴我們:質量越大、半徑越小的恆星,核心被壓得越緊、越熱,而核融合反應率對溫度極度敏感,這就埋下了「大質量恆星燒得又快又猛」的種子。
用一條方程式估出太陽核心有多熱
把流體靜力平衡再往前推一步,配合理想氣體狀態方程 $P = \dfrac{\rho k_B T}{\mu m_H}$($\mu$ 是平均分子量、$m_H$ 是氫原子質量、$k_B$ 是波茲曼常數),我們可以反過來估計核心溫度。將前面的 $P_c \sim GM^2/R^4$ 與 $\rho \sim M/R^3$ 代入:
$$k_B T_c \sim \frac{G M \mu m_H}{R}$$
這正是一個深刻的物理直覺的數學版本:重力把氣體壓縮做的功,轉化成了粒子的熱運動。 對太陽代入(取 $\mu \approx 0.6$,對應電離的氫氦混合氣體),可以估出 $T_c$ 約為一千萬度的量級,與標準太陽模型的 $1.5\times10^7\ \mathrm{K}$ 相符。
這個結果是恆星物理的轉捩點:唯有質量足夠大的天體,重力壓縮才能讓核心突破約 $10^7\ \mathrm{K}$,點燃氫核融合——這正是恆星與棕矮星(brown dwarf)、行星的物理分界線(約 $0.08\,M_\odot$)。質量決定核心溫度,核心溫度決定能否發光、發多亮,一切都從那條力平衡的式子流淌出來。
質光關係:為什麼亮一千倍的星,壽命卻短得多
恆星結構的四條方程式(流體靜力平衡、質量連續、能量產生、能量傳輸)聯立起來,會自然導出一條極重要的經驗與理論皆吻合的關係——質光關係(mass–luminosity relation)。對主序帶(main sequence)上的恆星,光度 $L$ 與質量 $M$ 大致呈
$$L \propto M^{\alpha}, \qquad \alpha \approx 3.5$$
指數 $\alpha$ 在不同質量區間會在 3 到 4 之間變動,但 3.5 是個好用的代表值。它的意思非常驚人:一顆質量是太陽 10 倍的恆星,光度不是太陽的 10 倍,而是約 $10^{3.5} \approx 3000$ 倍。
現在把這個結果和「燃料總量」放在一起。恆星可用的核燃料正比於質量 $M$,而它「燒錢的速度」正比於光度 $L$。因此主序壽命
$$t_{\text{MS}} \propto \frac{M}{L} \propto \frac{M}{M^{3.5}} = M^{-2.5}$$
這就是恆星生命中最反直覺、也最重要的一條定律:質量越大的恆星,活得越短。 一顆 $30\,M_\odot$ 的藍色巨星只能燃燒幾百萬年就壯烈死去;而一顆 $0.3\,M_\odot$ 的紅矮星(red dwarf)可以節儉地燃燒上兆年——遠遠超過目前宇宙 138 億年的年齡。換句話說,宇宙至今還沒有任何一顆紅矮星走完它的一生。
動手算一下:兩顆恆星,誰先謝幕?
讓我們把這套邏輯用在具體數字上。已知太陽的主序壽命約 $t_\odot \approx 10^{10}$ 年(一百億年)。請估計一顆 $15\,M_\odot$ 恆星與一顆 $0.5\,M_\odot$ 紅矮星的主序壽命。
我們用更精確的比例式,把光度也展開:
$$t_{\text{MS}} = t_\odot \cdot \frac{M/M_\odot}{L/L_\odot} = t_\odot \cdot \left(\frac{M}{M_\odot}\right)^{1-3.5} = t_\odot \cdot \left(\frac{M}{M_\odot}\right)^{-2.5}$$
$15\,M_\odot$ 的恆星:
$$t_{\text{MS}} = 10^{10} \times 15^{-2.5}\ \text{年}$$
計算 $15^{2.5} = 15^2 \times 15^{0.5} = 225 \times 3.87 \approx 871$,所以
$$t_{\text{MS}} \approx \frac{10^{10}}{871} \approx 1.1\times10^{7}\ \text{年}$$
約一千一百萬年——以宇宙尺度看不過是一瞬。
$0.5\,M_\odot$ 的紅矮星:
$$t_{\text{MS}} = 10^{10} \times 0.5^{-2.5}\ \text{年}$$
$0.5^{-2.5} = 2^{2.5} = 4 \times 1.41 \approx 5.66$,所以
$$t_{\text{MS}} \approx 5.7\times10^{10}\ \text{年}$$
約 570 億年。(真實的紅矮星因為內部完全對流、能把核心外的氫也送進去燒,實際壽命還更長,可達數兆年;簡單模型已經抓到「遠比太陽長壽」的本質。)
同一片星雲中誕生的兄弟姊妹,僅僅因為質量差了 30 倍,壽命就差了四個數量級。這正是「質量即命運」最赤裸的展現。
赫羅圖:恆星一生的「相空間地圖」
如果說結構方程式是恆星演化的「運動定律」,那麼赫羅圖就是它的「相空間」——一張把每顆恆星畫成一個點的地圖。橫軸是表面溫度(由左向右遞減,這是歷史慣例),縱軸是光度。一顆恆星在某個瞬間的狀態,就是圖上的一個點;它的一生,就是這個點走過的一條軌跡(evolutionary track)。
恆星的表面溫度、光度與半徑三者並不獨立,而被史蒂芬–波茲曼定律(Stefan–Boltzmann law) 綁在一起:
$$L = 4\pi R^2 \sigma T_{\text{eff}}^4$$
其中 $\sigma$ 是史蒂芬–波茲曼常數、$T_{\text{eff}}$ 是有效溫度。這條式子讓赫羅圖變得可被「解讀」:在圖上,沿著左下到右上的對角線是等半徑線。一顆星若在圖上往右上方移動(變冷卻又變亮),唯一的可能就是它膨脹了——這正是恆星離開主序、走向紅巨星(red giant)時發生的事。當核心的氫耗盡、核融合移到外圍的殼層,外層急速膨脹、表面冷卻變紅,但因半徑暴增,總光度反而飆升,於是恆星在赫羅圖上爬上「巨星支(giant branch)」。
理解這一點,赫羅圖就從一張「分類圖表」升級成一部「演化劇本」。主序帶不是一個分類抽屜,而是恆星花費約 90% 生命的「燃氫平衡態」聚集處;恆星演化的所有戲劇性轉折——膨脹、收縮、脈動、拋殼——都對應赫羅圖上一段可被星球結構模型計算出的軌跡。
看一個例子:太陽未來的軌跡
我們可以具體預測太陽在赫羅圖上的未來路線。此刻太陽是一顆穩定的主序星,位於圖中央偏下。約再過 50 億年,核心氫耗盡,核心收縮升溫、氫殼層燃燒啟動,外層膨脹——太陽會沿軌跡往右上方爬,膨脹成一顆紅巨星,半徑可能脹大到吞沒水星、金星,光度增為現在的上千倍,但表面溫度降到約 3000 K(偏紅)。
接著核心氦點燃(氦閃, helium flash),太陽短暫退回到較穩定的「水平分支(horizontal branch)」;氦也燒完後再度膨脹爬上漸近巨星支(asymptotic giant branch, AGB)。最後,太陽會把外層氣體溫和地拋散成行星狀星雲(planetary nebula),裸露的核心則收縮成一顆熾熱緻密的白矮星(white dwarf),在赫羅圖上落到左下角——又小、又熱、又暗——此後緩慢冷卻數百億年。整條軌跡,今天就能用電腦演化模型畫出來。
重點回顧
- 流體靜力平衡($\frac{dP}{dr} = -\frac{GM\rho}{r^2}$)是理解恆星的起點:重力與壓力梯度的平衡,決定了核心的壓力與溫度,並由此決定能否點燃核融合(恆星與棕矮星的分界約 $0.08\,M_\odot$)。
- 由量級估算可得核心溫度 $k_B T_c \sim \frac{GM\mu m_H}{R}$,本質是「重力壓縮做的功變成熱」,太陽核心因此達約 $1.5\times10^7\ \mathrm{K}$。
- 質光關係 $L\propto M^{3.5}$ 是結構方程式的自然結果;它直接導出主序壽命 $t_{\text{MS}}\propto M^{-2.5}$——質量越大壽命越短,差 30 倍質量可差四個數量級壽命。
- 赫羅圖是恆星一生的相空間地圖,配合史蒂芬–波茲曼定律 $L = 4\pi R^2\sigma T_{\text{eff}}^4$,圖上的位移可直接翻譯成恆星半徑與結構的變化(往右上 = 膨脹成巨星)。
- 「質量即命運」不是觀測巧合,而是少數幾條物理原理層層演繹出的必然——這正是恆星演化能被計算、被預測的根本原因。
深入探討(研究所視角)
把恆星物理嚴格化,需要聯立四條恆星結構方程式:流體靜力平衡、質量連續方程 $\frac{dM}{dr}=4\pi r^2\rho$、能量產生方程 $\frac{dL}{dr}=4\pi r^2\rho\,\varepsilon$($\varepsilon$ 為單位質量產能率),以及能量傳輸方程(輻射傳輸為 $\frac{dT}{dr}=-\frac{3\kappa\rho L}{16\pi a c r^2 T^3}$,$\kappa$ 為不透明度;當輻射梯度過陡則改由對流主導)。配合狀態方程 $P(\rho,T,X_i)$、不透明度 $\kappa(\rho,T,X_i)$ 與核反應率 $\varepsilon(\rho,T,X_i)$ 三組「微觀物理輸入」,再給定邊界條件,就能解出整顆恆星的結構。隨時間更新化學組成 $X_i$,便得到演化軌跡——這正是現代恆星演化碼(如 MESA)的骨架。
一個值得深思的理論結果是 Vogt–Russell 定理:在化學組成均勻的假設下,一顆恆星的結構幾乎完全由它的質量與組成決定。這是「質量即命運」的嚴格表述,也解釋了為何主序是一條窄帶而非散亂的雲。但這定理只是近似——真實恆星因核反應持續改變內部組成、因自轉、磁場、對流超射(convective overshooting)、質量流失與雙星交互作用而偏離單參數描述,這些「二階效應」正是當前恆星天文物理研究的前沿。
另一個關鍵概念是簡併壓力(degeneracy pressure) 接手的時刻。當核心密度高到電子遵守庖立不相容原理而產生與溫度幾乎無關的簡併壓力時,恆星的演化邏輯被改寫:氦閃之所以是「閃」,正因簡併核心的點火不會像理想氣體那樣靠膨脹自我調節,導致反應率失控暴衝直到簡併被解除。而錢德拉塞卡極限(Chandrasekhar limit, $\approx 1.4\,M_\odot$) 則標誌電子簡併壓力的上限——超過它,白矮星無法支撐,命運便交給中子簡併壓力(中子星)或徹底失敗(黑洞)。從理想氣體到簡併物質的這場「壓力來源接力」,才是恆星終局分歧的真正物理開關,也是把恆星天文與量子統計力學、廣義相對論縫合在一起的迷人交界。