事件视界望远镜合作组织(Event Horizon Telescope Collaboration, EHTC)利用毫米波干涉测量技术合成了地球大小的有效孔径,并发布了位于射电星系Messier 87(M87)核心的黑洞阴影的首张重建图像,标志着黑洞事件视界附近吸积动力学直接可视化新纪元的开启。该图像呈现为环状结构,代表围绕星系核中超大质量黑洞M87*的引力透镜效应和红移效应下的辐射。然而,应用不同的数据分析与成像启发式方法导致了对M87*精细尺度特征的分歧性表征,部分外部分析甚至对源的 fundamental 环状形态提出了质疑。这些困难源于甚长基线干涉测量(Very Long Baseline Interferometry, VLBI)数据固有的稀疏性,使得干涉图像重建成为一个约束不足的问题。此外,解对仪器校准、后处理算法和成像超参数高度敏感。本研究将研究人员独立开发的成像算法应用于M87*数据集,该算法不依赖规定的形态学先验,且建立在不受仪器校准误差影响的观测量基础之上。研究人员重建了黑洞阴影图像,明确再现了环状形态,其角直径为42.2 ± 1.5 μas,在假定标称标度参数下对应黑洞质量为6.5 ± 0.5 × 109 M⊙,观测倾角约为20°(相对于接近喷流方向)。
研究背景方面,M87*超大质量黑洞位于距离D = 16.8 Mpc处,其质量约为M
BH = (3–7) × 10
9 M
⊙,为研究视界尺度广义相对论物理提供了独特窗口。根据广义相对论预言,其在事件视界尺度上的外观应由中心阴影主导,对于非旋转史瓦西黑洞,阴影特征由临界碰撞参数描述,该参数与黑洞质量成正比,即3√3GM
BH/(c
2D),其中G为引力常数,c为光速。要观测M87*黑洞阴影的细节,需要前所未有的微角秒(μas)级角分辨率,大致相当于1000公里外人发丝宽度所张的角度。实现这一目标需要结合毫米波段(λ ≈ 1 mm)和与地球直径相当的有效孔径。事件视界望远镜(EHT)通过甚长基线干涉测量(VLBI)技术,关联全球独立射电天文台网络的信号,实现了地球大小孔径的合成。然而,由此产生的虚拟孔径仍大量未填充,缺乏填充孔径的灵敏度和重建保真度。数据稀疏性问题因大陆尺度上快速、基于站点的 атмосферный 相位和振幅波动的校准需求而进一步加剧。这些固有挑战使得VLBI和EHT数据集的重建成为一个不适定问题,对仪器校准和后处理启发式方法高度敏感。现有算法通常使用手动调谐和独立加权的数据项与正则化来偏置解 towards 合理形态,可能无意中将重建导向偏好解,导致M87*精细尺度特征(如表面亮度分布和热点定位)的分歧性表征。使用传统方法(包括基于CLEAN的算法或结合自校准的混合映射),独立分析表明重建结果对方法学和超参数选择敏感,部分独立分析甚至展示了没有清晰环状结构的重建结果。为解决围绕M87*结果的争议并超越方法学不确定性,研究人员从传统方法转向更自主和客观的框架。该方法的独立性主要基于使用稳健的干涉观测量——闭合不变量(closure invariants),其捕获的源形态信息不受波前和仪器校准误差的影响。闭合不变量对绝对位置和绝对流量标度不敏感,但這些并非固有形态参数。与传统的闭合相位和闭合振幅不同,研究人员使用的闭合不变量形式提供了一个统一框架,携带等效信息同时消除了文献中常见的异质处理需求。通过利用闭合不变量携带的丰富校准无关的纯形态信息,可以绕过需要敏感校准调谐和观察者驱动先验的传统技术。
研究人员开展了以下研究:开发并应用了生成式深度学习图像重建与闭合项(Generative Deep learning Image Reconstruction with Closure Terms, GenDIReCT)算法,该算法在非天文图像上训练,设计为仅使用闭合不变量生成数据驱动的图像重建。GenDIReCT使用闭合不变量作为条件来调节潜在扩散模型(latent diffusion model)。该算法已在合成和真实数据上进行了广泛测试,证明能够在无需手动干预、规定性先验或参数调谐的情况下一致地恢复多样形态结构。研究首次使用GenDIReCT算法通过闭合不变量对M87*进行成像,使用公开可用的1.3毫米EHT测量数据,涵盖2017年四个观测历元(4月5日、6日、10日、11日)、2018年两个历元(4月21日、25日)和2021年两个历元(4月13日、18日)。设置了多个视场用于重建(128×128、225×225和450×450 μas
2),每个数据集和视场的组合是一个独立模型,并作为最大源尺寸的间接先验。对于每个模型,将大量生成式重建聚合为单一代表中值图像,并与使用公开eht-imaging流程重现的EHTC基准重建进行比较。还从每个历元所有视场的所有重建的逐像素中值导出合成图像,对齐到具有最佳数据保真度得分的单一重建。
关键技术方法方面,本研究主要采用以下方法:(1)闭合不变量(closure invariants)形式体系,从传统闭合相位和闭合振幅推广而来,提供校准无关的干涉观测量;(2)GenDIReCT算法,一种基于潜在扩散模型的条件生成式深度学习图像重建方法,使用闭合不变量作为唯一数据项;(3)贝叶斯几何建模,使用xsringauss模型和DYNESTY嵌套采样算法对闭合不变量进行前向建模;(4)图像域特征提取方法,用于定量测量环直径、宽度、方向、方位不对称性和中心亮度分数等参数。数据来源为2017年、2018年和2021年EHT观测的公开网络校准数据集,经EHT-HOPS流程处理。
研究结果方面,首先,所有独立观测历元的数据集均令人信服地产生约40 μas尺度的紧凑环状结构。在每年最高质量数据集(2017年4月11日、2018年4月21日、2021年4月18日)上,中值重建为紧凑的圆形化环,具有一致的方向和直径,独立于重建视场窗口,并在南方或西南方向增亮。环的外观与从接近喷流方向观测到的约20°倾角最为一致。其次,源尺寸(约40 μas)跨所有视场的一致性强烈确认了紧凑源。在最受约束的128×128 μas
2视场下,中值重建图像与EHTC基准重建展现出最高的图像域相似度得分,且数据 adherence 指标性能与较大视场具有竞争力,支持了M87*喷流延展辐射不受EHT数据集约束的先验预期。在某些质量较低的数据集上,环状形态可能变得不太清晰,这可能是孔径覆盖不足或低质量数据的伪影。当考虑更高质量数据集时,对视场的敏感性消失,这与基于传统混合映射方法的预测相反。该工作因此产生了对M87*黑洞阴影图像环状几何特征的稳健、独立验证。第三,形态学分析方面,通过测量重建图像的环特征并与已发表参数比较,定量展示了与EHTC图像重建的一致性。测量了环直径、宽度、方向(从北逆时针测量)、方位不对称性和中心亮度分数。GenDIReCT和EHTC成像方法之间的环方向和直径保持相对稳定。2017年至2018年观测之间环方向的演化值得注意,归因于黑洞吸积的湍流性质。2018年和2021年观测到的环方向与M87*的大尺度前向喷流更为对齐。环宽度和中心亮度分数受方法学差异影响较大,可能不反映固有源属性。方位不对称参数在不同日期和方法之间不稳定,表明对孔径覆盖变化的敏感性。第四,通过互补的贝叶斯前向建模框架进一步交叉验证测量环参数,使用具有解析可见度的原始几何形状叠加来拟合观测闭合不变量。使用DYNESTY嵌套采样算法对xsringauss模型进行参数估计。所有历元中,从前向几何建模测得的平均直径和标准差为d = 42.3 ± 1.1 μas,与非参数图像重建的平均d = 42.2 ± 1.5 μas兼容。第五,从测量环直径推断黑洞引力半径及质量需要独立校准成像方法。使用测量环直径、距离测量值16.8
−0.7+0.8 Mpc和引力半径的标称标度参数,估计M87*的黑洞质量为M
BH = (6.5 ± 0.5) × 10
9 M
⊙。独立测量的环直径和推断的黑洞质量与先前的恒星动力学质量测量以及EHTC图像域特征提取结果高度一致。第六,黑洞阴影的外观还对黑洞自旋a
*和观测倾角敏感。将合成M87*重建与IPOLE算法产生的230 GHz合成图像进行比较,该算法在广义相对论磁流体动力学(General Relativistic Magnetohydrodynamic, GRMHD)模拟的流体数据上沿相对论光子轨迹积分协变辐射转移方程。使用归一化互相关相似度得分ρ
NX进行定量分析,识别出对倾角的明显敏感性。最高相似度得分倾向于非零黑洞自旋,但分析未约束自旋大小。最佳得分倾角约为20°,与从43 GHz VLBI观测推断的基准喷流视角17.2° ± 3.3°大体一致。
讨论部分总结了以下要点:与EHTC使用的公开成像流程中的正则化最大似然方法不同,该方法不需要手动指定图像初始化。当模型在不同视场上训练时,模型从其多样化训练库中学习隐式统计先验,以输入干涉观测量(即闭合不变量)为条件。跨不同视场向稳定、单一环结构的收敛表明模型成功地从数据本身推断固有源几何,而非在手动选择的初始化约束内运行。GenDIReCT的自主性和直接从校准无关观测量成像的能力代表了在具有挑战性的数据稀疏VLBI体系中向客观盲成像迈出的显著一步。通过这种方法在多个数据集上明确恢复M87*的环状结构,提供了对已发表EHTC结果的严格交叉验证。研究结论翻译如下:研究人员应用独立的、无形态学先验且基于校准无关观测量的成像算法于M87*数据集,明确再现了环状形态,角直径为42.2 ± 1.5 μas,对应黑洞质量为6.5 ± 0.5 × 10
9 M
⊙,观测倾角约为20°。该结果与先前恒星动力学质量测量及EHTC图像域特征提取高度一致,提供了对M87*黑洞阴影环状几何的稳健独立验证。该方法的自主性和直接从校准无关观测量成像的能力代表了在数据稀疏VLBI体系中向客观盲成像的重要进展。未来空间VLBI的视界尺度成像有潜力解析M87*的光子环,从而实现黑洞自旋的直接精确测定,目前自旋尚未被图像域分析所约束。