跳转到内容
← 返回知识库
恒星物理约 21 分钟阅读

造父变星与变星

造父变星变星周光关系宇宙距离哈勃常数不稳定带

概述

造父变星(Cepheid variable stars,或 Cepheids)是一类以周期性脉动为特征的超巨星,其光度随时间规律性起伏,周期从约 $1$ 天到约 $100$ 天不等。它们的重要性远超一颗闪烁的星星:1908年,亨利埃塔·勒维特(Henrietta Swan Leavitt)发现造父变星存在严格的周光关系——亮度越高、脉动周期越长——这一关系成为宇宙距离测量的基础,直接导致了河外星系的发现,推动了现代宇宙学的诞生。时至今日,造父变星仍是校准宇宙距离阶梯、测量哈勃常数的关键工具,置身于"哈勃张力"这一当代宇宙学核心争议的中心。

从观测宇宙学角度看,造父变星的价值不在于把单颗星的内部演化故事讲完,而在于它把可测的周期翻译成绝对光度,从而把视星等差变成距离模数。距离阶梯的每一级都依赖这种可转移的标准烛光:近处用几何方法钉死零点,远处用同一类天体把尺子接长。因此本文重心放在阶梯衔接、零点校准与 H0H_0 系统误差,而不是重述大质量恒星演化全程。

变星的物理机制:κ 机制

恒星为什么会脉动?恒星内部的热核反应是稳定的,不会自然产生脉动——必须存在某种驱动机制将能量转化为相干振荡。对于造父变星,主导机制是κ机制(kappa mechanism,κ\kappa 是不透明度的符号):

  1. 当恒星包层被压缩时,氦的部分电离区(位于约 10,00010{,}000–40,00040{,}000 K 的温度范围)变得更不透明(κ\kappa 增大),吸收辐射而非让其通过
  2. 积聚的辐射热量使包层受热膨胀,如同受压气体推开活塞
  3. 膨胀中包层降温,不透明度降低,辐射得以逃逸,驱动减弱
  4. 重力将包层拉回,压缩再次开始,循环重复

这一过程的关键在于:部分电离区的不透明度对温度和密度的响应方式,使压缩时不透明度增大(而非减小)——这与正常气体相反,因此能积聚能量驱动振荡,而非耗散能量阻尼振荡。天文学家将赫罗图上造父变星等脉动变星存在的区域称为不稳定带(instability strip),位于约 Teff≈5,000T_{eff} \approx 5{,}000–8,0008{,}000 K 的垂直带状区域,贯穿从主序(δ\delta 盾型变星)到白矮星前身(ZZ Ceti 型变星)的多种脉动变星。

对测距而言,不稳定带给出的是样本选择边界:只有落在带内、光变幅度与周期可可靠测定的星,才适合进入周光关系拟合。带外天体即使亮,也不能随便当作造父变星烛光使用。

驱动层的深度还决定了光变曲线的形状:基模造父变星常呈现快速上升、缓慢下降的锯齿状轮廓,傅立叶分解得到的相位差与振幅比可作为模式判别量。对测距流水线来说,这些参数不是装饰——它们用来把混入样本的泛音星、混合模式星或非造父变星剔除,避免把错误的周光零点带进宿主星系距离。不稳定带的蓝缘与红缘还会随金属丰度微移:金属更富时带略红,金属更贫时带略蓝,因此跨星系比较时必须把颜色—周期选择与金属丰度改正绑在同一套程序里。

勒维特定律:周光关系

亨利埃塔·勒维特(Henrietta Swan Leavitt,1868–1921)是哈佛天文台的"人工计算员"(computers)——一个主要由女性组成的团队,负责系统分析天文底片。1908年,勒维特在小麦哲伦云中识别出 1,7771{,}777 颗变星(其中约 47 颗被她归为造父变星),并注意到较亮的变星周期更长:星等越亮,周期越长。

