雙星與變星（進階）：質量函數、振動模態與距離階梯的當代危機

一條速度曲線藏著一個你永遠看不到的伴星

入門篇裡，我們用牛頓版克卜勒第三定律「秤」出一對目視聯星(visual binary)的質量。但那有一個奢侈的前提：你得能看見兩顆星，量出它們各自繞質心的軌道。問題是，宇宙中絕大多數聯星都太遠、靠得太近，望遠鏡裡只是一個無法分辨的光點。你連第二顆星都看不見，要怎麼稱它的體重？

更尖銳的問題是：天文學家後來在這類「看不見的伴星」裡，發現了第一批恆星級黑洞。一個你永遠看不到的天體，怎麼證明它的質量大到只能是黑洞？答案不在影像裡，而藏在一條徑向速度曲線(radial velocity curve)的形狀與振幅之中。這篇進階篇，就從入門篇略過的「光譜聯星到底怎麼解」出發，一路深入到脈動變星的振動模態、聯星演化的災難性結局，以及現代巡天如何把天空變成一座永不關閉的變星實驗室。

質量函數：從都卜勒擺動逆推看不見的伴星

入門篇提到光譜聯星(spectroscopic binary)時，只說「譜線會週期性偏移」。現在我們把這句話變成可以計算的物理。

假設一顆我們看得見的恆星（質量 $M_1$）繞著系統質心運動，它的譜線因都卜勒效應前後偏移，讓我們測出它沿視線方向的速度隨時間變化。若軌道接近圓形，這條速度曲線是一條正弦波，其振幅 $K_1$ 就是這顆星徑向速度的最大值。

關鍵在於：徑向速度只反映速度沿視線的分量。如果軌道平面相對視線傾斜了一個角度 $i$（傾角，inclination，$i=90^\circ$ 表示軌道側對我們、看得最清楚），那麼我們量到的振幅其實是真實軌道速度乘上 $\sin i$。把克卜勒運動與質心關係合併，可以推導出聯星研究的核心方程——質量函數(mass function)：

$$ f(M) = \frac{(M_2 \sin i)^3}{(M_1 + M_2)^2} = \frac{P\,K_1^3}{2\pi G}\,(1-e^2)^{3/2} $$

右邊全都是可觀測量：軌道週期 $P$、速度振幅 $K_1$、軌道離心率 $e$。換句話說，只要量出一顆看得見的星的速度曲線，不需要看見伴星，就能算出左邊這個只含 $M_1$、$M_2$、$i$ 的組合。

這個式子既強大又惱人。強大，是因為它把「看不見的伴星」拉進了可計算的範圍；惱人，是因為單靠一條速度曲線，$M_2$、$i$ 兩個未知數綁在一起，無法分離——這就是聯星天文學最根本的傾角簡併(inclination degeneracy)。

看一個例子：質量函數設下的質量下限

質量函數最漂亮的用途，是它能設下一個伴星質量的絕對下限。注意到 $\sin i \le 1$，而且 $(M_1+M_2)^2 \ge M_2^2$，於是：

$$ f(M) = \frac{(M_2 \sin i)^3}{(M_1+M_2)^2} \le \frac{M_2^3}{M_2^2} = M_2 $$

也就是說，質量函數的數值本身，就是伴星質量的最小可能值。

舉個具體情境：我們觀測到一個系統，看得見的恆星速度曲線給出 $P = 5.6$ 天、$K_1 = 470\ \mathrm{km/s}$、軌道近乎圓形（$e\approx 0$）。把這些代進去（換算單位後），算出 $f(M) \approx 6\ M_\odot$。

這代表什麼？不論傾角是多少、不論看得見那顆星多重，伴星至少有 $6\ M_\odot$。如果光譜還告訴我們看得見的那顆只是普通的 B 型恆星、而我們在那個位置看不到任何發光天體——一個質量超過 $6\ M_\odot$、卻完全不發光的緻密天體，候選名單上幾乎只剩黑洞。歷史上著名的天鵝座 X-1(Cygnus X-1)正是循這條邏輯被確認為黑洞候選者。質量函數，就是天文學家「稱黑洞」的第一把尺。

