跳转到内容
← 返回知识库
宇宙学观测约 19 分钟阅读

甚长基线干涉测量

VLBI射电干涉角分辨率事件视界望远镜地球自转国际参考系

概述

甚长基线干涉测量(Very Long Baseline Interferometry,VLBI)是一种将分布在地球各地(乃至地球轨道上)的射电天线联合起来,形成等效口径等于天线间距的超大型干涉仪的技术。VLBI 的角分辨率可达微角秒(μ\muas)量级,是目前人类能够实现的最高天文分辨率——远超任何光学望远镜。它同时也是测量板块运动、地球自转和建立国际天球参考系的重要大地测量工具。

把地球直径当成口径,分辨率可以进到微角秒,但各站之间没有共用的本振电缆。这与甚大阵一类的连线干涉不同:后者能在本地把信号相干相加,VLBI 只能先把电压时间序列连同氢频标时间戳写进磁盘,再在相关机上做事后互相关。分辨率由最长基线决定,灵敏度由天线总面积和系统温度决定,两者不能互相代替。

因此 VLBI 的“地球尺度望远镜”是稀疏孔径。台站一少,空间频率覆盖就出现空洞,图像必须从可见度重建,而不是在焦平面上直接曝光。后文的 VLBA、EVN 与 事件视界望远镜,都是在这条约束下工作。

干涉测量原理:时延与互相关

两台相距基线长度 $B$ 的射电天线接收同一天体信号时,由于光程差的存在,信号到达两台天线有时间延迟 τ\tau:

τ=Bcos⁡θc\tau = \frac{B \cos\theta}{c}

其中 θ\theta 为天线基线方向与天体方向的夹角,$c$ 为光速。通过对两路信号进行互相关(记录到高精度时频标准——氢频标——后在相关机上离线计算),可以得到干涉条纹,干涉仪的角分辨率约为:

θres≈λB\theta_{\rm res} \approx \frac{\lambda}{B}

VLBI 基线 $B$ 最大可达地球直径(约 12700 km),在 1.3 mm 波长(事件视界望远镜工作波段)的衍射极限约为 20 μ20 \, \muas(微角秒),约为人眼分辨率的一千万分之一。

几何时延是科学量:它把基线矢量投影到视线方向上,源的位置稍有变化,τ\tau 就变。实际记录里还叠着钟差、大气与仪器时延,相关机必须在时延—时延率平面上搜相关峰,找到条纹以后再把这些干扰项一并解出。互相关给出的是复可见度,不是一张现成的像。

两路信号要对齐到远小于一个波长的等效光程,条纹才会积起来。氢频标把各站的时间轴钉在同一套物理秒上,互相关才有意义。基线越长、波长越短,对时钟和大气的要求越苛刻,这也是毫米波全球阵比厘米波阵更难的原因。

VLBI 的技术特点

独立时钟与氢频标

VLBI 各站不能实时连接传输数据(基线过长),信号必须记录在磁盘阵列上后运回相关中心。因此每台天线必须配备极高精度的独立时频标准——氢原子频标(hydrogen maser),日积累误差约 10−1410^{-14}–10−1510^{-15} 秒,确保各站记录的时间戳足够一致以进行事后相关。

氢原子频标锁定在氢基态超精细跃迁附近,为的是提供可独立运行的时间基准,而不是为了观测 21 厘米线。日稳定度到文中给出的量级,相干积分才不会因为相位漂移而把自己抹平。钟一散相,基线再长也换不回条纹——这与把天线面积做大是两回事。

各站还要记下气象和仪器相位校准音,供事后改正。VLBI 把“同时观测”定义成:在各自的钟面上,对同一天区写下可对齐的电压序列。

相关处理

收集各站磁盘后,在专用软件相关机(如 SFXC、DiFX)上进行互相关,搜索特定延迟和延迟率下的相关峰(干涉条纹)。现代 VLBI 数据率约 32 Gbps 至 256 Gbps,一次 EHT 全球实验可产生约 3500 TB 数据(EHTS,2019 年实验)。

现代生产相关大量使用软件相关机。DiFX 是分布式 FX 结构,在普通计算集群上把各站记录变可见度,已用于 VLBA、大地测量和毫米波实验(Deller 等,2011,PASP,123,275)。欧洲网一侧常用 JIVE 的 SFXC。相关中心要的不是实时图像,而是在延迟空间里把条纹找出来。