由于小麦哲伦云中所有恒星到地球的距离近似相等(约 62 kpc),视亮度差异直接反映绝对光度差异。1912 年,勒维特挑选其中周期数据可靠的 $25$ 颗造父变星,将这一关系量化为线性规律(Leavitt & Pickering 1912)。短讯与相对台阶见 勒维特 1912 周光关系。现代校准的造父变星周光关系(Leavitt Law)形式为(以 $V$ 波段为例):

MV=−2.81log⁡10(P/day)−1.43±0.05(经金属丰度修正)M_V = -2.81 \log_{10}(P/\text{day}) - 1.43 \pm 0.05 \quad \text{(经金属丰度修正)}

(Freedman et al. 2001,The Astrophysical Journal,553,47;基于哈勃空间望远镜大麦哲伦云造父变星校准)近红外($H$ 波段和 WHW_H 带)的周光关系斜率更陡、离散更小,尘埃消光影响更弱,是现代测量首选(Monson & Pierce 2011,The Astrophysical Journal Supplement Series,193,12)。

观测程序通常分三步:先测周期与相位平均星等;再做消光与金属丰度修正;最后把校正后的绝对星等接到宿主星系距离。任何一步的系统偏差都会乘性地抬升或压低整条阶梯上方的全部距离,因此周光关系不是"多一个经验拟合",而是局部宇宙膨胀率的标尺原点之一。

近红外与中红外周光关系之所以更受青睐,是因为尘埃消光随波长下降,相位平均星等对周期的离散也更小。光学波段仍不可少:它提供颜色信息,用来估计 reddening 与温度,并把近红外绝对星等接到历史上以光学为主的样本。实践中常见做法是构造 Wesenheit 指数——用两个波段的线性组合近似消光无关量——再拟合周光斜率与零点。斜率若随金属丰度或宿主环境漂移,零点校准就不能只靠一个锚点星系;这正是盖亚银河系样本与 LMC 几何距离被同时使用的原因。

造父变星的分类

造父变星并非单一类型,而是包含几个不同的子群:

经典造父变星(Classical Cepheids,I 型)

经典造父变星(Cepheid Ia 或 I 型)是大质量、年轻的超巨星(质量 $4$–20 M⊙20\,M_\odot,年龄 10710^7–10810^8 yr),位于薄盘和旋臂中。它们的脉动模式主要是径向基模(fundamental mode)。

造父座 δ\delta(δ\delta Cephei)是原型星,周期 $5.366$ 天,亮度在 $V$ 波段从 $3.48$ 到 $4.37$ 等之间变化,距离约 $273$ pc(Gaia 视差测量,Benedict et al. 2007)。

W Vir 型(II 型造父变星)

II 型造父变星(Cepheid II 或 Population II Cepheids)是低金属丰度的老年恒星(质量约 $0.5$–0.6 M⊙0.6\,M_\odot,年龄 >1010> 10^{10} yr),主要出现在球状星团和晕中。关键区别:在相同周期下,II 型造父变星比经典(I 型)造父变星暗约 $1.5$–$2$ 等。

历史上,混淆 I 型和 II 型造父变星导致了严重的距离标度错误。1952 年,沃尔特·巴德(Walter Baade)发现仙女座大星云(M31)中的造父变星是 I 型而非 II 型,由此将 M31 的距离修正为原来的约两倍,宇宙年龄估计也相应修正(Baade 1952)。

这一事件成为距离阶梯方法论的经典教训:同一名称不等于同一烛光。阶梯衔接处必须先证明两端样本属于同一观测类别(周期范围、星族、消光处理一致),再谈提高统计精度;否则再多的数据也只是把错的尺子量得更细。

RR Lyrae 型

RR Lyrae 型变星是脉动的老年水平支恒星(质量约 $0.6$–0.8 M⊙0.8\,M_\odot),周期约 $0.2$–$1$ 天。它们同样遵循近似的光度-周期关系(光度几乎一致,MV≈0.5M_V \approx 0.5 mag,因此也可用作标准烛光),但只能用于银河系内和近邻星系(比造父变星更暗约 $2$–$3$ 等)。