要徹底解出 $M_2$ 與 $i$，得補上額外資訊：若這同時是食聯星，食的存在本身就要求 $i$ 接近 $90^\circ$；若能看到雙線光譜（兩顆星的譜線都看得見，稱雙譜光譜聯星 double-lined SB2），則 $K_1/K_2 = M_2/M_1$ 可直接給出質量比。資訊一層層補上，簡併才一層層解開。

變星動物園：不只造父變星

入門篇把造父變星(Cepheid)當主角，因為它撐起了宇宙距離階梯。但脈動變星(pulsating variable)是一個龐大的家族，住在赫羅圖(Hertzsprung–Russell diagram)上一條斜貫的不穩定帶(instability strip)裡。不同質量、不同演化階段的恆星穿過這條帶時，會以不同方式脈動：

• 造父變星：高光度黃色超巨星，週期 1–100 天，週光關係的主角，用於測量遙遠星系。
• 天琴座 RR 型變星(RR Lyrae)：低質量、老年的水平分支星，週期短（約 0.2–1 天）、光度幾乎固定在 $M_V \approx +0.6$ 附近。因為又老又遍布球狀星團與銀暈，是測量銀河系結構與近鄰距離的利器。
• 盾牌座 δ 型變星(δ Scuti)：不穩定帶下端的中等質量星，週期僅數小時，常同時激發多個振動模態。
• 芻藁型長週期變星(Mira variables)：演化到漸近巨星分支(AGB)的紅巨星，週期長達數百天、振幅可達好幾個星等，外層鬆散、強烈脈動並拋出物質。

把這些放在一起看，會浮現一個統一的圖像：脈動是恆星結構失穩的自然反應，而它停在赫羅圖的哪個位置、由什麼物理驅動，決定了它脈動的週期與振幅。

脈動的物理：駐波、基頻與泛音

入門篇用 $\kappa$ 機制解釋了「為什麼脈動會持續」（不透明度像閥門，把能量一存一放）。進階篇要問另一個問題：脈動的週期是由什麼決定的？

把恆星想成一個會共鳴的球體。脈動本質上是一種徑向聲波駐波(radial standing acoustic wave)——壓力擾動以聲速在恆星內部來回傳播，形成穩定的振動模態，就像管風琴的管子有基音與泛音一樣。最簡單的基頻模態(fundamental mode)，整顆星同步膨脹收縮，中心是唯一的節點；更高階的泛音(overtone)則在恆星內部出現額外的節點面，週期較短。

脈動週期與恆星平均密度的關係（入門篇曾給出 $P\sqrt{\bar\rho}\approx Q$）可以用一個更直觀的方式理解：聲波橫越恆星一趟的時間，大致就是脈動週期的尺度。聲速正比於 $\sqrt{T}$、而恆星越蓬鬆（密度越低）路徑越長，於是低密度的巨星週期長、緻密的星週期短。這就是為什麼週期可以反推恆星的物理狀態。

動手算一下：用脈動「秤」出恆星密度

脈動關係 $P\sqrt{\bar\rho} = Q$ 提供一個迷人的副產品：只要測到週期，就能估出恆星的平均密度。取造父變星典型的脈動常數 $Q \approx 0.04$ 天（以恆星平均密度相對太陽的比值 $\bar\rho/\bar\rho_\odot$ 為單位來表示）：

$$ P\,\sqrt{\bar\rho/\bar\rho_\odot} = Q \quad\Rightarrow\quad \bar\rho/\bar\rho_\odot = \left(\frac{Q}{P}\right)^2 $$

對一顆週期 $P = 10$ 天的造父變星：

$$ \bar\rho/\bar\rho_\odot = \left(\frac{0.04}{10}\right)^2 = (4\times10^{-3})^2 = 1.6\times10^{-5} $$

太陽平均密度約 $1.41\ \mathrm{g/cm^3}$，所以這顆造父變星的平均密度只有約 $2\times10^{-5}\ \mathrm{g/cm^3}$——比地球海平面空氣（約 $1.2\times10^{-3}\ \mathrm{g/cm^3}$）還稀薄百倍。這正具體說明了造父變星是何等蓬鬆膨脹的超巨星：質量是太陽的好幾倍，體積卻撐到密度比空氣還低。脈動，等於把恆星內部的密度資訊「廣播」了出來。

