LST方向移动平均与校准天空模型精度对宇宙黎明21 cm功率谱测量的影响：SHINONOME II实测与端到端模拟研究

《Publications of the Astronomical Society of Japan》：SHINONOME II: The role of LST-direction moving average and required sky model accuracy

　　

探测宇宙黎明（Cosmic Dawn）的21 cm功率谱（power spectrum）需要极高精度的校准（calibration）和对前景辐射（foreground emission）的显著抑制。本研究分析Murchison Widefield Array（MWA）超低频段（73.6–104.3 MHz）的漂移扫描（drift-scan）观测，并从真实数据表明增益自相关函数（gain autocorrelation functions）的过零时标（zero-crossing timescales）约为～2–3 h，而Allan偏差（Allan-deviation）极小值表明增益幅度（gain amplitude）的特征稳定时标约为～45 min、增益相位（gain phase）约为1.3 h。利用模拟，研究人员进一步评估代表性目录误差——位置误差（positional errors）、流量标度误差（flux-scale errors）和目录不完备性（catalog incompleteness）——如何在校准增益中诱发频率依赖结构（frequency-dependent structures），并在再电离窗口（Epoch of Reionization window, EoR window）内产生系统性功率。结果表明，校准天空目录（calibration sky catalog）所需精度取决于所假设的21 cm信号幅度。对于EDGES所提示的强吸收场景，即使存在～0.1%的流量标度误差和角秒级位置精度，信号仍可被探测。然而，对于标准场景，研究人员发现～0.01%的流量标度精度和～0.′′5的位置精度足以达到所需动态范围（dynamic range）。研究人员还发现，对于强或弱标准21 cm信号场景，源完备性（source completeness）低至1 mJy仍不足以达到必要动态范围，尽管对于EDGES启发的奇异场景（EDGES-inspired exotic scenario）已经足够。受Allan偏差分析所指示的增益稳定时标启发，研究人员在时间上应用增益移动平均（gain moving average）（～50 min），将增益的标准差降低30%–60%，并部分放宽对校准目录的要求。这些结果表明，校准目录的精度是低频21 cm宇宙学的一个关键限制因素，并表明未来的宇宙黎明观测将需要具有更高精度和深度的下一代目录。

研究背景：红移21 cm线来自宇宙黎明和再电离（Epoch of Reionization, EoR）时期的中性氢，可探测首批发光天体的形成以及星际介质的热历史和电离历史。EDGES-low曾报告红移z≈17.8、幅度约500 mK的深吸收特征，但该强信号难以与标准天体物理模型预测相协调，并引发多种非标准或“奇异”场景；同时其可靠性受到系统误差质疑，SARAS3则报告未探测到该特征。因此，独立且互补的测量十分必要，尤其是超越天区平均的空间21 cm涨落测量。低频干涉阵列如MWA、LOFAR、PAPER、GMRT、HERA以及未来SKA-Low正在开展21 cm功率谱统计测量，但尚未实现探测。宇宙黎明观测尤其困难：天体物理前景比预期21 cm信号亮四到五个数量级，干涉仪的色响应将前景污染扩展到(k_⊥,k_∥)平面的楔形区域，仪器色度和校准误差可把前景功率散射到EoR窗口。低频天空目录因波束更宽、角分辨率降低和电离层畸变更不准确；既有校准精度要求多基于约100 MHz以上EoR频段，是否适用于宇宙黎明低频仍不清楚。冗余基线校准虽可利用规则基线从数据中导出增益和真实可见度，但仍需绝对定标并重新依赖天空模型。真实观测中增益并不稳定，若未正确建模会引入额外谱涨落。为此，SHINONOME项目利用MWA Phase II超低频段漂移扫描数据研究电离层、射频干扰、可见度异常值、前景行为及校准与抑制技术。本文聚焦时间增益变化和天空目录精度对校准性能的作用，旨在从真实数据测量增益特征时标，并用端到端模拟评估目录不完美如何传播到校准误差和21 cm功率谱，从而给出宇宙黎明可靠探测所需的天空模型精度与时间稳定性定量要求。其意义在于确认校准目录精度是低频21 cm宇宙学的关键限制因素，并为未来观测的校准策略和下一代目录建设提供依据。

