彗星、小行星與矮行星：非引力力、共振與撞擊物理

一顆彗星為什麼會「無緣無故」改變軌道？

你在入門篇已經知道彗星來自遙遠的歐特雲（Oort Cloud），會在接近太陽時長出壯麗的彗尾。但這裡有個令天文學家頭痛了兩個世紀的謎題：當我們用萬有引力把太陽、各大行星全部算進去，精準預測一顆週期彗星下次回歸的時刻，它卻常常早到或遲到好幾天。1986 年哈雷彗星（1P/Halley）回歸時，純引力模型的誤差就達到約四天。

四天聽起來不多，但在天體力學裡，這是一個巨大的、系統性的偏差——它不是隨機誤差，而是每一圈都朝同一個方向累積。是什麼力量，能在真空中持續推著一顆幾公里大的冰球？

這篇進階文章不再重述「彗星是髒雪球」這類基礎，而是要帶你進入小天體動力學的核心：非引力力（non-gravitational forces）、共振與混沌、撞擊物理，以及小行星如何被陽光「悄悄搬運」。這些機制不只解釋了上面的謎題，更是當代行星防禦（planetary defense）與太空採礦的理論基石。

非引力力：彗星自帶的火箭引擎

謎底其實藏在彗星的本質裡。當彗核（nucleus）接近太陽，表面的冰昇華（sublimation，固態直接變氣態）噴出氣體與塵埃。根據牛頓第三定律，噴氣會給彗核一個反作用推力——彗星等於自帶了一具不對稱、會隨日距變化的天然火箭引擎。

天文學家用一組經驗公式描述這個加速度。Marsden、Sekanina 與 Yeomans 在 1973 年提出的標準模型，把非引力加速度分解到三個方向：徑向（radial，沿日彗連線）、橫向（transverse，沿運動方向）、法向（normal，垂直軌道面）：

$$ \vec{a}_{\rm ng} = \left(A_1 \hat{r} + A_2 \hat{t} + A_3 \hat{n}\right) g(r) $$

其中 $A_1, A_2, A_3$ 是各方向的強度係數，而關鍵是那個隨日距 $r$ 變化的標度函數 $g(r)$：

$$ g(r) = \alpha \left(\frac{r}{r_0}\right)^{-m} \left[1 + \left(\frac{r}{r_0}\right)^{n}\right]^{-k} $$

這個函數是針對水冰昇華率量身打造的：在 $r_0 \approx 2.808\,\mathrm{AU}$（水冰昇華開始顯著的日距）內側，昇華旺盛、推力強；到了更遠處則急遽衰減。其中 $m=2.15,\ n=5.093,\ k=4.6142$ 是擬合出來的常數。

真正讓彗星「遲到」的，是那個橫向分量 $A_2$。為什麼？因為彗核會自轉，而表面被加熱後昇華有熱滯後（thermal lag）——就像午後兩三點才是一天最熱的時候，彗核「正午」過後才噴得最猛。這使噴氣方向偏離正對太陽的方向，產生一個沿軌道切線的淨推力。如果這個推力與運動方向相反，彗星就被持續「煞車」，軌道半長軸略增、週期變長——於是一圈又一圈地遲到。哈雷彗星正是如此。

這不是什麼神祕力量，而是自轉 × 熱慣性 × 昇華三者耦合的必然結果。

亞爾科夫斯基效應：陽光也能搬動一座山

如果說彗星的非引力力來自昇華噴氣，那對於不太揮發的小行星，難道就純粹乖乖聽引力的話嗎？並非如此。有一個更精緻、更安靜的力量在長期作用——亞爾科夫斯基效應（Yarkovsky effect）。

想像一顆自轉的小行星。它向陽面被曬熱，但因為熱慣性，這份熱量不會立刻輻射出去，而是被自轉「帶」到了下午面甚至昏側才以紅外線輻射出去。光子帶有動量，這份延遲的熱輻射就在偏離正午的方向上給了小行星一個微小但持續的反衝力。

這個加速度極其微弱，量級大約是：

