[论文解读] On the suspected timing-offset-induced calibration error in the Wilkinson microwave anisotropy probe time-ordered data

Boudewijn F. Roukema|arXiv (Cornell University)|Jul 29, 2010

Geophysics and Gravity Measurements参考文献 14被引用 13

一句话总结

本文研究了WMAP时间有序数据中航天器姿态与通量密度时间戳之间疑似-25.6 ms的定时偏移，若在校准过程中未修正此偏移，将导致人工四极矩样信号及像素方差增加。通过在HEALPix分辨率N_side=8下对像素方差进行最小化分析，研究发现最优定时偏移为：K/Ka波段-38±8 ms，Q波段-27±3 ms，V波段-43±8 ms，W波段-47±194 ms，该结果在99.999%显著性水平上否定了WMAP合作组所采用的0±1.7 ms偏移，并表明宇宙微波背景四极矩被高估了64–94%。

ABSTRACT

In the time-ordered data (TOD) files of the WMAP CMB observations, there is an undocumented timing offset of -25.6 ms between the spacecraft attitude and radio flux density timestamps. If the offset induced an error during calibration of the raw TOD, then estimates of the WMAP CMB quadrupole might be substantially in error. A timing error during calibration would not only induce an artificial quadrupole-like signal in the mean sky map, it would also add variance per pixel. This variance would be present in the calibrated TOD. Low-resolution map-making as a function of timing offset should show a minimum variance for the correct timing offset. Three years of the calibrated, filtered WMAP 3-year TOD are compiled into sky maps at HEALPix resolution N_side=8, individually for each of the K, Ka, Q, V and W band differencing assemblies (DA's), as a function of timing offset. The median per map of the temperature fluctuation variance per pixel is calculated and minimised against timing offset. Minima are clearly present. The timing offsets that minimise the median variance are -38 \pm 8 ms (K, Ka), -27 \pm 3 ms (Q), -43 \pm 8 ms (V), and -47 \pm 194 ms (W), i.e. an average of -30 \pm 3 ms, where the WMAP collaboration's preferred offset is 0 \pm 1.7 ms. A non-parametric bootstrap analysis rejects the latter at a significance of 99.999%. The hypothesis of a -25.6 ms offset, suggested by Liu, Xiong & Li from the TOD file timing offset, is consistent with these minima. It is difficult to avoid the conclusion that the WMAP calibrated TOD and inferred maps are wrongly calibrated. CMB quadrupole estimates (3/pi)C_2 based on the incorrectly calibrated TOD are overestimated by roughly 64 \pm 6% (KQ85 mask) to 94 \pm 10% (KQ75 mask). Ideally, the WMAP map-making pipelines should be redone starting from the uncalibrated TOD and using the -25.6 ms timing offset correction.

研究动机与目标

• 调查WMAP TOD文件中航天器姿态与通量密度时间戳之间未记录的-25.6 ms定时偏移是否引发校准误差。
• 确定该偏移是否可能在校准后的天图中产生人工四极矩样信号并增加像素方差。
• 检验通过在不同定时偏移下对像素方差进行最小化处理，能否识别出真实的校准偏移。
• 评估此误差对CMB四极矩振幅估计的影响，特别是针对标准天区掩膜（KQ85与KQ75）。

提出的方法

• 对WMAP三年校准后滤波的时间有序数据（TOD）在HEALPix分辨率N_side=8下进行重映射，针对K、Ka、Q、V和W波段的差分组装（DAs），将定时偏移从-100 ms至+100 ms进行扫描。
• 对每个偏移值，生成温度涨落δT及其平方值(δT)²的天图，以计算像素方差。
• 计算每幅图中像素方差的中位数，并将其作为定时偏移的函数进行最小化，以识别最优校准偏移。
• 采用非参数自助法（bootstrap）分析，评估方差最小值相对于WMAP合作组首选偏移0±1.7 ms的统计显著性。
• 通过使用高度迭代的天图和替代拟合技术，检验方法的稳健性，并评估与四极矩或前景污染相关的潜在偏差。
• 在所有五个频段之间比较结果，并评估对观测到的方差最小值的其他可能解释。

实验结果

研究问题

• RQ1 WMAP TOD文件中航天器姿态与通量密度时间戳之间存在-25.6 ms的定时偏移，是否会导致校准误差并增加天图中的像素方差？
• RQ2 在多个频段中，何种定时偏移可使校准后WMAP天图的像素方差最小化？
• RQ3 WMAP合作组所采用的首选定时偏移0±1.7 ms是否与方差最小化偏移在统计上一致？
• RQ4 由于此校准误差，CMB四极矩振幅在KQ85与KQ75天区掩膜下被高估了多少？
• RQ5 观测到的方差最小值是否可能为统计巧合或与四极矩信号相关的偏差，而非真实的校准错误？

主要发现

• 使中位数像素方差最小化的定时偏移为：K与Ka波段-38±8 ms，Q波段-27±3 ms，V波段-43±8 ms，W波段-47±194 ms。
• 这些最优偏移的平均值为-30±3 ms，与WMAP合作组所采用的0±1.7 ms偏移在99.999%显著性水平上不一致。
• Liu等人（2010a）通过文件时间分析提出的-25.6 ms偏移，与所有波段中方差最小化偏移一致。
• 由于校准错误，基于错误校准TOD的CMB四极矩估计值被高估了64±6%（KQ85掩膜）至94±10%（KQ75掩膜）。
• 方差最小化方法在不同拟合技术、地图迭代水平和频段下均表现稳健，且对结果的其他解释被发现为推测性或与数据不一致。
• 本研究结论为：校准步骤可能使用了错误的时间戳，引入了系统性误差，且未在地图生成过程中修正，导致CMB四极矩被显著高估。