视向速度法

概述#

视向速度法（radial velocity method，又称多普勒测速法）是最早成功探测系外行星的方法，也是目前唯一能在不依赖凌星的情况下直接测量行星质量的常规手段。其原理基于牛顿力学中最基本的事实：一颗行星并非绕静止恒星运转，而是行星与恒星共同绕系统质量中心（重心）运动。当恒星随行星拉拽而摇摆时，其径向速度（沿视线方向的速度分量）会发生周期性变化，从而导致恒星光谱谱线发生多普勒频移。

这个频移的幅度极小。太阳因木星而产生的最大径向速度变化仅约12 m/s，相当于每小时43公里——一个勤快的自行车手的速度。如果是地球引起的，太阳的速度变化只有约0.1 m/s。探测如此微小的速度变化，需要分辨率极高的光谱仪和极其稳定的观测环境，这项技术的发展经历了数十年的仪器革命。

物理原理#

多普勒频移#

当光源以径向速度 $v_r$ 运动时，观测到的波长 $\lambda_{obs}$ 与静止波长 $\lambda_0$ 之间的关系：

$$\frac{\Delta\lambda}{\lambda_0} = \frac{v_r}{c}$$

对于10 m/s的速度变化，可见光（约550 nm）的波长偏移约 $1.8 \times 10^{-5}$ nm，即约0.00002纳米。探测这样的偏移需要光谱仪的波长分辨率达到约0.001 nm量级，并且长达数月到数年的漂移控制在亚毫秒量级。

径向速度半幅#

由轨道力学，行星对恒星的径向速度扰动半幅（$K$）为：

$$K = \left(\frac{2\pi G}{P}\right)^{1/3} \frac{M_p \sin i}{(M_* + M_p)^{2/3}} \frac{1}{\sqrt{1-e^2}}$$

对低质量行星近似为：

$$K \approx 28.4 \text{ m/s} \left(\frac{P}{1 \text{ yr}}\right)^{-1/3} \left(\frac{M_p \sin i}{M_J}\right) \left(\frac{M_*}{M_\odot}\right)^{-2/3}$$

关键之处：方法测量的是 $M_p \sin i$，其中 $i$ 是轨道倾角（相对视线方向），而非真实质量 $M_p$。只有当凌星信号确认 $i \approx 90^\circ$ 时，才能确定真实质量。

仪器技术：精度革命#

第一代：碘气池校正#

1990年代，马西（Geoff Marcy）和巴特勒（Paul Butler）等人将碘气泡（iodine cell）技术推向实用，在目标恒星光穿过碘气泡时叠加已知的碘分子吸收谱，以其精确已知的谱线位置作为内部波长参考。这将精度从约数百米每秒提升至约2–3 m/s，使探测热木星成为可能。

第二代：光纤稳定高分辨率光谱仪#

HARPS（High Accuracy Radial velocity Planet Searcher，安装于ESO 3.6米望远镜，2003年投入使用）利用真空密封的温控光谱仪腔体和单模光纤耦合，将系统精度推至约1 m/s，能够探测到数个地球质量的行星（超级地球，如Gl 876 d）。

第三代：亚米每秒精度#

ESPRESSO（ESO VLT，2018年）目标精度约0.1 m/s，NEID（WIYN，2021年）约0.3 m/s，目标是探测真正的类地行星（需要约0.09 m/s精度才能探测日-地系统的类似物）。极大的技术挑战来自恒星本身的"噪声"——恒星活动（stellar jitter）：黑子、耀斑、对流颗粒引起的表面径向速度抖动约0.5–5 m/s，往往掩盖小质量行星信号。如何将仪器噪声和恒星噪声分离，是当前研究的主战场。

视向速度法的独特价值#

两种方法相辅相成：对于有凌星信号的行星，RV提供质量，凌星提供半径，两者结合可得行星密度，进而推断内部成分。

重要历史发现#

• 飞马座51 b（51 Peg b）（1995年，Mayor & Queloz）：第一颗围绕类日恒星的系外行星，RV半幅 $K \approx 55$ m/s，轨道周期4.23天，行星质量 $\geq 0.46 M_J$。这一发现改变了整个行星科学领域，两人荣获2019年诺贝尔物理学奖。
• Gl 876 系统（1998年）：第一颗围绕M型红矮星的行星，含有轨道共振（1:2）的双行星系统，是多行星系统轨道动力学研究的经典样本。
• HD 40307 系统（2009年）：三颗超级地球（最小约4 $M_\oplus$）的系统，HARPS RV数据，示范了对小质量行星的探测能力。

