跳转到内容

凌星法

凌星法系外行星探测开普勒TESS光变曲线

概述

凌星法（transit method，又称凌星光度测量法）是目前发现系外行星最高效的方法，通过精密测量恒星亮度随时间的变化，探测行星从恒星圆面前方经过时造成的周期性亮度下降来推断行星的存在和大小。

该方法的逻辑简洁而优雅：一颗行星挡住宿主恒星一小部分光线，精度足够的光度计就能探测到这个"凹陷"。开普勒空间望远镜（2009–2018年）利用此方法在约15万颗恒星中发现了超过2600颗确认的系外行星，以及数千个候选体，彻底改变了系外行星科学的面貌。截至2024年，所有系外行星探测方法中，凌星法贡献了约四分之三（约75%）的已确认行星，是迄今最高产的探测方法（NASA系外行星档案数据）。

基本原理

凌星光变曲线

当行星从恒星圆面前通过时，恒星接收到的光通量下降，下降幅度（凌星深度，transit depth）由行星与恒星半径之比决定：

δ=(RpR∗)2\delta = \left(\frac{R_p}{R_*}\right)^2

例如，木星半径约为太阳的1/10，因此热木星造成的光变约1%；地球大小的行星造成的下降约0.01%（100 ppm）——这需要极高精度的光度测量，是开普勒卫星的关键设计指标（光度精度约20 ppm/6小时积分）。

凌星持续时间

凌星持续时间 TdurT_{dur} 取决于轨道周期、恒星半径和轨道倾角：

Tdur≈R∗πaPT_{dur} \approx \frac{R_*}{\pi a} P

（简化中心凌星情形），其中 $a$ 为轨道半长轴，$P$ 为轨道周期。典型热木星凌星持续2–4小时，地球尺寸的类地行星在宜居带凌星持续约13小时（以太阳-地球参数为例）。

几何概率

一颗行星正好从恒星圆面前通过需要特定的轨道倾角。凌星概率：

ptransit≈R∗ap_{transit} \approx \frac{R_*}{a}

对于0.05 AU的热木星：p≈10%p \approx 10\%；对于1 AU的类地行星：p≈0.5%p \approx 0.5\%；对于宜居带的行星绕K型恒星：约1%。这意味着凌星法存在强烈的几何选择效应，偏向于发现近星轨道和大半径行星。

凌星法能告诉我们什么

从凌星光变曲线中，天文学家可以提取：

  • 行星半径（RpR_p）：直接由凌星深度 δ\delta 和恒星半径 R∗R_*（由恒星物理模型确定）得出
  • 轨道周期（$P$）：由相邻凌星的时间间隔确定
  • 轨道倾角（$i$）：由凌星的"盒形"形状（入射角/出射角斜率）限定
  • 轨道半长轴（$a$）：由开普勒第三定律 a3∝P2M∗a^3 \propto P^2 M_* 结合恒星质量
  • 大气成分（传输光谱）：在不同波长下凌星深度不同，反映大气中各分子的吸收截面

凌星法本身不能直接给出行星质量——需要配合视向速度法（RV）或凌星时变（TTV）方法才能确定密度和成分。

凌星法的数学细节：超越简单模型

实际的凌星光变曲线建模涉及多种物理效应：

临边昏暗（limb darkening）：恒星并非均匀发光圆盘——中心更亮，边缘（临边）因视线穿过更厚大气层而较暗。精确的凌星拟合需要用二次或更高阶的临边昏暗模型：

I(μ)=I0[1−u1(1−μ)−u2(1−μ)2]I(\mu) = I_0 [1 - u_1(1-\mu) - u_2(1-\mu)^2]

其中 μ=cos⁡θ\mu = \cos\theta 为视线方向与法向的余弦，u1u_1、u2u_2 为临边昏暗系数（依赖恒星温度和观测波段）。若不正确处理临边昏暗，会导致行星半径被低估（甚至偏差达几个百分比）。

影响行星半径测量的系统误差：

  • 恒星半径不确定性（通常5–10%）是行星半径不确定性的主要来源，通过盖亚视差+恒星光谱测定可改善
  • 非对称凌星形状（egress/ingress不对称）暗示行星大气层延伸或恒星黑子交叉
  • 多行星系统中的引力摄动导致凌星时刻偏移（TTV），需要同时拟合所有行星参数

传输光谱：大气的化学指纹

当行星处于凌星位置时，部分恒星光穿过行星大气的临边，被大气中的分子吸收，留下波长相关的"透射光谱"。这一技术——传输光谱学（transmission spectroscopy）——是表征系外行星大气的核心工具：

  • 哈勃WFC3在HD 209458 b、WASP-39 b等热木星大气中探测到水汽（H₂O）
  • JWST 2022年的早期释放科学结果在WASP-39 b中首次明确探测到CO₂、SO₂（光化学产物）、Na、K
  • 在TRAPPIST-1系统的岩石行星上，JWST正在尝试探测是否存在大气（TRAPPIST-1 b的结果暗示没有明显大气）

主要空间任务

任务机构运行时间特点
开普勒（Kepler）NASA2009–2018约15万颗目标星，4年连续观测，发现统计规律
K2NASA2014–2018开普勒后续，黄道面观测，多视场
TESSNASA2018–至今全天巡测，主攻周期短于27天行星，重点离地球近的亮星
CoRoTESA/CNES2006–2013首批空间凌星探测，发现CoRoT-7 b（首颗岩石系外行星候选）
CHEOPSESA2019–至今定点精密测量已知系外行星系统，获取精确半径
PLATOESA计划2026寻找类日恒星宜居带类地行星，光度精度目标极高

误判与假阳性