数据介质从磁带演进到磁盘模块。一部分欧洲实验可以用光纤把原始电压送回相关中心,缩短周转;全球毫米波阵的站点——尤其是南极——仍常常只能等季节、把硬盘运走。后文 EHT 的数据量,就是这条记录链的极端情形。

主要 VLBI 网络

超长基线阵列(VLBA)

美国超长基线阵列(VLBA)由 10 台 25 m 天线组成,分布从夏威夷到维尔京群岛,最大基线约 8611 km,主要运行在 1.4–43 GHz。VLBA 在射电星系喷流结构、活动星系核的厘米尺度成像和脉冲星天体测量中有重要贡献(Napier,1994,VLBA 科学文件)。

VLBA 是为 VLBI 专门建造的阵列,而不是把现成望远镜临时联网。十台 25 米天线从夏威夷莫纳克亚到维尔京群岛圣克罗伊,由新墨西哥索科罗的阵列运行中心遥控,1993 年前后进入常规使用(Napier 等,1994)。每站自备氢频标与记录系统,观测时并不把中频信号用电缆送到同一座相关机。

厘米到短厘米波段上,它适合把活动星系核的喷流核成分和脉冲星位置测到毫角秒以内。后文的天体测量与 MOJAVE 监测,大量依赖这套固定几何。它不是 单口径射电望远镜 的放大版:分辨率来自最长基线,不是来自某一面 25 米天线本身。

欧洲 VLBI 网络(EVN)

EVN 是由欧洲、亚洲和南非的天线组成的协同网络,最大基线约 10000 km,灵敏度比 VLBA 高(包含 100 m 埃弗尔斯伯格和 64 m 洛弗尔天线),适合高灵敏度的 FRB 精确定位和致密射电源成像。

EVN 的逻辑相反:它把欧洲、亚洲和南非已经存在的大天线组成协同网,相关由荷兰德温厄洛的联合研究所 JIVE 完成。100 米埃弗尔斯伯格与 64 米洛弗尔把收集面积抬上去,所以对暗弱源更敏感。最大基线仍是地球尺度,但 uv 覆盖随参与台站而变,不是每年固定的十站图案。

对 快速射电暴 重复源,高灵敏度加上相位参考,可以把爆发位置钉到毫角秒,从而对准宿主星系里的恒星形成区。这与 VLBA 的成图专长互补:一个更重灵敏度,一个更重全年可调度的固定基线。

事件视界望远镜(EHT)

EHT 是在 1.3 mm 和 0.87 mm 工作的全球 VLBI 网络,包含夏威夷、智利、亚利桑那、法国、南极等 8 个站点(2017 年观测时)。2019 年发布 M87 黑洞阴影图像(Event Horizon Telescope Collaboration,2019,ApJL 875,L1),2022 年发布人马座 A 图像,两项成果均为人类天文史上的里程碑。EHT 持续扩展台站(包括加入 NOEMA 和 ALMA 更多天线),目标是实现对黑洞周围相对论性等离子体的动态成像("黑洞电影")。

2017 年那次全球毫米波观测用了八个台站,包括智利的 ALMA 与 APEX、西班牙 IRAM 30 米、夏威夷的 JCMT 与亚毫米阵、墨西哥大型毫米波望远镜、亚利桑那亚毫米波望远镜,以及南极点望远镜。工作波长 1.3 毫米,对应大约 230 GHz。把地球直径代入衍射公式,就得到文首已经给出的微角秒量级极限;具体成像还受天气、相位校准和正则化先验约束。

2019 年发布的是 M87 中心的阴影图像,2022 年是人马座 A\*。两张图都是从可见度重建的,不是相机底片。台站极少,所以“看到的环”必须对算法做交叉检验。EHT 后来加入 NOEMA 等站,目标是把稀疏采样变得稍密一点,而不是改写这条几何极限。

VLBI 在大地测量与基础物理中的应用

国际天球参考系(ICRF)

VLBI 通过精确测量遥远类星体的相对位置,建立了国际天球参考系(ICRF)——地球和天文观测的基础坐标系。ICRF3(2018 年)包含 4536 个致密射电源,参考方向不确定度约 30 μ30 \, \muas(Charlot 等,2020,A&A 644,A159)。

选遥远致密射电源,是因为它们在天球上几乎不显示可测的周年视差,适合当方向基准。活动星系核的喷流会在毫角秒上缓慢移动,所以参考架要不断剔除结构变化大的源,ICRF3 才把不确定度压到文中给出的微角秒量级(Charlot 等,2020)。光学星表要和它对齐,才能把空间卫星的姿态接到同一套惯性方向上。