球状星团因含大量 RR Lyrae 型变星而成为校准这类标准烛光的天然实验室(Cacciari & Clementini 2003,Stellar Candles for the Extragalactic Distance Scale)。在距离阶梯中,RR Lyrae 多用于银河系内与近邻矮星系的局部锚点,与造父变星的河外宿主路径互补,而不是互相替代。

$\delta$ 盾型(Delta Scuti)变星

δ\delta 盾型变星是主序上 A–F 型恒星的径向和非径向脉动,周期约 $0.01$–$0.3$ 天,振幅较小(约 $0.003$–$0.9$ mag)。开普勒卫星发现了数千颗,正在用于星震学精确测定恒星内部结构(Murphy et al. 2019,Monthly Notices of the Royal Astronomical Society,485,2380)。

长周期变星(Mira 型)

米拉型变星(Mira variables)是 AGB 巨星的大振幅脉动(振幅 $> 2.5$ mag,周期 $100$–1,0001{,}000 天),是晚期恒星演化的重要阶段,与大量尘粒驱动的质量损失密切相关。鲸鱼座 $o$(Mira,刍藁增二)是原型星,周期约 $332$ 天,距离约 $92$ pc。

哈勃常数测量与哈勃张力

造父变星是当代宇宙学最核心争议的主角之一——哈勃张力(Hubble tension)。

距离阶梯路径

宇宙距离阶梯(cosmic distance ladder)从造父变星出发:

视差⏟盖亚→造父变星⏟大麦哲伦云/银河系→遥远星系造父变星⏟HST/JWST→Ia型超新星⏟宇宙学→H0\underbrace{\text{视差}}_{\text{盖亚}} \rightarrow \underbrace{\text{造父变星}}_{\text{大麦哲伦云/银河系}} \rightarrow \underbrace{\text{遥远星系造父变星}}_{\text{HST/JWST}} \rightarrow \underbrace{\text{Ia型超新星}}_{\text{宇宙学}} \rightarrow H_0

关键在于"传递":几何锚点(盖亚视差、LMC 食双星)把周光关系零点钉死;造父变星再把零点带到能同时看见 Ia 型超新星的宿主星系;超新星样本才把距离延伸到 Hubble flow。任何一级的零点漂移,都会在最终 H0H_0 上表现为整体偏移,而不是局部噪声。

2022 年,Riess 等人(SH0ES 团队)利用 HST 校准的 42 个宿主星系造父变星,结合 Ia 型超新星距离,得到 H0=73.04±1.04H_0 = 73.04 \pm 1.04 km/s/Mpc(Riess et al. 2022,The Astrophysical Journal Letters,934,L7)。

CMB 路径的张力

独立地,Planck 卫星通过宇宙微波背景的声学峰测量,在 Λ\LambdaCDM 框架下给出 H0=67.4±0.5H_0 = 67.4 \pm 0.5 km/s/Mpc(Planck 2018)。两者之间约 5σ5\sigma 的差异——"哈勃张力"——是当代宇宙学最严峻的挑战之一。可能的原因包括:

  • 造父变星测量的系统误差(尘埃消光、金属丰度效应、周光关系的非线性)
  • 声学峰分析中的系统误差
  • 超越 Λ\LambdaCDM 的新物理(早期暗能量、超流暗物质等)

距离阶梯路径与 CMB 路径几乎不共享系统误差源:前者依赖晚宇宙标准烛光与红移—距离拟合,后者依赖早期声学尺度在 Λ\LambdaCDM 下的外推。正因为正交,张力才尖锐——它强迫观测者把造父变星零点、超新星量纲与早期物理假设拆开检验。

2024 年,詹姆斯·韦伯太空望远镜(JWST)开始独立校准近红外造父变星(Freedman et al. 2024,The Astrophysical Journal,961,L16),初步结果给出 H0≈69H_0 \approx 69–$72$ km/s/Mpc,暗示 HST 造父变星可能存在一定系统误差,但哈勃张力尚未完全解决。

观测宇宙学争论因此高度依赖"同一天体、不同望远镜、不同波段"的交叉检验:若差异随宿主星系环境、拥挤场改正或近红外消光处理系统变化,更像阶梯系统误差;若在多种独立烛光路径上仍稳定存在,才更值得认真考虑早期宇宙新物理。造父变星位于这场检验的枢纽,而不是旁证。

