克卜勒定律的深層機制：從橢圓的物理來源到宇宙天平

如果克卜勒第三定律「失準」了，問題出在哪？

你大概已經很熟悉克卜勒三大定律：行星走橢圓、面積速度守恆、週期平方正比於半長軸立方。但這裡有個讓很多人卡住的問題：當天文學家把月球繞地球、木星繞太陽、雙星互繞的資料一起畫進 $T^2 \propto a^3$ 的對數圖時，會發現這些點並不全落在同一條線上——每個系統有自己的斜率截距。如果克卜勒定律是「定律」，為什麼它對不同系統會給出不同的常數？

答案藏在克卜勒自己看不到的地方。克卜勒是從第谷（Tycho Brahe）的觀測資料歸納出這些規律的，他並不知道背後的原因。真正讓這三條經驗法則升級成「可推導、可修正、可推廣」的物理定律的，是牛頓。這篇進階篇，我們不再重述三大定律的內容，而是要回答一個更深的問題：這些定律從哪裡來、在哪裡會破、又被推廣到多遠？

從平方反比力到橢圓：第一定律其實是個「巧合」

我們先釐清一件常被誤會的事：橢圓軌道不是萬有引力的必然結果，而是「束縛態」的特例。

牛頓萬有引力給出的是平方反比的中心力：

$$ \vec{F} = -\frac{G M m}{r^2}\,\hat{r} $$

在這個力場中解運動方程式，得到的軌道是圓錐曲線（conic sections）的通式：

$$ r(\theta) = \frac{p}{1 + e\cos\theta} $$

其中 $p$ 是半通徑（semi-latus rectum），$e$ 是離心率（eccentricity）。離心率決定了軌道是哪一種圓錐曲線：

• $e = 0$：圓
• $0 < e < 1$：橢圓（克卜勒研究的行星）
• $e = 1$：拋物線（恰好逃逸）
• $e > 1$：雙曲線（飛掠後一去不回）

換句話說，克卜勒第一定律之所以成立，是因為太陽系的行星恰好都處於束縛態（$e<1$）。2017 年穿越太陽系的星際天體 ʻOumuamua（1I/2017 U1）離心率約 $e \approx 1.2$，走的就是雙曲線——它從來不「繞」太陽，只是被太陽的重力場偏折一下就離開了。克卜勒第一定律對它完全不適用。

這個觀點很重要：第一定律不是普世的，它是「能量為負」這個條件的幾何表現。 系統的總力學能決定了命運：

$$ E = \frac{1}{2}mv^2 - \frac{GMm}{r} $$

$E<0$ 是橢圓，$E=0$ 是拋物線，$E>0$ 是雙曲線。離心率與能量、角動量的關係是：

$$ e = \sqrt{1 + \frac{2EL^2}{G^2 M^2 m^3}} $$

你可以直接從這條式子讀出來：當 $E<0$，根號內小於 1，$e<1$，軌道閉合。這就是橢圓的物理來源。

第二定律的真身：角動量守恆

克卜勒第二定律（等面積定律）常被當成行星的特殊性質，但它其實是任何中心力下都成立的——力是不是平方反比根本無所謂，只要力永遠指向同一個中心點即可。

理由是角動量守恆。中心力對中心點不產生力矩（$\vec{\tau} = \vec{r}\times\vec{F} = 0$，因為 $\vec{F}$ 與 $\vec{r}$ 平行），所以角動量 $\vec{L} = \vec{r}\times m\vec{v}$ 是常數。面積速度（areal velocity）正好是角動量的一半除以質量：

$$ \frac{dA}{dt} = \frac{1}{2}r^2\dot\theta = \frac{L}{2m} = \text{常數} $$

這意味著克卜勒第二定律比第一、第三定律「更基本」。即使有一天我們發現引力不是嚴格的平方反比（例如某些修正重力理論），只要力仍是中心力，等面積定律依然成立。這也是為什麼第二定律在三條裡最不容易被觀測推翻。

實務上，第二定律告訴我們行星在近日點（perihelion）跑得快、遠日點（aphelion）跑得慢。地球近日點在一月初、遠日點在七月初，這就是為什麼北半球的冬季（地球跑得快）比夏季短了約 4 天——北半球冬天，地球正好衝過近日點附近。

