科學通報 | LAMOST助力解開熱木星“從何而來”又“向何而去”之謎

1995年, 飛馬座51b的發現 [1] 徹底打破了人類對行星系統的固有認知. 這顆質量與木星相當, 卻在短短4.2天內就繞母星旋轉一周的“熱木星”, 不僅為其發現者贏得了諾貝爾獎, 也給天文學家提出了一個持續三十年的核心挑戰: 這類“叛逆”的巨行星究竟從何而來?

依據基于太陽系建立的行星形成經典理論(如核吸積模型), 氣態巨行星理應在原行星盤外圍的“雪線”之外形成, 那里豐富的冰質物質才能堆積出足夠大的核心, 進而吸附大量氣體. 然而, 熱木星卻運行在緊貼恒星的熾熱軌道上, 這意味著它們必然在形成后經歷了軌道遷移, 從遙遠的寒冷地帶“搬遷”至此.

為解釋這一遷徙過程, 學者提出了多種競爭性的物理機制 [2] : 有的認為它們在行星系統的“嬰兒期”通過原行星盤遷移(disk migration), 在數百萬年內即完成 [3] ; 有的認為它們在“青春期”通過行星間的引力散射(planet-planet scattering)或Kozai-Lidov機制, 在幾千萬到幾億年之內到達; 還有觀點認為, 它們可能經歷了一種歷時數億甚至數十億年的“長期混沌(secular chaos)”過程 [4] , 才步履蹣跚地抵達現在的軌道. 然而, 何種機制占主導? 不同機制的定量貢獻如何? 這些關鍵問題長期懸而未決.

熱木星的奧秘不僅在于其“從何而來”的起源, 還關乎其“向何而去”的最終命運. 由于距離宿主恒星極近, 熱木星與恒星之間會發生強烈的潮汐相互作用. 這種潮汐作用主要體現為兩個方面: 一是恒星對行星內部造成的潮汐, 即行星潮汐; 二是行星對恒星內部造成的潮汐, 即恒星潮汐.

行星潮汐會推動行星的自轉與公轉逐漸趨于同步, 而恒星潮汐則促使恒星的自轉與其公轉同步. 由于行星質量相對較小, 其軌道同步通常在演化早期即可實現; 而恒星質量極大, 其軌道同步往往難以達成, 因此恒星潮汐會一直持續作用. 在熱木星系統長期演化的大部分時間里, 恒星潮汐占據主導地位.

在典型情況下, 恒星的自轉速度慢于行星的公轉速度. 因此, 恒星潮汐的主要效應表現為: 行星公轉“拖拽”著恒星自轉加速, 使恒星獲得角動量; 與此同時, 行星公轉逐漸失去角動量, 導致其軌道持續收縮. 這一機制與地月系統的潮汐作用原理相同, 但情況相反: 地球自轉速度快于月球公轉, 地球自轉減速并將角動量轉移給月球, 使得月球公轉軌道擴張, 月球逐漸遠離地球.

在恒星潮汐長期持續的作用下, 行星的軌道角動量會不斷耗散, 導致軌道逐漸收縮, 最終可能被恒星吞噬. 這一過程的快慢取決于一個關鍵物理參數——恒星的潮汐品質因子( Q ). Q 值越小, 表示潮汐耗散越有效, 軌道衰減越快. 然而, 學界對 Q 的估計存在巨大分歧, 不同研究給出的數值跨越多達5個數量級 [5] . 一個突出的矛盾是: 目前唯一被確認為軌道正在收縮的熱木星WASP-12b [6] , 其觀測數據暗示log10( Q )≈5.1; 而先前基于熱木星整體分布統計的研究則傾向于log10( Q )=6~9. 前者預言應有較多正處于收縮中的熱木星, 與當前僅發現一例的事實不符; 后者則預言軌道收縮極慢, 難以被當前觀測探測到. 這一矛盾凸顯了單一案例與群體統計之間的脫節.