交叉检验还包括换一类标准烛光:红巨星支先端(TRGB)、表面亮度涨落、水脉泽几何距离等路径不经过造父变星周光关系,却同样接到 Ia 型超新星或 Hubble flow。若造父路径与 TRGB 路径在同一批宿主上系统性分歧,优先怀疑的是拥挤改正、消光律或样本选择;若多条晚宇宙路径彼此一致却仍与 CMB 外推冲突,争论才更偏向早期物理或宇宙学模型扩展。JWST 的贡献不只是"看得更远",而是在近红外把造父变星从拥挤的盘面星场中分得更开,从而单独报价拥挤偏差这一项。

为什么造父变星如此重要

造父变星的故事是科学史上最精彩的章节之一。一个约 10,00010{,}000 K 的恒星表层不断收缩和膨胀的物理过程,竟然成为测量宇宙大小的标尺。哈勃(Edwin Hubble)于 1924 年用威尔逊山望远镜在仙女座"星云"中找到造父变星,计算出其距离约 $100$ 万光年(远超银河系大小),一举证明宇宙中存在其他星系——这是 20 世纪最重大的天文发现之一。

观测宇宙学此后长期把造父变星放在阶梯中腰:上接几何视差,下接超新星宇宙学。角色稳定,争议却转向零点与系统误差——这说明方法已成熟到足以支撑高精度 H0H_0 声明,也成熟到任何百分之一级偏差都会被放大成宇宙学危机。

勒维特的发现来自对大量底片的系统分析,她却未能活到哈勃用她的工具改变宇宙观的那一天(1921 年去世)。1925 年,诺贝尔奖委员会曾联系哈佛天文台,希望提名勒维特,却被告知她已逝世。诺贝尔奖从不追授,但历史记住了她。

就方法史而言,周光关系把"变星监测"从分类学提升为宇宙学测距:周期是时间量,几乎不受距离本身影响;一旦零点被几何方法固定,远处样本只需重复同一套测光与消光程序。正因如此,现代争论才集中在零点、金属丰度与拥挤改正,而不是质疑"周期可否观测"。

跨域连接

  • 振动与简正模式:径向脉动的周期本质是声波穿越恒星的驻波周期,因而随平均密度下降而变长;再叠加"光度越高半径越大",就得到周期与光度的单调关系。周光关系因此是驻波条件的推论,不是经验巧合——基模与泛音是同一根弦的不同谐波,斜率必然不同。
  • 卡诺循环与热机:自激振荡要求热量在压缩相被吸收、在膨胀相被释放,κ 机制正是用不透明度的反常温度依赖凑出这个相位。推论:不稳定带之所以窄,是因为氦部分电离区必须落在包层的特定深度——太热则太浅,太冷则被对流短路。
  • 测量等价性:跨群体比较要求工具在各群体中含义相同。同周期下两类造父变星差约一点五到两等,把它们当同一把尺子用,距离就整体错一倍——历史上确实这样发生过。每一级阶梯都必须先证明"两端是同一类天体",这比提高精度更优先。
  • 女性主义:周光关系出自哈佛天文台由女性组成的"计算员"团队对上千张底片的系统比对。把发现建立在批量重复劳动上,与这些劳动者被排除在观测和理论岗位之外,是同一件事的两面;勒维特的署名与提名境遇正是这种错位的后果。
  • 哈勃张力:周光关系的零点、金属丰度修正与尘埃消光都是乘性系统误差,会整体平移阶梯上方的全部距离,因而直接平移哈勃常数;微波背景那条路径完全不经过它们。若差异源于造父变星,就应随宿主环境系统变化——这可检验,但目前尚无定论。

银河系内造父变星的精确校准

