甚长基线干涉测量
概述
甚长基线干涉测量(Very Long Baseline Interferometry,VLBI)是一种将分布在地球各地(乃至地球轨道上)的射电天线联合起来,形成等效口径等于天线间距的超大型干涉仪的技术。VLBI 的角分辨率可达微角秒(as)量级,是目前人类能够实现的最高天文分辨率——远超任何光学望远镜。它同时也是测量板块运动、地球自转和建立国际天球参考系的重要大地测量工具。
把地球直径当成口径,分辨率可以进到微角秒,但各站之间没有共用的本振电缆。这与甚大阵一类的连线干涉不同:后者能在本地把信号相干相加,VLBI 只能先把电压时间序列连同氢频标时间戳写进磁盘,再在相关机上做事后互相关。分辨率由最长基线决定,灵敏度由天线总面积和系统温度决定,两者不能互相代替。
因此 VLBI 的“地球尺度望远镜”是稀疏孔径。台站一少,空间频率覆盖就出现空洞,图像必须从可见度重建,而不是在焦平面上直接曝光。后文的 VLBA、EVN 与 事件视界望远镜,都是在这条约束下工作。
干涉测量原理:时延与互相关
两台相距基线长度 $B$ 的射电天线接收同一天体信号时,由于光程差的存在,信号到达两台天线有时间延迟 :
其中 为天线基线方向与天体方向的夹角,$c$ 为光速。通过对两路信号进行互相关(记录到高精度时频标准——氢频标——后在相关机上离线计算),可以得到干涉条纹,干涉仪的角分辨率约为:
VLBI 基线 $B$ 最大可达地球直径(约 12700 km),在 1.3 mm 波长(事件视界望远镜工作波段)的衍射极限约为 as(微角秒),约为人眼分辨率的一千万分之一。
几何时延是科学量:它把基线矢量投影到视线方向上,源的位置稍有变化, 就变。实际记录里还叠着钟差、大气与仪器时延,相关机必须在时延—时延率平面上搜相关峰,找到条纹以后再把这些干扰项一并解出。互相关给出的是复可见度,不是一张现成的像。
两路信号要对齐到远小于一个波长的等效光程,条纹才会积起来。氢频标把各站的时间轴钉在同一套物理秒上,互相关才有意义。基线越长、波长越短,对时钟和大气的要求越苛刻,这也是毫米波全球阵比厘米波阵更难的原因。
VLBI 的技术特点
独立时钟与氢频标
VLBI 各站不能实时连接传输数据(基线过长),信号必须记录在磁盘阵列上后运回相关中心。因此每台天线必须配备极高精度的独立时频标准——氢原子频标(hydrogen maser),日积累误差约 – 秒,确保各站记录的时间戳足够一致以进行事后相关。
氢原子频标锁定在氢基态超精细跃迁附近,为的是提供可独立运行的时间基准,而不是为了观测 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 个致密射电源,参考方向不确定度约 as(Charlot 等,2020,A&A 644,A159)。
选遥远致密射电源,是因为它们在天球上几乎不显示可测的周年视差,适合当方向基准。活动星系核的喷流会在毫角秒上缓慢移动,所以参考架要不断剔除结构变化大的源,ICRF3 才把不确定度压到文中给出的微角秒量级(Charlot 等,2020)。光学星表要和它对齐,才能把空间卫星的姿态接到同一套惯性方向上。
板块运动与地球自转
各站 VLBI 测量的射电源方向差异积累,可以精确测定各台站的相对位置变化,从而确定板块漂移速率(精度约 1 mm/yr)和地球自转不规则性(UT1-UTC 精度约 s)。VLBI 数据是国际地球自转与参考系服务(IERS)维护全球时间和参考系的核心输入。
同一套时延,既是天体的方向,也是台站的坐标。把源的位置当已知,解出来的就是基线矢量随时间的变化,年速率到毫米。国际地球自转与参考系服务把 VLBI 的 UT1 与台站坐标,和卫星测地、月球激光测距并在一起维护。天文与 大地测量 在这里用的是同一类相关峰。
广义相对论检验
太阳引力对星光(和射电波)的偏折(Shapiro 延迟)可以通过精确测量太阳附近类星体位置来检验。VLBI 测量与广义相对论预言的相符程度已达到 量级(Fomalont & Sramek,2009)。
太阳附近的射电源过境时,引力造成的偏折和夏皮罗时延必须从几何模型里扣掉,否则微角秒天体测量会系统偏移。把偏折系数留作自由参数,拟合值与广义相对论预言的相对偏差已被压到文中的 量级(Fomalont & Sramek,2009)。这是把 检验 做成日常改正之后,再反过来读残差。
VLBI 在脉冲星天体测量与引力检验中的应用
VLBI 对脉冲星的高精度位置和自行测量(通常精度约 $10$–as)是天体测量领域最重要的应用之一。通过对多颗毫秒脉冲星的长期 VLBI 监测,可以:
- 测量脉冲星的三角视差(精度约 $0.1$–),获得精确距离,从而校准脉冲星光度和色散量-距离关系(DM 校准,Deller 等,2016,ApJ 817,16)。
- 测量脉冲星的自行(精度约 $0.01$–),追踪超新星遗迹喷出动力学。
- 通过毫秒脉冲星与伴星双星系统的 Shapiro 时延测量(用于约束脉冲星质量)结合 VLBI 自行,对 Shapiro 时延的预言进行广义相对论检验,精度可达 量级。
VLBI 与脉冲星计时阵(PTA,Pulsar Timing Array)的结合特别有价值:PTA 通过高精度监测数十颗毫秒脉冲星的到达时间,寻找纳赫兹引力波背景(纳赫兹 GWB)在脉冲星时序中的共同相关签名。2023 年,北美纳赫兹引力波天文台(NANOGrav)、欧洲脉冲星计时阵(EPTA)、印度脉冲星计时阵(InPTA)和帕克斯脉冲星计时阵(PPTA)联合公布了约 $3$– 的纳赫兹 GWB 证据(NANOGrav Collaboration,2023,ApJL 951,L8),可能来自超大质量黑洞双星并合的引力波背景。VLBI 为 PTA 提供了脉冲星精确距离,帮助建立更好的脉冲星噪声模型。
脉冲星计时给出的是到达时间,VLBI 给出的是天球位置与自行。两者合在一起,视差距离才不单靠色散量模型。中子星与脉冲星 的光度、超新星踢射速度,都对距离尺度敏感。脉冲星计时阵 要在纳赫兹引力波分析里建立噪声模型,也需要这些几何距离,而不是把色散量直接当成尺子。
VLBI 与活动星系核喷流的超光速运动
射电星系和耀变体(blazars)的 VLBI 分辨成像是研究 AGN 喷流动力学的核心工具。在多历元 VLBI 图像中,喷流成分(knots)的表观横向运动速度常常超过光速——这种"超光速运动"是相对论喷流以接近光速指向观测者时产生的几何投影效应,而非真正违反相对论。表观速度:
对于 和小观测角 , 可达数十甚至更大。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)已证明技术可行,角分辨率达约 as。未来的 ngVLA(下一代 VLA)和 EHT 扩展将进一步提升分辨率和灵敏度。
在 0.87 mm 波长工作的 EHT 下一阶段目标分辨率可能突破 as,能够对黑洞事件视界附近的磁场结构成像,检验广义相对论在强场极限下的预言。
空间 VLBI 把 $B$ 拉到地球轨道以外,公式仍是 ,难处在星上氢频标、下链带宽和与地面站的相位传递。RadioAstron 已经证明条纹可以在三十多万公里基线上找到。下一步若要做黑洞视界的电影,缺的往往是短基线对大尺度结构的填充,而不只是把最长基线再加长一点。
跨域连接
- 事件视界望远镜:把基线拉到地球直径,毫米波的衍射极限就进到微角秒,刚好够分辨近邻超大质量黑洞的视界尺度。代价是台站极少、频域采样极稀——所以那张"黑洞照片"是带正则化的重建结果,而不是直接曝光的图像。
- 分布式系统:基线太长,各站无法实时互联,数据只能先写盘再运回相关中心。于是系统退化为没有共享时钟的分布式测量:每站靠氢原子频标独立守时,日漂移量决定了可相干积分的时长。时钟一散相,再多灵敏度也换不回条纹。
- 最优化:可见度样本远少于待求的像素,反演天生欠定,解不唯一。可行做法是加先验——稀疏性、正性、最大熵——把病态问题正则化。推论:结论必须对先验的选择做敏感性检验,否则"看到的结构"可能来自算法而不是来自天体。
- 大地测量与大地水准面:反过来用同一批数据:把遥远类星体当作方向固定的参考点,测出的时延差就变成台站之间的基线变化,年积累量到毫米级。这使干涉阵成为板块运动与地球自转不均匀性的基准仪器,天文与大地测量共用一套解算。
- 广义相对论的实验检验:太阳引力造成的射电波偏折与信号时延,在微角秒精度上是必须逐项扣除的常规改正,而非待检验的效应。这个身份翻转本身就是一种检验:把偏折系数当自由参数拟合,它与相对论预言的偏离已被压到极小量级。
参考文献
- 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
- Event Horizon Telescope Collaboration (2019). First M87 Event Horizon Telescope Results. I-VI. ApJL, 875, L1–L6.
- Charlot, P. et al. (2020). The third realization of the International Celestial Reference Frame by very long baseline interferometry. A&A, 644, A159.
- Reid, M. J., & Honma, M. (2014). Microarcsecond Radio Astrometry. ARA&A, 52, 339–372. doi:10.1146/annurev-astro-081913-040006
- 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
- 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.
- Napier, P. J. et al. (1994). The Very Long Baseline Array. Proceedings of the IEEE, 82, 658–672.
延伸阅读
- NRAO VLBA 官方页面:https://public.nrao.edu/telescopes/vlba/