$$ a_{\rm Yark} \sim \frac{F_\odot}{c}\,\frac{\pi R^2}{m} \sim \frac{L_\odot}{4\pi r^2 c}\cdot \frac{3}{4\rho R} $$

對一顆公里級小行星，$a_{\rm Yark}$ 約為 $10^{-13}\,\mathrm{m/s^2}$ 量級——比地表重力小了十幾個數量級。但天體力學的精髓在於時間：在數百萬年的尺度上，這點微力足以讓小行星的軌道半長軸漂移數十萬公里。

亞爾科夫斯基效應有兩種模式： - 日週亞爾科夫斯基（diurnal）：由自轉造成，效果方向取決於自轉是順行還是逆行。順行自轉的小行星軌道向外漂移，逆行則向內。 - 季節亞爾科夫斯基（seasonal）：由公轉造成的季節性溫差驅動，恆使軌道向內收縮。

這個效應正是為什麼小行星帶能持續「補貨」近地小行星（NEA）的關鍵——主帶小行星被亞爾科夫斯基緩慢推向附近的共振區，再被共振踢出來。而它的「旋轉版」表親 YORP 效應（Yarkovsky–O'Keefe–Radzievskii–Paddack），則因為小行星形狀不規則，使熱輻射產生淨力矩，能改變自轉速率與自轉軸方向，甚至把小行星轉到解體、形成雙小行星系統。

OSIRIS-REx 探測器在貝努（Bennu）上實測到的亞爾科夫斯基漂移約為每年 $284\,\mathrm{m}$——這正是把貝努列為潛在威脅天體、並需要精算其 2182 年撞擊機率時，不能忽略的因素。

共振與混沌：小行星帶的「禁區」是怎麼來的

把主帶小行星的軌道半長軸畫成直方圖，你會看到幾道明顯的空隙（Kirkwood gaps）。這些空隙的位置不是隨機的，而是精準對應到與木星的平均運動共振（mean-motion resonance, MMR）。

所謂 $p:q$ 共振，是指小行星公轉 $p$ 圈時木星剛好公轉 $q$ 圈。以 3:1 共振為例，由克卜勒第三定律 $a \propto T^{2/3}$：

$$ a_{3:1} = a_{\rm Jup}\left(\frac{1}{3}\right)^{2/3} \approx 5.2 \times 0.481 \approx 2.50\,\mathrm{AU} $$

這恰好就是其中一道最顯著的 Kirkwood gap。為什麼共振會「清空」軌道？因為在共振處，木星的引力擾動週期性地在同一相位累積，而非互相抵消。小行星的偏心率（eccentricity）被一點一點推高，最終高到讓它的軌道穿越火星甚至地球，於是被踢出共振區或墜入內太陽系。

看一個例子：3:1 共振如何製造近地小行星

數值模擬顯示，一顆掉進 3:1 共振的小行星，其偏心率可以在短短數十萬到一百萬年內從接近圓形（$e \approx 0.1$）被激發到 $e > 0.6$。當 $e$ 夠大，近日點距離 $q = a(1-e)$ 就會降到 1 AU 以內：

$$ q = 2.50 \times (1 - 0.6) = 1.0\,\mathrm{AU} $$

此時這顆小行星就成了會穿越地球軌道的近地天體。這正是內太陽系小行星與隕石的主要「產線」之一。Jack Wisdom 在 1980 年代用混沌動力學首次定量解釋了 3:1 共振的這種行為，證明小行星帶不是一池靜水，而是一個由共振主導、充滿混沌通道的動態系統。

混沌（chaos）在這裡有精確意義：相鄰軌道會指數發散，發散速率由李雅普諾夫指數（Lyapunov exponent）$\lambda$ 刻畫，特徵時間 $t_{\rm Lyap} = 1/\lambda$。對共振區小行星，$t_{\rm Lyap}$ 可短至數千年——意味著超過這個時間尺度，個別軌道的長期預測在原則上就不可能，只能談統計分布。

矮行星與流體靜力平衡：那道「夠圓」的門檻

入門篇告訴你矮行星（dwarf planet）的定義包含「達到流體靜力平衡（hydrostatic equilibrium）而呈近球形」。但這道門檻的物理究竟是什麼？

