銀河系進階：用恆星軌道倒帶一座星系的成長史

如果星空是一座檔案館，恆星的軌道就是館藏

入門篇告訴你銀河系是一座直徑十萬光年、繞著中心緩慢旋轉的恆星島嶼，太陽以每秒約 $220$ 公里在它的「郊區」公轉。這幅圖像很美，卻是一張靜態快照。現在請換一個問題：如果你能精準量出每一顆恆星此刻的位置與三維速度，你能不能反過來「倒帶」這座島嶼一百多億年的成長史？

這正是 2010 年代之後天文學的革命。歐洲太空總署的蓋亞衛星(Gaia) 已經測出超過十億顆恆星的位置、距離與運動，把銀河系從一張模糊的照片變成一部可以暫停、快轉、回放的影片。要讀懂這部影片，我們需要的不再只是「盤、棒、暈」的形狀描述，而是一整套銀河動力學(galactic dynamics) 的語言——本篇就帶你進入這套語言，看恆星的軌道如何洩漏星系的質量、歷史與未來。

微分轉動與奧爾特常數：恆星不像唱片那樣轉

入門篇說銀河系「旋轉」，但這個旋轉和黑膠唱片完全不同。唱片是剛體轉動(rigid-body rotation)，所有點的角速度 $\Omega$ 相同；銀河系卻是微分轉動(differential rotation)——不同半徑的恆星以不同的角速度繞行。在太陽附近，軌道速度 $v(R)$ 大致保持平坦（這正是暗物質的證據），於是角速度 $\Omega(R)=v/R$ 隨半徑遞減：內側恆星轉得比我們快，會「超車」；外側恆星轉得比我們慢，會「落後」。

要量化太陽附近這種剪切，荷蘭天文學家奧爾特(Jan Oort) 在 1927 年引入兩個局域常數。對鄰近恆星的視向速度與自行做傅立葉展開，會發現它們隨銀經 $\ell$ 呈正弦變化，係數正是奧爾特常數：

$$A = \frac{1}{2}\left(\frac{v}{R} - \frac{dv}{dR}\right), \qquad B = -\frac{1}{2}\left(\frac{v}{R} + \frac{dv}{dR}\right)$$

其中 $A$ 量度剪切(shear)，$B$ 量度渦度(vorticity)。它們有兩個漂亮的物理意義：

• $A - B = v/R = \Omega$，即太陽繞銀心的角速度；
• $A + B = -\,dv/dR$，即轉動曲線在太陽處的斜率。

蓋亞給出的現代值約為 $A \approx 15.3\ \mathrm{km\,s^{-1}\,kpc^{-1}}$、$B \approx -11.9\ \mathrm{km\,s^{-1}\,kpc^{-1}}$。由 $A+B \approx 3.4 > 0$（但很小）可知，轉動曲線在太陽附近幾乎平坦、略微下降——這是「我們只用太陽鄰域的恆星就能秤出整個銀河剪切」的經典範例。注意：奧爾特常數只描述局域線性近似，距離太陽幾百秒差距以外就需要完整的轉動曲線模型。

本輪近似：恆星其實在跳一支橢圓小舞

如果軌道都是完美圓形，那談「恆星檔案」就沒戲了，因為圓軌道記不住歷史。真實的恆星軌道是近圓但帶擾動的。處理這種軌道有一個極優雅的工具——本輪近似(epicyclic approximation)。

想像一顆恆星的軌道半徑在某個導引半徑 $R_0$ 附近小幅振盪。把運動分解成「導引中心以 $\Omega(R_0)$ 繞銀心轉」加上「恆星繞導引中心畫一個小橢圓」。這個小橢圓的徑向振盪頻率稱為本輪頻率(epicyclic frequency) $\kappa$：

$$\kappa^2 = R\,\frac{d\Omega^2}{dR} + 4\Omega^2 = 2\Omega\left(2\Omega + R\frac{d\Omega}{dR}\right)$$

用奧爾特常數可改寫得非常簡潔：

$$\kappa^2 = -4B(A - B) = -4B\,\Omega$$

對平坦轉動曲線（$v$ 為常數），可算出 $\kappa = \sqrt{2}\,\Omega$。這代表恆星每繞銀心一圈，會在徑向「進出」約 $\sqrt{2}$ 次——本輪不會閉合，於是軌道整體呈現玫瑰花形(rosette) 而非閉合橢圓。太陽自己就在做這支舞：它一邊繞銀心、一邊相對盤面上下振盪，垂直方向的振盪頻率 $\nu$ 通常比 $\kappa$ 還高，所以太陽每個銀河年會穿越銀盤中面好幾次。

