The impact of interpolation in high-resolution spectroscopy -- The overlooked role of interpolation in radial velocity extraction

技術要約：高解像度分光法における補間（Interpolation）の影響

問題提起
地球型系外惑星の検出と特性評価には、cm s$^{-1}$レベルの視線速度（RV）精度が求められる。ESPRESSOやEXPRESといった最先端の分光器は高い計器安定性を実現しているが、データ解析パイプライン自体に起因する偽の信号が発生することがある。具体的には、テンプレートに基づくRV抽出手法では、高SNRのテンプレートを構築し、個々の観測データをそれと比較するために、ステラスペクトルを共通の波長グリッド上に補間する必要がある。本論文は、このプロセスにおける補間アルゴリズムの選択によって導入される系統的なバイアスについて調査している。著者らは、補間はスペクトル線の特定の形状を暗黙のうちに仮定していると主張している。もしこれらの仮定が成立しない場合、このプロセスはフラックスの残差、ラインの非対称性、および系統的なRVバイアスを引き起こす可能性がある。これは、バリセン地球半径速度（BERV）の変化がピクセルサイズに対して小さい場合に特に顕著となる。

手法
本研究は、合成データの生成と実際のESPRESSO観測の分析を組み合わせた多角的なアプローチを採用している。

1. 合成データの生成:

   • ガウスプロファイル: 著者らは、補間の影響を分離するためにガウス関数を用いて合成データセットを生成した。単一のガウス線と、「フォレスト（林）」状の線（$\epsilon$ IndiおよびHD40307でARESによって検出された実線のパラメータから導出されたパラメータを合算したもの）をシミュレートした。
   • サンプリングの変動: BERVの運動と計器のドリフトの効果を模倣するため、合成スペクトルを0から500 m s$^{-1}$（ESPRESSOの平均ピクセルサイズをカバー）の範囲でドップラーシフトさせた。その後、元の波長グリッドへとフラックスを補間した。
   • ノイズの注入: 決定論的な系統的効果と確率的な系統的効果の間の破壊的な干渉を評価するため、ノイズのない条件下と、ポアソンノイズを注入した条件（SNR=10からSNR=78の範囲）の両方でシミュレーションを行った。

2. 補間アルゴリズム:
   本研究では、s-BART（Template matchingを用いたRVのための準ベイズ的アプローチ）パイプライン内で実装されている5つの補間戦略を比較した：

   • 3次スプライン（デフォルト）。
   • 2次スプライン。
   • 部分区分3次エルミート補間多項式（PCHIP）。
   • Akima。
   • Floater-Hormann。

3. 実観測:
   異なる観測戦略を持つ4つの恒星のESPRESSO観測に対してs-BARTパイプラインを適用した：

   • 高頻度、単一夜観測: アステロセイズモロジー・キャンペーンに使用されているHD40307 (K2) および $\epsilon$ Indi (K6)。
   • 複数夜／長期観測: 系外惑星の検出と特性評価に使用されているProxima Centauri (M5) および HD39091 (G0)。
     RVは異なる補間アルゴリズムを用いて抽出され、デフォルトの3次スプラインの結果と比較された。

主な結果

• フラックス残差と非対称性: 補間は、特にスペクトル線の中心がピクセルノードの間に位置する場合、スペクトル線のフラックスに残差を導入する。これらの残差の大きさは、サンプリング密度（FWHMあたりのピクセル数）および使用されるアルゴリズムに依存する。3次スプラインおよびFloater-Hormann戦略は、単一ガウス線において一般に100 ppm未満の残差を示した。しかし、ラインプロファイルの非対称性が観察されており、これはサンプリングが不十分なライン（低SNRまたは広いピクセル間隔）においてより顕著となる。
• シミュレーションによるRVバイアス（単一線および線のフォレスト）:
  • ノイズのない、完全に対称なガウスシナリオにおいて、補間は$\sim$mm s$^{-1}$のピーク・トゥ・ピーク振幅を持つ系統的なRV信号を導入した。この信号は、BERVの変動と相関した正弦波状の形状を示した。
  • 低SNRシナリオ（SNR=10）では、この系統的なバイアスの振幅が大幅に増加し、$\sim$20–40 m s$^{-1}$のピーク・トゥ・ピーク値に達した。
  • バイアスの振幅はSNRと逆相関している。SNRが増加するにつれて、系統的な信号は減少し、mm s$^{-1}$の領域へと収束していく。
  • 本研究は、SNRに対するRV応答における3つのレジームを特定した：低SNRおよび高SNRにおける正弦波挙動、そして中間SNRにおけるほぼ平坦なプロファイルである。これは、決定論的な系統的効果と確率的な系統的効果の間の破壊的な干渉によるものと考えられる。
• 実観測（単一夜）:
  • 小さなBERV変化（ピクセルサイズと同程度）を伴う高頻度観測において、補間アルゴリズムの選択はRV時系列に大きな影響を与える。
  • 低SNRのケース（HD40307）では、異なるアルゴリズムが$\sim$15–25 m s$^{-1}$のピーク・トゥ・ピーク残差を生じさせた。対称なシミュレーションで良好な結果を示したFloater-Hormannアルゴリズムは、実際には実データにおいて系統的信号を増大させた。
  • 高SNRのケース（$\epsilon$ Indi）では、アルゴリズム間の差異はより小さくなり、ピーク・トゥ・ピーク残差は約1 m s$^{-1}$であり、1-$\sigma$の不確かさの範囲内に収まった。
• 実観測（複数夜）:
  • 大きなBERVウィンドウにわたる観測（Proxima CentauriおよびHD39091）では、補間によって引き起こされる系統的信号は大幅に減衰する。
  • 補間アルゴリズムの選択によるピーク・トゥ・ピークのRV散らつきは、20 cm s$^{-1}$の上限値であることが判明しており、これは測定の1-$\sigma$不確かさよりも十分に低い値である。

意義と結論
本論文は、スペクトルの補間は、特に以下の特定の科学ケースにおいて、RV抽出における無視できない系統誤差の源であることを結論付けている：

1. 高頻度、低SNRキャンペーン: 短いベースラインの集中的な観測（例：アステロセイズモロジーや高周波モニタリング）を伴い、BERV変化が小さい科学ケースは、m s$^{-1}$レベルに達する補間誘発バイアスに対して最も脆弱である。
2. 幅広い影響: 長期的な高SNRの系外惑星サーベイ（BERVの範囲が大きい場合）ではバイアスは最小化されるものの、著者らは、スペクトルが共通のグリッド上に配置されるすべての科学ケース（系外惑星の検出だけでなく、大気研究、アステロセイズモロジー分析、および宇宙論的な赤方偏移の決定を含む）において、補間アーティファクトの影響を受けると主張している。
3. アルゴリズムへの依存性: 普遍的に最適とされる単一の補間アルゴリズムは存在しない。その性能は、SNR、サンプリング密度、およびスペクトル線の固有の形状に大きく依存する。本研究は、補間手法に内在するライン形状に関する仮定が、フラックスの過小評価およびライン中心の系統的なバイアスを招く可能性があることを強調している。

著者らは、理想的な条件下でのmm s$^{-1}$レベルのバイアスは現在の計器精度を下回っているものの、低SNRかつ高頻度のシナリオで観察されるm s$^{-1}$レベルのバイアスは、将来の高精度RV解析において考慮すべき重要な系統的誤差であることを強調している。