無論是追溯起源還是預測最終演化命運, 時間(即宿主恒星的年齡)都是最關鍵的因素. 不同遷移機制作用于迥異的時間尺度, 而恒星潮汐作用的效率則決定了熱木星“墜落”的時標. 因此, 精確繪制熱木星出現頻率隨宿主恒星年齡演變的圖譜, 成為破解雙重謎題的關鍵.

然而, 精確測定單顆恒星的年齡本身就是一個難題. 傳統恒星演化模型對系外行星宿主恒星給出的年齡估計往往誤差較大(中值誤差可達50% [7] ). 我們另辟蹊徑, 采用恒星運動學方法來估算恒星群體的平均年齡. 其原理在于, 恒星在銀河系引力勢中的運動速度彌散會隨該星群年齡的增長而增加. 依托我國大天區面積多目標光纖光譜天文望遠鏡(LAMOST, 也稱郭守敬望遠鏡)提供的大規模、均勻恒星光譜數據, 結合Gaia衛星的高精度天體測量信息, 我們成功重新標定了速度彌散與年齡的關系, 將恒星群體年齡的估算誤差顯著降低至約15% [8] , 為在時間維度上研究行星演化奠定了堅實基礎.

基于我們新的年齡估計, 我們在2023年的工作 [9] 中發現熱木星的數量隨恒星變老而整體減少, 并成功解釋了觀測上的兩個長期謎題: (1) 視向速度法樣本和凌星法樣本中熱木星出現率不一致; (2) 1999年哈勃空間望遠鏡未能在球狀星團中發現熱木星. 這兩個謎題的根源均在于未充分考慮熱木星出現率隨恒星年齡演變的效應. 具體而言, 第一個謎題是由于兩種樣本中的恒星年齡存在顯著系統差異所致; 在扣除該年齡系統差后, 兩種方法得到的熱木星出現率趨于一致. 第二個謎題則源于1999年預測球狀星團中熱木星數量(預估為17顆)時, 未考慮熱木星會因球狀星團中恒星的極高年齡(通常超過百億年)而顯著減少. 我們基于該年齡效應重新評估了當年的觀測數據, 預計可探測熱木星約為0.8顆, 與當時實際“零發現”的結果吻合.

2025年這次最新研究 [10] , 我們更進一步: 從“定性”走向“定量”, 從“描述現象”到“鎖定理論”. 通過對LAMOST和Gaia等數據的分析, 我們繪制出了更精確的“熱木星年齡圖譜”. 一個突破性的發現是: 熱木星數量衰減并非勻速, 而是在恒星年齡約20億年處出現一個明顯的“拐點”: 前期減少得慢, 后期減少得快( 圖1(a ) ). 這揭示了熱木星家族由兩大“族群”構成. “早期到達者”在恒星誕生后數千萬年內就迅速就位( 圖1(b) ). “晚期到達者”不緊不慢, 在數千萬年乃至數十億年的時間里才陸續就位( 圖1(c) ). 這兩大族群共同作用, 塑造了先緩后急的“分段”衰減曲線.

圖 1 熱木星出現率的分段衰減與多機制形成示意圖. (a) 總樣本的演化規律; (b) 早期到達族群; (c) 晚期到達族群; (d) 兩個族群軌道構型示意. 詳見文獻 [10]

“晚期到達”族群還展現出與“長期混沌”機制高度吻合的特征: 它們更傾向于環繞金屬豐度較高的恒星, 且其軌道面與恒星自轉軸往往存在較大夾角. 這為“長期混沌”是其主要形成通道提供了強有力的旁證.

通過將觀測圖譜與包含“早期”和“晚期”形成機制的混合物理模型進行精細擬合, 我們首次實現了對熱木星起源的定量分解: 約62%的熱木星通過盤遷移、散射等較快機制形成(“早期到達者”), 而約38%則源自“長期混沌”這一慢速通道(“晚期到達者”). 同時, 我們首次從群體演化角度, 將類太陽恒星的潮汐品質因子約束在log10( Q )=5.4~6.1. 基于這一結果所預測的、當前應處于可探測軌道衰減階段的熱木星數量, 也與實際觀測到的候選體數目吻合, 從而彌合了以往研究中的矛盾.