這支「小舞」之所以重要，是因為它定義了共振(resonance)。當恆星的本輪頻率與旋臂或中央棒的圖樣速度 $\Omega_p$ 形成簡單整數比時，能量會持續累積，軌道被劇烈改變。最重要的兩個是：

• 共轉共振(corotation, CR)：$\Omega(R) = \Omega_p$，恆星與圖樣同步，這裡是旋臂的「不增不減」分界；
• 林布拉共振(Lindblad resonances)：$\Omega \pm \kappa/m = \Omega_p$（$m$ 為旋臂條數），內、外林布拉共振界定了密度波能穩定傳播的範圍。

入門篇講的「密度波理論」其實就建立在這套共振結構上：旋臂只能存在於內外林布拉共振之間，這也回答了「為什麼旋臂不會越捲越緊」的更深層原因。

動手算一下：太陽附近的本輪週期

我們來估算太陽完成一次徑向本輪振盪要多久。已知太陽的軌道角速度

$$\Omega = \frac{v}{R} = \frac{220\ \mathrm{km/s}}{8.2\ \mathrm{kpc}}.$$

先把 $8.2\ \mathrm{kpc}$ 換成公里（$1\ \mathrm{kpc} \approx 3.086\times10^{16}\ \mathrm{km}$）：

$$R \approx 8.2 \times 3.086\times10^{16} \approx 2.53\times10^{17}\ \mathrm{km}.$$

於是

$$\Omega \approx \frac{220}{2.53\times10^{17}} \approx 8.7\times10^{-16}\ \mathrm{s^{-1}}.$$

對平坦轉動曲線，$\kappa = \sqrt{2}\,\Omega \approx 1.23\times10^{-15}\ \mathrm{s^{-1}}$。本輪週期

$$T_\kappa = \frac{2\pi}{\kappa} \approx \frac{6.28}{1.23\times10^{-15}} \approx 5.1\times10^{15}\ \mathrm{s}.$$

換算成年（$1$ 年 $\approx 3.16\times10^7$ 秒）：

$$T_\kappa \approx \frac{5.1\times10^{15}}{3.16\times10^7} \approx 1.6\times10^{8}\ \mathrm{年}.$$

也就是約 $1.6$ 億年——比一個銀河年（約 $2.3$ 億年）短。兩個週期的比值 $T_{\text{銀河年}}/T_\kappa \approx \sqrt{2}$，正好對應 $\kappa = \sqrt{2}\,\Omega$。從一個速度、一個半徑，我們就推出了太陽在盤面裡「呼吸」的節奏，這就是動力學估算的力量。

蓋亞時代的銀河考古學：被消化的小星系

入門篇提到「銀河考古學」是透過恆星化學與運動重建歷史。蓋亞把這門學問推到了全新高度，因為它讓我們能在積分運動(integrals of motion) 的空間裡找出「同源」的恆星群。

關鍵概念是：恆星在盤面或暈中飄盪了上百億年，原本的空間位置早就被微分轉動「攪拌」均勻，再也看不出它們曾經屬於同一個系統。但有些物理量在平滑、近似不變的星系位能裡會被保留——例如能量 $E$ 與角動量 $L_z$。一群來自同一個被併吞的矮星系的恆星，即使如今散布全天，仍會在 $(E, L_z)$ 或所謂「作用量空間(action space)」裡聚成一團。這些殘骸就叫恆星流(stellar streams)。

2018 年蓋亞第二批資料公布後，研究者立刻發現銀暈裡有一大群恆星呈現高度徑向、幾乎不轉動的軌道，金屬量又自成一族。這被認證為一次重大併合事件的殘骸，命名為蓋亞–恩克拉多斯–香腸(Gaia-Enceladus / Gaia-Sausage)——名字裡的「香腸」正是因為這群恆星在速度空間 $(v_R, v_\phi)$ 裡拉成一條香腸形狀。它代表約 $80$ 至 $100$ 億年前，一個質量達數十億倍太陽的矮星系撞進年輕的銀河系，把當時的盤面加熱、打厚，並貢獻了今日銀暈相當大一部分的恆星。

這個發現徹底改寫了我們對銀河系的理解：它不是孤獨地、平滑地長大，而是靠不斷併吞小星系「滾雪球」累積質量——這正是「階層式結構形成(hierarchical structure formation)」在我們家門口的直接證據。今天還能看到的人馬座矮星系流(Sagittarius stream)、GD-1 流等，都是這場持續至今的「銀河系吃飯」現場。

