跳到主要内容
Sandbox Physics
宇宙距离阶梯怎样获得尺度

L1 · 宇宙距离阶梯

从反平方亮度到恒星距离标定

知道恒星本征光度,就能把实测流量换成距离。先从一个球面开始:同样的发光功率传到远处,会分散在更大面积上。星等把这种反平方关系写成加法,距离便藏在视星等与绝对星等的差里。

同样的功率铺到球面上

设光度 L 为指定波段里每秒发出的能量,单位瓦;流量 F 为同一波段每平方米每秒收到的能量,单位瓦每平方米。先假定光源向各方向均匀发光,空间透明且满足欧几里得几何,距离 d 用米。每个包围光源的球面收到相同总功率,除以球面面积就得到流量。

距离加倍,球面积变成四倍,每平方米收到的流量降到四分之一。近处还可以从角视差得到独立距离:一天文单位的基线,在一秒差距处张开一角秒。对带噪实测视差求距离时还要处理误差;这里先写理想几何关系。

F=L4πd2,d [pc]=1ϖ [arcsec]F=\frac{L}{4\pi d^2},\qquad d\,[\mathrm{pc}]=\frac{1}{\varpi\,[\mathrm{arcsec}]}
光度和流量必须指同一波段。视差使用角半振幅与一天文单位基线;整个轨道的两天文单位直径,对应角位移的峰峰值。

对应阅读:Riess et al. · Reid

相减星等,得到距离模数

星等越小,光越亮。对零点相同的两次流量测量,星等差等于流量比的常用对数乘负二点五。绝对星等 M 指同一光源放在十秒差距处、消光已去掉时的视星等。这里的视星等 m 也采用相同波段约定,并作相应改正。

把距离 d 处的流量除以十秒差距处的流量,再代入星等差。距离平方带来的系数二与二点五相乘,变成五。得到的距离模数,是同一星系全部恒星共同增加的一项。对宇宙学远源,要用光度距离代替欧氏距离,波段 K 改正也应纳入亮度模型。

m1−m2=−2.5log⁡10 ⁣F1F2,μ≡m−M=5log⁡10 ⁣d10 pcm_1-m_2=-2.5\log_{10}\!\frac{F_1}{F_2},\qquad \mu\equiv m-M=5\log_{10}\!\frac{d}{10\,\mathrm{pc}}
流量比没有量纲。距离增大时,视星等和距离模数都增大。
d [Mpc]=10(μ−25)/5d\,[\mathrm{Mpc}]=10^{(\mu-25)/5}
一百万秒差距等于一百万个秒差距,和十秒差距参考长度相比,换成这一单位后常数为二十五。

对应阅读:Riess et al.

用脉动周期选择平均光度

造父星的周期能预示平均光度,不过还要处理恒星族与观测改正。周期 P 用天,取对数后以十天为中心。这样,零点就是十天造父星的平均绝对星等。斜率为负,表示周期越长,恒星越亮。

Wesenheit 星等把几个实测波段按固定颜色系数组合起来,以减小消光影响。它是组合出的亮度量,不是一片新滤镜。剩余的消光律差异、恒星拥挤和金属丰度仍可能影响结果。设 Z 为相对参考值的金属丰度对数,单位 dex;系数 gamma 描述它对应的亮度变化。

mW=μ+M10W+b[log⁡10 ⁣(Pday)−1]+γZm^W=\mu+M^W_{10}+b\left[\log_{10}\!\left(\frac{P}{\mathrm{day}}\right)-1\right]+\gamma Z
视 Wesenheit 星等、距离模数和十天零点都用星等。斜率 b 的单位是每个周期 dex 对应的星等,gamma 为每个金属丰度 dex 对应的星等。随机散布与协方差在后面加入。

对应阅读:Riess et al. · Yuan et al.

用一行真实数据完整算一次

取 M101 编号54672的恒星,也就是原矩阵第0行,周期为6.853491天。发布矩阵先减掉了每 dex 负3.285星等的基准周期项,存下来的23.15228653是预处理量,不是原始 H 波段星等。加回基准周期项后,Wesenheit 值为23.69131604。它与另列的 F160W 星等24.04是不同量。

为了练习,暂取十天零点负5.80、斜率负3.30,并在这道算例里忽略金属丰度。周期偏移为负0.16408813 dex,算出的绝对星等为负5.25850917。用23.69131604减去它,得到模数28.94982521,距离约6.165 Mpc。这里先给定标定值以展示计算;后面的八参数拟合会同时估计它们。