該研究是“行星的空間分布和年齡演化”(Planets Across Space and Time, PAST)系列研究的第七項核心成果(PAST-7). 該計劃旨在繪制行星系統的“時空演化圖譜”, 追溯不同行星族群的完整“生命史” [ 11 , 12 ] . 該工作通過對“年齡”這一關鍵維度的深度挖掘, 首次構建了一個能自洽統一解釋熱木星“起源-遷移-潮汐演化”全過程的物理框架.

展望未來, 隨著TESS、GAIA、PLATO以及我國系外地球2.0(ET 2.0)搜尋等項目 [ 13 ~ 15 ] 的推進, 系外行星研究將邁入“十萬級樣本”時代. 這不僅能讓人們對熱木星進行更多維、更系統的研究, 還有望揭示更多行星族群(如超級地球、亞海王星)的形成與演化歷史. 通過將日益豐富的系外行星統計數據與我們的太陽系進行深度比較, 人類終將更清晰地回答那個根本性問題: 我們的太陽系在宇宙中是否特殊? 如果特殊, 它又是如何形成并演變成今天這般模樣的.

參考文獻

[1] Mayor M, Queloz D. A Jupiter-mass companion to a solar-type star . Nature , 1995 , 378: 355 -359

[2] Dawson R I, Johnson J A. Origins of hot Jupiters . Annu Rev Astron Astrophys , 2018 , 56: 175 -221

[3] Lin D N C, Bodenheimer P, Richardson D C. Orbital migration of the planetary companion of 51 Pegasi to its present location . Nature , 1996 , 380: 606 -607

[4] Wu Y, Lithwick Y. Secular chaos and the production of hot Jupiters . Astrophys J , 2011 , 735: 109

[5] Ogilvie G I. Tidal dissipation in stars and giant planets . Annu Rev Astron Astrophys , 2014 , 52: 171 -210

[6] Yee S W, Winn J N, Knutson H A, et al. The orbit of WASP-12b is decaying . Astronphys J Lett , 2020 , 888: L5

[7] Berger T A, Huber D, van Saders J L, et al. The Gaia–Kepler stellar properties catalog. I. Homogeneous fundamental properties for 186,301 Kepler stars . Astron J , 2020 , 159: 280

[8] Chen D C, Xie J W, Zhou J L, et al. Planets Across Space and Time (PAST). I. Characterizing the memberships of Galactic components and stellar ages: revisiting the kinematic methods and applying to planet host stars . Astrophys J , 2021 , 909: 115

[9] Chen D C, Xie J W, Zhou J L, et al. The evolution of hot Jupiters revealed by the age distribution of their host stars . Proc Natl Acad Sci USA , 2023 , 120: e2304179120

[10] Chen D C, Xie J W, Zhou J L, et al. The origin and tidal evolution of hot Jupiters constrained by a broken age–frequency relation . Nat Astron , 2026 , 10: 92 -104

[11] Yang J Y, Chen D C, Xie J W, et al. Planets Across Space and Time (PAST). IV. The occurrence and architecture of Kepler planetary systems as a function of kinematic age revealed by the LAMOST–Gaia–Kepler sample . Astron J , 2023 , 166: 243

[12] Tu P W, Xie J W, Chen D C, et al. Age dependence of the occurrence and architecture of ultra-short-period planet systems . Nat Astron , 2025 , 9: 995 -1006

[13] Kunimoto M, Winn J, Ricker G R, et al. Predicting the exoplanet yield of the TESS prime and extended missions through years 1–7 . Astron J , 2022 , 163: 290

[14] Rauer H, Aerts C, Cabrera J, et al. The PLATO mission . Exp Astron , 2025 , 59: 26

[15] Ge J, Chen W, Chen Y, et al. Search for a Second Earth – the Earth 2.0 (ET) Space Mission (in Chinese) . Chin J Space Sci , 2024 , 44: 400 -424