基于Huygens波前合成的高分辨率成像大气射电望远镜模拟及其在宇宙伽马射线与中微子观测中的应用

《Astroparticle Physics》：Simulating an imaging atmospheric radio telescope to observe cosmic gamma rays and cosmic neutrinos

　　

对能量高于10¹²电子伏特（TeV）范围的宇宙伽马射线的观测，往往是对理论预测受到极端条件挑战的事件的直接探测手段。由于高能宇宙伽马射线稀少，需要面积约为十几个足球场大小（10⁵ m²）或更大的收集面积，才能在短时间内积累足够的计数统计量。成像大气Cherenkov望远镜（IACT）能够在如此大的面积上探测宇宙伽马射线，目前在宇宙粒子的类型、能量和方向重建能力方面表现最佳。然而，每年仅有约1000小时的黑暗晴朗夜晚，成像大气Cherenkov望远镜的占空比（duty cycle）仅为约12%。为克服这一限制，研究人员提出了一种新型成像大气射电望远镜（imaging atmospheric radio telescope）。通过观测大气簇射在10 GHz频段发出的射电辐射，成像大气射电望远镜可以全天候（24/7）观测，不受昼夜循环影响，且几乎不受天气条件限制。其创新之处在于对大气簇射射电辐射进行高分辨率成像，有望获得与成像大气Cherenkov望远镜类似的高分辨率图像。研究人员提出了一种利用波动力学而非传统射线追踪来模拟射电望远镜中图像形成过程的技术。成像大气射电望远镜可将高能伽马射线观测的占空比提高八倍以上。此外，成像大气射电望远镜还可能将所谓“掠地”（Earth-skimming）技术对宇宙中微子的观测能力提高八倍。

研究背景方面，对宇宙伽马射线的观测为理解宇宙中最剧烈的过程提供了独特视角。伽马射线常产生于宇宙尺度的瞬变灾变事件中，最适合用能够及时积累统计量的探测器进行观测。此外，大多数伽马射线发射呈现幂律能谱，高能事件远较低能事件稀少。综合来看，在TeV及更高能区观测伽马射线天空，需要具有大收集面积的探测器。目前，宇宙伽马射线可直接在太空探测，也可借助地球大气间接探测。太空探测器虽然重建能力强、对伽马射线的纯度极高，但尽管成本高昂，其收集面积仅约1 m²，限制了对高能天空的探测能力。间接大气探测方法中，成像大气Cherenkov望远镜（IACT）以高分辨率（约0.1°）对大气簇射成像，能量阈值可低至约10^-2 TeV，是继太空探测器之后能收集最纯净宇宙伽马射线样本的技术，但其占空比仅约12%。大气射电技术则具有接近100%的占空比，此前的射电成像原型实验（AMBER、CROME、MIDAS）已在1–10 GHz频段对大气簇射射电辐射进行方向分箱成像，但能量阈值较高（10⁴–10⁶ TeV），不利于伽马射线观测。鉴于成像在Cherenkov技术中的强大能力，研究人员提出探索具有足够角分辨率的高分辨率成像大气射电望远镜的前景，以受益于图像空间中强大的重建技术。研究旨在通过实现一个现实可行的系统来探索高分辨率成像大气射电望远镜的前景。

研究人员开展了成像大气射电望远镜的模拟研究，提出了基于Huygens原理的图像形成方法，并设计了三种不同尺寸的望远镜（??rome、??edium、??arge）进行系统评估。研究得出了以下结论：Huygens图像形成方法能够自然再现波动光学（Airy斑）和几何光学（像差、畸变、景深变窄）的预测特征；能够正确分离同时到达的多个平面波信号；能够正确保持大气簇射的方向信息并再现粒子依赖的图像形态；频率依赖成像可区分簇射几何结构；能量阈值估计在1–2 GHz频段约为70–100 TeV，并随频率升高而增加。这些结果证明了利用大气射电望远镜对广延大气簇射进行成像的可行性，为未来纳入噪声、探测器性能和定量重建精度的研究奠定了基础。其重要意义在于：成像大气射电望远镜有望将高能伽马射线观测的占空比提高八倍以上，并可能将掠地技术对宇宙中微子的观测能力提高八倍，从而显著拓展高能宇宙粒子观测的科学 reach。

在技术方法方面，研究人员主要采用了以下关键技术方法：利用CORSIKA-CoREAS程序模拟广延大气簇射的射电辐射电场输出（观测海拔110 m，磁场配置参照CROME站点）；采用基于Huygens-Fresnel原理的图像形成方法，将望远镜镜面和相机屏幕表示为Huygens散射中心点云，通过传播复电磁场并叠加干涉、衍射、相位延迟和部分相干效应来模拟成像；通过平面波包能量校准确定全局能量守恒因子；利用理想带通滤波器进行频率分辨成像；基于系统温度（??_sys=50、70、100 K）估算背景噪声等效电场并推算能量阈值。样本队列来源为CORSIKA-CoREAS模拟的伽马射线、质子和铁核引发的大气簇射，未涉及实际实验队列。

