白矮星行星大气中的气溶胶与碳氢化合物

《Nature》：Aerosols and hydrocarbons in the atmosphere of a white dwarf planet

　　

大多数恒星，包括太阳，终将在未来演化为红巨星，随后成为白矮星。近期已有多颗候选行星被确认围绕白矮星运行，表明行星能够在恒星主序后阶段完整存活。然而，关于主序后行星的大气成分知之甚少，迄今大气探测最为演化的凌星行星仅围绕亚巨星运行。本研究报告了对白矮星行星WD 1856 b的大气探测，该探测通过詹姆斯·韦伯空间望远镜（JWST)近红外光谱仪（NIRSpec)棱镜模式（PRISM）的透射光谱实现。0.5–5.0 μm光谱揭示了碳氢化合物（几率比167:1–5,377:1，其中CH₄以17:1–30:1为优）、气溶胶（2 × 10⁵:1–2 × 10⁶:1）以及行星夜侧热辐射（2 × 10⁶³:1–2 × 10⁷³:1）的存在。光谱分析将WD 1856 b的质量约束为4.3–10.9 M_J，发现其大气富含碳（CH₄丰度约7%），有效温度超过预期行星平衡温度（390–412 K 对比 160 K）。基于冷却模型，这些结果表明WD 1856 b在白矮星阶段3.0–5.5 Gyr时经历了一次与迁移相关的再加热事件，这与主序后潮汐演化为现今0.02 au圆轨道的演化图景一致。研究结果为了解围绕类太阳质量恒星运行的巨行星的最终命运提供了一个窗口。

**研究背景与意义**

大多数恒星，包括太阳，终将演化为红巨星并最终成为白矮星。近年来，多颗候选行星被确认围绕白矮星运行，证明行星能够在恒星主序后阶段完整存活。然而，关于主序后行星的大气成分，科学界了解甚少。迄今大气探测最为演化的凌星行星仅围绕亚巨星运行，尚未有围绕白矮星行星的系统性大气研究。白矮星行星系统为理解行星系统的最终演化命运提供了独特实验室，但受限于宿主恒星亮度较低（如WD 1856的J星等为15.677），高精度光谱观测面临挑战。此外，WD 1856 b的近距离轨道（0.02 au/1.4天）需要其在恒星红巨星阶段后发生轨道演化以避免被吞噬，但具体的演化机制——是公共包层演化还是高偏心率迁移——长期存在争议，而行星质量的不确定性是区分这两种情景的关键瓶颈。

研究人员利用JWST NIRSpec PRISM对白矮星行星WD 1856 b开展透射光谱观测，首次实现了对白矮星行星大气的详细探测。该研究发表于《Nature》，不仅为理解类太阳质量恒星系统中巨行星的最终命运提供了直接观测证据，还揭示了主序后行星大气成分、热演化及轨道迁移机制的深层物理图像。

**关键技术方法**

样本来源于白矮星行星WD 1856 b的JWST第2358号客座观测程序（PI: Ryan MacDonald），观测时长1.98小时，覆盖0.6–5.5 μm波段，平均分辨本领λ/Δλ ≈ 100。数据处理采用两种独立管线：FIREFLy和Juniper，分别提取光谱凌星光变曲线并进行拟合。光谱建模使用辐射传输与反演代码POSEIDON，针对掠入射凌星几何和夜侧热辐射效应进行定制开发。反演分析涵盖多种碳、氧、氮、硫、磷载体气体及不平衡示踪分子，采用贝叶斯框架（PyMultiNest）进行大气参数反演，并通过Mie散射模型探索气溶胶特性。热历史重建基于ATMO2020和Sonora Bobcat次恒星天体冷却模型，结合MIST演化模型和MESA白矮星冷却模型约束系统年龄。

**研究结果**

**白矮星行星WD 1856 b的大气探测**：研究人员通过JWST NIRSpec PRISM在2023年4月27日观测了WD 1856 b的一次凌星事件。两种独立数据管线均获得高信噪比光谱凌星光变曲线，发现在4.19–4.96 μm近红外波段比0.55–1.71 μm可见光波段的凌星深度浅3%，这一差异的物理机制被确认为行星夜侧热辐射对凌星深度的稀释效应。