未来展望#

极大型望远镜（ELT，口径39米，预计2028年）配合下一代稳定光谱仪（ANDES）计划将RV精度推至亚厘米每秒量级，可能直接探测Alpha Centauri系统中的类地行星（若存在）。但恒星活动带来的基本噪声下限估计约为数十厘米每秒，要突破这一极限需要太阳物理层面对恒星表面速度场的深入理解。

视向速度法的哲学：倾听恒星的"心跳"#

视向速度法有一种深刻的间接性：我们从未直接看到被探测的行星，我们看到的是行星对恒星施加的引力影响——恒星在引力拉拽下发生的微小摇摆，通过光谱中极微小的多普勒偏移呈现给我们。

这种间接性并不意味着弱势。事实上，我们通过恒星光谱"感受"到行星的存在，在某种意义上类似于通过一个人的动作感受到另一个人的存在，而非直接目击。对于天文学家而言，这种间接探测提供了行星质量的最直接测量——这正是凌星法不能独立完成的事情。

RV法在1995–2010年间主导了系外行星探测，催生了一大批独立发现：在开普勒时代之前，约500颗已知系外行星中约90%来自视向速度法。这些早期发现的行星大多是热木星（因其信号强大，易于探测），但随着HARPS精度的提升，越来越多超级地球（如55 Cancri e、Gl 581系统）被揭示，预示着一个全新的多样性时代。

近邻星系外行星系统#

利用RV法，天文学家对太阳周围10秒差距内的大量恒星进行了系统性测量：

Alpha Centauri系统#

距太阳最近的恒星系统（约4.37光年）由三颗恒星组成：Alpha Cen A（G型，类太阳）、Alpha Cen B（K型）和Proxima Cen（M5.5型，距地球4.24光年）。

• Proxima Cen b（2016年，Anglada-Escudé et al.，Nature）：RV探测到最小质量约1.3 $M_\oplus$、轨道周期11.2天（位于宜居带）的信号，引发巨大关注。然而，Proxima Cen是极活跃的耀星，其RV jitter约1–2 m/s，与行星信号同量级，争议持续。2020年ESPRESSO对更精确数据的分析再次确认该信号（Suárez Mascareño et al. 2020）。
• Proxima Cen d（2022年）：ESPRESSO探测到最小质量约0.26 $M_\oplus$（火星到地球之间）、周期5.1天的超小质量行星候选。若确认，将是迄今RV探测的最小质量系外行星之一。

巴纳德星#

距地球约6光年（第二近的单星），2018年Ribas等人报告一颗候选行星信号（周期233天，最小质量约3.2 $M_\oplus$），但2021年后续分析认为该信号可能是恒星活动的假阳性。这个案例深刻说明了分离恒星活动与行星信号的挑战。

恒星活动：视向速度法的主要噪声来源#

恒星活动（stellar activity）是限制RV精度最难克服的障碍。以下机制产生虚假或叠加的径向速度信号：

• 黑子和耀斑：低温黑子减少恒星表面的辐射，使光谱谱线形状不对称，产生视向速度"抖动"，幅度约1–10 m/s（类太阳恒星）到数十 m/s（活跃M矮星）
• 对流颗粒（granulation）：恒星表面对流单元的上升和下降气流在分钟到小时时标上产生约1 m/s的随机速度扰动
• 磁活动周期：类似太阳11年活动周期，使恒星平均RV以数年为周期变化数十厘米每秒

应对方法：

1. 多指标同步监测：Ca II H&K谱线、H-alpha线的变化与磁活动相关，可与RV信号一起拟合，分离恒星活动贡献
2. 高采样率均化：短时白噪声（granulation）可通过在数小时内多次观测取平均来压低
3. 机器学习/GP回归：高斯过程（Gaussian Process）回归被用于同时拟合行星信号和恒星活动的协方差结构（Haywood et al. 2014）

自相关与贝叶斯分析：现代RV数据处理#

