Owl-z: a Bayesian tool to select z \geq 7 quasars

「Owl-z: $z \gtrsim 7$ のクエーサーを選択するためのベイズ統計ツール」技術要約

問題提起
高赤方偏移（$z \gtrsim 7$）のクエーサーの特定は、初期宇宙、超巨大ブラックホールの起源、および再電離の時代を理解する上で極めて重要である。しかし、広視野サーベイにおけるこれらの天体の測光選択は、その暗さと、混入する集団による深刻なコンタミネーション（汚染）のために困難を伴う。具体的には、Lyman-$\alpha$ ブレイクに起因する高-$z$ クエーサーの赤い色は、褐色矮星（スペクトル型 M, L, T）や中間赤方偏移の早期型銀河（$1 \lesssim z \lesssim 2$）によって模倣される。従来のカラーカット法は、特にクエーサーとコンタミナントのカラー空間が大きく重なる低赤方偏移域（$7 < z < 8$）において、これらの集団を区別することに苦慮することが多い。Euclidのような大規模サーベイの到来に伴い、コストのかかる分光フォローアップの要求を最小限に抑えつつ、高-$z$ クエーサー候補を選択するための、効率的かつ最適化された統計ツールの必要性が高まっている。

手法
本論文では、あるソースが以下の3つのモデル化された集団のいずれかに属する事後確率を算出するために設計された、オープンソースのベイズ・コードである Owl-z を紹介する：高-$z$ クエーサー（$7 \le z \le 12$）、中間赤方偏移の早期型銀河（$1 < z < 2$）、または褐色矮星（MLT型星、スペクトル型 M6–T9）。

• ベイズ形式: 本コードは、観測された測光データ $D$ が与えられたときのソースの型 $t$ に対する事後確率 $P(t|D)$ を計算するためにベイズの定理を用いる。これは、各集団の事前面密度と、ソースの型が与えられた場合のデータの統合尤度（エビデンス）によって駆動される。
• 集団モデリング:
  • クエーサー: 分光エネルギー分布（SED）テンプレートの集合（平滑化された合成スペクトルおよびパラメトリックモデル）と、Willott et al. (2010) および Matsuoka et al. (2023) から導出された二重冪関数状の光度関数（LF）を用いてモデル化される。
  • 褐色矮星: 観測されたスペクトルライブラリ（LHS および SpeX プリズム）と WISE 測光を用いてモデル化される。空間分布は、EWS の深さに対する薄い円盤近似に限定された、指数関数的な垂直・水平分布を持つ薄い円盤としてモデル化され、表面密度は銀河緯度に依存する。
  • コンタミナント銀河: 早期型銀河の SED（SSP モデル）と、Schechter 関数でフィッティングされた E-COSMOS カタログ（Euclidised COSMOS2022）から導出された光度関数を用いてモデル化される。
• 尤度計算: データの尤度は、ガウス型の測光誤差および非検出（強制測光または上限値を使用）を考慮し、モデルのフラックスを中心とする多変量密度の混合として計算される。
• 実装: Python で記述され（高速化のために NumPy, SciPy, Astropy, および Cython を利用）、測光データと銀河座標を入力として受け取り、各集団の確率と最良適合パラメータ（赤方偏移、等級、スペクトル型）を出力する。

主な貢献

1. Owl-z の開発: Euclid Wide Survey (EWS) に特化して調整されているが、他のサーベイ（例：LSST, Roman）にも適応可能な、柔軟かつ高速な公開ベイズツール。
2. 包括的なコンタミナント・モデリング: 褐色矮星の集団（銀河緯度の影響および薄い円盤／厚い円盤の近似を含む）と、主要なコンタミナントとして特定された中間赤方偏移の早期型銀河の両方を明示的にモデリングしている。
3. 実データによる検証: 様々なサーベイ（DECaLS, Pan-STARRS, HSC, UKIDSS, WISE）から得られた、$z > 7$ で分光的に確認されている既知のすべてのクエーサーの再特定に成功しており、異なる測光データセットに対するコードの汎用性を実証している。
4. 性能指標: 選択効率を評価するために、完備率（True Positives / (True Positives + False Negatives)）と純度（True Positives / (True Positives + False Positives)）を組み合わせたグローバルな性能指標である F値（F-measure） を導入した。

結果

• 再特定: Owl-z は、$z > 7$ において分光的に確認されているテスト対象の全8個のクエーサーに対して $P_q > 0.99$ の確率を達成し、これらを MLT 矮星のコンタミナントから正しく区別した。
• 完備率と純度（シミュレーション）: シミュレーションされた EWS データ（$7 \le z \le 12$, $18 \le H_E \le 24.5$）を用いた結果：
  • 明るいソース: $H_E \approx 22$ より明るいソースについては、Owl-z は赤方偏移に大きく依存することなく、高い完備率と純度を伴う完全な性能（高い F 値）に達する。
  • 暗いソース: 低い信号対雑音比（SNR）およびコンタミナントとの混同の増加により、より暗いソース（$H_E > 22$）では性能が低下する。
  • 赤方偏移依存性: 後期型矮星との混同のため、$7 < z < 8$ では完備率が低くなる。$z \approx 10$ では、早期型銀河との混同により完備率の低下が見られる。
  • 測光赤方偏移: 本コードは正確な測光赤方偏移推定値（$z_{out}$）を提供し、$H_E < 24$ において真の赤方偏移に対して 0.97 のピアソン相関係数を示し、外れ値の割合も低い（$\Delta z > 0.1(1+z)$ に対して $\sim 2.5\%$）。
• 閾値の最適化: 本研究は、選択確率の閾値（$\zeta$）を事後処理で調整できることを示している。明るいソースで完備率を最大化するには低い閾値（$\zeta = 0.1$）で十分であるが、暗いソースで早期型銀河のコンタミネーションに対して純度を維持するには、より高い閾値（$\zeta \approx 0.9$）が必要となる。
• 堅牢性: コードの性能は EWS のフットプリント全体で比較的均一であり、銀河座標や SNR の違いによる変動は僅かである。シミュレーションにおける厚い円盤の褐色矮星の導入は、EWS の深さにおける F 値への影響は無視できる程度であり、このサーベイにおける薄い円盤近似の妥当性を裏付けている。

意義
本論文は、Owl-z が Euclid やその他の広視野サーベイにおける高-$z$ クエーサーの選択を最適化するための強力な公開ツールであることを主張している。クエーサーをコンタミナントから区別するための堅牢な手法を提供することで、測光および分光による確認キャンペーンの最適化を促進する。等級と赤方偏移に基づいて選択閾値を調整できる能力により、観測者は、より多くの候補を見つけること（完備率）と、その信頼性を確保すること（純度）のトレードオフを調整することができる。著者らは、本ツールは Euclid 用に調整されているものの、そのベイズ的な枠組みと適応性により、次世代の大規模宇宙論サーベイの時代における汎用性の高いソリューションであると強調している。