科學(xué)通報(bào) | 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ìn)而吸附大量氣體. 然而, 熱木星卻運(yùn)行在緊貼恒星的熾熱軌道上, 這意味著它們必然在形成后經(jīng)歷了軌道遷移, 從遙遠(yuǎn)的寒冷地帶“搬遷”至此.

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

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

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

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

在恒星潮汐長期持續(xù)的作用下, 行星的軌道角動量會不斷耗散, 導(dǎo)致軌道逐漸收縮, 最終可能被恒星吞噬. 這一過程的快慢取決于一個關(guān)鍵物理參數(shù)——恒星的潮汐品質(zhì)因子( Q ). Q 值越小, 表示潮汐耗散越有效, 軌道衰減越快. 然而, 學(xué)界對 Q 的估計(jì)存在巨大分歧, 不同研究給出的數(shù)值跨越多達(dá)5個數(shù)量級 [5] . 一個突出的矛盾是: 目前唯一被確認(rèn)為軌道正在收縮的熱木星WASP-12b [6] , 其觀測數(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)前觀測探測到. 這一矛盾凸顯了單一案例與群體統(tǒng)計(jì)之間的脫節(jié).

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

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

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

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

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

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

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

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

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

參考文獻(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