被牛頓改寫的第三定律：那個「常數」根本不是常數

這是本篇的核心。克卜勒原版的第三定律寫成：

$$ \frac{T^2}{a^3} = k \quad (\text{對繞同一顆太陽的所有行星，} k \text{相同}) $$

克卜勒以為 $k$ 是宇宙常數。牛頓用萬有引力一推導，發現真相是：

$$ T^2 = \frac{4\pi^2}{G(M+m)}\,a^3 $$

關鍵差別在分母的 $(M+m)$。原版第三定律假設行星質量 $m$ 可忽略（$m \ll M$），所以 $k \approx 4\pi^2/(GM)$ 只跟中心天體質量有關。一旦兩個天體質量相當（雙星、行星與大衛星），$(M+m)$ 就不能省略，「常數」就會隨系統改變——這正是開頭那張對數圖上不同系統落在不同線上的原因。

這個修正把克卜勒第三定律從「描述太陽系的經驗規則」升級成全宇宙最強大的稱重工具。只要量到繞行週期 $T$ 和軌道半長軸 $a$，就能反推系統總質量 $M+m$。天文學家測黑洞質量、測系外行星、測星系暗物質，骨子裡用的都是這條被牛頓改寫的第三定律。

動手算一下：用第三定律稱出太陽的質量

我們用地球的軌道資料來秤太陽。已知：

• 地球公轉週期 $T = 1\ \text{年} = 3.156\times10^7\ \text{s}$
• 軌道半長軸 $a = 1\ \text{AU} = 1.496\times10^{11}\ \text{m}$
• 重力常數 $G = 6.674\times10^{-11}\ \text{N·m}^2/\text{kg}^2$

因為地球質量遠小於太陽（$m_\oplus/M_\odot \approx 3\times10^{-6}$），可近似 $M+m \approx M_\odot$。把第三定律反解：

$$ M_\odot = \frac{4\pi^2 a^3}{G T^2} $$

代入：

$$ M_\odot = \frac{4\pi^2 (1.496\times10^{11})^3}{(6.674\times10^{-11})(3.156\times10^7)^2} $$

分子：$4\pi^2 \times 3.348\times10^{33} \approx 1.322\times10^{35}$

分母：$6.674\times10^{-11} \times 9.96\times10^{14} \approx 6.647\times10^{4}$

$$ M_\odot \approx \frac{1.322\times10^{35}}{6.647\times10^{4}} \approx 1.99\times10^{30}\ \text{kg} $$

這正是太陽的質量。注意我們完全沒有去太陽上量任何東西——只靠地球繞太陽的兩個幾何量，就秤出了一顆 $2\times10^{30}$ 公斤的恆星。這就是克卜勒第三定律在天文學裡無可取代的威力：運動就是天平。

如果你想感受尺度：$1.99\times10^{30}$ 公斤大約是地球質量的 33 萬倍。而這整顆太陽，在第三定律眼中只是公式分母裡的一個數字。

第三定律的推廣：質心、約化質量與系外行星

當兩個天體質量相當時，它們並不是一個繞另一個轉，而是兩者一起繞共同質心（barycenter）轉。質心位置由槓桿原理決定：

$$ M\,r_1 = m\,r_2, \qquad r_1 + r_2 = a $$

質量大的那顆離質心近、畫小圈；質量小的離質心遠、畫大圈。地球與月球的質心其實在地球內部（距地心約 4670 公里，仍在地表以下），所以看起來像月球繞地球；但太陽與木星的質心落在太陽表面之外，太陽其實在繞著一個太陽外的點微微擺動。

這個「恆星擺動」正是早期偵測系外行星的核心原理——徑向速度法（radial velocity / Doppler method）。行星雖然看不見，但它拉著母恆星繞質心擺動，使恆星光譜週期性紅移／藍移。從擺動的週期和振幅，配合改寫後的第三定律，就能反推行星的最小質量。