关键技术方法（≤250字）：研究人员使用2023年9月5日MWA Phase II紧凑构型超低频漂移扫描真实数据（73.6–104.3 MHz，z≈12.6–18.3，40 kHz、2 s，5 min分文件，连续约8 h，LST=?3至+5 h）作为增益变化参考，并以同一观测设置开展端到端模拟；模拟用Hyperdrive生成可见度，时间8 s、频率80 kHz，每块37个时间步（296 s），连续10块共2960 s（约49 min），最后一块用于功率谱分析，其余用于移动平均，分析EoR0天区。校准目录取MWA EoR源表（2023年11月7日版），分别引入位置偏移、流量标度扰动和流量阈值截断，并在LST方向对幅度与相位做窗口N=10（约49.3 min）的滑动平均。

研究结果：
3.1 观测增益（Observed gain）：通过真实MWA漂移扫描增益解的ACF和Allan偏差分析，研究人员发现增益幅度和相位的ACF首次过零时标分别约为～2–3 h和～2 h；Allan偏差中值极小值分别出现在τ?2.6×10³ s（～45 min）和τ?4.7×10³ s（～1.3 h）。8 h内中值增益幅度相对观测开始逐渐增加，相位稳定在零附近且无明显系统趋势；增益幅度标准差σ_g在低频更大，并随LST/RA变化，在EoR0附近达到峰值、在RA≈50°的Fornax A附近最小。这些结果为选择约50 min的LST方向移动平均窗口提供经验依据。
3.2 位置偏移导致的增益误差（Gain error from positional offsets）：通过给目录源赤经赤纬独立加入0.′′5、5′′、50′′高斯随机偏移并逐个天线单元计算频率轴增益幅度标准差和相位圆标准差，研究人员发现位置误差增大时幅度和相位散射系统性增大，表明较大偏移诱发更强频率依赖结构；移动平均可抑制带内高频增益涨落。
3.3 流量密度误差导致的增益误差（Gain error from flux-density errors）：通过令每个源流量按均值零、标准差σ=0.01%、0.1%、1%、10%的高斯因子扰动，研究人员发现幅度和相位散射随流量误差增大而增大，说明目录流量不确定度直接表现为增益解误差。
3.4 目录不完备导致的增益误差（Gain error from catalog incompleteness）：通过按束斑衰减后视流量阈值1、10、100 mJy构造不完备目录，研究人员发现阈值越高、目录越浅，幅度和相位散射越大；该非线性增长可由缺失可见度幅度随阈值升高而增长解释，并与线性化天空校准理论一致。附录以缺失/全模型可见度角功率谱比值R_vis作定量比较，阈值1、10、100 mJy归一化后R_vis为1、2.27、8.40，未平滑增益幅度散射为1、2.33、8.39、相位散射为1、2.30、8.55，吻合约3%。
3.5 位置误差下的延迟功率谱（Delay power spectrum from positional errors）：通过在k_⊥=0.01 h Mpc^?1代表模式上比较源减除前后延迟功率谱，研究人员发现位置偏移增强k_∥≈0峰并引入向高延迟模扩展的过剩功率；50′′偏移在k_∥～0.2–0.3的残留功率超过10⁸。源减除后，无误差目录几乎消除残留，5′′偏移在k_∥≈0仍约10⁵，50′′超过10⁹；移动平均不显著改变整体谱形，但可降低50′′情形EoR窗口内功率。结论是目录天体测量精度对抑制EoR分析系统污染至关重要。
3.6 流量密度误差下的延迟功率谱（Delay power spectrum from flux-density errors）：通过将流量扰动从0.01%变化到10%并比较源减除前后谱，研究人员发现?0.01%误差与无误差情形几乎相同，0.1%及以上使k_∥≈0功率增加并向高延迟模扩展；源减除后百分比级误差产生明显残留，而N=10的LST方向移动平均在EoR窗口大致降低一个量级。结论是亚百分流量不确定度影响较小，百分比级误差导致显著前景泄漏，移动平均不能完全补偿大流量误差。
3.7 目录不完备下的延迟功率谱（Delay power spectrum from catalog incompleteness）：通过应用1 μJy、1 mJy、10 mJy、100 mJy阈值并比较源减除前后谱，研究人员发现目录越浅，k_∥≈0峰越强并扩展至高延迟模；100 mJy最浅情形峰达10¹²–10¹³。源减除后，1 μJy完全目录几乎消除EoR窗口功率，1 mJy阈值在k_∥≈0.2残留约10⁸，100 mJy约10⁹；移动平均可改善最多约一个量级，但1 mJy仍不足强弱标准21 cm场景。结论是延迟功率谱保真度强烈依赖目录深度，10–100 mJy级不完备性无法仅靠LST方向移动平均弥补。

