Circumbinary disk formation through AGB star winds

技術要約：AGB星の風による周連星円盤の形成

問題提起
周連星円盤は、ポスト漸近巨星分枝（post-AGB）段階にある連星系の周囲で頻繁に観測されているが、その形成および長期的な進化を支配するメカニズムは依然として十分に理解されていない。これまでの研究では、共通外層フェーズや、半分離系におけるロッシュローブ・オーバーフロー（RLOF）を介した質量移動が探索されてきたが、風と連星の相互作用が長期にわたって円盤を生成する具体的な役割については、十分に特徴付けられていない。既存のシミュレーションは、多くの場合100軌道周期未満の実行に留まっており、円盤の形態の長期的（セキュラー）な進化を捉えるには限界がある。さらに、異なる質量移動レジーム（具体的にはボンディ・ホイル（BH）降着と風ロッシュローブ・オーバーフロー（WRLOF）の間）の遷移や、連星パラメータ（離心率および質量比）が結果として生じる円盤構造にどのように影響するかについては、系統的な調査が必要である。

手法
著者らは、AGB星の風と連星の相互作用をモデル化するために、phantomコードを用いた粒子法（SPH）シミュレーションを用いている。系はAG_B星（$M_1$）と伴星（$M_2$）で構成され、それぞれ定義された降着半径（$R_{s,1} = 1.0$ AU、$R_{s,2} = 0.4$ AU）を持つシンク粒子として扱われる。シミュレーションはAC Her系のパラメータに基づいているが、注入される粒子の比エネルギーに応じてスケールフリーとなるよう調整されている。

主な手法の特徴は以下の通りである：

• 風の注入: ガス粒子はAGB星から$R_1 = 1.1$ AUの地点で注入される。本研究では、比エネルギー比$\eta$によって定義される2つの注入レジームを調査している：
  • $\eta > 1$: 高速の風であり、すべての粒子は非束縛状態となる（BHレジーム）。
  • $\eta < 0$: 低速の風であり、粒子は束縛状態となり、WRLOFへとつながる。
  • $0 \leq \eta \leq 1$: 中間レジーム。
• 注入幾何学: 一様分布とガウス分布（赤道方向に偏った分布）の両方がテストされており、詳細な円盤研究においては、中平面付近の解像度を高めるためにガウスモデルが採用されている。
• 熱力学: 主に局所的な等温状態（$H/R = 0.1$）が用いられ、熱的仮定を検証するためにバイナリ等温状態を用いた比較シミュレーションも行われている。
• 粘性: 有効なシャクラ・スニャイエ・ビスコシティ・パラメータとして$\alpha \sim 0.1$が課されているが、本研究では、課された粘性よりも物理的な構造（衝撃波やスパイラル）が角運動量輸送を支配しているかどうかを明示的に分析している。
• 期間: シミュレーションは600連星軌道周期（$\sim4500$年）を超えて実行されており、先行研究の期間を大幅に上回っている。

主な結果
シミュレーションの結果、流出の形態は、風の速度と伴星の質量の相互作用によって決定されることが明らかになり、これらは比エネルギー比$\eta$によって分類される：

1. ボンディ・ホイル・レジーム（高速風／低質量伴星）: $\eta \gg 1$のとき、風は大部分が非束縛状態のままである。伴星が風を重力的に集中させ、持続的なスパイラル密度パターンを形成する。ガス流は半径方向の速度に支配されており、物質は安定した周連星円盤を形成することなく、連星系から直接脱出する傾向がある。
2. WRLOFレジーム（低速風／高質量伴星）: $\eta < 0$（または中間）のとき、低速で高密度な風がドナーのロッシュローブを満たす。伴星の周囲に周単星円盤が形成され、これにより実効的な捕獲断面積が大幅に増加する。この円盤は効率的な角運動量輸送を促進し、物質が外側のラグランジュ点（$L_2$）を通じて脱出し、一貫した周連星円盤へと蓄積することを可能にする。
3. 連星離心率（$e_b$）の影響:
   • 円または軽度の離心（$e_b = 0, 0.2$）: 円盤は小さな半径（$\sim10$ AU）で形成され、外側に向かって拡大する。円盤は中程度の離心率（$e_d \approx 0.2$）を維持するが、粘性散逸と連星のトルクが競合することにより、時間の経過とともに徐々に円形化していく。
   • 高い離心（$e_b = 0.4$）: 伴星が近日点においてAGB星に接近するため、大量の物質が直接、高離心率の周連星軌道へと注入される。これにより、大きな長半径（$\sim30$ AU）を持ち、高い離心率（$e_d \approx 0.6$）を有する円盤が形成される。
4. 質量比（$q$）の影響: 質量比が低い（伴星が軽い）ほど、風の捕獲および円盤形成の効率は低下する。結果として得られる周連星円盤は密度が低く成長も緩やかであるが、それでも拡張した構造を形成することができる。注目すべきは、最も低い$q$の場合において周伴星円盤が解像されていない場合でも、拡張したほぼ円形の周連星円盤が発達していることであり、これは広範な円盤形成に高い初期密度が厳密に必要ではないことを示唆している。
5. 角運動量輸送: 有効応力（$\alpha_{eff}$）の分析により、周連星円盤内の角運動量輸送は、連星と風の相互作用によって生成されるスパイラル構造と衝撃波によって支配されていることが示された。これらの物理的な応力は、課された人工粘性に匹敵するかそれを上回っており、これは課された粘性が円盤の進化を支配していないことを意味している。

意義と主張
本論文は、風と連星の相互作用、特にWRLOFレジームにおける相互作用が、共通外層フェーズを必要とせずに、AGBおよびポストAGB系で観測される多様な周連星円盤の形態を自然に生成できると主張している。本研究は以下を実証している：

• スパイラル構造から安定した円盤への遷移を捉え、円盤の離心率と広がりの長期的（セキュラー）な進化を観察するためには、長期的なシミュレーションが必要である。
• 連星パラメータ（離心率および質量比）が円盤の半径方向の広がりと離心率を敏感に制御しており、これが観測される多様なポストAGB円盤構造の潜在的な説明となる。
• 周連星円盤の角運動量は完全に伴星によって供給されるため、自然に共面的な円盤形成へと導かれる。

著者らは、自らのモデルは共面的な円盤を説明できる一方で、AC Herのような系で見られる極方向に整列した周連星円盤については説明できておらず、そのような極方向の角運動量の起源は未解決の課題として残されていることを明示的に述べている。また、熱力学の簡略化（等温仮定）やダストと放射の結合の無視が詳細な進化に影響を与える可能性があることも認めているが、グローバルな円盤形成メカニズム自体は堅牢であると予想される。