paired: A Statistical Framework for Detecting Stellar Binarity with Gaia RVs. I. Sensitivity to Unresolved Binaries

技術要約：paired: Gaia RVSを用いた恒星連星検出のための統計的枠組み I. 分解されていない連星への感度

問題提起
Gaia データリリース3（DR3）には、3,300万以上の天体に関する視線速度（RV）の測定値が含まれている。このデータセットは、恒星の多重性とそれが系外惑星の人口統計に与える影響を大規模に研究するためのユニークな機会を提供する一方で、重大な制限が存在する。それは、大多数のソースにおいて、個々のエポックごとのRV測定値が公開されていないことである。その代わりに、ユーザーは中央値のRV、観測数、および中央値の推定誤差（$\epsilon$ または radial_velocity_error）にのみアクセスできる。この誤差項は、観測値の時系列の標準偏差から導出されるものであり、連星性の情報をコード化している可能性があるが、その解釈のための標準化された大規模な枠組みが欠けている。トランジット観測における感度解析は確立されているが、RVノイズを介したGaiaの「分解されていない（unresolved）」連星に対する感度の特性評価は、出現率の研究におけるバイアスを回避するために必要である。さらに、Gaiaは一部の明るい恒星（約600万天体）に対して特定の「ノイズ」キーワード（例：rv_chisqu_pvalue, rv_renormalised_gof）を提供しているが、これらは radial_velocity_error はあるものの公開されたp値を持たない残りの約2,300万天体には利用できない。

手法
著者らは、GaiaのRV観測量をフォワードモデリングし、空間的に分解されていない恒星伴星の尤度を抽出するために設計されたオープンソースの統計的枠組みである paired を導入する。コアとなる手法は、単一の星のRV測定値は、測定不確かさによって支配されるガウス分布に従うはずであるという前提に基づいている。「過剰な」ノイズは、潜在的な連星性を示す。

本フレームワークは以下のステップを通じて動作する：

1. データ駆動型ノイズモデリング： 観測されたGaiaカタログを、色・マグニチュード空間（$BP-RP$ vs. $G$ マグニチュード）のビンに分割する。各ビンについて、「典型的な」観測あたりのRV測定不確かさを推定する。ベースラインとなる「単一星」のノイズモデルを歪ませるノイズの多いソースによる汚染を防ぐため、反復的なクリッピングプロセスを用いてこれらを除外する。
2. p値の計算： 各天体に対し、そのソースが単一の星であるという帰無仮説の下でp値を計算する。これは、システムが真に単一であった場合に、測定されたRV誤差が観測される確率を定量化する。
   • 誤差がサンプル分散に基づく明るいターゲット（$GRVS < 12$）では、p値は閉形式で計算される。
   • より暗いターゲット（$GRVS > 12$）では、誤差が相互相関関数の鋭さから導出されるが、既知のカタログに対して検証された通り、同じ式が経験的に適用される。
3. 注入および回収分析（Injection-and-Recovery Analysis）： 感度を特性化するために、著者らは半長軸（$a$）および質量比（$q$）において一様対数間隔を持つ100万個の合成連星集団を生成した。これらのシステムをGaiaのノイズモデルに注入し、$q$、$a$、および恒星のマグニチュードの関数としての「完全性（completeness）」の検出能を決定した。
4. 半振幅の推論： 有意な過剰ノイズ（$p < 0.01$）を持つソースに対し、観測されたRV分散を利用して、ケプラー運動を仮定したRV半振幅（$K$）を確率的に制約する。

主な貢献