研究结果方面：（1）Huygens图像形成结果：通过与波动光学远场近似（Airy极限）和几何光学射线追踪对比，发现Huygens方法在光轴上接近Airy斑，在离轴处收敛于射线追踪曲线，中间角度区域同时体现衍射和像差的贡献，这是仅有完整Huygens方法才能捕捉的区域；三种望远镜均再现了能量守恒趋势，??edium和??arge的点扩散函数接近Airy斑极限（PSF_b≈1.0），畸变分别为约4%和2%，??rome畸变约28%。（2）同时观测两个平面波：??edium和??arge望远镜均能清晰区分两个同时到达的平面波包，在能量和频率上均可分辨，表明Huygens图像形成可正确处理大气簇射模拟中多粒子相互作用发射的叠加。（3）离焦观测：??arge望远镜离焦图像中可清晰看到镜面环状结构，但发现弥散的环状伪影和能量涨落。（4）伪影研究：当平面波来自视场外时，大多数相机屏幕散射中心包含约等于平面波能量除以镜面散射中心数??的伪影，伪影模式的出现显著依赖于镜面散射中心的排列方式，方形网格排列即使加入随机性也会产生显著伪影（鬼影）。（5）大气簇射图像验证：不同来向的1 TeV伽马射线引发簇射的相机图像随来向系统性位移，验证了Huygens图像形成正确保持方向信息。（6）Cherenkov锥处粒子依赖图像：在距簇射轴约112 m处，伽马射线、质子和铁核引发的簇射图像呈现特征差异——给定初级粒子类型，信号幅度随能量成比例增加；固定能量下，信号强度随核电荷增加而降低；高能下不同粒子类型的图像形态差异减弱，产生能量与成分之间的部分简并。（7）频率依赖成像：1–2 GHz频段发射在大部分簇射中相干，图像宽阔平滑主要反映总能量；2–4 GHz频段大尺度相干性降低，图像更局域化并揭示簇射方向；4–8 GHz频段相干性仅维持在小区域，发射沿投影簇射轴集中；8–12 GHz频段相干性大部分丧失，锐利分辨簇射核心但幅度下降，体现空间分辨率与可探测性之间的权衡。（8）背景噪声与能量阈值估计：最低阈值在1–2 GHz频段约为70–100 TeV（??_sys=50–100 K），2–4 GHz约0.2–0.4 PeV，4–8 GHz约几PeV，8–12 GHz约数十PeV；拟合幂律指数??在低频接近1并向高频递减。（9）计算性能：??arge望远镜配置下，望远镜响应计算约12分钟且与能量无关，而CORSIKA-CoREAS大气簇射模拟随能量显著增加（1 TeV约4分钟，10 TeV约35分钟，50 TeV近3小时），端到端总运行时间从1 TeV约16分钟增至50 TeV约3.1小时。

讨论部分总结：Huygens图像形成方法仍有改进空间，存在与预期的偏差，能量守恒较为微妙，仿真分辨率与计算量之间的平衡具有挑战性。但该方法自然再现了实际望远镜的Airy斑、像差、畸变和景深变窄等特征和限制，同时平面波观测表明复杂大气簇射成像可能可行。模拟射电图像证明Huygens图像形成方法在整个模拟链中保持了广延大气簇射的关键特征，伽马射线和强子引发簇射图像与IACT已知图像相似，表明已建立的IACT重建技术可能可迁移至射电成像领域。频率依赖成像进一步实现了簇射几何的解缠，凸显了多波段射电观测的潜力。能量阈值结果呈现强频率依赖性，但仅为基于单通道的估计，未包含天线响应且未优化带宽。展望部分指出：纳入真实接收链及其噪声贡献、扩展至不同天顶角簇射、覆盖更广泛能量和粒子类型以系统研究伽马-强子分离、探索不同望远镜设计、系统研究望远镜相对Cherenkov环的位置、探索掠地tau中微子探测（PeV能区）等是未来研究方向。研究结论翻译：模拟的广延大气簇射射电图像表明，Huygens图像形成方法在整个模拟链中保持了广延大气簇射的关键特征。方向信息和粒子依赖的形态被正确再现，为理想无噪声条件下成像大气射电望远镜提供了清晰的概念验证。伽马射线和强子引发簇射图像与成像大气Cherenkov望远镜已知图像相似，表明已建立的IACT重建技术可能可迁移至射电成像领域。频率依赖成像进一步实现了簇射几何的解缠，凸显了多波段射电观测用于详细大气簇射研究的潜力。能量阈值初步结果呈现强频率依赖性，范围从70 TeV到数十PeV。这些值应仅视为估计，因其基于单通道，未包含天线响应且未进行带宽优化。这些结果证明了利用大气射电望远镜对广延大气簇射进行成像的前景，并推动了未来纳入噪声、探测器性能和定量重建精度的研究。