當一顆星同時激發多個模態（例如 δ Scuti 或脈動白矮星），不同模態深入恆星的不同深度，攜帶不同層次的內部資訊。解讀這些頻率、反推恆星內部結構的學問，叫做星震學(asteroseismology)——這是現代恆星物理最熱門的前沿之一，我們留到最後再談。

聯星演化的災難劇場

入門篇的研究所視角提到了洛希瓣(Roche lobe)溢流與質量轉移。這裡我們把鏡頭拉遠，看質量轉移如何把一對平凡的恆星，推向宇宙中最劇烈的事件。

當較重的星先演化膨脹、填滿洛希瓣，物質流向伴星。若轉移溫和而穩定，系統可能變成一個 X 射線聯星——物質落入緻密伴星的吸積盤、釋放重力位能成為 X 射線。但若轉移失控、或兩顆星都膨脹到彼此的物質混在一起，系統會進入共同包層(common envelope)階段：兩顆恆星的核心一起埋在同一團氣體裡，像兩顆攪在一鍋濃湯裡的櫻桃。摩擦使兩核心快速向內螺旋、軌道急遽收縮，最後要嘛甩掉整個包層、留下一對靠得極近的緻密核心，要嘛直接併合。

這個「軌道收縮」的後果極為深遠。兩顆緻密殘骸（白矮星、中子星或黑洞）若被推得夠近，重力波輻射會持續帶走軌道能量，使它們在數百萬到數十億年內緩慢盤旋接近，最終以接近光速的速度相撞併合。2017 年的 GW170817——兩顆中子星併合——正是這樣一個系統的終曲，它同時被重力波偵測器與電磁望遠鏡捕捉到，開啟了多信使天文學(multi-messenger astronomy)的時代，也直接證實了宇宙中約一半的金、鉑等重元素，是在這類中子星併合中鍛造出來的。

換句話說，聯星演化不是恆星物理的一個邊角，而是一條把「兩顆普通恆星」與「重力波、伽瑪射線暴、宇宙重元素來源」串起來的主線。你手上那枚金戒指，很可能來自某對遠古聯星的死亡之舞。

巡天時代：把整片天空變成變星實驗室

入門篇的測距邏輯，假設你已經「找到」一顆造父變星並長期盯著它。但要怎麼從幾億顆星裡，自動篩出會變亮變暗的那一小撮？

這正是現代時域天文學(time-domain astronomy)的核心任務。像 OGLE、Gaia、ZTF(Zwicky Transient Facility) 這類巡天，反覆掃描同一片天空，為每顆星建立一條長期的亮度曲線(light curve)，再用演算法自動偵測週期性與爆發。一個關鍵工具是Lomb–Scargle 週期圖——它能在觀測時間點不規則（天氣、白天、季節都會打斷觀測）的資料裡，找出隱藏的週期訊號，這比均勻取樣的傅立葉分析困難得多，卻是真實天文資料的常態。

這套方法把變星研究從「一次盯一顆」升級成「一次處理數十億顆」。它的副產品改寫了多個領域：系外行星的凌日法(transit method)本質上就是在找一種特殊的「食」；而把數萬顆 RR Lyrae 的距離與運動拼起來，天文學家得以重建銀河系的三維結構與並吞矮星系的化石痕跡。變星，從一顆顆孤立的奇景，變成了測繪整個星系的統計工具。

重點回顧