• paired フレームワーク： 物理的な連星特性を予測される「ノイズ」特性へとマッピングするコミュニティツールであり、3,000万のRVSソースの全サンプルに対して、分解されていない連星の特定を可能にする。
• 感度の特性化： 分解されていない連星に対するGaia DR3の完全性の限界を定義する、厳密な注入および回収分析。本研究では、明るい星（$m_G < 7$）において、本手法は数AUまでの褐色矮星、および3 AUにおいて質量比 $q=0.1$ の連星に対して99%の信頼度で感度を持つことを示している。
• カタログ生成： p値と半振幅の制約を含む約3,000万個の恒星のカタログを作成。これには、DR3で rv_chisqu_pvalue が公開されていない約2,300万のソースも含まれ、事実上、より暗い集団に対しても連星診断を拡張している。
• ベンチマーキング： Gaia自身の公開製品（約600万の明るいソース）および独立した地上観測カタログ（Kepler/TESS食連星、SB9分光連星、APOGEE）に対する包括的な検証。

結果

• 検出率： 3,000万のソースのうち、約14.4%が過剰なRVノイズ（$p < 0.001$）を示す。この割合はヘルツスプルング・ラッセル（HR）図全体で変化し、若い主系列星では20%に上昇し、明るい進化星ではさらに高くなる。
• Gaia DR3との比較： 公開された rv_chisqu_pvalue を持つ約600万のソースについて、paired は強い一致（90%のソースで0.1以内）を示している。偏差は主に、経験的なノイズ事前分布が十分に堅牢ではない、色・マグニチュード図内の人口が希薄な領域で発生する。
• 既知の連星の回収：
  • 食連星： paired は既知のKeplerおよびTESS食連星の約74%を回収しており、これはエッジオン（真横を向いた）の短周期系がRVノイズを介して最も検出されやすいという期待と一致している。
  • 分光連星： 本フレームワークは、SB9カタログからの既知の単一線分光連星（SB1）の約83%を回収しており、この特定のサンプルに対しては、RVで選択された連星の幾何学的配置（エッジオン）により、アストロメトリ（位置天文学）のみの手法（RUWE）を上回る性能を示す。
  • Gaia軌道解： paired は、Gaiaの nss_two_body_orbit テーブルにある約30万のソースの約75%を特定しており、特にSB1としてフラグが立てられたものについては98%の回収率を示している。
• 半振幅の回収： $p < 0.01$ のソースについて、推論された半振幅は、高い信号対雑音比（S/N）において地上観測（APOGEE, SB9）と1-$\sigma$以内で一致する。精度は信号がGaiaのノイズフロア（$\sim$1 km/s）に近づくにつれて低下し、暗いソースについては、既知のGaia RV較正の問題と一致するわずかな系統的オフセット（$\sim$500 m/s）が認められる。
• 偽陽性： 著者らは、「過剰な」ノイズが非連星ソース（恒星活動：脈動、対流、高速回転など）からも生じ得ることを認めている。活動的なM型矮星や高速回転星はより頻繁にフラグが立てられており、これは paired のp値が、連星性と固有の恒星変動の両方を反映していることを示唆している。

意義と主張
本論文は、paired を、特に系外惑星科学や恒星集団研究のための、Gaiaコミュニティにとって必要な補助ツールとして位置づけている。著者らは、paired が以下のことを可能にすると主張している：

1. 集団レベルの解析： 地上でのフォローアップ観測に伴う選択バイアスを克服し、惑星を保有する星と保有しない星の間での連星率の一様な比較を可能にする。
2. 感度の推定： 分解されていない連星の影響を補正するために必要な「完全性」曲線を提供する。
3. 拡張された診断： DR3において特定の rv_chisqu_pvalue 製品を持たない大部分のGaiaソースに対して、連星のフラグを提供する。

著者らは「過剰な」ノイズの解釈について慎重であり、本フレームワークが統計的に有意な偏差を特定するものの、これらの偏差が排他的に連星のシグネチャであるわけではないことを明示している。彼らは、このツールが追加データなしに個々のシステムを確定的に確認するためのものではなく、候補を特定し、連星集団に対する統計的な制約を提供するためのものであることを強調している。本研究は、連星性が惑星の形成と人口統計に与える影響に関する将来の研究の基礎となるものである。