關鍵在於天體自身的重力是否能壓垮其物質強度。小天體靠固體強度（material strength）維持不規則形狀；一旦質量夠大，中心壓力超過岩石或冰的屈服強度，物質就會像流體一樣流動、填平凸起，趨於球形。

我們可以做個量級估計。一個半徑 $R$、密度 $\rho$ 的天體，其中心壓力約為：

$$ P_c \sim \frac{2}{3}\pi G \rho^2 R^2 $$

當 $P_c$ 超過材料屈服強度 $Y$（岩石約 $10^8\,\mathrm{Pa}$，冰約 $10^6$–$10^7\,\mathrm{Pa}$），天體就會塌縮成球。

動手算一下：冰質天體的「變圓臨界半徑」

取冰的密度 $\rho \approx 1000\,\mathrm{kg/m^3}$、屈服強度 $Y \approx 1\times10^7\,\mathrm{Pa}$，求中心壓力剛好等於 $Y$ 的臨界半徑 $R$：

$$ \frac{2}{3}\pi G \rho^2 R^2 = Y \;\Rightarrow\; R = \sqrt{\frac{3Y}{2\pi G \rho^2}} $$

代入數值（$G = 6.67\times10^{-11}$）：

$$ R = \sqrt{\frac{3 \times 10^7}{2\pi \times 6.67\times10^{-11} \times (10^3)^2}} = \sqrt{\frac{3\times10^7}{4.19\times10^{-4}}} \approx \sqrt{7.16\times10^{10}} \approx 2.7\times10^{5}\,\mathrm{m} $$

也就是約 270 公里。這個量級和觀測驚人地吻合：冰質天體大約在直徑 400 公里以上開始普遍呈球形（如土衛一 Mimas 直徑約 396 公里就已接近球形），而岩質天體因強度更高，門檻則更大（小行星灶神星 Vesta 直徑約 525 公里仍未完全達到平衡）。這解釋了為什麼冰質的矮行星普遍比岩質天體更早變圓——冰「軟」，更容易被自身重力捏圓。

當然，現實更複雜：自轉會讓天體呈扁球（鬩神星 Eris、妊神星 Haumea 因快速自轉而明顯拉長），「流體靜力平衡」其實是引力與離心力共同決定的平衡形狀。

撞擊物理：當動能取代化學

最後，談一個小天體研究中最暴力也最實用的面向——高速撞擊（hypervelocity impact）。這是隕石坑、行星防禦與小行星演化的共同核心。

小天體撞擊的相對速度動輒每秒數公里到數十公里，遠超岩石中的聲速（約 $5\,\mathrm{km/s}$）。在這種速度下，撞擊不再是「撞凹一個坑」這麼簡單，而是激震波（shock wave）把動能瞬間轉成熱與壓力，岩石被汽化、熔融、噴出。一次撞擊釋放的動能為：

$$ E_k = \frac{1}{2}m v^2 = \frac{1}{2}\cdot\frac{4}{3}\pi\rho R^3 \cdot v^2 $$

看一個例子：一顆 50 公尺小行星的能量

取 $R = 25\,\mathrm{m}$、$\rho = 3000\,\mathrm{kg/m^3}$、$v = 18\,\mathrm{km/s}$：

$$ m = \frac{4}{3}\pi (25)^3 \times 3000 \approx 1.96\times10^8\,\mathrm{kg} $$

$$ E_k = \frac{1}{2}\times 1.96\times10^8 \times (1.8\times10^4)^2 \approx 3.2\times10^{16}\,\mathrm{J} $$

換算成 TNT 當量（1 噸 TNT $= 4.184\times10^9\,\mathrm{J}$）：

$$ \frac{3.2\times10^{16}}{4.184\times10^9} \approx 7.6\times10^{6}\,\mathrm{tonne} \approx 7.6\,\mathrm{Mt} $$