看一個例子：厚盤是怎麼來的？

銀盤其實分兩層：薄盤(thin disk)厚約 $300$ 秒差距，由年輕、富金屬、軌道近圓的恆星組成；外面還包著一層厚盤(thick disk)，厚約 $1000$ 秒差距，恆星較老、較貧金屬、垂直速度彌散更大。

厚盤的成因長期是個謎。蓋亞–恩克拉多斯併合提供了一個有力答案：當那個矮星系以高速撞入時，它的引力擾動把原本薄薄的「古盤」加熱(dynamical heating)——恆星的垂直振盪振幅被泵大，盤就被「打厚」了。換句話說，厚盤可能是一場古老撞車事故的傷疤。

我們可以用化學標記驗證這個故事。厚盤恆星普遍有較高的 $[\alpha/\mathrm{Fe}]$ 比值（$\alpha$ 元素如氧、鎂來自快速的核心塌縮超新星，鐵則大量來自較慢的 Ia 型超新星）。高 $[\alpha/\mathrm{Fe}]$ 意味著這些恆星形成於早期、快速的恆星爆發階段，鐵還來不及被 Ia 型超新星大量釋放——時間戳記與「古盤被併合加熱」的劇本完全吻合。化學與動力學兩條線索互相印證，這就是現代銀河考古學的標準作法。

重返銀心：S 星、GRAVITY 與時空的彎曲

入門篇介紹了人馬座 A*(Sgr A*) 是約 $400$ 萬倍太陽質量的超大質量黑洞。進階地看，銀心是檢驗強重力下廣義相對論的天然實驗室，而關鍵主角是那群繞它狂奔的 S 星(S-stars)。

最著名的 S2 星軌道週期僅約 $16$ 年，近心點時距黑洞約 $120$ 個天文單位，速度飆到光速的百分之二以上。2018 年起，歐洲南方天文台用四台望遠鏡干涉合成的 GRAVITY 儀器，把 S2 通過近心點的過程量得無比精準，看到了兩個純粹來自廣義相對論、牛頓力學無法解釋的效應：

• 重力紅移(gravitational redshift)：S2 在近心點時，它發出的光要爬出黑洞的深重力井，波長被拉長。觀測到的紅移量與廣義相對論預言完全一致。
• 史瓦西進動(Schwarzschild precession)：S2 的橢圓軌道不會閉合，而是每繞一圈近心點方向就轉動一個小角度（正如水星近日點進動的星系尺度版本）。蓋亞與 GRAVITY 證實了這個進動，方向與大小都符合理論。

近年 GRAVITY 還在 Sgr A* 周圍偵測到閃焰(flares)——熱斑(hot spot)以接近光速繞行最內穩定圓軌道(ISCO)附近，在數十分鐘內畫出一個小圈。這幾乎是直接「看見」物質在黑洞事件視界外緣打轉。配合 2022 年事件視界望遠鏡(EHT)拍下的 Sgr A* 陰影影像，我們對自家黑洞的理解已經從「間接秤重」進化到「直接成像 + 時空檢驗」。

這裡值得再強調一個常被誇大的迷思：Sgr A* 雖然巨大，卻是個極度安靜的黑洞。它目前的吸積率非常低，亮度遠不及活躍星系核(AGN)。我們之所以能看到 S 星而不被輻射淹沒，正是因為它「沒在大吃」——銀心是一頭沉睡的巨獸，不是入門想像中那種貪婪的吸塵器。

重點回顧

• 銀河系是微分轉動而非剛體轉動：不同半徑角速度不同，太陽附近的剪切由奧爾特常數 $A$、$B$ 量化，且 $A-B=\Omega$、$A+B=-dv/dR$。
• 本輪近似把近圓軌道拆成「導引中心 + 小橢圓」，本輪頻率 $\kappa=\sqrt{-4B\Omega}$；平坦轉動曲線下 $\kappa=\sqrt{2}\,\Omega$，軌道呈不閉合的玫瑰花形。
• 旋臂與中央棒透過共轉共振與林布拉共振塑造盤面結構，也解釋密度波為何能穩定存在。
• 蓋亞在積分運動空間中辨識出恆星流，揭露銀河系靠併吞矮星系（如蓋亞–恩克拉多斯）成長，厚盤可能是古老併合加熱的傷疤。
• 銀心的 S 星 + GRAVITY 觀測直接檢驗了重力紅移與史瓦西進動；Sgr A* 是一頭安靜、低吸積率的超大質量黑洞，並非貪婪的吸塵器。

