Runaway OB stars within 1 kpc of the Sun

技術要約：太陽近傍1 kpc以内の逃走OB型星

問題提起
逃走星（Runaway stars）は、その誕生環境から放出された高速度の天体であり、星団の運動学的履歴や恒星放出を駆動するメカニズムに関する重要な知見を提供する。伝統的に、観測のしやすさからOB型星にこの用語が適用されることが多いが、逃走星の割合（OB型星のうち逃走星が占める割合）の推定値は、歴史的に約1.5%から50%超まで幅広い分散を示してきた。この不一致は、主に基礎となる恒星サンプル（多くの場合、等級制限や分光学的バイアスがある）の差異や、逃走星の定義および選択に使用される多様な手法に起因している。先行研究では、これらの相違を解決するための統一された、かつ体積完備な（volume-complete）逃走星のセンサスを提供することに苦慮してきた。

手法
本研究では、Quintana et al. (2025) によるカタログから導出された、太陽近傍1 kpc以内のOB型星の、新たに発表された体積完備なサンプルを利用している。最終的なサンプルは、40個のO型星と24,488個のB型星で構成される。解析は以下の手順で行われる：

1. 銀河回転モデリング： 著者らは、マルコフ連鎖モンテカルロ（MCMC）法を用いて、サンプルの観測された固有運動に対して銀河回転モデルを適合させる。このモデルは、銀河半径および銀道経度の関数として銀河中心速度を扱い、銀河半径に対する方位角速度（$V_{\phi}$）および動径方向速度（$V_R$）の一次微分をパラメータ化する。また、銀河中心に対する太陽の運動も考慮に入れている。
2. 残留速度の算出： モデルによって予測される速度を観測された速度から差し引くことで、銀河面（$\ell$）方向および垂直（$b$）方向の特異速度を導出する。
3. 統計的特性化： これらの残留速度の分布をマクスウェル・ボルツマン分布に適合させ、非逃走星集団の運動学的分散を特徴付ける。
4. 逃走星の選択基準： 逃走星を特定するために、2つの異なる手法を採用する：
   • 手法1（固定閾値）： 固定された2次元特異速度閾値 $23 \text{ km s}^{-1}$ を適用する。この値は、観測された確率密度とマクスウェル・ボルツマン分布の交点に基づいて選ばれており、ある星が逃走星である可能性が（確率的に）高くなる速度を表している。
   • 手法2（正規化分散）： 不確かさで正規化された特異速度が有意性の閾値（$E > 1$）を超える場合、その星を逃走星として識別する。これは、測定誤差に対するサンプルの分散に対して、有意に高い速度を持つ星を選択することを意味する。
5. モンテカルロ・シミュレーション： 測定の不確かさやモデルのパラメータ変動を考慮するため、500回のシミュレーションを行うモンテカルロ実験を実行する。ある星が逃走星である確率は、シミュレーション全体を通じてその星が選択基準を超える頻度に基づいて計算される。
6. 3次元校正： 2次元の結果は、APOGEEおよびLAMOSTから得られた、動径速度が利用可能な3,884個の星のサブセットを用いて校正され、3次元特異速度を計算し、2次元の閾値を3次元のマクスウェル・ボルツマン分布に対して検証する。

主な結果

• 逃走星の割合： 固定速度閾値（手法1）を用いると、O型星の逃走星の割合は 17.5% ($^{+0.1}_{-2.5}$) であり、B型星では 7.0% ($\pm 0.1$) であることが判明した。手法2では、O型星で10.0%、B型星で2.5%というより低い割合が得られた。
• 速度分布： サンプルの残留速度は、パラメータ $b \approx 11 \text{ km s}^{-1}$ (2D) および $b \approx 14 \text{ km s}^{-1}$ (3D) を持つマクスウェル・ボルツマン分布によって良好に記述される。3次元解析により、2次元の閾値 $23 \text{ km s}^{-1}$ が、約 $28 \text{ km s}^{-1}$ の3次元閾値に対応しており、完全な3次元的取り扱いのための妥当な近似であることが確認された。
• 個別の候補天体： 本研究は、固有運動に基づき、いくつかの著名な候補天体（例：$\mu$ Col, AE Aur, $\zeta$ Pup, 68 Cyg, $\lambda$ Cep）の逃走星としての状態を裏付けている。また、動径速度のみに基づいて逃走星とみなされていた星（例：$\xi$ Per, X Per）が、これらの基準の下では、逃走星として分類するには不十分な2次元固有運動を持っていることを特定し、特定の天体に対する2次元のみの解析の限界を浮き彫りにした。
• 空間分布： 逃走星、特にB型星は、一般的なOB型集団と比較して、より拡散した空間分布を示す。O型星の逃走星は、その年齢の若さから、銀河面により近い位置に留まっている。逃走星の垂直分布は、一般的なOB型集団と比較して「薄い」円盤成分の重みが減少しており、これは逃走星がより高い銀河緯度で見つかりやすいことを示している。

意義および主張
著者らは、本研究の主要な意義は、1 kpc以内のOB型星の体積完備なサンプルを用いたことにあると主張している。等級制限や分光学的バイアスのあるサンプル（非埋没環境において観測しやすいため、逃走星を過大評価する可能性がある）に依存していた従来の研究とは異なり、この体積完備なアプローチは選択バイアスを軽減している。

結果は、O型星がB型星よりも有意に高い逃走星の割合を持つことを示しており、これは先行文献でも示唆されていた傾向であるが、今回の完全な局所サンプルによって定量化された。著者らは、O型星の高い逃走星の割合は、それらの放出の支配的なメカニズムとして**力学的放出シナリオ（Dynamical Ejection Scenario: DES）**を支持していることを示唆しており、一方でB型星の低い割合は、**連星超新星シナリオ（Binary Supernova Scenario: BSS）**の寄与が大きい可能性を示唆している。本論文は、導出された割合は一部の先行研究と一致しているものの、文献に見られる相違は、天体固有の天文学的な差異ではなく、主に基礎となるサンプルと選択手法に起因するものであると結論付けている。本研究は、これらの星をその誕生星団へと遡る将来の研究のための強固な基準を提供するものである。