• 質量函數 $f(M) = P K_1^3 (1-e^2)^{3/2}/(2\pi G)$ 只用可觀測的週期、速度振幅與離心率，就能逼出看不見伴星的質量；其數值本身即為伴星質量的下限，是辨認黑洞候選者的關鍵工具。
• 單條徑向速度曲線存在傾角簡併（$M_2$ 與 $\sin i$ 綁在一起），需靠食聯星（$i\approx90^\circ$）或雙譜光譜（質量比）等額外資訊解開。
• 脈動變星是一個大家族——造父變星、RR Lyrae、δ Scuti、Mira——都住在赫羅圖的不穩定帶上，由演化階段與質量決定脈動性質。
• 脈動是恆星內部的聲波駐波，有基頻與泛音；$P\sqrt{\bar\rho}\approx Q$ 讓我們能由週期反推恆星平均密度，也是星震學的基礎。
• 緊密聯星經由質量轉移、共同包層、重力波盤旋，最終可導向 X 射線聯星、Ia 型超新星與中子星併合，是重力波與宇宙重元素的源頭。

深入探討（研究所視角）

星震學：聽見恆星的內部

恆星脈動遠不只「整顆星一起膨脹收縮」這麼簡單。一般而言，恆星的振動可分解為一組三維的簡正模態(normal modes)，每個模態以三個量子數標記：徑向階數 $n$、球諧函數的角度數 $\ell$ 與 $m$。其中以壓力為回復力的稱 p 模態(pressure modes)，主要探測恆星外層；以浮力為回復力的稱 g 模態(gravity modes)，深入恆星核心。

對類太陽恆星，外層對流會隨機激發數千個低振幅的 p 模態，形成一系列近乎等間距的頻率峰，這個間距（大頻率間隔 $\Delta\nu$）與恆星平均密度直接相關：

$$ \Delta\nu \propto \sqrt{\bar\rho} \propto \sqrt{M/R^3} $$

而頻率峰的包絡最大值位置 $\nu_{\max}$ 則與表面重力和溫度有關。把這兩個全域震波參數搭配，便能在不依賴聯星、不依賴模型的情況下，獨立給出恆星的質量與半徑——這就是著名的星震標度關係(asteroseismic scaling relations)。克卜勒(Kepler)與 TESS 太空望遠鏡的高精度測光，已用這套方法測定了數萬顆紅巨星的質量與年齡，為「銀河系考古學(galactic archaeology)」提供了關鍵的年齡標尺。一顆恆星的振動頻譜，等於它的內部結構指紋。

非線性脈動與週期倍化

入門篇的 $\kappa$ 機制是線性圖像：擾動小、振盪規則。但真實的大振幅脈動星（如 RV Tauri 型、某些 BL Herculis 型造父變星）會展現非線性行為——脈動振幅大到模態之間發生共振耦合，出現週期倍化(period doubling)：相鄰的兩次脈動高低交替，週期實質上翻倍。連續的倍化序列正是混沌動力學裡通往確定性混沌(deterministic chaos)的經典路徑，少數變星的亮度曲線已被證實具有低維混沌的特徵。脈動變星因此意外成為非線性動力學的天然實驗室：恆星不只是測距工具，它的振盪本身就是一個值得研究的非線性系統。

哈伯張力與距離階梯的當代危機

入門篇結尾提到造父變星與 Ia 型超新星構成宇宙距離階梯。當代最受矚目的問題，是這條階梯量出的哈伯常數 $H_0$ 與宇宙微波背景(CMB)推出的數值之間的哈伯張力(Hubble tension)：以 SH0ES 團隊為代表、用造父變星校準的局域測量給出約 $73\ \mathrm{km\,s^{-1}\,Mpc^{-1}}$，而普朗克(Planck)衛星由早期宇宙模型外推得約 $67$，兩者差距已達 $5\sigma$ 以上，遠超過彼此的誤差棒。

這個張力之所以難以化解，正因為兩端的測量都被反覆檢驗：造父變星端必須處理金屬量對週光關係的影響、星際塵埃紅化、以及擁擠星場中的光度混染。為了交叉驗證造父變星，天文學家近年大力發展獨立的距離指標，例如尖端紅巨星法(TRGB，利用紅巨星分支頂端的固定光度）與富碳米拉變星。若多種彼此獨立的局域方法持續指向偏高的 $H_0$，那麼問題或許不在測量，而在我們的宇宙學標準模型 $\Lambda$CDM 本身——是否需要新的物理（如早期暗能量）。一條源自「恆星眨眼」的測距傳統，如今正站在挑戰整個宇宙學框架的最前線。