**碳氢化合物的探测**：反演分析确认WD 1856 b大气中至少存在一种碳氢化合物（Juniper数据4.5σ/ln? = 8.59；FIREFLy数据3.6σ/ln? = 5.12）。主要贡献者为CH₄（Juniper数据3.1σ/ln? = 3.40；FIREFLy数据2.9σ/ln? = 2.87），C₂H₆也对最佳拟合模型有贡献。CH₄在1.75 μm、2.3 μm和3.3 μm处产生三条显著吸收带，反演得到的CH₄丰度为2–20%，对应约100×太阳金属丰度的大气金属增丰（log(C/H) = 2.24^+0.33_-0.40，FIREFLy数据）。掩蔽3.1–3.5 μm数据后的反演仍给出3.0σ的模型偏好，证实了CH₄探测的可靠性。最佳拟合模型中4.3 μm处存在PH₃的潜在贡献，可能指示类似木星上层大气的不平衡垂直混合，但受限于信噪比和模型简并性尚不能确认。3.4 μm处可归因于C₂H₆的特征与3.3 μm附近CH₄吸收部分简并，因此C₂H₆未被显著探测。综合所有碳氢化合物（CH₄、C₂H₂、C₂H₄、C₂H₆）的证据在统计上达到约4σ显著性。NH₃、H₂O等其他化学组分未被探测到。

**气溶胶与夜侧热辐射的探测**：除碳氢化合物外，研究还探测到来自观测者朝向夜侧半球的热辐射（Juniper数据18.匪夷所思的是，研究人员在分析中发现了5σ/ln? = 14.5；FIREFLy数据17.3σ/ln? = 146）和气溶胶（Juniper数据5.7σ/ln? = 14.5；FIREFLy数据5.3σ/ln? = 12.3）。夜侧半球约400 K的热辐射导致光谱在3.5 μm以上出现显著的向下倾斜。由于透明窗口中的热辐射来自更深、更热的大气层，反演要求约10 mbar以下存在光学厚的云 deck以解释CH₄吸收带之间缺乏"负尖峰"的现象。大多数热辐射源自100 mbar附近的云顶温度，这与MIRI测光在近四分相附近测得约190 K热过剩形成潜在对比，可能指示不同半球或时间变异的云顶压力。气溶胶的额外证据来自1 μm以下比纯净H₂瑞利散射更显著的散射斜率。Mie散射模型未能确定气溶胶的具体化学成分，但水富有机ExoHaze（400K变体）、KCl和NH₃冰分别作为最佳拟合的霾/烟炱、T-Y型褐矮星云和冰类气溶胶被识别。

**行星质量与温度约束**：研究首次约束了WD 1856 b的质量为4.3–10.9 M_J（两类数据降采样1σ范围，或FIREFLy数据6.7^+2.8_-2.4 M_J，Juniper数据7.8^+3.1_-2.7 M_J）。该质量约束源于CH₄吸收特征与夜侧大气发射压力积分光学深度之间的标高权衡。行星有效温度（T_eff = 390–412 K，或FIREFLy数据400⁺⁶_-10 K，Juniper数据405⁺⁷_-11 K）显著高于平衡温度（160 K）和约10 Gyr系统年龄下演化巨行星预期温度（?100 K）。

**热演化历史与轨道迁移**：基于次恒星天体冷却模型的热历史重建显示，WD 1856 b的再加热发生在渐近巨星支（AGB）阶段结束后3.0–5.5 Gyr（FIREFLy数据4.2^+1.0_-1.2 Gyr，Juniper数据4.6^+0.9_-1.0 Gyr）。更保守的2σ下限表明再加热至少在AGB阶段后1.4 Gyr（FIREFLy）或2.1 Gyr（Juniper）发生。AGB阶段和后期AGB/前白矮星阶段分别仅持续约2 Myr和0.1 Myr，因此再加热时间不与公共包层演化相符。这一结论对于FIREFLy数据具有>2σ显著性，对于Juniper数据具有>3σ显著性。WD 1856 b最可能经历了高偏心率迁移至当前轨道，推断的再加热事件对应潮汐圆化过程。

**讨论与结论**

WD 1856 b作为首颗被详细表征的白矮星凌星行星，其推断的热演化历史表明高偏心率迁移是巨行星在恒星主序后的合理命运。反演得到的CH₄丰度与海王星深层大气（4%）相当，这需要行星H₂包层从挥发性富集物质中获得显著碳增丰，无论是迁移前吸积还是迁移后吸积。高大气金属度（?100×太阳）增强了气溶胶产生，与透射光谱中短波长散射斜率的探测一致。

研究人员排除了多种替代解释：纯被动冷却在约10 Gyr系统年龄下需要约24 M_J的质量，与观测矛盾；氘核聚变虽可能在2σ质量范围内发生，但其持续时间远短于系统年龄且预测光度远高于观测值；潮汐加热虽在原理上可维持观测有效温度，但对应的轨道偏心率阻尼时标仅约0.075 Gyr，使研究人员处于高偏心率迁移最后阶段的概率极低；白矮星辐照的影响亦被评估为次要因素（? ≈ 25）。

WD 1856 b的观测证明，围绕白矮星的行星光谱学为确定恒星死亡后行星系统的命运提供了全新机遇，开启了理解行星系统终极演化的观测窗口。