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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

參考文獻(xiàn)

[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