Estimating stellar metallicities from Gaia DR3 XP data using LAMOST DR10

技術要約：LAMOST DR10を用いたGaia DR3 XPデータからの恒星金属度の推定

問題提起
Gaiaミッションは10億個以上の恒星に対して天体物理学的パラメータを提供しているが、Gaia DR3のGSP-Photモジュールから導出された金属度（[M/H]）は、色、等級、および消光に依存する顕著な系統的バイアスを抱えている。これまでの取り組みとして、Gaiaの低解像度ブルー・レッドフォトメーター（BP/RPまたは「XP」）スペクトルをAPOGEEなどの分光サーベイにアンカー（基準付け）して金属度を再導出する試みがなされてきたが、特に金属欠乏側の裾野や塵の多い円盤領域をより良く制約するために、LAMOSTの広範な光学領域を活用した、AFGK型星のための均質かつ高精度な金属度カタログへのニーズが依然として存在する。

手法
著者らは、LAMOST DR10の高精度な分光ラベルを用いて学習させた、Gaia DR3の入力のみから恒星の鉄存在量（[Fe/H]）を推定する機械学習アプローチを提示する。

• 学習データ： 本研究では、精査されたLAMOST DR10の低解像度AFGKカタログのサブセットを利用している。高S/N比および十分に制約された大気パラメータを確保するために品質カットを適用し、結果として約120万個の星からなる学習セットを作成した。ラベルは、範囲 $-2.5 \le [Fe/H] \le 1.0$ に制限された、繰り返し測定の逆分散重み付き平均である。
• 入力特徴量： モデルは、Gaia DR3からのみ導出された139個の特徴量を入力とする。これには以下が含まれる：
  • 110個のXPスペクトル係数（BP 55個、RP 55個）。これらは、スペクトルの形状を明るさから分離するため、標準的なフラックスレベル（$G=15$）に再スケーリングされている。
  • 広帯域測光および色（$G, G_{BP}, G_{RP}$）。
  • アストロメトリおよび品質指標（年周視差、固有運動、RUWE、フラックス過剰）。
  • 光度および運動学的プロキシ（疑似絶対等級、減少固有運動）。
  • 消光項（GSP-Photによる塵の推定値 $A_G$ および $E(BP-RP)$）。これは、事前分布を注入することを避けるため、GSP-Photの大気パラメータを明示的に除外している。
  • 銀緯（$b$）としての弱い位置指標。
• モデルアーキテクチャ： 勾配ブースティング決定木回帰器（XGBoost）を採用した。学習には、検証セットの性能に基づく早期停止（early stopping）を伴う二乗誤差目的関数を用いた。
• バイアス緩和： 限界等級に制限されたフィールド星のサンプルに特有の「色－金属度」の事前分布（これにより、化学的に均質な散開星団において偽の金属度傾向が生じる）に対処するため、著者らはサンプル重み付けスキームを実装した。これは、カラー・マグニチュード図（CMD）のセル内で星を重み付けすることで、金属リッチな集団と金属欠乏な集団の表現をバランスさせ、CMDの位置と金属度の間の学習された相関を弱めるものである。
• 不確実性推定： 星ごとの不確実性は、データのランダムな再サンプリングに基づいて訓練された5つのXGBoostモデルのアンサンブルの標準偏差から導出される。

主な結果
モデルは、保持された（held-out）LAMOSTデータ、外部サーベイ（GALAH DR4およびAPOGEE DR17）、および46個の散開星団を用いて検証された。

1. 保持されたデータに対する性能： LAMOSTテストセットにおいて、モデルは平均絶対誤差（MAE）0.052 dexおよび$R^2$ 0.94を達成し、バイアスは無視できる程度であった。これは、同一の星に対するネイティブなGaia DR3 GSP-Photの金属度（MAE = 0.242 dex）と比較して、約5倍の改善を示している。
2. 外部検証： 推定値は独立したサーベイによく適合する：
   • GALAH DR4： MAE = 0.066 dex。
   • APOGEE DR17： MAE = 0.068 dex。
   • 残差は有効温度、表面重力、または色に対して強い傾向を示さないが、GALAHとの比較において、極端に金属欠乏な星（$[Fe/H] < -1.0$）で平均値へと引き寄せられる「シュリンケージ（収縮）」現象が見られる。これは、この領域の予測値が下限値として扱われるべきであることを示唆している。
3. 散開星団： 高分解能の参照金属度を持つ26個の散開星団について、モデルの推定値と分光値の中央値の差は0.041 dexであった。これは、GSP-Photの推定値（0.248 dex）および以前のAPOGEE訓練済みXGBoostモデル（0.067 dex）の両方を上回る性能である。
4. 銀河構造：
   • 垂直勾配： モデルは $-0.237$ dex kpc$^{-1}$ の垂直金属度勾配を回収しており、これは従来のG型矮星の研究と一致している。
   • 半径方向勾配： 運動学的に選択された薄い円盤のサンプルに適用したところ、モデルは折れ曲がった半径方向勾配を明らかにした。内側は $+0.119 \pm 0.018$ dex kpc$^{-1}$、外側は $-0.058 \pm 0.002$ dex kpc$^{-1}$ の傾きを持ち、5.9 kpc付近で折れ曲がっている。これは近年の高分解能分光による知見と一致する。
   • 星団勾配： 散開星団の半径方向勾配は $-0.066$ dex kpc$^{-1}$ と測定され、APOGEEベースの測定値と一致している。

限界と系統誤差
著者らは、散開星団内の色に応じた予測金属度における系統的な「ハンプ（隆起）」を特定した（$G_{BP}-G_{RP} \approx 1.0$ 付近でピーク）。化学的に均質な星団においても存在するこの特徴（振幅の中央値は約0.07 dex）は、温度と金属度が縮退している低解像度（$R \sim 50$）の分光測光の固有の限界に起因するとされている。重み付けスキームはこの影響を軽減したが、完全に排除するには至らなかった。さらに、散開星団内の1 kpcを超える星の予測については、金属度の系統的な過小評価が見られる。

意義
本論文は、Gaia DR3のXPデータが、広範なLAMOST DR10の分光ラベルにXGBoostを介してアンカーされることで、AFGK型星に対して0.05から0.07 dexの精度で恒星金属度を生成できることを示している。この精度は、銀河円盤の化学的研究（複雑な半径方向勾配の回収や散開星団の特性評価を含む）を行うのに十分である。著者らは公開カタログと訓練済みモデルを提供しており、これはネイティブなGaia GSP-Phot製品よりも大幅な改善であり、Gaia DR4のリリースまでの間、他のXPベースの金属度推定に対する競争力のある代替手段となる。本研究は、モデルはフィールド星や広範な集団研究には堅牢であるが、ユーザーは化学的に均質な集団（狭い色の範囲を持つもの）を分析する際に、残留する色依存の系統誤差を考慮する必要があることを強調している。