随着盖亚(Gaia)卫星提供高精度视差数据,银河系内约 2,5002{,}500 颗经典造父变星的距离已被精确测量到 $1$–5%5\% 的精度(Ripepi et al. 2023,Astronomy & Astrophysics,674,A17),极大地改善了周光关系的零点校准:

  • 视差几何校准:对于距离 $< 5$ kpc 的造父变星,盖亚视差可以直接定标周光关系,消除了对大麦哲伦云(LMC)距离中间步骤的依赖
  • 大麦哲伦云几何距离:DEBs(分离双食双星)在 LMC 中提供了几何距离 49.59±0.09±0.5449.59 \pm 0.09 \pm 0.54 kpc(Pietrzyński et al. 2019,Nature,567,200),结合 LMC 中的造父变星,进一步锚定周光关系的零点
  • 金属丰度修正:不同金属丰度的造父变星在周光关系上有系统偏差,约 $-0.2$–$-0.4$ mag/dex,在比较不同环境的造父变星时必须修正

这些改进使哈勃常数的系统误差在 $2030$ 年代有望降至 1%1\% 以下,届时哈勃张力是否仍然存在将有明确答案。

对距离阶梯而言,盖亚与 LMC 几何锚点的意义是缩短依赖链:更少中间假设,意味着更少乘性系统误差的藏身之处。与此同时,宿主星系样本的选择效应、造父变星视场拥挤改正、以及不同金属丰度环境之间的零点迁移,仍是观测侧必须单独报价的不确定度项——它们属于测距程序,而不是恒星演化故事的附录。

造父变星的脉动模式

除基模(fundamental mode)脉动外,造父变星还可以以泛音模(overtone modes)脉动:

  • 第一泛音造父变星(First Overtone Cepheids):脉动振幅较小、周期较短,在赫罗图上的不稳定带内位于较蓝的区域。周光关系斜率与基模不同,需分开处理
  • 拍造父变星(Beat Cepheids):同时以两种模式脉动,光变曲线显示复杂的拍频图案,可提供星震学约束

对距离阶梯实践而言,模式识别直接影响该用哪条周光关系:把泛音造父变星误当作基模,周期—光度零点会错位,距离随之整体偏移。因此测光分类与傅立叶参数、颜色—周期图上的位置,都是观测宇宙学流程的前置质控,而不是可有可无的恒星物理细节。

拍造父变星同时给出两个周期,周期比接近恒星内部声速结构所允许的特征值,因而能约束平均密度与包层结构,反过来检验周光关系背后的物理假设是否自洽。对测距样本,拍星数量虽少,却是"斜率是否该随质量或金属丰度微调"的活体实验室:若理论预测的周期比与观测不符,就不应把单一线性周光关系无条件外推到所有宿主环境。现代测光巡天(如 OGLE、Gaia)把大量短周期样本按模式分开发布,正是为了让阶梯校准不再把不同物理状态的星混进同一条拟合。

参考文献

  • Leavitt, H.S. & Pickering, E.C. (1912). Periods of 25 Variable Stars in the Small Magellanic Cloud. Harvard College Observatory Circular, 173, 1–3.
  • Freedman, W.L. et al. (2001). Final Results from the Hubble Space Telescope Key Project. The Astrophysical Journal, 553, 47–72.
  • Riess, A.G. et al. (2022). A Comprehensive Measurement of the Local Value of the Hubble Constant. The Astrophysical Journal Letters, 934, L7.
  • Freedman, W.L. et al. (2024). Status Report on the Chicago-Carnegie Hubble Program (CCHP): Three Independent Estimates of the Hubble Constant Using JWST. The Astrophysical Journal, 961, L16.
  • Planck Collaboration (2020). Planck 2018 results. VI. Cosmological parameters. Astronomy & Astrophysics, 641, A6.

延伸阅读

  • Catelan, M. & Smith, H.A. (2015). Pulsating Stars. Wiley-VCH. — 脉动变星物理的全面参考书
  • Pietrzyński, G. et al. (2019). A distance to the Large Magellanic Cloud that is precise to one per cent. Nature, 567, 200–203. — LMC 几何距离的里程碑测量
  • Di Valentino, E. et al. (2021). In the realm of the Hubble tension — a review of solutions. Classical and Quantum Gravity, 38, 153001. — 哈勃张力及可能解决方案的综述