凌星法最大的挑战是假阳性（false positives）：许多凌星信号并非来自行星，而是：

  • 背景双星：行星圆面后方恰好有一对食双星，其周期性遮掩模拟行星凌星
  • 近处双星（grazing eclipsing binary）：双星擦边凌星的光变深度可能与行星相当
  • 双星绕主星（hierarchical triple）：远处双星绕近处恒星运行

统计上，约10%的开普勒候选体最终为假阳性。因此大多数行星需要地面望远镜的视向速度跟踪（确认行星质量量级）或高分辨率成像（排除背景双星）才能"确认"。

统计发现：宇宙行星的普遍性

凌星法最深远的贡献不是某颗具体行星，而是对行星频率（occurrence rate）的大统计揭示。开普勒数据分析显示：

  • 约90%以上的类FGK恒星在轨道周期小于400天内至少有一颗行星
  • 超级地球和迷你海王星（1–4 R⊕R_\oplus，周期<100天）的出现率约为每颗恒星0.5–1颗
  • 地球大小行星（0.8–1.25 R⊕R_\oplus）在宜居带内的出现率约为0.1–0.4（即每5–10颗类日恒星有一颗宜居带类地行星候选）
  • M型矮星周围的岩石行星出现率更高，约为类日恒星的2–3倍

这意味着银河系中潜在类地宜居带行星数量在数十亿颗量级——这是迄今系外行星科学给出的最具震撼力的统计结论。

次级凌星与相位曲线

除了行星从恒星前方通过的"主凌星"外，当行星移动到恒星背后时发生次级凌星（secondary eclipse）。通过测量次级凌星前后恒星亮度的微小变化，可直接测量行星的热辐射（在红外波段）或反射率（在可见光波段）。

结合整个轨道上恒星+行星的亮度随时间变化（相位曲线，phase curve），可以反演行星的东西侧温度分布图，即所谓的行星天气图。斯皮策在多颗热木星（如HD 189733 b、WASP-43 b）上绘制的相位曲线是了解其大气环流的核心依据；JWST正将这一技术推广到更小、更冷的行星。

开普勒遗产与统计宇宙观

开普勒任务（2009–2018年）是系外行星科学的转折点。除了发现数千颗具体行星外，其最深远的贡献是建立了行星统计学：

开普勒测量了约15万颗恒星，对其中大约2600颗行星进行了统计学意义上的"纠正"（corrected for detection efficiency），得到了无偏的行星频率分布：

  • 小于约4 R⊕R_\oplus 的行星在轨道周期<200天内的每颗恒星平均出现率约1–2颗
  • 行星在1–2 R⊕R_\oplus 和 2–4 R⊕R_\oplus 两个峰之间存在富尔顿间隙（见超级地球与迷你海王星条目）
  • 类地行星（0.8–1.25 R⊕R_\oplus）在宜居带内的出现率估计约为每7–10颗FGK型恒星1颗——以银河系约10¹⁰颗FGK恒星计算，银河系中可能有约10⁹颗宜居带类地行星候选

这些数字深刻改变了天文学界的集体认知：行星不是异常，而是恒星演化的自然产物；地球可能不是宇宙中唯一的岩石行星，甚至可能不是罕见的。

跨域连接

  • 赫罗图：凌星深度只给出半径比 Rp/R∗R_p/R_*，所以行星半径的误差从来不会比恒星半径的误差更小。恒星半径由光度与有效温度经赫罗图定位反推，正文指出这项百分之几到一成的不确定度正是行星半径误差的主要来源。推论很直接：改进恒星参数比延长测光时间更能压低行星半径误差。
  • 概率：轨道朝向若各向同性，能看到凌星的概率约为 R∗/aR_*/a——这是球面上的一个测度，而不是技术缺陷。由此得到本方法的偏倚：短周期、大半径的行星被系统性放大。要从被探测样本反推真实出现率必须除回这个概率，而概率越小除法的方差越大，一个天文单位处的地球类似物因此始终难以约束。
  • 机器学习：正文给出约一成开普勒候选体最终是假阳性，分类器要做的就是把盒形凌星与擦边食双星的 V 形凹陷分开。但训练集来自已确认行星，而确认流程本身偏向亮星与短周期，分类器会继承这层偏倚。所以判断产出质量要看注入-回收测试，而不能只看分类准确率。
  • 光谱学：不同波长上凌星深度不同，是因为大气分子的吸收截面不同——这是把测光变成化学分析的那一步。可检验的推论是：同一分子的多个吸收带必须同时出现，相对强度还受截面比约束；只在单一波段冒出来的"信号"更可能是仪器系统误差。
  • 万有引力与开普勒定律：周期本身只是钟表读数，要变成轨道距离必须用 a3∝P2M∗a^3 \propto P^2 M_* 换算。这一步把"任务持续多久"直接翻译成"能触及多远的轨道"：至少要看到两三次凌星，周期上限就被观测跨度锁死，于是宽轨道行星不是难探测，而是被时间排除在外。

参考文献

  • Charbonneau, D. et al. (2000). "Detection of Planetary Transits Across a Sun-like Star." The Astrophysical Journal Letters, 529, L45–L48.（首次观测到系外行星凌星HD 209458 b）
  • Borucki, W. J. et al. (2010). "Kepler Planet-Detection Mission: Introduction and First Results." Science, 327, 977–980.（开普勒任务正式结果）
  • Seager, S., & Mallén-Ornelas, G. (2003). "A Unique Solution of Planet and Star Parameters from an Extrasolar Planet Transit Light Curve." The Astrophysical Journal, 585, 1038–1055.（凌星参数分析的理论基础）
  • NASA系外行星档案可实时访问：exoplanetarchive.ipac.caltech.edu

数据来源：NASA系外行星档案（2024）；开普勒任务DR25数据集（MAST）；TESS任务数据。

正文中引用本条 · 1