約 760 萬噸 TNT，相當於數百顆廣島原子彈——而這只是一顆 50 公尺的小石頭。1908 年通古斯事件（Tunguska）的能量量級正在此範圍。這就是為什麼動能撞擊器（如 2022 年 NASA 的 DART 任務撞擊雙星小行星迪莫弗斯 Dimorphos）能成為偏轉小行星的可行手段：不需炸藥，純靠動量轉移就足以改變一顆小行星的軌道。

重點回顧

• 非引力力讓彗星自帶火箭引擎：自轉 + 熱滯後使昇華噴氣偏離正午方向，橫向分量 $A_2$ 系統性地改變週期，這就是哈雷彗星「遲到」的物理。
• 亞爾科夫斯基效應讓陽光能在百萬年尺度搬動小行星軌道，是主帶補貨近地小行星的關鍵；其旋轉版 YORP 效應能改變自轉甚至使小行星解體。
• 共振與混沌塑造了 Kirkwood gaps：與木星的平均運動共振（如 3:1 在 2.50 AU）激發偏心率，把小行星送進內太陽系，過程由混沌動力學主導。
• 矮行星的球形源於重力壓垮材料強度，冰質天體臨界半徑約 270 公里、岩質更大，解釋了「冰天體更早變圓」。
• 高速撞擊把動能轉為激震波，一顆 50 公尺小行星即達數 Mt 當量；DART 證明純動量轉移可偏轉小行星。

深入探討（研究所視角）

若你想再往深處走，以下是幾條當代研究的前沿線索：

1. 昇華驅動力的物理建模。 Marsden 經驗公式雖好用，但它把彗核當成均勻黑箱。現代研究改用熱物理模型（thermophysical model），在彗核表面網格上解一維熱傳導方程式 $\rho c_p \partial T/\partial t = \partial(\kappa\,\partial T/\partial x)/\partial x$，配合各冰種（$\mathrm{H_2O}$、$\mathrm{CO}$、$\mathrm{CO_2}$）的克勞修斯–克拉佩龍昇華曲線，逐網格算出噴氣通量與動量，再積分得到全核淨推力。Rosetta 任務對 67P/Churyumov–Gerasimenko 的就地觀測，讓這類模型首次能用真實的形狀、孔隙率與活躍區分布來檢驗，發現噴氣有強烈的季節性與緯度依賴。

2. YORP 與雙小行星形成。 YORP 力矩 $\tau_{\rm YORP} \propto L_\odot R^3/(c\,r^2)$ 對形狀極其敏感，理論上能把小行星自轉週期推到自轉解體極限（赤道離心加速度等於自重力時，週期約 2.2 小時）。一旦逼近此極限，鬆散堆積體（rubble pile）會從赤道甩出物質，形成衛星或雙系統——龍宮（Ryugu）與貝努那種「陀螺形」赤道隆起正是此機制的化石證據。這把熱物理、自轉動力學與顆粒物質力學耦合在一起，是目前活躍的交叉領域。

3. 共振重疊與 Chirikov 判據。 為什麼共振區是混沌的？Chirikov 共振重疊判據（resonance overlap criterion）給出定量答案：當相鄰共振的「寬度」在相空間中重疊，規則軌道就崩解為混沌海。這把小行星帶的結構、行星遷移（如 Nice 模型中海王星向外遷移掃過共振）與系外行星系統的穩定性，統一在同一套哈密頓微擾理論之下。

4. 撞擊的標度律與行星防禦。 DART 撞擊迪莫弗斯後，其公轉週期縮短了約 33 分鐘，遠超純動量守恆的預期——因為撞擊噴出的碎屑（ejecta）反衝提供了額外動量，效率因子 $\beta$ 達到約 3–4。這個 $\beta$ 取決於目標的孔隙率與內聚力，是當前撞擊標度律（impact scaling laws）研究的核心參數，也直接決定未來真正要偏轉一顆威脅小行星時，需要多大的撞擊器。這裡，小天體科學從純粹的好奇，變成了攸關文明存續的工程問題。

這些線索共同指向一個訊息：太陽系的小天體不是冰冷靜止的化石，而是被陽光、自轉、共振與碰撞持續雕琢的動態系統。理解它們的機制，既是讀懂太陽系四十六億年歷史的鑰匙，也是守護地球未來的科學基礎。