量保存或计算数值
周期6.853491 day6.853491\,\mathrm{day}
预处理观测23.1522865323.15228653
恢复的 Wesenheit 值23.6913160423.69131604
假定标定所得模数与距离28.94982521,6.165 Mpc28.94982521,\quad6.165\,\mathrm{Mpc}

对应阅读:PantheonPlusSH0ES DataRelease · Riess et al.

几何距离怎样定下零点

NGC4258 的水脉泽有可测角位置、多普勒速度和加速度,可以放入盘的几何距离模型。2019年结果为7.576 Mpc,统计、系统误差分别为0.082和0.076 Mpc。Riess 等采用参考模数29.398,不确定度约0.032星等;选定矩阵保留发布时的精确锚约束。

距离已知后,把造父星关系中的距离项移到另一边,就能拟合光度零点。红巨星支尖端的标定也是同一道理,只是改用年老恒星族亮度分布的尖端。两者恒星物理和选样不同,几何起点却可能相同。比较时要保留这项共同依赖。

对应阅读:Reid · Riess et al.

年老恒星族提供另一把标定尺

红巨星支尖端利用年老红巨星亮度分布的一个清楚上边缘。标定取决于波段、颜色和恒星族。Reid 等讨论的 F814W 约定里,NGC4258 锚给出尖端绝对星等约负4.01±0.044.01\pm0.04。只要把恒星分辨得足够清楚,即使年轻造父星稀少,也可能测到这个边缘。

做一道计算例题:固定上述颜色、恒星族标定,采用参考模数29.398。把尖端绝对星等与模数相加,得到预言的消光修正视星等25.388。这个数由给定标定算出,两项采用同一测光约定。换一种恒星指标后,距离锚仍可以是同一个。

mtip,0=Mtip+μ=−4.01+29.398=25.388m_{\mathrm{tip},0}=M_{\mathrm{tip}}+\mu=-4.01+29.398=25.388
下标零表示已作消光改正,数值例题固定使用来源中的标定约定。

对应阅读:Reid · Riess et al.

亲手观察

检查距离标定的方向

  1. 保持“距离阶梯”模式,把采用的锚模数增加0.10。先预测各个推断距离往哪边变,再看数值。
  2. 解释为什么采用的距离改变后,实测造父星周期与亮度仍保持不变。

练习与答案

恒星放到两倍距离,流量比与星等变化各是多少?

流量降到四分之一。星等增加二点五乘以四的常用对数,即1.50515星等。正的星等变化对应变暗,符号相符。

标定让采用的绝对星等亮0.10星等,视星等不变。推断距离如何变化?

变亮表示绝对星等减小0.10。距离模数增加0.10,距离乘以十的零点零二次方,约1.04713。

有人把发布矩阵第0行当成 F160W 星等,又拟合完整的负3.3周期斜率。错在哪里?

预处理值已经采用 Wesenheit 颜色组合,并减掉基准周期斜率。应拟合斜率修正量,或先恢复原周期项后再拟合完整斜率。还要区分 F160W 单波段亮度和 Wesenheit 组合亮度。

来源与计算范围

  1. Riess et al., A Comprehensive Measurement of the Local Value of the Hubble Constant, arXiv:2112.04510 v 3 (2022)

    第2.1节与第4节定义联合模型和标定;第5节、表5、图14—16给出完整拟合与相关性,第6—7节讨论变体和比较。课程所述论文结论固定在这一版本。

  2. PantheonPlusSH0ES DataRelease, commit c447f0fea703fcd0fff57de5000947b5ca81286b

    ally、alll、allc FITS 数组包含3492项预处理观测、47个参数列及完整协方差。公开教程保留其中87项观测与涉及的全部非零模型列。table2.README 明确提醒,不能只拿表2的对角误差重新拟合。

  3. Reid, Pesce & Riess, An Improved Distance to NGC4258 and its Implications for the Hubble Constant, arXiv:1908.05625 v 2 (2019)

    水脉泽几何距离为7.576 Mpc,统计误差0.082 Mpc、系统误差0.076 Mpc。论文也用此锚标定红巨星支尖端,并讨论起点、终点相同时的比较。

  4. Yuan et al., A New Cepheid Distance to the Maser-Host Galaxy NGC4258, arXiv:2203.06681 v 1 (2022)

    第2—4节介绍 HST 观测、造父星选择与周期—光度拟合。发现星表与近红外标定样本大小不同;发布矩阵含443颗 NGC4258 造父星。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策