完整的兩體問題透過約化質量（reduced mass） $\mu = \frac{Mm}{M+m}$，可以化簡成「一個質量為 $\mu$ 的虛擬粒子繞固定中心運動」的單體問題，而那個虛擬粒子的軌道週期，正好滿足牛頓版第三定律。這是把克卜勒從太陽系推廣到任意兩體系統的數學鑰匙。

重點回顧

• 第一定律是束縛態的特例：橢圓只在總能量 $E<0$ 時出現；$E\ge0$ 給出拋物線或雙曲線，像 ʻOumuamua 那樣有去無回。離心率由能量與角動量共同決定：$e=\sqrt{1+2EL^2/(G^2M^2m^3)}$。
• 第二定律的本質是角動量守恆：它對任何中心力都成立，不限於平方反比，因此是三條裡最「穩固」的一條。
• 牛頓把第三定律的常數改成 $4\pi^2/[G(M+m)]$：原版漏掉了 $(M+m)$，只在 $m\ll M$ 時才近似成立；質量相當時「常數」就因系統而異。
• 第三定律是宇宙天平：只靠週期與半長軸兩個幾何量，就能秤出太陽、黑洞、系外行星系統的總質量，無需接觸天體本身。
• 兩體要繞共同質心：質量相當時不是一繞一，而是同繞質心；恆星的擺動正是徑向速度法偵測系外行星的依據。

深入探討（研究所視角）

水星近日點進動與廣義相對論。 克卜勒—牛頓框架預測束縛軌道是封閉的橢圓——行星每圈都回到同一個近日點。但水星的近日點實際上每世紀額外多進動約 $43''$（角秒），這是牛頓力學連把所有已知行星攝動都算進去後仍無法解釋的殘差。這個微小偏差最終由愛因斯坦的廣義相對論（General Relativity）解決：在彎曲時空中，引力不再是嚴格的 $1/r^2$，有效位能多了一項 $\propto 1/r^3$ 的修正，使軌道變成緩慢旋轉的玫瑰花瓣（rosette）而非封閉橢圓。進動率為：

$$ \Delta\varphi = \frac{6\pi G M}{c^2 a (1-e^2)}\ \text{每軌道} $$

這正是克卜勒第一定律「封閉橢圓」失效之處——它是牛頓引力的精確結果，而牛頓引力本身只是相對論在弱場低速下的近似。

Bertrand 定理：為何偏偏是平方反比？ 1873 年 Bertrand 證明，在所有中心力中，只有兩種力場能讓所有束縛軌道都封閉：平方反比力（$\propto 1/r^2$，引力）與線性回復力（$\propto r$，三維諧振子）。任何其他冪次的中心力，軌道一般都會進動、不封閉。這解釋了為什麼克卜勒能觀測到「乾淨的封閉橢圓」——宇宙的引力恰好落在這兩種特殊力場之一上。這也說明水星進動所揭示的，是廣義相對論對純平方反比的微小偏離。

Laplace–Runge–Lenz 向量與隱藏對稱性。 平方反比力場有一個額外的守恆量——Laplace–Runge–Lenz 向量 $\vec{A} = \vec{p}\times\vec{L} - GMm^2\hat{r}$，它永遠指向近日點方向且大小不變，這正是軌道封閉、近日點不動的數學保證。從群論看，這對應克卜勒問題隱藏的 $SO(4)$ 對稱性（高於一般中心力只有的 $SO(3)$ 旋轉對稱）。當引力偏離純平方反比（如相對論修正），這個向量不再守恆，近日點便開始進動。這套對稱性分析後來在量子力學中對氫原子能階的「意外簡併」給出了同構的解釋——克卜勒問題與氫原子，數學上是同一個故事。

多體系統與混沌。 牛頓兩體問題有解析解，但三體以上一般沒有封閉形式解（Poincaré）。太陽系實際上是弱混沌系統：克卜勒橢圓只是零階近似，行星間的長期攝動（secular perturbation）使軌道要素在數百萬年尺度上緩慢演化。現代研究用 Laplace–Lagrange 攝動理論與數值積分追蹤這些演化，發現太陽系內行星軌道的 Lyapunov 時間約僅 500 萬年——意味著克卜勒的「精準鐘錶宇宙」在足夠長的時間尺度上，其實藏著不可預測性。