A constrained linear model for continuum normalization of stellar spectra

問題提起
連続体正規化（Continuum normalization）は、恒星パラメータや化学組成をスペクトル線の強度から推定するための極めて重要な前処理ステップである。しかし、現在の慣習は主観的であったり、信号対雑音比（S/N）によって一貫性が欠けていたり、あるいは連続体正規化を恒星パラメータ推定とは独立したステップとして扱うことでバイアスが生じたりしている。反復的なシグマ・クリッピングや「連続体ピクセル」の手動選択といった従来の手法は、金属に富む星、低S/Nのデータ、および連続体が存在しないスペクトル領域（例：低温の星の青色波長域）において苦戦する。連続体と恒星パラメータの同時モデリングは理論的には理想的であるが、計算コストが高くなったり、非凸最適化を伴ったりすることが多い。既存のデータ駆動型モデルは、大規模な訓練セットを必要とし、不安定になる場合がある。

手法
著者らは、恒星の吸収、テラリック（地球大気）透過、および連続体－インストゥルメント応答を同時に適合させる制約付き線形モデルを提案している。この手法の核となるのは、理論的な連続体正規化済みスペクトルのグリッドに対して適用される非負行列因子分解（NMF）である。

1. モデルの定式化: 観測されたフラックスは、恒星成分（$f$）、テラリック成分（$g$）、および連続体－インストゥルメント応答（$h$）の3つの要素の積としてモデル化される。本モデルは、この乗法関係を線形和に変換するために、対数変換された空間で動作する：$Y = \log f + \log g + \log h$。
2. 恒星およびテラリック吸収: 恒星成分とテラリック成分の両方は、理論スペクトルのグリッドを2つの非負行列（基底ベクトルと重み）に因子分解することによってモデル化される。
   • 星のグリッド（BOSZ）は、スパース性と非負性を確保するために負の対数を用いて変換され、$K$個の基底ベクトルに因子分解される。
   • テラリック・グリッド（PypeItプロジェクトによるもの）も同様に、$L$個の基底ベクトルに因子分解される。
   • 決定的なことに、非負制約により、再構成された正規化フラックスは厳密に1以下（$f \le 1$）となることが保証される。これにより、吸収特徴と連続体が物理的に分離される。すなわち、恒星モデルの予測値を超えるフラックスは、すべて連続体成分に帰属される。
3. 連続体モデリング: 連続体－インストゥルメント応答は、フーリエ基底（正弦波および余弦関数）を用いてモデル化され、滑らかな変化に対する柔軟かつ線形な表現を提供する。
4. 推論: 推論問題は、恒星およびテラリックの重みに対する非負制約を伴う重み付き最小二乗最適化へと簡略化される（$\alpha \ge 0, \beta \ge 0$）。モデルは線形であり、制約は凸であるため、解は安定しており、高速で、恒星パラメータの初期値も必要としない。

主な貢献

• 制約付き線形フレームワーク: 本論文は、連続体正規化を、摂動パラメータ（テラリック）のフィッティングと単一の凸最適化問題の中で統合する手法を紹介している。
• NMFによる物理的制約: 基底ベクトルと重みに非負性を強制することで、モデルは恒星成分が吸収のみを表すことを保証する。これにより、連続体モデルがスペクトル特徴を「吸収」してしまうことを防ぎ、正規化されたフラックスが物理的に境界内（$\le 1$）に留まるようにする。
• 効率性とスケーラビリティ: 本手法は、推論時に高価なスペクトル合成コードや非線形エミュレータを回避する。訓練は理論グリッドに対して一度行われるだけであり、高分解能スペクトル（例：HARPS）に対する推論には、わずか数コア秒しか要さない。
• 拡張性: モデルの線形的性質により、最適化の凸性を崩すことなく、視線速度のシフトや回転広がりといった他の線形な摂動を組み込むことが可能である。

結果
著者らは、ESO/HARPSの高分解能エシェル分光スペクトルを用いて本手法を検証した：

• 実験1（スペクトル型）: モデルは、B型からM型までの星を正常にフィッティングした。これには、訓練グリッドの温度範囲を超えたB4型星（HD 45418）も含まれる。連続体の推定値は、進化段階（赤色巨星を含む）に関わらず一貫していた。
• 実験 2 (S/Nの一貫性): 幅広いS/N比（10から200）にわたる$\alpha$-Centauri Aの繰り返し観測を用い、著者らは連続体正規化における高い一貫性を実証した。
  • S/N > 100 の場合、正規化フラックスの95パーセンタイルにおける標準偏差は 0.22% であった。
  • S/N > 30 の場合、散らつきは 0.46% であった。
  • これらの結果は、本手法が、従来の手法では失敗したり手動調整を必要としたりするような低S/N比においても、一貫した擬似連続体を生成することを示している。

意義と主張
本論文は、この制約付き線形モデルが、主観的な連続体正規化に代わる、堅牢で再現可能かつ物理的根拠に基づいた選択肢を提供すると主張している。非負制約を通じて連続体を吸収特徴から区別することで、本手法は、連続体ピクセルが存在しないスペクトル領域においても、「真の」連続体推定を提供することを保証する。著者らは、本モデルは連続体正規化のために設計されており、必ずしも精密な恒星パラメータ推論のためのものではない（基礎となる物理学の欠陥を隠蔽する可能性があるため）ことを強調しているが、摂動を同時に、かつ様々なS/N比において一貫してフィッティングできる能力は、大規模な分光サーベイにとって優れたツールとなる。本手法は、スタックされたスペクトルに適用された従来の最先端手法（例：RASSINE）よりも、1桁低いS/N比においてサブパーセントの一貫性を達成している。