[論文レビュー] 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

ひとこと要約

本稿では、個々の光子到着時刻を用いた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線およびガンマ線観測におけるスパarsely、タイムタグ付き光子データを対象とした周期信号検出のための統計的手法を開発すること。
• フーリエ変換やエポックフォールディングといった標準的手法の性能を著しく低下させる、観測の不連続性（ギャップ）という課題に対処すること。
• 一様な位相分布の帰無仮説の下でのクイパー統計量の信頼性の高い有意水準のテストを提供すること。
• 観測ギャップや非一様なカバレッジを考慮した周期探索における、有効な試行周波数の数の推定を行うこと。
• 全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 クイパー検定は、スパarsely、タイムタグ付き光子データに加え、顕著なギャップを含む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 $ の選択にかかわらず、堅牢である。これはエポックフォールディング法に比べて顕著な利点である。