Synthesis imaging with a lunar orbit array: I. global sky map and its systematics

技术摘要：月球轨道阵列的合成成像：I. 全天图及其系统误差

问题陈述
低于 30 MHz 的地面射电天文学受到电离层吸收和地面射电频率干扰（RFI）的严重限制。发现长波天体（DSL）任务（亦称为“鸿蒙”任务）提议使用一个绕月运行的线性卫星阵列来克服这些限制。虽然此类阵列收集的可见度数据与天空图呈线性相关，但重建全球天空图面临着不同于地面干涉测量技术的特定算法挑战。这些挑战包括由于月球遮挡部分天空而导致的 $uvw$ 覆盖不全、解决病态反演问题所需的正则化，以及由天空有限像素化引起的系统误差。具体而言，本文研究了亚像素噪声（即在像素尺度以下波动的噪声）如果在构建波束矩阵时未得到妥善处理，将如何诱发显著的混叠效应。

方法论
作者采用数值模拟来评估 DSL 阵列配置下的全球图重建算法。

• 天空模型： 为了解决现有全球天空模型（如 GSM、SSM）缺乏小尺度功率的问题，作者构建了一个模拟全天图，结合了高分辨率银河系弥散辐射模型（基于增强型 ULSA 模型）和源自 GLEAM 巡天计数的点源种群。该图生成于 $N_{side} = 2048$，并降级至 $N_{side} = 1024$ 作为模拟输入。
• 阵列配置： 模拟了一个由 8 颗子卫星组成的圆形月球轨道阵列（高度 300 km，倾角 30°），其周期为 1.3 年。卫星采用“呼吸”策略，通过定期压缩和拉伸其构型来改变基线长度。基线分布经过优化，范围从 100 m 到 100 km。
• 信号建模： 可见度数据被建模为天空强度的线性映射，并纳入了热噪声、有限时间分辨率和有限带宽效应。由于延迟相对于基线长度的关系，月球反射被认为对相干性影响可以忽略不计，尽管它略微增加了源的通量。
• 重建算法： 作者使用 Tikhornov 正则化求解线性反演问题。他们比较了两种构建波束矩阵（$B$）的方法：
  1. 像素中心法： 在每个像素的中心评估波束响应。
  2. 像素平均法： 在更高分辨率下评估响应，并向下采样至目标分辨率以平滑亚像素变化。
• 指标： 重建质量通过输入图与重建图之间的相关系数 ($\rho_l$)、球面谐波空间的信噪比 (SNR) 以及热噪声和不完美波束效应误差协方差矩阵的迹来进行量化。

主要贡献与结果

1. 通过像素平均缓解混叠： 研究表明，若天真地使用像素中心法构建波束矩阵，当包含长于 Nyquist 采样极限（$b > \lambda/\theta_p$）的基线时，会导致严重的混叠伪影。这些伪影表现为显著的畸变和相关系数的降低。作者展示了像素平均法能有效缓解这种混叠，从而允许在不降低图像质量的前提下包含更长的基线。
2. 基线选择： 分析显示，包含高达 $b < 2b_p$（其中 $b_p$ 是对应于像素尺寸的基线）的基线对于重建已足够。包含更长的基线（$b > 2b_p$）对最终图像的提升微乎其微，因为平均过程会将其贡献衰减至接近于零。
3. 正则化参数权衡： 本文研究了正则化参数 ($\epsilon$) 的选择。在抑制热噪声（有利于较大的 $\epsilon$）与最小化波束缺陷误差（有利于较小的 $\epsilon$）之间存在权衡。最优 $\epsilon$ 位于两条误差曲线的交点处，尽管作者指出这取决于具体情况。
4. 极区系统误差： 由于轨道倾角导致的短投影基线缺失，重建结果在极区表现出系统性的负偏差（变暗）。这种效应在较高频率（例如 10 MHz）下更为显著，因为在波长单位下，投影基线实际上变得更长。使用先验图可以缓解这种功率损失，但会引入对先验准确性的敏感性。
5. 频率依赖性： 与预期中较低频率因更高系统温度而产生更多噪声的结论相反，模拟显示在 3 MHz 时，强源的信噪比实际上优于 10 MHz。这归功于更强的源具有更高的亮度温度，以及重建过程中包含了更多的基线（由于波长更长）。

意义与主张
本文声称为 DSL 任务的合成成像提供了一个稳健的算法框架。其主要意义在于识别并解决了亚像素混叠问题，这是月球轨道阵列成像中一个关键的系统误差，且不同于标准的地面技术。通过验证像素平均法，作者证明了尽管存在独特的 $uvw$ 覆盖和轨道约束，仍可以重建高质量的全球天空图。作者明确指出，其结果是基于理想化模型的；他们承认诸如仪器定标、卫星位置/取向误差以及月球反射效应等系统误差在本文中被忽略了，这些将在后续论文中进行讨论。作者断言，只要卫星阵列满足其设计参数，所呈现的结论对于此类系统误差是稳健的。