现代RV数据分析远超简单的Lomb-Scargle周期图。高精度数据中多颗行星的信号混叠、恒星活动噪声和仪器系统误差共同存在，需要更复杂的统计方法：

贝叶斯马尔可夫链蒙特卡洛（MCMC）：对行星轨道参数（周期、半幅K、偏心率等）的后验概率分布进行全采样，而非单纯找最优拟合，给出参数的置信区间和参数之间的相关性。对于多行星系统，参数空间维度高达20维以上，需要嵌套采样（nested sampling）等专门算法。

高斯过程（GP）建模恒星活动：将恒星活动的RV贡献建模为具有特定协方差核的高斯过程（如准周期性核，与恒星自转周期相关），同时与行星轨道信号共同拟合。这使原本被恒星噪声掩盖的小质量行星信号能够被提取出来（Haywood et al. 2014；Rajpaul et al. 2015）。

谱线形态分析：用谱线的等值宽度（bisector span）和FWHM随时间变化监测恒星活动，活动引起的RV变化通常伴随bisector span的相关变化，而行星引起的仅有RV变化——这是区分行星和星斑的重要诊断工具。

凌星时变（TTV）：对RV的补充与竞争#

对于开普勒发现的多行星系统，行星间的引力摄动会导致凌星时刻偏离等间隔，形成凌星时变（Transit Timing Variations, TTV）信号。TTV：

• 可探测质量低至几分之一地球质量的行星（对近轨道共振系统）
• 直接给出行星质量而无 $\sin i$ 简并
• 但需要长时间连续凌星观测，且信号解读复杂

对于不凌星的行星，RV是唯一的质量测量手段；对于同时有凌星的系统，RV+凌星组合仍是最可靠的方法，TTV则提供了独立交叉验证。

跨域连接#

• 天体测量与盖亚卫星：本方法只测视线方向那一个速度分量，因此拿到的永远是 $M_p \sin i$，一个质量下限。天体测量测的是天球面上另外两个分量，两者拼起来才定出真实倾角与真实质量。推论是：对同一系统，天体测量若测不到摆动，等于给真实质量加了个上限，与视向速度的下限一起把质量夹住。
• 多普勒效应：整套灵敏度都压在 $\Delta\lambda/\lambda = v_r/c$ 这一个比值上，正文给出十米每秒对应可见光波段极微小的波长偏移。由此可见方法的物理边界：分母是光速，想提高精度只能加长时间基线、增加参与拟合的谱线数目，而不能改变这个比例关系本身。
• 随机过程：黑子与对流颗粒造成的抖动不是白噪声，而是带协方差结构的随机过程，其时标与恒星自转周期绑定。把它建成准周期核的高斯过程，就能与严格周期的开普勒信号分开。但当行星周期接近自转周期或其谐波时，二者原则上不可分辨——正文里巴纳德星的争议正是这一情形。
• 蒙特卡洛方法：多行星系统的参数空间可达二十维以上，真正要的不是最优拟合点而是后验分布：周期别名、偏心率与半幅之间高度相关，最优点常常并不典型。嵌套采样还给出模型证据，让"到底几颗行星"成为可比较的假设，而不是靠肉眼看周期图定夺。
• 光谱学：第一代精度革命靠碘气池——让星光先穿过碘分子，用位置由分子结构固定的吸收谱当作叠在同一张底片上的波长尺。要点在于参考尺与目标光走同一条光路，因此仪器漂移对两者同等作用，相减即抵消；这也解释了为什么后来改用真空恒温腔体。

参考文献#

• Mayor, M., & Queloz, D. (1995). "A Jupiter-mass companion to a solar-type star." Nature, 378, 355–359.
• Pepe, F. et al. (2021). "ESPRESSO at the VLT." Astronomy & Astrophysics, 645, A96.
• Fischer, D. A. et al. (2016). "State of the Field: Extreme Precision Radial Velocities." Publications of the Astronomical Society of the Pacific, 128, 066001.
• Mayor, M. et al. (2009). "The HARPS search for southern extra-solar planets." Astronomy & Astrophysics, 507, 487–494.（超级地球批量发现的代表性论文）

数据来源：HARPS/ESPRESSO观测计划（ESO Archive）；NASA系外行星档案；上述引用文献。