"Perfect" spectra for ESO's HARPS spectrograph

技術要約：ESOのHARPS分光器における「完璧な」スペクトル

問題提起
天文学用分光器における標準的なスペクトル抽出は、1980年代に開発された「最適抽出（optimal extraction）」という手法に依存している。この手法は、スペクトルのトレースが検出器のグリッドと完全に一致しており、かつ点広がり関数（PSF）が分散方向（$x$）とクロス分散方向（$y$）で分離可能（独立）であるという仮定に基づいている。実際には、HARPSのような高分解能エシェル分光器において、この仮定が満たされることは稀である。PSFが非分離である場合、標準的な抽出法はフラックス推定、波長較正、および分解能に誤差を導入する。これらの誤差は、基本物理定数の研究や系外惑星の大気観測など、ラインの形状や中心の正確な測定を必要とする精密科学にとって致命的となる。

手法：スペクトル・パーフェクショニズム（SP）
著者らは、高分解能エシェル分光器に対して「スペクトル・パーフェクショニズム（SP）」を初めて適用した事例である「Harperfect」を提示している。SPは、較正行列（$A$）とノイズ共分散行列（$N$）を用いて、2次元検出器データ（$p$）から固有の1次元スペクトル（$\hat{f}$）を再構成するフォワードモデリング・アプローチである。モデルは $p = Af + n$ と定義される。

主な手法構成要素は以下の通りである：

1. 較正行列の構築： 複雑な装置モデルに依存していた従来のSPの実装とは異なり、Harperfectは、レーザー周波数コム（LFC）の観測から行列 $A$ を完全に構築する。LFCは既知の周波数を持つ約10,000本のラインを提供する。HARPSでは個々のLFCラインを分解できないため（ラインが装置の分解能よりも広いため）、観測された各ラインは、特定の検出器位置および波長におけるピクセル積分されたPSFの直接的なサンプルとして機能する。
2. PSFモデリング： 著者らは、21,148本の特定されたLFCラインに対し、7つの自由パラメータ（振幅、中心、標準偏差 $\sigma_x, \sigma_y$、回転角 $\theta$、およびオフセット）を持つ2変数（2D）ガウス関数を用いてフィッティングを行った。これにより、検出器全体にわたるPSF形状の変化のマッピングが可能となった。
3. 波長およびトレース較正： 既知のLFC周波数とフィッティングされたライン中心を用いることで、著者らはガウス過程回帰を行い、2次元波長較正およびエシェル次数の幾何学的中心のトレースを行った。
4. 分解能行列とデコンボリューション： フォワードモデルの直接逆転に伴う「共分散のリンギング（反相関ノイズ）」を回避するため、著者らは分解能行列 $R$ を適用する。彼らは、$C^{-1}$（逆共分散行列）の対称行列平方根から $R$ を導出し、フラックス保存を保証し、単純な自己相関に伴う $\sqrt{2}$ の広がりペナルティなしに、ネイティブな装置分解能を保持するようにしている。
5. 計算実装： 計算の複雑さを管理するため、パイプラインは全フレームを同時に処理するのではなく、検出器の局所的な「リボン」（高さ約60ピクセル）を処理する。線形逆転には疎行列のチョレスキー分解を利用している。

主な貢献

• LFCベースの較正： 本論文は、従来の弧状放電ランプやフラットフィールドに頼ることなく、LFC観測を用いて完全な較正行列（PSF形状、波長、および次数トレース）を構築できることを示した。
• PSFマッピング： 著者らは、前例のない詳細さ（約15物理ピクセルごとにサンプリング）で、HARPS検出器全体のPSF変化をマッピングした。PSFの非分離性を相関係数 $\rho$ を用いて定量化したところ、弱く非分離（$|\rho| < 0.1$）ではあるものの、検出器全体で測定可能な値であることが判明した。
• 次数トレースの精度： LFCラインの中心から導出された次数トレースは、タングステンランプの観測から導出された標準的なHARPSデータ削減ソフトウェア（DRS）のトレースと比較して、次数の視覚的な幾何学的中心により良く一致することが示された。
• ソフトウェア実装： 任意の位置でのモックPSFの補間生成を含む、SPをHARPSデータに適用可能なPythonパイプライン「Harperfect」のリリース。

結果
本手法は、クエーサーHE0515−4414に対する52.5時間のHARPSデータに適用された。

• 信号対雑音比（S/N）： HarperfectのスペクトルのS/Nは、標準的なHARPS DRSの最適抽出よりも一貫して1%から5%低い。著者らは、これがSPフレームワークの失敗ではなく、現在のPSFモデルの限界に起因すると考えている。2変数ガウスモデルは、PSFの非対称性や背景の変化を完全には捉えきれず、それが較正行列 $A$ の不完全な不確かさ推定につながっている。
• ラインプロファイル： 6つの金属吸収遷移の視覚的比較において、HarperfectとDRSのスペクトルはほとんどのラインで密接に一致している。しかし、2つの遷移（Fe II $\lambda$2600 および Mg II $\lambda$2804）については、Harperfectのプロファイルがラインコア付近で視覚的に狭くなっているように見えた。著者らは、定量的なフォークト・プロファイル・フィッティングなしでは、これは予備的な視覚的観察であり、確認された分解能の向上ではないと注意を促している。
• 較正精度： LFCベースの波長較正は、DRSによる独立したThArランプ較正と250 m s$^{-1}$ 以内で一致している。残差には、HARPS検出器の既知のピクセルサイズ異常に対応する、512ピクセルごとの周期的な構造（$\pm 25$ m s$^{-1}$）が見られる。

意義と主張
本論文は、非分離なPSFに対して標準的な抽出法では保証できない特性である、厳密な分解能行列と独立したサンプルを提供することで、Harperfectが「完璧な」スペクトル抽出に向けた重要な一歩であると主張している。

著者らは、即時的な性能向上に関する主張については控えめである。彼らは、現在のS/Nのペナルティが簡略化されたPSFモデル（2変数ガウス）の結果であることを明示しており、完全なスペクトル・パーフェクショニズムの恩恵は、より正確なPSFモデリング、優れた背景処理、およびピクセルサイズ異常の補償が実現して初めて得られるものであるとしている。本研究の主要な価値は、概念実証（プルーフ・オブ・コンセプト）として示されている。すなわち、LFCを用いて必要な較正行列を構築できることの実現可能性を検証し、基本定数や系外惑星大気の測定を向上させるための枠組みを確立した点にある。著者らは、現在の処理時間（フルスペクトルで約2時間）が標準的な抽出よりも遅いものの、精密科学のために計算コストを支払う価値があるほど、スペクトルの忠実度向上のポテンシャルは高いと述べている。