White dwarf planets in star clusters: gravitational scattering versus mass-loss effects

技術要約：星団における白色矮星惑星：重力散乱と質量放出効果

問題提起
白色矮星（WD）は、惑星系の長期的な生存性と進化を理解するためのユニークな実験室として機能する。主系列星以降の質量放出が惑星の軌道特性を変化させ、デブリ円盤を撹乱すること（これにより大気の汚染が引き起こされる）は既知であるが、誕生時の星形成領域による環境の影響については、依然として制約が不十分である。具体的には、高密度な星形成領域内での力学的な相互作用（重力散乱や恒星のフライバイなど）が、恒星の質量放出によって駆動される軌道進化とどのように競合し、あるいはそれを修飾するのかが明らかになっていない。これまでの研究では、クラスター内の力学的相互作用、あるいは孤立した状態での恒星進化のいずれかをモデル化してきたが、星形成領域の形成からその前駆星が白色矮星になるまでの過程を、N体動力学進化と恒星進化の両面から自己整合的に結合させた研究はほとんど存在しない。

手法
著者らは、Starlab環境内のkiraインテグレータを用いて、星形成領域のN体シミュレーションを行った。シミュレーションは1 Gyr（10億年）間実行された。これは、最も質量の大きい恒星（$>2 M_\odot$）が主系列を離れ、白色矮星になるのに十分なタイムスケールである。

• 初期条件: シミュレーションには、Maschberger (2013) の初期質量関数（IMF）から抽出された$N_\star = 1000$個の恒星を用いた。初期質量 $M_\star \leq 2.5 M_\odot$ のすべての恒星に対して、単一の木星質量を持つ惑星を割り当てた。初期軌道長半径（$a_p$）は、一様分布（0.1–50 au）または特定の値（0.1, 1, 5, 30 au）のデルタ関数からサンプリングされた。初期離心率および軌道傾斜角はゼロに設定された。
• クラスター力学: 2種類の初期クラスター構成をテストした：
  1. フラクタル分布: 高度に構造化された束縛のある星形成領域を表現するため、ボックス・フラクタル法を用いて次元 $D=1.6$ で生成された。初期半径は、1 pc（高密度、$\sim 10^4 M_\odot \text{pc}^{-3}$）および 5 pc（低密度、$\sim 10^2 M_\odot \text{pc}^{-3}$）に設定された。
  2. プラマー球: 半質量半径が 0.1 pc である、滑らかで中心集中型の分布。
• 恒星進化: 恒星および連星の進化は、N体ステップタイムと同期したSeBaルックアップテーブルを用いて実装された。コードは、巨星段階における質量放出を追跡した。恒星半径が惑星の軌道長半径を超えた場合、惑星は飲み込まれた（engulfed）とみなされた。なお、本シミュレーションには、原始的な恒星連星や外部の銀河潮汐場は含まれていない。
• シミュレーション・セット: 初期条件の確率性を評価するため、各セットに対して10回の実現（realization）を実行した。本研究は、白色矮星へと進化する前駆星（初期質量 1–2.5 $M_\odot$）を回る惑星の運命に焦点を当てた。

主な結果

1. 質量放出の支配性: 生存する惑星の軌道進化の主要な要因は、力学的散乱ではなく、恒星の質量放出である。高密度な環境は、散乱によって最大50%の惑星系を破壊または大幅に変形させるが、残りの系については質量放出の効果がダイナミクスを支配する。この挙動は、誕生領域の初期の恒星密度にはほとんど依存しない。
2. 軌道長半径の拡大: 飲み込みを生き延び、束縛状態を維持した惑星については、質量放出により軌道長半径が著しく増大する。断熱近似において、最終的な軌道長半径は最終的な恒星質量に反比例する（$a_f \propto M_{initial}/M_{final}$）。前駆星は通常50–80%の質量を失うため、惑星の軌道は2〜5倍に拡大する。この拡大は、最終的な軌道長半径分布において顕著な特徴（初期分布が0.1–50 auの場合、40–200 au付近にピークを持つ）を生み出す。これは、100 Myr時点での全惑星（主系列星を回るものを含む）の分布には見られない特徴である。
3. 離心率の変化: 全惑星の30%に対し、白色矮星の惑星の約40%が離心率の上昇を経験する。極めて重要な点は、これらの離心率の変化は、白色矮星相に移行する「前」に発生した力学的遭遇によってほぼ完全に引き起こされていることである。一度恒星が白色矮星になると、質量放出のみでは、$\sim 100$ au以内の惑星に対して無視できる程度の離心率変化しか誘発しない。
4. 捕獲および三重連星系:
   • 捕獲された惑星: 最終的な白色矮星惑星系の約10%は、捕獲された惑星（自由浮遊惑星が恒星に捕獲されたもの）で構成されている。これらの系は通常、大きな軌道長半径（$>100$ au）と高い離心率（$e > 0.5$）を持ち、熱的分布に似ている。
   • 三重連星系: 本シミュレーションは、白色矮星、惑星、および第三天体（恒星または捕獲された惑星）を含む長寿命の三重連星系を生成した。これらの構成は、PSR B1620-26 (AB) bのような系の類似系となり得る。
5. 破壊率: 高密度領域では、惑星系の破壊率が高い。白色矮星の前駆星を回る惑星の場合、1 Gyr後の束縛された生存惑星の割合は、全恒星の約71%に対し、約52%まで低下する。これは、より質量の大きい恒星がより高い相互作用率を経験することに起因する。

意義および結論
本論文は、誕生クラスターにおける白色矮星惑星のダイナミカルな環境が、生存している大多数の惑星の軌道特性に対しては、限定的な影響しか与えないと結論付けている。具体的には：

• 生存している系の約80%について、現在の軌道構成は、誕生クラスターでの散乱ではなく、主系列以降の質量放出によって決定されている。
• 残りの約20%については、力学的遭遇が軌道長半径を大きく変えている。
• 白色矮星の周囲に見られる広離心率・高離心率の惑星（例：WD 0806-661 b）は、クラスター内での捕獲起源を示唆している可能性があり、これは銀河系内のフィールドよりもクラスター内でより効率的に起こるプロセスである。
• 逆に、非常に近距離にある惑星（例：WD 1856+534 b）は、質量放出による拡大のみではなく、クラスターの初期段階で発生した散乱メカニズム（他の惑星や恒星との相互作用を含む可能性あり）とその後の潮汐による切り詰めを経て、現在の軌道に至ったと考えられる。

著者らは、誕生クラスターから出現する白色矮星の巨大惑星（1–100 au）の個体群合成を提供しており、これはナンシー・グレース・ローマン宇宙望遠鏡、ガイア、JWATといった将来の観測施設に対する制約を与えるものである。本研究は、誕生環境が破壊と捕獲を通じて個体群を形成する一方で、生存する系の根本的な軌道拡大は、恒星進化の強固な帰結であることを強調している。