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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

參考文獻

[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