板块运动与地球自转

各站 VLBI 测量的射电源方向差异积累,可以精确测定各台站的相对位置变化,从而确定板块漂移速率(精度约 1 mm/yr)和地球自转不规则性(UT1-UTC 精度约 1 μ1 \, \mus)。VLBI 数据是国际地球自转与参考系服务(IERS)维护全球时间和参考系的核心输入。

同一套时延,既是天体的方向,也是台站的坐标。把源的位置当已知,解出来的就是基线矢量随时间的变化,年速率到毫米。国际地球自转与参考系服务把 VLBI 的 UT1 与台站坐标,和卫星测地、月球激光测距并在一起维护。天文与 大地测量 在这里用的是同一类相关峰。

广义相对论检验

太阳引力对星光(和射电波)的偏折(Shapiro 延迟)可以通过精确测量太阳附近类星体位置来检验。VLBI 测量与广义相对论预言的相符程度已达到 10−410^{-4} 量级(Fomalont & Sramek,2009)。

太阳附近的射电源过境时,引力造成的偏折和夏皮罗时延必须从几何模型里扣掉,否则微角秒天体测量会系统偏移。把偏折系数留作自由参数,拟合值与广义相对论预言的相对偏差已被压到文中的 10−410^{-4} 量级(Fomalont & Sramek,2009)。这是把 检验 做成日常改正之后,再反过来读残差。

VLBI 在脉冲星天体测量与引力检验中的应用

VLBI 对脉冲星的高精度位置和自行测量(通常精度约 $10$–100 μ100 \, \muas)是天体测量领域最重要的应用之一。通过对多颗毫秒脉冲星的长期 VLBI 监测,可以:

  • 测量脉冲星的三角视差(精度约 $0.1$–1 mas1 \, \text{mas}),获得精确距离,从而校准脉冲星光度和色散量-距离关系(DM 校准,Deller 等,2016,ApJ 817,16)。
  • 测量脉冲星的自行(精度约 $0.01$–1 mas/yr1 \, \text{mas/yr}),追踪超新星遗迹喷出动力学。
  • 通过毫秒脉冲星与伴星双星系统的 Shapiro 时延测量(用于约束脉冲星质量)结合 VLBI 自行,对 Shapiro 时延的预言进行广义相对论检验,精度可达 10−310^{-3} 量级。

VLBI 与脉冲星计时阵(PTA,Pulsar Timing Array)的结合特别有价值:PTA 通过高精度监测数十颗毫秒脉冲星的到达时间,寻找纳赫兹引力波背景(纳赫兹 GWB)在脉冲星时序中的共同相关签名。2023 年,北美纳赫兹引力波天文台(NANOGrav)、欧洲脉冲星计时阵(EPTA)、印度脉冲星计时阵(InPTA)和帕克斯脉冲星计时阵(PPTA)联合公布了约 $3$–4σ4\sigma 的纳赫兹 GWB 证据(NANOGrav Collaboration,2023,ApJL 951,L8),可能来自超大质量黑洞双星并合的引力波背景。VLBI 为 PTA 提供了脉冲星精确距离,帮助建立更好的脉冲星噪声模型。

脉冲星计时给出的是到达时间,VLBI 给出的是天球位置与自行。两者合在一起,视差距离才不单靠色散量模型。中子星与脉冲星 的光度、超新星踢射速度,都对距离尺度敏感。脉冲星计时阵 要在纳赫兹引力波分析里建立噪声模型,也需要这些几何距离,而不是把色散量直接当成尺子。

VLBI 与活动星系核喷流的超光速运动

射电星系和耀变体(blazars)的 VLBI 分辨成像是研究 AGN 喷流动力学的核心工具。在多历元 VLBI 图像中,喷流成分(knots)的表观横向运动速度常常超过光速——这种"超光速运动"是相对论喷流以接近光速指向观测者时产生的几何投影效应,而非真正违反相对论。表观速度:

βapp=βsin⁡θ1−βcos⁡θ\beta_{\rm app} = \frac{\beta \sin\theta}{1 - \beta \cos\theta}

