[論文レビュー] A Detection Threshold in the Amplitude Spectra Calculated from TESS Time-Series Data

A. S. Baran, C. Koen|arXiv (Cornell University)|Jun 17, 2021

Spectroscopy and Chemometric Analyses被引用数 26

ひとこと要約

論文は、TESSの時系列データの振幅スペクトルの偽陽性確率(FAP)を0.1%としたときのS/N検出閾値を推定し、データカバレッジとカデンスの関数として経験的S/Nを導出し、他のFAPやデータ長の閾値を調整する式を提供する。

ABSTRACT

We present results of time-series data simulation. We aimed at estimating the threshold used for detecting signals in amplitude spectra, calculated from simulating TESS photometry of up to one year duration. We selected the threshold at a false alarm probability FAP=0.1% and derived S/N ratios between 4.6 and 5.7 depending on the data cadence and coverage. We also provide a formula to estimate the threshold for any FAP adopted and a given number of data points. Our result confirms that, to avoid spurious detection, space-based photometry may require substantially higher S/N than that typically being employed for ground-based data.

研究の動機と目的

• 宇宙ベースの時系列解析において、地上データよりも偽陽性ピークの発生率が高いため、頑健な検出閾値の必要性を動機づける。
• TESSの超短、短、長カデンスで、1–13観測セクターに跨るFAP = 0.1%に対応するS/N閾値を推定する。
• 任意のFAPとデータ点数に対して閾値を推定する実用的な式を提供する。
• 宇宙ベースの光度測定は、偽検出を回避するために地上データより高いS/Nを必要とすることを示す。

提案手法

• 超短(20 s)、短(120 s)、長(1800 s)カデンスにおけるGaussianノイズのTESS時系列データを1–13セクターにわたってシミュレーションする。
• 高速フーリエ変換を用いて振幅スペクトルを計算し、スペクトルを中央値ベースの標準化で標準化する。
• 指定された偽陽性確率(FAP)を用いてスペクトル最大値として検出閾値をモデル化し、セクター数 Ns の関数としてS/N比を導出する。
• 標準化 W および U を用いたスペクトル最大値の累積分布関数の解析表現を導出し、シミュレーション結果と比較する。
• Ns（Eq. 3）に対するS/Nの閉形式の関係を提供し、周波数範囲を制限した場合の有効データ点数N_effの調整について議論する。
• FAP = 0.1% に跨るカデンスと Ns に対するS/N閾値を表として示す。

実験結果

リサーチクエスチョン

• RQ1 TESSの時系列データについて、FAP = 0.1% の検出を主張するためには振幅スペクトルでどのS/N閾値を超える必要があるか？
• RQ2 カデンス（超短、短、長）と観測セクター数の増加によって必要なS/N閾値はどのように変化するか？
• RQ3 異なるFAPや関心の周波数範囲を制限する場合、閾値をどう調整できるか？
• RQ4 データギャップやカデンス観測は偽陽性閾値にどの程度影響するか？
• RQ5 宇宙ベースの光度測定で地上データと比較して偽検出を避けるための実用的な指針は何か？

主な発見

Number of sectors	S/N (ultra-short 20 s)	S/N (short 120 s)	S/N (long 1800 s)
1	5.292	5.037	4.615
2	5.392	5.124	4.731
3	5.433	5.194	4.790
4	5.468	5.238	4.835
5	5.497	5.266	4.877
6	5.538	5.287	4.890
7	5.555	5.321	4.924
8	5.575	5.330	4.947
9	5.585	5.342	4.962
10	5.604	5.362	4.981
11	5.622	5.369	4.998
12	5.628	5.391	5.002
13	5.629	5.397	5.013

• FAP = 0.1% の場合、S/N閾値はセクター数 Ns によって超短で約5.3–5.6、短で約5.0–5.4、長カデンスで約4.6–5.0 の範囲になる（Ns に依存）。
• S/N閾値はセクター数が増え、カデンスが高いデータほど大きくなる。これは宇宙ベースの光度測定が地上の閾値（S/N = 4程度）より高いS/Nを必要とすることを示している。
• S/N閾値は Ns に対する対数フィット（式3）で予測でき、導出されたCDFs（付録）を用いて異なるFAPへ調整できる。
• Nyquist区間の分の周波数範囲に制限して周波数範囲を狭めると、有効データ点数N_eff = N_p fとなり、特定のFAPに対するS/N閾値は低下する。
• セクターのデータにギャップがあると、同じ点数でギャップなしの場合よりS/N閾値が最大約0.25上昇する可能性があり、長いギャップでは別名（エイリアシング）が生じ得る。
• 論文は3つのカデンスに対するNsのS/N閾値の表（Table 1）を提供している。