讨论部分总结：移动平均沿LST方向通常将增益解散射抑制30%–60%，原因是目录不完美引起的校准误差随时间变化，对连续LST区间平均可抵消大部分变化；频率方向平滑如多项式拟合或正则化也可考虑，但过度平滑有去除真实天体物理谱特征甚至21 cm信号的风险。探测所需动态范围由DR=P_fg(avg)/P₂₁给出，模拟采用前景平均功率P_fg(avg)?4.6×10¹⁵ mK² h^?3 Mpc³；在|k_∥|?0.2 h Mpc^?1附近考虑奇异、强、弱三种21 cm场景，功率分别取2.5×10⁸、5.0×10⁵、2.2×10⁴ mK² h^?3 Mpc³，对应所需DR为1.6×10⁷、9.2×10⁹、2.1×10¹¹。位置误差0.′′5时DR～10¹²，足以满足弱信号约10¹¹要求；5′′时即使强信号也不足，但增益平滑后5′′仍可达强信号要求；50′′时DR～10⁸，仅奇异场景可探测。流量误差～0.01%时DR～10¹³，足以弱信号；～0.1%时DR～10¹¹，不足弱信号但平滑后可；～1%时DR～10⁹，不足强信号；～10%时DR～10⁷，所有场景均不足，但平滑后可探测奇异场景。目录不完备至1 mJy时DR仅～10⁸，平滑后～10⁹，仍不足强场景；10 mJy足以探测EDGES启发的奇异场景，100 mJy接近下边界而不稳健；外推表明探测强信号需100 μJy目录深度并结合增益平滑。与当前目录比较，GLEAM-X覆盖72–231 MHz，170–231 MHz平均灵敏度1.5 mJy beam^?1，亮源位置好于1′′，流量标度不确定度约2%（内部）和2%–5%（外部），72–80 MHz RMS约10 mJy beam^?1，5.8 mJy和10.2 mJy处完备性分别约50%和90%；其位置精度满足要求，但2%–5%流量标度不满足强场景?0.1%要求，且5 mJy处仅约50%完备、低频?80 MHz灵敏度和流量标度可能退化。TGSS在150 MHz约25 mJy处50%完备、位置好于2′′但不覆盖全南天；LoTSS在0.8 mJy beam^?1约90%完备、位置0.′′2但主要覆盖北天，不能用于MWA校准；MWA宇宙黎明测量需要SKA-Low（最大基线75 km）等下一代射电源巡天。

研究结论部分翻译：SHINONOME项目利用MWA Phase II超低频段（73.6–104.3 MHz）漂移扫描识别宇宙黎明频率系统效应、改进校准天空理解并检验强谱涨落模型。本研究首先分析真实MWA观测增益变化：增益幅度和相位的ACF首次过零时标分别约～2–3 h和～2 h，Allan偏差极小值分别约～45 min和1.3 h。随后系统研究校准天空目录不完美对宇宙黎明21 cm功率谱测量的影响，端到端模拟纳入位置偏移、流量密度不确定度和目录不完备性三类误差。三类误差均在天线增益解中诱发频率依赖涨落，并在延迟谱中产生显著残留功率；角秒级位置误差、百分比级流量密度不准确和高于mJy级的不完备性均可产生与预期宇宙黎明21 cm信号相当或更强的污染。稳健探测需相对前景达到极高动态范围（～9.2×10⁹–2.1×10¹¹），意味着校准目录须具备亚角秒天体测量精度、亚百分流量精度和亚mJy级完备性。与现有低频巡天比较，GLEAM-X在带宽覆盖和均匀性方面代表当前先进水平，但流量密度精度仍在百分之几水平，深度也未在mJy级完全完备；现有目录尚不满足宇宙黎明功率谱分析所需严格精度，亟需下一代更深、更准确目录。研究人员采用与Allan偏差所指增益幅度稳定时标一致的～50 min平滑窗口，使增益解统计散射通常降低30%–60%；该方法缓解目录误差引起的随机涨落并部分放宽目录精度要求，但不能补偿不准确天空模型引入的系统谱结构。总体而言，校准目录精度是21 cm功率谱分析系统污染的主导因素之一，现有目录尚未达到必要指标；LST方向移动平均是有用的缓解策略，但必须与改进的天空目录结合。这些发现为即将开展的MWA超低频观测和未来宇宙黎明实验的校准策略提供重要指导。未来工作将聚焦于构建这些频率上位置和流量密度精度更高的天空目录，并确定可靠分析所需的目录完备性阈值。