深入探討（研究所視角）

從碰撞較少系統到 Jeans 方程

要嚴格描述銀河系這樣的恆星系統，核心觀念是它是一個碰撞較少系統(collisionless system)：恆星之間兩兩近距離散射的弛豫時標遠超過宇宙年齡，所以個別恆星感受到的是平滑的整體位能 $\Phi(\mathbf{x},t)$，而非鄰星的逐一拉扯。系統的狀態由相空間的分布函數(distribution function) $f(\mathbf{x},\mathbf{v},t)$ 描述，其演化遵守無碰撞玻爾茲曼方程(collisionless Boltzmann equation)：

$$\frac{\partial f}{\partial t} + \mathbf{v}\cdot\nabla f - \nabla\Phi\cdot\frac{\partial f}{\partial \mathbf{v}} = 0,$$

並與位能透過帕松方程(Poisson equation) $\nabla^2\Phi = 4\pi G\rho$ 耦合。直接解 $f$ 通常不可行，於是我們對速度取矩(moments)，得到流體力學形式的 Jeans 方程。對穩態、軸對稱的盤，垂直方向的 Jeans 方程給出一個極實用的關係：

$$\frac{1}{\nu}\frac{\partial(\nu\,\sigma_z^2)}{\partial z} = -\frac{\partial \Phi}{\partial z} = -4\pi G\,\Sigma(z)\ \text{(近似)},$$

其中 $\nu$ 是示蹤恆星的數密度、$\sigma_z$ 是垂直速度彌散、$\Sigma$ 是盤的面密度。這就是經典的 Oort 局域質量(Oort limit) 問題：用太陽鄰域恆星的垂直運動「秤」出當地總質量密度，再扣掉可見恆星與氣體，就能檢驗是否需要局域暗物質。蓋亞讓這類測量精度空前提升，目前結果顯示太陽鄰域的盤暗物質貢獻很小——這與「暗物質呈球形暈分布、在盤面局域不顯著」的圖像一致。

作用量–角度變數與本輪的嚴格版本

本文用的本輪近似只是線性化結果。研究所層級的工具是作用量–角度變數(action–angle variables) $(\mathbf{J},\boldsymbol{\theta})$。對可積位能，三個作用量 $(J_R, L_z, J_z)$ 是嚴格的運動積分，分別對應徑向、方位與垂直方向的軌道「振幅」，而角度 $\boldsymbol{\theta}$ 隨時間線性增長。作用量的威力在於它們是絕熱不變量(adiabatic invariants)：當星系位能緩慢演化（例如盤逐漸增厚、或暗物質暈緩慢沉降）時，作用量近似守恆。這讓恆星流在作用量空間裡比在能量–角動量空間更乾淨地聚集，是現代恆星流搜尋（如用 AGAMA、galpy 等動力學程式庫）的數學基礎。本輪頻率 $\kappa$、垂直頻率 $\nu$、方位頻率 $\Omega$ 正是作用量空間裡的三個基頻 $\partial H/\partial J_i$，共振條件 $\ell\,\Omega + m\,\kappa + n\,\nu = 0$ 統一描述了棒、旋臂與衛星星系激發的所有共振現象。

動力摩擦與併合的命運

為什麼被併吞的矮星系最終會沉入銀心、被瓦解成恆星流？答案是動力摩擦(dynamical friction)。當一個質量 $M$ 的天體穿過背景恆星海時，會在身後激起一道密度尾跡(wake)，這道過密區反過來拉住它，造成減速。Chandrasekhar 給出的經典公式為：

$$\frac{d\mathbf{v}_M}{dt} = -\frac{4\pi G^2 M\,\rho\,\ln\Lambda}{v_M^3}\,\mathbf{v}_M,$$

其中 $\ln\Lambda$ 為庫倫對數，$\rho$ 為背景密度。減速正比於 $M$，所以越重的衛星沉得越快——這解釋了為何大質量併合（如蓋亞–恩克拉多斯）能快速把能量灌進古盤、製造厚盤，而小質量流（如 GD-1）能在暈裡飄盪數十億年仍保持細長。同一個動力摩擦機制也驅動著超大質量黑洞在星系併合後彼此沉降、最終合併，連接到重力波天文學的最前沿。把這些拼在一起，銀河系就不再是一張靜態的恆星地圖，而是一部仍在書寫、可被動力學方程逐頁解讀的演化史。