对于 β→1\beta \to 1 和小观测角 θ\theta,βapp\beta_{\rm app} 可达数十甚至更大。VLBA 的 MOJAVE 程序(Monitoring Of Jets in Active galactic nuclei with VLBA Experiments)系统性地监测约 300 个 AGN 的厘米波喷流,测量了其运动学和磁场结构(Lister 等,2016,AJ 152,12)。这些数据是约束喷流速度、洛伦兹因子和 AGN 功率的唯一直接手段。

表观超光速不是信号在真空里超过 $c$,而是接近光速的喷流以小角度指向观测者时的投影。MOJAVE 用 VLBA 在厘米波多年监测同一批核,为的是把运动学从单历元的“看起来很快”变成可重复的速度场(Lister 等,2016)。没有这条地球基线,喷流核成分会糊成一个点。

未来方向:空间 VLBI 与亚毫米 VLBI

空间 VLBI 将天线放入轨道以获得超越地球尺度的基线。VSOP(1997—2003,日本 HALCA 卫星)和 RadioAstron(2011—2019,俄罗斯,最大基线 350000 km)已证明技术可行,角分辨率达约 7 μ7 \, \muas。未来的 ngVLA(下一代 VLA)和 EHT 扩展将进一步提升分辨率和灵敏度。

在 0.87 mm 波长工作的 EHT 下一阶段目标分辨率可能突破 10 μ10 \, \muas,能够对黑洞事件视界附近的磁场结构成像,检验广义相对论在强场极限下的预言。

空间 VLBI 把 $B$ 拉到地球轨道以外,公式仍是 λ/B\lambda/B,难处在星上氢频标、下链带宽和与地面站的相位传递。RadioAstron 已经证明条纹可以在三十多万公里基线上找到。下一步若要做黑洞视界的电影,缺的往往是短基线对大尺度结构的填充,而不只是把最长基线再加长一点。

跨域连接

  • 事件视界望远镜:把基线拉到地球直径,毫米波的衍射极限就进到微角秒,刚好够分辨近邻超大质量黑洞的视界尺度。代价是台站极少、频域采样极稀——所以那张"黑洞照片"是带正则化的重建结果,而不是直接曝光的图像。
  • 分布式系统:基线太长,各站无法实时互联,数据只能先写盘再运回相关中心。于是系统退化为没有共享时钟的分布式测量:每站靠氢原子频标独立守时,日漂移量决定了可相干积分的时长。时钟一散相,再多灵敏度也换不回条纹。
  • 最优化:可见度样本远少于待求的像素,反演天生欠定,解不唯一。可行做法是加先验——稀疏性、正性、最大熵——把病态问题正则化。推论:结论必须对先验的选择做敏感性检验,否则"看到的结构"可能来自算法而不是来自天体。
  • 大地测量与大地水准面:反过来用同一批数据:把遥远类星体当作方向固定的参考点,测出的时延差就变成台站之间的基线变化,年积累量到毫米级。这使干涉阵成为板块运动与地球自转不均匀性的基准仪器,天文与大地测量共用一套解算。
  • 广义相对论的实验检验:太阳引力造成的射电波偏折与信号时延,在微角秒精度上是必须逐项扣除的常规改正,而非待检验的效应。这个身份翻转本身就是一种检验:把偏折系数当自由参数拟合,它与相对论预言的偏离已被压到极小量级。

参考文献

  1. Thompson, A. R., Moran, J. M., & Swenson, G. W. (2017). Interferometry and Synthesis in Radio Astronomy (3rd ed.). Springer Open. doi:10.1007/978-3-319-44431-4
  2. Event Horizon Telescope Collaboration (2019). First M87 Event Horizon Telescope Results. I-VI. ApJL, 875, L1–L6.
  3. Charlot, P. et al. (2020). The third realization of the International Celestial Reference Frame by very long baseline interferometry. A&A, 644, A159.
  4. Reid, M. J., & Honma, M. (2014). Microarcsecond Radio Astrometry. ARA&A, 52, 339–372. doi:10.1146/annurev-astro-081913-040006
  5. Deller, A. T. et al. (2011). DiFX-2: A More Flexible, Efficient, Robust, and Powerful Software Correlator. PASP, 123, 275–287. doi:10.1086/658907
  6. Event Horizon Telescope Collaboration (2022). First Sagittarius A\* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way. ApJL, 930, L12.
  7. Napier, P. J. et al. (1994). The Very Long Baseline Array. Proceedings of the IEEE, 82, 658–672.

延伸阅读

  1. NRAO VLBA 官方页面:https://public.nrao.edu/telescopes/vlba/