The Effects of Instrumental Deadtime on NICER Timing Products

技術要約：NICERのタイミング製品における計器デッドタイムの影響

問題提起
中性子星内部組成探査機（NICER）のX線タイミング計器（XTI）は、前身のRXTEやXMM-Newtonを凌駕する、前例のないタイミング分解能（~300 ns）と有効面積を提供している。しかし、X線連星における高いカウントレートは、計器によるデッドタイム（検出器が新しいイベントを処理できなくなる期間）を導入する。RXTEやNuSTARといった過去のミッションでは、デッドタイムに起因するポアソンノイズの歪みを考慮するために、特定の周波数依存の補正をパワースペクトルに適用する必要があったが、確立されたこれらの手法がNICERのマルチ検出器アーキテクチャにおいて必要かつ適切であるかどうかは依然として不明である。具体的には、著者らは、NICERのデータにおいて、フーリエ解析の前にライトカーブ（光度曲線）をデッドタイムで補正する必要があるのか、あるいは、補正されていないライトカーブから周波数依存のノイズを減算するという標準的なアプローチが依然として十分であるのかを調査している。

手法
著者らは、低、中、高のカウントレート領域（~31 から ~7000 counts/s の範囲）を代表するように選択された、X線連星 GX 33 9-4 の9つの観測データを分析した。研究では、テレメトリ飽和、南大西洋異常（SAA）通過、およびバックグラウンドイベント（アンダーシュート、オーバーシュート、強制トリガー）を除外した、標準的なNICER処理パイプラインを用いてグッドタイムインターバル（GTI）を生成した。

分析は3つの段階で進められた：

1. 検出器レベルのデッドタイム特性評価： 著者らは、各焦点面モジュール（FPM）単位でデッドタイムを調査した。52個の稼働中の検出器にわたり、記録されたイベントあたりの平均デッドタイム、およびX線イベントあたりの平均デッドタイム（スクリーニング後）を算出した。彼らは、デッドタイムが全検出器間で一貫しているという仮説を検証するため、定数をフィッティングし、簡約カイ二乗統計量を評価した。
2. ライトカーブの補正： カスタム補正がライトカーブに適用された。各タイムビンについて、イベントごとの記録に基づき、デッドタイムがビンの境界を跨ぐイベントを考慮しながら、総デッドタイムを合算した。観測カウントレート ($r_{Obs}$) は、以下の式を用いて「真の」レート ($r_{True}$) に補正された：$r_{True} = r_{Obs} / (1 - \tau_{Dead,prop})$。ここで、$\tau_{Dead,prop}$ はデッドタイムによって失われたタイムビンの割合である。
3. フーリエ領域での比較： 補正されていないライトカーブと、デッドタイム補正済みのライトカーブの両方からパワースペクトルが生成された。
   • 補正されていないライトカーブについては、Zhangら(1995)によるvan der Klis (1989)の定式化（論文の式5）の修正を用いた、周波数依存のポアソンノイズレベルを減算した。
   • 補正されたライトカーブについては、平坦なポアソンノイズレベル ($2/r_e$) を減算した。
   • 得られたパワースペクトルとノイズレベルを9つの観測全体で比較し、系統的な偏差を特定した。

主な結果

• 検出器の一貫性： 9つの観測すべてにおいて、記録されたイベントあたりの平均デッドタイムは約22 $\mu$sでほぼ一定であることが判明した。デッドタイム分布が全検出器間で一定であるという統計的フィッティングの結果、ほとんどの観測において簡約カイ二乗値は0.1未満であり、特定の検出器が著しく大きなデッドタイムイベントの影響を受けていることはないことが示された。
• 比例デッドタイム： ライトカーブのビンにおける分数デッドタイムは、観測カウントレートに対して線形に増加する。解析された最高カウントレート（7000 counts/s）においても、比例デッドタイムは0.6%未満（最大観測値 ~0.006）であった。著者らは、検出器のデッドタイムが有意なレベル（例：1%）に達するには、カウントレートが23,000 counts/sに近づく必要があるが、その領域ではテレメトリのデッドタイム（バッファ飽和）によって既に観測が断片化されているであろうと指摘している。
• パワースペクトルの比較： デッドタイム補正されたライトカーブから得られたパワースペクトルと、周波数依存ノイズ減算を用いた未補正のライトカーブから得られたパワースペクトルを比較したところ、差異は軽微であった。
  • 周波数 $f \lesssim 1$ Hz において、スペクトルはほぼ同一であった。
  • 高周波数域では偏差が生じたが、一般に再ビン化されたスペクトルの誤差範囲内に留まっていた。
  • 偏差は負の値であり、デッドタイム補正されたライトカーブからのパワーは、周波数依存補正されたスペクトルよりもわずかに高かったが、その差の大きさは、検討されたカウントレートにおいては統計的に有意ではなかった。
• ノイズモデリング： 本研究では、NICERの周波数依存ノイズ項（式5における余弦および正弦依存項）は、検出器数が多いこと（50個対5個）および大きなイベントに対するデッドタイムが小さいことから、RXTEと比較して著しく小さいことが判明した。その結果、「理想的な」平坦なポアソンノイズ ($2/r_e$) は、従来の計器よりもNICERにとって実際のノイズフロアのより良い近似となっている。

意義と主張
本論文は、現在NICERでアクセス可能なカウントレート（~7000 counts/sまで）において、計器の検出器デッドタイムの影響は、解析前にライトカーブに対して複雑な補正を行うほどの重要性を持っていないと結論付けている。

著者らは以下のように主張している：

1. 標準的な慣行は依然として有効である： RXTEやその他の検出器に適用されているデッドタイムの取り扱い方法（具体的には、補正されていないライトカーブから周波数依存ノイズを減算すること）は、依然としてNICERのデータに対しても適切である。
2. 事前補正の必要性はない： フーリエ製品を作成する前に、デッドタイム補正をライトカーブに対して行う必要はない。未校正のイベントファイルを処理してデッドタイム補正済みライトカーブを作成するための計算コストとデータストレージ要件は、結果として得られるタイミング製品に対して実質的な利益をもたらさない。
3. 検出器の均一性： 本計器の機能的な検出器間に重大な不一致はなく、個別の検出器のデッドタイム重み付けなしにアンサンブルデータを使用することが妥当であることを示している。

本研究は、NICERのデータが標準的なパイプラインを用いて処理され、歴史的な計器から得られたものと同等のフーリエ領域の製品を生成できるという信頼を提供しており、これによりX線連星進化の長期モニタリングが容易になる。また、著者らは、基礎となるスペクトル形状を推論するための機械学習アプローチが登場しつつある一方で、現在の分析は、NICERのタイミング解析における伝統的な手法の継続的な使用を支持していると述べている。