Drifters on the edge of town: $λ$ Boötis stars in clusters

技術要約：街の端を漂う者たち：星団におけるλ Boötis型星

問題提起
λ Boötis（λ Boo）型星は、軽元素（C, N, O, S）については太陽と同程度の組成を示すが、鉄族元素が著しく欠乏している化学的特異性を持つA型星である。主要な理論的説明は、これらの星が原始的な金属欠乏ガスを星間物質（ISM）から降着しているという仮説に基づいている。しかし、このシナリオは重大な課題に直面している。最近、星団内またはその周辺部でλ Boo型星が発見されたことは、そのような環境における中間質量および大質量星（>5 M⊙）からの強力な電離放射線が、原始的なISMガスの降着を妨げたり、プロセスに必要な星周円盤を剥ぎ取ったりするという予想と矛盾するためである。さらに、λ Boo型の組成を生み出すために必要な降着率（10⁻¹⁴ to 10⁻⁹ M⊙ yr⁻¹）と、白色矮星の観測から推測されるはるかに低い降着率（<10⁻¹⁷ M⊙ yr⁻¹）との間には乖離が存在する。本研究が取り組む具体的な問題は、進化する星団におけるA型星の力学的履歴が、高放射線量の星団コアから脱出し、原始ガスを降着し、その後再び戻ってくることを可能にするのかどうかであり、それによって星団環境におけるλ Boo型星の存在とISM降着仮説を両立させられるかという点である。

手法
著者らは、純粋な重力N体シミュレーションを用いて、100 Myrにわたる星形成領域の形成と力学的進化をモデル化した。

• シミュレーションの設定: 全星数（$N_*$）を150、500、1000、3000とする4組のシミュレーションを実施した。初期条件には、ボックス・フラクタル法（フラクタル次元 $D=1.6$）によって生成された空間的および運動学的サブ構造、およびビリアル比 $\alpha_{vir} = 0.3$ を含み、これにより初期崩壊と中心集中型の星団形成を導いた。星質量は、上界50 M⊙のMaschberger (2013) の初期質量関数（IMF）に従ってサンプリングされた。原始連星は除外された。
• 後処理解析: シミュレーションは、以下の項目を決定するために事後解析された。
  1. 降着ポテンシャル: A型星を追跡し、それらがヤコビ（潮汐）半径（$r_J$）を超えて移動した時期を特定した（$r_J$ は、原始的なISMガスと遭遇する可能性のある境界として定義）。$r_J$ の外側にいるとき、静止したISM雲（$\rho = 10$ cm⁻³）に対する星の速度に基づき、ボンディ・ホイル・リトルトン（BHL）降着率を計算した。
  2. 放射環境: 大質量星（>5 M⊙ および >20 M⊙）からの距離に基づき、各A型星に照射される遠紫外（FUV）および極端紫外（EUV）フラックスを計算した。
  3. 円盤の生存: FRIEDグリッドの光蒸発モデルおよび解析的な処方箋を用い、A型星の星周円盤の質量損失をモデル化した。円盤の全質量損失が初期円盤質量（星質量の10%と仮定）を超えた場合、その円盤は「破壊された」とみなした。星がλ Boo型の形成候補となるのは、円盤を保持したままISM物質を降着した場合のみとした。
  4. 時間的制約: Type Ia 超新星によるISMの富化を考慮し、40 Myr 以後の降着を制限する制約を設けて計算を繰り返した。

主な結果

• 力学的履歴: すべてのシミュレーション実現例において、大多数（90–98%）のA型星が、100 Myrの進化過程のどこかでヤコビ半径を超えて移動した。これは、星団コアからの放出、あるいは「漂流（drifting）」、すなわち星が潮汐半径の外側へ移動した後に再び内側領域へと戻ってくることによって起こった。
• 降着率: ヤコビ半径を超えて移動したA型星は、$10^{-12}$ から $10^{-5}$ M⊙ yr⁻¹ の範囲のBHL降着率を達成した。高速で放出された星は低い降着率（$10^{-12}–10^{-10}$ M⊙ yr⁻¹）を示したが、星団中心に向かう軌道を描き、外側へ漂流しては戻ってきた星は高い降着率（$10^{-9}–10^{-5}$ M⊙ yr⁻¹）を示し、これはλ Boo型の組成を説明するために必要な範囲と重なった。
• 円盤の生存 vs 放射: 星周円盤の生存は、星団内の総星数（およびそれに伴う電離源の数）に強く依存していた。
  • 低質量領域（$N_* = 150$）では、ヤコビ半径の外側へ出たA型星の約55%が10 Myr以上円盤を保持していた。
  • 高質量領域（$N_* = 3000$）では、強烈なFUV/EUV場により、円盤を保持できたA型星は0%であり、これらの環境におけるλ Boo型形成は事実上不可能であった。
  • 中間質量領域（$N_* = 500, 1000$）では、一部のA型星（それぞれ8.8%および4.8%）が降着しながら円盤を保持していた。
• 連星系: 本研究では、星形成領域の解体中に捕獲を通じて形成され得る、A型星を含む広い連星系（>100 au）が存在することを見出した。決定的なことに、星は原始的なISM物質を降着してλ Boo型の組成を発達させた後に、伴星を捕獲する（あるいはその逆）ことが可能である。このメカニズムにより、構成成分が著しく異なる化学組成を持つ連星系が生じる。シミュレーションされたこのような連星の分離分布は、観測されているλ Boo型星を含む広い連星（例：HD 87304 + CD-33 6615B）と交差している。

意義および主張
本論文は、星団内またはその近傍にλ Boo型星が存在することが、必ずしもISM降着シナリオを否定するものではないことを示すと主張している。著者らは以下の通り論じている：

1. 力学的脱出は十分である: A型星は、原始ガスを降着するために永久に放出される必要はない。彼らは潮汐半径の外側へ漂流し、降着し、そして星団の中心へと戻ってくることができる。
2. 放射閾値: 降着シナリオの実行可能性は、星団の質量と強く相関している。低質量星団（Theia 139やHSC 1640のような、最近λ Boo型が見つかった場所）は、円盤の生存を許容するほど弱い放射場を持っているが、高質量星団はそうではない。
3. 連星の診断: 構成成分間で組成が異なる連星系の観測は、降着仮説を無効にするものではない。そのような系は、降着イベントの後に力学的に形成され得る。

著者らは、降着率の絶対的な妥当性については、白色矮星の限界値との間の緊張関係が未解決であることを指摘し、控えめな姿勢を保っている。結論として、彼らの力学的モデルは星団におけるλ Boo型星の形成の「可能性」を支持するものであるが、特定の降着率や、降着円盤と放射場の物理的な相互作用を検証するためには、さらなる研究が必要であるとしている。