[论文解读] Searching for periods in X-ray observations using Kuiper's test. Application to the ROSAT PSPC archive

S. Paltani|ArXiv.org|Mar 8, 2004

Astrophysics and Cosmic Phenomena参考文献 25被引用 23

一句话总结

本文提出库珀检验（Kuiper’s test）作为一种稳健方法，用于利用单个光子到达时间检测X射线和伽马射线观测中的周期性信号，有效处理不连续观测并避免间隙带来的污染。应用于ROSAT PSPC档案数据时，识别出28个新的周期源候选体，其在长周期搜索中优于传统方法，尤其在间隙严重影响其他检验性能时表现更优。

ABSTRACT

We use Kuiper's test to detect periodicities in X-ray and gamma-ray observations. Like Rayleigh's test, it uses the individual photon arrival times, and is therefore well suited to the analysis of faint sources. Our method makes it possible to take into account the discontinuities in the observation, and to completely get rid of the contamination that results from them. This makes it particularly adapted to the search of periods long compared to the total observation duration. We propose a semi-analytical approach to determine the effective number of trial frequencies when searching for unknown periods over a frequency range. This approach can be easily adapted to other tests. We show that, using Kuiper's test, we can recover periods in frequency domains where other tests are completely confused by contamination. We finally search the entire ROSAT Position-Sensitive Proportional Counter (PSPC) archive for long periods, and find 28 new periodic-source candidates.

研究动机与目标

• 开发一种用于检测X射线和伽马射线观测中稀疏、时间标记光子数据周期信号的统计方法。
• 解决标准周期检测方法（如傅里叶变换和折叠法）在不连续观测（间隙）下性能下降的挑战。
• 为库珀统计量在均匀相位分布零假设下的显著性提供可靠检验。
• 估计周期搜索中有效试验频率数，考虑观测间隙和非均匀覆盖的影响。
• 将该方法应用于ROSAT PSPC档案全集，以识别新的周期X射线源，特别是长周期源。

提出的方法

• 使用库珀检验，一种基于经验累积分布函数的非参数统计检验，评估折叠于候选周期下的光子相位均匀性。
• 通过使用测试周期 $ P_0 = 1/f_0 $ 将光子到达时间相位折叠，将时间转换为分数相位 $ \psi_i(f_0) \in [0,1) $，并应用于单个光子到达时间。
• 采用渐近公式（公式3）计算库珀统计量 $ V^U $ 的显著性，对小样本量进行校正，以避免低估假阳性概率。
• 引入校正因子 $ R $ 以估计有效试验频率数，考虑观测间隙和非均匀相位覆盖，该因子通过蒙特卡洛模拟推导得出。
• 使用半解析方法计算 $ R $ 作为光子数 $ N $、试验频率数 $ N_F $ 和有效观测时间与总时长比值的函数，发现 $ R $ 主要依赖于 $ N $。
• 通过 $ 10^9 $ 次蒙特卡洛模拟验证显著性检验和校正因子，覆盖多种样本大小和概率阈值，确认其在不同随机数生成器下的稳健性。

实验结果

研究问题

• RQ1 库珀检验能否在具有稀疏时间标记光子数据和显著间隙的X射线观测中可靠检测周期性信号？
• RQ2 在间隙污染观测条件下，库珀检验与雷利检验、傅里叶方法及折叠法相比，在灵敏度和假阳性控制方面表现如何？
• RQ3 在具有非均匀覆盖和中断的观测中搜索未知周期时，有效试验频率数是多少？
• RQ4 对于在暗弱X射线源中常见的小样本量，如何准确估计库珀检验的显著性水平？
• RQ5 利用该方法在ROSAT PSPC档案中可发现哪些新的周期X射线源，特别是长周期源？

主要发现

• 库珀检验在其他方法（如傅里叶变换和折叠法）因间隙严重污染的频率范围内成功检测到周期性信号。
• 渐近显著性公式（公式3）能可靠估计假阳性概率（FPP），FPP的高估是可接受的，且不会导致假阴性。
• 对于样本量 $ N \geq 100 $，渐近公式估计的FPP值与真实值相差不超过40%；对于 $ N < 100 $，高估程度更高但仍具保守性。
• 校正因子 $ R $（用于调整有效试验频率数）主要依赖于光子数 $ N $，当 $ k = 20 $ 时平均值为0.335，且对 $ N_F $ 和有效时间比值的依赖性较弱。
• 该方法已应用于ROSAT PSPC档案全集，成功识别出28个新的周期源候选体，尤其在长周期信号检测中表现优异。
• 该方法对相位偏移的不变性确保了其鲁棒性，不受参考时间 $ t_0 $ 任意选择的影响，这是其相对于折叠法的关键优势。